Search

Single-cell quantification of RNAs is important for understanding cellular heterogeneity and gene regulation, yet current approaches suffer from low sensitivity for individual transcripts, limiting their utility for many applications. Here we present Hybridization of Probes to RNA for sequencing (HyPR-seq), a method to sensitively quantify the expression of hundreds of chosen genes in single cells. HyPR-seq involves hybridizing DNA probes to RNA, distributing cells into nanoliter droplets, amplifying the probes with PCR, and sequencing the amplicons to quantify the expression of chosen genes. HyPR-seq achieves high sensitivity for individual transcripts, detects nonpolyadenylated and low-abundance transcripts, and can profile more than 100,000 single cells. We demonstrate how HyPR-seq can profile the effects of CRISPR perturbations in pooled screens, detect time-resolved changes in gene expression via measurements of gene introns, and detect rare transcripts and quantify cell-type frequencies in tissue using low-abundance marker genes. By directing sequencing power to genes of interest and sensitively quantifying individual transcripts, HyPR-seq reduces costs by up to 100-fold compared to whole-transcriptome single-cell RNA-sequencing, making HyPR-seq a powerful method for targeted RNA profiling in single cells.

Characterizing the interactions that SARS-CoV-2 viral RNAs make with host cell proteins during infection can improve our understanding of viral RNA functions and the host innate immune response. Using RNA antisense purification and mass spectrometry, we identified up to 104 human proteins that directly and specifically bind to SARS-CoV-2 RNAs in infected human cells. We integrated the SARS-CoV-2 RNA interactome with changes in proteome abundance induced by viral infection and linked interactome proteins to cellular pathways relevant to SARS-CoV-2 infections. We demonstrated by genetic perturbation that cellular nucleic acid-binding protein (CNBP) and La-related protein 1 (LARP1), two of the most strongly enriched viral RNA binders, restrict SARS-CoV-2 replication in infected cells and provide a global map of their direct RNA contact sites. Pharmacological inhibition of three other RNA interactome members, PPIA, ATP1A1, and the ARP2/3 complex, reduced viral replication in two human cell lines. The identification of host dependency factors and defence strategies as presented in this work will improve the design of targeted therapeutics against SARS-CoV-2., Interactions between SARS-CoV-2 viral RNAs and host cell proteins during infection are evaluated to improve our understanding of viral RNA functions and the host innate immune response.

Structural variants in more than 17,000 human genomes are mapped and characterized using whole-genome sequencing, showing how this type of variation contributes to rare deleterious coding and noncoding alleles. A key goal of whole-genome sequencing for studies of human genetics is to interrogate all forms of variation, including single-nucleotide variants, small insertion or deletion (indel) variants and structural variants. However, tools and resources for the study of structural variants have lagged behind those for smaller variants. Here we used a scalable pipeline(1)to map and characterize structural variants in 17,795 deeply sequenced human genomes. We publicly release site-frequency data to create the largest, to our knowledge, whole-genome-sequencing-based structural variant resource so far. On average, individuals carry 2.9 rare structural variants that alter coding regions; these variants affect the dosage or structure of 4.2 genes and account for 4.0-11.2% of rare high-impact coding alleles. Using a computational model, we estimate that structural variants account for 17.2% of rare alleles genome-wide, with predicted deleterious effects that are equivalent to loss-of-function coding alleles; approximately 90% of such structural variants are noncoding deletions (mean 19.1 per genome). We report 158,991 ultra-rare structural variants and show that 2% of individuals carry ultra-rare megabase-scale structural variants, nearly half of which are balanced or complex rearrangements. Finally, we infer the dosage sensitivity of genes and noncoding elements, and reveal trends that relate to element class and conservation. This work will help to guide the analysis and interpretation of structural variants in the era of whole-genome sequencing.

Genome-wide association studies have associated thousands of genetic variants with complex traits and diseases, but pinpointing the causal variant(s) among those in tight linkage disequilibrium with each associated variant remains a major challenge. Here, we use seven experimental assays to characterize all common variants at the multiple disease-associated TNFAIP3 locus in five disease-relevant immune cell lines, based on a set of features related to regulatory potential. Trait/disease-associated variants are enriched among SNPs prioritized based on either: (1) residing within CRISPRi-sensitive regulatory regions, or (2) localizing in a chromatin accessible region while displaying allele-specific reporter activity. Of the 15 trait/disease-associated haplotypes at TNFAIP3, 9 have at least one variant meeting one or both of these criteria, 5 of which are further supported by genetic fine-mapping. Our work provides a comprehensive strategy to characterize genetic variation at important disease-associated loci, and aids in the effort to identify trait causal genetic variants., While genome-wide association studies have yielded thousands of trait-associated loci, identifying causal variants remains challenging. Here, the authors perform seven genomics assays in various cell types to prioritize genetic variants in the TNFAIP3 locus, and report high-priority variants within disease-associated haplotypes.

Gene regulation in the human genome is controlled by distal enhancers that activate specific nearby promoters. One model for the specificity of enhancer-promoter regulation is that different promoters might have sequence-encoded preferences for distinct classes of enhancers, for example mediated by interacting sets of transcription factors or cofactors. This “biochemical compatibility” model has been supported by observations at individual human promoters and by genome-wide measurements in Drosophila. However, the degree to which human enhancers and promoters are intrinsically compatible or specific has not been systematically measured, and how their activities combine to control RNA expression remains unclear. To address these questions, we designed a high-throughput reporter assay called enhancer x promoter (ExP) STARR-seq and applied it to examine the combinatorial compatibilities of 1,000 enhancer and 1,000 promoter sequences in human K562 cells. We identify a simple logic for enhancer-promoter compatibility – virtually all enhancers activated all promoters by similar amounts, and intrinsic enhancer and promoter activities combine multiplicatively to determine RNA output (R2=0.82). In addition, two classes of enhancers and promoters showed subtle preferential effects. Promoters of housekeeping genes contained built-in activating sequences, corresponding to motifs for factors such as GABPA and YY1, that correlated with both stronger autonomous promoter activity and enhancer activity, and weaker responsiveness to distal enhancers. Promoters of context-specific genes lacked these motifs and showed stronger responsiveness to enhancers. Together, this systematic assessment of enhancer-promoter compatibility suggests a multiplicative model tuned by enhancer and promoter class to control gene transcription in the human genome.

SUMMARYGenetic variation at the SLC16A11 locus contributes to the disproportionate impact of type 2 diabetes (T2D) on Latino populations. We recently demonstrated that T2D risk variants reduce SLC16A11 liver expression and function of MCT11, the monocarboxylate transporter encoded by the SLC16A11 gene. Here, we show that SLC16A11 expression within the liver is primarily localized to the low oxygen pericentral region, and that T2D risk variants disrupt oxygen-regulated SLC16A11 expression in human hepatocytes. Under physiologic oxygen conditions, MCT11 deficiency alters hepatocyte glucose metabolism, resulting in elevated intracellular lactate and a metabolic shift toward triacylglycerol accumulation. We also demonstrate an impact of Mct11 deficiency on glucose and lipid metabolism in Slc16a11 knockout mice, which display physiological changes that are observed in individuals with T2D. Our findings provide mechanistic insight into how SLC16A11 disruption impacts hepatic energy metabolism and T2D risk, and highlight MCT11-mediated regulation of lactate levels as a potential therapeutic target.

New antibiotics are needed to combat rising levels of resistance, with new Mycobacterium tuberculosis (Mtb) drugs having the highest priority. However, conventional whole-cell and biochemical antibiotic screens have failed. Here we develop a strategy termed PROSPECT (primary screening of strains to prioritize expanded chemistry and targets), in which we screen compounds against pools of strains depleted of essential bacterial targets. We engineered strains that target 474 essential Mtb genes and screened pools of 100–150 strains against activity-enriched and unbiased compound libraries, probing more than 8.5 million chemical–genetic interactions. Primary screens identified over tenfold more hits than screening wild-type Mtb alone, with chemical–genetic interactions providing immediate, direct target insights. We identified over 40 compounds that target DNA gyrase, the cell wall, tryptophan, folate biosynthesis and RNA polymerase, as well as inhibitors that target EfpA. Chemical optimization yielded EfpA inhibitors with potent wild-type activity, thus demonstrating the ability of PROSPECT to yield inhibitors against targets that would have eluded conventional drug discovery. A high-throughput chemical–genetic screening approach for the discovery of targets and chemicals to treat Mycobacterium tuberculosis yields tenfold more hit compounds than conventional whole-cell screening methods.

Genome-wide association studies (GWAS) have identified thousands of noncoding loci that are associated with human diseases and complex traits, each of which could reveal insights into the mechanisms of disease1. Many of the underlying causal variants may affect enhancers2,3, but we lack accurate maps of enhancers and their target genes to interpret such variants. We recently developed the activity-by-contact (ABC) model to predict which enhancers regulate which genes and validated the model using CRISPR perturbations in several cell types4. Here we apply this ABC model to create enhancer–gene maps in 131 human cell types and tissues, and use these maps to interpret the functions of GWAS variants. Across 72 diseases and complex traits, ABC links 5,036 GWAS signals to 2,249 unique genes, including a class of 577 genes that appear to influence multiple phenotypes through variants in enhancers that act in different cell types. In inflammatory bowel disease (IBD), causal variants are enriched in predicted enhancers by more than 20-fold in particular cell types such as dendritic cells, and ABC achieves higher precision than other regulatory methods at connecting noncoding variants to target genes. These variant-to-function maps reveal an enhancer that contains an IBD risk variant and that regulates the expression of PPIF to alter the membrane potential of mitochondria in macrophages. Our study reveals principles of genome regulation, identifies genes that affect IBD and provides a resource and generalizable strategy to connect risk variants of common diseases to their molecular and cellular functions. Mapping enhancer regulation across human cell types and tissues illuminates genome function and provides a resource to connect risk variants for common diseases to their molecular and cellular functions.

EBioMedicine 72, 103588 (2021). doi:10.1016/j.ebiom.2021.103588, Published by Elsevier, Amsterdam [u.a.]

Summary Skeletal and glycemic traits have shared etiology, but the underlying genetic factors remain largely unknown. To identify genetic loci that may have pleiotropic effects, we studied Genome-wide association studies (GWASs) for bone mineral density and glycemic traits and identified a bivariate risk locus at 3q21. Using sequence and epigenetic modeling, we prioritized an adenylate cyclase 5 (ADCY5) intronic causal variant, rs56371916. This SNP changes the binding affinity of SREBP1 and leads to differential ADCY5 gene expression, altering the chromatin landscape from poised to repressed. These alterations result in bone- and type 2 diabetes-relevant cell-autonomous changes in lipid metabolism in osteoblasts and adipocytes. We validated our findings by directly manipulating the regulator SREBP1, the target gene ADCY5, and the variant rs56371916, which together imply a novel link between fatty acid oxidation and osteoblast differentiation. Our work, by systematic functional dissection of pleiotropic GWAS loci, represents a framework to uncover biological mechanisms affecting pleiotropic traits., Graphical abstract, Highlights • 3q21 variants are associated with type 2 diabetes and increased bone mineral density • 3q21 variants affect chromatin accessibility in mesenchymal cells • rs56371916 is a 3q21 causal variant and ADCY5 its target in adipocytes and osteoblasts • ADCY5 and rs56371916 affect lipid oxidation processes in adipocytes and osteoblasts, Nasa Sinnott-Armstrong and colleagues identify a pleiotropic risk locus at 3q21 that is associated with type 2 diabetes (T2D) and greater bone mineral density (BMD) and its associated cell-autonomous mechanisms in adipocytes and osteoblasts. Together, these findings provide a possible explanation for the perplexing finding that individuals with T2D have higher BMD but greater susceptibility to bone fracture.

Genome-wide association studies (GWAS) are a valuable tool for understanding the biology of complex traits, but the associations found rarely point directly to causal genes. Here, we introduce a new method to identify the causal genes by integrating GWAS summary statistics with gene expression, biological pathway, and predicted protein-protein interaction data. We further propose an approach that effectively leverages both polygenic and locus-specific genetic signals by combining results across multiple gene prioritization methods, increasing confidence in prioritized genes. Using a large set of gold standard genes to evaluate our approach, we prioritize 8,402 unique gene-trait pairs with greater than 75% estimated precision across 113 complex traits and diseases, including known genes such as SORT1 for LDL cholesterol, SMIM1 for red blood cell count, and DRD2 for schizophrenia, as well as novel genes such as TTC39B for cholelithiasis. Our results demonstrate that a polygenic approach is a powerful tool for gene prioritization and, in combination with locus-specific signal, improves upon existing methods.

Genome-wide association studies have now identified tens of thousands of noncoding loci associated with human diseases and complex traits, each of which could reveal insights into biological mechanisms of disease. Many of the underlying causal variants are thought to affect enhancers, but we have lacked genome-wide maps of enhancer-gene regulation to interpret such variants. We previously developed the Activity-by-Contact (ABC) Model to predict enhancer-gene connections and demonstrated that it can accurately predict the results of CRISPR perturbations across several cell types. Here, we apply this ABC Model to create enhancer-gene maps in 131 cell types and tissues, and use these maps to interpret the functions of fine-mapped GWAS variants. For inflammatory bowel disease (IBD), causal variants are >20-fold enriched in enhancers in particular cell types, and ABC outperforms other regulatory methods at connecting noncoding variants to target genes. Across 72 diseases and complex traits, ABC links 5,036 GWAS signals to 2,249 unique genes, including a class of 577 genes that appear to influence multiple phenotypes via variants in enhancers that act in different cell types. Guided by these variant-to-function maps, we show that an enhancer containing an IBD risk variant regulates the expression ofPPIFto tune mitochondrial membrane potential. Together, our study reveals insights into principles of genome regulation, illuminates mechanisms that influence IBD, and demonstrates a generalizable strategy to connect common disease risk variants to their molecular and cellular functions.

Genetic variation can predispose to disease both through (i) monogenic risk variants that disrupt a physiologic pathway with large effect on disease and (ii) polygenic risk that involves many variants of small effect in different pathways. Few studies have explored the interplay between monogenic and polygenic risk. Here, we study 80,928 individuals to examine whether polygenic background can modify penetrance of disease in tier 1 genomic conditions — familial hypercholesterolemia, hereditary breast and ovarian cancer, and Lynch syndrome. Among carriers of a monogenic risk variant, we estimate substantial gradients in disease risk based on polygenic background — the probability of disease by age 75 years ranged from 17% to 78% for coronary artery disease, 13% to 76% for breast cancer, and 11% to 80% for colon cancer. We propose that accounting for polygenic background is likely to increase accuracy of risk estimation for individuals who inherit a monogenic risk variant., Genetic variation predisposes to disease via monogenic and polygenic risk variants. Here, the authors assess the interplay between these types of variation on disease penetrance in 80,928 individuals. In carriers of monogenic variants, they show that disease risk is a gradient influenced by polygenic background.

SARS-CoV-2 infections pose a global threat to human health and an unprecedented research challenge. Among the most urgent tasks is obtaining a detailed understanding of the molecular interactions that facilitate viral replication or contribute to host defense mechanisms in infected cells. While SARS-CoV-2 co-opts cellular factors for viral translation and genome replication, a comprehensive map of the host cell proteome in direct contact with viral RNA has not been elucidated. Here, we use RNA antisense purification and mass spectrometry (RAP-MS) to obtain an unbiased and quantitative picture of the human proteome that directly binds the SARS-CoV-2 RNA in infected human cells. We discover known host factors required for coronavirus replication, regulators of RNA metabolism and host defense pathways, along with dozens of potential drug targets among direct SARS-CoV-2 binders. We further integrate the SARS-CoV-2 RNA interactome with proteome dynamics induced by viral infection, linking interactome proteins to the emerging biology of SARS-CoV-2 infections. Validating RAP-MS, we show that CNBP, a regulator of proinflammatory cytokines, directly engages the SARS-CoV-2 RNA. Supporting the functional relevance of identified interactors, we show that the interferon-induced protein RYDEN suppresses SARS-CoV-2 ribosomal frameshifting and demonstrate that inhibition of SARS-CoV-2-bound proteins is sufficient to manipulate viral replication. The SARS-CoV-2 RNA interactome provides an unprecedented molecular perspective on SARS-CoV-2 infections and enables the systematic dissection of host dependency factors and host defense strategies, a crucial prerequisite for designing novel therapeutic strategies.

Single-cell quantification of RNAs is important for understanding cellular heterogeneity and gene regulation, yet current approaches suffer from low sensitivity for individual transcripts, limiting their utility for many applications. Here we present Hybridization of Probes to RNA for sequencing (HyPR-seq), a method to sensitively quantify the expression of up to 100 chosen genes in single cells. HyPR-seq involves hybridizing DNA probes to RNA, distributing cells into nanoliter droplets, amplifying the probes with PCR, and sequencing the amplicons to quantify the expression of chosen genes. HyPR-seq achieves high sensitivity for individual transcripts, detects nonpolyadenylated and low-abundance transcripts, and can profile more than 100,000 single cells. We demonstrate how HyPR-seq can profile the effects of CRISPR perturbations in pooled screens, detect time-resolved changes in gene expression via measurements of gene introns, and detect rare transcripts and quantify cell type frequencies in tissue using low-abundance marker genes. By directing sequencing power to genes of interest and sensitively quantifying individual transcripts, HyPR-seq reduces costs by up to 100-fold compared to whole-transcriptome scRNA-seq, making HyPR-seq a powerful method for targeted RNA profiling in single cells.

Structural variants (SVs) rearrange large segments of DNA1 and can have profound consequences in evolution and human disease2,3. As national biobanks, disease-association studies, and clinical genetic testing have grown increasingly reliant on genome sequencing, population references such as the Genome Aggregation Database (gnomAD)4 have become integral in the interpretation of single-nucleotide variants (SNVs)5. However, there are no reference maps of SVs from high-coverage genome sequencing comparable to those for SNVs. Here we present a reference of sequence-resolved SVs constructed from 14,891 genomes across diverse global populations (54% non-European) in gnomAD. We discovered a rich and complex landscape of 433,371 SVs, from which we estimate that SVs are responsible for 25–29% of all rare protein-truncating events per genome. We found strong correlations between natural selection against damaging SNVs and rare SVs that disrupt or duplicate protein-coding sequence, which suggests that genes that are highly intolerant to loss-of-function are also sensitive to increased dosage6. We also uncovered modest selection against noncoding SVs in cis-regulatory elements, although selection against protein-truncating SVs was stronger than all noncoding effects. Finally, we identified very large (over one megabase), rare SVs in 3.9% of samples, and estimate that 0.13% of individuals may carry an SV that meets the existing criteria for clinically important incidental findings7. This SV resource is freely distributed via the gnomAD browser8 and will have broad utility in population genetics, disease-association studies, and diagnostic screening., A large empirical assessment of sequence-resolved structural variants from 14,891 genomes across diverse global populations in the Genome Aggregation Database (gnomAD) provides a reference map for disease-association studies, population genetics, and diagnostic screening.

Delivery mode is the variable with the greatest influence on the infant gut microbiome composition in the first few months of life. Children born by Cesarean section (C-section) lack species from the Bacteroides genus in their gut microbial community, and this difference can be detectable until 6-18 months of age. One hypothesis is that these differences stem from lack of exposure to the maternal vaginal microbiome, as children born by C-section do not pass through the birth canal; however, Bacteroides species are not common members of the vaginal microbiome, thus this explanation seems inadequate. Here, we set out to re-evaluate this hypothesis by collecting rectal and vaginal samples before delivery from 73 mothers with paired stool from their infants in the first two weeks of life. We compared microbial profiles of infants born by planned, pre-labor C-section to those born by emergent, post-labor surgery (where the child was in the birth canal, but eventually delivered through an abdominal incision), and found no significant differences in the microbiome between these two groups. Both groups showed the characteristic signature lack of Bacteroides species, despite their difference in exposure to the birth canal. Surprisingly, this signature was only evident in samples from week two of life, but not in the first week. Children born by C-section often had high abundance of Bacteroides in their first few days of life, but these were not stable colonizers of the infant gut, as they were not detectable by week two. Finally, we used metagenomic sequencing to compare microbial strains in maternal vaginal and rectal samples and samples from their infants; we found evidence for mother-to-child transmission of rectal rather than vaginal strains. These results challenge birth canal exposure as the dominant factor in infant gut microbiome establishment and implicate colonization efficiency rather than exposure as a dictating factor of the newborn gut microbiome composition.

Summary Mode of delivery strongly influences the early infant gut microbiome. Children born by cesarean section (C-section) lack Bacteroides species until 6–18 months of age. One hypothesis is that these differences stem from lack of exposure to the maternal vaginal microbiome. Here, we re-evaluate this hypothesis by comparing the microbial profiles of 75 infants born vaginally or by planned versus emergent C-section. Multiple children born by C-section have a high abundance of Bacteroides in their first few days of life, but at 2 weeks, both C-section groups lack Bacteroides (primarily according to 16S sequencing), despite their difference in exposure to the birth canal. Finally, a comparison of microbial strain profiles between infants and maternal vaginal or rectal samples finds evidence for mother-to-child transmission of rectal rather than vaginal strains. These results suggest differences in colonization stability as an important factor in infant gut microbiome composition rather than birth canal exposure., Graphical Abstract, Highlights Week 1 gut microbiota does not differ between infants born vaginally versus C-section Week 2 gut microbiota of C-section infants lacks Bacteroides Microbiota of infants born by C-section after labor resembles scheduled C-section Bacterial strains in infants match maternal rectal rather than vaginal strains, Mitchell et al. compare early-life infant gut microbiota by delivery mode, suggesting early colonization by Bacteroides regardless of delivery mode, but loss of Bacteroides by 2 weeks in C-section-delivered infants, whether or not exposed to the vagina in labor. Infant strains matched maternal rectal rather vaginal strains.

The Zoonomia Project is investigating the genomics of shared and specialized traits in eutherian mammals. Here we provide genome assemblies for 131 species, of which all but 9 are previously uncharacterized, and describe a whole-genome alignment of 240 species of considerable phylogenetic diversity, comprising representatives from more than 80% of mammalian families. We find that regions of reduced genetic diversity are more abundant in species at a high risk of extinction, discern signals of evolutionary selection at high resolution and provide insights from individual reference genomes. By prioritizing phylogenetic diversity and making data available quickly and without restriction, the Zoonomia Project aims to support biological discovery, medical research and the conservation of biodiversity., This project was funded by NIH NHGRI R01HG008742 (E.K.K., B.B., D.P.G., R.S., J.T.-M., J.J., H.J.N., B.P. and J. Armstrong), Swedish Research Council Distinguished Professor Award (K.L.-T., V.D.M., E.M. and J.R.S.M.), Swedish Research Council grant 2018-05973 (K.L.-T.), Knut and Alice Wallenberg Foundation (K.L.-T., V.D.M., E.M. and J.R.S.M.), Uppsala University (K.L.-T., V.D.M., E.M., J.R.S.M., J.J., J. Alfoldi and L.G.), Broad Institute Next10 (L.G.), Gladstone Institutes (K.S.P.), NIH NHGRI 5R01HG002939 (A.F.A.S. and R.H.), NIH NHGRI 5U24HG010136 (A.F.A.S. and R.H.), NIH NHGRI 5R01HG010485 (B.P. and M.D.), NIH NHGRI 2U41HG007234 (B.P., M.D. and J. Armstrong), NIH NIA 5PO1AG047200 (V.N.G.), NIH NIA 1UH2AG064706 (V.N.G.), BFU2017-86471-P MINECO/FEDER, UE (T.M.-B.), Secretaria d’Universitats i Recerca and CERCA Programme del Departament d’Economia i Coneixement de la Generalitat de Catalunya GRC 2017 SGR 880 (T.M.-B.), Howard Hughes International Early Career (T.M.-B.), European Research Council Horizon 2020 no. 864203 (T.M.-B.), Obra Social ‘La Caixa’ (T.M.-B.), BBSRC BBS/E/T/000PR9818, BBS/E/T/ 000PR9783 (W.H. and W.N.), BBSRC Core Strategic Programme Grant BB/P016774/1 (W.H., W.N. and F.D.), Sir Henry Dale Fellowship 200517/Z/16/Z jointly funded by the Wellcome Trust and the Royal Society (N.R.C.), FJCI-2016-29558 MICINN (D.J.), Prince Albert II Foundation of Monaco and Canada, Global Genome Initiative, Smithsonian Institution (M.N.), European Research Council Research Grant ERC-2012-StG311000 (E.C.T.), Irish Research Council Laureate Award (E.C.T.), UK Medical Research Council MR/P026028/1 (W.H. and W.N.), National Science Foundation DEB-1457735 (M.S.S.), National Science Foundation DEB-1753760 (W.J.M.), National Science Foundation IOS-2029774 (E.K.K. and D.P.G.), Robert and Rosabel Osborne Endowment (H.A.L. and J.D.), Swedish Research Council, FORMAS 221-2012-1531 (J.R.S.M.), NSF RoL: FELS: EAGER: DEB 1838283 (D.A.R.) and Academy of Finland grant to Center of Excellence in Tumor Genetics Research no. 312042 (T.K. and J.T.).

The Alzheimer’s Disease Sequencing Project (ADSP) undertook whole exome sequencing in 5,740 late-onset Alzheimer disease (AD) cases and 5,096 cognitively normal controls primarily of European ancestry (EA), among whom 218 cases and 177 controls were Caribbean Hispanic (CH). An age-, sex- and APOE based risk score and family history were used to select cases most likely to harbor novel AD risk variants and controls least likely to develop AD by age 85 years. We tested ~1.5 million single nucleotide variants (SNVs) and 50,000 insertion-deletion polymorphisms (indels) for association to AD, using multiple models considering individual variants as well as gene-based tests aggregating rare, predicted functional, and loss of function variants. Sixteen single variants and 19 genes that met criteria for significant or suggestive associations after multiple-testing correction were evaluated for replication in four independent samples; three with whole exome sequencing (2,778 cases, 7,262 controls) and one with genome-wide genotyping imputed to the Haplotype Reference Consortium panel (9,343 cases, 11,527 controls). The top findings in the discovery sample were also followed-up in the ADSP whole-genome sequenced family-based dataset (197 members of 42 EA families and 501 members of 157 CH families). We identified novel and predicted functional genetic variants in genes previously associated with AD. We also detected associations in three novel genes: IGHG3 (p = 9.8 × 10−7), an immunoglobulin gene whose antibodies interact with β-amyloid, a long non-coding RNA AC099552.4 (p = 1.2 × 10−7), and a zinc-finger protein ZNF655 (gene-based p = 5.0 × 10−6). The latter two suggest an important role for transcriptional regulation in AD pathogenesis.

Large-scale deep-coverage whole-genome sequencing (WGS) is now feasible and offers potential advantages for locus discovery. We perform WGS in 16,324 participants from four ancestries at mean depth >29X and analyze genotypes with four quantitative traits—plasma total cholesterol, low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol, and triglycerides. Common variant association yields known loci except for few variants previously poorly imputed. Rare coding variant association yields known Mendelian dyslipidemia genes but rare non-coding variant association detects no signals. A high 2M-SNP LDL-C polygenic score (top 5th percentile) confers similar effect size to a monogenic mutation (~30 mg/dl higher for each); however, among those with severe hypercholesterolemia, 23% have a high polygenic score and only 2% carry a monogenic mutation. At these sample sizes and for these phenotypes, the incremental value of WGS for discovery is limited but WGS permits simultaneous assessment of monogenic and polygenic models to severe hypercholesterolemia., Common genetic variants associated with plasma lipids have been extensively studied for a better understanding of common diseases. Here, the authors use whole-genome sequencing of 16,324 individuals to analyze rare variant associations and to determine their monogenic and polygenic contribution to lipid traits.

RNA plays important and diverse roles in biology, but molecular tools to manipulate and measure RNA are limited. For example, RNA interference (RNAi)1-3 can efficiently knockdown RNAs, but it is prone to off-target effects4, and visualizing RNAs typically relies on the introduction of exogenous tags5. Here, we demonstrate that the class 2 type VI6,7 RNA-guided RNA-targeting CRISPR-Cas effector Cas13a8 (previously known as C2c2) can be engineered for mammalian cell RNA knockdown and binding. After initial screening of fifteen orthologs in E. coli, we identified Cas13a from Leptotrichia wadei (LwaCas13a) as the most effective. LwaCas13a can be heterologously expressed in mammalian and plant cells for targeted knockdown of either reporter or endogenous transcripts. We demonstrate that LwaCas13a is capable of providing comparable levels of knockdown as RNAi, but with dramatically improved specificity. Moreover, catalytically inactive LwaCas13a maintains targeted RNA binding, allowing for programmable tracking of transcripts in live cells. Our results establish CRISPR-Cas13a as a flexible platform for RNA targeting with wide applicability for studying RNA in mammalian cells.

Hi-C for mosquito genomes Most genomes sequenced today are determined through the generation of short sequenced bits of DNA that are computationally pieced together like a jigsaw puzzle. This has resulted in the need for funds and additional data to fill in gaps in order to fully assemble the many chromosomes that make up a eukaryotic genome. Dudchenko et al. used the Hi-C method, which measures the distance between contact points within and between chromosomes for scaffold validation, together with correction and ordering to more completely determine the arrangement of short sequencing reads for genome mapping. They validated their approach through the de novo generation of a complete human genome. A comparative analysis of mosquito genomes was made possible by improving the Culex quinquefasciatus genome assembly and generating the genome of Aedes aegypti , the vector of Zika virus. Science , this issue p. 92

Age is the dominant risk factor for most chronic human diseases; yet the mechanisms by which aging confers this risk are largely unknown.1 Recently, the age-related acquisition of somatic mutations in regenerating hematopoietic stem cell populations was associated with both hematologic cancer incidence2–4 and coronary heart disease prevalence.5 Somatic mutations with leukemogenic potential may confer selective cellular advantages leading to clonal expansion, a phenomenon termed ‘Clonal Hematopoiesis of Indeterminate Potential’ (CHIP).6 Simultaneous germline and somatic whole genome sequence analysis now provides the opportunity to identify root causes of CHIP. Here, we analyze high-coverage whole genome sequences from 97,691 participants of diverse ancestries in the NHLBI TOPMed program and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid, and inflammatory traits specific to different CHIP genes. Association of a genome-wide set of germline genetic variants identified three genetic loci associated with CHIP status, including one locus at TET2 that was African ancestry specific. In silico-informed in vitro evaluation of the TET2 germline locus identified a causal variant that disrupts a TET2 distal enhancer. Aggregates of rare germline loss-of-function variants in CHEK2, a DNA damage repair gene, predisposed to CHIP acquisition. Overall, we observe that germline genetic variation altering hematopoietic stem cell function and the fidelity of DNA-damage repair increase the likelihood of somatic mutations leading to CHIP.

© 2020, The Author(s), under exclusive licence to Springer Nature Limited. Age is the dominant risk factor for most chronic human diseases, but the mechanisms through which ageing confers this risk are largely unknown1. The age-related acquisition of somatic mutations that lead to clonal expansion in regenerating haematopoietic stem cell populations has recently been associated with both haematological cancer2–4 and coronary heart disease5—this phenomenon is termed clonal haematopoiesis of indeterminate potential (CHIP)6. Simultaneous analyses of germline and somatic whole-genome sequences provide the opportunity to identify root causes of CHIP. Here we analyse high-coverage whole-genome sequences from 97,691 participants of diverse ancestries in the National Heart, Lung, and Blood Institute Trans-omics for Precision Medicine (TOPMed) programme, and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid and inflammatory traits that are specific to different CHIP driver genes. Association of a genome-wide set of germline genetic variants enabled the identification of three genetic loci associated with CHIP status, including one locus at TET2 that was specific to individuals of African ancestry. In silico-informed in vitro evaluation of the TET2 germline locus enabled the identification of a causal variant that disrupts a TET2 distal enhancer, resulting in increased self-renewal of haematopoietic stem cells. Overall, we observe that germline genetic variation shapes haematopoietic stem cell function, leading to CHIP through mechanisms that are specific to clonal haematopoiesis as well as shared mechanisms that lead to somatic mutations across tissues.

© 2020, The Author(s), under exclusive licence to Springer Nature America, Inc. Cancer genomes contain large numbers of somatic mutations but few of these mutations drive tumor development. Current approaches either identify driver genes on the basis of mutational recurrence or approximate the functional consequences of nonsynonymous mutations by using bioinformatic scores. Passenger mutations are enriched in characteristic nucleotide contexts, whereas driver mutations occur in functional positions, which are not necessarily surrounded by a particular nucleotide context. We observed that mutations in contexts that deviate from the characteristic contexts around passenger mutations provide a signal in favor of driver genes. We therefore developed a method that combines this feature with the signals traditionally used for driver-gene identification. We applied our method to whole-exome sequencing data from 11,873 tumor–normal pairs and identified 460 driver genes that clustered into 21 cancer-related pathways. Our study provides a resource of driver genes across 28 tumor types with additional driver genes identified according to mutations in unusual nucleotide contexts.

Severe obesity is a rapidly growing global health threat. Although often attributed to unhealthy lifestyle choices or environmental factors, obesity is known to be heritable and highly polygenic; the majority of inherited susceptibility is related to the cumulative effect of many common DNA variants. Here we derive and validate a new polygenic predictor comprised of 2.1 million common variants to quantify this susceptibility and test this predictor in more than 300,000 individuals ranging from middle age to birth. Among middle-aged adults, we observe a 13-kg gradient in weight and a 25-fold gradient in risk of severe obesity across polygenic score deciles. In a longitudinal birth cohort, we note minimal differences in birthweight across score deciles, but a significant gradient emerged in early childhood and reached 12 kg by 18 years of age. This new approach to quantify inherited susceptibility to obesity affords new opportunities for clinical prevention and mechanistic assessment. © 2019 Author(s), National Human Genome Research Institute (1K08HG0101), Wellcome Trust (202802/Z/16/Z), University of Bristol NIHR Biomedical Research Centre (S- BRC-1215-20011), National Human Genome Research Institute (HG008895), National Heart, Lung, and Blood Institute (NHLBI) HHSN268201300025C, National Heart, Lung, and Blood Institute (NHLBI) HHSN268201300026C, National Heart, Lung, and Blood Institute (NHLBI) HHSN268201300027C, National Heart, Lung, and Blood Institute (NHLBI) HHSN268201300028C, National Heart, Lung, and Blood Institute (NHLBI) HHSN268201300029C, National Heart, Lung, and Blood Institute (NHLBI) HHSN268200900041C, National Institute on Aging (AG0005), NHLBI (AG0005), National Human Genome Research Institute (U01-HG004729), National Human Genome Research Institute (U01-HG04424), National Human Genome Research Institute (U01-HG004446), Wellcome (102215/2/13/2)

Summary: Human genetic variants in SLC16A11 are associated with increased risk of type 2 diabetes (T2D). We previously identified two distinct mechanisms through which co-inherited T2D-risk coding and non-coding variants disrupt SLC16A11 expression and activity, thus implicating reduced SLC16A11 function as the disease-relevant direction of effect. In a recent publication, Zhao et al. (2019a) argue that human SLC16A11 coding variants confer gain of function, basing their conclusions on phenotypic changes observed following overexpression of mutant murine Slc16a11. However, data necessary to demonstrate gain-of-function activity are not reported. Furthermore, several fundamental flaws in their experimental system—including inaccurate modeling of the human variant haplotype and expression conditions that are not physiologically relevant—prevent conclusions about T2D-risk variant effects on human physiology. This Matters Arising paper is in response to Zhao et al. (2019a), published in Cell Reports. See also the response by Zhao et al. (2019b) in this issue of Cell Reports. : Hoch et al. discuss the analysis of coding variants of SLC16A11 in mice and humans in light of a recent Cell Reports publication. This Matters Arising paper is in response to Zhao et al. (2019a), published in Cell Reports. See also the response by Zhao et al. (2019b), published in this issue. Keywords: SLC16A11, type 2 diabetes, T2D

Natural killer (NK) cells distinguish cancer cells from healthy cells using an array of germline-encoded receptors that interact with ligands expressed on target cells. A balance of inhibitory and activating signals transduced by these receptors regulate NK cell function to provide anti-tumor immunity while maintaining self-tolerance. However, knowledge of the spectrum of factors regulating NK-cell-mediated cytotoxicity, including the contribution of specific ligands and regulatory mechanisms for their expression on tumor cells, remains incomplete. Here, we apply a genome-wide loss-of-function screen in tumor cells using CRISPR/Cas9 technology to identify the factors that promote NK-cell cytotoxicity towards tumor cells. We established the drivers of tumor-cell sensitivity towards NK-cell attack (TuSeNKA) screening approach using the chronic myeloid leukemia (CML) cell line, K562. Interestingly, we identified B7H6, the ligand for the activating NK cell receptor NKp30, as the single factor whose loss resulted in increased resistance of K562 cells towards NK cells. Our study shows that combination of CRISPR-based genetic screens with NK-cell cytotoxicity assays is a valuable tool for identifying functionally relevant NK cell-tumor cell interactions, paving the way for further investigations that unravel the complexity of signals that promote NK-cell recognition of transformed cells and develop therapies that target these modes of tumor-cell killing.

Mammalian genomes harbor millions of noncoding elements called enhancers that quantitatively regulate gene expression, but it remains unclear which enhancers regulate which genes. Here we describe an experimental approach, based on CRISPR interference, RNA FISH, and flow cytometry (CRISPRi-FlowFISH), to perturb enhancers in the genome, and apply it to test >3,000 potential regulatory enhancer-gene connections across multiple genomic loci. A simple equation based on a mechanistic model for enhancer function performed remarkably well at predicting the complex patterns of regulatory connections we observe in our CRISPR dataset. This Activity-by-Contact (ABC) model involves multiplying measures of enhancer activity and enhancer-promoter 3D contacts, and can predict enhancer-gene connections in a given cell type based on chromatin state maps. Together, CRISPRi-FlowFISH and the ABC model provide a systematic approach to map and predict which enhancers regulate which genes, and will help to interpret the functions of the thousands of disease risk variants in the noncoding genome.

Enhancer elements in the human genome control how genes are expressed in specific cell types and harbor thousands of genetic variants that influence risk for common diseases1-4. Yet, we still do not know how enhancers regulate specific genes, and we lack general rules to predict enhancer-gene connections across cell types5,6. We developed an experimental approach, CRISPRi-FlowFISH, to perturb enhancers in the genome, and we applied it to test >3,500 potential enhancer-gene connections for 30 genes. We found that a simple activity-by-contact model substantially outperformed previous methods at predicting the complex connections in our CRISPR dataset. This activity-by-contact model allows us to construct genome-wide maps of enhancer-gene connections in a given cell type, on the basis of chromatin state measurements. Together, CRISPRi-FlowFISH and the activity-by-contact model provide a systematic approach to map and predict which enhancers regulate which genes, and will help to interpret the functions of the thousands of disease risk variants in the noncoding genome.

Sequencing-based studies have identified novel risk genes for rare, severe epilepsies and revealed a role of rare deleterious variation in common epilepsies. To identify the shared and distinct ultra-rare genetic risk factors for rare and common epilepsies, we performed a whole-exome sequencing (WES) analysis of 9,170 epilepsy-affected individuals and 8,364 controls of European ancestry. We focused on three phenotypic groups; the rare but severe developmental and epileptic encephalopathies (DEE), and the commoner phenotypes of genetic generalized epilepsy (GGE) and non-acquired focal epilepsy (NAFE). We observed that compared to controls, individuals with any type of epilepsy carried an excess of ultra-rare, deleterious variants in constrained genes and in genes previously associated with epilepsy, with the strongest enrichment seen in DEE and the least in NAFE. Moreover, we found that inhibitory GABAA receptor genes were enriched for missense variants across all three classes of epilepsy, while no enrichment was seen in excitatory receptor genes. The larger gene groups for the GABAergic pathway or cation channels also showed a significant mutational burden in DEE and GGE. Although no single gene surpassed exome-wide significance among individuals with GGE or NAFE, highly constrained genes and genes encoding ion channels were among the top associations, including CACNA1G, EEF1A2, and GABRG2 for GGE and LGI1, TRIM3, and GABRG2 for NAFE. Our study confirms a convergence in the genetics of common and rare epilepsies associated with ultra-rare coding variation and highlights a ubiquitous role for GABAergic inhibition in epilepsy etiology in the largest epilepsy WES study to date.

A key goal of whole genome sequencing (WGS) for human genetics studies is to interrogate all forms of variation, including single nucleotide variants (SNV), small insertion/deletion (indel) variants and structural variants (SV). However, tools and resources for the study of SV have lagged behind those for smaller variants. Here, we used a cloud-based pipeline to map and characterize SV in 17,795 deeply sequenced human genomes from common disease trait mapping studies. We publicly release site-frequency information to create the largest WGS-based SV resource to date. On average, individuals carry 2.9 rare SVs that alter coding regions, which affect the dosage or structure of 4.2 genes and account for 4.0-11.2% of rare high-impact coding alleles. Based on a computational model, we estimate that SVs account for 17.2% of rare alleles genome-wide whose predicted deleterious effects are equivalent to loss-of-function (LoF) coding alleles; ~90% of such SVs are non-coding deletions (mean 19.1 per genome). We report 158,991 ultra-rare SVs and show that ~2% of individuals carry ultra-rare megabase-scale SVs, nearly half of which are balanced and/or complex rearrangements. Finally, we exploit this resource to infer the dosage sensitivity of genes and non-coding elements, revealing strong trends related to regulatory element class, conservation and cell-type specificity. This work will help guide SV analysis and interpretation in the era of WGS.

In this protocol, pooled sgRNA libraries targeting thousands of genes are computationally designed, generated using microarray-based synthesis techniques, and packaged into lentiviral particles. Target cells of interest are transduced with the lentiviral sgRNA pools to generate a collection of knockout mutants-via Cas9-mediated genomic cleavage-and screened for a phenotype of interest. The relative abundance of each mutant in the population can be monitored over time through high-throughput sequencing of the integrated sgRNA expression cassettes. Using this technique, we outline strategies for the identification of cancer driver genes and genes mediating drug response.

Many cancer genomes contain large numbers of somatic mutations, but few of these mutations drive tumor development. Current approaches to identify cancer driver genes are largely based on mutational recurrence, i.e. they search for genes with an increased number of nonsynonymous mutations relative to the local background mutation rate. Multiple studies have noted that the sensitivity of recurrence-based methods is limited in tumors with high background mutation rates, because passenger mutations dilute their statistical power. Here, we observe that passenger mutations tend to occur in characteristic nucleotide sequence contexts, while driver mutations follow a different distribution pattern determined by the location of functionally relevant genomic positions along the protein-coding sequence. To discover new cancer genes, we searched for genes with an excess of mutations in unusual nucleotide contexts that deviate from the characteristic context around passenger mutations. By applying this statistical framework to whole-exome sequencing data from 12,004 tumors, we discovered a long tail of novel candidate cancer genes with mutation frequencies as low as 1% and functional supporting evidence. Our results show that considering both the number and the nucleotide context around mutations helps identify novel cancer driver genes, particularly in tumors with high background mutation rates.

Diamond-Blackfan anemia (DBA) is a rare bone marrow failure disorder that affects 1 in 100,000 to 200,000 live births and has been associated with mutations in components of the ribosome. In order to characterize the genetic landscape of this genetically heterogeneous disorder, we recruited a cohort of 472 individuals with a clinical diagnosis of DBA and performed whole exome sequencing (WES). Overall, we identified rare and predicted damaging mutations in likely causal genes for 78% of individuals. The majority of mutations were singletons, absent from population databases, predicted to cause loss of function, and in one of 19 previously reported genes encoding for a diverse set of ribosomal proteins (RPs). Using WES exon coverage estimates, we were able to identify and validate 31 deletions in DBA associated genes. We also observed an enrichment for extended splice site mutations and validated the diverse effects of these mutations using RNA sequencing in patientderived cell lines. Leveraging the size of our cohort, we observed several robust genotype-phenotype associations with congenital abnormalities and treatment outcomes. In addition to comprehensively identifying mutations in known genes, we further identified rare mutations in 7 previously unreported RP genes that may cause DBA. We also identified several distinct disorders that appear to phenocopy DBA, including 9 individuals with biallelicCECR1mutations that result in deficiency of ADA2. However, no new genes were identified at exome-wide significance, suggesting that there are no unidentified genes containing mutations readily identified by WES that explain > 5% of DBA cases. Overall, this comprehensive report should not only inform clinical practice for DBA patients, but also the design and analysis of future rare variant studies for heterogeneous Mendelian disorders.

A Correction to this paper has been published: https://doi.org/10.1038/s41586-020-03176-6.

The genetic architecture of common traits, including the number, frequency, and effect sizes of inherited variants that contribute to individual risk, has been long debated. Genome-wide association studies have identified scores of common variants associated with type 2 diabetes, but in aggregate, these explain only a fraction of the heritability of this disease. Here, to test the hypothesis that lower-frequency variants explain much of the remainder, the GoT2D and T2D-GENES consortia performed whole-genome sequencing in 2,657 European individuals with and without diabetes, and exome sequencing in 12,940 individuals from five ancestry groups. To increase statistical power, we expanded the sample size via genotyping and imputation in a further 111,548 subjects. Variants associated with type 2 diabetes after sequencing were overwhelmingly common and most fell within regions previously identified by genome-wide association studies. Comprehensive enumeration of sequence variation is necessary to identify functional alleles that provide important clues to disease pathophysiology, but large-scale sequencing does not support the idea that lower-frequency variants have a major role in predisposition to type 2 diabetes.

Mucin-1 kidney disease, previously described as medullary cystic kidney disease type 1 (MCKD1, OMIM 174000), is an autosomal dominant tubulointerstitial kidney disease recently shown to be caused by a single-base insertion within the variable number tandem repeat region of the MUC1 gene. Because of variable age of disease onset and often subtle signs and symptoms, clinical diagnosis of mucin-1 kidney disease and differentiation from other forms of hereditary kidney disease have been difficult. The causal insertion resides in a variable number tandem repeat region with high GC content, which has made detection by standard next-generation sequencing impossible to date. The inherently difficult nature of this mutation required an alternative method for routine detection and clinical diagnosis of the disease. We therefore developed and validated a mass spectrometry–based probe extension assay with a series of internal controls to detect the insertion event using 24 previously characterized positive samples from patients with mucin-1 kidney disease and 24 control samples known to be wild type for the variant. Validation results indicate an accurate and reliable test for clinically establishing the molecular diagnosis of mucin-1 kidney disease with 100% sensitivity and specificity across 275 tests called., Carlos Slim Foundation. Slim Initiative for Genomic Medicine

Although studies have identified hundreds of loci associated with human traits and diseases, pinpointing causal alleles remains difficult, particularly for non-coding variants. To address this challenge, we adapted the massively parallel reporter assay (MPRA) to identify variants that directly modulate gene expression. We applied it to 32,373 variants from 3,642 cis-expression quantitative trait loci and control regions. Detection by MPRA was strongly correlated with measures of regulatory function. We demonstrate MPRA’s capabilities for pinpointing causal alleles, using it to identify 842 variants showing differential expression between alleles, including 53 well-annotated variants associated with diseases and traits. We investigated one in detail, a risk allele for ankylosing spondylitis, and provide direct evidence of a non-coding variant that alters expression of the prostaglandin EP4 receptor. These results create a resource of concrete leads and illustrate the promise of this approach for comprehensively interrogating how non-coding polymorphism shapes human biology., National Institutes of Health (U.S.) (grant DP2OD006514), National Institutes of Health (U.S.) (grant K99HG0081), National Institutes of Health (U.S.) (grant R01HG006785)

Three years ago, scientists reported that CRISPR technology can enable precise and efficient genome editing in living eukaryotic cells. Since then, the method has taken the scientific community by storm, with thousands of labs using it for applications from biomedicine to agriculture. Yet, the preceding 20-year journey--the discovery of a strange microbial repeat sequence; its recognition as an adaptive immune system; its biological characterization; and its repurposing for genome engineering--remains little known. This Perspective aims to fill in this backstory--the history of ideas and the stories of pioneers--and draw lessons about the remarkable ecosystem underlying scientific discovery.

Background: Approximately 7% of American adults have severe hypercholesterolemia (untreated low-density lipoprotein [LDL] cholesterol ≥190 mg/dl), which may be due to familial hypercholeste...

A typo in the 'Reviewer information' section of this Letter was corrected online.

New antibiotics are needed to combat rising resistance, with new Mycobacterium tuberculosis (Mtb) drugs of highest priority. Conventional whole-cell and biochemical antibiotic screens have failed. We developed a novel strategy termed PROSPECT (PRimary screening Of Strains to Prioritize Expanded Chemistry and Targets) in which we screen compounds against pools of strains depleted for essential bacterial targets. We engineered strains targeting 474 Mtb essential genes and screened pools of 100-150 strains against activity-enriched and unbiased compounds libraries, measuring > 8.5-million chemical-genetic interactions. Primary screens identified >10-fold more hits than screening wild-type Mtb alone, with chemical-genetic interactions providing immediate, direct target insight. We identified > 40 novel compounds targeting DNA gyrase, cell wall, tryptophan, folate biosynthesis, and RNA polymerase, as well as inhibitors of a novel target EfpA. Chemical optimization yielded EfpA inhibitors with potent wild-type activity, thus demonstrating PROSPECT’s ability to yield inhibitors against novel targets which would have eluded conventional drug discovery.

Significance Gene expression is controlled by sequence-specific transcription factors (TFs), which bind to regulatory sequences in DNA. The degree to which the arrangement of motif sites within regulatory elements determines their function remains unclear. Here, we show that the positional distribution of TF motif sites within nucleosome-depleted regions of DNA fall into six distinct classes. These patterns are highly consistent across cell types and bring together factors that have similar functional and binding properties. Furthermore, the position of motif sites appears to be related to their known functions. Our results suggest that TFs play distinct roles in forming a functional enhancer, facilitated by their position within a regulatory sequence., Gene expression is controlled by sequence-specific transcription factors (TFs), which bind to regulatory sequences in DNA. TF binding occurs in nucleosome-depleted regions of DNA (NDRs), which generally encompass regions with lengths similar to those protected by nucleosomes. However, less is known about where within these regions specific TFs tend to be found. Here, we characterize the positional bias of inferred binding sites for 103 TFs within ∼500,000 NDRs across 47 cell types. We find that distinct classes of TFs display different binding preferences: Some tend to have binding sites toward the edges, some toward the center, and some at other positions within the NDR. These patterns are highly consistent across cell types, suggesting that they may reflect TF-specific intrinsic structural or functional characteristics. In particular, TF classes with binding sites at NDR edges are enriched for those known to interact with histones and chromatin remodelers, whereas TFs with central enrichment interact with other TFs and cofactors such as p300. Our results suggest distinct regiospecific binding patterns and functions of TF classes within enhancers.

Bacterial community acquisition in the infant gut impacts immune education and disease susceptibility. We compared bacterial strains across and within families in a prospective birth cohort of 44 infants and their mothers, sampled longitudinally in the first months of each child's life. We identified mother-to-child bacterial transmission events and describe the incidence of family-specific antibiotic resistance genes. We observed two inheritance patterns across multiple species, where often the mother's dominant strain is transmitted to the child, but occasionally her secondary strains colonize the infant gut. In families where the secondary strain of B. uniformis was inherited, a starch utilization gene cluster that was absent in the mother's dominant strain was identified in the child, suggesting the selective advantage of a mother's secondary strain in the infant gut. Our findings reveal mother-to-child bacterial transmission events at high resolution and give insights into early colonization of the infant gut. Using longitudinal metagenomic sequencing from 44 mother/child pairs, Yassour et al. characterized mother-to-child strain transmission patterns. While mothers' dominant strains were often inherited, nondominant secondary strain transmissions were also observed. Microbial functional analysis reveals that inherited maternal secondary strains may have a selective advantage to colonize infant guts., National Institutes of Health (Grant 1DP3DK094338–01)

Understanding the molecular programs that guide differentiation during development is a major challenge. Here, we introduce Waddington-OT, an approach for studying developmental time courses to infer ancestor-descendant fates and model the regulatory programs that underlie them. We apply the method to reconstruct the landscape of reprogramming from 315,000 single-cell RNA sequencing (scRNA-seq) profiles, collected at half-day intervals across 18 days. The results reveal a wider range of developmental programs than previously characterized. Cells gradually adopt either a terminal stromal state or a mesenchymal-to-epithelial transition state. The latter gives rise to populations related to pluripotent, extra-embryonic, and neural cells, with each harboring multiple finer subpopulations. The analysis predicts transcription factors and paracrine signals that affect fates and experiments validate that the TF Obox6 and the cytokine GDF9 enhance reprogramming efficiency. Our approach sheds light on the process and outcome of reprogramming and provides a framework applicable to diverse temporal processes in biology., National Institutes of Health (U.S.) (grant HD045022), National Institutes of Health (U.S.) (grant R01 MH104610-15), National Institutes of Health (U.S.) (grant R01NS088538)

Blood cell formation is classically thought to occur through a hierarchical differentiation process, although recent studies have shown that lineage commitment may occur earlier in hematopoietic stem and progenitor cells (HSPCs). The relevance to human blood diseases and the underlying regulation of these refined models remain poorly understood. By studying a genetic blood disorder, Diamond-Blackfan anemia (DBA), where the majority of mutations affect ribosomal proteins and the erythroid lineage is selectively perturbed, we are able to gain mechanistic insight into how lineage commitment is programmed normally and disrupted in disease. We show that in DBA, the pool of available ribosomes is limited, while ribosome composition remains constant. Surprisingly, this global reduction in ribosome levels more profoundly alters translation of a select subset of transcripts. We show how the reduced translation of select transcripts in HSPCs can impair erythroid lineage commitment, illuminating a regulatory role for ribosome levels in cellular differentiation. A global reduction in ribosome levels in Diamond-Blackfan anemia profoundly alters translation of a select subset of transcripts, thereby impeding erythroid lineage commitment. Keyword: hematopoiesis; lineage commitment; GATA1; genetics; ribosome; translation; Diamond-Blackman anemia; erythropoiesis, National Institutes of Health (U.S.) (Grants R01 DK103794), National Institutes of Health (U.S.) (Grant R33 HL120791)