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]

Tinospora cordifolia is an important medicinal plant in ayurvedic system and modern medicine in India and traditionally applied for the treatment of inflammation, allergies and various other ailments. There is no genomic information available for T. cordifolia. In the present study, first time reported complete cp genome of T. cordifolia using high throughput ion torrent genome machine with ion torrent server. The cp genome of T. cordifolia was 150 945 bp in length with GC content of 38.7%, displays a quadripartite structure with two IR regions (IRa: 22 375 bp and IRb: 22 380 bp) that are separated by the LSC region (84 899 bp) and the SSC region (21 291 bp). There are 123 gene, including 85 protein-coding genes, 29 tRNA genes, seven rRNA genes and two pseudogenes (ycf1, rps16). Phylogenetic analysis based on whole cp genome revealed that T. cordifolia is as positioned near to the Menispermum dauricum within the Menispermaceae family. Total 27 SSRs were detected inthe plastomes which include di-, tri-, tetra- penta- and hexanucleotieds repeats. Further, compared with cp genome of other species of Menispermaceae it showed ten unique SSRs. Total twenty-five hypervariable regions loci were found in genome. These data could serve as DNA barcodes for species identification. [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.

Oreomecon nudicaulis, commonly known as mountain poppy, is a significant perennial herb. In 2022, the species O. nudicaulis, which was previously classified under the genus Papaver, was reclassified within the genus Oreomecon. Nevertheless, the phylogenetic status and chloroplast genome within the genus Oreomecon have not yet been reported. This study elucidates the chloroplast genome sequence and structural features of O. nudicaulis and explores its evolutionary relationships within Papaveraceae. Using Illumina sequencing technology, the chloroplast genome of O. nudicaulis was sequenced, assembled, and annotated. The results indicate that the chloroplast genome of O. nudicaulis exhibits a typical circular quadripartite structure. The chloroplast genome is 153,903 bp in length, with a GC content of 38.87%, containing 84 protein-coding genes, 8 rRNA genes, 38 tRNA genes, and 2 pseudogenes. The genome encodes 25,815 codons, with leucine (Leu) being the most abundant codon, and the most frequently used codon is AUU. Additionally, 129 microsatellite markers were identified, with mononucleotide repeats being the most abundant (53.49%). Our phylogenetic analysis revealed that O. nudicaulis has a relatively close relationship with the genus Meconopsis within the Papaveraceae family. The phylogenetic analysis supported the taxonomic status of O. nudicaulis, as it did not form a clade with other Papaver species, consistent with the revised taxonomy of Papaveraceae. This is the first report of a phylogenomic study of the complete chloroplast genome in the genus Oreomecon, which is a significant genus worldwide. This analysis of the O. nudicaulis chloroplast genome provides a theoretical basis for research on genetic diversity, molecular marker development, and species identification, enriching genetic information and supporting the evolutionary relationships among Papaveraceae. [ABSTRACT FROM AUTHOR]

Protein evolution is guided by structural, functional, and dynamical constraints ensuring organismal viability. Pseudogenes are genomic sequences identified in many eukaryotes that lack translational activity due to sequence degradation and thus over time have undergone "devolution." Previously pseudogenized genes sometimes regain their protein-coding function, suggesting they may still encode robust folding energy landscapes despite multiple mutations. We study both the physical folding landscapes of protein sequences corresponding to human pseudogenes using the Associative Memory, Water Mediated, Structure and Energy Model, and the evolutionary energy landscapes obtained using direct coupling analysis (DCA) on their parent protein families. We found that generally mutations that have occurred in pseudogene sequences have disrupted their native global network of stabilizing residue interactions, making it harder for them to fold if they were translated. In some cases, however, energetic frustration has apparently decreased when the functional constraints were removed. We analyzed this unexpected situation for Cyclophilin A, Profilin-1, and Small Ubiquitin-like Modifier 2 Protein. Our analysis reveals that when such mutations in the pseudogene ultimately stabilize folding, at the same time, they likely alter the pseudogenes' former biological activity, as estimated by DCA. We localize most of these stabilizing mutations generally to normally frustrated regions required for binding to other partners. [ABSTRACT FROM AUTHOR]

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

Objective: To investigate the expression of cutaneous T cell lymphoma associated antigen 15 (CTAGE15) and its association with clinical prognosis of esophageal squamous cell carcinoma (ESCC) patients, as well as its effect on the proliferation, migration, and invasion abilities of ESCC cells. Methods: Kaplan-Meier analysis was used to analyze the effect of CTAGE15 on the prognosis of ESCC patients in public database. The mRNA expression levels of CTAGE15 in esophageal epithelial cells and ESCC cells were measured by qRT⁃PCR. Furthermore, the expression of CTAGE15 was detected in cancer and adjacent normal tissues from 120 ESCC patients with a complete follow⁃up information. COX proportional survival risk model was enrolled to analyze the survival risk factors of ESCC patients, and multivariate correlation analysis was used for the analysis of the association between CTAGE15 expression and the clinicopathology of ESCC patients. Kaplan⁃Meier analysis was employed to analyze the effect of CTAGE15 on the prognosis of ESCC patients. The knockdown of CTAGE15 in ESCC cells was achieved by using small interfering RNA (siRNA), and the changes in cell proliferation, migration, and invasion abilities of ESCC cells were assessed. Results: The public database showed that the CTAGE15 expression was significantly associated with overall survival (P < 0.05) and disease⁃specific survival (P < 0.05) in ESCC patients, and a higher GTAGE15 expression was correlated with a better survival. The qRT-PCR results revealed significantly high expression of CTAGE15 in ESCC tissues and cells, compared with normal esophageal epithelial tissues and cells. In addition, the expression of CTAGE15 was significantly associated with T stage and N stage of ESCC patients, and a higher GTAGE15 expression was associated with an earlier T stage and N stage. Knockdown of CTAGE15 enhanced the proliferation, migration, and invasion abilities of ESCC cells (P < 0.05). Conclusion: CTAGE15 expression is positively associated with the prognosis of ESCC patients, and it inhibits the proliferation, migration and invasion abilities of ESCC cells. [ABSTRACT FROM AUTHOR]

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]

The mitochondrial genomes possess genes of five complexes involved in respiration. Mitochondrial gene content greatly varies across eukaryotes. Frequent gene losses have been reported in the mitogenomes of halophytes. In this study, we report, assemble, and annotate complete mitochondrial genome of a halophyte Pandanus odorifer (Forssk.) Kuntze (Pandanaceae). It is a commercially important species for its essential oil. It grows naturally in coastal regions and can withstand winds as well as salt sprays. The mitogenome of P. odorifer is circular with length of 330,962 bp. The genome contains 33 protein coding genes, 3 rRNAs, and 6 tRNAs. The analysis revealed that the important respiratory genes ccmC, cox1, nad1, nad4, and nad9 were pseudogenized in mitogenome but nuclear encoded copies of cox1, nad1, and nad9 were recovered. The phylogenetic analysis placed P. odorifer sister to P. tectorius. Further, mitogenome was compared with the other monocot mitogenomes for their gene contents. Thus, the study provides first complete mitochondrial genome for the family Pandanaceae and order Pandanales. [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]

The genes ECERIFERUM1 (CER1) and ECERIFERUM3 (CER3) encode the biosynthesis of alkane waxes, a key component of the plant cuticle. To study the evolution of CER1 and CER3 in a highly diverse group of eucalypts, we performed a genome-wide survey using recently released genome assemblies of 28 Myrtaceae species, with 22 species from the main eucalypt lineage and 6 non-eucalypt Myrtaceae tree species. We manually annotated 250 genes and pseudogenes, identifying a near-ubiquitous single copy of CER3 and 2 to 10 CER1 gene copies per Myrtaceae species. Phylogenetic analysis suggested that copy number variation in eucalypts is due to multiple tandem duplication events, both ancient (shared by all Myrtaceae species studied) and relatively recent (present only in eucalypts). Inter-chromosomal translocations were discovered for both CER1 and CER3, along with recurrent loss of often the same CER1 introns in the WAX2 domain, the domain that is essential for wax production. Despite the varied environments occupied by the eucalypt species in this study, we did not find statistically significant associations between intra-genic structural changes or CER1 copy number and aspects of the environment they occupy (including aridity). The challenge is now to explain the species-specific evolutionary histories that contributed to the observed variation in CER1 and the extent to which it may contribute to the adaptability of eucalypts. [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]

Comprehensive in silico analysis aimed to identify the genes involved in the biosafety and probiotic features of Limosilactobacillus fermentum J23 was performed. The strain had a single circular chromosome of 2,004,882 base pairs (bp), with a Guanine + Cytosine (G + C) content of 51.74%. The chromosome contained 2,001 genes predicted as protein-coding sequences (CDS), 86 pseudogenes, 56 tRNAs, five rRNAs, and three ncRNAs. The in silico analysis demonstrated the absence of risk-related genes corresponding to transferable/acquirable antibiotic resistance, virulence factors, toxic metabolites, and detrimental metabolic/enzymatic activities. The strain was predicted as nonpathogenic to humans and the presence of a variety of genes potentially related to its probiotic profile were revealed, including those of carbohydrate metabolism and transport, adhesion and aggregation ability, stress/adaptive response, CRISPR/Cas adaptive immunity, and immunomodulation. These findings highlight the probiotic potential of L. fermentum J23 and provide genomic clues of its biosafety to advance its future applications. [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]

Due to its faster evolution rate compared to nuclear genes, haploid mitochondrial DNA (mtDNA) is a promising species identification tool. This has led to its significant use in taxonomic and phylogenetic studies. While mtDNA is subject to selective constraints that prevent the accumulation of deleterious mutations, the prevalence of nuclear-mitochondrial fragments (known as NUMTs or pseudogenes) in mammals has complicated the use of mtDNA for taxonomy. In the present study, a pseudogene of the mitochondrial cytochrome b (Cytb) was detected in Meriones persicus. This pseudogene differed from its mitochondrial counterpart at 235 out of 1140 sites and is characterized by frame-shift mutations, indels, and accumulation of non-synonymous substitutions. It is the first report of the Cytb pseudogene in a jird, highlighting the risk of misidentifying NUMTs as authentic mtDNA and the importance of addressing this potential pitfall in taxonomic studies. [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]

Killer cell immunoglobulin-like receptor (KIR) molecules on the surface of Natural Killer cells HLA class I by interacting directly with their molecules, they cause inhibition or triggering of cytokine synthesis. The KIR gene family, plays an important role in the regulation of the immune response, is located in the 19q13.4 region in humans and consists of 15 homologous genes. Inhibitory KIR (iKIR) genes are 2DL1, 2DL2, 2DL3, 2DL4, 2DL5, 3DL1, 3DL2, 3DL3, activating KIR (aKIR) genes are KIR2DS1, KIR2DS2, KIR2DS3, KIR2DS4, KIR2DS5, KIR3DS1; two KIR pseudogenes are KIR2DP1 and KIR3DP1. There is not enough data on the KIR gene profile of the Turkish population. It was aimed to determine the genotype distrubition of KIR2DS1, KIR2DS2, KIR2DS3, KIR2DS4, KIR2DS5, KIR3DS1, KIR2DL1-5, KIR3DL1-3) genes in the Southeastern Anatolia Region (SAR). In the present study 200 unrelated healthy individuals were included. Polymerase-chain reaction with sequence-specific primers method (PCR-SSP) were performed. 3DL3 had the highest (99.5%) frequency and the lowest frequency among the inhibitory KIR genes was obtained for 2DL2 (56.5%). The highest frequency among the activator KIR genes was obtained for 2DS4 (93.5%) and the lowest was for 2DS5 (28%). It can be concluded that the SAR population has similar KIR gene features to the Caucasians. But studies with more extensive sample groups are needed to understand the relationship between KIR gene distribution and similarities between populations. [ABSTRACT FROM AUTHOR]

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]

Tandemly arranged repetitive regions (repeats) that encode 5S rRNA (5S rDNA) are an indispensable component of eukaryotic genomes. Typically, 5S rDNA repeats within a genome are very similar due to the concerted nature of the evolution of this type of repeats. Each 5S rDNA repeat consists of an evolutionarily conserved coding sequence (CDS) and a variable intergenic spacer (IGS). 5S rDNA is a popular model for studying the molecular evolution of repetitive sequences, and the high rate of IGS mutations determines its wide use in phylogenetic analysis of closely related taxa. Nevertheless, 5S rDNA remains unexplored for many groups of higher plants, especially for the Plumbaginaceae family. Some taxa of this family are endemic to southern Ukraine and listed in the Red Book. However, their taxonomic status is controversial, and its clarification requires the use of molecular phylogenetic methods. In this work, we examined the molecular organization of 5S rDNA for representatives of four genera of the tribe Limonieae, the largest in the family Plumbaginaceae. It was shown that the CDS of 5S rDNA of representatives of the genera Limonium, Armeria, and Ceratolimon possess single mutations that do not affect the formation of the secondary structure of 5S rRNA. In contrast, in the genomes of Goniolimon species, in addition to functionally normal 5S rDNA repeats, numerous pseudogenes were found that do not evolve in a concerted manner and contain numerous mutations in the CDS that disrupt the secondary structure of 5S rRNA. A significant phylogenetic distance between representatives of the subgenera Pteroclados and Limonium of the genus Limonium indicates that Pteroclados can be considered a separate genus. The high rate of molecular evolution makes 5S rDNA IGS a convenient tool for the reconstruction of phylogenetic relationships within the studied genera of the tribe Limonieae and the barcoding of Ukrainian endemics of the genus Limonium. [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.

Purpose: Hepatocellular carcinoma (HCC) remains a highly deadly malignant tumor with high recurrence and metastasis rates. Cancer stem cells (CSCs) are involved in tumor metastasis, recurrence, and resistance to drugs, which have attracted widespread attention in recent years. Research has shown that pseudogenes may regulate stemness to promote the progression of HCC, but its specific mechanisms and impact on prognosis remain unclear. Methods: In this study, clinical prognosis information of HCC was first downloaded from The Cancer Genome Atlas (TCGA) database. Then we calculated the mRNA expression-based stemness index (mRNAsi) of HCC. We also screened the differentially expressed pseudogene (DEPs) and conducted univariate Cox regression analysis to investigate their effect on the prognosis of HCC. Further, genomic mutation frequency analysis and weighted gene co-expression network analysis (WGCNA) were performed to compare the role of pseudogenes and stemness in promoting the progression of HCC. Finally, we conducted the correlation analysis to examine the potential mechanism of pseudogenes regulating stemness to promote the progression of HCC and detected the possible pathways through the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Results: Herein, we revealed that the high stemness of HCC correlated with an unfavorable prognosis. We obtained 31 up-regulated and 8 down-regulated DEPs in HCC and screened CTB-63M22.1, a poor prognostic indicator of HCC. In addition, CTB-63M22.1 had a mutation frequency similar to mRNAsi and acted in a module similar to that of mRNAsi on HCC. We then screened two RNA-binding proteins (RBPs) LIN28B and NOP56 with the highest correlation with stemness. We also discovered that they were primarily enriched in the biological process as examples of mitotic nuclear division and cell cycle. Conclusions: Collectively, these results revealed that pseudogenes CTB-63M22.1 may regulate cancer stemness by regulating RBPs, suggesting that CTB-63M22.1 may serve as an innovative therapeutic target and a reliable prognostic marker for HCC. [ABSTRACT FROM AUTHOR]

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]

Analysis of publicly available whole‐genome sequence data from the Human Pangenome Project and the 1000 Genomes Project has identified a DNA segment of approximately 60 kb in the major histocompatibility complex (MHC) between HLA‐W and HLA‐J that is present in some MHC haplotypes but not others. This DNA segment is largely repeat element‐rich but includes the pseudogene HLA‐Y, thus pinpointing the location of this pseudogene, and a new HLA class I sequence we have called HLA‐OLI. HLA‐OLI clusters phylogenetically with the HLA class I pseudogenes, HLA‐P and HLA‐W, and appears to have a similar genetic structure. The availability of whole‐genome sequence data from diverse populations enables a detailed characterization of the MHC at the population level and will have implications for understanding MHC disease associations and the non‐HLA MHC factors that impact unrelated hematopoietic cell transplant outcomes. [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]