The plastomes of Lepismium cruciforme (Vell.) Miq and Schlumbergera truncata (Haw.) Moran reveal tribe-specific rearrangements and the first loss of the trnT-GGU gene in Cactaceae.

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. METHODS AND 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.

The primate-specific presence of interferon regulatory factor-5 pseudogene 1.

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.

High-throughput fitness experiments reveal specific vulnerabilities of human-adapted Salmonella during stress and infection.

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.

Pseudogene OCT4-pg5 upregulates OCT4B expression to promote bladder cancer progression by competing with miR-145-5p.

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.

Improved Genetic Characterization of Congenital Adrenal Hyperplasia by Long-Read Sequencing Compared with Multiplex Ligation-Dependent Probe Amplification Plus Sanger Sequencing.

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.

The annotation of GBA1 has been concealed by its protein-coding pseudogene GBAP1.

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.

Probiotic functional gene explorations in the genome of Limosilactobacillus fermentum GD5MG.

Limosilactobacillus fermentum is an isolate obtained from oral gingival samples of healthy human individuals. The whole genome of Lb. fermentum GD5MG is composed of a circular DNA molecule containing 1,834,134 bp and exhibits a GC content of 52.80 %. The sequencing effort produced 38.6 million reads, each 150 bp in length, resulting in a sequencing depth of 2912.48x. Our examination unveiled a total of 1961 protein-coding genes, 27 rRNA genes, 24 tRNA genes, 3 non-coding RNA genes, and 63 pseudogenes with the use of gene annotations in NCBI Prokaryotic Genome Annotation tool. RAST revealed 1863 coding genes distributed across 209 subsystems, with a predominant involvement in amino acid, carbohydrate, and protein metabolism. Phylogenetic analysis infers that the Lb. fermentum GD5MG shares 281 gene clusters. Furthermore, the genome features showed a single CRISPR locus of 45 bp in length. Three genes associated with adhesion ability (strA, dltD, and dltA) and 26 genes related to acid tolerance, digestive enzyme secretion, and bile salt resistance were identified. Numerous genes associated with oral probiotic properties, comprising adhesion, acid and bile salt tolerance, oxidative stress tolerance, and sugar metabolism, were identified in the genome. Our findings shed light on the genomic characteristics of Lb. fermentum GD5MG, which are probable probiotics with functional benefits in humans.

Pseudogenes in plasmid genomes reveal past transitions in plasmid mobility.

Evidence for gene non-functionalization due to mutational processes is found in genomes in the form of pseudogenes. Pseudogenes are known to be rare in prokaryote chromosomes, with the exception of lineages that underwent an extreme genome reduction (e.g. obligatory symbionts). Much less is known about the frequency of pseudogenes in prokaryotic plasmids; those are genetic elements that can transfer between cells and may encode beneficial traits for their host. Non-functionalization of plasmid-encoded genes may alter the plasmid characteristics, e.g. mobility, or their effect on the host. Analyzing 10 832 prokaryotic genomes, we find that plasmid genomes are characterized by threefold-higher pseudogene density compared to chromosomes. The majority of plasmid pseudogenes correspond to deteriorated transposable elements. A detailed analysis of enterobacterial plasmids furthermore reveals frequent gene non-functionalization events associated with the loss of plasmid self-transmissibility. Reconstructing the evolution of closely related plasmids reveals that non-functionalization of the conjugation machinery led to the emergence of non-mobilizable plasmid types. Examples are virulence plasmids in Escherichia and Salmonella. Our study highlights non-functionalization of core plasmid mobility functions as one route for the evolution of domesticated plasmids. Pseudogenes in plasmids supply insights into past transitions in plasmid mobility that are akin to transitions in bacterial lifestyle.

Loss to gain: pseudogenes in microorganisms, focusing on eubacteria, and their biological significance.

Pseudogenes are defined as "non-functional" copies of corresponding parent genes. The cognition of pseudogenes continues to be refreshed through accumulating and updating research findings. Previous studies have predominantly focused on mammals, but pseudogenes have received relatively less attention in the field of microbiology. Given the increasing recognition on the importance of pseudogenes, in this review, we focus on several aspects of microorganism pseudogenes, including their classification and characteristics, their generation and fate, their identification, their abundance and distribution, their impact on virulence, their ability to recombine with functional genes, the extent to which some pseudogenes are transcribed and translated, and the relationship between pseudogenes and viruses. By summarizing and organizing the latest research progress, this review will provide a comprehensive perspective and improved understanding on pseudogenes in microorganisms. KEY POINTS: • Concept, classification and characteristics, identification and databases, content, and distribution of microbial pseudogenes are presented. • How pseudogenization contribute to pathogen virulence is highlighted. • Pseudogenes with potential functions in microorganisms are discussed.

Pseudogene CSPG4P12 inhibits colorectal cancer progression by attenuating epithelial-mesenchymal transition.

Colorectal cancer is one of the most common malignant cancers. Pseudogenes have been identified as oncogenes or tumor suppressor genes in the development of various cancers. However, the function of pseudogene CSPG4P12 in colorectal cancer remains unclear. Therefore, the aim of this study was to investigate the potential role of CSPG4P12 in colorectal cancer and explore the possible underlying mechanism. The difference of CSPG4P12 expression between colorectal cancer tissues and adjacent normal tissues was analyzed using the online Gene Expression Profiling Interactive Analysis 2 (GEPIA2) database. Cell viability and colony formation assays were conducted to evaluate cell viability. Transwell and wound healing assays were performed to assess cell migration and invasion capacities. Western blot was used to measure the expression levels of epithelial-mesenchymal transition-related proteins. Colorectal cancer tissues had lower CSPG4P12 expression than adjacent normal tissues. The overexpression of CSPG4P12 inhibited cell proliferation, invasion, and migration in colorectal cancer cells. Overexpressed CSPG4P12 promoted the expression of E-cadherin, whereas it inhibited the expression of vimentin, N-cadherin, and MMP9. These findings suggested that CSPG4P12 inhibits colorectal cancer development and may serve as a new potential target for colorectal cancer.

Conservation of a Chromosome 8 Inversion and Exon Mutations Confirm Common Gulonolactone Oxidase Gene Evolution Among Primates, Including H. Neanderthalensis.

Ascorbic acid functions as an antioxidant and facilitates other biochemical processes such as collagen triple helix formation, and iron uptake by cells. Animals which endogenously produce ascorbic acid have a functional gulonolactone oxidase gene (GULO); however, humans have a GULO pseudogene (GULOP) and depend on dietary ascorbic acid. In this study, the conservation of GULOP sequences in the primate haplorhini suborder were investigated and compared to the GULO sequences belonging to the primates strepsirrhini suborder. Phylogenetic analysis suggested that the conserved GULOP exons in the haplorhini primates experienced a high rate of mutations following the haplorhini/strepsirrhini divergence. This high mutation rate has decreased during the evolution of the haplorhini primates. Additionally, indels of the haplorhini GULOP sequences were conserved across the suborder. A separate analysis for GULO sequences and well-conserved GULOP sequences focusing on placental mammals identified an in-frame GULO sequence in the Brazilian guinea pig, and a potential GULOP sequence in the pika. Similar to haplorhini primates, the guinea pig and lagomorph species have experienced a high substitution rate when compared to the mammals used in this study. A shared synteny to examine the conservation of local genes near GULO/GULOP identified a conserved inversion around the GULO/GULOP locus between the haplorhini and strepsirrhini primates. Fischer's exact test did not support an association between GULOP and the chromosomal inversion. Mauve alignment showed that the inversion of the length of the syntenic block that the GULO/GULOP genes belonged to was variable. However, there were frequent rearrangements around ~ 2 million base pairs adjacent to GULOP involving the KIF13B and MSRA genes. These data may suggest that genes acquiring deleterious mutations in the coding sequence may respond to these deleterious mutations with rapid substitution rates.

Structural Features and Physiological Associations of Human 14-3-3ζ Pseudogenes.

There are about 14,000 pseudogenes that are mutated or truncated sequences resembling functional parent genes. About two-thirds of pseudogenes are processed, while others are duplicated. Although initially thought dead, emerging studies indicate they have functional and regulatory roles. We study 14-3-3ζ, an adaptor protein that regulates cytokine signaling and inflammatory diseases, including rheumatoid arthritis, cancer, and neurological disorders. To understand how 14-3-3ζ (gene symbol YWHAZ) performs diverse functions, we examined the human genome and identified nine YWHAZ pseudogenes spread across many chromosomes. Unlike the 32 kb exon-to-exon sequence in YWHAZ, all pseudogenes are much shorter and lack introns. Out of six, four YWHAZ exons are highly conserved, but the untranslated region (UTR) shows significant diversity. The putative amino acid sequence of pseudogenes is 78-97% homologous, resulting in striking structural similarities with the parent protein. The OMIM and Decipher database searches revealed chromosomal loci containing pseudogenes are associated with human diseases that overlap with the parent gene. To the best of our knowledge, this is the first report on pseudogenes of the 14-3-3 family protein and their implications for human health. This bioinformatics-based study introduces a new insight into the complexity of 14-3-3ζ's functions in biology.

Many purported pseudogenes in bacterial genomes are bona fide genes.

BACKGROUND: Microbial genomes are largely comprised of protein coding sequences, yet some genomes contain many pseudogenes caused by frameshifts or internal stop codons. These pseudogenes are believed to result from gene degradation during evolution but could also be technical artifacts of genome sequencing or assembly. RESULTS: Using a combination of observational and experimental data, we show that many putative pseudogenes are attributable to errors that are incorporated into genomes during assembly. Within 126,564 publicly available genomes, we observed that nearly identical genomes often substantially differed in pseudogene counts. Causal inference implicated assembler, sequencing platform, and coverage as likely causative factors. Reassembly of genomes from raw reads confirmed that each variable affects the number of putative pseudogenes in an assembly. Furthermore, simulated sequencing reads corroborated our observations that the quality and quantity of raw data can significantly impact the number of pseudogenes in an assembler dependent fashion. The number of unexpected pseudogenes due to internal stops was highly correlated (R2 = 0.96) with average nucleotide identity to the ground truth genome, implying relative pseudogene counts can be used as a proxy for overall assembly correctness. Applying our method to assemblies in RefSeq resulted in rejection of 3.6% of assemblies due to significantly elevated pseudogene counts. Reassembly from real reads obtained from high coverage genomes showed considerable variability in spurious pseudogenes beyond that observed with simulated reads, reinforcing the finding that high coverage is necessary to mitigate assembly errors. CONCLUSIONS: Collectively, these results demonstrate that many pseudogenes in microbial genome assemblies are actually genes. Our results suggest that high read coverage is required for correct assembly and indicate an inflated number of pseudogenes due to internal stops is indicative of poor overall assembly quality.

Evidence That Aquaporin 11 (AQP11) in the Spiny Dogfish (Squalus acanthias) May Represent a Pseudogene.

Various attempts to amplify an AQP11 cDNA from tissues of the spiny dogfish (Squalus acanthias) were made. Two pairs of deoxy-inosine-containing degenerate primers were designed based on conserved amino acid sequences from an AQP11 alignment. These primers yielded some faint bands from gill cDNA that were sequenced. Blast searches with the sequences showed they were not AQP11. An elasmobranch AQP11 nucleotide sequence alignment was produced to identify conserved regions to make further degenerate primers. One primer pair produced a short 148 bp fragment showing particularly strong amplification in gill and intestine. It was sequenced and represented a piece of the AQP11 gene. However, as the fragment may have resulted from contaminating genomic DNA (in total RNA used to make cDNA), 5' and 3' RACE were performed to amplify the two ends of the putative cDNA. Furthermore, 5' and 3' RACE amplifications depend on the presence of a 5' cap nucleotide and a poly A tail, respectively on the putative AQP11 mRNA. Hence, successful amplification was only possible from cDNA and not genomic DNA. Nested RACE amplifications were performed using gill and intestinal RACE cDNA, but none of the DNA fragments sequenced were AQP11. Consequently, the spiny dogfish AQP11 gene may represent a pseudogene.

UBDP1 pseudogene and UBD network competitively bind miR­6072 to promote glioma progression.

Increasing evidence suggests that pseudogenes play crucial roles in various cancers, yet their functions and regulatory mechanisms in glioma pathogenesis remain enigmatic. In the present study, a novel pseudogene was identified, UBDP1, which is significantly upregulated in glioblastoma and positively correlated with the expression of its parent gene, UBD. Additionally, high levels of these paired genes are linked with a poor prognosis for patients. In the present study, clinical samples were collected followed by various analyses including microarray for long non­coding RNAs, reverse transcription­quantitative PCR, fluorescence in situ hybridization and western blotting. Cell lines were authenticated and cultured then subjected to various assays for proliferation, migration, and invasion to investigate the molecular mechanisms. Bioinformatic tools identified miRNA targets, and luciferase reporter assays validated these interactions. A tumor xenograft model in mice was used for in vivo studies. In vitro and in vivo studies have demonstrated that UBDP1, localized in the cytoplasm, functions as a tumor­promoting factor influencing cell proliferation, migration, invasion and tumor growth. Mechanistic investigations have indicated that UBDP1 exerts its oncogenic effects by decoying miR­6072 from UBD mRNA, thus forming a competitive endogenous RNA network, which results in the enhanced oncogenic activity of UBD. The present findings offered new insights into the role of pseudogenes in glioma progression, suggesting that targeting the UBDP1/miR­6072/UBD network may serve as a potential therapeutic strategy for glioma patients.

Pathological roles of miRNAs and pseudogene-derived lncRNAs in human cancers, and their comparison as prognosis/diagnosis biomarkers.

This review examines and compares the diagnostic and prognostic capabilities of miRNAs and lncRNAs derived from pseudogenes in cancer patients. Additionally, it delves into their roles in cancer pathogenesis. Both miRNAs and pseudogene-derived lncRNAs have undergone thorough investigation as remarkably sensitive and specific cancer biomarkers, offering significant potential for cancer detection and monitoring. . Extensive research is essential to gain a complete understanding of the precise roles these non-coding RNAs play in cancer, allowing the development of novel targeted therapies and biomarkers for improved cancer detection and treatment approaches.

Construction and analysis of pseudogene-related ceRNA network in breast cancer.

Breast cancer (BC) is one of the leading causes of cancer-related deaths in women. The present study explored the potential role of pseudogenes in BC via construction and analysis of a competing endogenous RNA (ceRNA) network through a three-step process. First, we screened differentially expressed genes in nine BC datasets. Then the gene-pseudogenes pairs (nine hub genes) were selected according to the functional enrichment and correlation analysis. Second, the candidate hub genes and interacting miRNAs were used to construct the ceRNA network. Further analysis of the ceRNA network revealed a crucial ceRNA module with two genes-pseudogene pairs and two miRNAs. The in-depth analysis identified the GBP1/hsa-miR-30d-5p/GBP1P1 axis as a potential tumorigenic axis in BC patients. In the third step, the GBP1/hsa-miR-30d-5p/GBP1P1 axis expression level was assessed in 40 tumor/normal BC patients and MCF-7 cell lines. The expression of GBP1 and GBP1P1 was significantly higher in the tumor compared to the normal tissue. However, the expression of hsa-miR-30d-5p was lower in tumor samples. Then, we introduced the GBP1P1 pseudogene into the MCF-7 cell line to evaluate its effect on GBP1 and hsa-miR-30d-5p expression. As expected, the GBP1 level increased while the hsa-miR-30d-5p level decreased in the GBP1P1-overexprsssing cell line. In addition, the oncogenic properties of MCF-7 (cell viability, clonogenicity, and migration) were improved after GBP1P1 overexpression. In conclusion, we report a ceRNA network that may provide new insight into the role of pseudogenes in BC development.

Prediction of Tumor Microenvironment Characteristics and Treatment Response in Lung Squamous Cell Carcinoma by Pseudogene OR7E47P-related Immune Genes.

OBJECTIVE: Pseudogenes are initially regarded as nonfunctional genomic sequences, but some pseudogenes regulate tumor initiation and progression by interacting with other genes to modulate their transcriptional activities. Olfactory receptor family 7 subfamily E member 47 pseudogene (OR7E47P) is expressed broadly in lung tissues and has been identified as a positive regulator in the tumor microenvironment (TME) of lung adenocarcinoma (LUAD). This study aimed to elucidate the correlation between OR7E47P and tumor immunity in lung squamous cell carcinoma (LUSC). METHODS: Clinical and molecular information from The Cancer Genome Atlas (TCGA) LUSC cohort was used to identify OR7E47P-related immune genes (ORIGs) by weighted gene correlation network analysis (WGCNA). Based on the ORIGs, 2 OR7E47P clusters were identified using non-negative matrix factorization (NMF) clustering, and the stability of the clustering was tested by an extreme gradient boosting classifier (XGBoost). LASSO-Cox and stepwise regressions were applied to further select prognostic ORIGs and to construct a predictive model (ORPScore) for immunotherapy. The Botling cohorts and 8 immunotherapy cohorts (the Samstein, Braun, Jung, Gide, IMvigor210, Lauss, Van Allen, and Cho cohorts) were included as independent validation cohorts. RESULTS: OR7E47P expression was positively correlated with immune cell infiltration and enrichment of immune-related pathways in LUSC. A total of 57 ORIGs were identified to classify the patients into 2 OR7E47P clusters (Cluster 1 and Cluster 2) with distinct immune, mutation, and stromal programs. Compared to Cluster 1, Cluster 2 had more infiltration by immune and stromal cells, lower mutation rates of driver genes, and higher expression of immune-related proteins. The clustering performed well in the internal and 5 external validation cohorts. Based on the 7 ORIGs (HOPX, STX2, WFS, DUSP22, SLFN13, GGCT, and CCSER2), the ORPScore was constructed to predict the prognosis and the treatment response. In addition, the ORPScore was a better prognostic factor and correlated positively with the immunotherapeutic response in cancer patients. The area under the curve values ranged from 0.584 to 0.805 in the 6 independent immunotherapy cohorts. CONCLUSION: Our study suggests a significant correlation between OR7E47P and TME modulation in LUSC. ORIGs can be applied to molecularly stratify patients, and the ORPScore may serve as a biomarker for clinical decision-making regarding individualized prognostication and immunotherapy.

Comprehensive Identification of Mitochondrial Pseudogenes (NUMTs) in the Human Telomere-to-Telomere Reference Genome.

Practices related to mitochondrial research have long been hindered by the presence of mitochondrial pseudogenes within the nuclear genome (NUMTs). Even though partially assembled human reference genomes like hg38 have included NUMTs compilation, the exhaustive NUMTs within the only complete reference genome (T2T-CHR13) remain unknown. Here, we comprehensively identified the fixed NUMTs within the reference genome using human pan-mitogenome (HPMT) from GeneBank. The inclusion of HPMT serves the purpose of establishing an authentic mitochondrial DNA (mtDNA) mutational spectrum for the identification of NUMTs, distinguishing it from the polymorphic variations found in NUMTs. Using HPMT, we identified approximately 10% of additional NUMTs in three human reference genomes under stricter thresholds. And we also observed an approximate 6% increase in NUMTs in T2T-CHR13 compared to hg38, including NUMTs on the short arms of chromosomes 13, 14, and 15 that were not assembled previously. Furthermore, alignments based on 20-mer from mtDNA suggested the presence of more mtDNA-like short segments within the nuclear genome, which should be avoided for short amplicon or cell free mtDNA detection. Finally, through the assay of transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) on cell lines before and after mtDNA elimination, we concluded that NUMTs have a minimal impact on bulk ATAC-seq, even though 16% of sequencing data originated from mtDNA.

Tracking the Diversity and Chromosomal Distribution of the Olfactory Receptor Gene Repertoires of Three Anurans Species.

Olfaction is a crucial capability for most vertebrates and is realized through olfactory receptors in the nasal cavity. The enormous diversity of olfactory receptors has been created by gene duplication, following a birth-and-death model of evolution. The olfactory receptor genes of the amphibians have received relatively little attention up to now, although recent studies have increased the number of species for which data are available. This study analyzed the diversity and chromosomal distribution of the OR genes of three anuran species (Engystomops pustulosus, Bufo bufo and Hymenochirus boettgeri). The OR genes were identified through searches for homologies, and sequence filtering and alignment using bioinformatic tools and scripts. A high diversity of OR genes was found in all three species, ranging from 917 in B. bufo to 1194 in H. boettgeri, and a total of 2076 OR genes in E. pustulosus. Six OR groups were recognized using an evolutionary gene tree analysis. While E. pustulosus has one of the highest numbers of genes of the gamma group (which detect airborne odorants) yet recorded in an anuran, B. bufo presented the smallest number of pseudogene sequences ever identified, with no pseudogenes in either the beta or epsilon groups. Although H. boettgeri shares many morphological adaptations for an aquatic lifestyle with Xenopus, and presented a similar number of genes related to the detection of water-soluble odorants, it had comparatively far fewer genes related to the detection of airborne odorants. This study is the first to describe the complete OR repertoire of the three study species and represents an important contribution to the understanding of the evolution and function of the sense of smell in vertebrates.