BACKGROUND: Recent studies have revealed atypical features in the plastomes of the family Cactaceae, the largest lineage of succulent species adapted to arid and semi-arid regions. Most plastomes sequenced to date are from short-globose and cylindrical cacti, while little is known about plastomes of epiphytic cacti. Published cactus plastomes reveal reduction and complete loss of IRs, loss of genes, pseudogenization, and even degeneration of tRNA structures. Aiming to contribute with new insights into the plastid evolution of Cactaceae, particularly within the tribe Rhipsalideae, we de novo assembled and analyzed the plastomes of Lepismium cruciforme and Schlumbergera truncata, two South American epiphytic cacti.
RESULTS: Our data reveal many gene losses in both plastomes and the first loss of functionality of the trnT-GGU gene in Cactaceae. The trnT-GGU is a pseudogene in L. cruciforme plastome and appears to be degenerating in the tribe Rhipsalideae. Although the plastome structure is conserved among the species of the tribe Rhipsalideae, with tribe-specific rearrangements, we mapped around 200 simple sequence repeats and identified nine nucleotide polymorphism hotspots, useful to improve the phylogenetic resolutions of the Rhipsalideae. Furthermore, our analysis indicated high gene divergence and rapid evolution of RNA editing sites in plastid protein-coding genes in Cactaceae.
CONCLUSIONS: Our findings show that some characteristics of the Rhipsalideae tribe are conserved, such as plastome structure with IRs containing only the ycf2 and two tRNA genes, structural degeneration of the trnT-GGU gene and ndh complex, and lastly, pseudogenization of rpl33 and rpl23 genes, both plastid translation-related genes.

Interferon regulatory factor 5 (IRF5) is a key transcription factor in inflammatory and immune responses, with its dysregulation linked to autoimmune diseases. Using bioinformatic approaches, including Basic Local Alignment Search Tool (BLAST) for sequence similarity searches, BLAST-Like Alignment Tool (BLAT) for genome-wide alignments, and several phylogenetics software, such as Multiple Alignment using Fast Fourier Transform (MAFFT), for phylogenetic analyses, we characterized the structure, origin, and evolutionary history of the human IRF5 pseudogene 1 (IRF5P1). Our analyses reveal that IRF5P1 is a chimeric processed pseudogene containing sequences derived from multiple sources, including IRF5-like sequences from disparate organisms. We find that IRF5P1 is specific to higher primates, likely originating through an ancient retroviral integration event approximately 60 million years ago. Interestingly, IRF5P1 resides within the triple QxxK/R motif-containing (TRIQK) gene, and its antisense strand is predominantly expressed as part of the TRIQK pre-messenger RNA (mRNA). Analysis of publicly available RNA-seq data suggests potential expression of antisense IRF5P1 RNA. We hypothesize that this antisense RNA may regulate IRF5 expression through complementary binding to IRF5 mRNA, with human genetic variants potentially modulating this interaction. The conservation of IRF5P1 in the primate lineage suggests its positive effects on primate evolution and innate immunity. This study highlights the importance of investigating pseudogenes and their potential regulatory roles in shaping lineage-specific immune adaptations.

BACKGROUND: China is the hotspot of global freshwater crab diversity, but their wild populations are facing severe pressures associated with anthropogenic factors, necessitating the need to map their taxonomic and genetic diversity and design conservation policies.
RESULTS: Herein, we sequenced the mitochondrial genome of a Chinese freshwater crab species Bottapotamon fukienense, and found that it is fragmented into two chromosomes. We confirmed that fragmentation was not limited to a single specimen or population. Chromosome 1 comprised 15,111 base pairs (bp) and there were 26 genes and one pseudogene (pseudo-nad1) encoded on it. Chromosome 2 comprised 8,173 bp and there were 12 genes and two pseudogenes (pseudo-trnL2 and pseudo-rrnL) encoded on it. Combined, they comprise the largest mitogenome (23,284 bp) among the Potamidae. Bottapotamon was the only genus in the Potamidae dataset exhibiting rearrangements of protein-coding genes. Bottapotamon fukienense exhibited average rates of sequence evolution in the dataset and did not differ in selection pressures from the remaining Potamidae.
CONCLUSIONS: This is the first experimentally confirmed fragmentation of a mitogenome in crustaceans. While the mitogenome of B. fukienense exhibited multiple signs of elevated mitogenomic architecture evolution rates, including the exceptionally large size, duplicated genes, pseudogenisation, rearrangements of protein-coding genes, and fragmentation, there is no evidence that this is matched by elevated sequence evolutionary rates or changes in selection pressures.

BACKGROUND: Association testing between molecular phenotypes and genomic variants can help to understand how genotype affects phenotype. RNA sequencing provides access to molecular phenotypes such as gene expression and alternative splicing while DNA sequencing or microarray genotyping are the prevailing options to obtain genomic variants.
RESULTS: We genotype variants for 74 male Braunvieh cattle from both DNA (~ 13-fold coverage) and deep total RNA sequencing from testis, vas deferens, and epididymis tissue (~ 250 million reads per tissue). We show that RNA sequencing can be used to identify approximately 40% of variants (7-10 million) called from DNA sequencing, with over 80% precision. Within highly expressed coding regions, over 92% of expected variants were called with nearly 98% precision. Allele-specific expression and putative post-transcriptional modifications negatively impact variant genotyping accuracy from RNA sequencing and contribute to RNA-DNA differences. Variants called from RNA sequencing detect roughly 75% of eGenes identified using variants called from DNA sequencing, demonstrating a nearly 2-fold enrichment of eQTL variants. We observe a moderate-to-strong correlation in nominal association p-values (Spearman ρ2 ~ 0.6), although only 9% of eGenes have the same top associated variant.
CONCLUSIONS: We find hundreds of thousands of RNA-DNA differences in variants called from RNA and DNA sequencing on the same individuals. We identify several highly significant eQTL when using RNA sequencing variant genotypes which are not found with DNA sequencing variant genotypes, suggesting that using RNA sequencing variant genotypes for association testing results in an increased number of false positives. Our findings demonstrate that caution must be exercised beyond filtering for variant quality or imputation accuracy when analysing or imputing variants called from RNA sequencing.

Recently, substantial evidence has demonstrated that pseudogene-derived long noncoding RNAs (lncRNAs) as regulatory RNAs have been implicated in basic physiological processes and disease development through multiple modes of functional interaction with DNA, RNA, and proteins. Here, we report an important role for GBP1P1, the pseudogene of guanylate-binding protein 1, in regulating influenza A virus (IAV) replication in A549 cells. GBP1P1 was dramatically upregulated after IAV infection, which is controlled by JAK/STAT signaling. Functionally, ectopic expression of GBP1P1 in A549 cells resulted in significant suppression of IAV replication. Conversely, silencing GBP1P1 facilitated IAV replication and virus production, suggesting that GBP1P1 is one of the interferon-inducible antiviral effectors. Mechanistically, GBP1P1 is localized in the cytoplasm and functions as a sponge to trap DHX9 (DExH-box helicase 9), which subsequently restricts IAV replication. Together, these studies demonstrate that GBP1P1 plays an important role in antagonizing IAV replication.IMPORTANCELong noncoding RNAs (lncRNAs) are extensively expressed in mammalian cells and play a crucial role as regulators in various biological processes. A growing body of evidence suggests that host-encoded lncRNAs are important regulators involved in host-virus interactions. Here, we define a novel function of GBP1P1 as a decoy to compete with viral mRNAs for DHX9 binding. We demonstrate that GBP1P1 induction by IAV is mediated by JAK/STAT activation. In addition, GBP1P1 has the ability to inhibit IAV replication. Importantly, we reveal that GBP1P1 acts as a decoy to bind and titrate DHX9 away from viral mRNAs, thereby attenuating virus production. This study provides new insight into the role of a previously uncharacterized GBP1P1, a pseudogene-derived lncRNA, in the host antiviral process and a further understanding of the complex GBP network.

Genetic analysis of congenital adrenal hyperplasia (CAH) has been challenging because of high homology between CYP21A2 and its pseudogene CYP21A1P. This study aimed to evaluate the clinical utility of long-read sequencing (LRS) in diagnosis of CAH attributable to 21-hydroxylase deficiency by comparing with multiplex ligation-dependent probe amplification plus Sanger sequencing. In this retrospective study, 69 samples, including 49 probands from 47 families with high-risk of CAH, were enrolled and blindly subjected to detection of CAH by LRS. The genotype results were compared with control methods, and discordant samples were validated by additional Sanger sequencing. LRS successfully identified biallelic variants of CYP21A2 in the 39 probands diagnosed as having CAH. The remaining 10 probands were not patients with CAH. Additionally, LRS directly identified two pathogenic single-nucleotide variations (SNVs; c.293-13C/A>G and c.955C>T) in the presence of interference caused by nearby insertions/deletions (indels). The cis-trans configuration of two or more SNVs and indels identified in 18 samples was directly determined by LRS without family analysis. Eight CYP21A1P/A2 or TNXA/B deletion chimeras, composed of five subtypes, were identified; and the junction sites were precisely determined. Moreover, LRS determined the exact genotype in two probands who had three heterozygous SNVs/indels and duplication, which could not be clarified by control methods. These findings highlight that LRS could assist in more accurate genotype imputation and more precise CAH diagnosis.

Mutations in GBA1 cause Gaucher disease and are the most important genetic risk factor for Parkinson\'s disease. However, analysis of transcription at this locus is complicated by its highly homologous pseudogene, GBAP1. We show that >50% of short RNA-sequencing reads mapping to GBA1 also map to GBAP1. Thus, we used long-read RNA sequencing in the human brain, which allowed us to accurately quantify expression from both GBA1 and GBAP1. We discovered significant differences in expression compared to short-read data and identify currently unannotated transcripts of both GBA1 and GBAP1. These included protein-coding transcripts from both genes that were translated in human brain, but without the known lysosomal function-yet accounting for almost a third of transcription. Analyzing brain-specific cell types using long-read and single-nucleus RNA sequencing revealed region-specific variations in transcript expression. Overall, these findings suggest nonlysosomal roles for GBA1 and GBAP1 with implications for our understanding of the role of GBA1 in health and disease.

UNASSIGNED: Pseudogenes are sequences that have lost the ability to transcribe RNA molecules or encode truncated but possibly functional proteins. While they were once considered to be meaningless remnants of evolution, recent researches have shown that pseudogenes play important roles in various biological processes. However, the studies of pseudogenes in the silkworm, an important model organism, are limited and have focused on single or only a few specific genes.
UNASSIGNED: To fill these gaps, we present a systematic genome-wide studies of pseudogenes in the silkworm.
UNASSIGNED: We identified the pseudogenes in the silkworm using the silkworm genome assemblies, transcriptome, protein sequences from silkworm and its related species. Then we used transcriptome datasets from 832 RNA-seq analyses to construct spatio-temporal expression profiles for these pseudogenes. Additionally, we identified tissue-specifically expressed and differentially expressed pseudogenes to further understand their characteristics. Finally, the functional roles of pseudogenes as lncRNAs were systematically analyzed.
UNASSIGNED: We identified a total of 4410 pseudogenes, which were grouped into 4 groups, including duplications (DUPs), unitary pseudogenes (Unitary), processed pseudogenes (retropseudogenes, RETs), and fragments (FRAGs). The most of pseudogenes in the domestic silkworm were generated before the divergence of wild and domestic silkworm, however, the domestication may also involve in the accumulation of pseudogenes. These pseudogenes were clearly divided into 2 cluster, a highly expressed and a lowly expressed, and the posterior silk gland was the tissue with the most tissue-specific pseudogenes (199), implying these pseudogenes may be involved in the development and function of silkgland. We identified 3299 lncRNAs in these pseudogenes, and the target genes of these lncRNAs in silkworm pseudogenes were enriched in the egg formation and olfactory function.
UNASSIGNED: This study replenishes the genome annotations for silkworm, provide valuable insights into the biological roles of pseudogenes. It will also contribute to our understanding of the complex gene regulatory networks in the silkworm and will potentially have implications for other organisms as well.

Bladder cancer (BC) is one of the most common malignant neoplasms worldwide. Competing endogenous RNA (ceRNA) networks may identify potential biomarkers associated with the progression and prognosis of BC. The OCT4-pg5/miR-145-5p/OCT4B ceRNA network was found to be related to the progression and prognosis of BC. OCT4-pg5 expression was significantly higher in BC cell lines than in normal bladder cells, with OCT4-pg5 expression correlating with OCT4B expression and advanced tumor grade. Overexpression of OCT4-pg5 and OCT4B promoted the proliferation and invasion of BC cells, whereas miR-145-5p suppressed these activities. The 3\' untranslated region (3\'UTR) of OCT4-pg5 competed for miR-145-5p, thereby increasing OCT4B expression. In addition, OCT4-pg5 promoted epithelial-mesenchymal transition (EMT) by activating the Wnt/β-catenin pathway and upregulating the expression of matrix metalloproteinases (MMPs) 2 and 9 as well as the transcription factors zinc finger E-box binding homeobox (ZEB) 1 and 2. Elevated expression of OCT4-pg5 and OCT4B reduced the sensitivity of BC cells to cisplatin by reducing apoptosis and increasing the proportion of cells in G1. The OCT4-pg5/miR-145-5p/OCT4B axis promotes the progression of BC by inducing EMT via the Wnt/β-catenin pathway and enhances cisplatin resistance. This axis may represent a therapeutic target in patients with BC.

Salmonella enterica is comprised of genetically distinct \'serovars\' that together provide an intriguing model for exploring the genetic basis of pathogen evolution. Although the genomes of numerous Salmonella isolates with broad variations in host range and human disease manifestations have been sequenced, the functional links between genetic and phenotypic differences among these serovars remain poorly understood. Here, we conduct high-throughput functional genomics on both generalist (Typhimurium) and human-restricted (Typhi and Paratyphi A) Salmonella at unprecedented scale in the study of this enteric pathogen. Using a comprehensive systems biology approach, we identify gene networks with serovar-specific fitness effects across 25 host-associated stresses encountered at key stages of human infection. By experimentally perturbing these networks, we characterize previously undescribed pseudogenes in human-adapted Salmonella. Overall, this work highlights specific vulnerabilities encoded within human-restricted Salmonella that are linked to the degradation of their genomes, shedding light into the evolution of this enteric pathogen.