Characterizing the Major Structural Variant Alleles of the Human Genome.

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

CpG Island Hypermethylation Mediated by DNMT3A Is a Consequence of AML Progression.

DNMT3A mutations occur in ∼25% of acute myeloid leukemia (AML) patients. The most common mutation, DNMT3AR882H, has dominant negative activity that reduces DNA methylation activity by ∼80% in vitro. To understand the contribution of DNMT3A-dependent methylation to leukemogenesis, we performed whole-genome bisulfite sequencing of primary leukemic and non-leukemic cells in patients with or without DNMT3AR882 mutations. Non-leukemic hematopoietic cells with DNMT3AR882H displayed focal methylation loss, suggesting that hypomethylation antedates AML. Although virtually all AMLs with wild-type DNMT3A displayed CpG island hypermethylation, this change was not associated with gene silencing and was essentially absent in AMLs with DNMT3AR882 mutations. Primary hematopoietic stem cells expanded with cytokines were hypermethylated in a DNMT3A-dependent manner, suggesting that hypermethylation may be a response to, rather than a cause of, cellular proliferation. Our findings suggest that hypomethylation is an initiating phenotype in AMLs with DNMT3AR882, while DNMT3A-dependent CpG island hypermethylation is a consequence of AML progression.

Sequencing the mouse Y chromosome reveals convergent gene acquisition and amplification on both sex chromosomes.

We sequenced the MSY (male-specific region of the Y chromosome) of the C57BL/6J strain of the laboratory mouse Mus musculus. In contrast to theories that Y chromosomes are heterochromatic and gene poor, the mouse MSY is 99.9% euchromatic and contains about 700 protein-coding genes. Only 2% of the MSY derives from the ancestral autosomes that gave rise to the mammalian sex chromosomes. Instead, all but 45 of the MSY's genes belong to three acquired, massively amplified gene families that have no homologs on primate MSYs but do have acquired, amplified homologs on the mouse X chromosome. The complete mouse MSY sequence brings to light dramatic forces in sex chromosome evolution: lineage-specific convergent acquisition and amplification of X-Y gene families, possibly fueled by antagonism between acquired X-Y homologs. The mouse MSY sequence presents opportunities for experimental studies of a sex-specific chromosome in its entirety, in a genetically tractable model organism.

Vaccination with Glycan-Modified HIV NFL Envelope Trimer-Liposomes Elicits Broadly Neutralizing Antibodies to Multiple Sites of Vulnerability.

The elicitation of broadly neutralizing antibodies (bNAbs) against the HIV-1 envelope glycoprotein (Env) trimer remains a major vaccine challenge. Most cross-conserved protein determinants are occluded by self-N-glycan shielding, limiting B cell recognition of the underlying polypeptide surface. The exceptions to the contiguous glycan shield include the conserved receptor CD4 binding site (CD4bs) and glycoprotein (gp)41 elements proximal to the furin cleavage site. Accordingly, we performed heterologous trimer-liposome prime:boosting in rabbits to drive B cells specific for cross-conserved sites. To preferentially expose the CD4bs to B cells, we eliminated proximal N-glycans while maintaining the native-like state of the cleavage-independent NFL trimers, followed by gradual N-glycan restoration coupled with heterologous boosting. This approach successfully elicited CD4bs-directed, cross-neutralizing Abs, including one targeting a unique glycan-protein epitope and a bNAb (87% breadth) directed to the gp120:gp41 interface, both resolved by high-resolution cryoelectron microscopy. This study provides proof-of-principle immunogenicity toward eliciting bNAbs by vaccination.

The next-generation sequencing revolution and its impact on genomics.

Genomics is a relatively new scientific discipline, having DNA sequencing as its core technology. As technology has improved the cost and scale of genome characterization over sequencing's 40-year history, the scope of inquiry has commensurately broadened. Massively parallel sequencing has proven revolutionary, shifting the paradigm of genomics to address biological questions at a genome-wide scale. Sequencing now empowers clinical diagnostics and other aspects of medical care, including disease risk, therapeutic identification, and prenatal testing. This Review explores the current state of genomics in the massively parallel sequencing era.

The somatic genomic landscape of glioblastoma.

We describe the landscape of somatic genomic alterations based on multidimensional and comprehensive characterization of more than 500 glioblastoma tumors (GBMs). We identify several novel mutated genes as well as complex rearrangements of signature receptors, including EGFR and PDGFRA. TERT promoter mutations are shown to correlate with elevated mRNA expression, supporting a role in telomerase reactivation. Correlative analyses confirm that the survival advantage of the proneural subtype is conferred by the G-CIMP phenotype, and MGMT DNA methylation may be a predictive biomarker for treatment response only in classical subtype GBM. Integrative analysis of genomic and proteomic profiles challenges the notion of therapeutic inhibition of a pathway as an alternative to inhibition of the target itself. These data will facilitate the discovery of therapeutic and diagnostic target candidates, the validation of research and clinical observations and the generation of unanticipated hypotheses that can advance our molecular understanding of this lethal cancer.

Chromothripsis and human disease: piecing together the shattering process.

The unprecedented resolution of high-throughput genomics has enabled the recent discovery of a phenomenon by which specific regions of the genome are shattered and then stitched together via a single devastating event, referred to as chromothripsis. Potential mechanisms governing this process are now emerging, with implications for our understanding of the role of genomic rearrangements in development and disease.

Evolution of human-specific neural SRGAP2 genes by incomplete segmental duplication.

Gene duplication is an important source of phenotypic change and adaptive evolution. We leverage a haploid hydatidiform mole to identify highly identical sequences missing from the reference genome, confirming that the cortical development gene Slit-Robo Rho GTPase-activating protein 2 (SRGAP2) duplicated three times exclusively in humans. We show that the promoter and first nine exons of SRGAP2 duplicated from 1q32.1 (SRGAP2A) to 1q21.1 (SRGAP2B) ∼3.4 million years ago (mya). Two larger duplications later copied SRGAP2B to chromosome 1p12 (SRGAP2C) and to proximal 1q21.1 (SRGAP2D) ∼2.4 and ∼1 mya, respectively. Sequence and expression analyses show that SRGAP2C is the most likely duplicate to encode a functional protein and is among the most fixed human-specific duplicate genes. Our data suggest a mechanism where incomplete duplication created a novel gene function-antagonizing parental SRGAP2 function-immediately "at birth" 2-3 mya, which is a time corresponding to the transition from Australopithecus to Homo and the beginning of neocortex expansion.

Genomic landscape of non-small cell lung cancer in smokers and never-smokers.

We report the results of whole-genome and transcriptome sequencing of tumor and adjacent normal tissue samples from 17 patients with non-small cell lung carcinoma (NSCLC). We identified 3,726 point mutations and more than 90 indels in the coding sequence, with an average mutation frequency more than 10-fold higher in smokers than in never-smokers. Novel alterations in genes involved in chromatin modification and DNA repair pathways were identified, along with DACH1, CFTR, RELN, ABCB5, and HGF. Deep digital sequencing revealed diverse clonality patterns in both never-smokers and smokers. All validated EFGR and KRAS mutations were present in the founder clones, suggesting possible roles in cancer initiation. Analysis revealed 14 fusions, including ROS1 and ALK, as well as novel metabolic enzymes. Cell-cycle and JAK-STAT pathways are significantly altered in lung cancer, along with perturbations in 54 genes that are potentially targetable with currently available drugs.

Multiple hominin dispersals into Southwest Asia over the past 400,000 years.

Pleistocene hominin dispersals out of, and back into, Africa necessarily involved traversing the diverse and often challenging environments of Southwest Asia1-4. Archaeological and palaeontological records from the Levantine woodland zone document major biological and cultural shifts, such as alternating occupations by Homo sapiens and Neanderthals. However, Late Quaternary cultural, biological and environmental records from the vast arid zone that constitutes most of Southwest Asia remain scarce, limiting regional-scale insights into changes in hominin demography and behaviour1,2,5. Here we report a series of dated palaeolake sequences, associated with stone tool assemblages and vertebrate fossils, from the Khall Amayshan 4 and Jubbah basins in the Nefud Desert. These findings, including the oldest dated hominin occupations in Arabia, reveal at least five hominin expansions into the Arabian interior, coinciding with brief 'green' windows of reduced aridity approximately 400, 300, 200, 130-75 and 55 thousand years ago. Each occupation phase is characterized by a distinct form of material culture, indicating colonization by diverse hominin groups, and a lack of long-term Southwest Asian population continuity. Within a general pattern of African and Eurasian hominin groups being separated by Pleistocene Saharo-Arabian aridity, our findings reveal the tempo and character of climatically modulated windows for dispersal and admixture.

Author Correction: Elephant shark genome provides unique insights into gnathostome evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Key gp120 Glycans Pose Roadblocks to the Rapid Development of VRC01-Class Antibodies in an HIV-1-Infected Chinese Donor.

VRC01-class antibodies neutralize diverse HIV-1 strains by targeting the conserved CD4-binding site. Despite extensive investigations, crucial events in the early stage of VRC01 development remain elusive. We demonstrated how VRC01-class antibodies emerged in a Chinese donor by antigen-specific single B cell sorting, structural and functional studies, and longitudinal antibody and virus repertoire analyses. A monoclonal antibody DRVIA7 with modest neutralizing breadth was isolated that displayed a subset of VRC01 signatures. X-ray and EM structures revealed a VRC01-like angle of approach, but less favorable interactions between the DRVIA7 light-chain CDR1 and the N terminus with N276 and V5 glycans of gp120. Although the DRVIA7 lineage was unable to acquire broad neutralization, longitudinal analysis revealed a repertoire-encoded VRC01 light-chain CDR3 signature and VRC01-like neutralizing heavy-chain precursors that rapidly matured within 2 years. Thus, light chain accommodation of the glycan shield should be taken into account in vaccine design targeting this conserved site of vulnerability.

A human genome structural variation sequencing resource reveals insights into mutational mechanisms.

Understanding the prevailing mutational mechanisms responsible for human genome structural variation requires uniformity in the discovery of allelic variants and precision in terms of breakpoint delineation. We develop a resource based on capillary end sequencing of 13.8 million fosmid clones from 17 human genomes and characterize the complete sequence of 1054 large structural variants corresponding to 589 deletions, 384 insertions, and 81 inversions. We analyze the 2081 breakpoint junctions and infer potential mechanism of origin. Three mechanisms account for the bulk of germline structural variation: microhomology-mediated processes involving short (2-20 bp) stretches of sequence (28%), nonallelic homologous recombination (22%), and L1 retrotransposition (19%). The high quality and long-range continuity of the sequence reveals more complex mutational mechanisms, including repeat-mediated inversions and gene conversion, that are most often missed by other methods, such as comparative genomic hybridization, single nucleotide polymorphism microarrays, and next-generation sequencing.

Author Correction: Exome sequencing of Finnish isolates enhances rare-variant association power.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

Exome sequencing of Finnish isolates enhances rare-variant association power.

Exome-sequencing studies have generally been underpowered to identify deleterious alleles with a large effect on complex traits as such alleles are mostly rare. Because the population of northern and eastern Finland has expanded considerably and in isolation following a series of bottlenecks, individuals of these populations have numerous deleterious alleles at a relatively high frequency. Here, using exome sequencing of nearly 20,000 individuals from these regions, we investigate the role of rare coding variants in clinically relevant quantitative cardiometabolic traits. Exome-wide association studies for 64 quantitative traits identified 26 newly associated deleterious alleles. Of these 26 alleles, 19 are either unique to or more than 20 times more frequent in Finnish individuals than in other Europeans and show geographical clustering comparable to Mendelian disease mutations that are characteristic of the Finnish population. We estimate that sequencing studies of populations without this unique history would require hundreds of thousands to millions of participants to achieve comparable association power.

The Impact of Inflammation on Thermal Hyperpnea: Relevance for Heat Stress and Febrile Seizures.

Extreme heat caused by climate change is increasing transmission of infectious diseases resulting in a sharp rise in heat-related illness and mortality. Understanding mechanistic link between heat, inflammation and disease is thus important for public health. Thermal hyperpnea, and consequent respiratory alkalosis is crucial in febrile seizures and convulsions induced by heat stress in humans. Here we address what causes thermal hyperpnea in neonates and how is it affected by inflammation. TRPV1, a heat-activated channel is sensitized by inflammation and modulates breathing, and thus may play a key role. To investigate whether inflammatory sensitization of TRPV1 modifies neonatal ventilatory responses to heat stress, leading to respiratory alkalosis and an increased susceptibility to hyperthermic seizures we treated neonatal rats with bacterial lipopolysaccharide, and breathing, arterial pH, in-vitro vagus nerve activity, and seizure susceptibility were assessed during heat stress in the presence or absence of a TRPV1 antagonist (AMG-9810) or shRNA-mediated TRPV1 suppression. Lipopolysaccharide-induced inflammatory preconditioning lowered the threshold temperature and latency of hyperthermic seizures. This was accompanied by increased tidal volume, minute ventilation, expired CO2, and arterial pH (alkalosis). Lipopolysaccharide exposure also elevated vagal spiking and intracellular calcium levels in response to hyperthermia. TRPV1 inhibition with AMG-9810 or shRNA reduced the lipopolysaccharide-induced susceptibility to hyperthermic seizures and altered the breathing pattern to fast shallow breaths (tachypnea), making each breath less efficient and restoring arterial pH. These results indicate that inflammation exacerbates thermal hyperpnea-induced respiratory alkalosis associated with increased susceptibility to hyperthermic seizures, primarily mediated by TRPV1 localized to vagus neurons. This article is open access and distributed under the terms of the Creative Commons Attribution Non-Commercial No Derivatives License 4.0 (http://creativecommons.org/licenses/by-nc-nd/4.0/).

Full-length isoform concatenation sequencing to resolve cancer transcriptome complexity.

BACKGROUND: Cancers exhibit complex transcriptomes with aberrant splicing that induces isoform-level differential expression compared to non-diseased tissues. Transcriptomic profiling using short-read sequencing has utility in providing a cost-effective approach for evaluating isoform expression, although short-read assembly displays limitations in the accurate inference of full-length transcripts. Long-read RNA sequencing (Iso-Seq), using the Pacific Biosciences (PacBio) platform, can overcome such limitations by providing full-length isoform sequence resolution which requires no read assembly and represents native expressed transcripts. A constraint of the Iso-Seq protocol is due to fewer reads output per instrument run, which, as an example, can consequently affect the detection of lowly expressed transcripts. To address these deficiencies, we developed a concatenation workflow, PacBio Full-Length Isoform Concatemer Sequencing (PB_FLIC-Seq), designed to increase the number of unique, sequenced PacBio long-reads thereby improving overall detection of unique isoforms. In addition, we anticipate that the increase in read depth will help improve the detection of moderate to low-level expressed isoforms. RESULTS: In sequencing a commercial reference (Spike-In RNA Variants; SIRV) with known isoform complexity we demonstrated a 3.4-fold increase in read output per run and improved SIRV recall when using the PB_FLIC-Seq method compared to the same samples processed with the Iso-Seq protocol. We applied this protocol to a translational cancer case, also demonstrating the utility of the PB_FLIC-Seq method for identifying differential full-length isoform expression in a pediatric diffuse midline glioma compared to its adjacent non-malignant tissue. Our data analysis revealed increased expression of extracellular matrix (ECM) genes within the tumor sample, including an isoform of the Secreted Protein Acidic and Cysteine Rich (SPARC) gene that was expressed 11,676-fold higher than in the adjacent non-malignant tissue. Finally, by using the PB_FLIC-Seq method, we detected several cancer-specific novel isoforms. CONCLUSION: This work describes a concatenation-based methodology for increasing the number of sequenced full-length isoform reads on the PacBio platform, yielding improved discovery of expressed isoforms. We applied this workflow to profile the transcriptome of a pediatric diffuse midline glioma and adjacent non-malignant tissue. Our findings of cancer-specific novel isoform expression further highlight the importance of long-read sequencing for characterization of complex tumor transcriptomes.

Sirt5-mediated lysine desuccinylation regulates oxidative stress adaptation in Magnaporthe oryzae during host intracellular infection.

Plant pathogenic fungi elaborate numerous detoxification strategies to suppress host reactive oxygen species (ROS), but their coordination is not well-understood. Here, we show that Sirt5-mediated protein desuccinylation in Magnaporthe oryzae is central to host ROS detoxification. SIRT5 encodes a desuccinylase important for virulence via adaptation to host oxidative stress. Quantitative proteomics analysis identified a large number of succinylated proteins targeted by Sirt5, most of which were mitochondrial proteins involved in oxidative phosphorylation, TCA cycle, and fatty acid oxidation. Deletion of SIRT5 resulted in hypersuccinylation of detoxification-related enzymes, and significant reduction in NADPH : NADP+ and GSH : GSSG ratios, disrupting redox balance and impeding invasive growth. Sirt5 desuccinylated thioredoxin Trx2 and glutathione peroxidase Hyr1 to activate their enzyme activity, likely by affecting proper folding. Altogether, this work demonstrates the importance of Sirt5-mediated desuccinylation in controlling fungal process required for detoxifying host ROS during M. oryzae infection.

Cerebral organoids containing an AUTS2 missense variant model microcephaly.

Variants in the AUTS2 gene are associated with a broad spectrum of neurological conditions characterized by intellectual disability, microcephaly, and congenital brain malformations. Here, we use a human cerebral organoid model to investigate the pathophysiology of a heterozygous de novo missense AUTS2 variant identified in a patient with multiple neurological impairments including primary microcephaly and profound intellectual disability. Proband cerebral organoids exhibit reduced growth, deficits in neural progenitor cell (NPC) proliferation and disrupted NPC polarity within ventricular zone-like regions compared to control cerebral organoids. We used CRISPR-Cas9-mediated gene editing to correct this variant and demonstrate rescue of impaired organoid growth and NPC proliferative deficits. Single-cell RNA sequencing revealed a marked reduction of G1/S transition gene expression and alterations in WNT-ß-catenin signalling within proband NPCs, uncovering a novel role for AUTS2 in NPCs during human cortical development. Collectively, these results underscore the value of cerebral organoids to investigate molecular mechanisms underlying AUTS2 syndrome.

Solvent effects on extractant conformational energetics in liquid-liquid extraction: a simulation study of molecular solvents and ionic liquids.

Extractant design in liquid-liquid extraction (LLE) is a research frontier of metal ion separations that typically focuses on the direct extractant-metal interactions. However, a more detailed understanding of energetic drivers of separations beyond primary metal coordination is often lacking, including the role of solvent in the extractant phase. In this work, we propose a new mechanism for enhancing metal-complexant energetics with nanostructured solvents. Using molecular dynamics simulations with umbrella sampling, we find that the organic solvent can reshape the energetics of the extractant's intramolecular conformational landscape. We calculate free energy profiles of different conformations of a representative bidentate extractant, n-octyl(phenyl)-N,N-diisobutyl carbamoyl methyl phosphinoxide (CMPO), in four different solvents: dodecane, tributyl phosphate (TBP), and dry and wet ionic liquid (IL) 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf2N]). By promoting reorganization of the extractant molecule into its binding conformation, our findings reveal how particular solvents can ameliorate this unfavorable step of the metal separation process. In particular, the charge alternating nanodomains formed in ILs substantially reduce the free energy penalty associated with extractant reorganization. Importantly, using alchemical free energy calculations, we find that this stabilization persists even when we explicitly include the extracted cation. These findings provide insight into the energetic drivers of metal ion separations and potentially suggest a new approach to designing effective separations using a molecular-level understanding of solvent effects.