Search

SLC22A10 is an orphan transporter with unknown substrates and function. The goal of this study is to elucidate its substrate specificity and functional characteristics. In contrast to orthologs from great apes, human SLC22A10, tagged with green fluorescent protein, is not expressed on the plasma membrane. Cells expressing great ape SLC22A10 orthologs exhibit significant accumulation of estradiol-17β-glucuronide, unlike those expressing human SLC22A10. Sequence alignments reveal a proline at position 220 in humans, which is a leucine in great apes. Replacing proline with leucine in SLC22A10-P220L restores plasma membrane localization and uptake function. Neanderthal and Denisovan genomes show proline at position 220, akin to modern humans, indicating functional loss during hominin evolution. Human SLC22A10 is a unitary pseudogene due to a fixed missense mutation, P220, while in great apes, its orthologs transport sex steroid conjugates. Characterizing SLC22A10 across species sheds light on its biological role, influencing organism development and steroid homeostasis.

Background: The Synotis (C. B. Clarke) C. Jeffrey & Y. L. Chen is an ecologically important genus of the tribe Senecioneae, family Asteraceae. Because most species of the genus bear similar morphology, traditional morphological identification methods are very difficult to discriminate them. Therefore, it is essential to develop a reliable and effective identification method for Synotis species. In this study, the complete chloroplast (cp.) genomes of four Synotis species, S. cavaleriei (H.Lév.) C. Jeffrey & Y.L. Chen, S. duclouxii (Dunn) C. Jeffrey & Y.L. Chen, S. nagensium (C.B. Clarke) C. Jeffrey & Y.L. Chen and S. erythropappa (Bureau & Franch.) C. Jeffrey & Y. L. Chen had been sequenced using next-generation sequencing technology and reported here. Results: These four cp. genomes exhibited a typical quadripartite structure and contained the large single-copy regions (LSC, 83,288 to 83,399 bp), the small single-copy regions (SSC, 18,262 to 18,287 bp), and the inverted repeat regions (IR, 24,837 to 24,842 bp). Each of the four cp. genomes encoded 134 genes, including 87 protein-coding genes, 37 tRNA genes, 8 rRNA genes, and 2 pseudogenes (ycf1 and rps19). The highly variable regions (trnC-GCA-petN, ccsA-psaC, trnE-UUC-rpoB, ycf1, ccsA and petN) may be used as potential molecular barcodes. The complete cp. genomes sequence of Synotis could be used as the potentially effective super-barcode to accurately identify Synotis species. Phylogenetic analysis demonstrated that the four Synotis species were clustered into a monophyletic group, and they were closed to the Senecio, Crassocephalum and Dendrosenecio in tribe Senecioneae. Conclusions: This study will be useful for further species identification, evolution, genetic diversity and phylogenetic studies within this genus Synotis and the tribe Senecioneae. [ABSTRACT FROM AUTHOR]

The genus Ainsliaea DC. is one of the major groups within the tribe Pertyeae (Asteraceae). It comprises several important Chinese medicinal species. However, the phylogenetic position has undergone a long process of exploration. The complete chloroplast (cp) genome sequences data has not been employed in species identification and phylogeny of Ainsliaea. In this study, the complete cp genomes of two Ainsliaea species (A. gracilis and A. henryi) were reported, followed by structural, comparative, and phylogenetic analyses within the tribe Peryteae. Both cp genomes displayed a typical quadripartite circular structure, with the LSC and SSC regions separated by the IR regions. The genomes were 152,959 (A. gracilis) and 152,805 (A. henryi) base pairs (bp) long, with a GC content of 37.6%. They were highly conserved, containing 134 genes, including 87 protein-coding genes, 37 tRNA genes, 8 rRNA genes, and 2 pseudogenes (rps19 and ycf1). Moreover, thirteen highly polymorphic regions (e.g., trnK-UUU, trnG-UCC, trnT-GGU, accD-psaI, and rpl22-rps19) were identified, indicating their potential as DNA barcodes. The phylogenetic analysis confirmed the placement of Ainsliaea in the tribe Pertyeae, revealing close relationships with the genera Myripnois and Pertya. In comparison with Ainsliaea, Myripnois was more closely related to Pertya. This study lays a theoretical foundation for future research on species identification, population genetics, resource conservation, and sustainable utilization within Ainsliaea and Pertyeae. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Azolla is the only plant with a co-evolving nitrogen-fixing (diazotrophic) cyanobacterial symbiont (cyanobiont), Nostoc azollae, resulting from whole-genome duplication (WGD) 80 million years ago in Azolla's ancestor. Additional genes from the WGD resulted in genetic, biochemical, and morphological changes in the plant that enabled the transmission of the cyanobiont to successive generations via its megaspores. The resulting permanent symbiosis and co-evolution led to the loss, downregulation, or conversion of non-essential genes to pseudogenes in the cyanobiont, changing it from a free-living organism to an obligate symbiont. The upregulation of other genes in the cyanobiont increased its atmospheric dinitrogen fixation and the provision of nitrogen-based products to the plant. As a result, Azolla can double its biomass in less than two days free-floating on fresh water and sequester large amounts of atmospheric CO2, giving it the potential to mitigate anthropogenic climate change through carbon capture and storage. Azolla's biomass can also provide local, low-cost food, biofertiliser, feed, and biofuel that are urgently needed as our population increases by a billion every twelve years. This paper integrates data from biology, genetics, geology, and palaeontology to identify the location, timing and mechanism for the acquisition of a co-evolving diazotrophic cyanobiont by Azolla's ancestor in the Late Cretaceous (Campanian) of North America. [ABSTRACT FROM AUTHOR]

Background: 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. Objective: To fill these gaps, we present a systematic genome-wide studies of pseudogenes in the silkworm. Methods: 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. Results: 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. Conclusions: 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. [ABSTRACT FROM AUTHOR]

Abstract 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.

Abstract Background Several studies suggested that cavefish populations of Astyanax mexicanus settled during the Late Pleistocene. This implies that the cavefish’s most conspicuous phenotypic changes, blindness and depigmentation, and more cryptic characters important for cave life, evolved rapidly. Results Using the published genomes of 47 Astyanax cavefish from la Cueva de El Pachón, El Sótano de la Tinaja, La Cueva Chica and El Sótano de Molino, we searched for putative loss-of-function mutations in previously defined sets of genes, i.e., vision, circadian clock and pigmentation genes. Putative non-functional alleles for four vision genes were identified. Then, we searched genome-wide for putative non-functional alleles in these four cave populations. Among 512 genes with segregating putative non-functional alleles in cavefish that are absent in surface fish, we found an enrichment in visual perception genes. Among cavefish populations, different levels of shared putative non-functional alleles were found. Using a subset of 12 genes for which putative loss-of-function mutations were found, we extend the analysis of shared pseudogenes to 11 cave populations. Using a subset of six genes for which putative loss-of-function mutations were found in the El Sótano del Toro population, where extensive hybridization with surface fish occurs, we found a correlation between the level of eye regression and the amount of putative non-functional alleles. Conclusions We confirm that very few putative non-functional alleles are present in a large set of vision genes, in accordance with the recent origin of Astyanax mexicanus cavefish. Furthermore, the genome-wide analysis indicates an enrichment of putative loss-of-function alleles in genes with vision-related GO-terms, suggesting that visual perception may be the function chiefly impacted by gene losses related to the shift from a surface to a cave environment. The geographic distribution of putative loss-of-function alleles newly suggests that cave populations from Sierra de Guatemala and Sierra de El Abra share a common origin, albeit followed by independent evolution for a long period. It also supports that populations from the Micos area have an independent origin. In El Sótano del Toro, the troglomorphic phenotype is maintained despite massive introgression of the surface genome.

In the dynamic landscape of hepatocellular carcinoma (HCC) or the liver cancer research, pseudogenes have emerged from the shadows of genetic obscurity to become central figures, significantly influencing the disease molecular development and clinical trajectory. This review explores a transformative shift in perspective, recognizing pseudogenes not as genetic remnants without function, but as critical regulators in the molecular underpinnings of HCC. Engaging in complex interactions such as microRNA sponging, gene expression modulation, and signaling pathway disruptions, pseudogenes orchestrate a part of the molecular complexity driving tumor genesis, progression, and drug resistance in the liver cancer. Their unique expression patterns in hepatoma tissues herald new opportunities for early HCC detection, offering insights into patient prognosis, and identifying novel targets for therapeutic intervention of this disease. Such advancements underscore the importance of pseudogenes in enriching our understanding and management of HCC, paving the way for more effective diagnostic strategies and targeted therapies in the ongoing battle against this challenging malignancy. [ABSTRACT FROM AUTHOR]

This document is an editorial published in the journal Frontiers in Plant Science. It discusses the analysis of plastomes, which are the genomes of plant chloroplasts, and related bioinformatics tools. The editorial emphasizes the importance of comprehensive studies on plastomes for plant phylogenetic reconstruction and molecular identification. It also mentions the challenges in assembling and analyzing plastomes with unique features and introduces a research topic on this subject, which includes ten research articles. The document suggests that future research should focus on studying plastomes in parasitic plants, saprophytic plants, and plants living in extreme environments. [Extracted from the article]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

This study aimed to understand the genetic and metabolic traits of a Lactiplantibacillus plantarum JS21 strain and its probiotic abilities through laboratory tests and computer analysis. L. plantarum JS21 was isolated from a traditional fermented food known as "Jiangshui" in Hanzhong city. In this research, the complete genetic makeup of JS21 was determined using Illumina and PacBio technologies. The JS21 genome consisted of a 3.423 Mb circular chromosome and five plasmids. It was found to contain 3023 protein-coding genes, 16 tRNA genes, 64 rRNA operons, 40 non-coding RNA genes, 264 pseudogenes, and six CRISPR array regions. The GC content of the genome was 44.53%. Additionally, the genome harbored three complete prophages. The evolutionary relationship and the genome collinearity of JS21 were compared with other L. plantarum strains. The resistance genes identified in JS21 were inherent. Enzyme genes involved in the Embden–Meyerhof–Parnas (EMP) and phosphoketolase (PK) pathways were detected, indicating potential for facultative heterofermentative pathways. JS21 possessed bacteriocins plnE/plnF genes and genes for polyketide and terpenoid assembly, possibly contributing to its antibacterial properties against Escherichia coli (ATCC 25922), Escherichia coli (K88), Staphylococcus aureus (CMCC 26003), and Listeria monocytogenes (CICC 21635). Furthermore, JS21 carried genes for Na+/H+ antiporters, F0F1 ATPase, and other stress resistance genes, which may account for its ability to withstand simulated conditions of the human gastrointestinal tract in vitro. The high hydrophobicity of its cell surface suggested the potential for intestinal colonization. Overall, L. plantarum JS21 exhibited probiotic traits as evidenced by laboratory experiments and computational analysis, suggesting its suitability as a dietary supplement. [ABSTRACT FROM AUTHOR]

Pantoea sp. strain MHSD4 is a bacterial endophyte isolated from the leaves of the medicinal plant Pellaea calomelanos. Here, we report on strain MHSD4 draft whole genome sequence and annotation. The draft genome size of Pantoea sp. strain MHSD4 is 4 647 677 bp with a G+C content of 54.2% and 41 contigs. The National Center for Biotechnology Information Prokaryotic Genome Annotation Pipeline tool predicted a total of 4395 genes inclusive of 4235 protein-coding genes, 87 total RNA genes, 14 non-coding (nc) RNAs and 70 tRNAs, and 73 pseudogenes. Biosynthesis pathways for naphthalene and anthracene degradation were identified. Putative genes involved in bioremediation such as copA, copD, cueO, cueR, glnGm, and trxC were identified. Putative genes involved in copper homeostasis and tolerance were identified which may suggest that Pantoea sp. strain MHSD4 has biotechnological potential for bioremediation of heavy metals. [ABSTRACT FROM AUTHOR]

The honey bee trypanosomatid parasite, Lotmaria passim, contains two genes that encode the flagellar calcium binding protein (FCaBP) through tandem duplication in its genome. FCaBPs localize in the flagellum and entire body membrane of L. passim through specific N-terminal sorting sequences. This finding suggests that this is an example of protein subcellular relocalization resulting from gene duplication, altering the intracellular localization of FCaBP. However, this phenomenon may not have occurred in Leishmania, as one or both of the duplicated genes have become pseudogenes. Multiple copies of the FCaBP gene are present in several Trypanosoma species and Leptomonas pyrrhocoris, indicating rapid evolution of this gene in trypanosomatid parasites. The N-terminal flagellar sorting sequence of L. passim FCaBP1 is in close proximity to the BBSome complex, while that of Trypanosoma brucei FCaBP does not direct GFP to the flagellum in L. passim. Deletion of the two FCaBP genes in L. passim affected growth and impaired flagellar morphogenesis and motility, but it did not impact host infection. Therefore, FCaBP represents a duplicated gene with a rapid evolutionary history that is essential for flagellar structure and function in a trypanosomatid parasite. Author summary: Protein subcellular relocalization (PSR) was proposed as a mechanism for the functional divergence and retention of duplicate genes. We explored this hypothesis using flagellar calcium binding protein (FCaBP) genes duplicated in many trypanosomatid parasites. FCaBPs localize in the flagellum and entire body membrane of honey bee trypanosomatid parasite, Lotmaria passim through specific N-terminal sorting sequences. However, this is not common with all trypanosomatid parasites since one or both of the duplicated genes have become pseudogenes in Leishmania. N-terminal flagellar sorting sequence of FCaBP1 interacts with BBSome complex, and thus they must have co-evolved in each trypanosomatid species. FCaBPs are essential for the normal growth, flagellar morphogenesis, and motility but not host infection of L. passim. Our findings demonstrate that duplicated FCaBP genes have undergone PSR by changing amino acids in the N-terminal end with some but not all trypanosomatid parasites. In the parasites with multiple FCaBP genes, they play critical roles for the growth control as well as flagellar structure and function. [ABSTRACT FROM AUTHOR]

The Tripterospermum, comprising 34 species, is a genus of Gentianaceae. Members of Tripterospermum are mostly perennial, entwined herbs with high medicinal value and rich in iridoids, xanthones, flavonoids, and triterpenes. However, our inadequate understanding of the differences in the plastid genome sequences of Tripterospermum species has severely hindered the study of their evolution and phylogeny. Therefore, we first analyzed the 86 Gentianae plastid genomes to explore the phylogenetic relationships within the Gentianae subfamily where Tripterospermum is located. Then, we analyzed six plastid genomes of Tripterospermum, including two newly sequenced plastid genomes and four previously published plastid genomes, to explore the plastid genomes' evolution and phylogenetic relationships in the genus Tripterospermum. The Tripterospermum plastomes have a quadripartite structure and are between 150,929 and 151,350 bp in size. The plastomes of Tripterospermum encoding 134 genes were detected, including 86 protein-coding genes (CDS), 37 transfer RNA (tRNA) genes, eight ribosomal RNA (rRNA) genes, and three pseudogenes (infA, rps19, and ycf1). The result of the comparison shows that the Tripterospermum plastomes are very conserved, with the total plastome GC content ranging from 37.70% to 37.79%. In repeat sequence analysis, the number of single nucleotide repeats (A/T) varies among the six Tripterospermum species, and the identified main long repeat types are forward and palindromic repeats. The degree of conservation is higher at the SC/IR boundary. The regions with the highest divergence in the CDS and the intergenic region (IGS) are psaI and rrn4.5-rrn5, respectively. The average pi of the CDS and the IGS are only 0.071% and 0.232%, respectively, indicating that the Tripterospermum plastomes are highly conserved. Phylogenetic analysis indicated that Gentianinae is divided into two clades, with Tripterospermum as a sister to Sinogeniana. Phylogenetic trees based on CDS and CDS + IGS combined matrices have strong support in Tripterospermum. These findings contribute to the elucidation of the plastid genome evolution of Tripterospermum and provide a foundation for further exploration and resource utilization within this genus. [ABSTRACT FROM AUTHOR]

Several species in the oomycete genus Peronosclerospora cause downy mildew on maize and can result in significant yield losses in Asia. Bio-surveillance of these pathogens is a high priority to prevent epidemics on maize in the United States and consequent damage to the US economy. The unresolved taxonomy and dearth of molecular resources for Peronosclerospora spp. hinder these efforts. P. sorghi is a pathogen of sorghum and maize with a global distribution, for which limited diversity has been detected in the southern USA. We characterized the genome, transcriptome, and mitogenome of an isolate, representing the US pathotype 6. The highly homozygous genome was assembled using 10× Genomics linked reads and scaffolded using Hi-C into 13 chromosomes. The total assembled length was 303.2 Mb, larger than any other oomycete previously assembled. The mitogenome was 38 kb, similar in size to other oomycetes, although it had a unique gene order. Nearly 20,000 genes were annotated in the nuclear genome, more than described for other downy mildew causing oomycetes. The 13 chromosomes of P. sorghi were highly syntenic with the 17 chromosomes of Peronospora effusa with conserved centromeric regions and distinct chromosomal fusions. The increased assembly size and gene count of P. sorghi is due to extensive retrotransposition, resulting in putative pseudogenization. Ancestral genes had higher transcript abundance and were enriched for differential expression. This study provides foundational resources for analysis of Peronosclerospora and comparisons to other oomycete genera. Further genomic studies of global Peronosclerospora spp. will determine the suitability of the mitogenome, ancestral genes, and putative pseudogenes for marker development and taxonomic relationships.

Abstract Introduction Pseudogenes have been implicated for their role in regulating cellular differentiation and organismal development. However, their role in promoting cancer-associated differentiation has not been well-studied. This study explores the tumour landscape of oesophageal carcinoma to identify pseudogenes that may regulate events of differentiation to promote oncogenic transformation. Materials and method De-regulated differentiation-associated pseudogenes were identified using DeSeq2 followed by ‘InteractiVenn’ analysis to identify their expression pattern. Gene expression dependent and independent enrichment analyses were performed with GSEA and ShinyGO, respectively, followed by quantification of cellular reprogramming, extent of differentiation and pleiotropy using three unique metrics. Stage-specific gene regulatory networks using Bayesian Network Splitting Average were generated, followed by network topology analysis. MEME, STREME and Tomtom were employed to identify transcription factors and miRNAs that play a regulatory role downstream of pseudogenes to initiate cellular reprogramming and further promote oncogenic transformation. The patient samples were stratified based on the expression pattern of pseudogenes, followed by GSEA, mutation analysis and survival analysis using GSEA, MAF and ‘survminer’, respectively. Results Pseudogenes display a unique stage-wise expression pattern that characterizes stage II (SII) ESCA with a high rate of cellular reprogramming, degree of differentiation and pleiotropy. Gene regulatory network and associated topology indicate high robustness, thus validating high pleiotropy observed for SII. Pseudogene-regulated expression of SOX2, FEV, PRRX1 and TFAP2A in SII may modulate cellular reprogramming and promote oncogenesis. Additionally, patient stratification-based mutational analysis in SII signifies APOBEC3A (A3A) as a potential hallmark of homeostatic mutational events of reprogrammed cells which in addition to de-regulated APOBEC3G leads to distinct events of hypermutations. Further enrichment analysis for both cohorts revealed the critical role of combinatorial expression of pseudogenes in cellular reprogramming. Finally, survival analysis reveals distinct genes that promote poor prognosis in SII ESCA and patient-stratified cohorts, thus providing valuable prognostic bio-markers along with markers of differentiation and oncogenesis for distinct landscapes of pseudogene expression. Conclusion Pseudogenes associated with the events of differentiation potentially aid in the initiation of cellular reprogramming to facilitate oncogenic transformation, especially during SII ESCA. Despite a better overall survival of SII, patient stratification reveals combinatorial de-regulation of pseudogenes as a notable marker for a high degree of cellular differentiation with a unique mutational landscape.

Abstract Background Buruli ulcer (BU) disease, caused by Mycobacterium ulcerans (MU), and characterized by necrotic ulcers is still a health problem in Africa and Australia. The genome of the bacterium has several pseudogenes due to recent evolutionary events and environmental pressures. Pseudogenes are genetic elements regarded as nonessential in bacteria, however, they are less studied due to limited available tools to provide understanding of their evolution and roles in MU pathogenicity. Results This study developed a bioinformatic pipeline to profile the pseudogenomes of sequenced MU clinical isolates from different countries. One hundred and seventy-two MU genomes analyzed revealed that pseudogenomes of African strains corresponded to the two African lineages 1 and 2. Pseudogenomes were lineage and location specific and African lineage 1 was further divided into A and B. Lineage 2 had less relaxation in positive selection than lineage 1 which may signify different evolutionary points. Based on the Gil-Latorre model, African MU strains may be in the latter stages of evolutionary adaption and are adapting to an environment rich in metabolic resources with a lower temperature and decreased UV radiation. The environment fosters oxidative metabolism and MU may be less reliant on some secondary metabolites. In-house pseudogenomes from Ghana and Cote d’Ivoire were different from other African strains, however, they were identified as African strains. Conclusion Our bioinformatic pipeline provides pseudogenomic insights to complement other whole genome analyses, providing a better view of the evolution of the genome of MU and suggest an adaptation model which is important in understanding transmission. MU pseudogene profiles vary based on lineage and country, and an apparent reduction in insertion sequences used for the detection of MU which may adversely affect the sensitivity of diagnosis.

Different species of toothed whales (Odontoceti) exhibit a variety of tooth forms and enamel types. Some odontocetes have highly prismatic enamel with Hunter-Schreger bands, whereas enamel is vestigial or entirely lacking in other species. Different tooth forms and enamel types are associated with alternate feeding strategies that range from biting and grasping prey with teeth in most oceanic and river dolphins to the suction feeding of softer prey items without the use of teeth in many beaked whales. At the molecular level, previous studies have documented inactivating mutations in the enamel-specific genes of some odontocete species that lack complex enamel. At a broader scale, however, it is unclear whether enamel complexity across the full diversity of extant Odontoceti correlates with the relative strength of purifying selection on enamel-specific genes. Here, we employ sequence alignments for seven enamel-specific genes (ACP4, AMBN, AMELX, AMTN, ENAM, KLK4, MMP20) in 62 odontocete species that are representative of all extant families. The sequences for 33 odontocete species were obtained from databases, and sequences for the remaining 29 species were newly generated for this study. We screened these alignments for inactivating mutations (e.g., frameshift indels) and provide a comprehensive catalog of these mutations in species with one or more inactivated enamel genes. Inactivating mutations are rare in Delphinidae (oceanic dolphins) and Platanistidae/Inioidea (river dolphins) that have higher enamel complexity scores. By contrast, mutations are much more numerous in clades such as Monodontidae (narwhal, beluga), Ziphiidae (beaked whales), Physeteroidea (sperm whales), and Phocoenidae (porpoises) that are characterized by simpler enamel or even enamelless teeth. Further, several higher-level taxa (e.g., Hyperoodon, Kogiidae, Monodontidae) possess shared inactivating mutations in one or more enamel genes, which suggests loss of function of these genes in the common ancestor of each clade. We also performed selection (dN/dS) analyses on a concatenation of these genes and used linear regression and Spearman's rank-order correlation to test for correlations between enamel complexity and two different measures of selection intensity (# of inactivating mutations per million years, dN/dS values). Selection analyses revealed that relaxed purifying selection is especially prominent in physeteroids, monodontids, and phocoenids. Linear regressions and correlation analyses revealed a strong negative correlation between selective pressure (dN/dS values) and enamel complexity. Stronger purifying selection (low dN/dS) is found on branches with more complex enamel and weaker purifying selection (higher dN/dS) occurs on branches with less complex enamel or enamelless teeth. As odontocetes diversified into a variety of feeding modes, in particular, the suction capture of prey, a reduced reliance on the dentition for prey capture resulted in the relaxed selection of genes that are critical to enamel development. [ABSTRACT FROM AUTHOR]

Messenger RNA polyadenylation, the process wherein the primary RNA polymerase II transcript is cleaved and a poly(A) tract added, is a key step in the expression of genes in plants. Moreover, it is a point at which gene expression may be regulated by determining the functionality of the mature mRNA. Polyadenylation is mediated by a complex (the polyadenylation complex, or PAC) that consists of between 15 and 20 subunits. While the general functioning of these subunits may be inferred by extending paradigms established in welldeveloped eukaryotic models, much remains to be learned about the roles of individual subunits in the regulation of polyadenylation in plants. To gain further insight into this, we conducted a survey of variability in the plant PAC. For this, we drew upon a database of naturally-occurring variation in numerous geographic isolates of Arabidopsis thaliana. For a subset of genes encoding PAC subunits, the patterns of variability included the occurrence of premature stop codons in some Arabidopsis accessions. These and other observations lead us to conclude that some genes purported to encode PAC subunits in Arabidopsis are actually pseudogenes, and that others may encode proteins with dispensable functions in the plant. Many subunits of the PAC showed patterns of variability that were consistent with their roles as essential proteins in the cell. Several other PAC subunits exhibit patterns of variability consistent with selection for new or altered function. We propose that these latter subunits participate in regulatory interactions important for differential usage of poly(A) sites. [ABSTRACT FROM AUTHOR]

Background: Buruli ulcer (BU) disease, caused by Mycobacterium ulcerans (MU), and characterized by necrotic ulcers is still a health problem in Africa and Australia. The genome of the bacterium has several pseudogenes due to recent evolutionary events and environmental pressures. Pseudogenes are genetic elements regarded as nonessential in bacteria, however, they are less studied due to limited available tools to provide understanding of their evolution and roles in MU pathogenicity. Results: This study developed a bioinformatic pipeline to profile the pseudogenomes of sequenced MU clinical isolates from different countries. One hundred and seventy-two MU genomes analyzed revealed that pseudogenomes of African strains corresponded to the two African lineages 1 and 2. Pseudogenomes were lineage and location specific and African lineage 1 was further divided into A and B. Lineage 2 had less relaxation in positive selection than lineage 1 which may signify different evolutionary points. Based on the Gil-Latorre model, African MU strains may be in the latter stages of evolutionary adaption and are adapting to an environment rich in metabolic resources with a lower temperature and decreased UV radiation. The environment fosters oxidative metabolism and MU may be less reliant on some secondary metabolites. In-house pseudogenomes from Ghana and Cote d'Ivoire were different from other African strains, however, they were identified as African strains. Conclusion: Our bioinformatic pipeline provides pseudogenomic insights to complement other whole genome analyses, providing a better view of the evolution of the genome of MU and suggest an adaptation model which is important in understanding transmission. MU pseudogene profiles vary based on lineage and country, and an apparent reduction in insertion sequences used for the detection of MU which may adversely affect the sensitivity of diagnosis. [ABSTRACT FROM AUTHOR]

Background/aim: Chemical biological radiological nuclear threats are at an important point in the agenda of world health today, as they can cause mass deaths. B. mallei attracts attention as a potential biological warfare agent due to its features such as multidrug resistance, a rapid transmission mechanism via aerosol, the absence of a complete treatment protocol for the infection it causes, and the absence of an approved vaccine for protection against the bacteria. B. mallei suspect samples must be studied by experienced personnel in biosafety level III laboratories. B mallei is a difficult and troublesome pathogen to diagnose and many unknowns about B. mallei today. Therefore, the aim of the study was to determine the molecular differences and potential resistance genes of B mallei strains. Materials and methods: Determination of the molecular differences and potential resistance genes of B mallei strains with new bioinformatics approaches by comparatively examining the data of 29 B mallei strains, 10 of which were isolated from Türkiye, on the genome list of the National Biotechnology Information Center (NCBI). Results: According to the genome annotations of the origins, the origin containing the highest number of CDS which is 5172 was found as the 11th strain obtained in Türkiye in 1949. The origin with the highest number of pseudogenes was determined as 23,344 (China 7) origin. Two hundred and eighty-five pseudogenes found in this strain were obtained from a knee effusion in Myanmar. According to chromosome 2 data, B. mallei strain was determined as the most similar strain to ATCC 23344, line 11 with NCTC 10229 strain, and SAVP1 strain was determined as the least similar strain. When the antimicrobial resistance gene markers of the isolates included in the study were examined, amrA and amrB, qacG ade, Burkholderia pseudomallei Omp38 were found to be carrying. Conclusion: In terms of public health, it was thought that the data obtained as a result of our study about B mallei, which is defined as a biological weapon, is very valuable for creating treatment protocols to be applied to possible epidemics in the future. In addition, the available genetic epidemiological data of these strains belonging to a category that is dangerous to work with in a laboratory environment were reviewed. [ABSTRACT FROM AUTHOR]

The specific characteristics of k-mer words (2 ≤ k ≤ 11) regarding genomic distribution and evolutionary conservation were recently found. Among them are, in high abundance, words with a tandem repeat structure (repeat unit length of 1 bp to 3 bp). Furthermore, there seems to be a class of extremely short tandem repeats (≤12 bp), so far overlooked, that are non-random-distributed and, therefore, may play a crucial role in the functioning of the genome. In the following article, the positional distributions of these motifs we call super-short tandem repeats (SSTRs) were compared to other functional elements, like genes and retrotransposons. We found length- and sequence-dependent correlations between the local SSTR density and G+C content, and also between the density of SSTRs and genes, as well as correlations with retrotransposon density. In addition to many general interesting relations, we found that SINE Alu has a strong influence on the local SSTR density. Moreover, the observed connection of SSTR patterns to pseudogenes and -exons might imply a special role of SSTRs in gene expression. In summary, our findings support the idea of a special role and the functional relevance of SSTRs in the genome. [ABSTRACT FROM AUTHOR]

Objective: Although radiation workers are exposed to much lower doses of neutron-γ rays than those suffered in nuclear explosions and accidents, it does not mean that their health is not affected by radiation. Lower doses of radiation do not always cause morphological aberrations in chromosomes, so more sophisticated tests must be sought to specific alterations in the exposed cells. Our goal was to characterize the specific gene expression in lymphocytes from logging workers who were continuously exposed to low doses of neutron-γ radiation. We hypothesized that the combination of cell type-specific transcriptomes and open chromatin profiles would identify lymphocyte-specific gene alterations induced by long-term radiation with low-dose neutron^-rays and discover new regulatory pathways and transcriptional regulatory elements. Methods: Lymphocytes were extracted from workers who have been occupationally exposed to neutron-γ and workers unexposed to radiation in the same company. mRNA-seq and ATAC-seq (Assay for Transposase-Accessible Chromatin with high-throughput sequencing) were performed, followed integrative analysis to identify specific gene regulatory regions induced by neutron-γ radiation. A qPCR assay was then performed to verify the downregulation of RNA coding for ribosomal proteins and flow cytometry was used to detect ribosomal protein expression and cell cycle alterations. Results: We identified transcripts that were specifically induced by neutron-γ radiation and discovered differential open chromatin regions that correlated with these gene activation patterns. Notably, we observed a downward trend in the expression of both differentially expressed genes and open chromatin peaks. Our most significant finding was that the differential peak upregulated in ATAC-seq, while the differential gene was downregulated in the ribosome pathway. We confirmed that neutron-γ radiation leads to transcriptional inhibition by analyzing the most enriched promoters, examining RPS18 and RPS27A expression by qPCR, and analyzing protein-protein interactions of the differential genes. Ribosomal protein expression and cell cycle were also affected by neutron-γ as detected by flow cytometry. Conclusion: We have comprehensively analyzed the genetic landscape of human lymphocytes based on chromatin accessibility and transcript levels, enabling the identification of novel neutron-γ induced signature genes not previously known. By comparing fine-mapping of open chromatin and RNA reads, we have determined that neutron-γ specifically leads to downregulation of genes in the ribosome pathway, with pseudogenes potentially playing a crucial role. [ABSTRACT FROM AUTHOR]

Parnassius glacialis is a typical "Out of the QTP" alpine butterfly that originated on the Qinghai-Tibet Plateau (QTP) and dispersed into relatively low-altitude mountainous. Here we assemble a chromosome-level genome of P. glacialis and resequence 9 populations in order to explore the genome evolution and local adaptation of this species. These results indicated that the rapid accumulation and slow unequal recombination of transposable elements (TEs) contributed to the formation of its large genome. Several ribosomal gene families showed extensive expansion and selective evolution through transposon-mediated processed pseudogenes. Additionally, massive structural variations (SVs) of TEs affected the genetic differentiation of low-altitude populations. These low-altitude populations might have experienced a genetic bottleneck in the past and harbor genes with selective signatures which may be responsible for the potential adaptation to low-altitude environments. These results provide a foundation for understanding genome evolution and local adaptation for "Out of the QTP" of P. glacialis. The butterfly Parnassius glacialis experienced vast environmental transition when it dispersed out of the Qinghai-Tibet plateau. Here, the authors find that P. glacialis has an unusually large genome with rapid accumulation of transposable elements possibly facilitating its evolutionary adaptation. [ABSTRACT FROM AUTHOR]

"The Terpene Synthase Gene Family in Maize – A Clarification of Existing Community Nomenclature" is a document that explores the significance of terpenes in plants and other organisms, with a specific focus on terpene synthases (TPS). The authors analyze TPS genes in maize and identify approximately 40 TPS genes, including pseudogenes. They propose a standardized naming system for these genes to address gaps in the existing nomenclature. The document provides a comprehensive list of maize TPS genes and their corresponding identifiers, highlighting their roles in terpene biosynthesis. The authors also critique a previous study for its incomplete and inconsistent nomenclature. This document aims to provide clarity and consistency in the naming of maize TPS genes for future research. [Extracted from the article]

Most characterized metazoan mitochondrial genomes are compact and encode a small set of proteins that are essential for oxidative phosphorylation, as well as rRNA and tRNA for their expression. However, in rare cases, invertebrate taxa have additional open reading frames (ORFs) in their mtDNA sequences. Here, we sequenced and analyzed the mitochondrial genome of a polychaete worm, Polydora cf. ciliata , part of whose life cycle takes place in low-oxygen conditions. In the mitogenome, we found three "ORFan" regions (544, 1,060, and 427 bp) that have no resemblance to any standard metazoan mtDNA gene but lack stop codons in one of the reading frames. Similar regions are found in the mitochondrial genomes of three other Polydora species and Bocardiella hamata. All five species share the same gene order in their mitogenomes, which differ from that of other known Spionidae mitogenomes. By analyzing the ORFan sequences, we found that they are under purifying selection pressure and contain conservative regions. The codon adaptation indices (CAIs) of the ORFan genes were in the same range of values as the CAI of conventional protein-coding genes in corresponding mitochondrial genomes. The analysis of the P. cf. ciliata mitochondrial transcriptome showed that ORFan-544, ORFan-427, and a portion of the ORFan-1060 are transcribed. Together, this suggests that ORFan-544 and ORFan-427 encode functional proteins. It is likely that the ORFans originated when the Polydora/Bocardiella species complex separated from the rest of the Spionidae, and this event coincided with massive gene rearrangements in their mitochondrial genomes and tRNA-Met duplication. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Retrocopies are gene duplicates arising from reverse transcription of mature mRNA transcripts and their insertion back into the genome. While long being regarded as processed pseudogenes, more and more functional retrocopies have been discovered. How the stripped-down retrocopies recover expression capability and become functional paralogs continually intrigues evolutionary biologists. Here, we investigated the function and evolution of retrocopies in the context of 3D genome organization. By mapping retrocopy–parent pairs onto sequencing-based and imaging-based chromatin contact maps in human and mouse cell lines and onto Hi-C interaction maps in 5 other mammals, we found that retrocopies and their parental genes show a higher-than-expected interchromosomal colocalization frequency. The spatial interactions between retrocopies and parental genes occur frequently at loci in active subcompartments and near nuclear speckles. Accordingly, colocalized retrocopies are more actively transcribed and translated and are more evolutionarily conserved than noncolocalized ones. The active transcription of colocalized retrocopies may result from their permissive epigenetic environment and shared regulatory elements with parental genes. Population genetic analysis of retroposed gene copy number variants in human populations revealed that retrocopy insertions are not entirely random in regard to interchromosomal interactions and that colocalized retroposed gene copy number variants are more likely to reach high frequencies, suggesting that both insertion bias and natural selection contribute to the colocalization of retrocopy–parent pairs. Further dissection implies that reduced selection efficacy, rather than positive selection, contributes to the elevated allele frequency of colocalized retroposed gene copy number variants. Overall, our results hint a role of interchromosomal colocalization in the "resurrection" of initially neutral retrocopies. [ABSTRACT FROM AUTHOR]

ABSTRACTSalmonella enterica serovar Typhi and Paratyphi A are the cause of typhoid and paratyphoid fever in humans, which are systemic life-threatening illnesses. Both serovars are exclusively adapted to the human host, where they can cause life-long persistent infection. A distinct feature of these serovars is the presence of a relatively high number of degraded coding sequences coding for metabolic pathways, most likely a consequence of their adaptation to a single host. As a result of convergent evolution, these serovars shared many of the degraded coding sequences although often affecting different genes in the same metabolic pathway. However, there are several coding sequences that appear intact in one serovar while clearly degraded in the other, suggesting differences in their metabolic capabilities. Here, we examined the functionality of metabolic pathways that appear intact in S. Typhi but that show clear signs of degradation in S. Paratyphi A. We found that, in all cases, the existence of single amino acid substitutions in S. Typhi metabolic enzymes, transporters, or transcription regulators resulted in the inactivation of these metabolic pathways. Thus, the inability of S. Typhi to metabolize Glucose-6-Phosphate or 3-phosphoglyceric acid is due to the silencing of the expression of the genes encoding the transporters for these compounds due to point mutations in the transcriptional regulatory proteins. In contrast, its inability to utilize glucarate or galactarate is due to the presence of point mutations in the transporter and enzymes necessary for the metabolism of these sugars. These studies provide additional support for the concept of adaptive convergent evolution of these two human-adapted S. enterica serovars and highlight a limitation of bioinformatic approaches to predict metabolic capabilities.IMPORTANCESalmonella enterica serovar Typhi and Paratyphi A are the cause of typhoid and paratyphoid fever in humans, which are systemic life-threatening illnesses. Both serovars can only infect the human host, where they can cause life-long persistent infection. Because of their adaptation to the human host, these bacterial pathogens have changed their metabolism, leading to the loss of their ability to utilize certain nutrients. In this study we examined the functionality of metabolic pathways that appear intact in S. Typhi but that show clear signs of degradation in S. Paratyphi A. We found that, in all cases, the existence of single amino acid substitutions in S. Typhi metabolic enzymes, transporters, or transcription regulators resulted in the inactivation of these metabolic pathways. These studies provide additional support for the concept of adaptive convergent evolution of these two human-adapted S. enterica serovars.

The Ga1 locus controls cross-incompatibility between field corn and popcorn. The Ga1-S haplotype contains 2 types of pectin methylesterase (PME) genes, ZmPme3 and several copies of ZmGa1P that are expressed in silk and pollen, respectively. The ga1 haplotype contains nonfunctional tandem repeat sequences related to ZmPme3 and ZmGa1P. This haplotype can cross-pollinate freely and is widely present in field corn. The primary objective of this study is to characterize the repeat sequences from a diverse collection of maize and teosinte lines and use this information to understand the evolution of the Ga1 locus. First, we characterized the complexity of the Ga1 genome region in high-quality maize genome assemblies that led to their categorization into 5 groups based on the number and type of PME-like sequences found at this region. Second, we studied duplication events that led to the ga1 and Ga1-S repeats using maximum likelihood phylogenetic reconstruction. Divergence estimates of the ga1 haplotype suggest that the duplication events occurred more than 600 KYA whereas those in Ga1-S occurred at 3 time points, i.e. >600, ~260, and ~100 KYA. These estimates suggest that the ga1 and Ga1-S tandem duplication events occurred independently. Finally, analysis of ZmPme3 and ZmGa1P homologs in Zea and Tripsacum genomes suggests that ga1 and Ga1-S repeats originated from an ancestral pair of PME genes that duplicated and diverged through 2 evolutionary branches prior to the domestication of maize. [ABSTRACT FROM AUTHOR]

Olfactory receptor (OR) genes represent the largest multigenic family in mammalian genomes and encode proteins that bind environmental odorant molecules. The OR repertoire is extremely variable among species and is subject to many gene duplications and losses, which have been linked to ecological adaptations in mammals. Although they have been studied on a broad taxonomic scale (i.e. placental), finer sampling has rarely been explored in order to better capture the mechanisms that drove the evolution of the OR repertoire. Among placental mammals, rodents are well-suited for this task, as they exhibit diverse life history traits, and genomic data are available for most major families and a diverse array of lifestyles. In this study, 53 rodent published genomes were mined for their OR subgenomes. We retrieved more than 85,000 functional and pseudogene OR sequences that were subsequently classified into phylogenetic clusters. Copy number variation among rodents is similar to that of other mammals. Using our OR counts along with comparative phylogenetic approaches, we demonstrated that ecological niches such as diet, period of activity, and a fossorial lifestyle strongly impacted the proportion of OR pseudogenes. Within the OR subgenome, phylogenetic inertia was the main factor explaining the relative variations of the 13 OR gene families. However, a striking exception was a convergent 10-fold expansion of the OR family 14 among the phylogenetically divergent subterranean mole-rat lineages belonging to Bathyergidae and Spalacidae families. This study illustrates how the diversity of the OR repertoire has evolved among rodents, both shaped by selective forces stemming from species life history traits and neutral evolution along the rodent phylogeny. [ABSTRACT FROM AUTHOR]

Small nuclear RNAs (snRNAs) are structural and functional cores of the spliceosome. In metazoan genomes, each snRNA has multiple copies/variants, up to hundreds in mammals. However, the expressions and functions of each copy/variant in one organism have not been systematically studied. Focus on U1 snRNA genes, we investigated all five copies in Drosophila melanogaster using two series of constructed strains. Analyses of transgenic flies that each have a U1 promoter-driven gfp revealed that U1:21D is the major and ubiquitously expressed copy, and the other four copies have specificities in developmental stages and tissues. Mutant strains that each have a precisely deleted copy of U1-gene exhibited various extents of defects in fly morphology or mobility, especially deletion of U1:82Eb. Interestingly, splicing was changed at limited levels in the deletion strains, while large amounts of differentially-expressed genes and alternative polyadenylation events were identified, showing preferences in the down-regulation of genes with 1–2 introns and selection of proximal sites for 3'-end polyadenylation. In vitro assays suggested that Drosophila U1 variants pulled down fewer SmD2 proteins compared to the canonical U1. This study demonstrates that all five U1-genes in Drosophila have physiological functions in development and play regulatory roles in transcription and 3'-end formation. Author summary: Small nuclear RNAs (snRNAs) are the core of the spliceosome, which catalyzes the removal of introns. In metazoan genomes, each snRNA has multiple copies with variant sequences, up to hundreds in mammals, in which many singular copies/variants have shown function and are related to human diseases. However, expressions and functions of all snRNA copies in one organism have remained unknown. U1 snRNA is critical for early intron recognition and spliceosome assembly. The Drosophila genome has five U1 genes that generate three U1 variants. Here, we report that the U1:21D is ubiquitously high in all tissues and developmental stages, 95Ca is high in the salivary gland and male accessory gland, and 82Eb, 95Cb, and 95Cc are expressed in the ovary and embryo. Deletions of each singular U1-gene cause various developmental defects, demonstrating that none of them are pseudogenes. The absence of a U1-gene results in more serious defects in transcription and 3'-end formation than in pre-mRNA splicing, showing large amounts of expression-changed genes and alternative polyadenylation changes but much fewer changes in alternative splicing. Thus, our findings advance the understanding of in vivo functions of all the U1 snRNA genes in one organism. [ABSTRACT FROM AUTHOR]

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 [ABSTRACT FROM AUTHOR]

In the current study, we sequenced the complete plastome of Monstera deliciosa Liebm. (1849), an attractive foliage plant. The total length of the plastome of M. deliciosa is 163,499 bp and it consists of three distinct regions: a large single-copy (90,092 bp), a small single-copy (21,737 bp), and a pair of inverted repeats (IRs, 25,835 bp). The overall GC content is 36.2%, and the genome contains 110 functional genes, excluding pseudogenes. These functional genes encompass 77 protein-coding genes, 29 transfer RNA genes, and four ribosomal RNA genes. Notably, both the infA and rpl23 genes have been identified as pseudogenes. Phylogenetic analysis based on 14 representative plastomes from seven subfamilies indicates that Monsteroideae is monophyletic and sister to Pothoideae. Furthermore, M. deliciosa and M. adansonii were shown to share a recent common ancestor, the finding for which is supported by a strong bootstrap value. The sequenced plastome of M. deliciosa can serve as a valuable resource for establishing phylogenetic relationships and enhancing species identification within the genus Monstera. In addition, it can facilitate investigations into the genetic characteristics of this plant. [ABSTRACT FROM AUTHOR]

Luisia hancockii Rolfe 1896 is an epiphytic orchid species. In our present study, the whole chloroplast genome sequence of L. hancockii was de novo assembled by using high-throughput Illumina reads, and phylogenetic analysis was conducted within species of subtribe Aeridinae. The complete chloroplast genome sequence of L. hancockii was 146,243 bp in length, with a typical quadripartite structure, and its large single-copy, small single-copy, and inverted repeat were 84,441 bp, 11,412 bp, and 25,195 bp long, respectively. The GC content of the whole chloroplast genome was 36.6%, while the GC contents of LSC, SSC, and IR were 33.8%, 27.5%, and 43.3%, respectively. The chloroplast genome consisted of 129 genes, including 74 protein-coding genes, eight rRNAs, 38 tRNAs, and nine pseudogenes. Phylogenic tree was generated using the best model GTR + R, and the results showed that L. hancockii was sister to Holcoglossum and Vanda species, with a support of 100%. [ABSTRACT FROM AUTHOR]

The Elsholtzieae, comprising ca. 7 genera and 70 species, is a small tribe of Lamiaceae (mint family). Members of Elsholtzieae are of high medicinal, aromatic, culinary, and ornamentals value. Despite the rich diversity and value of Elsholtzieae, few molecular markers or plastomes are available for phylogenetics. In the present study, we employed high-throughput sequencing to assemble two Mosla plastomes, M. dianthera and M. scabra, for the first time, and compared with other plastomes of Elsholtzieae. The plastomes of Elsholtzieae exhibited a quadripartite structure, ranging in size from 148,288 bp to 152,602 bp. Excepting the absence of the pseudogene rps19 in Elsholtzia densa, the exhaustive tally revealed the presence of 132 genes (113 unique genes). Among these, 85 protein-coding genes (CDS), 37 tRNA genes, 8 rRNA genes, and 2 pseudogenes (rps19 and ycf1) were annotated. Comparative analyses showed that the plastomes of these species have minor variations at the gene level. Notably, the E. eriostchya plastid genome exhibited increased GC content regions in the LSC and SSC, resulting in an increased overall GC content of the entire plastid genome. The E. densa plastid genome displayed modified boundaries due to inverted repeat (IR) contraction. The sequences of CDS and intergenic regions (IGS) with elevated variability were identified as potential molecular markers for taxonomic inquiries within Elsholtzieae. Phylogenetic analysis indicated that four genera formed monophyletic entities, with Mosla and Perilla forming a sister clade. This clade was, in turn, sister to Collinsonia, collectively forming a sister group to Elsholtzia. Both CDS, and CDS + IGS could construct a phylogenetic tree with stronger support. These findings facilitate species identification and DNA barcoding investigations in Elsholtzieae and provide a foundation for further exploration and resource utilization within this tribe. [ABSTRACT FROM AUTHOR]

As an endosymbiont, Wolbachia exerts significant effects on the host, including on reproduction, immunity, and metabolism. However, the study of Wolbachia in Thysanopteran insects, such as flower thrips Frankliniella intonsa, remains limited. Here, we assembled a gap-free looped genome assembly of Wolbachia strain wFI in a length of 1,463,884 bp (GC content 33.80%), using Nanopore long reads and Illumina short reads. The annotation of wFI identified a total of 1838 protein-coding genes (including 85 pseudogenes), 3 ribosomal RNAs (rRNAs), 35 transfer RNAs (tRNAs), and 1 transfer-messenger RNA (tmRNA). Beyond this basic description, we identified mobile genetic elements, such as prophage and insertion sequences (ISs), which make up 17% of the entire wFI genome, as well as genes involved in riboflavin and biotin synthesis and metabolism. This research lays the foundation for understanding the nutritional mutualism between Wolbachia and flower thrips. It also serves as a valuable resource for future studies delving into the intricate interactions between Wolbachia and its host. [ABSTRACT FROM AUTHOR]

Anaplasma capra is an emerging tickborne human pathogen initially recognized in China in 2015; it has been reported in ticks and in a wide range of domestic and wild animals worldwide. We describe whole-genome sequences of 2 A. capra strains from metagenomic sequencing of purified erythrocytes from infected goats in China. The genome of A. capra was the smallest among members of the genus Anaplasma. The genomes of the 2 A. capra strains contained comparable G+C content and numbers of pseudogenes with intraerythrocytic Anaplasma species. The 2 A. capra strains had 54 unique genes. The prevalence of A. capra was high among goats in the 2 endemic areas. Phylogenetic analyses revealed that the A. capra strains detected in this study were basically classified into 2 subclusters with those previously detected in Asia. Our findings clarify details of the genomic characteristics of A. capra and shed light on its genetic diversity. [ABSTRACT FROM AUTHOR]

The complete chloroplast genomes of red currant cultivar 'Belaya Potapenko' and gooseberry cultivar 'Nekrasovskij' were sequenced and assembled for the first time. The plastomes are 157,802 bp and 157,559 bp in length for Ribes rubrum and R. uva-crispa, respectively. The R. rubrum cp genome is 243 b.p. longer. It has one more protein-coding gene ycf1, which is pseudogenized in the R. uva-crispa cp genome. In total, 56 and 54 simple sequence repeats (SSRs) were identified within the assembled plastid genomes. The SSR content of plastid genomes was assessed for the 18 Saxifragales species. Phylogeny inference based on plastome data of 18 Saxifragales revealed that all Ribes species are clustered together on the phylogenetic tree, though R. fasciculatum seems to be the most distant from the other analyzed Ribes species. The position of taxa inside the Ribes genus clade does not support the concept of its division into five subgenera. All Ribes species share approximately the same set of protein-coding genes in their plastome sequences. There was multiple independent pseudogenization of the ycf1 gene within the Ribes genus as well as other Saxifragales taxa. Negative selection was observed for most of the genes in both the Ribes group and Saxifragales. A positive selection ratio was observed only inside the Ribes group for the ycf4 and clpP genes. Together with positive selection signatures, pseudogenization events of ycfs genes perhaps reflect that these genes' evolution was important for Ribes' adaptation. Thus, our study provides genomic resources and valuable reference for marker development, and makes some clarifications of the phylogenomics of the Ribes genus. [ABSTRACT FROM AUTHOR]

The idea that proteins are the main determining factors in the functioning of cells and organisms, and their dysfunctions are the first cause of pathologies, has been predominant in biology and biomedicine until recently. This protein-centered view was too simplistic and failed to explain the physiological and pathological complexity of the cell. About 80% of the human genome is dynamically and pervasively transcribed, mostly as non-protein-coding RNAs (ncRNAs), which competitively interact with each other and with coding RNAs generating a complex RNA network regulating RNA processing, stability, and translation and, accordingly, fine-tuning the gene expression of the cells. Qualitative and quantitative dysregulations of RNA-RNA interaction networks are strongly involved in the onset and progression of many pathologies, including cancers and degenerative diseases. This review will summarize the RNA species involved in the competitive endogenous RNA network, their mechanisms of action, and involvement in pathological phenotypes. Moreover, it will give an overview of the most advanced experimental and computational methods to dissect and rebuild RNA networks.

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.

ABSTRACT Pseudogenes, once considered "junk DNA" based on the incorrect assumption that the absence of full coding potential means a complete lack of functionality, have recently become a subject of significant interest in the scientific community. Concurrently, it is widely assumed that bacterial genomes are compact and have a high density of coding genes with little room for non-coding genes, including pseudogenes. A key aspect of genome annotation is the correct identification of genes and the distinction between coding genes and pseudogenes, as it directly impacts functional and comparative genomics studies. In this study, we analyzed the genomic data of 4,699 strains of the bacterium Pseudomonas aeruginosa (P. aeruginosa) as they exhibit high variability in the number of annotated pseudogenes. In particular, we looked for correlations between the number of pseudogenes and other genomic and meta-features of the strains. We identified clusters of orthologous genes and pseudogenes and compared cluster size distributions and length homogeneity within clusters. We then mapped and examined orthology relationships between genes and pseudogenes. Additionally, we generated a phylogenetic tree of the strains and found that phylogenetically related strains are more homogeneous in the number of pseudogenes and share a significant amount of pseudogenes. Finally, we delved into clusters of orthologous genes and pseudogenes and quantified their phylogenetic neighborhood, classifying pseudogenes into evolutionary preserved pseudogenes, mis-annotated pseudogenes, or pseudogenes formed by failed horizontal transfer events. This in-depth study provides important insights that can be incorporated into pseudogene annotation pipelines in the future. IMPORTANCE Accurate annotation of genes and pseudogenes is vital for comparative genomics analysis. Recent studies have shown that bacterial pseudogenes have an important role in regulatory processes and can provide insight into the evolutionary history of homologous genes or the genome as a whole. Due to pseudogenes’ nature as non-functional genes, there is no commonly accepted definition of a pseudogene, which poses difficulties in verifying the annotation through experimental methods and resolving discrepancies among different annotation techniques. Our study introduces an in-depth analysis of annotated genes and pseudogenes and insights that can be incorporated into improved pseudogene annotation pipelines in the future.

Gibberellin (GAs) plays an important regulatory role in the development and growth of pineapple (Ananas comosus (L.) Merr.). Bioinformatics was used to confirm the differential expression of GA2 gibberellin oxidase gene AcGA2oxs in the pineapple genome, which laid the foundation for exploring its role in pineapple. In this study, 42 GA2ox genes (AcGA2oxs) were identified in the pineapple genome, named from AcGA2ox1 to AcGA2ox42, and divided into four groups according to phylogenetic analysis. We also analyzed the gene structure, conserved motifs and chromosome localization of AcGA2oxs. AcGA2oxs within the same group had similar gene structure and motifs composition. Collinear analysis and cis-element analysis provided the basis for understanding the evolution and function of GA2ox genes in pineapple. In addition, we selected different tissue parts to analyze the expression profile of AcGA2oxs, and the results show that 41 genes were expressed, except for AcGA2ox18. AcGA2ox18 may not be expressed in these sites or may be pseudogenes. qRT-PCR (real-time fluorescence quantitative PCR) was used to detect the relative expression levels of the GA2ox gene family under different concentrations of GA3 treatment, and it was found that AcGA2ox gene expression was upregulated in different degrees under GA3 treatment. These results provide useful information for further study on the evolution and function of the GA2ox family in pineapple. [ABSTRACT FROM AUTHOR]

Complex regions in the human genome such as repeat motifs, pseudogenes and structural (SVs) and copy number variations (CNVs) present ongoing challenges to accurate genetic analysis, particularly for short-read Next-Generation- Sequencing (NGS) technologies. One such region is the highly polymorphic CYP2D loci, containing CYP2D6, a clinically relevant pharmacogene contributing to the metabolism of >20% of common drugs, and two highly similar pseudogenes, CYP2D7 and CYP2D8. Multiple complex SVs, including CYP2D6/CYP2D7-derived hybrid genes are known to occur in different configurations and frequencies across populations and are difficult to detect and characterize accurately. This can lead to incorrect enzyme activity assignment and impact drug dosing recommendations, often disproportionally affecting underrepresented populations. To improve CYP2D6 genotyping accuracy, we developed a PCR-free CRISPR-Cas9 based enrichment method for targeted long-read sequencing that fully characterizes the entire CYP2D6- CYP2D7-CYP2D8 loci. Clinically relevant sample types, including blood, saliva, and liver tissue were sequenced, generating high coverage sets of continuous single molecule reads spanning the entire targeted region of up to 52 kb, regardless of SV present (n = 9). This allowed for fully phased dissection of the entire loci structure, including breakpoints, to accurately resolve complex CYP2D6 diplotypes with a single assay. Additionally, we identified three novel CYP2D6 suballeles, and fully characterized 17 CYP2D7 and 18 CYP2D8 unique haplotypes. This method for CYP2D6 genotyping has the potential to significantly improve accurate clinical phenotyping to inform drug therapy and can be adapted to overcome testing limitations of other clinically challenging genomic regions. [ABSTRACT FROM AUTHOR]

Objective. Skin cutaneous melanoma (SKCM) is a highly lethal malignancy that poses a significant threat to human health. Recent research has shown that competing endogenous RNA (ceRNA) regulatory networks play a critical role in the development and progression of various types of cancer, including SKCM. The objective of this study is to investigate the ceRNA regulatory network associated with the transmembrane protein semaphorin 6A (SEMA6A) and identify the underlying molecular mechanisms involved in SKCM. Methods. Expression profiles of four RNAs, including pseudogenes, long non-coding RNAs, microRNAs, and mRNAs were obtained from The Cancer Genome Atlas database. The analysis was completed by bioinformatics methods, and the expression levels of the selected genes were verified by cell experiments. Results. Bioinformatics analysis revealed that the LINC00511–hsa-miR-625-5p–SEMA6A ceRNA network was associated with SKCM prognosis. Furthermore, immune infiltration analysis indicated that the LINC00511–hsa-miR-625-5p–SEMA6A axis may have an impact on changes in the tumor immune microenvironment of SKCM. Conclusion. The LINC00511–hsa-miR-625-5p–SEMA6A axis could be a promising therapeutic target and a prognostic biomarker for SKCM. [ABSTRACT FROM AUTHOR]

Recent advances in long read technologies not only enable large consortia to aim to sequence all eukaryotes on Earth, but they also allow individual laboratories to sequence their species of interest with relatively low investment. Long read technologies embody the promise of overcoming scaffolding problems associated with repeats and low complexity sequences, but the number of contigs often far exceeds the number of chromosomes and they may contain many insertion and deletion errors around homopolymer tracts. To overcome these issues, we have implemented the ILRA pipeline to correct long read-based assemblies. Contigs are first reordered, renamed, merged, circularized, or filtered if erroneous or contaminated. Illumina short reads are used subsequently to correct homopolymer errors. We successfully tested our approach by improving the genome sequences of Homo sapiens , Trypanosoma brucei , and Leptosphaeria spp. and by generating four novel Plasmodium falciparum assemblies from field samples. We found that correcting homopolymer tracts reduced the number of genes incorrectly annotated as pseudogenes, but an iterative approach seems to be required to correct more sequencing errors. In summary, we describe and benchmark the performance of our new tool, which improved the quality of novel long read assemblies up to 1 Gbp. The pipeline is available at GitHub: https://github.com/ThomasDOtto/ILRA. [ABSTRACT FROM AUTHOR]

Background: Voltage-dependent anion selective channels (VDACs) are the most abundant mitochondrial outer membrane proteins, encoded in mammals by three genes, VDAC1, 2 and 3, mostly ubiquitously expressed. As 'mitochondrial gatekeepers', VDACs control organelle and cell metabolism and are involved in many diseases. Despite the presence of numerous VDAC pseudogenes in the human genome, their significance and possible role in VDAC protein expression has not yet been considered. Results: We investigated the relevance of processed pseudogenes of human VDAC genes, both in physiological and in pathological contexts. Using high-throughput tools and querying many genomic and transcriptomic databases, we show that some VDAC pseudogenes are transcribed in specific tissues and pathological contexts. The obtained experimental data confirm an association of the VDAC1P8 pseudogene with acute myeloid leukemia (AML). Conclusions: Our in-silico comparative analysis between the VDAC1 gene and its VDAC1P8 pseudogene, together with experimental data produced in AML cellular models, indicate a specific over-expression of the VDAC1P8 pseudogene in AML, correlated with a downregulation of the parental VDAC1 gene. [ABSTRACT FROM AUTHOR]

The vast majority of the human genome is non-coding. There is a diversity of non-coding features, some of which have functional importance. Although the non-coding regions constitute the majority of the genome, they remain understudied, and for a long time, these regions have been referred to as junk DNA. Pseudogenes are one of these features. A pseudogene is a non-functional copy of a protein-coding gene. Pseudogenes may arise through a variety of genetic mechanisms. Processed pseudogenes are formed through reverse transcription of mRNA by LINE elements, after which the cDNA is integrated into the genome. Processed pseudogenes are known to be variable across populations; however, the variability and distribution remains unknown. Herein, we apply a custom-designed processed pseudogene pipeline on the whole genome sequencing data of 3,500 individuals; 2,500 individuals from the thousand genomes dataset, as well as 1,000 Swedish individuals. Through these analyses, we discover over 3,000 pseudogenes missing from the GRCh38 reference. Utilising our pipeline, we position 74% of the detected processed pseudogenes—allowing for analyses of formation. Notably, we find that common structural variant callers, such as Delly, classify the processed pseudogenes as deletion events, which are later predicted to be truncating variants. By compiling lists of non-reference processed pseudogenes and their frequencies, we find a great variability of pseudogenes; indicating that non-reference processed pseudogenes may be useful for DNA testing and as population-specific markers. In summary, our findings highlight a great diversity of processed pseudogenes, that processed pseudogenes are actively formed in the human genome; and that our pipeline may be used to reduce false positive structural variation caused by the misalignment and subsequent misclassification of non-reference processed pseudogenes. [ABSTRACT FROM AUTHOR]