SinEx DB: a database for single exon coding sequences in mammalian genomes

Mutational Bias and Natural Selection Driving the Synonymous Codon Usage of Single-Exon Genes in Rice (Oryza sativa L.)

The relative abundance of single-exon genes (SEGs) in higher plants is perplexing. Uncovering the synonymous codon usage pattern of SEGs will benefit for further understanding their underlying evolutionary mechanism in plants. Using internal correspondence analysis (ICA), we reveal a significant difference in synonymous codon usage between SEGs and multiple-exon genes (MEGs) in rice. But the effect is weak, accounting for only 2.61% of the total codon usage variability. SEGs and MEGs contain remarkably different base compositions, and are under clearly differential selective constraints, with the former having higher GC content, and evolving relatively faster during evolution. In the group of SEGs, the variability in synonymous codon usage among genes is partially due to the variations in GC content, gene function, and gene expression level, which accounts for 22.03%, 5.99%, and 3.32% of the total codon usage variability, respectively. Therefore, mutational bias and natural selection should work on affecting the synonymous codon usage of SEGs in rice. These findings may deepen our knowledge for the mechanisms of origination, differentiation and regulation of SEGs in plants.

Identification of Potential Long Non-Coding RNA Candidates that Contribute to Triple-Negative Breast Cancer in Humans through Computational Approach

Breast cancer (BC) is the most frequent malignancy identified in adult females, resulting in enormous financial losses worldwide. Owing to the heterogeneity as well as various molecular subtypes, the molecular pathways underlying carcinogenesis in various forms of BC are distinct. Therefore, the advancement of alternative therapy is required to combat the ailment. Recent analyses propose that long non-coding RNAs (lncRNAs) perform an essential function in controlling immune response, and therefore, may provide essential information about the disorder. However, their function in patients with triple-negative BC (TNBC) has not been explored in detail. Here, we analyzed the changes in the genomic expression of messenger RNA (mRNA) and lncRNA in standard control in response to cancer metastasis using publicly available single-cell RNA-Seq data. We identified a total of 197 potentially novel lncRNAs in TNBC patients of which 86 were differentially upregulated and 111 were differentially downregulated. In addition, among the 909 candidate lncRNA transcripts, 19 were significantly differentially expressed (DE) of which three were upregulated and 16 were downregulated. On the other hand, 1901 mRNA transcripts were significantly DE of which 1110 were upregulated and 791 were downregulated by TNBCs subtypes. The Gene Ontology (GO) analyses showed that some of the host genes were enriched in various biological, molecular, and cellular functions. The Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis showed that some of the genes were involved in only one pathway of prostate cancer. The lncRNA-miRNA-gene network analysis showed that the lncRNAs TCONS_00076394 and TCONS_00051377 interacted with breast cancer-related micro RNAs (miRNAs) and the host genes of these lncRNAs were also functionally related to breast cancer. Thus, this study provides novel lncRNAs as potential biomarkers for the therapeutic intervention of this cancer subtype.

Evolutionary Perspective and Expression Analysis of Intronless Genes Highlight the Conservation of Their Regulatory Role

The structure of eukaryotic genes is generally a combination of exons interrupted by intragenic non-coding DNA regions (introns) removed by RNA splicing to generate the mature mRNA. A fraction of genes, however, comprise a single coding exon with introns in their untranslated regions or are intronless genes (IGs), lacking introns entirely. The latter code for essential proteins involved in development, growth, and cell proliferation and their expression has been proposed to be highly specialized for neuro-specific functions and linked to cancer, neuropathies, and developmental disorders. The abundant presence of introns in eukaryotic genomes is pivotal for the precise control of gene expression. Notwithstanding, IGs exempting splicing events entail a higher transcriptional fidelity, making them even more valuable for regulatory roles. This work aimed to infer the functional role and evolutionary history of IGs centered on the mouse genome. IGs consist of a subgroup of genes with one exon including coding genes, non-coding genes, and pseudogenes, which conform approximately 6% of a total of 21,527 genes. To understand their prevalence, biological relevance, and evolution, we identified and studied 1,116 IG functional proteins validating their differential expression in transcriptomic data of embryonic mouse telencephalon. Our results showed that overall expression levels of IGs are lower than those of MEGs. However, strongly up-regulated IGs include transcription factors (TFs) such as the class 3 of POU (HMG Box), Neurog1, Olig1, and BHLHe22, BHLHe23, among other essential genes including the β-cluster of protocadherins. Most striking was the finding that IG-encoded BHLH TFs fit the criteria to be classified as microproteins. Finally, predicted protein orthologs in other six genomes confirmed high conservation of IGs associated with regulating neural processes and with chromatin organization and epigenetic regulation in Vertebrata. Moreover, this study highlights that IGs are essential modulators of regulatory processes, such as the Wnt signaling pathway and biological processes as pivotal as sensory organ developing at a transcriptional and post-translational level. Overall, our results suggest that IG proteins have specialized, prevalent, and unique biological roles and that functional divergence between IGs and MEGs is likely to be the result of specific evolutionary constraints.

Genome-wide analysis of European sea bass provides insights into the evolution and functions of single-exon genes

Several studies have attempted to understand the origin and evolution of single-exon genes (SEGs) in eukaryotic organisms, including fishes, but few have examined the functional and evolutionary relationships between SEGs and multiple-exon gene (MEG) paralogs, in particular the conservation of promoter regions. Given that SEGs originate via the reverse transcription of mRNA from a "parental" MEGs, such comparisons may enable identifying evolutionarily-related SEG/MEG paralogs, which might fulfill equivalent physiological functions. Here, the relationship of SEG proportion with MEG count, gene density, intron count, and chromosome size was assessed for the genome of the European sea bass, Dicentrarchus labrax. Then, SEGs with an MEG parent were identified, and promoter sequences of SEG/MEG paralogs were compared, to identify highly conserved functional motifs. The results revealed a total count of 1,585 (8.3% of total genes) SEGs in the European sea bass genome, which was correlated with MEG count but not with gene density. The significant correlation of SEG content with the number of MEGs suggests that SEGs were continuously and independently generated over evolutionary time following species divergence through retrotranscription events, followed by tandem duplications. Functional annotation showed that the majority of SEGs are functional, as is evident from their expression in RNA-seq data used to support homology-based genome annotation. Differences in 5'UTR and 3'UTR lengths between SEG/MEG paralogs observed in this study may contribute to gene expression divergence between them and therefore lead to the emergence of new SEG functions. The comparison of nonsynonymous to synonymous changes (Ka/Ks) between SEG/MEG parents showed that 74 of them are under positive selection (Ka/Ks > 1; p = .0447). An additional fifteen SEGs with an MEG parent have a common promoter, which implies that they are under the influence of common regulatory networks.

SinEx DB 2.0 update 2020: database for eukaryotic single-exon coding sequences

Single-exon coding sequences (CDSs), also known as ‘single-exon genes’ (SEGs), are defined as nuclear, protein-coding genes that lack introns in their CDSs. They have been studied not only to determine their origin and evolution but also because their expression has been linked to several types of human cancers and neurological/developmental disorders, and many exhibit tissue-specific transcription. We developed SinEx DB that houses DNA and protein sequence information of SEGs from 10 mammalian genomes including human. SinEx DB includes their functional predictions (KOG (euKaryotic Orthologous Groups)) and the relative distribution of these functions within species. Here, we report SinEx 2.0, a major update of SinEx DB that includes information of the occurrence, distribution and functional prediction of SEGs from 60 completely sequenced eukaryotic genomes, representing animals, fungi, protists and plants. The information is stored in a relational database built with MySQL Server 5.7, and the complete dataset of SEG sequences and their GO (Gene Ontology) functional assignations are available for downloading. SinEx DB 2.0 was built with a novel pipeline that helps disambiguate single-exon isoforms from SEGs. SinEx DB 2.0 is the largest available database for SEGs and provides a rich source of information for advancing our understanding of the evolution, function of SEGs and their associations with disorders including cancers and neurological and developmental diseases.

Database URL:http://v2.sinex.cl/

The evolution of the huntingtin-associated protein 40 (HAP40) in conjunction with huntingtin

BACKGROUND

The huntingtin-associated protein 40 (HAP40) abundantly interacts with huntingtin (HTT), the protein that is altered in Huntington's disease (HD). Therefore, we analysed the evolution of HAP40 and its interaction with HTT.

RESULTS

We found that in amniotes HAP40 is encoded by a single-exon gene, whereas in all other organisms it is expressed from multi-exon genes. HAP40 co-occurs with HTT in unikonts, including filastereans such as Capsaspora owczarzaki and the amoebozoan Dictyostelium discoideum, but both proteins are absent from fungi. Outside unikonts, a few species, such as the free-living amoeboflagellate Naegleria gruberi, contain putative HTT and HAP40 orthologs. Biochemically we show that the interaction between HTT and HAP40 extends to fish, and bioinformatic analyses provide evidence for evolutionary conservation of this interaction. The closest homologue of HAP40 in current protein databases is the family of soluble N-ethylmaleimide-sensitive factor attachment proteins (SNAPs).

CONCLUSION

Our results indicate that the transition from a multi-exon to a single-exon gene appears to have taken place by retroposition during the divergence of amphibians and amniotes, followed by the loss of the parental multi-exon gene. Furthermore, it appears that the two proteins probably originated at the root of eukaryotes. Conservation of the interaction between HAP40 and HTT and their likely coevolution strongly indicate functional importance of this interaction.

[Virus-host coevolution: Endogenous RNA viral elements as pseudogenes]

RNA viruses do not need to take the form of DNAs, and RNAs alone complete their replication cycles. On the other hand, since the 1970s, it has been known that DNA fragments derived from RNA viruses can be detected in RNA virus-infected cells. Furthermore, in this decade, it has become clear that the eukaryotic genomes contain genetic sequences derived from non-retroviral RNA viruses. The DNA sequences derived from these RNA viruses are thought to be generatedby using a transposable mechanism of retrotransposon, such as LINE-1. Many endogenous RNA viral sequences are formed by the same mechanism as processed pseudogenes in eukaryotic cells, but the significance of the production of RNA viral "pseudogenes " in infected cells has not been elucidated. We have discovered endogenous bornavirus-like elements (EBLs), which derived from a negative-sense, single-stranded RNA virus, Bornaviruses, and have studied the evolution and function of EBLs in host animals. The analysis of EBLs provides us a clue to unravel the history of host-RNA virus coexistence. In this review, I overview about the function of endogenous RNA virus sequences, especially EBLs in mammalian genomes, and discuss the significance of endogenization of RNA viruses as viral pseudogenes in evolution.

Genetic basis of functional variability in adhesion G protein-coupled receptors

The enormous sizes of adhesion G protein-coupled receptors (aGPCRs) go along with complex genomic exon-intron architectures giving rise to multiple mRNA variants. There is a need for a comprehensive catalog of aGPCR variants for proper evaluation of the complex functions of aGPCRs found in structural, in vitro and animal model studies. We used an established bioinformatics pipeline to extract, quantify and visualize mRNA variants of aGPCRs from deeply sequenced transcriptomes. Data analysis showed that aGPCRs have multiple transcription start sites even within introns and that tissue-specific splicing is frequent. On average, 19 significantly expressed transcript variants are derived from a given aGPCR gene. The domain architecture of the N terminus encoded by transcript variants often differs and N termini without or with an incomplete seven-helix transmembrane anchor as well as separate seven-helix transmembrane domains are frequently derived from aGPCR genes. Experimental analyses of selected aGPCR transcript variants revealed marked functional differences. Our analysis has an impact on a rational design of aGPCR constructs for structural analyses and gene-deficient mouse lines and provides new support for independent functions of both, the large N terminus and the transmembrane domain of aGPCRs.

Small, but surprisingly repetitive genomes: transposon expansion and not polyploidy has driven a doubling in genome size in a metazoan species complex

BACKGROUND

The causes and consequences of genome size variation across Eukaryotes, which spans five orders of magnitude, have been hotly debated since before the advent of genome sequencing. Previous studies have mostly examined variation among larger taxonomic units (e.g., orders, or genera), while comparisons among closely related species are rare. Rotifers of the Brachionus plicatilis species complex exhibit a seven-fold variation in genome size and thus represent a unique opportunity to study such changes on a relatively short evolutionary timescale. Here, we sequenced and analysed the genomes of four species of this complex with nuclear DNA contents spanning 110-422 Mbp. To establish the likely mechanisms of genome size change, we analysed both sequencing read libraries and assemblies for signatures of polyploidy and repetitive element content. We also compared these genomes to that of B. calyciflorus, the closest relative with a sequenced genome (293 Mbp nuclear DNA content).

RESULTS

Despite the very large differences in genome size, we saw no evidence of ploidy level changes across the B. plicatilis complex. However, repetitive element content explained a large portion of genome size variation (at least 54%). The species with the largest genome, B. asplanchnoidis, has a strikingly high 44% repetitive element content, while the smaller B. plicatilis genomes contain between 14 and 25% repetitive elements. According to our analyses, the B. calyciflorus genome contains 39% repetitive elements, which is substantially higher than previously reported (21%), and suggests that high repetitive element load could be widespread in monogonont rotifers.

CONCLUSIONS

Even though the genome sizes of these species are at the low end of the metazoan spectrum, their genomes contain substantial amounts of repetitive elements. Polyploidy does not appear to play a role in genome size variations in these species, and these variations can be mostly explained by changes in repetitive element content. This contradicts the naïve expectation that small genomes are streamlined, or less complex, and that large variations in nuclear DNA content between closely related species are due to polyploidy.

Analysis of long non-coding RNA and mRNA expression in bovine macrophages brings up novel aspects of Mycobacterium avium subspecies paratuberculosis infections

Paratuberculosis is a major disease in cattle that severely affects animal welfare and causes huge economic losses worldwide. Development of alternative diagnostic methods is of urgent need to control the disease. Recent studies suggest that long non-coding RNAs (lncRNAs) play a crucial role in regulating immune function and may confer valuable information about the disease. However, their role has not yet been investigated in cattle with respect to infection towards Paratuberculosis. Therefore, we investigated the alteration in genomic expression profiles of mRNA and lncRNA in bovine macrophages in response to Paratuberculosis infection using RNA-Seq. We identified 397 potentially novel lncRNA candidates in macrophages of which 38 were differentially regulated by the infection. A total of 820 coding genes were also significantly altered by the infection. Co-expression analysis of lncRNAs and their neighbouring coding genes suggest regulatory functions of lncRNAs in pathways related to immune response. For example, this included protein coding genes such as TNIP3, TNFAIP3 and NF-κB2 that play a role in NF-κB2 signalling, a pathway associated with immune response. This study advances our understanding of lncRNA roles during Paratuberculosis infection.

Evolutionary, structural and functional analysis of the caleosin/peroxygenase gene family in the Fungi

BACKGROUND

Caleosin/peroxygenases, CLO/PXG, (designated PF05042 in Pfam) are a group of genes/proteins with anomalous distributions in eukaryotic taxa. We have previously characterised CLO/PXGs in the Viridiplantae. The aim of this study was to investigate the evolution and functions of the CLO/PXGs in the Fungi and other non-plant clades and to elucidate the overall origin of this gene family.

RESULTS

CLO/PXG-like genes are distributed across the full range of fungal groups from the basal clades, Cryptomycota and Microsporidia, to the largest and most complex Dikarya species. However, the genes were only present in 243 out of 844 analysed fungal genomes. CLO/PXG-like genes have been retained in many pathogenic or parasitic fungi that have undergone considerable genomic and structural simplification, indicating that they have important functions in these species. Structural and functional analyses demonstrate that CLO/PXGs are multifunctional proteins closely related to similar proteins found in all major taxa of the Chlorophyte Division of the Viridiplantae. Transcriptome and physiological data show that fungal CLO/PXG-like genes have complex patterns of developmental and tissue-specific expression and are upregulated in response to a range of biotic and abiotic stresses as well as participating in key metabolic and developmental processes such as lipid metabolism, signalling, reproduction and pathogenesis. Biochemical data also reveal that the Aspergillus flavus CLO/PXG has specific functions in sporulation and aflatoxin production as well as playing roles in lipid droplet function.

CONCLUSIONS

In contrast to plants, CLO/PXGs only occur in about 30% of sequenced fungal genomes but are present in all major taxa. Fungal CLO/PXGs have similar but not identical roles to those in plants, including stress-related oxylipin signalling, lipid metabolism, reproduction and pathogenesis. While the presence of CLO/PXG orthologs in all plant genomes sequenced to date would suggest that they have core housekeeping functions in plants, the selective loss of CLO/PXGs in many fungal genomes suggests more restricted functions in fungi as accessory genes useful in particular environments or niches. We suggest an ancient origin of CLO/PXG-like genes in the 'last eukaryotic common ancestor' (LECA) and their subsequent loss in ancestors of the Metazoa, after the latter had diverged from the ancestral fungal lineage.

The Reprimo Gene Family: A Novel Gene Lineage in Gastric Cancer with Tumor Suppressive Properties

The reprimo (RPRM) gene family is a group of single exon genes present exclusively within the vertebrate lineage. Two out of three members of this family are present in humans: RPRM and RPRM-Like (RPRML). RPRM induces cell cycle arrest at G2/M in response to p53 expression. Loss-of-expression of RPRM is related to increased cell proliferation and growth in gastric cancer. This evidence suggests that RPRM has tumor suppressive properties. However, the molecular mechanisms and signaling partners by which RPRM exerts its functions remain unknown. Moreover, scarce studies have attempted to characterize RPRML, and its functionality is unclear. Herein, we highlight the role of the RPRM gene family in gastric carcinogenesis, as well as its potential applications in clinical settings. In addition, we summarize the current knowledge on the phylogeny and expression patterns of this family of genes in embryonic zebrafish and adult humans. Strikingly, in both species, RPRM is expressed primarily in the digestive tract, blood vessels and central nervous system, supporting the use of zebrafish for further functional characterization of RPRM. Finally, drawing on embryonic and adult expression patterns, we address the potential relevance of RPRM and RPRML in cancer. Active investigation or analytical research in the coming years should contribute to novel translational applications of this poorly understood gene family as potential biomarkers and development of novel cancer therapies.

Phylogenetic and genomic analyses of the ribosomal oxygenases Riox1 (No66) and Riox2 (Mina53) provide new insights into their evolution

BACKGROUND

Translation of specific mRNAs can be highly regulated in different cells, tissues or under pathological conditions. Ribosome heterogeneity can originate from variable expression or post-translational modifications of ribosomal proteins. The ribosomal oxygenases RIOX1 (NO66) and RIOX2 (MINA53) modify ribosomal proteins by histidine hydroxylation. A similar mechanism is present in prokaryotes. Thus, ribosome hydroxylation may be a well-conserved regulatory mechanism with implications in disease and development. However, little is known about the evolutionary history of Riox1 and Riox2 genes and their encoded proteins across eukaryotic taxa.

RESULTS

In this study, we have analysed Riox1 and Riox2 orthologous genes from 49 metazoen species and have constructed phylogenomic trees for both genes. Our genomic and phylogenetic analyses revealed that Arthropoda, Annelida, Nematoda and Mollusca lack the Riox2 gene, although in the early phylum Cnidaria both genes, Riox1 and Riox2, are present and expressed. Riox1 is an intronless single-exon-gene in several species, including humans. In contrast to Riox2, Riox1 is ubiquitously present throughout the animal kingdom suggesting that Riox1 is the phylogenetically older gene from which Riox2 has evolved. Both proteins have maintained a unique protein architecture with conservation of active sites within the JmjC domains, a dimerization domain, and a winged-helix domain. In addition, Riox1 proteins possess a unique N-terminal extension domain. Immunofluorescence analyses in Hela cells and in Hydra vulgaris identified a nucleolar localisation signal within the extended N-terminal domain of human RIOX1 and an altered subnuclear localisation for the Hydra Riox2.

CONCLUSIONS

Conserved active site residues and uniform protein domain architecture suggest a consistent enzymatic activity within the Riox orthologs throughout evolution. However, differences in genomic architecture, like single exon genes and alterations in subnuclear localisation, as described for Hydra, point towards adaption mechanisms that may correlate with taxa- or species-specific requirements. The diversification of Riox1/Riox2 gene structures throughout evolution suggest that functional requirements in expression of protein isoforms and/or subcellular localisation of proteins may have evolved by adaptation to lifestyle.

Improved ontology for eukaryotic single-exon coding sequences in biological databases

Efficient extraction of knowledge from biological data requires the development of structured vocabularies to unambiguously define biological terms. This paper proposes descriptions and definitions to disambiguate the term 'single-exon gene'. Eukaryotic Single-Exon Genes (SEGs) have been defined as genes that do not have introns in their protein coding sequences. They have been studied not only to determine their origin and evolution but also because their expression has been linked to several types of human cancer and neurological/developmental disorders and many exhibit tissue-specific transcription. Unfortunately, the term 'SEGs' is rife with ambiguity, leading to biological misinterpretations. In the classic definition, no distinction is made between SEGs that harbor introns in their untranslated regions (UTRs) versus those without. This distinction is important to make because the presence of introns in UTRs affects transcriptional regulation and post-transcriptional processing of the mRNA. In addition, recent whole-transcriptome shotgun sequencing has led to the discovery of many examples of single-exon mRNAs that arise from alternative splicing of multi-exon genes, these single-exon isoforms are being confused with SEGs despite their clearly different origin. The increasing expansion of RNA-seq datasets makes it imperative to distinguish the different SEG types before annotation errors become indelibly propagated in biological databases. This paper develops a structured vocabulary for their disambiguation, allowing a major reassessment of their evolutionary trajectories, regulation, RNA processing and transport, and provides the opportunity to improve the detection of gene associations with disorders including cancers, neurological and developmental diseases.

Protein-Coding Genes' Retrocopies and Their Functions

Transposable elements, often considered to be not important for survival, significantly contribute to the evolution of transcriptomes, promoters, and proteomes. Reverse transcriptase, encoded by some transposable elements, can be used in trans to produce a DNA copy of any RNA molecule in the cell. The retrotransposition of protein-coding genes requires the presence of reverse transcriptase, which could be delivered by either non-long terminal repeat (non-LTR) or LTR transposons. The majority of these copies are in a state of “relaxed” selection and remain “dormant” because they are lacking regulatory regions; however, many become functional. In the course of evolution, they may undergo subfunctionalization, neofunctionalization, or replace their progenitors. Functional retrocopies (retrogenes) can encode proteins, novel or similar to those encoded by their progenitors, can be used as alternative exons or create chimeric transcripts, and can also be involved in transcriptional interference and participate in the epigenetic regulation of parental gene expression. They can also act in trans as natural antisense transcripts, microRNA (miRNA) sponges, or a source of various small RNAs. Moreover, many retrocopies of protein-coding genes are linked to human diseases, especially various types of cancer.