The loud minority: Transcriptionally active HIV-1-infected cells survive, proliferate, and persist.

The shock-and-kill strategy reactivates HIV-1 latent reservoir for immune clearance. Einkauf et al. found that some HIV-1-infected cells that persist and proliferate have transcriptionally active HIV-1 in permissive chromatin. Silent proviruses in repressive chromatin resist reactivation. Understanding HIV-1-chromatin interactions and how transcriptionally active HIV-1-infected cells survive is a pressing need.

Large-Scale Exome Sequencing Study Implicates Both Developmental and Functional Changes in the Neurobiology of Autism.

We present the largest exome sequencing study of autism spectrum disorder (ASD) to date (n = 35,584 total samples, 11,986 with ASD). Using an enhanced analytical framework to integrate de novo and case-control rare variation, we identify 102 risk genes at a false discovery rate of 0.1 or less. Of these genes, 49 show higher frequencies of disruptive de novo variants in individuals ascertained to have severe neurodevelopmental delay, whereas 53 show higher frequencies in individuals ascertained to have ASD; comparing ASD cases with mutations in these groups reveals phenotypic differences. Expressed early in brain development, most risk genes have roles in regulation of gene expression or neuronal communication (i.e., mutations effect neurodevelopmental and neurophysiological changes), and 13 fall within loci recurrently hit by copy number variants. In cells from the human cortex, expression of risk genes is enriched in excitatory and inhibitory neuronal lineages, consistent with multiple paths to an excitatory-inhibitory imbalance underlying ASD.

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.

Single-cell epigenetic, transcriptional, and protein profiling of latent and active HIV-1 reservoir revealed that IKZF3 promotes HIV-1 persistence.

Understanding how HIV-1-infected cells proliferate and persist is key to HIV-1 eradication, but the heterogeneity and rarity of HIV-1-infected cells hamper mechanistic interrogations. Here, we used single-cell DOGMA-seq to simultaneously capture transcription factor accessibility, transcriptome, surface proteins, HIV-1 DNA, and HIV-1 RNA in memory CD4+ T cells from six people living with HIV-1 during viremia and after suppressive antiretroviral therapy. We identified increased transcription factor accessibility in latent HIV-1-infected cells (RORC) and transcriptionally active HIV-1-infected cells (interferon regulatory transcription factor [IRF] and activator protein 1 [AP-1]). A proliferation program (IKZF3, IL21, BIRC5, and MKI67 co-expression) promoted the survival of transcriptionally active HIV-1-infected cells. Both latent and transcriptionally active HIV-1-infected cells had increased IKZF3 (Aiolos) expression. Distinct epigenetic programs drove the heterogeneous cellular states of HIV-1-infected cells: IRF:activation, Eomes:cytotoxic effector differentiation, AP-1:migration, and cell death. Our study revealed the single-cell epigenetic, transcriptional, and protein states of latent and transcriptionally active HIV-1-infected cells and cellular programs promoting HIV-1 persistence.

The tumor-derived cytokine Chi3l1 induces neutrophil extracellular traps that promote T cell exclusion in triple-negative breast cancer.

In triple-negative breast cancer (TNBC), stromal restriction of CD8+ T cells associates with poor clinical outcomes and lack of responsiveness to immune-checkpoint blockade (ICB). To identify mediators of T cell stromal restriction, we profiled murine breast tumors lacking the transcription factor Stat3, which is commonly hyperactive in breast cancers and promotes an immunosuppressive tumor microenvironment. Expression of the cytokine Chi3l1 was decreased in Stat3-/- tumors. CHI3L1 expression was elevated in human TNBCs and other solid tumors exhibiting T cell stromal restriction. Chi3l1 ablation in the polyoma virus middle T (PyMT) breast cancer model generated an anti-tumor immune response and delayed mammary tumor onset. These effects were associated with increased T cell tumor infiltration and improved response to ICB. Mechanistically, Chi3l1 promoted neutrophil recruitment and neutrophil extracellular trap formation, which blocked T cell infiltration. Our findings provide insight into the mechanism underlying stromal restriction of CD8+ T cells and suggest that targeting Chi3l1 may promote anti-tumor immunity in various tumor types.

Loss of CD4+ T cell-intrinsic arginase 1 accelerates Th1 response kinetics and reduces lung pathology during influenza infection.

Arginase 1 (Arg1), the enzyme catalyzing the conversion of arginine to ornithine, is a hallmark of IL-10-producing immunoregulatory M2 macrophages. However, its expression in T cells is disputed. Here, we demonstrate that induction of Arg1 expression is a key feature of lung CD4+ T cells during mouse in vivo influenza infection. Conditional ablation of Arg1 in CD4+ T cells accelerated both virus-specific T helper 1 (Th1) effector responses and its resolution, resulting in efficient viral clearance and reduced lung pathology. Using unbiased transcriptomics and metabolomics, we found that Arg1-deficiency was distinct from Arg2-deficiency and caused altered glutamine metabolism. Rebalancing this perturbed glutamine flux normalized the cellular Th1 response. CD4+ T cells from rare ARG1-deficient patients or CRISPR-Cas9-mediated ARG1-deletion in healthy donor cells phenocopied the murine cellular phenotype. Collectively, CD4+ T cell-intrinsic Arg1 functions as an unexpected rheostat regulating the kinetics of the mammalian Th1 lifecycle with implications for Th1-associated tissue pathologies.

Single-cell multiomics reveals persistence of HIV-1 in expanded cytotoxic T cell clones.

Understanding the drivers and markers of clonally expanding HIV-1-infected CD4+ T cells is essential for HIV-1 eradication. We used single-cell ECCITE-seq, which captures surface protein expression, cellular transcriptome, HIV-1 RNA, and TCR sequences within the same single cell to track clonal expansion dynamics in longitudinally archived samples from six HIV-1-infected individuals (during viremia and after suppressive antiretroviral therapy) and two uninfected individuals, in unstimulated conditions and after CMV and HIV-1 antigen stimulation. Despite antiretroviral therapy, persistent antigen and TNF responses shaped T cell clonal expansion. HIV-1 resided in Th1-polarized, antigen-responding T cells expressing BCL2 and SERPINB9 that may resist cell death. HIV-1 RNA+ T cell clones were larger in clone size, established during viremia, persistent after viral suppression, and enriched in GZMB+ cytotoxic effector memory Th1 cells. Targeting HIV-1-infected cytotoxic CD4+ T cells and drivers of clonal expansion provides another direction for HIV-1 eradication.

Regioselective aliphatic C-H functionalization using frustrated radical pairs.

Frustrated Lewis pairs (FLPs) are well documented for the activation of small molecules such as dihydrogen and carbon dioxide1-4. Although canonical FLP chemistry is heterolytic in nature, recent work has shown that certain FLPs can undergo single-electron transfer to afford radical pairs5. Owing to steric encumbrance and/or weak bonding association, these radicals do not annihilate one another, and they have thus been named frustrated radical pairs (FRPs). Notable preliminary results suggest that FRPs may be useful reagents in chemical synthesis6-8, although their applications remain limited. Here we demonstrate that the functionalization of C(sp3)-H bonds can be accomplished using a class of FRPs generated from disilazide donors and an N-oxoammonium acceptor. Together, these species undergo single-electron transfer to generate a transient and persistent radical pair capable of cleaving unactivated C-H bonds to furnish aminoxylated products. By tuning the structure of the donor, it is possible to control regioselectivity and tailor reactivity towards tertiary, secondary or primary C-H bonds. Mechanistic studies lend strong support for the formation and involvement of radical pairs in the target reaction.

Toward effective probiotics for autism and other neurodevelopmental disorders.

Hsaio and colleagues link gut microbes to autism spectrum disorders (ASD) in a mouse model. They show that ASD symptoms are triggered by compositional and structural shifts of microbes and associated metabolites, but symptoms are relieved by a Bacteroides fragilis probiotic. Thus probiotics may provide therapeutic strategies for neurodevelopmental disorders.

Whole-genome sequencing reveals host factors underlying critical COVID-19.

Critical COVID-19 is caused by immune-mediated inflammatory lung injury. Host genetic variation influences the development of illness requiring critical care1 or hospitalization2-4 after infection with SARS-CoV-2. The GenOMICC (Genetics of Mortality in Critical Care) study enables the comparison of genomes from individuals who are critically ill with those of population controls to find underlying disease mechanisms. Here we use whole-genome sequencing in 7,491 critically ill individuals compared with 48,400 controls to discover and replicate 23 independent variants that significantly predispose to critical COVID-19. We identify 16 new independent associations, including variants within genes that are involved in interferon signalling (IL10RB and PLSCR1), leucocyte differentiation (BCL11A) and blood-type antigen secretor status (FUT2). Using transcriptome-wide association and colocalization to infer the effect of gene expression on disease severity, we find evidence that implicates multiple genes-including reduced expression of a membrane flippase (ATP11A), and increased expression of a mucin (MUC1)-in critical disease. Mendelian randomization provides evidence in support of causal roles for myeloid cell adhesion molecules (SELE, ICAM5 and CD209) and the coagulation factor F8, all of which are potentially druggable targets. Our results are broadly consistent with a multi-component model of COVID-19 pathophysiology, in which at least two distinct mechanisms can predispose to life-threatening disease: failure to control viral replication; or an enhanced tendency towards pulmonary inflammation and intravascular coagulation. We show that comparison between cases of critical illness and population controls is highly efficient for the detection of therapeutically relevant mechanisms of disease.

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.

Ribosome decision graphs for the representation of eukaryotic RNA translation complexity.

The application of ribosome profiling has revealed an unexpected abundance of translation in addition to that responsible for the synthesis of previously annotated protein-coding regions. Multiple short sequences have been found to be translated within single RNA molecules, within both annotated protein-coding and noncoding regions. The biological significance of this translation is a matter of intensive investigation. However, current schematic or annotation-based representations of mRNA translation generally do not account for the apparent multitude of translated regions within the same molecules. They also do not take into account the stochasticity of the process that allows alternative translations of the same RNA molecules by different ribosomes. There is a need for formal representations of mRNA complexity that would enable the analysis of quantitative information on translation and more accurate models for predicting the phenotypic effects of genetic variants affecting translation. To address this, we developed a conceptually novel abstraction that we term ribosome decision graphs (RDGs). RDGs represent translation as multiple ribosome paths through untranslated and translated mRNA segments. We termed the latter "translons." Nondeterministic events, such as initiation, reinitiation, selenocysteine insertion, or ribosomal frameshifting, are then represented as branching points. This representation allows for an adequate representation of eukaryotic translation complexity and focuses on locations critical for translation regulation. We show how RDGs can be used for depicting translated regions and for analyzing genetic variation and quantitative genome-wide data on translation for characterization of regulatory modulators of translation.

Time is of the essence: The importance of considering biological rhythms in an increasingly polluted world.

Biological rhythms have a crucial role in shaping the biology and ecology of organisms. Light pollution is known to disrupt these rhythms, and evidence is emerging that chemical pollutants can cause similar disruption. Conversely, biological rhythms can influence the effects and toxicity of chemicals. Thus, by drawing insights from the extensive study of biological rhythms in biomedical and light pollution research, we can greatly improve our understanding of chemical pollution. This Essay advocates for the integration of biological rhythmicity into chemical pollution research to gain a more comprehensive understanding of how chemical pollutants affect wildlife and ecosystems. Despite historical barriers, recent experimental and technological advancements now facilitate the integration of biological rhythms into ecotoxicology, offering unprecedented, high-resolution data across spatiotemporal scales. Recognizing the importance of biological rhythms will be essential for understanding, predicting, and mitigating the complex ecological repercussions of chemical pollution.

Elevated NSD3 histone methylation activity drives squamous cell lung cancer.

Amplification of chromosomal region 8p11-12 is a common genetic alteration that has been implicated in the aetiology of lung squamous cell carcinoma (LUSC)1-3. The FGFR1 gene is the main candidate driver of tumorigenesis within this region4. However, clinical trials evaluating FGFR1 inhibition as a targeted therapy have been unsuccessful5. Here we identify the histone H3 lysine 36 (H3K36) methyltransferase NSD3, the gene for which is located in the 8p11-12 amplicon, as a key regulator of LUSC tumorigenesis. In contrast to other 8p11-12 candidate LUSC drivers, increased expression of NSD3 correlated strongly with its gene amplification. Ablation of NSD3, but not of FGFR1, attenuated tumour growth and extended survival in a mouse model of LUSC. We identify an LUSC-associated variant NSD3(T1232A) that shows increased catalytic activity for dimethylation of H3K36 (H3K36me2) in vitro and in vivo. Structural dynamic analyses revealed that the T1232A substitution elicited localized mobility changes throughout the catalytic domain of NSD3 to relieve auto-inhibition and to increase accessibility of the H3 substrate. Expression of NSD3(T1232A) in vivo accelerated tumorigenesis and decreased overall survival in mouse models of LUSC. Pathological generation of H3K36me2 by NSD3(T1232A) reprograms the chromatin landscape to promote oncogenic gene expression signatures. Furthermore, NSD3, in a manner dependent on its catalytic activity, promoted transformation in human tracheobronchial cells and growth of xenografted human LUSC cell lines with amplification of 8p11-12. Depletion of NSD3 in patient-derived xenografts from primary LUSCs containing NSD3 amplification or the NSD3(T1232A)-encoding variant attenuated neoplastic growth in mice. Finally, NSD3-regulated LUSC-derived xenografts were hypersensitive to bromodomain inhibition. Thus, our work identifies NSD3 as a principal 8p11-12 amplicon-associated oncogenic driver in LUSC, and suggests that NSD3-dependency renders LUSC therapeutically vulnerable to bromodomain inhibition.

Exome sequencing and analysis of 454,787 UK Biobank participants.

A major goal in human genetics is to use natural variation to understand the phenotypic consequences of altering each protein-coding gene in the genome. Here we used exome sequencing1 to explore protein-altering variants and their consequences in 454,787 participants in the UK Biobank study2. We identified 12 million coding variants, including around 1 million loss-of-function and around 1.8 million deleterious missense variants. When these were tested for association with 3,994 health-related traits, we found 564 genes with trait associations at P ≤ 2.18 × 10-11. Rare variant associations were enriched in loci from genome-wide association studies (GWAS), but most (91%) were independent of common variant signals. We discovered several risk-increasing associations with traits related to liver disease, eye disease and cancer, among others, as well as risk-lowering associations for hypertension (SLC9A3R2), diabetes (MAP3K15, FAM234A) and asthma (SLC27A3). Six genes were associated with brain imaging phenotypes, including two involved in neural development (GBE1, PLD1). Of the signals available and powered for replication in an independent cohort, 81% were confirmed; furthermore, association signals were generally consistent across individuals of European, Asian and African ancestry. We illustrate the ability of exome sequencing to identify gene-trait associations, elucidate gene function and pinpoint effector genes that underlie GWAS signals at scale.

Novel breast cancer susceptibility loci under linkage peaks identified in African ancestry consortia.

BACKGROUND: Expansion of genome-wide association studies across population groups is needed to improve our understanding of shared and unique genetic contributions to breast cancer. We performed association and replication studies guided by a priori linkage findings from African ancestry (AA) relative pairs. METHODS: We performed fixed-effect inverse-variance weighted meta-analysis under three significant AA breast cancer linkage peaks (3q26-27, 12q22-23, and 16q21-22) in 9241 AA cases and 10 193 AA controls. We examined associations with overall breast cancer as well as estrogen receptor (ER)-positive and negative subtypes (193,132 SNPs). We replicated associations in the African-ancestry Breast Cancer Genetic Consortium (AABCG). RESULTS: In AA women, we identified two associations on chr12q for overall breast cancer (rs1420647, OR = 1.15, p = 2.50×10-6; rs12322371, OR = 1.14, p = 3.15×10-6), and one for ER-negative breast cancer (rs77006600, OR = 1.67, p = 3.51×10-6). On chr3, we identified two associations with ER-negative disease (rs184090918, OR = 3.70, p = 1.23×10-5; rs76959804, OR = 3.57, p = 1.77×10-5) and on chr16q we identified an association with ER-negative disease (rs34147411, OR = 1.62, p = 8.82×10-6). In the replication study, the chr3 associations were significant and effect sizes were larger (rs184090918, OR: 6.66, 95% CI: 1.43, 31.01; rs76959804, OR: 5.24, 95% CI: 1.70, 16.16). CONCLUSION: The two chr3 SNPs are upstream to open chromatin ENSR00000710716, a regulatory feature that is actively regulated in mammary tissues, providing evidence that variants in this chr3 region may have a regulatory role in our target organ. Our study provides support for breast cancer variant discovery using prioritization based on linkage evidence.

Integration site-dependent HIV-1 promoter activity shapes host chromatin conformation.

HIV-1 integration introduces ectopic transcription factor binding sites into host chromatin. We postulate that the integrated provirus serves as an ectopic enhancer that recruits additional transcription factors to the integration locus, increases chromatin accessibility, changes 3D chromatin interactions, and enhances both retroviral and host gene expression. We used four well-characterized HIV-1-infected cell line clones having unique integration sites and low to high levels of HIV-1 expression. Using single-cell DOGMA-seq, which captured the heterogeneity of HIV-1 expression and host chromatin accessibility, we found that HIV-1 transcription correlated with HIV-1 accessibility and host chromatin accessibility. HIV-1 integration increased local host chromatin accessibility within an ∼5- to 30-kb distance. CRISPRa- and CRISPRi-mediated HIV-1 promoter activation and inhibition confirmed integration site-dependent HIV-1-driven changes of host chromatin accessibility. HIV-1 did not drive chromatin confirmation changes at the genomic level (by Hi-C) or the enhancer connectome (by H3K27ac HiChIP). Using 4C-seq to interrogate HIV-1-chromatin interactions, we found that HIV-1 interacted with host chromatin ∼100-300 kb from the integration site. By identifying chromatin regions having both increased transcription factor activity (by ATAC-seq) and HIV-1-chromatin interaction (by 4C-seq), we identified enrichment of ETS, RUNT, and ZNF-family transcription factor binding that may mediate HIV-1-host chromatin interactions. Our study has found that HIV-1 promoter activity increases host chromatin accessibility, and HIV-1 interacted with host chromatin within the existing chromatin boundaries in an integration site-dependent manner.

Revisiting Cardiac Biology in the Era of Single Cell and Spatial Omics.

Throughout our lifetime, each beat of the heart requires the coordinated action of multiple cardiac cell types. Understanding cardiac cell biology, its intricate microenvironments, and the mechanisms that govern their function in health and disease are crucial to designing novel therapeutical and behavioral interventions. Recent advances in single-cell and spatial omics technologies have significantly propelled this understanding, offering novel insights into the cellular diversity and function and the complex interactions of cardiac tissue. This review provides a comprehensive overview of the cellular landscape of the heart, bridging the gap between suspension-based and emerging in situ approaches, focusing on the experimental and computational challenges, comparative analyses of mouse and human cardiac systems, and the rising contextualization of cardiac cells within their niches. As we explore the heart at this unprecedented resolution, integrating insights from both mouse and human studies will pave the way for novel diagnostic tools and therapeutic interventions, ultimately improving outcomes for patients with cardiovascular diseases.

Transcript expression-aware annotation improves rare variant interpretation.

The acceleration of DNA sequencing in samples from patients and population studies has resulted in extensive catalogues of human genetic variation, but the interpretation of rare genetic variants remains problematic. A notable example of this challenge is the existence of disruptive variants in dosage-sensitive disease genes, even in apparently healthy individuals. Here, by manual curation of putative loss-of-function (pLoF) variants in haploinsufficient disease genes in the Genome Aggregation Database (gnomAD)1, we show that one explanation for this paradox involves alternative splicing of mRNA, which allows exons of a gene to be expressed at varying levels across different cell types. Currently, no existing annotation tool systematically incorporates information about exon expression into the interpretation of variants. We develop a transcript-level annotation metric known as the 'proportion expressed across transcripts', which quantifies isoform expression for variants. We calculate this metric using 11,706 tissue samples from the Genotype Tissue Expression (GTEx) project2 and show that it can differentiate between weakly and highly evolutionarily conserved exons, a proxy for functional importance. We demonstrate that expression-based annotation selectively filters 22.8% of falsely annotated pLoF variants found in haploinsufficient disease genes in gnomAD, while removing less than 4% of high-confidence pathogenic variants in the same genes. Finally, we apply our expression filter to the analysis of de novo variants in patients with autism spectrum disorder and intellectual disability or developmental disorders to show that pLoF variants in weakly expressed regions have similar effect sizes to those of synonymous variants, whereas pLoF variants in highly expressed exons are most strongly enriched among cases. Our annotation is fast, flexible and generalizable, making it possible for any variant file to be annotated with any isoform expression dataset, and will be valuable for the genetic diagnosis of rare diseases, the analysis of rare variant burden in complex disorders, and the curation and prioritization of variants in recall-by-genotype studies.

The mutational constraint spectrum quantified from variation in 141,456 humans.

Genetic variants that inactivate protein-coding genes are a powerful source of information about the phenotypic consequences of gene disruption: genes that are crucial for the function of an organism will be depleted of such variants in natural populations, whereas non-essential genes will tolerate their accumulation. However, predicted loss-of-function variants are enriched for annotation errors, and tend to be found at extremely low frequencies, so their analysis requires careful variant annotation and very large sample sizes1. Here we describe the aggregation of 125,748 exomes and 15,708 genomes from human sequencing studies into the Genome Aggregation Database (gnomAD). We identify 443,769 high-confidence predicted loss-of-function variants in this cohort after filtering for artefacts caused by sequencing and annotation errors. Using an improved model of human mutation rates, we classify human protein-coding genes along a spectrum that represents tolerance to inactivation, validate this classification using data from model organisms and engineered human cells, and show that it can be used to improve the power of gene discovery for both common and rare diseases.