Common and rare variant associations with clonal haematopoiesis phenotypes.

Clonal haematopoiesis involves the expansion of certain blood cell lineages and has been associated with ageing and adverse health outcomes1-5. Here we use exome sequence data on 628,388 individuals to identify 40,208 carriers of clonal haematopoiesis of indeterminate potential (CHIP). Using genome-wide and exome-wide association analyses, we identify 24 loci (21 of which are novel) where germline genetic variation influences predisposition to CHIP, including missense variants in the lymphocytic antigen coding gene LY75, which are associated with reduced incidence of CHIP. We also identify novel rare variant associations with clonal haematopoiesis and telomere length. Analysis of 5,041 health traits from the UK Biobank (UKB) found relationships between CHIP and severe COVID-19 outcomes, cardiovascular disease, haematologic traits, malignancy, smoking, obesity, infection and all-cause mortality. Longitudinal and Mendelian randomization analyses revealed that CHIP is associated with solid cancers, including non-melanoma skin cancer and lung cancer, and that CHIP linked to DNMT3A is associated with the subsequent development of myeloid but not lymphoid leukaemias. Additionally, contrary to previous findings from the initial 50,000 UKB exomes6, our results in the full sample do not support a role for IL-6 inhibition in reducing the risk of cardiovascular disease among CHIP carriers. Our findings demonstrate that CHIP represents a complex set of heterogeneous phenotypes with shared and unique germline genetic causes and varied clinical implications.

Pan-ancestry exome-wide association analyses of COVID-19 outcomes in 586,157 individuals.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19), a respiratory illness that can result in hospitalization or death. We used exome sequence data to investigate associations between rare genetic variants and seven COVID-19 outcomes in 586,157 individuals, including 20,952 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome wide or when specifically focusing on (1) 13 interferon pathway genes in which rare deleterious variants have been reported in individuals with severe COVID-19, (2) 281 genes located in susceptibility loci identified by the COVID-19 Host Genetics Initiative, or (3) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, and results are publicly available through the Regeneron Genetics Center COVID-19 Results Browser.

Exome sequencing and characterization of 49,960 individuals in the UK Biobank.

The UK Biobank is a prospective study of 502,543 individuals, combining extensive phenotypic and genotypic data with streamlined access for researchers around the world1. Here we describe the release of exome-sequence data for the first 49,960 study participants, revealing approximately 4 million coding variants (of which around 98.6% have a frequency of less than 1%). The data include 198,269 autosomal predicted loss-of-function (LOF) variants, a more than 14-fold increase compared to the imputed sequence. Nearly all genes (more than 97%) had at least one carrier with a LOF variant, and most genes (more than 69%) had at least ten carriers with a LOF variant. We illustrate the power of characterizing LOF variants in this population through association analyses across 1,730 phenotypes. In addition to replicating established associations, we found novel LOF variants with large effects on disease traits, including PIEZO1 on varicose veins, COL6A1 on corneal resistance, MEPE on bone density, and IQGAP2 and GMPR on blood cell traits. We further demonstrate the value of exome sequencing by surveying the prevalence of pathogenic variants of clinical importance, and show that 2% of this population has a medically actionable variant. Furthermore, we characterize the penetrance of cancer in carriers of pathogenic BRCA1 and BRCA2 variants. Exome sequences from the first 49,960 participants highlight the promise of genome sequencing in large population-based studies and are now accessible to the scientific community.

The Loss of TBK1 Kinase Activity in Motor Neurons or in All Cell Types Differentially Impacts ALS Disease Progression in SOD1 Mice.

DNA sequence variants in the TBK1 gene associate with or cause sporadic or familial amyotrophic lateral sclerosis (ALS). Here we show that mice bearing human ALS-associated TBK1 missense loss-of-function mutations, or mice in which the Tbk1 gene is selectively deleted in motor neurons, do not display a neurodegenerative disease phenotype. However, loss of TBK1 function in motor neurons of the SOD1G93A mouse model of ALS impairs autophagy, increases SOD1 aggregation, and accelerates early disease onset without affecting lifespan. By contrast, point mutations that decrease TBK1 kinase activity in all cells also accelerate disease onset but extend the lifespan of SOD1 mice. This difference correlates with the failure to activate high levels of expression of interferon-inducible genes in glia. We conclude that loss of TBK1 kinase activity impacts ALS disease progression through distinct pathways in different spinal cord cell types and further implicate the importance of glia in neurodegeneration.

Sour Sensing from the Tongue to the Brain.

The ability to sense sour provides an important sensory signal to prevent the ingestion of unripe, spoiled, or fermented foods. Taste and somatosensory receptors in the oral cavity trigger aversive behaviors in response to acid stimuli. Here, we show that the ion channel Otopetrin-1, a proton-selective channel normally involved in the sensation of gravity in the vestibular system, is essential for sour sensing in the taste system. We demonstrate that knockout of Otop1 eliminates acid responses from sour-sensing taste receptor cells (TRCs). In addition, we show that mice engineered to express otopetrin-1 in sweet TRCs have sweet cells that also respond to sour stimuli. Next, we genetically identified the taste ganglion neurons mediating each of the five basic taste qualities and demonstrate that sour taste uses its own dedicated labeled line from TRCs in the tongue to finely tuned taste neurons in the brain to trigger aversive behaviors.

Antisense lncRNA Transcription Mediates DNA Demethylation to Drive Stochastic Protocadherin α Promoter Choice.

Stochastic activation of clustered Protocadherin (Pcdh) α, ß, and γ genes generates a cell-surface identity code in individual neurons that functions in neural circuit assembly. Here, we show that Pcdhα gene choice involves the activation of an antisense promoter located in the first exon of each Pcdhα alternate gene. Transcription of an antisense long noncoding RNA (lncRNA) from this antisense promoter extends through the sense promoter, leading to DNA demethylation of the CTCF binding sites proximal to each promoter. Demethylation-dependent CTCF binding to both promoters facilitates cohesin-mediated DNA looping with a distal enhancer (HS5-1), locking in the transcriptional state of the chosen Pcdhα gene. Uncoupling DNA demethylation from antisense transcription by Tet3 overexpression in mouse olfactory neurons promotes CTCF binding to all Pcdhα promoters, resulting in proximity-biased DNA looping of the HS5-1 enhancer. Thus, antisense transcription-mediated promoter demethylation functions as a mechanism for distance-independent enhancer/promoter DNA looping to ensure stochastic Pcdhα promoter choice.

Pcdhαc2 is required for axonal tiling and assembly of serotonergic circuitries in mice.

Serotonergic neurons project their axons pervasively throughout the brain and innervate various target fields in a space-filling manner, leading to tiled arrangements of their axon terminals to allow optimal allocation of serotonin among target neurons. Here we show that conditional deletion of the mouse protocadherin α (Pcdhα) gene cluster in serotonergic neurons disrupts local axonal tiling and global assembly of serotonergic circuitries and results in depression-like behaviors. Genetic dissection and expression profiling revealed that this role is specifically mediated by Pcdhαc2, which is the only Pcdhα isoform expressed in serotonergic neurons. We conclude that, in contrast to neurite self-avoidance, which requires single-cell identity mediated by Pcdh diversity, a single cell-type identity mediated by the common C-type Pcdh isoform is required for axonal tiling and assembly of serotonergic circuitries.

Multicluster Pcdh diversity is required for mouse olfactory neural circuit assembly.

The vertebrate clustered protocadherin (Pcdh) cell surface proteins are encoded by three closely linked gene clusters (Pcdhα, Pcdhß, and Pcdhγ). Here, we show that all three gene clusters functionally cooperate to provide individual mouse olfactory sensory neurons (OSNs) with the cell surface diversity required for their assembly into distinct glomeruli in the olfactory bulb. Although deletion of individual Pcdh clusters had subtle phenotypic consequences, the loss of all three clusters (tricluster deletion) led to a severe axonal arborization defect and loss of self-avoidance. By contrast, when endogenous Pcdh diversity is overridden by the expression of a single-tricluster gene repertoire (α and ß and γ), OSN axons fail to converge to form glomeruli, likely owing to contact-mediated repulsion between axons expressing identical combinations of Pcdh isoforms.

hnRNP U protein is required for normal pre-mRNA splicing and postnatal heart development and function.

We report that mice lacking the heterogeneous nuclear ribonucleoprotein U (hnRNP U) in the heart develop lethal dilated cardiomyopathy and display numerous defects in cardiac pre-mRNA splicing. Mutant hearts have disorganized cardiomyocytes, impaired contractility, and abnormal excitation-contraction coupling activities. RNA-seq analyses of Hnrnpu mutant hearts revealed extensive defects in alternative splicing of pre-mRNAs encoding proteins known to be critical for normal heart development and function, including Titin and calcium/calmodulin-dependent protein kinase II delta (Camk2d). Loss of hnRNP U expression in cardiomyocytes also leads to aberrant splicing of the pre-mRNA encoding the excitation-contraction coupling component Junctin. We found that the protein product of an alternatively spliced Junctin isoform is N-glycosylated at a specific asparagine site that is required for interactions with specific protein partners. Our findings provide conclusive evidence for the essential role of hnRNP U in heart development and function and in the regulation of alternative splicing.

An RNA-sequencing transcriptome and splicing database of glia, neurons, and vascular cells of the cerebral cortex.

The major cell classes of the brain differ in their developmental processes, metabolism, signaling, and function. To better understand the functions and interactions of the cell types that comprise these classes, we acutely purified representative populations of neurons, astrocytes, oligodendrocyte precursor cells, newly formed oligodendrocytes, myelinating oligodendrocytes, microglia, endothelial cells, and pericytes from mouse cerebral cortex. We generated a transcriptome database for these eight cell types by RNA sequencing and used a sensitive algorithm to detect alternative splicing events in each cell type. Bioinformatic analyses identified thousands of new cell type-enriched genes and splicing isoforms that will provide novel markers for cell identification, tools for genetic manipulation, and insights into the biology of the brain. For example, our data provide clues as to how neurons and astrocytes differ in their ability to dynamically regulate glycolytic flux and lactate generation attributable to unique splicing of PKM2, the gene encoding the glycolytic enzyme pyruvate kinase. This dataset will provide a powerful new resource for understanding the development and function of the brain. To ensure the widespread distribution of these datasets, we have created a user-friendly website (http://web.stanford.edu/group/barres_lab/brain_rnaseq.html) that provides a platform for analyzing and comparing transciption and alternative splicing profiles for various cell classes in the brain.

Global effects of the CSR-1 RNA interference pathway on the transcriptional landscape.

Argonaute proteins and their small RNA cofactors short interfering RNAs are known to inhibit gene expression at the transcriptional and post-transcriptional levels. In Caenorhabditis elegans, the Argonaute CSR-1 binds thousands of endogenous siRNAs (endo-siRNAs) that are antisense to germline transcripts. However, its role in gene expression regulation remains controversial. Here we used genome-wide profiling of nascent RNA transcripts and found that the CSR-1 RNA interference pathway promoted sense-oriented RNA polymerase II transcription. Moreover, a loss of CSR-1 function resulted in global increase in antisense transcription and ectopic transcription of silent chromatin domains, which led to reduced chromatin incorporation of centromere-specific histone H3. On the basis of these findings, we propose that the CSR-1 pathway helps maintain the directionality of active transcription, thereby propagating the distinction between transcriptionally active and silent genomic regions.

A neurodegeneration-specific gene-expression signature of acutely isolated microglia from an amyotrophic lateral sclerosis mouse model.

Microglia are resident immune cells of the CNS that are activated by infection, neuronal injury, and inflammation. Here, we utilize flow cytometry and deep RNA sequencing of acutely isolated spinal cord microglia to define their activation in vivo. Analysis of resting microglia identified 29 genes that distinguish microglia from other CNS cells and peripheral macrophages/monocytes. We then analyzed molecular changes in microglia during neurodegenerative disease activation using the SOD1(G93A) mouse model of amyotrophic lateral sclerosis (ALS). We found that SOD1(G93A) microglia are not derived from infiltrating monocytes, and that both potentially neuroprotective and toxic factors, including Alzheimer's disease genes, are concurrently upregulated. Mutant microglia differed from SOD1(WT), lipopolysaccharide-activated microglia, and M1/M2 macrophages, defining an ALS-specific phenotype. Concurrent messenger RNA/fluorescence-activated cell sorting analysis revealed posttranscriptional regulation of microglia surface receptors and T cell-associated changes in the transcriptome. These results provide insights into microglia biology and establish a resource for future studies of neuroinflammation.

Intricate interplay between astrocytes and motor neurons in ALS.

ALS results from the selective and progressive degeneration of motor neurons. Although the underlying disease mechanisms remain unknown, glial cells have been implicated in ALS disease progression. Here, we examine the effects of glial cell/motor neuron interactions on gene expression using the hSOD1(G93A) (the G93A allele of the human superoxide dismutase gene) mouse model of ALS. We detect striking cell autonomous and nonautonomous changes in gene expression in cocultured motor neurons and glia, revealing that the two cell types profoundly affect each other. In addition, we found a remarkable concordance between the cell culture data and expression profiles of whole spinal cords and acutely isolated spinal cord cells during disease progression in the G93A mouse model, providing validation of the cell culture approach. Bioinformatics analyses identified changes in the expression of specific genes and signaling pathways that may contribute to motor neuron degeneration in ALS, among which are TGF-ß signaling pathways.

Accelerated high-yield generation of limb-innervating motor neurons from human stem cells.

Human pluripotent stem cells are a promising source of differentiated cells for developmental studies, cell transplantation, disease modeling, and drug testing. However, their widespread use even for intensely studied cell types like spinal motor neurons is hindered by the long duration and low yields of existing protocols for in vitro differentiation and by the molecular heterogeneity of the populations generated. We report a combination of small molecules that within 3 weeks induce motor neurons at up to 50% abundance and with defined subtype identities of relevance to neurodegenerative disease. Despite their accelerated differentiation, motor neurons expressed combinations of HB9, ISL1, and column-specific markers that mirror those observed in vivo in human embryonic spinal cord. They also exhibited spontaneous and induced activity, and projected axons toward muscles when grafted into developing chick spinal cord. Strikingly, this novel protocol preferentially generates motor neurons expressing markers of limb-innervating lateral motor column motor neurons (FOXP1(+)/LHX3(-)). Access to high-yield cultures of human limb-innervating motor neuron subtypes will facilitate in-depth study of motor neuron subtype-specific properties, disease modeling, and development of large-scale cell-based screening assays.

Transcription factor binding predictions using TRAP for the analysis of ChIP-seq data and regulatory SNPs.

The transcription factor affinity prediction (TRAP) method calculates the affinity of transcription factors for DNA sequences on the basis of a biophysical model. This method has proven to be useful for several applications, including for determining the putative target genes of a given factor. This protocol covers two other applications: (i) determining which transcription factors have the highest affinity in a set of sequences (illustrated with chromatin immunoprecipitation-sequencing (ChIP-seq) peaks), and (ii) finding which factor is the most affected by a regulatory single-nucleotide polymorphism. The protocol describes how to use the TRAP web tools to address these questions, and it also presents a way to run TRAP on random control sequences to better estimate the significance of the results. All of the tools are fully available online and do not need any additional installation. The complete protocol takes about 45 min, but each individual tool runs in a few minutes.

Combinatorial binding in human and mouse embryonic stem cells identifies conserved enhancers active in early embryonic development.

Transcription factors are proteins that regulate gene expression by binding to cis-regulatory sequences such as promoters and enhancers. In embryonic stem (ES) cells, binding of the transcription factors OCT4, SOX2 and NANOG is essential to maintain the capacity of the cells to differentiate into any cell type of the developing embryo. It is known that transcription factors interact to regulate gene expression. In this study we show that combinatorial binding is strongly associated with co-localization of the transcriptional co-activator Mediator, H3K27ac and increased expression of nearby genes in embryonic stem cells. We observe that the same loci bound by Oct4, Nanog and Sox2 in ES cells frequently drive expression in early embryonic development. Comparison of mouse and human ES cells shows that less than 5% of individual binding events for OCT4, SOX2 and NANOG are shared between species. In contrast, about 15% of combinatorial binding events and even between 53% and 63% of combinatorial binding events at enhancers active in early development are conserved. Our analysis suggests that the combination of OCT4, SOX2 and NANOG binding is critical for transcription in ES cells and likely plays an important role for embryogenesis by binding at conserved early developmental enhancers. Our data suggests that the fast evolutionary rewiring of regulatory networks mainly affects individual binding events, whereas "gene regulatory hotspots" which are bound by multiple factors and active in multiple tissues throughout early development are under stronger evolutionary constraints.

PASTAA: identifying transcription factors associated with sets of co-regulated genes.

MOTIVATION: A major challenge in regulatory genomics is the identification of associations between functional categories of genes (e.g. tissues, metabolic pathways) and their regulating transcription factors (TFs). While, for a limited number of categories, the regulating TFs are already known, still for many functional categories the responsible factors remain to be elucidated. RESULTS: We put forward a novel method (PASTAA) for detecting transcriptions factors associated with functional categories, which utilizes the prediction of binding affinities of a TF to promoters. This binding strength information is compared to the likelihood of membership of the corresponding genes in the functional category under study. Coherence between the two ranked datasets is seen as an indicator of association between a TF and the category. PASTAA is applied primarily to the determination of TFs driving tissue-specific expression. We show that PASTAA is capable of recovering many TFs acting tissue specifically and, in addition, provides novel associations so far not detected by alternative methods. The application of PASTAA to detect TFs involved in the regulation of tissue-specific gene expression revealed a remarkable number of experimentally supported associations. The validated success for various datasets implies that PASTAA can directly be applied for the detection of TFs associated with newly derived gene sets. AVAILABILITY: The PASTAA source code as well as a corresponding web interface is freely available at http://trap.molgen.mpg.de.

A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome.

The functional complexity of the human transcriptome is not yet fully elucidated. We report a high-throughput sequence of the human transcriptome from a human embryonic kidney and a B cell line. We used shotgun sequencing of transcripts to generate randomly distributed reads. Of these, 50% mapped to unique genomic locations, of which 80% corresponded to known exons. We found that 66% of the polyadenylated transcriptome mapped to known genes and 34% to nonannotated genomic regions. On the basis of known transcripts, RNA-Seq can detect 25% more genes than can microarrays. A global survey of messenger RNA splicing events identified 94,241 splice junctions (4096 of which were previously unidentified) and showed that exon skipping is the most prevalent form of alternative splicing.

The RNA transport element RTE is essential for IAP LTR-retrotransposon mobility.

We previously identified an RNA transport element (RTE) present at a high copy number in the mouse genome. Here, we show that a related element, RTE-D, is part of a mobile LTR-retrotransposon, which belongs to a family of intracisternal A-particle related elements (IAP). We demonstrate that RTE-D is essential for the mobility of the retrotransposon and it can be substituted by other known RNA export signals. RTE-deficient IAP transcripts are retained in the nucleus, while the RTE-containing transcripts accumulate in the cytoplasm allowing Gag protein expression. RTE-D acts as a posttranscriptional control element in a heterologous reporter mRNA and is activated by the cellular RNA binding protein 15 (RBM15), as reported for the previously described RTE. We identified a complex family of RTE-containing IAPs in mouse and mapped the active RTE-D-containing IAPs to the Mmr10 group of LTR-retrotransposons. These data reveal that, despite a complex evolutionary history, retroelements and retroviruses share the dependency on posttranscriptional regulation.