PLEKHS1 drives PI3Ks and remodels pathway homeostasis in PTEN-null prostate.

The PIP3/PI3K network is a central regulator of metabolism and is frequently activated in cancer, commonly by loss of the PIP3/PI(3,4)P2 phosphatase, PTEN. Despite huge research investment, the drivers of the PI3K network in normal tissues and how they adapt to overactivation are unclear. We find that in healthy mouse prostate PI3K activity is driven by RTK/IRS signaling and constrained by pathway feedback. In the absence of PTEN, the network is dramatically remodeled. A poorly understood YXXM- and PIP3/PI(3,4)P2-binding PH domain-containing adaptor, PLEKHS1, became the dominant activator and was required to sustain PIP3, AKT phosphorylation, and growth in PTEN-null prostate. This was because PLEKHS1 evaded pathway-feedback and experienced enhanced PI3K- and Src-family kinase-dependent phosphorylation of Y258XXM, eliciting PI3K activation. hPLEKHS1 mRNA and activating Y419 phosphorylation of hSrc correlated with PI3K pathway activity in human prostate cancers. We propose that in PTEN-null cells receptor-independent, Src-dependent tyrosine phosphorylation of PLEKHS1 creates positive feedback that escapes homeostasis, drives PIP3 signaling, and supports tumor progression.

FOXG1-Dependent Dysregulation of GABA/Glutamate Neuron Differentiation in Autism Spectrum Disorders.

Autism spectrum disorder (ASD) is a disorder of brain development. Most cases lack a clear etiology or genetic basis, and the difficulty of re-enacting human brain development has precluded understanding of ASD pathophysiology. Here we use three-dimensional neural cultures (organoids) derived from induced pluripotent stem cells (iPSCs) to investigate neurodevelopmental alterations in individuals with severe idiopathic ASD. While no known underlying genomic mutation could be identified, transcriptome and gene network analyses revealed upregulation of genes involved in cell proliferation, neuronal differentiation, and synaptic assembly. ASD-derived organoids exhibit an accelerated cell cycle and overproduction of GABAergic inhibitory neurons. Using RNA interference, we show that overexpression of the transcription factor FOXG1 is responsible for the overproduction of GABAergic neurons. Altered expression of gene network modules and FOXG1 are positively correlated with symptom severity. Our data suggest that a shift toward GABAergic neuron fate caused by FOXG1 is a developmental precursor of ASD.

A genomic mutational constraint map using variation in 76,156 human genomes.

The depletion of disruptive variation caused by purifying natural selection (constraint) has been widely used to investigate protein-coding genes underlying human disorders1-4, but attempts to assess constraint for non-protein-coding regions have proved more difficult. Here we aggregate, process and release a dataset of 76,156 human genomes from the Genome Aggregation Database (gnomAD)-the largest public open-access human genome allele frequency reference dataset-and use it to build a genomic constraint map for the whole genome (genomic non-coding constraint of haploinsufficient variation (Gnocchi)). We present a refined mutational model that incorporates local sequence context and regional genomic features to detect depletions of variation. As expected, the average constraint for protein-coding sequences is stronger than that for non-coding regions. Within the non-coding genome, constrained regions are enriched for known regulatory elements and variants that are implicated in complex human diseases and traits, facilitating the triangulation of biological annotation, disease association and natural selection to non-coding DNA analysis. More constrained regulatory elements tend to regulate more constrained protein-coding genes, which in turn suggests that non-coding constraint can aid the identification of constrained genes that are as yet unrecognized by current gene constraint metrics. We demonstrate that this genome-wide constraint map improves the identification and interpretation of functional human genetic variation.

Non-canonical autophagy drives alternative ATG8 conjugation to phosphatidylserine.

Autophagy is a fundamental catabolic process that uses a unique post-translational modification, the conjugation of ATG8 protein to phosphatidylethanolamine (PE). ATG8 lipidation also occurs during non-canonical autophagy, a parallel pathway involving conjugation of ATG8 to single membranes (CASM) at endolysosomal compartments, with key functions in immunity, vision, and neurobiology. It is widely assumed that CASM involves the same conjugation of ATG8 to PE, but this has not been formally tested. Here, we discover that all ATG8s can also undergo alternative lipidation to phosphatidylserine (PS) during CASM, induced pharmacologically, by LC3-associated phagocytosis or influenza A virus infection, in mammalian cells. Importantly, ATG8-PS and ATG8-PE adducts are differentially delipidated by the ATG4 family and bear different cellular dynamics, indicating significant molecular distinctions. These results provide important insights into autophagy signaling, revealing an alternative form of the hallmark ATG8 lipidation event. Furthermore, ATG8-PS provides a specific "molecular signature" for the non-canonical autophagy pathway.

Cooperativity and rapid evolution of cobound transcription factors in closely related mammals.

To mechanistically characterize the microevolutionary processes active in altering transcription factor (TF) binding among closely related mammals, we compared the genome-wide binding of three tissue-specific TFs that control liver gene expression in six rodents. Despite an overall fast turnover of TF binding locations between species, we identified thousands of TF regions of highly constrained TF binding intensity. Although individual mutations in bound sequence motifs can influence TF binding, most binding differences occur in the absence of nearby sequence variations. Instead, combinatorial binding was found to be significant for genetic and evolutionary stability; cobound TFs tend to disappear in concert and were sensitive to genetic knockout of partner TFs. The large, qualitative differences in genomic regions bound between closely related mammals, when contrasted with the smaller, quantitative TF binding differences among Drosophila species, illustrate how genome structure and population genetics together shape regulatory evolution.

Clonally expanded B cells in multiple sclerosis bind EBV EBNA1 and GlialCAM.

Multiple sclerosis (MS) is a heterogenous autoimmune disease in which autoreactive lymphocytes attack the myelin sheath of the central nervous system. B lymphocytes in the cerebrospinal fluid (CSF) of patients with MS contribute to inflammation and secrete oligoclonal immunoglobulins1,2. Epstein-Barr virus (EBV) infection has been epidemiologically linked to MS, but its pathological role remains unclear3. Here we demonstrate high-affinity molecular mimicry between the EBV transcription factor EBV nuclear antigen 1 (EBNA1) and the central nervous system protein glial cell adhesion molecule (GlialCAM) and provide structural and in vivo functional evidence for its relevance. A cross-reactive CSF-derived antibody was initially identified by single-cell sequencing of the paired-chain B cell repertoire of MS blood and CSF, followed by protein microarray-based testing of recombinantly expressed CSF-derived antibodies against MS-associated viruses. Sequence analysis, affinity measurements and the crystal structure of the EBNA1-peptide epitope in complex with the autoreactive Fab fragment enabled tracking of the development of the naive EBNA1-restricted antibody to a mature EBNA1-GlialCAM cross-reactive antibody. Molecular mimicry is facilitated by a post-translational modification of GlialCAM. EBNA1 immunization exacerbates disease in a mouse model of MS, and anti-EBNA1 and anti-GlialCAM antibodies are prevalent in patients with MS. Our results provide a mechanistic link for the association between MS and EBV and could guide the development of new MS therapies.

Candidate Cancer Driver Mutations in Distal Regulatory Elements and Long-Range Chromatin Interaction Networks.

A comprehensive catalog of cancer driver mutations is essential for understanding tumorigenesis and developing therapies. Exome-sequencing studies have mapped many protein-coding drivers, yet few non-coding drivers are known because genome-wide discovery is challenging. We developed a driver discovery method, ActiveDriverWGS, and analyzed 120,788 cis-regulatory modules (CRMs) across 1,844 whole tumor genomes from the ICGC-TCGA PCAWG project. We found 30 CRMs with enriched SNVs and indels (FDR < 0.05). These frequently mutated regulatory elements (FMREs) were ubiquitously active in human tissues, showed long-range chromatin interactions and mRNA abundance associations with target genes, and were enriched in motif-rewiring mutations and structural variants. Genomic deletion of one FMRE in human cells caused proliferative deficiencies and transcriptional deregulation of cancer genes CCNB1IP1, CDH1, and CDKN2B, validating observations in FMRE-mutated tumors. Pathway analysis revealed further sub-significant FMREs at cancer genes and processes, indicating an unexplored landscape of infrequent driver mutations in the non-coding genome.

A harmonized public resource of deeply sequenced diverse human genomes.

Underrepresented populations are often excluded from genomic studies owing in part to a lack of resources supporting their analyses. The 1000 Genomes Project (1kGP) and Human Genome Diversity Project (HGDP), which have recently been sequenced to high coverage, are valuable genomic resources because of the global diversity they capture and their open data sharing policies. Here, we harmonized a high-quality set of 4094 whole genomes from 80 populations in the HGDP and 1kGP with data from the Genome Aggregation Database (gnomAD) and identified over 153 million high-quality SNVs, indels, and SVs. We performed a detailed ancestry analysis of this cohort, characterizing population structure and patterns of admixture across populations, analyzing site frequency spectra, and measuring variant counts at global and subcontinental levels. We also show substantial added value from this data set compared with the prior versions of the component resources, typically combined via liftOver and variant intersection; for example, we catalog millions of new genetic variants, mostly rare, compared with previous releases. In addition to unrestricted individual-level public release, we provide detailed tutorials for conducting many of the most common quality-control steps and analyses with these data in a scalable cloud-computing environment and publicly release this new phased joint callset for use as a haplotype resource in phasing and imputation pipelines. This jointly called reference panel will serve as a key resource to support research of diverse ancestry populations.

Waves of retrotransposon expansion remodel genome organization and CTCF binding in multiple mammalian lineages.

CTCF-binding locations represent regulatory sequences that are highly constrained over the course of evolution. To gain insight into how these DNA elements are conserved and spread through the genome, we defined the full spectrum of CTCF-binding sites, including a 33/34-mer motif, and identified over five thousand highly conserved, robust, and tissue-independent CTCF-binding locations by comparing ChIP-seq data from six mammals. Our data indicate that activation of retroelements has produced species-specific expansions of CTCF binding in rodents, dogs, and opossum, which often functionally serve as chromatin and transcriptional insulators. We discovered fossilized repeat elements flanking deeply conserved CTCF-binding regions, indicating that similar retrotransposon expansions occurred hundreds of millions of years ago. Repeat-driven dispersal of CTCF binding is a fundamental, ancient, and still highly active mechanism of genome evolution in mammalian lineages.

Global absence and targeting of protective immune states in severe COVID-19.

Although infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has pleiotropic and systemic effects in some individuals1-3, many others experience milder symptoms. Here, to gain a more comprehensive understanding of the distinction between severe and mild phenotypes in the pathology of coronavirus disease 2019 (COVID-19) and its origins, we performed a whole-blood-preserving single-cell analysis protocol to integrate contributions from all major immune cell types of the blood-including neutrophils, monocytes, platelets, lymphocytes and the contents of the serum. Patients with mild COVID-19 exhibit a coordinated pattern of expression of interferon-stimulated genes (ISGs)3 across every cell population, whereas these ISG-expressing cells are systemically absent in patients with severe disease. Paradoxically, individuals with severe COVID-19 produce very high titres of anti-SARS-CoV-2 antibodies and have a lower viral load compared to individuals with mild disease. Examination of the serum from patients with severe COVID-19 shows that these patients uniquely produce antibodies that functionally block the production of the ISG-expressing cells associated with mild disease, by activating conserved signalling circuits that dampen cellular responses to interferons. Overzealous antibody responses pit the immune system against itself in many patients with COVID-19, and perhaps also in individuals with other viral infections. Our findings reveal potential targets for immunotherapies in patients with severe COVID-19 to re-engage viral defence.

Break Check: Transcription-Driven Topoisomerase II Collisions near Chromatin Loop Anchors Are Hotspots for DNA Damage and Translocations.

In this issue of Molecular Cell, Gothe, Maman et al. and Canela et al. demonstrate that type II topoisomerase (TOP2B/TOP2A)-mediated DNA breaks, and the oncogenic translocations they generate, are dependent on transcription and governed by chromatin architecture.

Interactome Rewiring Following Pharmacological Targeting of BET Bromodomains.

Targeting bromodomains (BRDs) of the bromo-and-extra-terminal (BET) family offers opportunities for therapeutic intervention in cancer and other diseases. Here, we profile the interactomes of BRD2, BRD3, BRD4, and BRDT following treatment with the pan-BET BRD inhibitor JQ1, revealing broad rewiring of the interaction landscape, with three distinct classes of behavior for the 603 unique interactors identified. A group of proteins associate in a JQ1-sensitive manner with BET BRDs through canonical and new binding modes, while two classes of extra-terminal (ET)-domain binding motifs mediate acetylation-independent interactions. Last, we identify an unexpected increase in several interactions following JQ1 treatment that define negative functions for BRD3 in the regulation of rRNA synthesis and potentially RNAPII-dependent gene expression that result in decreased cell proliferation. Together, our data highlight the contributions of BET protein modules to their interactomes allowing for a better understanding of pharmacological rewiring in response to JQ1.

Trio RNA sequencing in a cohort of medically complex children.

Genome sequencing (GS) is a powerful test for the diagnosis of rare genetic disorders. Although GS can enumerate most non-coding variation, determining which non-coding variants are disease-causing is challenging. RNA sequencing (RNA-seq) has emerged as an important tool to help address this issue, but its diagnostic utility remains understudied, and the added value of a trio design is unknown. We performed GS plus RNA-seq from blood using an automated clinical-grade high-throughput platform on 97 individuals from 39 families where the proband was a child with unexplained medical complexity. RNA-seq was an effective adjunct test when paired with GS. It enabled clarification of putative splice variants in three families, but it did not reveal variants not already identified by GS analysis. Trio RNA-seq decreased the number of candidates requiring manual review when filtering for de novo dominant disease-causing variants, allowing for the exclusion of 16% of gene-expression outliers and 27% of allele-specific-expression outliers. However, clear diagnostic benefit from the trio design was not observed. Blood-based RNA-seq can facilitate genome analysis in children with suspected undiagnosed genetic disease. In contrast to DNA sequencing, the clinical advantages of a trio RNA-seq design may be more limited.

The evolution of early hominin food production and sharing.

How did humans evolve from individualistic to collective foraging with sex differences in production and widespread sharing of plant and animal foods? While current evolutionary scenarios focus on meat, cooking, or grandparental subsidies, considerations of the economics of foraging for extracted plant foods (e.g., roots, tubers), inferred to be important for early hominins (∼6 to 2.5 mya), suggest that early hominins shared such foods with offspring and others. Here, we present a conceptual and mathematical model of early hominin food production and sharing, prior to the emergence of frequent hunting, cooking, and increased lifespan. We hypothesize that extracted plant foods were vulnerable to theft, and that male mate guarding protected females from food theft. We identify conditions favoring extractive foraging and food sharing across mating systems (i.e., monogamy, polygyny, promiscuity), and we assess which system maximizes female fitness with changes in the profitability of extractive foraging. Females extract foods and share them with males only when: i) extracting rather than collecting plant foods pays off energetically; and ii) males guard females. Males extract foods when they are sufficiently high in value, but share with females only under promiscuous mating and/or no mate guarding. These results suggest that if early hominins had mating systems with pair-bonds (monogamous or polygynous), then food sharing by adult females with unrelated adult males occurred before hunting, cooking, and extensive grandparenting. Such cooperation may have enabled early hominins to expand into more open, seasonal habitats, and provided a foundation for the subsequent evolution of human life histories.

T cell deletional tolerance restricts AQP4 but not MOG CNS autoimmunity.

Aquaporin-4 (AQP4)-specific Th17 cells are thought to have a central role in neuromyelitis optica (NMO) pathogenesis. When modeling NMO, only AQP4-reactive Th17 cells from AQP4-deficient (AQP4-/-), but not wild-type (WT) mice, caused CNS autoimmunity in recipient WT mice, indicating that a tightly regulated mechanism normally ensures tolerance to AQP4. Here, we found that pathogenic AQP4 T cell epitopes bind MHC II with exceptionally high affinity. Examination of T cell receptor (TCR) α/ß usage revealed that AQP4-specific T cells from AQP4-/- mice employed a distinct TCR repertoire and exhibited clonal expansion. Selective thymic AQP4 deficiency did not fully restore AQP4-reactive T cells, demonstrating that thymic negative selection alone did not account for AQP4-specific tolerance in WT mice. Indeed, AQP4-specific Th17 cells caused paralysis in recipient WT or B cell-deficient mice, which was followed by complete recovery that was associated with apoptosis of donor T cells. However, donor AQP4-reactive T cells survived and caused persistent paralysis in recipient mice deficient in both T and B cells or mice lacking T cells only. Thus, AQP4 CNS autoimmunity was limited by T cell-dependent deletion of AQP4-reactive T cells. In contrast, myelin oligodendrocyte glycoprotein (MOG)-specific T cells survived and caused sustained disease in WT mice. These findings underscore the importance of peripheral T cell deletional tolerance to AQP4, which may be relevant to understanding the balance of AQP4-reactive T cells in health and in NMO. T cell tolerance to AQP4, expressed in multiple tissues, is distinct from tolerance to MOG, an autoantigen restricted in its expression.

Critical assessment of variant prioritization methods for rare disease diagnosis within the rare genomes project.

BACKGROUND: A major obstacle faced by families with rare diseases is obtaining a genetic diagnosis. The average "diagnostic odyssey" lasts over five years and causal variants are identified in under 50%, even when capturing variants genome-wide. To aid in the interpretation and prioritization of the vast number of variants detected, computational methods are proliferating. Knowing which tools are most effective remains unclear. To evaluate the performance of computational methods, and to encourage innovation in method development, we designed a Critical Assessment of Genome Interpretation (CAGI) community challenge to place variant prioritization models head-to-head in a real-life clinical diagnostic setting. METHODS: We utilized genome sequencing (GS) data from families sequenced in the Rare Genomes Project (RGP), a direct-to-participant research study on the utility of GS for rare disease diagnosis and gene discovery. Challenge predictors were provided with a dataset of variant calls and phenotype terms from 175 RGP individuals (65 families), including 35 solved training set families with causal variants specified, and 30 unlabeled test set families (14 solved, 16 unsolved). We tasked teams to identify causal variants in as many families as possible. Predictors submitted variant predictions with estimated probability of causal relationship (EPCR) values. Model performance was determined by two metrics, a weighted score based on the rank position of causal variants, and the maximum F-measure, based on precision and recall of causal variants across all EPCR values. RESULTS: Sixteen teams submitted predictions from 52 models, some with manual review incorporated. Top performers recalled causal variants in up to 13 of 14 solved families within the top 5 ranked variants. Newly discovered diagnostic variants were returned to two previously unsolved families following confirmatory RNA sequencing, and two novel disease gene candidates were entered into Matchmaker Exchange. In one example, RNA sequencing demonstrated aberrant splicing due to a deep intronic indel in ASNS, identified in trans with a frameshift variant in an unsolved proband with phenotypes consistent with asparagine synthetase deficiency. CONCLUSIONS: Model methodology and performance was highly variable. Models weighing call quality, allele frequency, predicted deleteriousness, segregation, and phenotype were effective in identifying causal variants, and models open to phenotype expansion and non-coding variants were able to capture more difficult diagnoses and discover new diagnoses. Overall, computational models can significantly aid variant prioritization. For use in diagnostics, detailed review and conservative assessment of prioritized variants against established criteria is needed.

Conformational rigidity of cytochrome c'-α from a thermophile is associated with slow NO binding.

Cytochromes c'-α are nitric oxide (NO)-binding heme proteins derived from bacteria that can thrive in a wide range of temperature environments. Studies of mesophilic Alcaligenes xylosoxidans cytochrome c'-α (AxCP-α) have revealed an unusual NO-binding mechanism involving both heme faces, in which NO first binds to form a distal hexa-coordinate Fe(II)-NO (6cNO) intermediate and then displaces the proximal His to form a proximal penta-coordinate Fe(II)-NO (5cNO) final product. Here, we characterize a thermally stable cytochrome c'-α from thermophilic Hydrogenophilus thermoluteolus (PhCP-α) to understand how protein thermal stability affects NO binding. Electron paramagnetic and resonance Raman spectroscopies reveal the formation of a PhCP-α 5cNO product, with time-resolved (stopped-flow) UV-vis absorbance indicating the involvement of a 6cNO intermediate. Relative to AxCP-α, the rates of 6cNO and 5cNO formation in PhCP-α are ∼11- and ∼13-fold lower, respectively. Notably, x-ray crystal structures of PhCP-α in the presence and absence of NO suggest that the sluggish formation of the proximal 5cNO product results from conformational rigidity: the Arg-132 residue (adjacent to the proximal His ligand) is held in place by a salt bridge between Arg-75 and Glu-135 (an interaction not present in AxCP-α or a psychrophilic counterpart). Overall, our data provide fresh insights into structural factors controlling NO binding in heme proteins, including 5cNO complexes relevant to eukaryotic NO sensors.

Ventilator-induced Lung Injury Promotes Inflammation within the Pleural Cavity.

Mechanical ventilation contributes to the morbidity and mortality of patients in intensive care, likely through the exacerbation and dissemination of inflammation. Despite the proximity of the pleural cavity to the lungs and exposure to physical forces, little attention has been paid to its potential as an inflammatory source during ventilation. Here, we investigate the pleural cavity as a novel site of inflammation during ventilator-induced lung injury. Mice were subjected to low or high tidal volume ventilation strategies for up to 3 hours. Ventilation with a high tidal volume significantly increased cytokine and total protein levels in BAL and pleural lavage fluid. In contrast, acid aspiration, explored as an alternative model of injury, only promoted intraalveolar inflammation, with no effect on the pleural space. Resident pleural macrophages demonstrated enhanced activation after injurious ventilation, including upregulated ICAM-1 and IL-1ß expression, and the release of extracellular vesicles. In vivo ventilation and in vitro stretch of pleural mesothelial cells promoted ATP secretion, whereas purinergic receptor inhibition substantially attenuated extracellular vesicles and cytokine levels in the pleural space. Finally, labeled protein rapidly translocated from the pleural cavity into the circulation during high tidal volume ventilation, to a significantly greater extent than that of protein translocation from the alveolar space. Overall, we conclude that injurious ventilation induces pleural cavity inflammation mediated through purinergic pathway signaling and likely enhances the dissemination of mediators into the vasculature. This previously unidentified consequence of mechanical ventilation potentially implicates the pleural space as a focus of research and novel avenue for intervention in critical care.

Cardiac specification during gastrulation - The Yellow Brick Road leading to Tinman.

The question of how the heart develops, and the genetic networks governing this process have become intense areas of research over the past several decades. This research is propelled by classical developmental studies and potential clinical applications to understand and treat congenital conditions in which cardiac development is disrupted. Discovery of the tinman gene in Drosophila, and examination of its vertebrate homolog Nkx2.5, along with other core cardiac transcription factors has revealed how cardiac progenitor differentiation and maturation drives heart development. Careful observation of cardiac morphogenesis along with lineage tracing approaches indicated that cardiac progenitors can be divided into two broad classes of cells, namely the first and second heart fields, that contribute to the heart in two distinct waves of differentiation. Ample evidence suggests that the fate of individual cardiac progenitors is restricted to distinct cardiac structures quite early in development, well before the expression of canonical cardiac progenitor markers like Nkx2.5. Here we review the initial specification of cardiac progenitors, discuss evidence for the early patterning of cardiac progenitors during gastrulation, and consider how early gene expression programs and epigenetic patterns can direct their development. A complete understanding of when and how the developmental potential of cardiac progenitors is determined, and their potential plasticity, is of great interest developmentally and also has important implications for both the study of congenital heart disease and therapeutic approaches based on cardiac stem cell programming.