Deep sequencing of 10,000 human genomes

We report on the sequencing of 10,545 human genomes at 30×-40× coverage with an emphasis on quality metrics and novel variant and sequence discovery. We find that 84% of an individual human genome can be sequenced confidently. This high-confidence region includes 91.5% of exon sequence and 95.2% of known pathogenic variant positions. We present the distribution of over 150 million single-nucleotide variants in the coding and noncoding genome. Each newly sequenced genome contributes an average of 8,579 novel variants. In addition, each genome carries on average 0.7 Mb of sequence that is not found in the main build of the hg38 reference genome. The density of this catalog of variation allowed us to construct high-resolution profiles that define genomic sites that are highly intolerant of genetic variation. These results indicate that the data generated by deep genome sequencing is of the quality necessary for clinical use.

Role of a single noncoding nucleotide in the evolution of an epidemic African clade of Salmonella

Invasive nontyphoidal Salmonella disease is a major and previously neglected tropical disease responsible for an estimated ∼390,000 deaths per year in Africa, largely caused by a variant of Salmonella Typhimurium called ST313. Despite the availability of >100,000 Salmonella genomes, it has proven challenging to associate individual SNPs with pathogenic traits of this dangerous bacterium. Here, we used a transcriptomic strategy to identify a single-nucleotide change in a promoter region responsible for crucial phenotypic differences of African S. Typhimurium. Our findings show that a noncoding nucleotide of the bacterial genome can have a profound effect upon the pathogenesis of infectious disease.

Salmonella enterica serovar Typhimurium ST313 is a relatively newly emerged sequence type that is causing a devastating epidemic of bloodstream infections across sub-Saharan Africa. Analysis of hundreds of Salmonella genomes has revealed that ST313 is closely related to the ST19 group of S. Typhimurium that cause gastroenteritis across the world. The core genomes of ST313 and ST19 vary by only ∼1,000 SNPs. We hypothesized that the phenotypic differences that distinguish African Salmonella from ST19 are caused by certain SNPs that directly modulate the transcription of virulence genes. Here we identified 3,597 transcriptional start sites of the ST313 strain D23580, and searched for a gene-expression signature linked to pathogenesis of Salmonella. We identified a SNP in the promoter of the pgtE gene that caused high expression of the PgtE virulence factor in African S. Typhimurium, increased the degradation of the factor B component of human complement, contributed to serum resistance, and modulated virulence in the chicken infection model. We propose that high levels of PgtE expression by African S. Typhimurium ST313 promote bacterial survival and dissemination during human infection. Our finding of a functional role for an extragenic SNP shows that approaches used to deduce the evolution of virulence in bacterial pathogens should include a focus on noncoding regions of the genome.

High-Throughput Approaches to Pinpoint Function within the Noncoding Genome

The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas nuclease system is a powerful tool for genome editing, and its simple programmability has enabled high-throughput genetic and epigenetic studies. These high-throughput approaches offer investigators a toolkit for functional interrogation of not only protein-coding genes but also noncoding DNA. Historically, noncoding DNA has lacked the detailed characterization that has been applied to protein-coding genes in large part because there has not been a robust set of methodologies for perturbing these regions. Although the majority of high-throughput CRISPR screens have focused on the coding genome to date, an increasing number of CRISPR screens targeting noncoding genomic regions continue to emerge. Here, we review high-throughput CRISPR-based approaches to uncover and understand functional elements within the noncoding genome and discuss practical aspects of noncoding library design and screen analysis.

Coding of Glioblastoma Progression and Therapy Resistance through Long Noncoding RNAs

Glioblastoma is the most aggressive and lethal primary brain malignancy, with an average patient survival from diagnosis of 14 months. Glioblastoma also usually progresses as a more invasive phenotype after initial treatment. A major step forward in our understanding of the nature of glioblastoma was achieved with large-scale expression analysis. However, due to genomic complexity and heterogeneity, transcriptomics alone is not enough to define the glioblastoma "fingerprint", so epigenetic mechanisms are being examined, including the noncoding genome. On the basis of their tissue specificity, long noncoding RNAs (lncRNAs) are being explored as new diagnostic and therapeutic targets. In addition, growing evidence indicates that lncRNAs have various roles in resistance to glioblastoma therapies (e.g., MALAT1, H19) and in glioblastoma progression (e.g., CRNDE, HOTAIRM1, ASLNC22381, ASLNC20819). Investigations have also focused on the prognostic value of lncRNAs, as well as the definition of the molecular signatures of glioma, to provide more precise tumor classification. This review discusses the potential that lncRNAs hold for the development of novel diagnostic and, hopefully, therapeutic targets that can contribute to prolonged survival and improved quality of life for patients with glioblastoma.

Emerging roles of macrosatellite repeats in genome organization and disease development

Abundant repetitive DNA sequences are an enigmatic part of the human genome. Despite increasing evidence on the functionality of DNA repeats, their biologic role is still elusive and under frequent debate. Macrosatellites are the largest of the tandem DNA repeats, located on one or multiple chromosomes. The contribution of macrosatellites to genome regulation and human health was demonstrated for the D4Z4 macrosatellite repeat array on chromosome 4q35. Reduced copy number of D4Z4 repeats is associated with local euchromatinization and the onset of facioscapulohumeral muscular dystrophy. Although the role other macrosatellite families may play remains rather obscure, their diverse functionalities within the genome are being gradually revealed. In this review, we will outline structural and functional features of coding and noncoding macrosatellite repeats, and highlight recent findings that bring these sequences into the spotlight of genome organization and disease development.

Transcriptional enhancers at 40: evolution of a viral DNA element to nuclear architectural structures

Gene regulation by transcriptional enhancers is the dominant mechanism driving cell type- and signal-specific transcriptional diversity in metazoans. However, over four decades since the original discovery, how enhancers operate in the nuclear space remains largely enigmatic. Recent multidisciplinary efforts combining real-time imaging, genome sequencing, and biophysical strategies provide insightful but conflicting models of enhancer-mediated gene control. Here, we review the discovery and progress in enhancer biology, emphasizing the recent findings that acutely activated enhancers assemble regulatory machinery as mesoscale architectural structures with distinct physical properties. These findings help formulate novel models that explain several mysterious features of the assembly of transcriptional enhancers and the mechanisms of spatial control of gene expression.

Stabilizing selection on microsatellite allele length at arginine vasopressin 1a receptor and oxytocin receptor loci

The loci arginine vasopressin receptor 1a (avpr1a) and oxytocin receptor (oxtr) have evolutionarily conserved roles in vertebrate social and sexual behaviour. Allelic variation at a microsatellite locus in the 5′ regulatory region of these genes is associated with fitness in the bank vole Myodes glareolus. Given the low frequency of long and short alleles at these microsatellite loci in wild bank voles, we used breeding trials to determine whether selection acts against long and short alleles. Female bank voles with intermediate length avpr1a alleles had the highest probability of breeding, while male voles whose avpr1a alleles were very different in length had reduced probability of breeding. Moreover, there was a significant interaction between male and female oxtr genotypes, where potential breeding pairs with dissimilar length alleles had reduced probability of breeding. These data show how genetic variation at microsatellite loci associated with avpr1a and oxtr is associated with fitness, and highlight complex patterns of selection at these loci. More widely, these data show how stabilizing selection might act on allele length frequency distributions at gene-associated microsatellite loci.

Quantitative spatial and temporal assessment of regulatory element activity in zebrafish

Mutations or genetic variation in noncoding regions of the genome harbouring cis-regulatory elements (CREs), or enhancers, have been widely implicated in human disease and disease risk. However, our ability to assay the impact of these DNA sequence changes on enhancer activity is currently very limited because of the need to assay these elements in an appropriate biological context. Here, we describe a method for simultaneous quantitative assessment of the spatial and temporal activity of wild-type and disease-associated mutant human CRE alleles using live imaging in zebrafish embryonic development. We generated transgenic lines harbouring a dual-CRE dual-reporter cassette in a pre-defined neutral docking site in the zebrafish genome. The activity of each CRE allele is reported via expression of a specific fluorescent reporter, allowing simultaneous visualisation of where and when in development the wild-type allele is active and how this activity is altered by mutation.

tRNA-like Transcripts from the NEAT1-MALAT1 Genomic Region Critically Influence Human Innate Immunity and Macrophage Functions

The evolutionary conserved NEAT1-MALAT1 gene cluster generates large noncoding transcripts remaining nuclear, while tRNA-like transcripts (mascRNA, menRNA) enzymatically generated from these precursors translocate to the cytosol. Whereas functions have been assigned to the nuclear transcripts, data on biological functions of the small cytosolic transcripts are sparse. We previously found NEAT1^-/- and MALAT1^-/- mice to display massive atherosclerosis and vascular inflammation. Here, employing selective targeted disruption of menRNA or mascRNA, we investigate the tRNA-like molecules as critical components of innate immunity. CRISPR-generated human ΔmascRNA and ΔmenRNA monocytes/macrophages display defective innate immune sensing, loss of cytokine control, imbalance of growth/angiogenic factor expression impacting upon angiogenesis, and altered cell-cell interaction systems. Antiviral response, foam cell formation/oxLDL uptake, and M1/M2 polarization are defective in ΔmascRNA/ΔmenRNA macrophages, defining first biological functions of menRNA and describing new functions of mascRNA. menRNA and mascRNA represent novel components of innate immunity arising from the noncoding genome. They appear as prototypes of a new class of noncoding RNAs distinct from others (miRNAs, siRNAs) by biosynthetic pathway and intracellular kinetics. Their NEAT1-MALAT1 region of origin appears as archetype of a functionally highly integrated RNA processing system.

Fitness Landscape of the Fission Yeast Genome

The relationship between DNA sequence, biochemical function, and molecular evolution is relatively well-described for protein-coding regions of genomes, but far less clear in noncoding regions, particularly, in eukaryote genomes. In part, this is because we lack a complete description of the essential noncoding elements in a eukaryote genome. To contribute to this challenge, we used saturating transposon mutagenesis to interrogate the Schizosaccharomyces pombe genome. We generated 31 million transposon insertions, a theoretical coverage of 2.4 insertions per genomic site. We applied a five-state hidden Markov model (HMM) to distinguish insertion-depleted regions from insertion biases. Both raw insertion-density and HMM-defined fitness estimates showed significant quantitative relationships to gene knockout fitness, genetic diversity, divergence, and expected functional regions based on transcription and gene annotations. Through several analyses, we conclude that transposon insertions produced fitness effects in 66-90% of the genome, including substantial portions of the noncoding regions. Based on the HMM, we estimate that 10% of the insertion depleted sites in the genome showed no signal of conservation between species and were weakly transcribed, demonstrating limitations of comparative genomics and transcriptomics to detect functional units. In this species, 3'- and 5'-untranslated regions were the most prominent insertion-depleted regions that were not represented in measures of constraint from comparative genomics. We conclude that the combination of transposon mutagenesis, evolutionary, and biochemical data can provide new insights into the relationship between genome function and molecular evolution.

Ontogeny and heterogeneity of innate lymphoid cells and the noncoding genome

Since their discovery a decade ago, it has become evident that innate lymphoid cells (ILCs) play critical roles in protective immune responses against intracellular and extracellular pathogens but are also central regulators of epithelial barrier integrity and tissue homeostasis. ILCs populate almost every tissue in mammalian organisms; therefore, not surprisingly, dysregulation of their functions contributes to the development and progression of multiple inflammatory and metabolic diseases. Our knowledge of the transcriptional programs governing the development, differentiation, and functions of the different groups of ILCs has increased dramatically in the last ten years. However, with the advent of new technologies, an unprecedented level of heterogeneity, plasticity, and developmental complexity has started to be revealed. In this review, we highlight recent advances in our understanding of ILC development and their biological functions. In particular, we aim to emphasize how our increasing knowledge of the chromatin landscape and the noncoding genome of these innate lymphocytes is allowing us to better understand their development and functions in different contexts during homeostasis and inflammation. Moreover, we propose that the design of more refined genetic tools to study tissue-specific ILCs and their functions can be accomplished by leveraging our understanding of how specific noncoding elements of the genome regulate gene expression in ILCs.

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

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

Gauging the trends of pseudogenes in plants

Pseudogenes, the debilitated parts of ancient genes, were previously scrapped off as junk or discarded genes with no functional significance. Pseudogenes have come under scrutiny for their functionality, since recent studies have unveiled their importance in the regulation of their corresponding parent genes and various biological mechanisms. Despite the enormous occurrence of pseudogenes in plants, the lack of experimental validation has contributed toward their unresolved roles in gene regulation. Contrarily, most of the studies associated with gene regulation have been mainly reported for humans, mice, and other mammalian genomes. Consequently, in order to present a cumulative report on plant-based pseudogenes research, an attempt has been made to assemble multiple studies presenting the pseudogene classification, the prediction and the determination of comparative accuracies of various computational pipelines, and recent trends in analyzing their biological functions, and regulatory mechanisms. This review represents the classical, as well as the recent advances on pseudogene identification and their potential roles in transcriptional regulation, which could possibly invigorate the quality of genome annotation, evolutionary analysis, and complexity surrounding the regulatory pathways in plants. Thus, when the ambiguous boundary girdling the pseudogenes eventually recedes on account of their explicit orchestration role, research in flora would no longer saunter compared to that on fauna.

RetroFun-RVS: A Retrospective Family-Based Framework for Rare Variant Analysis Incorporating Functional Annotations

A large proportion of genetic variations involved in complex diseases are rare and located within noncoding regions, making the interpretation of underlying biological mechanisms a daunting task. Although technical and methodological progress has been made to annotate the genome, current disease-rare-variant association tests incorporating such annotations suffer from two major limitations. First, they are generally restricted to case-control designs of unrelated individuals, which often require tens or hundreds of thousands of individuals to achieve sufficient power. Second, they were not evaluated with region-based annotations needed to interpret the causal regulatory mechanisms. In this work, we propose RetroFun-RVS, a new retrospective family-based score test, incorporating functional annotations. A critical feature of the proposed method is to aggregate genotypes to compare against rare variant-sharing expectations among affected family members. Through extensive simulations, we have demonstrated that RetroFun-RVS integrating networks based on 3D genome contacts as functional annotations reach greater power over the region-wide test, other strategies to include subregions and competing methods. Also, the proposed framework shows robustness to non-informative annotations, maintaining its power when causal variants are spread across regions. Asymptotic p-values are susceptible to Type I error inflation when the number of families with rare variants is small, and a bootstrap procedure is recommended in these instances. Application of RetroFun-RVS is illustrated on whole genome sequence in the Eastern Quebec Schizophrenia and Bipolar Disorder Kindred Study with networks constructed from 3D contacts and epigenetic data on neurons. In summary, the integration of functional annotations corresponding to regions or networks with transcriptional impacts in rare variant tests appears promising to highlight regulatory mechanisms involved in complex diseases.

CWAS-Plus: estimating category-wide association of rare noncoding variation from whole-genome sequencing data with cell-type-specific functional data

Variants in cis-regulatory elements link the noncoding genome to human pathology; however, detailed analytic tools for understanding the association between cell-level brain pathology and noncoding variants are lacking. CWAS-Plus, adapted from a Python package for category-wide association testing (CWAS), enhances noncoding variant analysis by integrating both whole-genome sequencing (WGS) and user-provided functional data. With simplified parameter settings and an efficient multiple testing correction method, CWAS-Plus conducts the CWAS workflow 50 times faster than CWAS, making it more accessible and user-friendly for researchers. Here, we used a single-nuclei assay for transposase-accessible chromatin with sequencing to facilitate CWAS-guided noncoding variant analysis at cell-type-specific enhancers and promoters. Examining autism spectrum disorder WGS data (n = 7280), CWAS-Plus identified noncoding de novo variant associations in transcription factor binding sites within conserved loci. Independently, in Alzheimer's disease WGS data (n = 1087), CWAS-Plus detected rare noncoding variant associations in microglia-specific regulatory elements. These findings highlight CWAS-Plus's utility in genomic disorders and scalability for processing large-scale WGS data and in multiple-testing corrections. CWAS-Plus and its user manual are available at https://github.com/joonan-lab/cwas/ and https://cwas-plus.readthedocs.io/en/latest/, respectively.

CWAS-Plus: Estimating category-wide association of rare noncoding variation from whole-genome sequencing data with cell-type-specific functional data

Variants in cis-regulatory elements link the noncoding genome to human brain pathology; however, detailed analytic tools for understanding the association between cell-level brain pathology and noncoding variants are lacking. CWAS-Plus, adapted from a Python package for category-wide association testing (CWAS) employs both whole-genome sequencing and user-provided functional data to enhance noncoding variant analysis, with a faster and more efficient execution of the CWAS workflow. Here, we used single-nuclei assay for transposase-accessible chromatin with sequencing to facilitate CWAS-guided noncoding variant analysis at cell-type specific enhancers and promoters. Examining autism spectrum disorder whole-genome sequencing data (n = 7,280), CWAS-Plus identified noncoding de novo variant associations in transcription factor binding sites within conserved loci. Independently, in Alzheimer's disease whole-genome sequencing data (n = 1,087), CWAS-Plus detected rare noncoding variant associations in microglia-specific regulatory elements. These findings highlight CWAS-Plus's utility in genomic disorders and scalability for processing large-scale whole-genome sequencing data and in multiple-testing corrections. CWAS-Plus and its user manual are available at https://github.com/joonan-lab/cwas/ and https://cwas-plus.readthedocs.io/en/latest/, respectively.

Innate Immunity in Cardiovascular Diseases-Identification of Novel Molecular Players and Targets

During the past few years, unexpected developments have driven studies in the field of clinical immunology. One driver of immense impact was the outbreak of a pandemic caused by the novel virus SARS-CoV-2. Excellent recent reviews address diverse aspects of immunological re-search into cardiovascular diseases. Here, we specifically focus on selected studies taking advantage of advanced state-of-the-art molecular genetic methods ranging from genome-wide epi/transcriptome mapping and variant scanning to optogenetics and chemogenetics. First, we discuss the emerging clinical relevance of advanced diagnostics for cardiovascular diseases, including those associated with COVID-19-with a focus on the role of inflammation in cardiomyopathies and arrhythmias. Second, we consider newly identified immunological interactions at organ and system levels which affect cardiovascular pathogenesis. Thus, studies into immune influences arising from the intestinal system are moving towards therapeutic exploitation. Further, powerful new research tools have enabled novel insight into brain-immune system interactions at unprecedented resolution. This latter line of investigation emphasizes the strength of influence of emotional stress-acting through defined brain regions-upon viral and cardiovascular disorders. Several challenges need to be overcome before the full impact of these far-reaching new findings will hit the clinical arena.

Peeling the onion: how complex is the artemisinin resistance genetic trait of malaria parasites?

The genetics of Plasmodium as an intracellular, mostly haploid, sexually reproducing, eukaryotic organism with a complex life cycle, presents unprecedented challenges in studying drug resistance. This article summarizes current knowledge on the genetic basis of artemisinin resistance (AR) - a main component of current drug therapies for falciparum malaria. Although centered on nonsynonymous single-nucleotide polymorphisms (nsSNPs), we describe multifaceted resistance mechanisms as part of a complex, cumulative genetic trait that involves regulation of expression by a wide array of polymorphisms in noncoding regions. These genetic variations alter transcriptome profiles linked to Plasmodium's development and population dynamics, ultimately influencing the emergence and spread of the resistance.