Genetic and molecular architecture of complex traits.

Human genetics has emerged as one of the most dynamic areas of biology, with a broadening societal impact. In this review, we discuss recent achievements, ongoing efforts, and future challenges in the field. Advances in technology, statistical methods, and the growing scale of research efforts have all provided many insights into the processes that have given rise to the current patterns of genetic variation. Vast maps of genetic associations with human traits and diseases have allowed characterization of their genetic architecture. Finally, studies of molecular and cellular effects of genetic variants have provided insights into biological processes underlying disease. Many outstanding questions remain, but the field is well poised for groundbreaking discoveries as it increases the use of genetic data to understand both the history of our species and its applications to improve human health.

3D reconstruction of a gastrulating human embryo.

Progress in understanding early human development has been impeded by the scarcity of reference datasets from natural embryos, particularly those with spatial information during crucial stages like gastrulation. We conducted high-resolution spatial transcriptomics profiling on 38,562 spots from 62 transverse sections of an intact Carnegie stage (CS) 8 human embryo. From this spatial transcriptomic dataset, we constructed a 3D model of the CS8 embryo, in which a range of cell subtypes are identified, based on gene expression patterns and positional register, along the anterior-posterior, medial-lateral, and dorsal-ventral axis in the embryo. We further characterized the lineage trajectories of embryonic and extra-embryonic tissues and associated regulons and the regionalization of signaling centers and signaling activities that underpin lineage progression and tissue patterning during gastrulation. Collectively, the findings of this study provide insights into gastrulation and post-gastrulation development of the human embryo.

A bat MERS-like coronavirus circulates in pangolins and utilizes human DPP4 and host proteases for cell entry.

It is unknown whether pangolins, the most trafficked mammals, play a role in the zoonotic transmission of bat coronaviruses. We report the circulation of a novel MERS-like coronavirus in Malayan pangolins, named Manis javanica HKU4-related coronavirus (MjHKU4r-CoV). Among 86 animals, four tested positive by pan-CoV PCR, and seven tested seropositive (11 and 12.8%). Four nearly identical (99.9%) genome sequences were obtained, and one virus was isolated (MjHKU4r-CoV-1). This virus utilizes human dipeptidyl peptidase-4 (hDPP4) as a receptor and host proteases for cell infection, which is enhanced by a furin cleavage site that is absent in all known bat HKU4r-CoVs. The MjHKU4r-CoV-1 spike shows higher binding affinity for hDPP4, and MjHKU4r-CoV-1 has a wider host range than bat HKU4-CoV. MjHKU4r-CoV-1 is infectious and pathogenic in human airways and intestinal organs and in hDPP4-transgenic mice. Our study highlights the importance of pangolins as reservoir hosts of coronaviruses poised for human disease emergence.

HGT is widespread in insects and contributes to male courtship in lepidopterans.

Horizontal gene transfer (HGT) is an important evolutionary force shaping prokaryotic and eukaryotic genomes. HGT-acquired genes have been sporadically reported in insects, a lineage containing >50% of animals. We systematically examined HGT in 218 high-quality genomes of diverse insects and found that they acquired 1,410 genes exhibiting diverse functions, including many not previously reported, via 741 distinct transfers from non-metazoan donors. Lepidopterans had the highest average number of HGT-acquired genes. HGT-acquired genes containing introns exhibited substantially higher expression levels than genes lacking introns, suggesting that intron gains were likely involved in HGT adaptation. Lastly, we used the CRISPR-Cas9 system to edit the prevalent unreported gene LOC105383139, which was transferred into the last common ancestor of moths and butterflies. In diamondback moths, males lacking LOC105383139 courted females significantly less. We conclude that HGT has been a major contributor to insect adaptation.

Reconstruction of neocortex: Organelles, compartments, cells, circuits, and activity.

We assembled a semi-automated reconstruction of L2/3 mouse primary visual cortex from ∼250 × 140 × 90 µm3 of electron microscopic images, including pyramidal and non-pyramidal neurons, astrocytes, microglia, oligodendrocytes and precursors, pericytes, vasculature, nuclei, mitochondria, and synapses. Visual responses of a subset of pyramidal cells are included. The data are publicly available, along with tools for programmatic and three-dimensional interactive access. Brief vignettes illustrate the breadth of potential applications relating structure to function in cortical circuits and neuronal cell biology. Mitochondria and synapse organization are characterized as a function of path length from the soma. Pyramidal connectivity motif frequencies are predicted accurately using a configuration model of random graphs. Pyramidal cells receiving more connections from nearby cells exhibit stronger and more reliable visual responses. Sample code shows data access and analysis.

A single-cell atlas of chromatin accessibility in the human genome.

Current catalogs of regulatory sequences in the human genome are still incomplete and lack cell type resolution. To profile the activity of gene regulatory elements in diverse cell types and tissues in the human body, we applied single-cell chromatin accessibility assays to 30 adult human tissue types from multiple donors. We integrated these datasets with previous single-cell chromatin accessibility data from 15 fetal tissue types to reveal the status of open chromatin for ∼1.2 million candidate cis-regulatory elements (cCREs) in 222 distinct cell types comprised of >1.3 million nuclei. We used these chromatin accessibility maps to delineate cell-type-specificity of fetal and adult human cCREs and to systematically interpret the noncoding variants associated with complex human traits and diseases. This rich resource provides a foundation for the analysis of gene regulatory programs in human cell types across tissues, life stages, and organ systems.

The Allen Mouse Brain Common Coordinate Framework: A 3D Reference Atlas.

Recent large-scale collaborations are generating major surveys of cell types and connections in the mouse brain, collecting large amounts of data across modalities, spatial scales, and brain areas. Successful integration of these data requires a standard 3D reference atlas. Here, we present the Allen Mouse Brain Common Coordinate Framework (CCFv3) as such a resource. We constructed an average template brain at 10 µm voxel resolution by interpolating high resolution in-plane serial two-photon tomography images with 100 µm z-sampling from 1,675 young adult C57BL/6J mice. Then, using multimodal reference data, we parcellated the entire brain directly in 3D, labeling every voxel with a brain structure spanning 43 isocortical areas and their layers, 329 subcortical gray matter structures, 81 fiber tracts, and 8 ventricular structures. CCFv3 can be used to analyze, visualize, and integrate multimodal and multiscale datasets in 3D and is openly accessible (https://atlas.brain-map.org/).

Trans Effects on Gene Expression Can Drive Omnigenic Inheritance.

Early genome-wide association studies (GWASs) led to the surprising discovery that, for typical complex traits, most of the heritability is due to huge numbers of common variants with tiny effect sizes. Previously, we argued that new models are needed to understand these patterns. Here, we provide a formal model in which genetic contributions to complex traits are partitioned into direct effects from core genes and indirect effects from peripheral genes acting in trans. We propose that most heritability is driven by weak trans-eQTL SNPs, whose effects are mediated through peripheral genes to impact the expression of core genes. In particular, if the core genes for a trait tend to be co-regulated, then the effects of peripheral variation can be amplified such that nearly all of the genetic variance is driven by weak trans effects. Thus, our model proposes a framework for understanding key features of the architecture of complex traits.

Predicting Splicing from Primary Sequence with Deep Learning.

The splicing of pre-mRNAs into mature transcripts is remarkable for its precision, but the mechanisms by which the cellular machinery achieves such specificity are incompletely understood. Here, we describe a deep neural network that accurately predicts splice junctions from an arbitrary pre-mRNA transcript sequence, enabling precise prediction of noncoding genetic variants that cause cryptic splicing. Synonymous and intronic mutations with predicted splice-altering consequence validate at a high rate on RNA-seq and are strongly deleterious in the human population. De novo mutations with predicted splice-altering consequence are significantly enriched in patients with autism and intellectual disability compared to healthy controls and validate against RNA-seq in 21 out of 28 of these patients. We estimate that 9%-11% of pathogenic mutations in patients with rare genetic disorders are caused by this previously underappreciated class of disease variation.

Characterizing the Major Structural Variant Alleles of the Human Genome.

In order to provide a comprehensive resource for human structural variants (SVs), we generated long-read sequence data and analyzed SVs for fifteen human genomes. We sequence resolved 99,604 insertions, deletions, and inversions including 2,238 (1.6 Mbp) that are shared among all discovery genomes with an additional 13,053 (6.9 Mbp) present in the majority, indicating minor alleles or errors in the reference. Genotyping in 440 additional genomes confirms the most common SVs in unique euchromatin are now sequence resolved. We report a ninefold SV bias toward the last 5 Mbp of human chromosomes with nearly 55% of all VNTRs (variable number of tandem repeats) mapping to this portion of the genome. We identify SVs affecting coding and noncoding regulatory loci improving annotation and interpretation of functional variation. These data provide the framework to construct a canonical human reference and a resource for developing advanced representations capable of capturing allelic diversity.

A Compendium of Genetic Modifiers of Mitochondrial Dysfunction Reveals Intra-organelle Buffering.

Mitochondrial dysfunction is associated with a spectrum of human conditions, ranging from rare, inborn errors of metabolism to the aging process. To identify pathways that modify mitochondrial dysfunction, we performed genome-wide CRISPR screens in the presence of small-molecule mitochondrial inhibitors. We report a compendium of chemical-genetic interactions involving 191 distinct genetic modifiers, including 38 that are synthetic sick/lethal and 63 that are suppressors. Genes involved in glycolysis (PFKP), pentose phosphate pathway (G6PD), and defense against lipid peroxidation (GPX4) scored high as synthetic sick/lethal. A surprisingly large fraction of suppressors are pathway intrinsic and encode mitochondrial proteins. A striking example of such "intra-organelle" buffering is the alleviation of a chemical defect in complex V by simultaneous inhibition of complex I, which benefits cells by rebalancing redox cofactors, increasing reductive carboxylation, and promoting glycolysis. Perhaps paradoxically, certain forms of mitochondrial dysfunction may best be buffered with "second site" inhibitors to the organelle.

Deficiency of CBL and CBLB ubiquitin ligases leads to hyper T follicular helper cell responses and lupus by reducing BCL6 degradation.

Recent evidence reveals hyper T follicular helper (Tfh) cell responses in systemic lupus erythematosus (SLE); however, molecular mechanisms responsible for hyper Tfh cell responses and whether they cause SLE are unclear. We found that SLE patients downregulated both ubiquitin ligases, casitas B-lineage lymphoma (CBL) and CBLB (CBLs), in CD4+ T cells. T cell-specific CBLs-deficient mice developed hyper Tfh cell responses and SLE, whereas blockade of Tfh cell development in the mutant mice was sufficient to prevent SLE. ICOS was upregulated in SLE Tfh cells, whose signaling increased BCL6 by attenuating BCL6 degradation via chaperone-mediated autophagy (CMA). Conversely, CBLs restrained BCL6 expression by ubiquitinating ICOS. Blockade of BCL6 degradation was sufficient to enhance Tfh cell responses. Thus, the compromised expression of CBLs is a prevalent risk trait shared by SLE patients and causative to hyper Tfh cell responses and SLE. The ICOS-CBLs axis may be a target to treat SLE.

Multi-omics analysis of innate and adaptive responses to BCG vaccination reveals epigenetic cell states that predict trained immunity.

Immune responses are tightly regulated yet highly variable between individuals. To investigate human population variation of trained immunity, we immunized healthy individuals with Bacillus Calmette-Guérin (BCG). This live-attenuated vaccine induces not only an adaptive immune response against tuberculosis but also triggers innate immune activation and memory that are indicative of trained immunity. We established personal immune profiles and chromatin accessibility maps over a 90-day time course of BCG vaccination in 323 individuals. Our analysis uncovered genetic and epigenetic predictors of baseline immunity and immune response. BCG vaccination enhanced the innate immune response specifically in individuals with a dormant immune state at baseline, rather than providing a general boost of innate immunity. This study advances our understanding of BCG's heterologous immune-stimulatory effects and trained immunity in humans. Furthermore, it highlights the value of epigenetic cell states for connecting immune function with genotype and the environment.

Metabolic Induction of Trained Immunity through the Mevalonate Pathway.

Innate immune cells can develop long-term memory after stimulation by microbial products during infections or vaccinations. Here, we report that metabolic signals can induce trained immunity. Pharmacological and genetic experiments reveal that activation of the cholesterol synthesis pathway, but not the synthesis of cholesterol itself, is essential for training of myeloid cells. Rather, the metabolite mevalonate is the mediator of training via activation of IGF1-R and mTOR and subsequent histone modifications in inflammatory pathways. Statins, which block mevalonate generation, prevent trained immunity induction. Furthermore, monocytes of patients with hyper immunoglobulin D syndrome (HIDS), who are mevalonate kinase deficient and accumulate mevalonate, have a constitutive trained immunity phenotype at both immunological and epigenetic levels, which could explain the attacks of sterile inflammation that these patients experience. Unraveling the role of mevalonate in trained immunity contributes to our understanding of the pathophysiology of HIDS and identifies novel therapeutic targets for clinical conditions with excessive activation of trained immunity.

Western Diet Triggers NLRP3-Dependent Innate Immune Reprogramming.

Long-term epigenetic reprogramming of innate immune cells in response to microbes, also termed "trained immunity," causes prolonged altered cellular functionality to protect from secondary infections. Here, we investigated whether sterile triggers of inflammation induce trained immunity and thereby influence innate immune responses. Western diet (WD) feeding of Ldlr-/- mice induced systemic inflammation, which was undetectable in serum soon after mice were shifted back to a chow diet (CD). In contrast, myeloid cell responses toward innate stimuli remained broadly augmented. WD-induced transcriptomic and epigenomic reprogramming of myeloid progenitor cells led to increased proliferation and enhanced innate immune responses. Quantitative trait locus (QTL) analysis in human monocytes trained with oxidized low-density lipoprotein (oxLDL) and stimulated with lipopolysaccharide (LPS) suggested inflammasome-mediated trained immunity. Consistently, Nlrp3-/-/Ldlr-/- mice lacked WD-induced systemic inflammation, myeloid progenitor proliferation, and reprogramming. Hence, NLRP3 mediates trained immunity following WD and could thereby mediate the potentially deleterious effects of trained immunity in inflammatory diseases.

CRELD1 modulates homeostasis of the immune system in mice and humans.

CRELD1 is a pivotal factor for heart development, the function of which is unknown in adult life. We here provide evidence that CRELD1 is an important gatekeeper of immune system homeostasis. Exploiting expression variance in large human cohorts contrasting individuals with the lowest and highest CRELD1 expression levels revealed strong phenotypic, functional and transcriptional differences, including reduced CD4+ T cell numbers. These findings were validated in T cell-specific Creld1-deficient mice. Loss of Creld1 was associated with simultaneous overactivation and increased apoptosis, resulting in a net loss of T cells with age. Creld1 was transcriptionally and functionally linked to Wnt signaling. Collectively, gene expression variance in large human cohorts combined with murine genetic models, transcriptomics and functional testing defines CRELD1 as an important modulator of immune homeostasis.

An Expanded View of Complex Traits: From Polygenic to Omnigenic.

A central goal of genetics is to understand the links between genetic variation and disease. Intuitively, one might expect disease-causing variants to cluster into key pathways that drive disease etiology. But for complex traits, association signals tend to be spread across most of the genome-including near many genes without an obvious connection to disease. We propose that gene regulatory networks are sufficiently interconnected such that all genes expressed in disease-relevant cells are liable to affect the functions of core disease-related genes and that most heritability can be explained by effects on genes outside core pathways. We refer to this hypothesis as an "omnigenic" model.

TMCO1 Is an ER Ca(2+) Load-Activated Ca(2+) Channel.

Maintaining homeostasis of Ca(2+) stores in the endoplasmic reticulum (ER) is crucial for proper Ca(2+) signaling and key cellular functions. The Ca(2+)-release-activated Ca(2+) (CRAC) channel is responsible for Ca(2+) influx and refilling after store depletion, but how cells cope with excess Ca(2+) when ER stores are overloaded is unclear. We show that TMCO1 is an ER transmembrane protein that actively prevents Ca(2+) stores from overfilling, acting as what we term a "Ca(2+) load-activated Ca(2+) channel" or "CLAC" channel. TMCO1 undergoes reversible homotetramerization in response to ER Ca(2+) overloading and disassembly upon Ca(2+) depletion and forms a Ca(2+)-selective ion channel on giant liposomes. TMCO1 knockout mice reproduce the main clinical features of human cerebrofaciothoracic (CFT) dysplasia spectrum, a developmental disorder linked to TMCO1 dysfunction, and exhibit severe mishandling of ER Ca(2+) in cells. Our findings indicate that TMCO1 provides a protective mechanism to prevent overfilling of ER stores with Ca(2+) ions.

A Functional Genomics Approach to Understand Variation in Cytokine Production in Humans.

As part of the Human Functional Genomics Project, which aims to understand the factors that determine the variability of immune responses, we investigated genetic variants affecting cytokine production in response to ex vivo stimulation in two independent cohorts of 500 and 200 healthy individuals. We demonstrate a strong impact of genetic heritability on cytokine production capacity after challenge with bacterial, fungal, viral, and non-microbial stimuli. In addition to 17 novel genome-wide significant cytokine QTLs (cQTLs), our study provides a comprehensive picture of the genetic variants that influence six different cytokines in whole blood, blood mononuclear cells, and macrophages. Important biological pathways that contain cytokine QTLs map to pattern recognition receptors (TLR1-6-10 cluster), cytokine and complement inhibitors, and the kallikrein system. The cytokine QTLs show enrichment for monocyte-specific enhancers, are more often located in regions under positive selection, and are significantly enriched among SNPs associated with infections and immune-mediated diseases. PAPERCLIP.

Host and Environmental Factors Influencing Individual Human Cytokine Responses.

Differences in susceptibility to immune-mediated diseases are determined by variability in immune responses. In three studies within the Human Functional Genomics Project, we assessed the effect of environmental and non-genetic host factors of the genetic make-up of the host and of the intestinal microbiome on the cytokine responses in humans. We analyzed the association of these factors with circulating mediators and with six cytokines after stimulation with 19 bacterial, fungal, viral, and non-microbial metabolic stimuli in 534 healthy subjects. In this first study, we show a strong impact of non-genetic host factors (e.g., age and gender) on cytokine production and circulating mediators. Additionally, annual seasonality is found to be an important environmental factor influencing cytokine production. Alpha-1-antitrypsin concentrations partially mediate the seasonality of cytokine responses, whereas the effect of vitamin D levels is limited. The complete dataset has been made publicly available as a comprehensive resource for future studies. PAPERCLIP.