Perspectives on ENCODE

Transcriptional determinism and stochasticity contribute to the complexity of autism-associated SHANK family genes

Precision of transcription is critical because transcriptional dysregulation is disease causing. Traditional methods of transcriptional profiling are inadequate to elucidate the full spectrum of the transcriptome, particularly for longer and less abundant mRNAs. SHANK3 is one of the most common autism causative genes. Twenty-four Shank3-mutant animal lines have been developed for autism modeling. However, their preclinical validity has been questioned due to incomplete Shank3 transcript structure. We apply an integrative approach combining cDNA-capture and long-read sequencing to profile the SHANK3 transcriptome in humans and mice. We unexpectedly discover an extremely complex SHANK3 transcriptome. Specific SHANK3 transcripts are altered in Shank3-mutant mice and postmortem brain tissues from individuals with autism spectrum disorder. The enhanced SHANK3 transcriptome significantly improves the detection rate for potential deleterious variants from genomics studies of neuropsychiatric disorders. Our findings suggest that both deterministic and stochastic transcription of the genome is associated with SHANK family genes.

Big data and deep learning for RNA biology

The exponential growth of big data in RNA biology (RB) has led to the development of deep learning (DL) models that have driven crucial discoveries. As constantly evidenced by DL studies in other fields, the successful implementation of DL in RB depends heavily on the effective utilization of large-scale datasets from public databases. In achieving this goal, data encoding methods, learning algorithms, and techniques that align well with biological domain knowledge have played pivotal roles. In this review, we provide guiding principles for applying these DL concepts to various problems in RB by demonstrating successful examples and associated methodologies. We also discuss the remaining challenges in developing DL models for RB and suggest strategies to overcome these challenges. Overall, this review aims to illuminate the compelling potential of DL for RB and ways to apply this powerful technology to investigate the intriguing biology of RNA more effectively.

Development of DNA Insertion-specific Markers Based on the Intergenic Region of Oplegnathus punctatus Cdkn1/srsf3 for Sex Identification

Spotted knifejaw (Oplegnathus punctatus) is a marine economic fish with high food and ecological value, and its growth process has obvious male and female sexual dimorphism, with males growing significantly faster than females. However, the current sex identification technology is not yet mature, which will limit the growth rate of O. punctatus aquaculture and the efficiency of separate sex breeding, so the development of efficient sex molecular markers is imperative. This study identified a 926 bp DNA insertion fragment in the cdkn1/srsf3 intergenic region of O. punctatus males through whole-genome scanning, comparative genomics, and structural variant analysis. A pair of primers was designed based on the insertion information of the Y chromosome intergenic region in male individuals. Agarose gel electrophoresis revealed the amplification of two DNA fragments, 1118 bp and 192 bp, in male O. punctatus individuals. The 926 bp fragment was identified as the insertion in the intergenic region of cdkn1/srsf3 in males, while only a single 192 bp DNA fragment was amplified in females. The biological sex of the individuals identified in this manner was consistent with their known phenotypic sex. In this study, we developed a method to detect DNA insertion variants in the intergenic region of O. punctatus. Additionally, we introduced a new DNA marker for the rapid identification of the sex of O. punctatus, which enhances detection efficiency. The text has important reference significance and application value in sex identification, all-male breeding, and lineage selection. It provides new insights into the regulation of variation in the intergenic region of cdkn1/srsf3 genes and the study of RNA shearing.

An oncogenic enhancer promotes melanoma progression via regulating ETV4 expression

BACKGROUND

Enhancers are important gene regulatory elements that promote the expression of critical genes in development and disease. Aberrant enhancer can modulate cancer risk and activate oncogenes that lead to the occurrence of various cancers. However, the underlying mechanism of most enhancers in cancer remains unclear. Here, we aim to explore the function and mechanism of a crucial enhancer in melanoma.

METHODS

Multi-omics data were applied to identify an enhancer (enh17) involved in melanoma progression. To evaluate the function of enh17, CRISPR/Cas9 technology were applied to knockout enh17 in melanoma cell line A375. RNA-seq, ChIP-seq and Hi-C data analysis integrated with luciferase reporter assay were performed to identify the potential target gene of enh17. Functional experiments were conducted to further validate the function of the target gene ETV4. Multi-omics data integrated with CUT&Tag sequencing were performed to validate the binding profile of the inferred transcription factor STAT3.

RESULTS

An enhancer, named enh17 here, was found to be aberrantly activated and involved in melanoma progression. CRISPR/Cas9-mediated deletion of enh17 inhibited cell proliferation, migration, and tumor growth of melanoma both in vitro and in vivo. Mechanistically, we identified ETV4 as a target gene regulated by enh17, and functional experiments further support ETV4 as a target gene that is involved in cancer-associated phenotypes. In addition, STAT3 acts as a transcription factor binding with enh17 to regulate the transcription of ETV4.

CONCLUSIONS

Our findings revealed that enh17 plays an oncogenic role and promotes tumor progression in melanoma, and its transcriptional regulatory mechanisms were fully elucidated, which may open a promising window for melanoma prevention and treatment.

Multiomic profiling of transcription factor binding and function in human brain

Transcription factors (TFs) orchestrate gene expression programs crucial for brain function, but we lack detailed information about TF binding in human brain tissue. We generated a multiomic resource (ChIP-seq, ATAC-seq, RNA-seq, DNA methylation) on bulk tissues and sorted nuclei from several postmortem brain regions, including binding maps for more than 100 TFs. We demonstrate improved measurements of TF activity, including motif recognition and gene expression modeling, upon identification and removal of high TF occupancy regions. Further, predictive TF binding models demonstrate a bias for these high-occupancy sites. Neuronal TFs SATB2 and TBR1 bind unique regions depleted for such sites and promote neuronal gene expression. Binding sites for TFs, including TBR1 and PKNOX1, are enriched for risk variants associated with neuropsychiatric disorders, predominantly in neurons. This work, titled BrainTF, is a powerful resource for future studies seeking to understand the roles of specific TFs in regulating gene expression in the human brain.

Genome-wide Cas9-mediated screening of essential non-coding regulatory elements via libraries of paired single-guide RNAs

The functions of non-coding regulatory elements (NCREs), which constitute a major fraction of the human genome, have not been systematically studied. Here we report a method involving libraries of paired single-guide RNAs targeting both ends of an NCRE as a screening system for the Cas9-mediated deletion of thousands of NCREs genome-wide to study their functions in distinct biological contexts. By using K562 and 293T cell lines and human embryonic stem cells, we show that NCREs can have redundant functions, and that many ultra-conserved elements have silencer activity and play essential roles in cell growth and in cellular responses to drugs (notably, the ultra-conserved element PAX6_Tarzan may be critical for heart development, as removing it from human embryonic stem cells led to defects in cardiomyocyte differentiation). The high-throughput screen, which is compatible with single-cell sequencing, may allow for the identification of druggable NCREs.

Deconvolution of polygenic risk score in single cells unravels cellular and molecular heterogeneity of complex human diseases

Polygenic risk scores (PRSs) are commonly used for predicting an individual's genetic risk of complex diseases. Yet, their implication for disease pathogenesis remains largely limited. Here, we introduce scPRS, a geometric deep learning model that constructs single-cell-resolved PRS leveraging reference single-cell chromatin accessibility profiling data to enhance biological discovery as well as disease prediction. Real-world applications across multiple complex diseases, including type 2 diabetes (T2D), hypertrophic cardiomyopathy (HCM), and Alzheimer's disease (AD), showcase the superior prediction power of scPRS compared to traditional PRS methods. Importantly, scPRS not only predicts disease risk but also uncovers disease-relevant cells, such as hormone-high alpha and beta cells for T2D, cardiomyocytes and pericytes for HCM, and astrocytes, microglia and oligodendrocyte progenitor cells for AD. Facilitated by a layered multi-omic analysis, scPRS further identifies cell-type-specific genetic underpinnings, linking disease-associated genetic variants to gene regulation within corresponding cell types. We substantiate the disease relevance of scPRS-prioritized HCM genes and demonstrate that the suppression of these genes in HCM cardiomyocytes is rescued by Mavacamten treatment. Additionally, we establish a novel microglia-specific regulatory relationship between the AD risk variant rs7922621 and its target genes ANXA11 and TSPAN14. We further illustrate the detrimental effects of suppressing these two genes on microglia phagocytosis. Our work provides a multi-tasking, interpretable framework for precise disease prediction and systematic investigation of the genetic, cellular, and molecular basis of complex diseases, laying the methodological foundation for single-cell genetics.

Are non-protein coding RNAs junk or treasure?: An attempt to explain and reconcile opposing viewpoints of whether the human genome is mostly transcribed into non-functional or functional RNAs

The human genome project's lasting legacies are the emerging insights into human physiology and disease, and the ascendance of biology as the dominant science of the 21st century. Sequencing revealed that >90% of the human genome is not coding for proteins, as originally thought, but rather is overwhelmingly transcribed into non-protein coding, or non-coding, RNAs (ncRNAs). This discovery initially led to the hypothesis that most genomic DNA is "junk", a term still championed by some geneticists and evolutionary biologists. In contrast, molecular biologists and biochemists studying the vast number of transcripts produced from most of this genome "junk" often surmise that these ncRNAs have biological significance. What gives? This essay contrasts the two opposing, extant viewpoints, aiming to explain their bases, which arise from distinct reference frames of the underlying scientific disciplines. Finally, it aims to reconcile these divergent mindsets in hopes of stimulating synergy between scientific fields.

A systematic analyses of different bioinformatics pipelines for genomic data and its impact on deep learning models for chromatin loop prediction

Genomic data analysis has witnessed a surge in complexity and volume, primarily driven by the advent of high-throughput technologies. In particular, studying chromatin loops and structures has become pivotal in understanding gene regulation and genome organization. This systematic investigation explores the realm of specialized bioinformatics pipelines designed specifically for the analysis of chromatin loops and structures. Our investigation incorporates two protein (CTCF and Cohesin) factor-specific loop interaction datasets from six distinct pipelines, amassing a comprehensive collection of 36 diverse datasets. Through a meticulous review of existing literature, we offer a holistic perspective on the methodologies, tools and algorithms underpinning the analysis of this multifaceted genomic feature. We illuminate the vast array of approaches deployed, encompassing pivotal aspects such as data preparation pipeline, preprocessing, statistical features and modelling techniques. Beyond this, we rigorously assess the strengths and limitations inherent in these bioinformatics pipelines, shedding light on the interplay between data quality and the performance of deep learning models, ultimately advancing our comprehension of genomic intricacies.

Transcriptional Determinism and Stochasticity Contribute to the Complexity of Autism Associated SHANK Family Genes

Precision of transcription is critical because transcriptional dysregulation is disease causing. Traditional methods of transcriptional profiling are inadequate to elucidate the full spectrum of the transcriptome, particularly for longer and less abundant mRNAs. SHANK3 is one of the most common autism causative genes. Twenty-four Shank3 mutant animal lines have been developed for autism modeling. However, their preclinical validity has been questioned due to incomplete Shank3 transcript structure. We applied an integrative approach combining cDNA-capture and long-read sequencing to profile the SHANK3 transcriptome in human and mice. We unexpectedly discovered an extremely complex SHANK3 transcriptome. Specific SHANK3 transcripts were altered in Shank3 mutant mice and postmortem brains tissues from individuals with ASD. The enhanced SHANK3 transcriptome significantly improved the detection rate for potential deleterious variants from genomics studies of neuropsychiatric disorders. Our findings suggest the stochastic transcription of genome associated with SHANK family genes.

Epigenetic regulatory layers in the 3D nucleus

Nearly 7 decades have elapsed since Francis Crick introduced the central dogma of molecular biology, as part of his ideas on protein synthesis, setting the fundamental rules of sequence information transfer from DNA to RNAs and proteins. We have since learned that gene expression is finely tuned in time and space, due to the activities of RNAs and proteins on regulatory DNA elements, and through cell-type-specific three-dimensional conformations of the genome. Here, we review major advances in genome biology and discuss a set of ideas on gene regulation and highlight how various biomolecular assemblies lead to the formation of structural and regulatory features within the nucleus, with roles in transcriptional control. We conclude by suggesting further developments that will help capture the complex, dynamic, and often spatially restricted events that govern gene expression in mammalian cells.

HLA-B*27 and Ankylosing Spondylitis: 50 Years of Insights and Discoveries

PURPOSE OF REVIEW

To commemorate the 50th anniversary of the groundbreaking discovery of a remarkably strong association between HLA-B*27 and ankylosing spondylitis (AS).

RECENT FINDINGS

In addition to HLA-B*27, more than 116 other recognized genetic risk variants have been identified, while epigenetic factors largely remain unexplored in this context. Among patients with AS who carry the HLA-B*27 gene, clonally expanded CD8 + T cells can be found in their bloodstream and within inflamed tissues. Moreover, the α and β chain motifs of these T-cell receptors demonstrate a distinct affinity for certain self- and microbial-derived peptides, leading to an autoimmune response that ultimately results in the onset of the disease. These distinctive peptide-binding and presentation characteristics are a hallmark of the disease-associated HLA-B*27:05 subtype but are absent in HLA-B*27:09, a subtype not associated with the disease, differing by only a single amino acid. This discovery represents a significant advancement in unraveling the 50-year-old puzzle of how HLA-B*27 contributes to the development of AS. These findings will significantly accelerate the process of identifying peptides, both self- and microbial-derived, that instigate autoimmunity. This, in return, will pave the way for the development of more accurate and effective targeted treatments. Moreover, the discovery of improved biomarkers, in conjunction with the emerging technology of electric field molecular fingerprinting, has the potential to greatly bolster early diagnosis capabilities. A very recently published groundbreak paper underscores the remarkable effectiveness of targeting and eliminating disease-causing T cells in a HLA-B*27 patients with AS. This pivotal advancement not only signifies a paradigm shift but also bolsters the potential for preventing the disease in individuals carrying high-risk genetic variants.

RIBFIND2: Identifying rigid bodies in protein and nucleic acid structures

Molecular structures are often fitted into cryo-EM maps by flexible fitting. When this requires large conformational changes, identifying rigid bodies can help optimize the model-map fit. Tools for identifying rigid bodies in protein structures exist, however an equivalent for nucleic acid structures is lacking. With the increase in cryo-EM maps containing RNA and progress in RNA structure prediction, there is a need for such tools. We previously developed RIBFIND, a program for clustering protein secondary structures into rigid bodies. In RIBFIND2, this approach is extended to nucleic acid structures. RIBFIND2 can identify biologically relevant rigid bodies in important groups of complex RNA structures, capturing a wide range of dynamics, including large rigid-body movements. The usefulness of RIBFIND2-assigned rigid bodies in cryo-EM model refinement was demonstrated on three examples, with two conformations each: Group II Intron complexed IEP, Internal Ribosome Entry Site and the Processome, using cryo-EM maps at 2.7-5 Å resolution. A hierarchical refinement approach, performed on progressively smaller sets of RIBFIND2 rigid bodies, was clearly shown to have an advantage over classical all-atom refinement. RIBFIND2 is available via a web server with structure visualization and as a standalone tool.

Noncoding RNAs and Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes in Cardiac Arrhythmic Brugada Syndrome

Hundreds of thousands of people die each year as a result of sudden cardiac death, and many are due to heart rhythm disorders. One of the major causes of these arrhythmic events is Brugada syndrome, a cardiac channelopathy that results in abnormal cardiac conduction, severe life-threatening arrhythmias, and, on many occasions, death. This disorder has been associated with mutations and dysfunction of about two dozen genes; however, the majority of the patients do not have a definite cause for the diagnosis of Brugada Syndrome. The protein-coding genes represent only a very small fraction of the mammalian genome, and the majority of the noncoding regions of the genome are actively transcribed. Studies have shown that most of the loci associated with electrophysiological traits are located in noncoding regulatory regions and are expected to affect gene expression dosage and cardiac ion channel function. Noncoding RNAs serve an expanding number of regulatory and other functional roles within the cells, including but not limited to transcriptional, post-transcriptional, and epigenetic regulation. The major noncoding RNAs found in Brugada Syndrome include microRNAs; however, others such as long noncoding RNAs are also identified. They contribute to pathogenesis by interacting with ion channels and/or are detectable as clinical biomarkers. Stem cells have received significant attention in the recent past, and can be differentiated into many different cell types including those in the heart. In addition to contractile and relaxational properties, BrS-relevant electrophysiological phenotypes are also demonstrated in cardiomyocytes differentiated from stem cells induced from adult human cells. In this review, we discuss the current understanding of noncoding regions of the genome and their RNA biology in Brugada Syndrome. We also delve into the role of stem cells, especially human induced pluripotent stem cell-derived cardiac differentiated cells, in the investigation of Brugada syndrome in preclinical and clinical studies.

Interpreting non-coding disease-associated human variants using single-cell epigenomics

Genome-wide association studies (GWAS) have linked hundreds of thousands of sequence variants in the human genome to common traits and diseases. However, translating this knowledge into a mechanistic understanding of disease-relevant biology remains challenging, largely because such variants are predominantly in non-protein-coding sequences that still lack functional annotation at cell-type resolution. Recent advances in single-cell epigenomics assays have enabled the generation of cell type-, subtype- and state-resolved maps of the epigenome in heterogeneous human tissues. These maps have facilitated cell type-specific annotation of candidate cis-regulatory elements and their gene targets in the human genome, enhancing our ability to interpret the genetic basis of common traits and diseases.

Bioinformatics for wet-lab scientists: practical application in sequencing analysis

BACKGROUND

Genomics data is available to the scientific community after publication of research projects and can be investigated for a multitude of research questions. However, in many cases deposited data is only assessed and used for the initial publication, resulting in valuable resources not being exploited to their full depth. MAIN: A likely reason for this is that many wetlab-based researchers are not formally trained to apply bioinformatic tools and may therefore assume that they lack the necessary experience to do so themselves. In this article, we present a series of freely available, predominantly web-based platforms and bioinformatic tools that can be combined in analysis pipelines to interrogate different types of next-generation sequencing data. Additionally to the presented exemplary route, we also list a number of alternative tools that can be combined in a mix-and-match fashion. We place special emphasis on tools that can be followed and used correctly without extensive prior knowledge in programming. Such analysis pipelines can be applied to existing data downloaded from the public domain or be compared to the results of own experiments.

CONCLUSION

Integrating transcription factor binding to chromatin (ChIP-seq) with transcriptional output (RNA-seq) and chromatin accessibility (ATAC-seq) can not only assist to form a deeper understanding of the molecular interactions underlying transcriptional regulation but will also help establishing new hypotheses and pre-testing them in silico.

Cooperation between bHLH transcription factors and histones for DNA access

The basic helix-loop-helix (bHLH) family of transcription factors recognizes DNA motifs known as E-boxes (CANNTG) and includes 108 members1. Here we investigate how chromatinized E-boxes are engaged by two structurally diverse bHLH proteins: the proto-oncogene MYC-MAX and the circadian transcription factor CLOCK-BMAL1 (refs. 2,3). Both transcription factors bind to E-boxes preferentially near the nucleosomal entry-exit sites. Structural studies with engineered or native nucleosome sequences show that MYC-MAX or CLOCK-BMAL1 triggers the release of DNA from histones to gain access. Atop the H2A-H2B acidic patch4, the CLOCK-BMAL1 Per-Arnt-Sim (PAS) dimerization domains engage the histone octamer disc. Binding of tandem E-boxes5-7 at endogenous DNA sequences occurs through direct interactions between two CLOCK-BMAL1 protomers and histones and is important for circadian cycling. At internal E-boxes, the MYC-MAX leucine zipper can also interact with histones H2B and H3, and its binding is indirectly enhanced by OCT4 elsewhere on the nucleosome. The nucleosomal E-box position and the type of bHLH dimerization domain jointly determine the histone contact, the affinity and the degree of competition and cooperativity with other nucleosome-bound factors.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

A modERN resource: identification of Drosophila transcription factor candidate target genes using RNAi

Transcription factors (TFs) play a key role in development and in cellular responses to the environment by activating or repressing the transcription of target genes in precise spatial and temporal patterns. In order to develop a catalog of target genes of Drosophila melanogaster TFs, the modERN consortium systematically knocked down the expression of TFs using RNAi in whole embryos followed by RNA-seq. We generated data for 45 TFs which have 18 different DNA-binding domains and are expressed in 15 of the 16 organ systems. The range of inactivation of the targeted TFs by RNAi ranged from log2fold change -3.52 to +0.49. The TFs also showed remarkable heterogeneity in the numbers of candidate target genes identified, with some generating thousands of candidates and others only tens. We present detailed analysis from five experiments, including those for three TFs that have been the focus of previous functional studies (ERR, sens, and zfh2) and two previously uncharacterized TFs (sens-2 and CG32006), as well as short vignettes for selected additional experiments to illustrate the utility of this resource. The RNA-seq datasets are available through the ENCODE DCC (http://encodeproject.org) and the Sequence Read Archive (SRA). TF and target gene expression patterns can be found here: https://insitu.fruitfly.org. These studies provide data that facilitate scientific inquiries into the functions of individual TFs in key developmental, metabolic, defensive, and homeostatic regulatory pathways, as well as provide a broader perspective on how individual TFs work together in local networks during embryogenesis.

Leveraging Base Pair Mammalian Constraint to Understand Genetic Variation and Human Disease

Although thousands of genomic regions have been associated with heritable human diseases, attempts to elucidate biological mechanisms are impeded by a general inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function that is agnostic to cell type or disease mechanism. Here, single base phyloP scores from the whole genome alignment of 240 placental mammals identified 3.5% of the human genome as significantly constrained, and likely functional. We compared these scores to large-scale genome annotation, genome-wide association studies (GWAS), copy number variation, clinical genetics findings, and cancer data sets. Evolutionarily constrained positions are enriched for variants explaining common disease heritability (more than any other functional annotation). Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Leveraging transcriptomics for precision diagnosis: Lessons learned from cancer and sepsis

Diagnostics require precision and predictive ability to be clinically useful. Integration of multi-omic with clinical data is crucial to our understanding of disease pathogenesis and diagnosis. However, interpretation of overwhelming amounts of information at the individual level requires sophisticated computational tools for extraction of clinically meaningful outputs. Moreover, evolution of technical and analytical methods often outpaces standardisation strategies. RNA is the most dynamic component of all -omics technologies carrying an abundance of regulatory information that is least harnessed for use in clinical diagnostics. Gene expression-based tests capture genetic and non-genetic heterogeneity and have been implemented in certain diseases. For example patients with early breast cancer are spared toxic unnecessary treatments with scores based on the expression of a set of genes (e.g., Oncotype DX). The ability of transcriptomics to portray the transcriptional status at a moment in time has also been used in diagnosis of dynamic diseases such as sepsis. Gene expression profiles identify endotypes in sepsis patients with prognostic value and a potential to discriminate between viral and bacterial infection. The application of transcriptomics for patient stratification in clinical environments and clinical trials thus holds promise. In this review, we discuss the current clinical application in the fields of cancer and infection. We use these paradigms to highlight the impediments in identifying useful diagnostic and prognostic biomarkers and propose approaches to overcome them and aid efforts towards clinical implementation.

[The birth of a gene]

The availability of an extensive set of vertebrate genome sequences, together with large-scale transcriptome studies, has allowed the discovery of numerous non-canonical ORFs (usually quite short) with evidence of transcription, translation and functional involvement. Orthologs for these ORFs can be detected in many vertebrates, and the time of appearance of a functional mini-gene can be ascertained. Some of these have appeared quite recently in evolution and have already very specific expression patterns in humans.

Lack of GPR180 ameliorates hepatic lipid depot via downregulation of mTORC1 signaling

Our previous genome-wide association study to explore genetic loci associated with lean nonalcoholic fatty liver disease (NAFLD) in Japan suggested four candidate loci, which were mapped to chr6, chr7, chr12 and chr13. The present study aimed to identify the locus involved functionally in NAFLD around the association signal observed in chr13. Chromosome conformation capture assay and a database survey suggested the intermolecular interaction among DNA fragments in association signals with the adjacent four coding gene promoters. The four genes were further screened by knockdown (KD) in mice using shRNA delivered by an adeno-associated virus vector (AAV8), and KD of G protein-coupled receptor 180 (Gpr180) showed amelioration of hepatic lipid storage. Gpr180 knockout (KO) mice also showed ameliorated hepatic and plasma lipid levels without influencing glucose metabolism after high-fat diet intake. Transcriptome analyses showed downregulation of mTORC1 signaling and cholesterol homeostasis, which was confirmed by weakened phosphorylation of mTOR and decreased activated SREBP1 in Gpr180KO mice and a human hepatoma cell line (Huh7). AAV8-mediated hepatic rescue of GPR180 expression in KO mice showed recovery of plasma and hepatic lipid levels. In conclusion, ablation of GPR180 ameliorated plasma and hepatic lipid levels, which was mediated by downregulation of mTORC1 signaling.

ATM suppresses c-Myc overexpression in the mammary epithelium in response to estrogen

ATM gene mutation carriers are predisposed to estrogen-receptor-positive breast cancer (BC). ATM prevents BC oncogenesis by activating p53 in every cell; however, much remains unknown about tissue-specific oncogenesis after ATM loss. Here, we report that ATM controls the early transcriptional response to estrogens. This response depends on topoisomerase II (TOP2), which generates TOP2-DNA double-strand break (DSB) complexes and rejoins the breaks. When TOP2-mediated ligation fails, ATM facilitates DSB repair. After estrogen exposure, TOP2-dependent DSBs arise at the c-MYC enhancer in human BC cells, and their defective repair changes the activation profile of enhancers and induces the overexpression of many genes, including the c-MYC oncogene. CRISPR/Cas9 cleavage at the enhancer also causes c-MYC overexpression, indicating that this DSB causes c-MYC overexpression. Estrogen treatment induced c-Myc protein overexpression in mammary epithelial cells of ATM-deficient mice. In conclusion, ATM suppresses the c-Myc-driven proliferative effects of estrogens, possibly explaining such tissue-specific oncogenesis.

Chicken chromatin accessibility atlas accelerates epigenetic annotation of birds and gene fine-mapping associated with growth traits

The development of epigenetic maps, such as the ENCODE project in humans, provides resources for gene regulation studies and a reference for research of disease-related regulatory elements. However, epigenetic information, such as a bird-specific chromatin accessibility atlas, is currently lacking for the thousands of bird species currently described. The major genomic difference between birds and mammals is their shorter introns and intergenic distances, which seriously hinders the use of humans and mice as a reference for studying the function of important regulatory regions in birds. In this study, using chicken as a model bird species, we systematically compiled a chicken chromatin accessibility atlas using 53 Assay of Transposase Accessible Chromatin sequencing (ATAC-seq) samples across 11 tissues. An average of 50 796 open chromatin regions were identified per sample, cumulatively accounting for 20.36% of the chicken genome. Tissue specificity was largely reflected by differences in intergenic and intronic peaks, with specific functional regulation achieved by two mechanisms: recruitment of several sequence-specific transcription factors and direct regulation of adjacent functional genes. By integrating data from genome-wide association studies, our results suggest that chicken body weight is driven by different regulatory variants active in growth-relevant tissues. We propose CAB39L (active in the duodenum), RCBTB1 (muscle and liver), and novel long non-coding RNA ENSGALG00000053256 (bone) as candidate genes regulating chicken body weight. Overall, this study demonstrates the value of epigenetic data in fine-mapping functional variants and provides a compendium of resources for further research on the epigenetics and evolution of birds and mammals.

Antagonistic action of GPS2 and KDM1A at enhancers governs alternative macrophage activation by interleukin 4

The Th2 cytokine interleukin 4 (IL4) promotes macrophage differentiation into alternative subtypes and plays important roles in physiology, in metabolic and inflammatory diseases, in cancer and in tissue regeneration. While the regulatory transcription factor networks governing IL4 signaling are already well-characterized, it is currently less understood which transcriptional coregulators are involved and how they operate mechanistically. In this study, we discover that G protein pathway suppressor 2 (GPS2), a core subunit of the HDAC3 corepressor complex assembled by SMRT and NCOR, represses IL4-dependent enhancer activation in mouse macrophages. Our genome-wide and gene-specific characterization revealed that, instead of directly repressing STAT6, chromatin-bound GPS2 cooperates with SMRT and NCOR to antagonize enhancer activation by lysine demethylase 1A (KDM1A, LSD1). Mechanistically, corepressor depletion increased KDM1A recruitment to enhancers linked to IL4-induced genes, accompanied by demethylation of the repressive histone marks H3K9me2/3 without affecting H3K4me1/2, the classic KDM1A substrates for demethylation in other cellular contexts. This in turn caused enhancer and gene activation already in the absence of IL4/STAT6 and sensitized the STAT6-dependent IL4 responsiveness of macrophages. Thus, our work identified with the antagonistic action of a GPS2-containing corepressor complex and the lysine demethylase KDM1A a hitherto unknown epigenetic corepressor-coactivator switching mechanism that governs alternative macrophage activation.

The epigenetic regulation of cancer cell recovery from therapy exposure and its implications as a novel therapeutic strategy for preventing disease recurrence

The ultimate goal of cancer therapy is the elimination of disease from patients. Most directly, this occurs through therapy-induced cell death. Therapy-induced growth arrest can also be a desirable outcome, if prolonged. Unfortunately, therapy-induced growth arrest is rarely durable and the recovering cell population can contribute to cancer recurrence. Consequently, therapeutic strategies that eliminate residual cancer cells reduce opportunities for recurrence. Recovery can occur through diverse mechanisms including quiescence or diapause, exit from senescence, suppression of apoptosis, cytoprotective autophagy, and reductive divisions resulting from polyploidy. Epigenetic regulation of the genome represents a fundamental regulatory mechanism integral to cancer-specific biology, including the recovery from therapy. Epigenetic pathways are particularly attractive therapeutic targets because they are reversible, without changes in DNA, and are catalyzed by druggable enzymes. Previous use of epigenetic-targeting therapies in combination with cancer therapeutics has not been widely successful because of either unacceptable toxicity or limited efficacy. The use of epigenetic-targeting therapies after a significant interval following initial cancer therapy could potentially reduce the toxicity of combination strategies, and possibly exploit essential epigenetic states following therapy exposure. This review examines the feasibility of targeting epigenetic mechanisms using a sequential approach to eliminate residual therapy-arrested populations, that might possibly prevent recovery and disease recurrence.

Heritability enrichment in context-specific regulatory networks improves phenotype-relevant tissue identification

Systems genetics holds the promise to decipher complex traits by interpreting their associated SNPs through gene regulatory networks derived from comprehensive multi-omics data of cell types, tissues, and organs. Here, we propose SpecVar to integrate paired chromatin accessibility and gene expression data into context-specific regulatory network atlas and regulatory categories, conduct heritability enrichment analysis with genome-wide association studies (GWAS) summary statistics, identify relevant tissues, and estimate relevance correlation to depict common genetic factors acting in the shared regulatory networks between traits. Our method improves power upon existing approaches by associating SNPs with context-specific regulatory elements to assess heritability enrichments and by explicitly prioritizing gene regulations underlying relevant tissues. Ablation studies, independent data validation, and comparison experiments with existing methods on GWAS of six phenotypes show that SpecVar can improve heritability enrichment, accurately detect relevant tissues, and reveal causal regulations. Furthermore, SpecVar correlates the relevance patterns for pairs of phenotypes and better reveals shared SNP-associated regulations of phenotypes than existing methods. Studying GWAS of 206 phenotypes in UK Biobank demonstrates that SpecVar leverages the context-specific regulatory network atlas to prioritize phenotypes' relevant tissues and shared heritability for biological and therapeutic insights. SpecVar provides a powerful way to interpret SNPs via context-specific regulatory networks and is available at https://github.com/AMSSwanglab/SpecVar, copy archived at swh:1:rev:cf27438d3f8245c34c357ec5f077528e6befe829.

Chromatin structure in cancer

In the past decade, we have seen the emergence of sequence-based methods to understand chromosome organization. With the confluence of in situ approaches to capture information on looping, topological domains, and larger chromatin compartments, understanding chromatin-driven disease is becoming feasible. Excitingly, recent advances in single molecule imaging with capacity to reconstruct “bulk-cell” features of chromosome conformation have revealed cell-to-cell chromatin structural variation. The fundamental question motivating our analysis of the literature is, can altered chromatin structure drive tumorigenesis? As our community learns more about rare disease, including low mutational frequency cancers, understanding “chromatin-driven” pathology will illuminate the regulatory structures of the genome. We describe recent insights into altered genome architecture in human cancer, highlighting multiple pathways toward disruptions of chromatin structure, including structural variation, noncoding mutations, metabolism, and de novo mutations to architectural regulators themselves. Our analysis of the literature reveals that deregulation of genome structure is characteristic in distinct classes of chromatin-driven tumors. As we begin to integrate the findings from single cell imaging studies and chromatin structural sequencing, we will be able to understand the diversity of cells within a common diagnosis, and begin to define structure–function relationships of the misfolded genome.

An efficient five-lncRNA signature for lung adenocarcinoma prognosis, with AL606489.1 showing sexual dimorphism

Background: Lung adenocarcinoma (LUAD) is a sex-biased and easily metastatic malignant disease. A signature based on 5 long non-coding RNAs (lncRNAs) has been established to promote the overall survival (OS) prediction effect on LUAD.

Methods: The RNA expression profiles of LUAD patients were obtained from The Cancer Genome Atlas. OS-associated lncRNAs were identified based on the differential expression analysis between LUAD and normal samples followed by survival analysis, univariate and multivariate Cox proportional hazards regression analyses. OS-associated lncRNA with sex dimorphism was determined based on the analysis of expression between males and females. Functional enrichment analysis of the Gene Ontology (GO) terms and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways was performed to explore the possible mechanisms of 5-lncRNA signatures.

Results: A 5-lncRNA signature (composed of AC068228.1, SATB2-AS1, LINC01843, AC026355.1, and AL606489.1) was found to be effective in predicting high-risk LUAD patients as well as applicable to female and male subgroups and <65-year and ≥65-year age subgroups. The forecasted effect of the 5-lncRNA signature was more efficient and stable than the TNM stage and other clinical risk factors (such as sex and age). Functional enrichment analysis revealed that the mRNA co-expressed with these five OS-related lncRNAs was associated with RNA regulation within the nucleus. AL606489.1 demonstrated a sexual dimorphism that may be associated with microtubule activity.

Conclusion: Our 5-lncRNA signature could efficaciously predict the OS of LUAD patients. AL606489.1 demonstrated gender dimorphism, which provides a new direction for mechanistic studies on sexual dimorphism.

Unravelling the impact of aging on the human endothelial lncRNA transcriptome

The incidence and prevalence of cardiovascular disease is highest among the elderly. There is a need to further understand the mechanisms behind endothelial cell aging in order to achieve vascular rejuvenation and minimize the onset of age-related vascular diseases. Long non-coding RNAs (lncRNAs) have been proposed to regulate numerous processes in the human genome, yet their function in vascular aging and their therapeutic potential remain largely unknown. This is primarily because the majority of studies investigating the impact of aging on lncRNA expression heavily rely on in vitro studies based on replicative senescence. Here, using a unique collection of young and aged endothelial cells isolated from native human arteries, we sought to characterize the age-related alterations in lncRNA expression profiles. We were able to detect a total of 4463 lncRNAs expressed in the human endothelium from which ∼17% (798) were altered in advanced age. One of the most affected lncRNAs in aging was the primate-specific, Prostate Cancer Associated Transcript (PCAT) 14. In our follow up analysis, using single molecule RNA FISH, we showed that PCAT14 is relatively abundant, localized almost exclusively in the nucleus of young endothelial cells, and silenced in the aged endothelium. Functionally, our studies proposed that downregulation of PCAT14 alters endothelial cell transcription profile and cell functions including endothelial cell migration, sprouting and inflammatory responses in vitro. Taken together, our data highlight that endothelial cell aging correlates with altered expression of lncRNAs, which could impair the endothelial regenerative capacity and enhance inflammatory phenotypes.

From bugs to bedside: functional annotation of human genetic variation for neurological disorders using invertebrate models

The exponential accumulation of DNA sequencing data has opened new avenues for discovering the causative roles of single-nucleotide polymorphisms (SNPs) in neurological diseases. The opportunities emerging from this are staggering, yet only as good as our abilities to glean insights from this surplus of information. Whereas computational biology continues to improve with respect to predictions and molecular modeling, the differences between in silico and in vivo analysis remain substantial. Invertebrate in vivo model systems represent technically advanced, experimentally mature, high-throughput, efficient and cost-effective resources for investigating a disease. With a decades-long track record of enabling investigators to discern function from DNA, fly (Drosophila) and worm (Caenorhabditis elegans) models have never been better poised to serve as living engines of discovery. Both of these animals have already proven useful in the classification of genetic variants as either pathogenic or benign across a range of neurodevelopmental and neurodegenerative disorders-including autism spectrum disorders, ciliopathies, amyotrophic lateral sclerosis, Alzheimer's and Parkinson's disease. Pathogenic SNPs typically display distinctive phenotypes in functional assays when compared with null alleles and frequently lead to protein products with gain-of-function or partial loss-of-function properties that contribute to neurological disease pathogenesis. The utility of invertebrates is logically limited by overt differences in anatomical and physiological characteristics, and also the evolutionary distance in genome structure. Nevertheless, functional annotation of disease-SNPs using invertebrate models can expedite the process of assigning cellular and organismal consequences to mutations, ascertain insights into mechanisms of action, and accelerate therapeutic target discovery and drug development for neurological conditions.

The Genomics of Auditory Function and Disease

Current estimates suggest that nearly half a billion people worldwide are affected by hearing loss. Because of the major psychological, social, economic, and health ramifications, considerable efforts have been invested in identifying the genes and molecular pathways involved in hearing loss, whether genetic or environmental, to promote prevention, improve rehabilitation, and develop therapeutics. Genomic sequencing technologies have led to the discovery of genes associated with hearing loss. Studies of the transcriptome and epigenome of the inner ear have characterized key regulators and pathways involved in the development of the inner ear and have paved the way for their use in regenerative medicine. In parallel, the immense preclinical success of using viral vectors for gene delivery in animal models of hearing loss has motivated the industry to work on translating such approaches into the clinic. Here, we review the recent advances in the genomics of auditory function and dysfunction, from patient diagnostics to epigenetics and gene therapy.

Sodium-calcium exchanger-3 regulates pain "wind-up": From human psychophysics to spinal mechanisms

Repeated application of noxious stimuli leads to a progressively increased pain perception; this temporal summation is enhanced in and predictive of clinical pain disorders. Its electrophysiological correlate is "wind-up," in which dorsal horn spinal neurons increase their response to repeated nociceptor stimulation. To understand the genetic basis of temporal summation, we undertook a GWAS of wind-up in healthy human volunteers and found significant association with SLC8A3 encoding sodium-calcium exchanger type 3 (NCX3). NCX3 was expressed in mouse dorsal horn neurons, and mice lacking NCX3 showed normal, acute pain but hypersensitivity to the second phase of the formalin test and chronic constriction injury. Dorsal horn neurons lacking NCX3 showed increased intracellular calcium following repetitive stimulation, slowed calcium clearance, and increased wind-up. Moreover, virally mediated enhanced spinal expression of NCX3 reduced central sensitization. Our study highlights Ca2+ efflux as a pathway underlying temporal summation and persistent pain, which may be amenable to therapeutic targeting.

Repetitive DNA symmetry elements negatively regulate gene expression in embryonic stem cells

Transcription factor (TF) binding to genomic DNA elements constitutes one of the key mechanisms that regulates gene expression program in cells. Both consensus and nonconsensus DNA sequence elements influence the recognition specificity of TFs. Based on the analysis of experimentally determined c-Myc binding preferences to genomic DNA, here we statistically predict that certain repetitive, nonconsensus DNA symmetry elements can relatively reduce TF-DNA binding preferences. This is in contrast to a different set of repetitive, nonconsensus symmetry elements that can increase the strength of TF-DNA binding. Using c-Myc enhancer reporter system containing consensus motif flanked by nonconsensus sequences in embryonic stem cells, we directly demonstrate that the enrichment in such negatively regulating repetitive symmetry elements is sufficient to reduce the gene expression level compared with native genomic sequences. Negatively regulating repetitive symmetry elements around consensus c-Myc motif and DNA sequences containing consensus c-Myc motif flanked by entirely randomized sequences show similar expression baseline. A possible explanation for this observation is that rather than complete repression, negatively regulating repetitive symmetry elements play a regulatory role in fine-tuning the reduction of gene expression, most probably by binding TFs other than c-Myc.

Pan-African genome demonstrates how population-specific genome graphs improve high-throughput sequencing data analysis

Graph-based genome reference representations have seen significant development, motivated by the inadequacy of the current human genome reference to represent the diverse genetic information from different human populations and its inability to maintain the same level of accuracy for non-European ancestries. While there have been many efforts to develop computationally efficient graph-based toolkits for NGS read alignment and variant calling, methods to curate genomic variants and subsequently construct genome graphs remain an understudied problem that inevitably determines the effectiveness of the overall bioinformatics pipeline. In this study, we discuss obstacles encountered during graph construction and propose methods for sample selection based on population diversity, graph augmentation with structural variants and resolution of graph reference ambiguity caused by information overload. Moreover, we present the case for iteratively augmenting tailored genome graphs for targeted populations and demonstrate this approach on the whole-genome samples of African ancestry. Our results show that population-specific graphs, as more representative alternatives to linear or generic graph references, can achieve significantly lower read mapping errors and enhanced variant calling sensitivity, in addition to providing the improvements of joint variant calling without the need of computationally intensive post-processing steps.

DNA Methylation Analysis of Turner Syndrome BAV

Turner Syndrome (TS) is a rare cytogenetic disorder caused by the complete loss or structural variation of the second sex chromosome. The most common cause of early mortality in TS results from a high incidence of left-sided congenital heart defects, including bicuspid aortic valve (BAV), which occurs in about 30% of individuals with TS. BAV is also the most common congenital heart defect in the general population with a prevalence of 0.5–2%, with males being three-times more likely to have a BAV than females. TS is associated with genome-wide hypomethylation when compared to karyotypically normal males and females. Alterations in DNA methylation in primary aortic tissue are associated with BAV in euploid individuals. Here we show significant differences in DNA methylation patterns associated with BAV in TS found in peripheral blood by comparing TS BAV (n = 12), TS TAV (n = 13), and non-syndromic BAV (n = 6). When comparing TS with BAV to TS with no heart defects we identified a differentially methylated region encompassing the BAV-associated gene MYRF, and enrichment for binding sites of two known transcription factor contributors to BAV. When comparing TS with BAV to euploid women with BAV, we found significant overlapping enrichment for ChIP-seq transcription factor targets including genes in the NOTCH1 pathway, known for involvement in the etiology of non-syndromic BAV, and other genes that are essential regulators of heart valve development. Overall, these findings suggest that altered DNA methylation affecting key aortic valve development genes contributes to the greatly increased risk for BAV in TS.

SETD2: from chromatin modifier to multipronged regulator of the genome and beyond

Histone modifying enzymes play critical roles in many key cellular processes and are appealing proteins for targeting by small molecules in disease. However, while the functions of histone modifying enzymes are often linked to epigenetic regulation of the genome, an emerging theme is that these enzymes often also act by non-catalytic and/or non-epigenetic mechanisms. SETD2 (Set2 in yeast) is best known for associating with the transcription machinery and methylating histone H3 on lysine 36 (H3K36) during transcription. This well-characterized molecular function of SETD2 plays a role in fine-tuning transcription, maintaining chromatin integrity, and mRNA processing. Here we give an overview of the various molecular functions and mechanisms of regulation of H3K36 methylation by Set2/SETD2. These fundamental insights are important to understand SETD2’s role in disease, most notably in cancer in which SETD2 is frequently inactivated. SETD2 also methylates non-histone substrates such as α-tubulin which may promote genome stability and contribute to the tumor-suppressor function of SETD2. Thus, to understand its role in disease, it is important to understand and dissect the multiple roles of SETD2 within the cell. In this review we discuss how histone methylation by Set2/SETD2 has led the way in connecting histone modifications in active regions of the genome to chromatin functions and how SETD2 is leading the way to showing that we also have to look beyond histones to truly understand the physiological role of an ‘epigenetic’ writer enzyme in normal cells and in disease.

Challenges in Analyzing Functional Epigenetic Data in Perspective of Adolescent Psychiatric Health

The formative period of adolescence plays a crucial role in the development of skills and abilities for adulthood. Adolescents who are affected by mental health conditions are at risk of suicide and social and academic impairments. Gene–environment complementary contributions to the molecular mechanisms involved in psychiatric disorders have emphasized the need to analyze epigenetic marks such as DNA methylation (DNAm) and non-coding RNAs. However, the large and diverse bioinformatic and statistical methods, referring to the confounders of the statistical models, application of multiple-testing adjustment methods, questions regarding the correlation of DNAm across tissues, and sex-dependent differences in results, have raised challenges regarding the interpretation of the results. Based on the example of generalized anxiety disorder (GAD) and depressive disorder (MDD), we shed light on the current knowledge and usage of methodological tools in analyzing epigenetics. Statistical robustness is an essential prerequisite for a better understanding and interpretation of epigenetic modifications and helps to find novel targets for personalized therapeutics in psychiatric diseases.

An Amish founder population reveals rare-population genetic determinants of the human lipidome

Identifying the genetic determinants of inter-individual variation in lipid species (lipidome) may provide deeper understanding and additional insight into the mechanistic effect of complex lipidomic pathways in CVD risk and progression beyond simple traditional lipids. Previous studies have been largely population based and thus only powered to discover associations with common genetic variants. Founder populations represent a powerful resource to accelerate discovery of previously unknown biology associated with rare population alleles that have risen to higher frequency due to genetic drift. We performed a genome-wide association scan of 355 lipid species in 650 individuals from the Amish founder population including 127 lipid species not previously tested. To the best of our knowledge, we report for the first time the lipid species associated with two rare-population but Amish-enriched lipid variants: APOB_rs5742904 and APOC3_rs76353203. We also identified novel associations for 3 rare-population Amish-enriched loci with several sphingolipids and with proposed potential functional/causal variant in each locus including GLTPD2_rs536055318, CERS5_rs771033566, and AKNA_rs531892793. We replicated 7 previously known common loci including novel associations with two sterols: androstenediol with UGT locus and estriol with SLC22A8/A24 locus. Our results show the double power of founder populations and detailed lipidome to discover novel trait-associated variants.

A GWAS of 355 lipid species in the Old Order Amish founder population reveals associations between Amish-enriched loci and several sphingolipids.

The noncoding genome and hearing loss

The age of sequencing has provided unprecedented insights into the human genome. The coding region of the genome comprises nearly 20,000 genes, of which approximately 4000 are associated with human disease. Beyond the protein-coding genome, which accounts for only 3% of the genome, lies a vast pool of regulatory elements in the form of promoters, enhancers, RNA species, and other intricate elements. These features undoubtably influence human health and disease, and as a result, a great deal of effort is currently being invested in deciphering their identity and mechanism. While a paucity of material has caused a lag in identifying these elements in the inner ear, the emergence of technologies for dealing with a minimal number of cells now has the field working overtime to catch up. Studies on microRNAs (miRNAs), long non-coding RNAs (lncRNAs), methylation, histone modifications, and more are ongoing. A number of microRNAs and other noncoding elements are known to be associated with hearing impairment and there is promise that regulatory elements will serve as future tools and targets of therapeutics and diagnostics. This review covers the current state of the field and considers future directions for the noncoding genome and implications for hearing loss.

MicroRNAs and their delivery in diabetic fibrosis

The global prevalence of diabetes mellitus was estimated to be 463 million people in 2019 and is predicted to rise to 700 million by 2045. The associated financial and societal costs of this burgeoning epidemic demand an understanding of the pathology of this disease, and its complications, that will inform treatment to enable improved patient outcomes. Nearly two decades after the sequencing of the human genome, the significance of noncoding RNA expression is still being assessed. The family of functional noncoding RNAs known as microRNAs regulates the expression of most genes encoded by the human genome. Altered microRNA expression profiles have been observed both in diabetes and in diabetic complications. These transcripts therefore have significant potential and novelty as targets for therapy, therapeutic agents and biomarkers.

Why sequence all eukaryotes?

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.

OUP accepted manuscript

Mouse models have been engineered to reveal the biological mechanisms of human diseases based on an assumption. The assumption is that orthologous genes underlie conserved phenotypes across species. However, genetically modified mouse orthologs of human genes do not often recapitulate human disease phenotypes which might be due to the molecular evolution of phenotypic differences across species from the time of the last common ancestor. Here, we systematically investigated the evolutionary divergence of regulatory relationships between transcription factors (TFs) and target genes in functional modules, and found that the rewiring of gene regulatory networks (GRNs) contributes to the phenotypic discrepancies that occur between humans and mice. We confirmed that the rewired regulatory networks of orthologous genes contain a higher proportion of species-specific regulatory elements. Additionally, we verified that the divergence of target gene expression levels, which was triggered by network rewiring, could lead to phenotypic differences. Taken together, a careful consideration of evolutionary divergence in regulatory networks could be a novel strategy to understand the failure or success of mouse models to mimic human diseases. To help interpret mouse phenotypes in human disease studies, we provide quantitative comparisons of gene expression profiles on our website (http://sbi.postech.ac.kr/w/RN).

Transcriptional network orchestrating regional patterning of cortical progenitors

We uncovered a transcription factor (TF) network that regulates cortical regional patterning in radial glial stem cells. Screening the expression of hundreds of TFs in the developing mouse cortex identified 38 TFs that are expressed in gradients in the ventricular zone (VZ). We tested whether their cortical expression was altered in mutant mice with known patterning defects (Emx2, Nr2f1, and Pax6), which enabled us to define a cortical regionalization TF network (CRTFN). To identify genomic programming underlying this network, we performed TF ChIP-seq and chromatin-looping conformation to identify enhancer-gene interactions. To map enhancers involved in regional patterning of cortical progenitors, we performed assays for epigenomic marks and DNA accessibility in VZ cells purified from wild-type and patterning mutant mice. This integrated approach has identified a CRTFN and VZ enhancers involved in cortical regional patterning in the mouse.

Interrogating cell type-specific cooperation of transcriptional regulators in 3D chromatin

Context-specific activities of transcription regulators (TRs) in the nucleus modulate spatiotemporal gene expression precisely. Using the largest ChIP-seq data and chromatin loops in the human K562 cell line, we initially interrogated TR cooperation in 3D chromatin via a graphical model and revealed many known and novel TRs manipulating context-specific pathways. To explore TR cooperation across broad tissue/cell types, we systematically leveraged large-scale open chromatin profiles, computational footprinting, and high-resolution chromatin interactions to investigate tissue/cell type-specific TR cooperation. We first delineated a landscape of TR cooperation across 40 human tissue/cell types. Network modularity analyses uncovered the commonality and specificity of TR cooperation in different conditions. We also demonstrated that TR cooperation information can better interpret the disease-causal variants identified by genome-wide association studies and recapitulate cell states during neural development. Our study characterizes shared and unique patterns of TR cooperation associated with the cell type specificity of gene regulation in 3D chromatin.

•

Computational inference of transcriptional regulator (TR) cooperation in 3D chromatin

•

A landscape of 3D TR cooperation across 40 human tissue/cell types

•

TR cooperation can better interpret the disease-causal variants identified by GWAS

•

Cooperation of certain TRs shapes context-specific gene regulation in cell development