Recombination of repeat elements generates somatic complexity in human genomes.

Non-allelic recombination between homologous repetitive elements contributes to evolution and human genetic disorders. Here, we combine short- and long-DNA read sequencing of repeat elements with a new bioinformatics pipeline to show that somatic recombination of Alu and L1 elements is widespread in the human genome. Our analysis uncovers tissue-specific non-allelic homologous recombination hallmarks; moreover, we find that centromeres and cancer-associated genes are enriched for retroelements that may act as recombination hotspots. We compare recombination profiles in human-induced pluripotent stem cells and differentiated neurons and find that the neuron-specific recombination of repeat elements accompanies chromatin changes during cell-fate determination. Finally, we report that somatic recombination profiles are altered in Parkinson's and Alzheimer's disease, suggesting a link between retroelement recombination and genomic instability in neurodegeneration. This work highlights a significant contribution of the somatic recombination of repeat elements to genomic diversity in health and disease.

Long non-coding RNAs: definitions, functions, challenges and recommendations.

Genes specifying long non-coding RNAs (lncRNAs) occupy a large fraction of the genomes of complex organisms. The term 'lncRNAs' encompasses RNA polymerase I (Pol I), Pol II and Pol III transcribed RNAs, and RNAs from processed introns. The various functions of lncRNAs and their many isoforms and interleaved relationships with other genes make lncRNA classification and annotation difficult. Most lncRNAs evolve more rapidly than protein-coding sequences, are cell type specific and regulate many aspects of cell differentiation and development and other physiological processes. Many lncRNAs associate with chromatin-modifying complexes, are transcribed from enhancers and nucleate phase separation of nuclear condensates and domains, indicating an intimate link between lncRNA expression and the spatial control of gene expression during development. lncRNAs also have important roles in the cytoplasm and beyond, including in the regulation of translation, metabolism and signalling. lncRNAs often have a modular structure and are rich in repeats, which are increasingly being shown to be relevant to their function. In this Consensus Statement, we address the definition and nomenclature of lncRNAs and their conservation, expression, phenotypic visibility, structure and functions. We also discuss research challenges and provide recommendations to advance the understanding of the roles of lncRNAs in development, cell biology and disease.

The Secret Life of lncRNAs: Conserved, yet Not Conserved.

Guo and colleagues discover a new layer of complexity to the lncRNA evolution where positionally conserved lncRNAs in human ESCs are broadly spliced and exported to the cytoplasm contrary to their mouse counterpart that are predominantly unspliced and nuclear retained. Distinct processing leads to species-specific lncRNA function in pluripotency maintenance.

Gut microbial carbohydrate metabolism contributes to insulin resistance.

Insulin resistance is the primary pathophysiology underlying metabolic syndrome and type 2 diabetes1,2. Previous metagenomic studies have described the characteristics of gut microbiota and their roles in metabolizing major nutrients in insulin resistance3-9. In particular, carbohydrate metabolism of commensals has been proposed to contribute up to 10% of the host's overall energy extraction10, thereby playing a role in the pathogenesis of obesity and prediabetes3,4,6. Nevertheless, the underlying mechanism remains unclear. Here we investigate this relationship using a comprehensive multi-omics strategy in humans. We combine unbiased faecal metabolomics with metagenomics, host metabolomics and transcriptomics data to profile the involvement of the microbiome in insulin resistance. These data reveal that faecal carbohydrates, particularly host-accessible monosaccharides, are increased in individuals with insulin resistance and are associated with microbial carbohydrate metabolisms and host inflammatory cytokines. We identify gut bacteria associated with insulin resistance and insulin sensitivity that show a distinct pattern of carbohydrate metabolism, and demonstrate that insulin-sensitivity-associated bacteria ameliorate host phenotypes of insulin resistance in a mouse model. Our study, which provides a comprehensive view of the host-microorganism relationships in insulin resistance, reveals the impact of carbohydrate metabolism by microbiota, suggesting a potential therapeutic target for ameliorating insulin resistance.

The status of the human gene catalogue.

Scientists have been trying to identify every gene in the human genome since the initial draft was published in 2001. In the years since, much progress has been made in identifying protein-coding genes, currently estimated to number fewer than 20,000, with an ever-expanding number of distinct protein-coding isoforms. Here we review the status of the human gene catalogue and the efforts to complete it in recent years. Beside the ongoing annotation of protein-coding genes, their isoforms and pseudogenes, the invention of high-throughput RNA sequencing and other technological breakthroughs have led to a rapid growth in the number of reported non-coding RNA genes. For most of these non-coding RNAs, the functional relevance is currently unclear; we look at recent advances that offer paths forward to identifying their functions and towards eventually completing the human gene catalogue. Finally, we examine the need for a universal annotation standard that includes all medically significant genes and maintains their relationships with different reference genomes for the use of the human gene catalogue in clinical settings.

Endogenous retrotransposition activates oncogenic pathways in hepatocellular carcinoma.

LINE-1 (L1) retrotransposons are mobile genetic elements comprising ~17% of the human genome. New L1 insertions can profoundly alter gene function and cause disease, though their significance in cancer remains unclear. Here, we applied enhanced retrotransposon capture sequencing (RC-seq) to 19 hepatocellular carcinoma (HCC) genomes and elucidated two archetypal L1-mediated mechanisms enabling tumorigenesis. In the first example, 4/19 (21.1%) donors presented germline retrotransposition events in the tumor suppressor mutated in colorectal cancers (MCC). MCC expression was ablated in each case, enabling oncogenic ß-catenin/Wnt signaling. In the second example, suppression of tumorigenicity 18 (ST18) was activated by a tumor-specific L1 insertion. Experimental assays confirmed that the L1 interrupted a negative feedback loop by blocking ST18 repression of its enhancer. ST18 was also frequently amplified in HCC nodules from Mdr2(-/-) mice, supporting its assignment as a candidate liver oncogene. These proof-of-principle results substantiate L1-mediated retrotransposition as an important etiological factor in HCC.

Deep sequencing of short capped RNAs reveals novel families of noncoding RNAs.

In eukaryotes, capped RNAs include long transcripts such as messenger RNAs and long noncoding RNAs, as well as shorter transcripts such as spliceosomal RNAs, small nucleolar RNAs, and enhancer RNAs. Long capped transcripts can be profiled using cap analysis gene expression (CAGE) sequencing and other methods. Here, we describe a sequencing library preparation protocol for short capped RNAs, apply it to a differentiation time course of the human cell line THP-1, and systematically compare the landscape of short capped RNAs to that of long capped RNAs. Transcription initiation peaks associated with genes in the sense direction have a strong preference to produce either long or short capped RNAs, with one out of six peaks detected in the short capped RNA libraries only. Gene-associated short capped RNAs have highly specific 3' ends, typically overlapping splice sites. Enhancers also preferentially generate either short or long capped RNAs, with 10% of enhancers observed in the short capped RNA libraries only. Enhancers producing either short or long capped RNAs show enrichment for GWAS-associated disease SNPs. We conclude that deep sequencing of short capped RNAs reveals new families of noncoding RNAs and elucidates the diversity of transcripts generated at known and novel promoters and enhancers.

Piwil2 (Mili) sustains neurogenesis and prevents cellular senescence in the postnatal hippocampus.

Adult neural progenitor cells (aNPCs) ensure lifelong neurogenesis in the mammalian hippocampus. Proper regulation of aNPC fate has thus important implications for brain plasticity and healthy aging. Piwi proteins and the small noncoding RNAs interacting with them (piRNAs) have been proposed to control memory and anxiety, but the mechanism remains elusive. Here, we show that Piwil2 (Mili) is essential for proper neurogenesis in the postnatal mouse hippocampus. RNA sequencing of aNPCs and their differentiated progeny reveal that Mili and piRNAs are dynamically expressed in neurogenesis. Depletion of Mili and piRNAs in the adult hippocampus impairs aNPC differentiation toward a neural fate, induces senescence, and generates reactive glia. Transcripts modulated upon Mili depletion bear sequences complementary or homologous to piRNAs and include repetitive elements and mRNAs encoding essential proteins for proper neurogenesis. Our results provide evidence of a critical role for Mili in maintaining fitness and proper fate of aNPCs, underpinning a possible involvement of the piRNA pathway in brain plasticity and successful aging.

An atlas of combinatorial transcriptional regulation in mouse and man.

Combinatorial interactions among transcription factors are critical to directing tissue-specific gene expression. To build a global atlas of these combinations, we have screened for physical interactions among the majority of human and mouse DNA-binding transcription factors (TFs). The complete networks contain 762 human and 877 mouse interactions. Analysis of the networks reveals that highly connected TFs are broadly expressed across tissues, and that roughly half of the measured interactions are conserved between mouse and human. The data highlight the importance of TF combinations for determining cell fate, and they lead to the identification of a SMAD3/FLI1 complex expressed during development of immunity. The availability of large TF combinatorial networks in both human and mouse will provide many opportunities to study gene regulation, tissue differentiation, and mammalian evolution.

Antisense RNAs during early vertebrate development are divided in groups with distinct features.

Long noncoding RNAs or lncRNAs are a class of non-protein-coding RNAs that are >200 nt in length. Almost 50% of lncRNAs during zebrafish development are transcribed in an antisense direction to a protein-coding gene. However, the role of these natural antisense transcripts (NATs) during development remains enigmatic. To understand NATs in early vertebrate development, we took a computational biology approach and analyzed existing as well as novel data sets. Our analysis indicates that zebrafish NATs can be divided into two major classes based on their coexpression patterns with respect to the overlapping protein-coding genes. Group 1 NATs have characteristics similar to maternally deposited RNAs in that their levels decrease as development progresses. Group 1 NAT levels are negatively correlated with that of overlapping sense-strand protein-coding genes. Conversely, Group 2 NATs are coexpressed with overlapping protein-coding genes. In contrast to Group 1, which is enriched in genes involved in developmental pathways, Group 2 protein-coding genes are enriched in housekeeping functions. Group 1 NATs also show larger overlap and higher complementarity with the sense-strand mRNAs compared to other NATs. In addition, our transcriptomics data, quantifying RNA levels from cytoplasmic and nuclear compartments, indicates that Group 1 NATs are more abundant in the cytosol. Based on their expression pattern, cytosolic nature, and their higher complementarity to the overlapping developmental mRNAs, we speculate that Group 1 NATs function post-transcriptionally to silence spurious expression of developmental genes.

Embryonic LTR retrotransposons supply promoter modules to somatic tissues.

Long terminal repeat (LTR) retrotransposons are widely distributed across the human genome. They have accumulated through retroviral integration into germline DNA and are latent genetic modules. Active LTR promoters are observed in germline cells; however, little is known about the mechanisms underlying their active transcription in somatic tissues. Here, by integrating our previous transcriptome data set with publicly available data sets, we show that the LTR families MLT2A1 and MLT2A2 are primarily expressed in human four-cell and eight-cell embryos and are also activated in some adult somatic tissues, particularly pineal gland. Three MLT2A elements function as the promoters and first exons of the protein-coding genes ABCE1, COL5A1, and GALNT13 specifically in the pineal gland of humans but not in that of macaques, suggesting that the exaptation of these LTRs as promoters occurred during recent primate evolution. This analysis provides insight into the possible transition from germline insertion to somatic expression of LTR retrotransposons.

A field guide to cultivating computational biology.

Evolving in sync with the computation revolution over the past 30 years, computational biology has emerged as a mature scientific field. While the field has made major contributions toward improving scientific knowledge and human health, individual computational biology practitioners at various institutions often languish in career development. As optimistic biologists passionate about the future of our field, we propose solutions for both eager and reluctant individual scientists, institutions, publishers, funding agencies, and educators to fully embrace computational biology. We believe that in order to pave the way for the next generation of discoveries, we need to improve recognition for computational biologists and better align pathways of career success with pathways of scientific progress. With 10 outlined steps, we call on all adjacent fields to move away from the traditional individual, single-discipline investigator research model and embrace multidisciplinary, data-driven, team science.

Amphioxus functional genomics and the origins of vertebrate gene regulation.

Vertebrates have greatly elaborated the basic chordate body plan and evolved highly distinctive genomes that have been sculpted by two whole-genome duplications. Here we sequence the genome of the Mediterranean amphioxus (Branchiostoma lanceolatum) and characterize DNA methylation, chromatin accessibility, histone modifications and transcriptomes across multiple developmental stages and adult tissues to investigate the evolution of the regulation of the chordate genome. Comparisons with vertebrates identify an intermediate stage in the evolution of differentially methylated enhancers, and a high conservation of gene expression and its cis-regulatory logic between amphioxus and vertebrates that occurs maximally at an earlier mid-embryonic phylotypic period. We analyse regulatory evolution after whole-genome duplications, and find that-in vertebrates-over 80% of broadly expressed gene families with multiple paralogues derived from whole-genome duplications have members that restricted their ancestral expression, and underwent specialization rather than subfunctionalization. Counter-intuitively, paralogues that restricted their expression increased the complexity of their regulatory landscapes. These data pave the way for a better understanding of the regulatory principles that underlie key vertebrate innovations.

Comparative transcriptomics of primary cells in vertebrates.

Gene expression profiles in homologous tissues have been observed to be different between species, which may be due to differences between species in the gene expression program in each cell type, but may also reflect differences in cell type composition of each tissue in different species. Here, we compare expression profiles in matching primary cells in human, mouse, rat, dog, and chicken using Cap Analysis Gene Expression (CAGE) and short RNA (sRNA) sequencing data from FANTOM5. While we find that expression profiles of orthologous genes in different species are highly correlated across cell types, in each cell type many genes were differentially expressed between species. Expression of genes with products involved in transcription, RNA processing, and transcriptional regulation was more likely to be conserved, while expression of genes encoding proteins involved in intercellular communication was more likely to have diverged during evolution. Conservation of expression correlated positively with the evolutionary age of genes, suggesting that divergence in expression levels of genes critical for cell function was restricted during evolution. Motif activity analysis showed that both promoters and enhancers are activated by the same transcription factors in different species. An analysis of expression levels of mature miRNAs and of primary miRNAs identified by CAGE revealed that evolutionary old miRNAs are more likely to have conserved expression patterns than young miRNAs. We conclude that key aspects of the regulatory network are conserved, while differential expression of genes involved in cell-to-cell communication may contribute greatly to phenotypic differences between species.

Recounting the FANTOM CAGE-Associated Transcriptome.

Long noncoding RNAs (lncRNAs) have emerged as key coordinators of biological and cellular processes. Characterizing lncRNA expression across cells and tissues is key to understanding their role in determining phenotypes, including human diseases. We present here FC-R2, a comprehensive expression atlas across a broadly defined human transcriptome, inclusive of over 109,000 coding and noncoding genes, as described in the FANTOM CAGE-Associated Transcriptome (FANTOM-CAT) study. This atlas greatly extends the gene annotation used in the original recount2 resource. We demonstrate the utility of the FC-R2 atlas by reproducing key findings from published large studies and by generating new results across normal and diseased human samples. In particular, we (a) identify tissue-specific transcription profiles for distinct classes of coding and noncoding genes, (b) perform differential expression analysis across thirteen cancer types, identifying novel noncoding genes potentially involved in tumor pathogenesis and progression, and (c) confirm the prognostic value for several enhancer lncRNAs expression in cancer. Our resource is instrumental for the systematic molecular characterization of lncRNA by the FANTOM6 Consortium. In conclusion, comprised of over 70,000 samples, the FC-R2 atlas will empower other researchers to investigate functions and biological roles of both known coding genes and novel lncRNAs.

SCAFE: a software suite for analysis of transcribed cis-regulatory elements in single cells.

MOTIVATION: Cell type-specific activities of cis-regulatory elements (CRE) are central to understanding gene regulation and disease predisposition. Single-cell RNA 5'end sequencing (sc-end5-seq) captures the transcription start sites (TSS) which can be used as a proxy to measure the activity of transcribed CREs (tCREs). However, a substantial fraction of TSS identified from sc-end5-seq data may not be genuine due to various artifacts, hindering the use of sc-end5-seq for de novo discovery of tCREs. RESULTS: We developed SCAFE-Single-Cell Analysis of Five-prime Ends-a software suite that processes sc-end5-seq data to de novo identify TSS clusters based on multiple logistic regression. It annotates tCREs based on the identified TSS clusters and generates a tCRE-by-cell count matrix for downstream analyses. The software suite consists of a set of flexible tools that could either be run independently or as pre-configured workflows. AVAILABILITY AND IMPLEMENTATION: SCAFE is implemented in Perl and R. The source code and documentation are freely available for download under the MIT License from https://github.com/chung-lab/SCAFE. Docker images are available from https://hub.docker.com/r/cchon/scafe. The submitted software version and test data are archived at https://doi.org/10.5281/zenodo.7023163 and https://doi.org/10.5281/zenodo.7024060, respectively. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

An atlas of human long non-coding RNAs with accurate 5' ends.

Long non-coding RNAs (lncRNAs) are largely heterogeneous and functionally uncharacterized. Here, using FANTOM5 cap analysis of gene expression (CAGE) data, we integrate multiple transcript collections to generate a comprehensive atlas of 27,919 human lncRNA genes with high-confidence 5' ends and expression profiles across 1,829 samples from the major human primary cell types and tissues. Genomic and epigenomic classification of these lncRNAs reveals that most intergenic lncRNAs originate from enhancers rather than from promoters. Incorporating genetic and expression data, we show that lncRNAs overlapping trait-associated single nucleotide polymorphisms are specifically expressed in cell types relevant to the traits, implicating these lncRNAs in multiple diseases. We further demonstrate that lncRNAs overlapping expression quantitative trait loci (eQTL)-associated single nucleotide polymorphisms of messenger RNAs are co-expressed with the corresponding messenger RNAs, suggesting their potential roles in transcriptional regulation. Combining these findings with conservation data, we identify 19,175 potentially functional lncRNAs in the human genome.

Krimper Enforces an Antisense Bias on piRNA Pools by Binding AGO3 in the Drosophila Germline.

Piwi-interacting RNAs (piRNAs) suppress transposon activity in animal germ cells. In the Drosophila ovary, primary Aubergine (Aub)-bound antisense piRNAs initiate the ping-pong cycle to produce secondary AGO3-bound sense piRNAs. This increases the number of secondary Aub-bound antisense piRNAs that can act to destroy transposon mRNAs. Here we show that Krimper (Krimp), a Tudor-domain protein, directly interacts with piRNA-free AGO3 to promote symmetrical dimethylarginine (sDMA) modification, ensuring sense piRNA-loading onto sDMA-modified AGO3. In aub mutant ovaries, AGO3 associates with ping-pong signature piRNAs, suggesting AGO3's compatibility with primary piRNA loading. Krimp sequesters ectopically expressed AGO3 within Krimp bodies in cultured ovarian somatic cells (OSCs), in which only the primary piRNA pathway operates. Upon krimp-RNAi in OSCs, AGO3 loads with piRNAs, further showing the capacity of AGO3 for primary piRNA loading. We propose that Krimp enforces an antisense bias on piRNA pools by binding AGO3 and blocking its access to primary piRNAs.

FANTOM enters 20th year: expansion of transcriptomic atlases and functional annotation of non-coding RNAs.

The Functional ANnoTation Of the Mammalian genome (FANTOM) Consortium has continued to provide extensive resources in the pursuit of understanding the transcriptome, and transcriptional regulation, of mammalian genomes for the last 20 years. To share these resources with the research community, the FANTOM web-interfaces and databases are being regularly updated, enhanced and expanded with new data types. In recent years, the FANTOM Consortium's efforts have been mainly focused on creating new non-coding RNA datasets and resources. The existing FANTOM5 human and mouse miRNA atlas was supplemented with rat, dog, and chicken datasets. The sixth (latest) edition of the FANTOM project was launched to assess the function of human long non-coding RNAs (lncRNAs). From its creation until 2020, FANTOM6 has contributed to the research community a large dataset generated from the knock-down of 285 lncRNAs in human dermal fibroblasts; this is followed with extensive expression profiling and cellular phenotyping. Other updates to the FANTOM resource includes the reprocessing of the miRNA and promoter atlases of human, mouse and chicken with the latest reference genome assemblies. To facilitate the use and accessibility of all above resources we further enhanced FANTOM data viewers and web interfaces. The updated FANTOM web resource is publicly available at https://fantom.gsc.riken.jp/.

Dynamics of cardiomyocyte transcriptome and chromatin landscape demarcates key events of heart development.

Organogenesis involves dynamic regulation of gene transcription and complex multipathway interactions. Despite our knowledge of key factors regulating various steps of heart morphogenesis, considerable challenges in understanding its mechanism still exist because little is known about their downstream targets and interactive regulatory network. To better understand transcriptional regulatory mechanism driving heart development and the consequences of its disruption in vivo, we performed time-series analyses of the transcriptome and genome-wide chromatin accessibility in isolated cardiomyocytes (CMs) from wild-type zebrafish embryos at developmental stages corresponding to heart tube morphogenesis, looping, and maturation. We identified genetic regulatory modules driving crucial events of heart development that contained key cardiac TFs and are associated with open chromatin regions enriched for DNA sequence motifs belonging to the family of the corresponding TFs. Loss of function of cardiac TFs Gata5, Tbx5a, and Hand2 affected the cardiac regulatory networks and caused global changes in chromatin accessibility profile, indicating their role in heart development. Among regions with differential chromatin accessibility in mutants were highly conserved noncoding elements that represent putative enhancers driving heart development. The most prominent gene expression changes, which correlated with chromatin accessibility modifications within their proximal promoter regions, occurred between heart tube morphogenesis and looping, and were associated with metabolic shift and hematopoietic/cardiac fate switch during CM maturation. Our results revealed the dynamic regulatory landscape throughout heart development and identified interactive molecular networks driving key events of heart morphogenesis.