Evolutionary dynamics of coding and non-coding transcriptomes

Oligodendrocyte-specific expression of PSG8-AS1 suggests a role in myelination with prognostic value in oligodendroglioma

The segmentally duplicated Pregnancy-specific glycoprotein (PSG) locus on chromosome 19q13 may be one of the most rapidly evolving in the human genome. It comprises ten coding genes (PSG1-9, 11) and one predominantly non-coding gene (PSG10) that are expressed in the placenta and gut, in addition to several poorly characterized long non-coding RNAs. We report that long non-coding RNA PSG8-AS1 has an oligodendrocyte-specific expression pattern and is co-expressed with genes encoding key myelin constituents. PSG8-AS1 exhibits two peaks of expression during human brain development coinciding with the most active periods of oligodendrogenesis and myelination. PSG8-AS1 orthologs were found in the genomes of several primates but significant expression was found only in the human, suggesting a recent evolutionary origin of its proposed role in myelination. Additionally, because co-deletion of chromosomes 1p/19q is a genomic marker of oligodendroglioma, expression of PSG8-AS1 was examined in these tumors. PSG8-AS1 may be a promising diagnostic biomarker for glioma, with prognostic value in oligodendroglioma.

Targeting and engineering long non-coding RNAs for cancer therapy

RNA therapeutics (RNATx) aim to treat diseases, including cancer, by targeting or employing RNA molecules for therapeutic purposes. Amongst the most promising targets are long non-coding RNAs (lncRNAs), which regulate oncogenic molecular networks in a cell type-restricted manner. lncRNAs are distinct from protein-coding genes in important ways that increase their therapeutic potential yet also present hurdles to conventional clinical development. Advances in genome editing, oligonucleotide chemistry, multi-omics and RNA engineering are paving the way for efficient and cost-effective lncRNA-focused drug discovery pipelines. In this Review, we present the emerging field of lncRNA therapeutics for oncology, with emphasis on the unique strengths and challenges of lncRNAs within the broader RNATx framework. We outline the necessary steps for lncRNA therapeutics to deliver effective, durable, tolerable and personalized treatments for cancer.

FEVER: an interactive web-based resource for evolutionary transcriptomics across fishes

Teleost fish represent one of the largest and most diverse clades of vertebrates, which makes them great models in various research areas such as ecology and evolution. Recent sequencing endeavors provided high-quality genomes for species covering the main fish evolutionary lineages, opening up large-scale comparative genomics studies. However, transcriptomic data across fish species and organs are heterogenous and have not been integrated with newly sequenced genomes making gene expression quantification and comparative analyses particularly challenging. Thus, resources integrating genomic and transcriptomic data across fish species and organs are still lacking. Here, we present FEVER, a web-based resource allowing evolutionary transcriptomics across species and tissues. First, based on query genes FEVER reconstructs gene trees providing orthologous and paralogous relationships as well as their evolutionary dynamics across 13 species covering the major fish lineages, and 4 model species as evolutionary outgroups. Second, it provides unbiased gene expression across 11 tissues using up-to-date fish genomes. Finally, genomic and transcriptomic data are combined together allowing the exploration of gene expression evolution following speciation and duplication events. FEVER is freely accessible at https://fever.sk8.inrae.fr/.

Emerging roles of long non-coding RNAs in human epilepsy

Genome-scale biological studies conducted in the post-genomic era have revealed that two-thirds of human genes do not encode proteins. Most functional non-coding RNA transcripts in humans are products of long non-coding RNA (lncRNA) genes, an abundant but still poorly understood class of human genes. As a result of their fundamental and multitasking regulatory roles, lncRNAs are associated with a wide range of human diseases, including neurological disorders. Approximately 40% of lncRNAs are specifically expressed in the brain, and many of them exhibit distinct spatiotemporal patterns of expression. Comparative genomics approaches have determined that 65%-75% of human lncRNA genes are primate-specific and hence can be posited as a contributing potential cause of the higher-order complexity of primates, including human, brains relative to those of other mammals. Although lncRNAs present important mechanistic examples of epileptogenic functions, the human/primate specificity of lncRNAs questions their relevance in rodent models. Here, we present an in-depth review that supports the contention that human lncRNAs are direct contributors to the etiology and pathogenesis of human epilepsy, as a means to accelerate the integration of lncRNAs into clinical practice as potential diagnostic biomarkers and therapeutic targets. Meta-analytically, the major finding of our review is the commonality of lncRNAs in epilepsy and cancer pathogenesis through mitogen-activated protein kinase (MAPK)-related pathways. In addition, neuroinflammation may be a relevant part of the common pathophysiology of cancer and epilepsy. LncRNAs affect neuroinflammation-related signaling pathways such as nuclear factor kappa- light- chain- enhancer of activated B cells (NF-κB), Notch, and phosphatidylinositol 3- kinase/ protein kinase B (Akt) (PI3K/AKT), with the NF-κB pathway being the most common. Besides the controversy over lncRNA research in non-primate models, whether neuroinflammation is triggered by injury and/or central nervous system (CNS) toxicity during epilepsy modeling in animals or is a direct consequence of epilepsy pathophysiology needs to be considered meticulously in future studies.

Host adaptive radiation is associated with rapid virus diversification and cross-species transmission in African cichlid fishes

Adaptive radiations are generated through a complex interplay of biotic and abiotic factors. Although adaptive radiations have been widely studied in the context of animal and plant evolution, little is known about how they impact the evolution of the viruses that infect these hosts, which in turn may provide insights into the drivers of cross-species transmission and hence disease emergence. We examined how the rapid adaptive radiation of the cichlid fishes of African Lake Tanganyika over the last 10 million years has shaped the diversity and evolution of the viruses they carry. Through metatranscriptomic analysis of 2,242 RNA sequencing libraries, we identified 121 vertebrate-associated viruses among various tissue types that fell into 13 RNA and 4 DNA virus groups. Host-switching was commonplace, particularly within the Astroviridae, Metahepadnavirus, Nackednavirus, Picornaviridae, and Hepacivirus groups, occurring more frequently than in other fish communities. A time-calibrated phylogeny revealed that hepacivirus diversification was not constant throughout the cichlid radiation but accelerated 2-3 million years ago, coinciding with a period of rapid cichlid diversification and niche packing in Lake Tanganyika, thereby providing more closely related hosts for viral infection. These data depict a dynamic virus ecosystem within the cichlids of Lake Tanganyika, characterized by rapid virus diversification and frequent host jumping, and likely reflecting their close phylogenetic relationships that lower the barriers to cross-species virus transmission.

Identification of constrained sequence elements across 239 primate genomes

Noncoding DNA is central to our understanding of human gene regulation and complex diseases1,2, and measuring the evolutionary sequence constraint can establish the functional relevance of putative regulatory elements in the human genome3-9. Identifying the genomic elements that have become constrained specifically in primates has been hampered by the faster evolution of noncoding DNA compared to protein-coding DNA10, the relatively short timescales separating primate species11, and the previously limited availability of whole-genome sequences12. Here we construct a whole-genome alignment of 239 species, representing nearly half of all extant species in the primate order. Using this resource, we identified human regulatory elements that are under selective constraint across primates and other mammals at a 5% false discovery rate. We detected 111,318 DNase I hypersensitivity sites and 267,410 transcription factor binding sites that are constrained specifically in primates but not across other placental mammals and validate their cis-regulatory effects on gene expression. These regulatory elements are enriched for human genetic variants that affect gene expression and complex traits and diseases. Our results highlight the important role of recent evolution in regulatory sequence elements differentiating primates, including humans, from other placental mammals.

Population-level annotation of lncRNAs in Arabidopsis reveals extensive expression variation associated with transposable element-like silencing

Long noncoding RNAs (lncRNAs) are understudied and underannotated in plants. In mammals, lncRNA loci are nearly as ubiquitous as protein-coding genes, and their expression is highly variable between individuals of the same species. Using Arabidopsis thaliana as a model, we aimed to elucidate the true scope of lncRNA transcription across plants from different regions and study its natural variation. We used transcriptome deep sequencing data sets spanning hundreds of natural accessions and several developmental stages to create a population-wide annotation of lncRNAs, revealing thousands of previously unannotated lncRNA loci. While lncRNA transcription is ubiquitous in the genome, most loci appear to be actively silenced and their expression is extremely variable between natural accessions. This high expression variability is largely caused by the high variability of repressive chromatin levels at lncRNA loci. High variability was particularly common for intergenic lncRNAs (lincRNAs), where pieces of transposable elements (TEs) present in 50% of these lincRNA loci are associated with increased silencing and variation, and such lncRNAs tend to be targeted by the TE silencing machinery. We created a population-wide lncRNA annotation in Arabidopsis and improve our understanding of plant lncRNA genome biology, raising fundamental questions about what causes transcription and silencing across the genome.

Evolutionary analysis of conserved non-coding elements subsequent to whole-genome duplication in opium poppy

Whole-genome duplication (WGD) leads to the duplication of both coding and non-coding sequences within an organism's genome, providing an abundant supply of genetic material that can drive evolution, ultimately contributing to plant adaptation and speciation. Although non-coding sequences contain numerous regulatory elements, they have been understudied compared to coding sequences. In order to address this gap, we explored the evolutionary patterns of regulatory sequences, coding sequences and transcriptomes using conserved non-coding elements (CNEs) as regulatory element proxies following the recent WGD event in opium poppy (Papaver somniferum). Our results showed similar evolutionary patterns in subgenomes of regulatory and coding sequences. Specifically, the biased or unbiased retention of coding sequences reflected the same pattern as retention levels in regulatory sequences. Further, the divergence of gene expression patterns mediated by regulatory element variations occurred at a more rapid pace than that of gene coding sequences. However, gene losses were purportedly dependent on relaxed selection pressure in coding sequences. Specifically, the rapid evolution of tissue-specific benzylisoquinoline alkaloid production in P. somniferum was associated with regulatory element changes. The origin of a novel stem-specific ACR, which utilized ancestral cis-elements as templates, is likely to be linked to the evolutionary trajectory behind the transition of the PSMT1-CYP719A21 cluster from high levels of expression solely in P. rhoeas root tissue to its elevated expression in P. somniferum stem tissue. Our findings demonstrate that rapid regulatory element evolution can contribute to the emergence of new phenotypes and provide valuable insights into the high evolvability of regulatory elements.

Exploring the genetic landscape of nitrogen uptake in durum wheat: genome-wide characterization and expression profiling of NPF and NRT2 gene families

Nitrate uptake by plants primarily relies on two gene families: Nitrate transporter 1/peptide transporter (NPF) and Nitrate transporter 2 (NRT2). Here, we extensively characterized the NPF and NRT2 families in the durum wheat genome, revealing 211 NPF and 20 NRT2 genes. The two families share many Cis Regulatory Elements (CREs) and Transcription Factor binding sites, highlighting a partially overlapping regulatory system and suggesting a coordinated response for nitrate transport and utilization. Analyzing RNA-seq data from 9 tissues and 20 cultivars, we explored expression profiles and co-expression relationships of both gene families. We observed a strong correlation between nucleotide variation and gene expression within the NRT2 gene family, implicating a shared selection mechanism operating on both coding and regulatory regions. Furthermore, NPF genes showed highly tissue-specific expression profiles, while NRT2s were mainly divided in two co-expression modules, one expressed in roots (NAR2/NRT3 dependent) and the other induced in anthers and/ovaries during maturation. Our evidences confirmed that the majority of these genes were retained after small-scale duplication events, suggesting a neo- or sub-functionalization of many NPFs and NRT2s. Altogether, these findings indicate that the expansion of these gene families in durum wheat could provide valuable genetic variability useful to identify NUE-related and candidate genes for future breeding programs in the context of low-impact and sustainable agriculture.

The potential role of miRNA in regulating macrophage polarization

Macrophage polarization is a dynamic process determining the outcome of various physiological and pathological situations through inducing pro-inflammatory responses or resolving inflammation via exerting anti-inflammatory effects. The miRNAs are epigenetic regulators of different biologic pathways that target transcription factors and signaling molecules to promote macrophage phenotype transition and regulate immune responses. Modulating the macrophage activation, differentiation, and polarization by miRNAs is crucial for immune responses in response to microenvironmental signals and under various physiological and pathological conditions. In term of clinical significance, regulating macrophage polarization via miRNAs could be utilized for inflammation control. Also, understanding the role of miRNAs in macrophage polarization can provide insights into diagnostic strategies associated with dysregulated miRNAs and for developing macrophage-centered therapeutic methods. In this case, targeting miRNAs to further regulate of macrophage polarization may become an efficient strategy for treating immune-associated disorders. The current review investigated and categorized various miRNAs directly or indirectly involved in macrophage polarization by targeting different transcription factors and signaling pathways. In addition, prospects for regulating macrophage polarization via miRNA as a therapeutic choice that could be implicated in various pathological conditions, including cancer or inflammation-mediated injuries, were discussed.

Evolutionary Analyses of Gene Expression Divergence in Panicum hallii: Exploring Constitutive and Plastic Responses Using Reciprocal Transplants

The evolution of gene expression is thought to be an important mechanism of local adaptation and ecological speciation. Gene expression divergence occurs through the evolution of cis- polymorphisms and through more widespread effects driven by trans-regulatory factors. Here, we explore expression and sequence divergence in a large sample of Panicum hallii accessions encompassing the species range using a reciprocal transplantation experiment. We observed widespread genotype and transplant site drivers of expression divergence, with a limited number of genes exhibiting genotype-by-site interactions. We used a modified FST-QST outlier approach (QPC analysis) to detect local adaptation. We identified 514 genes with constitutive expression divergence above and beyond the levels expected under neutral processes. However, no plastic expression responses met our multiple testing correction as QPC outliers. Constitutive QPC outlier genes were involved in a number of developmental processes and responses to abiotic environments. Leveraging earlier expression quantitative trait loci results, we found a strong enrichment of expression divergence, including for QPC outliers, in genes previously identified with cis and cis-environment interactions but found no patterns related to trans-factors. Population genetic analyses detected elevated sequence divergence of promoters and coding sequence of constitutive expression outliers but little evidence for positive selection on these proteins. Our results are consistent with a hypothesis of cis-regulatory divergence as a primary driver of expression divergence in P. hallii.

GC-biased gene conversion drives accelerated evolution of ultraconserved elements in mammalian and avian genomes

Ultraconserved elements (UCEs) are the most conserved regions among the genomes of evolutionarily distant species and are thought to play critical biological functions. However, some UCEs rapidly evolved in specific lineages, and whether they contributed to adaptive evolution is still controversial. Here, using an increased number of sequenced genomes with high taxonomic coverage, we identified 2191 mammalian UCEs and 5938 avian UCEs from 95 mammal and 94 bird genomes, respectively. Our results show that these UCEs are functionally constrained and that their adjacent genes are prone to widespread expression with low expression diversity across tissues. Functional enrichment of mammalian and avian UCEs shows different trends indicating that UCEs may contribute to adaptive evolution of taxa. Focusing on lineage-specific accelerated evolution, we discover that the proportion of fast-evolving UCEs in nine mammalian and 10 avian test lineages range from 0.19% to 13.2%. Notably, up to 62.1% of fast-evolving UCEs in test lineages are much more likely to result from GC-biased gene conversion (gBGC). A single cervid-specific gBGC region embracing the uc.359 allele significantly alters the expression of Nova1 and other neural-related genes in the rat brain. Combined with the altered regulatory activity of ancient gBGC-induced fast-evolving UCEs in eutherians, our results provide evidence that synergy between gBGC and selection shaped lineage-specific substitution patterns, even in the most constrained regulatory elements. In summary, our results show that gBGC played an important role in facilitating lineage-specific accelerated evolution of UCEs, and further support the idea that a combination of multiple evolutionary forces shapes adaptive evolution.

De novo assembly and annotation of the singing mouse genome

BACKGROUND

Developing genomic resources for a diverse range of species is an important step towards understanding the mechanisms underlying complex traits. Specifically, organisms that exhibit unique and accessible phenotypes-of-interest allow researchers to address questions that may be ill-suited to traditional model organisms. We sequenced the genome and transcriptome of Alston's singing mouse (Scotinomys teguina), an emerging model for social cognition and vocal communication. In addition to producing advertisement songs used for mate attraction and male-male competition, these rodents are diurnal, live at high-altitudes, and are obligate insectivores, providing opportunities to explore diverse physiological, ecological, and evolutionary questions.

RESULTS

Using PromethION, Illumina, and PacBio sequencing, we produced an annotated genome and transcriptome, which were validated using gene expression and functional enrichment analyses. To assess the usefulness of our assemblies, we performed single nuclei sequencing on cells of the orofacial motor cortex, a brain region implicated in song coordination, identifying 12 cell types.

CONCLUSIONS

These resources will provide the opportunity to identify the molecular basis of complex traits in singing mice as well as to contribute data that can be used for large-scale comparative analyses.

Visual opsin gene expression evolution in the adaptive radiation of cichlid fishes of Lake Tanganyika

Tuning the visual sensory system to the ambient light is essential for survival in many animal species. This is often achieved through duplication, functional diversification, and/or differential expression of visual opsin genes. Here, we examined 753 new retinal transcriptomes from 112 species of cichlid fishes from Lake Tanganyika to unravel adaptive changes in gene expression at the macro-evolutionary and ecosystem level of one of the largest vertebrate adaptive radiations. We found that, across the radiation, all seven cone opsins-but not the rhodopsin-rank among the most differentially expressed genes in the retina, together with other vision-, circadian rhythm-, and hemoglobin-related genes. We propose two visual palettes characteristic of very shallow- and deep-water living species, respectively, and show that visual system adaptations along two major ecological axes, macro-habitat and diet, occur primarily via gene expression variation in a subset of cone opsin genes.

The molecular evolution of mammalian spermatogenesis

The testis is a key male reproductive organ that produces gametes through the process of spermatogenesis. Testis morphologies, sperm phenotypes, and the process of spermatogenesis evolve rapidly in mammals, presumably due to the evolutionary pressure on males to give rise to their own offspring. Here, we review studies illuminating the molecular evolution of the testis, in particular large-scale transcriptomic studies, which were based on bulk tissue samples and, more recently, individual cells. Together with various genomic and epigenomic data, these studies have unveiled the cellular source, molecular mechanisms, and evolutionary forces that underlie the rapid phenotypic evolution of the testis. They also revealed shared (ancestral) and species-specific spermatogenic gene expression programs. The insights and available data that have accumulated also provide a valuable resource for the investigation and treatment of male fertility disorders - a dramatically increasing problem in modern industrial societies.

LncRNA Functional Screening in Organismal Development

Controversy continues over the functional prevalence of long non-coding RNAs (lncRNAs) despite their being widely investigated in all kinds of cells and organisms. In animals, lncRNAs have aroused general interest from exponentially increasing transcriptomic repertoires reporting their highly tissue-specific and developmentally dynamic expression, and more importantly, from growing experimental evidence supporting their functionality in facilitating organogenesis and individual fitness. In mammalian testes, while a great multitude of lncRNA species are identified, only a minority of them have been shown to be useful, and even fewer have been demonstrated as true requirements for male fertility using knockout models to date. This noticeable gap is attributed to the virtual existence of a large number of junk lncRNAs, the lack of an ideal germline culture system, difficulty in loss-of-function interrogation, and limited screening strategies. Facing these challenges, in this review, we discuss lncRNA functionality in organismal development and especially in mouse testis, with a focus on lncRNAs with functional screening.

Roles of lncRNAs in brain development and pathogenesis: Emerging therapeutic opportunities

The human genome is pervasively transcribed, producing a majority of short and long noncoding RNAs (lncRNAs) that can influence cellular programs through a variety of transcriptional and post-transcriptional regulatory mechanisms. The brain houses the richest repertoire of long noncoding transcripts, which function at every stage during central nervous system development and homeostasis. An example of functionally relevant lncRNAs is species involved in spatiotemporal organization of gene expression in different brain regions, which play roles at the nuclear level and in transport, translation, and decay of other transcripts in specific neuronal sites. Research in the field has enabled identification of the contributions of specific lncRNAs to certain brain diseases, including Alzheimer's disease, Parkinson's disease, cancer, and neurodevelopmental disorders, resulting in notions of potential therapeutic strategies that target these RNAs to recover the normal phenotype. Here, we summarize the latest mechanistic findings associated with lncRNAs in the brain, focusing on their dysregulation in neurodevelopmental or neurodegenerative disorders, their use as biomarkers for central nervous system (CNS) diseases in vitro and in vivo, and their potential utility for therapeutic strategies.

Pervasive tandem duplications and convergent evolution shape coral genomes

BACKGROUND

Over the last decade, several coral genomes have been sequenced allowing a better understanding of these symbiotic organisms threatened by climate change. Scleractinian corals are reef builders and are central to coral reef ecosystems, providing habitat to a great diversity of species.

RESULTS

In the frame of the Tara Pacific expedition, we assemble two coral genomes, Porites lobata and Pocillopora cf. effusa, with vastly improved contiguity that allows us to study the functional organization of these genomes. We annotate their gene catalog and report a relatively higher gene number than that found in other public coral genome sequences, 43,000 and 32,000 genes, respectively. This finding is explained by a high number of tandemly duplicated genes, accounting for almost a third of the predicted genes. We show that these duplicated genes originate from multiple and distinct duplication events throughout the coral lineage. They contribute to the amplification of gene families, mostly related to the immune system and disease resistance, which we suggest to be functionally linked to coral host resilience.

CONCLUSIONS

At large, we show the importance of duplicated genes to inform the biology of reef-building corals and provide novel avenues to understand and screen for differences in stress resilience.

Structural analysis of MALAT1 long noncoding RNA in cells and in evolution

Although not canonically polyadenylated, the long noncoding RNA MALAT1 (metastasis-associated lung adenocarcinoma transcript 1) is stabilized by a highly conserved 76-nt triple helix structure on its 3' end. The entire MALAT1 transcript is over 8000 nt long in humans. The strongest structural conservation signal in MALAT1 (as measured by covariation of base pairs) is in the triple helix structure. Primary sequence analysis of covariation alone does not reveal the degree of structural conservation of the entire full-length transcript, however. Furthermore, RNA structure is often context dependent; RNA binding proteins that are differentially expressed in different cell types may alter structure. We investigate here the in-cell and cell-free structures of the full-length human and green monkey (Chlorocebus sabaeus) MALAT1 transcripts in multiple tissue-derived cell lines using SHAPE chemical probing. Our data reveal levels of uniform structural conservation in different cell lines, in cells and cell-free, and even between species, despite significant differences in primary sequence. The uniformity of the structural conservation across the entire transcript suggests that, despite seeing covariation signals only in the triple helix junction of the lncRNA, the rest of the transcript's structure is remarkably conserved, at least in primates and across multiple cell types and conditions.

Long non-coding RNAs: Fundamental regulators and emerging targets of cancer stem cells

Cancer is one of the leading causes of death worldwide, primarily due to the dearth of efficient therapies that result in long-lasting remission. This is especially true in cases of metastatic cancer where drug resistance causes the disease to recur after treatment. One of the factors contributing to drug resistance, metastasis, and aggressiveness of the cancer is cancer stem cells (CSCs) or tumor-initiating cells. As a result, CSCs have emerged as a potential target for drug development. In the present review, we have examined and highlighted the lncRNAs with their regulatory functions specific to CSCs. Moreover, we have discussed the difficulties and various methods involved in identifying lncRNAs that can play a particular role in regulating and maintaining CSCs. Interestingly, this review only focuses on those lncRNAs with strong functional evidence for CSC specificity and the mechanistic role that allows them to be CSC regulators and be the focus of CSC-specific drug development.

Splicing complexity as a pivotal feature of alternative exons in mammalian species

BACKGROUND

As a significant process of post-transcriptional gene expression regulation in eukaryotic cells, alternative splicing (AS) of exons greatly contributes to the complexity of the transcriptome and indirectly enriches the protein repertoires. A large number of studies have focused on the splicing inclusion of alternative exons and have revealed the roles of AS in organ development and maturation. Notably, AS takes place through a change in the relative abundance of the transcript isoforms produced by a single gene, meaning that exons can have complex splicing patterns. However, the commonly used percent spliced-in (Ψ) values only define the usage rate of exons, but lose information about the complexity of exons' linkage pattern. To date, the extent and functional consequence of splicing complexity of alternative exons in development and evolution is poorly understood.

RESULTS

By comparing splicing complexity of exons in six tissues (brain, cerebellum, heart, liver, kidney, and testis) from six mammalian species (human, chimpanzee, gorilla, macaque, mouse, opossum) and an outgroup species (chicken), we revealed that exons with high splicing complexity are prevalent in mammals and are closely related to features of genes. Using traditional machine learning and deep learning methods, we found that the splicing complexity of exons can be moderately predicted with features derived from exons, among which length of flanking exons and splicing strength of downstream/upstream splice sites are top predictors. Comparative analysis among human, chimpanzee, gorilla, macaque, and mouse revealed that, alternative exons tend to evolve to an increased level of splicing complexity and higher tissue specificity in splicing complexity. During organ development, not only developmentally regulated exons, but also 10-15% of non-developmentally regulated exons show dynamic splicing complexity.

CONCLUSIONS

Our analysis revealed that splicing complexity is an important metric to characterize the splicing dynamics of alternative exons during the development and evolution of mammals.

A hybrid approach to assess the structural impact of long noncoding RNA mutations uncovers key NEAT1 interactions in colorectal cancer

Long noncoding RNAs (lncRNAs) are emerging players in cancer and they entail potential as prognostic biomarkers or therapeutic targets. Earlier studies have identified somatic mutations in lncRNAs that are associated with tumor relapse after therapy, but the underlying mechanisms behind these associations remain unknown. Given the relevance of secondary structure for the function of some lncRNAs, some of these mutations may have a functional impact through structural disturbance. Here, we examined the potential structural and functional impact of a novel A > G point mutation in NEAT1 that has been recurrently observed in tumors of colorectal cancer patients experiencing relapse after treatment. Here, we used the nextPARS structural probing approach to provide first empirical evidence that this mutation alters NEAT1 structure. We further evaluated the potential effects of this structural alteration using computational tools and found that this mutation likely alters the binding propensities of several NEAT1-interacting miRNAs. Differential expression analysis on these miRNA networks shows upregulation of Vimentin, consistent with previous findings. We propose a hybrid pipeline that can be used to explore the potential functional effects of lncRNA somatic mutations.

Mechanisms of Immune-Related Long Non-Coding RNAs in Spleens of Mice Vaccinated with 23-Valent Pneumococcal Polysaccharide Vaccine (PPV23)

The 23-valent pneumococcal vaccine (PPV23) is a classical common vaccine used to prevent pneumococcal disease. In past decades, it was thought that vaccination with this vaccine induces humoral immunity, thereby reducing the disease associated with infection with 23 common serotypes of Streptococcus pneumoniae (Sp). However, for this polysaccharide vaccine, the mechanism of immune response at the transcriptional level has not been fully studied. To identify the lncRNAs (long noncoding RNAs) and mRNAs in spleens related to immunity after PPV23 vaccination in mice, high-throughput RNA sequencing of spleens between a PPV23 treatment group and a control group were performed and evaluated in this study. The RNA-seq results identified a total of 41,321 mRNAs and 34,375 lncRNAs, including 55 significantly differentially expressed (DE) mRNAs and 389 DE lncRNAs (p < 0.05) between the two groups. GO and KEGG annotation analysis indicated that the target genes of DE lncRNAs and DE mRNAs were related to T-cell costimulation, positive regulation of alpha–beta T-cell differentiation, the CD86 biosynthetic process, and the PI3K-Akt signaling pathway, indicating that the polysaccharide component antigens of PPV23 might activate a cellular immune response during the PPV23 immunization process. Moreover, we found that Trim35 (tripartite motif containing 35), a target gene of lncRNA MSTRG.9127, was involved in regulating immunity. Our study provides a catalog of lncRNAs and mRNAs associated with immune cells’ proliferation and differentiation, and they deserve further study to deepen the understanding of the biological processes in the regulation of PPV23 during humoral immunity and cellular immunity.

Evolution of gene regulatory network of C4 photosynthesis in the genus Flaveria reveals the evolutionary status of C3-C4 intermediate species

C₄ photosynthesis evolved from ancestral C₃ photosynthesis by recruiting pre-existing genes to fulfill new functions. The enzymes and transporters required for the C₄ metabolic pathway have been intensively studied and well documented; however, the transcription factors (TFs) that regulate these C₄ metabolic genes are not yet well understood. In particular, how the TF regulatory network of C₄ metabolic genes was rewired during the evolutionary process is unclear. Here, we constructed gene regulatory networks (GRNs) for four closely evolutionarily related species from the genus Flaveria, which represent four different evolutionary stages of C₄ photosynthesis: C₃ (F. robusta), type I C₃-C₄ (F. sonorensis), type II C₃-C₄ (F. ramosissima), and C₄ (F. trinervia). Our results show that more than half of the co-regulatory relationships between TFs and core C₄ metabolic genes are species specific. The counterparts of the C₄ genes in C₃ species were already co-regulated with photosynthesis-related genes, whereas the required TFs for C₄ photosynthesis were recruited later. The TFs involved in C₄ photosynthesis were widely recruited in the type I C₃-C₄ species; nevertheless, type II C₃-C₄ species showed a divergent GRN from C₄ species. In line with these findings, a ¹³CO₂ pulse-labeling experiment showed that the CO₂ initially fixed into C₄ acid was not directly released to the Calvin-Benson-Bassham cycle in the type II C₃-C₄ species. Therefore, our study uncovered dynamic changes in C₄ genes and TF co-regulation during the evolutionary process; furthermore, we showed that the metabolic pathway of the type II C₃-C₄ species F. ramosissima represents an alternative evolutionary solution to the ammonia imbalance in C₃-C₄ intermediate species.

Activity-Dependent Non-Coding RNA MAPK Interactome of the Human Epileptic Brain

The human brain has evolved to have extraordinary capabilities, enabling complex behaviors. The uniqueness of the human brain is increasingly posited to be due in part to the functions of primate-specific, including human-specific, long non-coding RNA (lncRNA) genes, systemically less conserved than protein-coding genes in evolution. Patients who have surgery for drug-resistant epilepsy are subjected to extensive electrical recordings of the brain tissue that is subsequently removed in order to treat their epilepsy. Precise localization of brain tissues with distinct electrical properties offers a rare opportunity to explore the effects of brain activity on gene expression. Here, we identified 231 co-regulated, activity-dependent lncRNAs within the human MAPK signaling cascade. Six lncRNAs, four of which were antisense to known protein-coding genes, were further examined because of their high expression and potential impact on the disease phenotype. Using a model of repeated depolarizations in human neuronal-like cells (Sh-SY5Y), we show that five out of six lncRNAs were electrical activity-dependent, with three of four antisense lncRNAs having reciprocal expression patterns relative to their protein-coding gene partners. Some were directly regulated by MAPK signaling, while others effectively downregulated the expression of the protein-coding genes encoded on the opposite strands of their genomic loci. These lncRNAs, therefore, likely contribute to highly evolved and primate-specific human brain regulatory functions that could be therapeutically modulated to treat epilepsy.

The molecular evolution of spermatogenesis across mammals

The testis produces gametes through spermatogenesis and evolves rapidly at both the morphological and molecular level in mammals^1-6, probably owing to the evolutionary pressure on males to be reproductively successful⁷. However, the molecular evolution of individual spermatogenic cell types across mammals remains largely uncharacterized. Here we report evolutionary analyses of single-nucleus transcriptome data for testes from 11 species that cover the three main mammalian lineages (eutherians, marsupials and monotremes) and birds (the evolutionary outgroup), and include seven primates. We find that the rapid evolution of the testis was driven by accelerated fixation rates of gene expression changes, amino acid substitutions and new genes in late spermatogenic stages, probably facilitated by reduced pleiotropic constraints, haploid selection and transcriptionally permissive chromatin. We identify temporal expression changes of individual genes across species and conserved expression programs controlling ancestral spermatogenic processes. Genes predominantly expressed in spermatogonia (germ cells fuelling spermatogenesis) and Sertoli (somatic support) cells accumulated on X chromosomes during evolution, presumably owing to male-beneficial selective forces. Further work identified transcriptomal differences between X- and Y-bearing spermatids and uncovered that meiotic sex-chromosome inactivation (MSCI) also occurs in monotremes and hence is common to mammalian sex-chromosome systems. Thus, the mechanism of meiotic silencing of unsynapsed chromatin, which underlies MSCI, is an ancestral mammalian feature. Our study illuminates the molecular evolution of spermatogenesis and associated selective forces, and provides a resource for investigating the biology of the testis across mammals.

A large-scale benchmark study of tools for the classification of protein-coding and non-coding RNAs

Identification of protein-coding and non-coding transcripts is paramount for understanding their biological roles. Computational approaches have been addressing this task for over a decade; however, generalized and high-performance models are still unreliable. This benchmark study assessed the performance of 24 tools producing >55 models on the datasets covering a wide range of species. We have collected 135 small and large transcriptomic datasets from existing studies for comparison and identified the potential bottlenecks hampering the performance of current tools. The key insights of this study include lack of standardized training sets, reliance on homogeneous training data, gradual changes in annotated data, lack of augmentation with homology searches, the presence of false positives and negatives in datasets and the lower performance of end-to-end deep learning models. We also derived a new dataset, RNAChallenge, from the benchmark considering hard instances that may include potential false alarms. The best and least well performing models under- and overfit the dataset, respectively, thereby serving a dual purpose. For computational approaches, it will be valuable to develop accurate and unbiased models. The identification of false alarms will be of interest for genome annotators, and experimental study of hard RNAs will help to untangle the complexity of the RNA world.

Evolution and function of developmentally dynamic pseudogenes in mammals

BACKGROUND

Pseudogenes are excellent markers for genome evolution, which are emerging as crucial regulators of development and disease, especially cancer. However, systematic functional characterization and evolution of pseudogenes remain largely unexplored.

RESULTS

To systematically characterize pseudogenes, we date the origin of human and mouse pseudogenes across vertebrates and observe a burst of pseudogene gain in these two lineages. Based on a hybrid sequencing dataset combining full-length PacBio sequencing, sample-matched Illumina sequencing, and public time-course transcriptome data, we observe that abundant mammalian pseudogenes could be transcribed, which contribute to the establishment of organ identity. Our analyses reveal that developmentally dynamic pseudogenes are evolutionarily conserved and show an increasing weight during development. Besides, they are involved in complex transcriptional and post-transcriptional modulation, exhibiting the signatures of functional enrichment. Coding potential evaluation suggests that 19% of human pseudogenes could be translated, thus serving as a new way for protein innovation. Moreover, pseudogenes carry disease-associated SNPs and conduce to cancer transcriptome perturbation.

CONCLUSIONS

Our discovery reveals an unexpectedly high abundance of mammalian pseudogenes that can be transcribed and translated, and these pseudogenes represent a novel regulatory layer. Our study also prioritizes developmentally dynamic pseudogenes with signatures of functional enrichment and provides a hybrid sequencing dataset for further unraveling their biological mechanisms in organ development and carcinogenesis in the future.

Signatures of Convergent Evolution and Natural Selection at the Alcohol Dehydrogenase Gene Region are Correlated with Agriculture in Ethnically Diverse Africans

The alcohol dehydrogenase (ADH) family of genes encodes enzymes that catalyze the metabolism of ethanol into acetaldehyde. Nucleotide variation in ADH genes can affect the catalytic properties of these enzymes and is associated with a variety of traits, including alcoholism and cancer. Some ADH variants, including the ADH1B*48His (rs1229984) mutation in the ADH1B gene, reduce the risk of alcoholism and are under positive selection in multiple human populations. The advent of Neolithic agriculture and associated increase in fermented foods and beverages is hypothesized to have been a selective force acting on such variants. However, this hypothesis has not been tested in populations outside of Asia. Here, we use genome-wide selection scans to show that the ADH gene region is enriched for variants showing strong signals of positive selection in multiple Afroasiatic-speaking, agriculturalist populations from Ethiopia, and that this signal is unique among sub-Saharan Africans. We also observe strong selection signals at putatively functional variants in nearby lipid metabolism genes, which may influence evolutionary dynamics at the ADH region. Finally, we show that haplotypes carrying these selected variants were introduced into Northeast Africa from a West-Eurasian source within the last ∼2,000 years and experienced positive selection following admixture. These selection signals are not evident in nearby, genetically similar populations that practice hunting/gathering or pastoralist subsistence lifestyles, supporting the hypothesis that the emergence of agriculture shapes patterns of selection at ADH genes. Together, these results enhance our understanding of how adaptations to diverse environments and diets have influenced the African genomic landscape.

Human-specific gene CT47 blocks PRMT5 degradation to lead to meiosis arrest

Exploring the functions of human-specific genes (HSGs) is challenging due to the lack of a tractable genetic model system. Testosterone is essential for maintaining human spermatogenesis and fertility, but the underlying mechanism is unclear. Here, we identified Cancer/Testis Antigen gene family 47 (CT47) as an essential regulator of human-specific spermatogenesis by stabilizing arginine methyltransferase 5 (PRMT5). A humanized mouse model revealed that CT47 functions to arrest spermatogenesis by interacting with and regulating CT47/PRMT5 accumulation in the nucleus during the leptotene/zygotene-to-pachytene transition of meiosis. We demonstrate that testosterone induces nuclear depletion of CT47/PRMT5 and rescues leptotene-arrested spermatocyte progression in humanized testes. Loss of CT47 in human embryonic stem cells (hESCs) by CRISPR/Cas9 led to an increase in haploid cells but blocked the testosterone-induced increase in haploid cells when hESCs were differentiated into haploid spermatogenic cells. Moreover, CT47 levels were decreased in nonobstructive azoospermia. Together, these results established CT47 as a crucial regulator of human spermatogenesis by preventing meiosis initiation before the testosterone surge.

Understanding of the various aspects of gene regulatory networks related to crop improvement

The hierarchical relationship between transcription factors, associated proteins, and their target genes is defined by a gene regulatory network (GRN). GRNs allow us to understand how the genotype and environment of a plant are incorporated to control the downstream physiological responses. During plant growth or environmental acclimatization, GRNs are diverse and can be differently regulated across tissue types and organs. An overview of recent advances in the development of GRN that speed up basic and applied plant research is given here. Furthermore, the overview of genome and transcriptome involving GRN research along with the exciting advancement and application are discussed. In addition, different approaches to GRN predictions were elucidated. In this review, we also describe the role of GRN in crop improvement, crop plant manipulation, stress responses, speed breeding and identifying genetic variations/locus. Finally, the challenges and prospects of GRN in plant biology are discussed.

The Long Non-Coding Antisense RNA JHDM1D-AS1 Regulates Inflammatory Responses in Human Monocytes

Monocytes are key players in innate immunity, with their ability to regulate inflammatory responses and combat invading pathogens. There is a growing body of evidence indicating that long non-coding RNA (lncRNA) participate in various cellular biological processes, including the innate immune response. The immunoregulatory properties of numerous lncRNAs discovered in monocytes remain largely unexplored. Here, by RNA sequencing, we identified a lncRNA JHDM1D-AS1, which was upregulated in blood monocytes obtained from patients with sepsis relative to healthy controls. JHDM1D-AS1 expression was induced in primary human monocytes exposed to Toll-like receptor ligands, such as lipopolysaccharide (LPS), or bacteria. The inducibility of JHDM1D-AS1 expression in monocytes depended, at least in part, on nuclear factor-κB activation. JHDM1D-AS1 knockdown experiments in human monocyte-derived macrophages revealed significantly enhanced expression of inflammatory mediators, before and after exposure to LPS, relative to control cells. Specifically, genes involved in inflammatory responses were upregulated (e.g., CXCL2, CXCL8, IL1RN, TREM1, TNF, and IL6), whereas genes involved in anti-inflammatory pathways were downregulated (e.g., SOCS1 and IL10RA). JHDM1D-AS1 overexpression in a pro-monocytic cell line revealed diminished pro-inflammatory responses subsequent to LPS challenge. Collectively, these findings identify JHDM1D-AS1 as a potential anti-inflammatory mediator induced in response to inflammatory stimuli.

Evolution of EGF-like and Zona pellucida domains containing shell matrix proteins in mollusks

Several types of shell matrix proteins (SMPs) have been identified in molluskan shells. Their diversity is the consequence of various molecular processes, including domain shuffling and gene duplication. However, the evolutionary origin of most SMPs remains unclear. In this study, we investigated the evolutionary process EGF-like and zona pellucida (ZP) domains containing SMPs. Two types of the proteins (EGF-like protein (EGFL) and EGF-like and ZP domains containing protein (EGFZP)) were found in the pearl oyster, Pinctada fucata. In contrast, only EGFZP was identified in the gastropods. Phylogenetic analysis and genomic arrangement studies showed that EGFL and EGFZP formed a clade in bivalves, and their encoding genes were localized in tandem repeats on the same scaffold. In P. fucata, EGFL genes were expressed in the outer part of mantle epithelial cells are related to the calcitic shell formation. However, in both P. fucata and the limpet Nipponacmea fuscoviridis, EGFZP genes were expressed in the inner part of the mantle epithelial cells are related to aragonitic shell formation. Furthermore, our analysis showed that in P. fucata, the ZP domain interacts with eight SMPs that have various functions in the nacreous shell mineralization. The data suggest that the ZP domain can interact with other SMPs, and EGFL evolution in pterimorph bivalves represents an example of neo-functionalization that involves the acquisition of a novel protein through gene duplication.

Detecting signatures of selection on gene expression

A substantial amount of phenotypic diversity results from changes in gene expression levels and patterns. Understanding how the transcriptome evolves is therefore a key priority in identifying mechanisms of adaptive change. However, in contrast to powerful models of sequence evolution, we lack a consensus model of gene expression evolution. Furthermore, recent work has shown that many of the comparative approaches used to study gene expression are subject to biases that can lead to false signatures of selection. Here we first outline the main approaches for describing expression evolution and their inherent biases. Next, we bridge the gap between the fields of phylogenetic comparative methods and transcriptomics to reinforce the main pitfalls of inferring selection on expression patterns and use simulation studies to show that shifts in tissue composition can heavily bias inferences of selection. We close by highlighting the multi-dimensional nature of transcriptional variation and identifying major unanswered questions in disentangling how selection acts on the transcriptome.

MicroRNAs and Efferocytosis: Implications for Diagnosis and Therapy

About 10-100 billion cells are generated in the human body in a day, and accordingly, 10-100 billion cells predominantly die for maintaining homeostasis. Dead cells generated by apoptosis are also rapidly engulfed by macrophages (Mθs) to be degraded. In case of the inefficient engulfment of apoptotic cells (ACs) via Mθs, they experience secondary necrosis and thus release intracellular materials, which display damage-associated molecular patterns (DAMPs) and result in diseases. Over the last decades, researchers have also reflected on the significant contribution of microRNAs (miRNAs) to autoimmune diseases through the regulation of Mθs functions. Moreover, miRNAs have shown intricate involvement with completely adjusting basic Mθs functions, such as phagocytosis, inflammation, efferocytosis, tumor promotion, and tissue repair. In this review, the mechanism of efferocytosis containing "Find-Me", "Eat-Me", and "Digest-Me" signals is summarized and the biogenesis of miRNAs is briefly described. Finally, the role of miRNAs in efferocytosis is discussed. It is concluded that miRNAs represent promising treatments and diagnostic targets in impaired phagocytic clearance, which leads to different diseases.

CircTHBS1 drives gastric cancer progression by increasing INHBA mRNA expression and stability in a ceRNA- and RBP-dependent manner

Circular RNAs (circRNAs) play vital regulatory roles in the progression of multiple cancers. In our study, transcriptome analysis and self-organizing maps (SOM) were applied to screen backbone circRNAs in gastric cancer (GC). Upon validation of the expression patterns of screened circRNAs, gain- and loss-of-function assays were performed in vitro and in vivo. Underlying mechanisms were investigated using RNA pull-down, luciferase reporter assay and RNA immunoprecipitation. The expression of circTHBS1 was significantly increased in GC and associated with poor prognosis. CircTHBS1 facilitated the malignant behavior and epithelial-to-mesenchymal transition of GC cells. Mechanistically, circTHBS1 sponged miR-204-5p to promote the expression of Inhibin Subunit Beta A (INHBA). Moreover, circTHBS1 could enhance the HuR-mediated mRNA stability of INHBA, which subsequently activated the TGF-β pathway. Our research identified circTHBS1 as an oncogenic circRNA that enhances GC malignancy by elevating INHBA expression, providing new insight and a feasible target for the diagnosis and treatment of GC.

Comparative Transcriptome Analyses of Gayal (Bos frontalis), Yak (Bos grunniens), and Cattle (Bos taurus) Reveal the High-Altitude Adaptation

Gayal and yak are well adapted to their local high-altitude environments, yet the transcriptional regulation difference of the plateau environment among them remains obscure. Herein, cross-tissue and cross-species comparative transcriptome analyses were performed for the six hypoxia-sensitive tissues from gayal, yak, and cattle. Gene expression profiles for all single-copy orthologous genes showed tissue-specific expression patterns. By differential expression analysis, we identified 3,020 and 1,995 differentially expressed genes (DEGs) in at least one tissue of gayal vs. cattle and yak vs. cattle, respectively. Notably, we found that the adaptability of the gayal to the alpine canyon environment is highly similar to the yak living in the Qinghai-Tibet Plateau, such as promoting red blood cell development, angiogenesis, reducing blood coagulation, immune system activation, and energy metabolism shifts from fatty acid β-oxidation to glycolysis. By further analyzing the common and unique DEGs in the six tissues, we also found that numerous expressed regulatory genes related to these functions are unique in the gayal and yak, which may play important roles in adapting to the corresponding high-altitude environment. Combined with WGCNA analysis, we found that UQCRC1 and COX5A are the shared differentially expressed hub genes related to the energy supply of myocardial contraction in the heart-related modules of gayal and yak, and CAPS is a shared differential hub gene among the hub genes of the lung-related module, which is related to pulmonary artery smooth muscle contraction. Additionally, EDN3 is the unique differentially expressed hub gene related to the tracheal epithelium and pulmonary vasoconstriction in the lung of gayal. CHRM2 is a unique differentially expressed hub gene that was identified in the heart of yak, which has an important role in the autonomous regulation of the heart. These results provide a basis for further understanding the complex transcriptome expression pattern and the regulatory mechanism of high-altitude domestication of gayal and yak.

Identification of Long Non-coding RNA Transcripts in Glycyrrhiza uralensis

Background:

Chinese liquorice (Glycyrrhiza uralensis), an important medicinal plant, contains various valuable secondary metabolites. Secondary metabolites biosynthesis is very tightly regulated; therefore, elucidation and manipulation of the biosynthetic pathways are of great interest. Recent studies have shown that lncRNAs play important regulatory roles in many biological processes, thus identification and modification of their expression is essential to metabolic pathways for biosynthesis of secondary metabolites.

Objectives:

In this study we attempted to identify non-coding RNA transcripts (lncRNAs) that may act as important regulators of diverse biological processes, including stress responses and developmental programs in Glycyrrhiza uralensis.

Materials and Methods:

Identification of potential lncRNAs in Chinese liquorice was performed using a bioinformatics pipeline from the available EST dataset of G. uralensis.

Results:

Bioinformatics analysis revealed that 1365 identical sequences in the range of 200 to 1286 base pair are putative lncRNAs. Only less than one percent of the predicted lncRNAs display sequence conservation with lncRNAs from other species. Moreover, 13 lncRNAs were detected as the potential precursors of 16 miRNAs. From this analysis, we also detected possible target genes of 16 known miRNA genes. The majority of the predicted miRNA target genes have important role in response to plant disease and a couple of them contribute to signalling and metabolic pathways.

Conclusion:

This study demonstrates the existence of lncRNAs in G. uralensis which has not been found before and provides valuable resources for further understanding and characterizing of lncRNAs and also a basis for additional investigation to reveal specific roles of lncRNAs in various biological processes and particularly in response to plant diseases.

Long-range promoter-enhancer contacts are conserved during evolution and contribute to gene expression robustness

Gene expression is regulated through complex molecular interactions, involving cis-acting elements that can be situated far away from their target genes. Data on long-range contacts between promoters and regulatory elements are rapidly accumulating. However, it remains unclear how these regulatory relationships evolve and how they contribute to the establishment of robust gene expression profiles. Here, we address these questions by comparing genome-wide maps of promoter-centered chromatin contacts in mouse and human. We show that there is significant evolutionary conservation of cis-regulatory landscapes, indicating that selective pressures act to preserve not only regulatory element sequences but also their chromatin contacts with target genes. The extent of evolutionary conservation is remarkable for long-range promoter–enhancer contacts, illustrating how the structure of regulatory landscapes constrains large-scale genome evolution. We show that the evolution of cis-regulatory landscapes, measured in terms of distal element sequences, synteny, or contacts with target genes, is significantly associated with gene expression evolution.

Processed pseudogenes: A substrate for evolutionary innovation: Retrotransposition contributes to genome evolution by propagating pseudogene sequences with rich regulatory potential throughout the genome

Processed pseudogenes may serve as a genetic reservoir for evolutionary innovation. Here, we argue that through the activity of long interspersed element-1 retrotransposons, processed pseudogenes disperse coding and noncoding sequences rich with regulatory potential throughout the human genome. While these sequences may appear to be non-functional, a lack of contemporary function does not prohibit future development of biological activity. Here, we discuss the dynamic evolution of certain processed pseudogenes into coding and noncoding genes and regulatory elements, and their implication in wide-ranging biological and pathological processes. Also see the video abstract here: https://youtu.be/iUY_mteVoPI.

Gene expression levels modulate germline mutation rates through the compound effects of transcription-coupled repair and damage

Of all mammalian organs, the testis has long been observed to have the most diverse gene expression profile. To account for this widespread gene expression, we have proposed a mechanism termed 'transcriptional scanning', which reduces germline mutation rates through transcription-coupled repair (TCR). Our hypothesis contrasts with an earlier observation that mutation rates are overall positively correlated with gene expression levels in yeast, implying that transcription is mutagenic due to effects dominated by transcription-coupled damage (TCD). Here we report evidence that the compound effects of both TCR and TCD during spermatogenesis modulate human germline mutation rates, with TCR dominating in most genes, thus supporting the transcriptional scanning hypothesis. Our analyses address potentially confounding factors, distinguish the differential mutagenic effects acting on the highly expressed genes and the low-to-moderately expressed genes, and resolve concerns relating to the validation of the results using a de novo mutation dataset. We also discuss the theoretical possibility of transcriptional scanning hypothesis from an evolutionary perspective. Together, these analyses support a model by which the coupling of transcription-coupled repair and damage establishes the pattern of germline mutation rates and provide an evolutionary explanation for widespread gene expression during spermatogenesis.

A putative de novo evolved gene required for spermatid chromatin condensation in Drosophila melanogaster

Comparative genomics has enabled the identification of genes that potentially evolved de novo from non-coding sequences. Many such genes are expressed in male reproductive tissues, but their functions remain poorly understood. To address this, we conducted a functional genetic screen of over 40 putative de novo genes with testis-enriched expression in Drosophila melanogaster and identified one gene, atlas, required for male fertility. Detailed genetic and cytological analyses showed that atlas is required for proper chromatin condensation during the final stages of spermatogenesis. Atlas protein is expressed in spermatid nuclei and facilitates the transition from histone- to protamine-based chromatin packaging. Complementary evolutionary analyses revealed the complex evolutionary history of atlas. The protein-coding portion of the gene likely arose at the base of the Drosophila genus on the X chromosome but was unlikely to be essential, as it was then lost in several independent lineages. Within the last ~15 million years, however, the gene moved to an autosome, where it fused with a conserved non-coding RNA and evolved a non-redundant role in male fertility. Altogether, this study provides insight into the integration of novel genes into biological processes, the links between genomic innovation and functional evolution, and the genetic control of a fundamental developmental process, gametogenesis.

The role of the multifaceted long non-coding RNAs: A nuclear-cytosolic interplay to regulate hyaluronan metabolism

In the extracellular matrix (ECM), the glycosaminoglycan (GAG) hyaluronan (HA) has different physiological roles favouring hydration, elasticity and cell survival. Three different isoforms of HA synthases (HAS1, 2, and 3) are responsible for the production of HA. In several pathologies the upregulation of HAS enzymes leads to an abnormal HA accumulation causing cell dedifferentiation, proliferation and migration thus favouring cancer progression, fibrosis and vascular wall thickening. An intriguing new player in HAS2 gene expression regulation and HA production is the long non-coding RNA (lncRNA) hyaluronan synthase 2 antisense 1 (HAS2-AS1). A significant part of mammalian genomes corresponds to genes that transcribe lncRNAs; they can regulate gene expression through several mechanisms, being involved not only in maintaining the normal homeostasis of cells and tissues, but also in the onset and progression of different diseases, as demonstrated by the increasing number of studies published through the last decades. HAS2-AS1 is no exception: it can be localized both in the nucleus and in the cytosol, regulating cancer cells as well as vascular smooth muscle cells behaviour.

•

Hyaluronan is a component of the extracellular matrix and is synthetised by three isoenzymes named HAS1, 2, and 3.

•

In several pathologies an upregulation of HAS2 leads to an abnormal accumulation of HA.

•

The long non-coding RNA is a new specific epigenetic regulator of HAS2.

•

In the nucleus HAS2-AS1 modulates chromatin structure around HAS2 promoter increasing transcription.

•

In the cytosol, HAS2-AS1 can interact with several miRNAs altering the expression of several genes as well as can stabilise HAS2 mRNA forming RNA: RNA duplex.

From Junk to Function: LncRNAs in CNS Health and Disease

Recent advances in RNA sequencing technologies helped to uncover the existence of tens of thousands of long non-coding RNAs (lncRNAs) that arise from the dark matter of the genome. These lncRNAs were originally thought to be transcriptional noise but an increasing number of studies demonstrate that these transcripts can modulate protein-coding gene expression by a wide variety of transcriptional and post-transcriptional mechanisms. The spatiotemporal regulation of lncRNA expression is particularly evident in the central nervous system, suggesting that they may directly contribute to specific brain processes, including neurogenesis and cellular homeostasis. Not surprisingly, lncRNAs are therefore gaining attention as putative novel therapeutic targets for disorders of the brain. In this review, we summarize the recent insights into the functions of lncRNAs in the brain, their role in neuronal maintenance, and their potential contribution to disease. We conclude this review by postulating how these RNA molecules can be targeted for the treatment of yet incurable neurological disorders.

Gene expression evolution in pattern-triggered immunity within Arabidopsis thaliana and across Brassicaceae species

Plants recognize surrounding microbes by sensing microbe-associated molecular patterns (MAMPs) to activate pattern-triggered immunity (PTI). Despite their significance for microbial control, the evolution of PTI responses remains largely uncharacterized. Here, by employing comparative transcriptomics of six Arabidopsis thaliana accessions and three additional Brassicaceae species to investigate PTI responses, we identified a set of genes that commonly respond to the MAMP flg22 and genes that exhibit species-specific expression signatures. Variation in flg22-triggered transcriptome responses across Brassicaceae species was incongruent with their phylogeny, while expression changes were strongly conserved within A. thaliana. We found the enrichment of WRKY transcription factor binding sites in the 5'-regulatory regions of conserved and species-specific responsive genes, linking the emergence of WRKY-binding sites with the evolution of gene expression patterns during PTI. Our findings advance our understanding of the evolution of the transcriptome during biotic stress.

Structural and Functional Characterization of a Testicular Long Non-coding RNA (4930463O16Rik) Identified in the Meiotic Arrest of the Mouse Topaz1 -/- Testes

Spermatogenesis involves coordinated processes, including meiosis, to produce functional gametes. We previously reported Topaz1 as a germ cell-specific gene highly conserved in vertebrates. Topaz1 knockout males are sterile with testes that lack haploid germ cells because of meiotic arrest after prophase I. To better characterize Topaz1^–/– testes, we used RNA-sequencing analyses at two different developmental stages (P16 and P18). The absence of TOPAZ1 disturbed the expression of genes involved in microtubule and/or cilium mobility, biological processes required for spermatogenesis. Moreover, a quarter of P18 dysregulated genes are long non-coding RNAs (lncRNAs), and three of them are testis-specific and located in spermatocytes, their expression starting between P11 and P15. The suppression of one of them, 4939463O16Rik, did not alter fertility although sperm parameters were disturbed and sperm concentration fell. The transcriptome of P18-4939463O16Rik^–/– testes was altered and the molecular pathways affected included microtubule-based processes, the regulation of cilium movement and spermatogenesis. The absence of TOPAZ1 protein or 4930463O16Rik produced the same enrichment clusters in mutant testes despite a contrasted phenotype on male fertility. In conclusion, although Topaz1 is essential for the meiosis in male germ cells and regulate the expression of numerous lncRNAs, these studies have identified a Topaz1 regulated lncRNA (4930463O16Rik) that is key for both sperm production and motility.

LncRNA BRCAT54 inhibits the tumorigenesis of non-small cell lung cancer by binding to RPS9 to transcriptionally regulate JAK-STAT and calcium pathway genes

OBJECTIVES

Increasing evidence suggest that long non-coding RNAs (lncRNAs) play critical roles in cancers. However, the expression pattern and underlying mechanisms of lncRNAs in non-small cell lung cancer (NSCLC) remain incompletely understood. This study aimed to elucidate the functions and molecular mechanisms of a certain lncRNA in NSCLC.

METHODS

LncRNA microarray was performed to identify differential expressed lncRNAs between pre- and postoperation plasma in NSCLC patients. The expression level of candidate lncRNA in NSCLC tissues, plasma and cells was determined by quantitative real-time PCR (qRT-PCR) and in situ hybridization. The functional roles of lncRNA were assessed in vitro and in vivo. Furthermore, RNA pull-down, RNA immunoprecipitation, microarray, qRT-PCR and rescue assays were conducted to explore the mechanism action of lncRNA in NSCLC cells.

RESULTS

We identified a novel lncRNA (BRCAT54), which was significantly upregulated in preoperative plasma, NSCLC tissues and NSCLC cells, and its higher expression was associated with better prognosis in patients with NSCLC. Overexpression of BRCAT54 inhibited proliferation, migration and activated apoptosis in NSCLC cells. Conversely, knockdown of BRCAT54 reversed the suppressive effects of BRCAT54. Moreover, overexpression of BRCAT54 repressed NSCLC cell growth in vivo. Mechanistically, BRCAT54 directly bound to RPS9. Knockdown of RPS9 substantially reversed the promoting effects of si-BRCAT54 on cell proliferation and enhanced the inhibitive effect of si-BRCAT54 on BRCAT54 expression. In addition, silencing of RPS9 activated JAK-STAT pathway and suppressed calcium signaling pathway gene expressions.

CONCLUSION

This study identified BRCAT54 as a tumor suppressor in NSCLC. Targeting the BRCAT54 and RPS9 feedback loop might be a novel therapeutic strategy for NSCLC.

The Role of microRNAs and Long Non-Coding RNAs in the Regulation of the Immune Response to Mycobacterium tuberculosis Infection

Non-coding RNAs have emerged as critical regulators of the immune response to infection. MicroRNAs (miRNAs) are small non-coding RNAs which regulate host defense mechanisms against viruses, bacteria and fungi. They are involved in the delicate interplay between Mycobacterium tuberculosis, the causative agent of tuberculosis (TB), and its host, which dictates the course of infection. Differential expression of miRNAs upon infection with M. tuberculosis, regulates host signaling pathways linked to inflammation, autophagy, apoptosis and polarization of macrophages. Experimental evidence suggests that virulent M. tuberculosis often utilize host miRNAs to promote pathogenicity by restricting host-mediated antibacterial signaling pathways. At the same time, host- induced miRNAs augment antibacterial processes such as autophagy, to limit bacterial proliferation. Targeting miRNAs is an emerging option for host-directed therapies. Recent studies have explored the role of long non-coding RNA (lncRNAs) in the regulation of the host response to mycobacterial infection. Among other functions, lncRNAs interact with chromatin remodelers to regulate gene expression and also function as miRNA sponges. In this review we attempt to summarize recent literature on how miRNAs and lncRNAs are differentially expressed during the course of M. tuberculosis infection, and how they influence the outcome of infection. We also discuss the potential use of non-coding RNAs as biomarkers of active and latent tuberculosis. Comprehensive understanding of the role of these non-coding RNAs is the first step towards developing RNA-based therapeutics and diagnostic tools for the treatment of TB.

A pig BodyMap transcriptome reveals diverse tissue physiologies and evolutionary dynamics of transcription

A comprehensive transcriptomic survey of pigs can provide a mechanistic understanding of tissue specialization processes underlying economically valuable traits and accelerate their use as a biomedical model. Here we characterize four transcript types (lncRNAs, TUCPs, miRNAs, and circRNAs) and protein-coding genes in 31 adult pig tissues and two cell lines. We uncover the transcriptomic variability among 47 skeletal muscles, and six adipose depots linked to their different origins, metabolism, cell composition, physical activity, and mitochondrial pathways. We perform comparative analysis of the transcriptomes of seven tissues from pigs and nine other vertebrates to reveal that evolutionary divergence in transcription potentially contributes to lineage-specific biology. Long-range promoter–enhancer interaction analysis in subcutaneous adipose tissues across species suggests evolutionarily stable transcription patterns likely attributable to redundant enhancers buffering gene expression patterns against perturbations, thereby conferring robustness during speciation. This study can facilitate adoption of the pig as a biomedical model for human biology and disease and uncovers the molecular bases of valuable traits.

A comprehensive transcriptomic survey of the pig could enable mechanistic understanding of tissue specialization and accelerate its use as a biomedical model. Here the authors characterize four distinct transcript types in 31 adult pig tissues to dissect their distinct structural and transcriptional features and uncover transcriptomic variability related to tissue physiology.

Annotation of snoRNA abundance across human tissues reveals complex snoRNA-host gene relationships

BACKGROUND

Small nucleolar RNAs (snoRNAs) are mid-size non-coding RNAs required for ribosomal RNA modification, implying a ubiquitous tissue distribution linked to ribosome synthesis. However, increasing numbers of studies identify extra-ribosomal roles of snoRNAs in modulating gene expression, suggesting more complex snoRNA abundance patterns. Therefore, there is a great need for mapping the snoRNome in different human tissues as the blueprint for snoRNA functions.

RESULTS

We used a low structure bias RNA-Seq approach to accurately quantify snoRNAs and compare them to the entire transcriptome in seven healthy human tissues (breast, ovary, prostate, testis, skeletal muscle, liver, and brain). We identify 475 expressed snoRNAs categorized in two abundance classes that differ significantly in their function, conservation level, and correlation with their host gene: 390 snoRNAs are uniformly expressed and 85 are enriched in the brain or reproductive tissues. Most tissue-enriched snoRNAs are embedded in lncRNAs and display strong correlation of abundance with them, whereas uniformly expressed snoRNAs are mostly embedded in protein-coding host genes and are mainly non- or anticorrelated with them. Fifty-nine percent of the non-correlated or anticorrelated protein-coding host gene/snoRNA pairs feature dual-initiation promoters, compared to only 16% of the correlated non-coding host gene/snoRNA pairs.

CONCLUSIONS

Our results demonstrate that snoRNAs are not a single homogeneous group of housekeeping genes but include highly regulated tissue-enriched RNAs. Indeed, our work indicates that the architecture of snoRNA host genes varies to uncouple the host and snoRNA expressions in order to meet the different snoRNA abundance levels and functional needs of human tissues.