Evolutionary conservation in noncoding genomic regions

The regulatory landscape of interacting RNA and protein pools in cellular homeostasis and cancer

Biological systems encompass intricate networks governed by RNA-protein interactions that play pivotal roles in cellular functions. RNA and proteins constituting 1.1% and 18% of the mammalian cell weight, respectively, orchestrate vital processes from genome organization to translation. To date, disentangling the functional fraction of the human genome has presented a major challenge, particularly for noncoding regions, yet recent discoveries have started to unveil a host of regulatory functions for noncoding RNAs (ncRNAs). While ncRNAs exist at different sizes, structures, degrees of evolutionary conservation and abundances within the cell, they partake in diverse roles either alone or in combination. However, certain ncRNA subtypes, including those that have been described or remain to be discovered, are poorly characterized given their heterogeneous nature. RNA activity is in most cases coordinated through interactions with RNA-binding proteins (RBPs). Extensive efforts are being made to accurately reconstruct RNA-RBP regulatory networks, which have provided unprecedented insight into cellular physiology and human disease. In this review, we provide a comprehensive view of RNAs and RBPs, focusing on how their interactions generate functional signals in living cells, particularly in the context of post-transcriptional regulatory processes and cancer.

T-2 toxin exposure induces ovarian damage in sows: lncRNA CUFF.253988.1 promotes cell apoptosis by inhibiting the SIRT3/PGC1α pathway

T-2 toxin, a mycotoxin found in foods and feeds, poses a threat to female reproductive health in both humans and animals. LncRNA CUFF.253988.1 (CUFF.253988.1), highly expressed in pigs, has an undisclosed regulatory role. This study aimed to establish a model of T-2 toxin-induced ovarian injury in sows, both in vivo and in vitro, and to explore the regulatory role and potential mechanisms of CUFF.253988.1. The results showed that feeding T-2 toxin-contaminated feed (1 mg/kg) induced ovarian follicle atresia and mitochondrial structural damage, accompanied by a significant upregulation of CUFF.253988.1 expression in the ovaries. Additionally, T-2 toxin inhibited the SIRT3/PGC1-α pathway associated with mitochondrial function. Moreover, T-2 toxin induced cell apoptosis by upregulating the expression of Cyt c, Bax, cleaved-caspase-9, and cleaved-caspase-3 proteins. In T-2 toxin-induced injury to the ovarian granulosa AVG-16 cells at concentrations of 10, 40 and 160 nM, not only were the previously mentioned effects observed, but also a decrease in mitochondrial membrane potential, ATP content, and an elevation in ROS levels. However, downregulating CUFF.253988.1 reversed T-2 toxin's inhibition of the SIRT3/PGC1-α pathway, alleviating mitochondrial dysfunction and reducing cell apoptosis. Notably, this may be attributed to the inhibition of T-2 toxin-induced enrichment of CUFF.253988.1 in mitochondria. In conclusion, CUFF.253988.1 plays a pivotal role in T-2 toxin-induced ovarian damage, operating through the inhibition of the SIRT3/PGC1-α pathway and promotion of cell apoptosis.

A systematic search for RNA structural switches across the human transcriptome

RNA structural switches are key regulators of gene expression in bacteria, but their characterization in Metazoa remains limited. Here, we present SwitchSeeker, a comprehensive computational and experimental approach for systematic identification of functional RNA structural switches. We applied SwitchSeeker to the human transcriptome and identified 245 putative RNA switches. To validate our approach, we characterized a previously unknown RNA switch in the 3' untranslated region of the RORC (RAR-related orphan receptor C) transcript. In vivo dimethyl sulfate (DMS) mutational profiling with sequencing (DMS-MaPseq), coupled with cryogenic electron microscopy, confirmed its existence as two alternative structural conformations. Furthermore, we used genome-scale CRISPR screens to identify trans factors that regulate gene expression through this RNA structural switch. We found that nonsense-mediated messenger RNA decay acts on this element in a conformation-specific manner. SwitchSeeker provides an unbiased, experimentally driven method for discovering RNA structural switches that shape the eukaryotic gene expression landscape.

The 9p21.3 coronary artery disease risk locus drives vascular smooth muscle cells to an osteochondrogenic state

Genome-wide association studies have identified common genetic variants at ~400 human genomic loci linked to coronary artery disease (CAD) susceptibility. Among these genomic regions, the most impactful is the 9p21.3 CAD risk locus, which spans a 60 kb gene desert and encompasses ~80 SNPs in high linkage disequilibrium. Despite nearly two decades since its discovery, the functional mechanism of this genomic region remains incompletely resolved. To investigate the transcriptional gene programs mediated by 9p21.3 risk locus, we applied a model of induced pluripotent stem cells (iPSCs) from risk and non-risk donors at 9p21.3, as well as isogenic lines with a full haplotype deletion. Upon differentiation to vascular smooth muscle cells (VSMC), single-cell transcriptomic profiling demonstrated iPSC-VSMC phenotypes resembling those from native human coronary arteries, confirming the robustness of this model. Remarkably, our analyses revealed that VSMCs harboring the 9p21.3 risk haplotype preferentially adopt an osteochondrogenic state. Importantly, we identified LIMCH1 and CRABP1 as signature genes critical for defining this transcriptional program. Our study provides new insights into the mechanism at the 9p21.3 risk locus and defines its role in driving a disease-prone transcriptional state in VSMCs.

Conserved Noncoding Elements Evolve Around the Same Genes Throughout Metazoan Evolution

Conserved noncoding elements (CNEs) are DNA sequences located outside of protein-coding genes that can remain under purifying selection for up to hundreds of millions of years. Studies in vertebrate genomes have revealed that most CNEs carry out regulatory functions. Notably, many of them are enhancers that control the expression of homeodomain transcription factors and other genes that play crucial roles in embryonic development. To further our knowledge of CNEs in other parts of the animal tree, we conducted a large-scale characterization of CNEs in more than 50 genomes from three of the main branches of the metazoan tree: Cnidaria, Mollusca, and Arthropoda. We identified hundreds of thousands of CNEs and reconstructed the temporal dynamics of their appearance in each lineage, as well as determining their spatial distribution across genomes. We show that CNEs evolve repeatedly around the same genes across the Metazoa, including around homeodomain genes and other transcription factors; they also evolve repeatedly around genes involved in neural development. We also show that transposons are a major source of CNEs, confirming previous observations from vertebrates and suggesting that they have played a major role in wiring developmental gene regulatory mechanisms since the dawn of animal evolution.

Differential Regulation of Male-Hormones-Related Enhancers Revealed by Chromatin Accessibility and Transcriptional Profiles in Pig Liver

Surgical castration can effectively avoid boar taint and improve pork quality by removing the synthesis of androstenone in the testis, thereby reducing its deposition in adipose tissue. The expression of genes involved in testis-derived hormone metabolism was altered following surgical castration, but the upstream regulatory factors and underlying mechanism remain unclear. In this study, we systematically profiled chromatin accessibility and transcriptional dynamics in liver tissue of castrated and intact full-sibling Yorkshire pigs. First, we identified 897 differentially expressed genes and 6864 differential accessible regions (DARs) using RNA- and ATAC-seq. By integrating the RNA- and ATAC-seq results, 227 genes were identified, and a significant positive correlation was revealed between differential gene expression and the ATAC-seq signal. We constructed a transcription factor regulatory network after motif analysis of DARs and identified a candidate transcription factor (TF) SP1 that targeted the HSD3B1 gene, which was responsible for the metabolism of androstenone. Subsequently, we annotated DARs by incorporating H3K27ac ChIP-seq data, marking 2234 typical enhancers and 245 super enhancers involved in the regulation of all testis-derived hormones. Among these, four typical enhancers associated with HSD3B1 were identified. Furthermore, an in-depth investigation was conducted on the androstenone-related enhancers, and an androstenone-related mutation was identified in a newfound candidatetypical enhancer (andEN) with dual-luciferase assays. These findings provide further insights into how enhancers function as links between phenotypic and non-coding area variations. The discovery of upstream TF and enhancers of HSD3B1 contributes to understanding the regulatory networks of androstenone metabolism and provides an important foundation for improving pork quality.

Human exons and introns classification using pre-trained Resnet-50 and GoogleNet models and 13-layers CNN model

BACKGROUND

Examining functions and characteristics of DNA sequences is a highly challenging task. When it comes to the human genome, which is made up of exons and introns, this task is more challenging. Human exons and introns contain millions to billions of nucleotides, which contributes to the complexity observed in this sequences. Considering how complicated the subject of genomics is, it is obvious that using signal processing techniques and deep learning tools to build a strong predictive model can be very helpful for the development of the research of the human genome.

RESULTS

After representing human exons and introns with color images using Frequency Chaos Game Representation, two pre-trained convolutional neural network models (Resnet-50 and GoogleNet) and a proposed CNN model having 13 hidden layers were used to classify our obtained images. We have reached a value of 92% for the accuracy rate for Resnet-50 model in about 7 h for the execution time, a value of 91.5% for the accuracy rate for the GoogleNet model in 2 h and a half for the execution time. For our proposed CNN model, we have reached 91.6% for the accuracy rate in 2 h and 37 min.

CONCLUSIONS

Our proposed CNN model is faster than the Resnet-50 model in terms of execution time. It was able to slightly exceed the GoogleNet model for the accuracy rate value.

Multiplex generation and single cell analysis of structural variants in a mammalian genome

The functional consequences of structural variants (SVs) in mammalian genomes are challenging to study. This is due to several factors, including: 1) their numerical paucity relative to other forms of standing genetic variation such as single nucleotide variants (SNVs) and short insertions or deletions (indels); 2) the fact that a single SV can involve and potentially impact the function of more than one gene and/or cis regulatory element; and 3) the relative immaturity of methods to generate and map SVs, either randomly or in targeted fashion, in in vitro or in vivo model systems. Towards addressing these challenges, we developed Genome-Shuffle-seq, a straightforward method that enables the multiplex generation and mapping of several major forms of SVs (deletions, inversions, translocations) throughout a mammalian genome. Genome-Shuffle-seq is based on the integration of "shuffle cassettes" to the genome, wherein each shuffle cassette contains components that facilitate its site-specific recombination (SSR) with other integrated shuffle cassettes (via Cre-loxP), its mapping to a specific genomic location (via T7-mediated in vitro transcription or IVT), and its identification in single-cell RNA-seq (scRNA-seq) data (via T7-mediated in situ transcription or IST). In this proof-of-concept, we apply Genome-Shuffle-seq to induce and map thousands of genomic SVs in mouse embryonic stem cells (mESCs) in a single experiment. Induced SVs are rapidly depleted from the cellular population over time, possibly due to Cre-mediated toxicity and/or negative selection on the rearrangements themselves. Leveraging T7 IST of barcodes whose positions are already mapped, we further demonstrate that we can efficiently genotype which SVs are present in association with each of many single cell transcriptomes in scRNA-seq data. Finally, preliminary evidence suggests our method may be a powerful means of generating extrachromosomal circular DNAs (ecDNAs). Looking forward, we anticipate that Genome-Shuffle-seq may be broadly useful for the systematic exploration of the functional consequences of SVs on gene expression, the chromatin landscape, and 3D nuclear architecture. We further anticipate potential uses for in vitro modeling of ecDNAs, as well as in paving the path to a minimal mammalian genome.

Exploration of the Noncoding Genome for Human-Specific Therapeutic Targets-Recent Insights at Molecular and Cellular Level

While it is well known that 98-99% of the human genome does not encode proteins, but are nevertheless transcriptionally active and give rise to a broad spectrum of noncoding RNAs [ncRNAs] with complex regulatory and structural functions, specific functions have so far been assigned to only a tiny fraction of all known transcripts. On the other hand, the striking observation of an overwhelmingly growing fraction of ncRNAs, in contrast to an only modest increase in the number of protein-coding genes, during evolution from simple organisms to humans, strongly suggests critical but so far essentially unexplored roles of the noncoding genome for human health and disease pathogenesis. Research into the vast realm of the noncoding genome during the past decades thus lead to a profoundly enhanced appreciation of the multi-level complexity of the human genome. Here, we address a few of the many huge remaining knowledge gaps and consider some newly emerging questions and concepts of research. We attempt to provide an up-to-date assessment of recent insights obtained by molecular and cell biological methods, and by the application of systems biology approaches. Specifically, we discuss current data regarding two topics of high current interest: (1) By which mechanisms could evolutionary recent ncRNAs with critical regulatory functions in a broad spectrum of cell types (neural, immune, cardiovascular) constitute novel therapeutic targets in human diseases? (2) Since noncoding genome evolution is causally linked to brain evolution, and given the profound interactions between brain and immune system, could human-specific brain-expressed ncRNAs play a direct or indirect (immune-mediated) role in human diseases? Synergistic with remarkable recent progress regarding delivery, efficacy, and safety of nucleic acid-based therapies, the ongoing large-scale exploration of the noncoding genome for human-specific therapeutic targets is encouraging to proceed with the development and clinical evaluation of novel therapeutic pathways suggested by these research fields.

Genome-wide mapping of regulatory variants for temperature- and salinity-adaptive genes reveals genetic basis of genotype-by-environment interaction in Crassostrea ariakensis

Regulatory variants in gene expression serve as bridges linking genetic variation and phenotypic plasticity. Environmental conditions typically influence the effects of regulatory variants on phenotypic plasticity; however, such genotype-by-environment interactions (G × E) are poorly understood. This study aimed to investigate the genetic basis of G × E in estuarine oyster (Crassostrea ariakensis), which is an important model animal for studying environmental adaption owing to its high plasticity and large intraspecific divergence. Genome-wide mapping of expression quantitative trait loci (eQTLs) for 23 environmental adaptive genes was performed for 256 estuarine oysters. We identified 1194 eQTL single nucleotide polymorphisms (eSNPs), including 433 cis-eSNPs in four genes and 722 trans-eSNPs in eight genes. The expression variation explanation of cis-eSNPs (9.95%) was significantly higher than that of trans-eSNPs (9.15%). We specifically showed cis- and trans-eSNPs with high linkage disequilibrium (LD) for Traf7, Slc6a5, Ggt, and Dap3. For example, we identified a cis-regulatory LD block containing 68 cis-eSNP and a trans-regulatory LD block, including 20 trans-eSNPs in Traf7. A high proportion (85%) of 40 vital eSNPs exhibited significant G × E effects. We identified crossing and nonparallel interactions of G × E, with the tag cis-eSNPs of Baat and Slc6a5 as representatives. Our results indicated that cis-eQTLs are highly conserved. This study provides insights into the understanding of adaptive evolutionary mechanisms and phenotypic response prediction to variable environments, as well as the genetic improvement for superior adaptive traits for genetic resource conservation and aquaculture.

Functional inference of long non-coding RNAs through exploration of highly conserved regions

Background: Long non-coding RNAs (lncRNAs), which are generally less functionally characterized or less annotated, evolve more rapidly than mRNAs and substantially possess fewer sequence conservation patterns than protein-coding genes across divergent species. People assume that the functional inference could be conducted on the evolutionarily conserved long non-coding RNAs as they are most likely to be functional. In the past decades, substantial progress has been made in discussions on the evolutionary conservation of non-coding genomic regions from multiple perspectives. However, understanding their conservation and the functions associated with sequence conservation in relation to further corresponding phenotypic variability or disorders still remains incomplete. Results: Accordingly, we determined a highly conserved region (HCR) to verify the sequence conservation among long non-coding RNAs and systematically profiled homologous long non-coding RNA clusters in humans and mice based on the detection of highly conserved regions. Moreover, according to homolog clustering, we explored the potential function inference via highly conserved regions on representative long non-coding RNAs. On lncRNA XACT, we investigated the potential functional competence between XACT and lncRNA XIST by recruiting miRNA-29a, regulating the downstream target genes. In addition, on lncRNA LINC00461, we examined the interaction relationship between LINC00461 and SND1. This interaction or association may be perturbed during the progression of glioma. In addition, we have constructed a website with user-friendly web interfaces for searching, analyzing, and downloading to present the homologous clusters of humans and mice. Conclusion: Collectively, homolog clustering via the highly conserved region definition and detection on long non-coding RNAs, as well as the functional explorations on representative sequences in our research, would provide new evidence for the potential function of long non-coding RNAs. Our results on the remarkable roles of long non-coding RNAs would presumably provide a new theoretical basis and candidate diagnostic indicators for tumors.

Molecular hallmarks of long non-coding RNAs in aging and its significant effect on aging-associated diseases

Aging is linked to the deterioration of many physical and cognitive abilities and is the leading risk factor for Alzheimer's disease. The growing aging population is a significant healthcare problem globally that researchers must investigate to better understand the underlying aging processes. Advances in microarrays and sequencing techniques have resulted in deeper analyses of diverse essential genomes (e.g., mouse, human, and rat) and their corresponding cell types, their organ-specific transcriptomes, and the tissue involved in aging. Traditional gene controllers such as DNA- and RNA-binding proteins significantly influence such programs, causing the need to sort out long non-coding RNAs, a new class of powerful gene regulatory elements. However, their functional significance in the aging process and senescence has yet to be investigated and identified. Several recent researchers have associated the initiation and development of senescence and aging in mammals with several well-reported and novel long non-coding RNAs. In this review article, we identified and analyzed the evolving functions of long non-coding RNAs in cellular processes, including cellular senescence, aging, and age-related pathogenesis, which are the major hallmarks of long non-coding RNAs in aging.

Recent Advances in Genetics and Genomics of Snub-Nosed Monkeys (Rhinopithecus) and Their Implications for Phylogeny, Conservation, and Adaptation

The snub-nosed monkey genus Rhinopithecus (Colobinae) comprises five species (Rhinopithecus roxellana, Rhinopithecus brelichi, Rhinopithecus bieti, Rhinopithecus strykeri, and Rhinopithecus avunculus). They are range-restricted species occurring only in small areas in China, Vietnam, and Myanmar. All extant species are listed as endangered or critically endangered by the International Union for Conservation of Nature (IUCN) Red List, all with decreasing populations. With the development of molecular genetics and the improvement and cost reduction in whole-genome sequencing, knowledge about evolutionary processes has improved largely in recent years. Here, we review recent major advances in snub-nosed monkey genetics and genomics and their impact on our understanding of the phylogeny, phylogeography, population genetic structure, landscape genetics, demographic history, and molecular mechanisms of adaptation to folivory and high altitudes in this primate genus. We further discuss future directions in this research field, in particular how genomic information can contribute to the conservation of snub-nosed monkeys.

Nucleotide Composition of Ultra-Conserved Elements Shows Excess of GpC and Depletion of GG and CC Dinucleotides

The public UCNEbase database, comprising 4273 human ultra-conserved noncoding elements (UCNEs), was thoroughly investigated with the aim to find any nucleotide signals or motifs that have made these DNA sequences practically unchanged over three hundred million years of evolution. Each UCNE comprises over 200 nucleotides and has at least 95% identity between humans and chickens. A total of 31,046 SNPs were found within the UCNE database. We demonstrated that every human has over 300 mutations within 4273 UCNEs. No association of UCNEs with non-coding RNAs, nor preference of a particular meiotic recombination rate within them were found. No sequence motifs associated with UCNEs nor their flanking regions have been found. However, we demonstrated that UCNEs have strong nucleotide and dinucleotide sequence abnormalities compared to genome averages. Specifically, UCNEs are depleted for CC and GG dinucleotides, while GC dinucleotides are in excess of 28%. Importantly, GC dinucleotides have extraordinarily strong stacking free-energy inside the DNA helix and unique resistance to dissociation. Based on the adjacent nucleotide stacking abnormalities within UCNEs, we conjecture that peculiarities in dinucleotide distribution within UCNEs may create unique 3D conformation and specificity to bind proteins. We also discuss the strange dynamics of multiple SNPs inside UCNEs and reasons why these sequences are extraordinarily conserved.

Consensus guidelines on the construct validity of rodent models of restless legs syndrome

Our understanding of the causes and natural course of restless legs syndrome (RLS) is incomplete. The lack of objective diagnostic biomarkers remains a challenge for clinical research and for the development of valid animal models. As a task force of preclinical and clinical scientists, we have previously defined face validity parameters for rodent models of RLS. In this article, we establish new guidelines for the construct validity of RLS rodent models. To do so, we first determined and agreed on the risk, and triggering factors and pathophysiological mechanisms that influence RLS expressivity. We then selected 20 items considered to have sufficient support in the literature, which we grouped by sex and genetic factors, iron-related mechanisms, electrophysiological mechanisms, dopaminergic mechanisms, exposure to medications active in the central nervous system, and others. These factors and biological mechanisms were then translated into rodent bioequivalents deemed to be most appropriate for a rodent model of RLS. We also identified parameters by which to assess and quantify these bioequivalents. Investigating these factors, both individually and in combination, will help to identify their specific roles in the expression of rodent RLS-like phenotypes, which should provide significant translational implications for the diagnosis and treatment of RLS.

A clinician’s handbook for using ctDNA throughout the patient journey

Background

The promise of precision cancer medicine presently centers around the genomic sequence of a patient’s tumor being translated into timely, actionable information to inform clinical care. The analysis of cell-free DNA from liquid biopsy, which contains circulating tumor DNA (ctDNA) in patients with cancer, has proven to be amenable to various settings in oncology. However, open questions surrounding the clinical validity and utility of plasma-based analyses have hindered widespread clinical adoption.

Main body

Owing to the rapid evolution of the field, studies supporting the use of ctDNA as a biomarker throughout a patient’s journey with cancer have accumulated in the last few years, warranting a review of the latest status for clinicians who may employ ctDNA in their precision oncology programs. In this work, we take a step back from the intricate coverage of detection approaches described extensively elsewhere and cover basic concepts around the practical implementation of next generation sequencing (NGS)-guided liquid biopsy. We compare relevant targeted and untargeted approaches to plasma DNA analysis, describe the latest evidence for clinical validity and utility, and highlight the value of genome-wide ctDNA analysis, particularly as it relates to early detection strategies and discovery applications harnessing the non-coding genome.

Conclusions

The maturation of liquid biopsy for clinical application will require interdisciplinary efforts to address current challenges. However, patients and clinicians alike may greatly benefit in the future from its incorporation into routine oncology care.