Explore the vast range of titles available.

Regulated transcription controls the diversity, developmental pathways and spatial organization of the hundreds of cell types that make up a mammal. Using single-molecule cDNA sequencing, we mapped transcription start sites (TSSs) and their usage in human and mouse primary cells, cell lines and tissues to produce a comprehensive overview of mammalian gene expression across the human body. We find that few genes are truly ‘housekeeping’, whereas many mammalian promoters are composite entities composed of several closely separated TSSs, with independent cell-type-specific expression profiles. TSSs specific to different cell types evolve at different rates, whereas promoters of broadly expressed genes are the most conserved. Promoter-based expression analysis reveals key transcription factors defining cell states and links them to binding-site motifs. The functions of identified novel transcripts can be predicted by coexpression and sample ontology enrichment analyses. The functional annotation of the mammalian genome 5 (FANTOM5) project provides comprehensive expression profiles and functional annotation of mammalian cell-type-specific transcriptomes with wide applications in biomedical research. A study from the FANTOM consortium using single-molecule cDNA sequencing of transcription start sites and their usage in human and mouse primary cells, cell lines and tissues reveals insights into the specificity and diversity of transcription patterns across different mammalian cell types. Mapping the human transcription FANTOM5 (standing for functional annotation of the mammalian genome 5) is the fifth major stage of a major international collaboration that aims to dissect the transcriptional regulatory networks that define every human cell type. Two Articles in this issue of Nature present some of the project's latest results. The first paper uses the FANTOM5 panel of tissue and primary cell samples to define an atlas of active, in vivo bidirectionally transcribed enhancers across the human body. These authors show that bidirectional capped RNAs are a signature feature of active enhancers and identify more than 40,000 enhancer candidates from over 800 human cell and tissue samples. The enhancer atlas is used to compare regulatory programs between different cell types and identify disease-associated regulatory SNPs, and will be a resource for studies on cell-type-specific enhancers. In the second paper, single-molecule sequencing is used to map human and mouse transcription start sites and their usage in a panel of distinct human and mouse primary cells, cell lines and tissues to produce the most comprehensive mammalian gene expression atlas to date. The data provide a plethora of insights into open reading frames and promoters across different cell types in addition to valuable annotation of mammalian cell-type-specific transcriptomes.

Every organism has a different set of genes essential for its viability. This indicates that an organism can become tolerant to the loss of an essential gene under certain circumstances during evolution, via the manifestation of ‘masked’ alternative mechanisms. In our quest to systematically uncover masked mechanisms in eukaryotic cells, we developed an extragenic suppressor screening method using haploid spores deleted of an essential gene in the fission yeast Schizosaccharomyces pombe. We screened for the ‘bypass’ suppressors of lethality of 92 randomly selected genes that are essential for viability in standard laboratory culture conditions. Remarkably, extragenic mutations bypassed the essentiality of as many as 20 genes (22%), 15 of which have not been previously reported. Half of the bypass-suppressible genes were involved in mitochondria function; we also identified multiple genes regulating RNA processing. 18 suppressible genes were conserved in the budding yeast Saccharomyces cerevisiae, but 13 of them were non-essential in that species. These trends suggest that essentiality bypass is not a rare event and that each organism may be endowed with secondary or backup mechanisms that can substitute for primary mechanisms in various biological processes. Furthermore, the robustness of our simple spore-based methodology paves the way for genome-scale screening.Key words: Schizosaccharomyces pombe, extragenic suppressor screening, bypass of essentiality (BOE), cut7 (kinesin-5), hul5 (E3 ubiquitin ligase)

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 sizes 1 . 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. A catalogue of predicted loss-of-function variants in 125,748 whole-exome and 15,708 whole-genome sequencing datasets from the Genome Aggregation Database (gnomAD) reveals the spectrum of mutational constraints that affect these human protein-coding genes.

Since the first half of the twentieth century, evolutionary theory has been dominated by the idea that mutations occur randomly with respect to their consequences 1 . Here we test this assumption with large surveys of de novo mutations in the plant Arabidopsis thaliana . In contrast to expectations, we find that mutations occur less often in functionally constrained regions of the genome—mutation frequency is reduced by half inside gene bodies and by two-thirds in essential genes. With independent genomic mutation datasets, including from the largest Arabidopsis mutation accumulation experiment conducted to date, we demonstrate that epigenomic and physical features explain over 90% of variance in the genome-wide pattern of mutation bias surrounding genes. Observed mutation frequencies around genes in turn accurately predict patterns of genetic polymorphisms in natural Arabidopsis accessions ( r  = 0.96). That mutation bias is the primary force behind patterns of sequence evolution around genes in natural accessions is supported by analyses of allele frequencies. Finally, we find that genes subject to stronger purifying selection have a lower mutation rate. We conclude that epigenome-associated mutation bias 2 reduces the occurrence of deleterious mutations in Arabidopsis , challenging the prevailing paradigm that mutation is a directionless force in evolution. Data on de novo mutations in Arabidopsis thaliana reveal that mutations do not occur randomly; instead, epigenome-associated mutation bias reduces the occurrence of deleterious mutations.

Salamanders serve as important tetrapod models for developmental, regeneration and evolutionary studies. An extensive molecular toolkit makes the Mexican axolotl ( Ambystoma mexicanum ) a key representative salamander for molecular investigations. Here we report the sequencing and assembly of the 32-gigabase-pair axolotl genome using an approach that combined long-read sequencing, optical mapping and development of a new genome assembler (MARVEL). We observed a size expansion of introns and intergenic regions, largely attributable to multiplication of long terminal repeat retroelements. We provide evidence that intron size in developmental genes is under constraint and that species-restricted genes may contribute to limb regeneration. The axolotl genome assembly does not contain the essential developmental gene Pax3 . However, mutation of the axolotl Pax3 paralogue Pax7 resulted in an axolotl phenotype that was similar to those seen in Pax3 −/− and Pax7 −/− mutant mice. The axolotl genome provides a rich biological resource for developmental and evolutionary studies. Sequencing and assembly of the 32-Gb genome of the Mexican axolotl reveals that it lacks the developmental gene Pax3 , which is essential in other vertebrates; the genome sequence could improve our understanding of the evolution of the axolotl’s remarkable regenerative capabilities. Axolotl genome sequence Elly Tanaka, Eugene Myers and colleagues report the genome sequence of the axolotl, a model organism for developmental, regeneration and evolutionary studies. To sequence and assemble this large and complex genome, the authors used a combination of long- and short-read sequencing, optical mapping and a new genome assembly pipeline, MARVEL, optimized for long-read sequencing of complex genomes. The genome assembly shows an expansion of long terminal repeat retroelements and the presence of a large HoxA cluster, but also a reduction in the number of Pax-family genes in the genome of this popular salamander.

Naturally occurring human genetic variants that are predicted to inactivate protein-coding genes provide an in vivo model of human gene inactivation that complements knockout studies in cells and model organisms. Here we report three key findings regarding the assessment of candidate drug targets using human loss-of-function variants. First, even essential genes, in which loss-of-function variants are not tolerated, can be highly successful as targets of inhibitory drugs. Second, in most genes, loss-of-function variants are sufficiently rare that genotype-based ascertainment of homozygous or compound heterozygous ‘knockout’ humans will await sample sizes that are approximately 1,000 times those presently available, unless recruitment focuses on consanguineous individuals. Third, automated variant annotation and filtering are powerful, but manual curation remains crucial for removing artefacts, and is a prerequisite for recall-by-genotype efforts. Our results provide a roadmap for human knockout studies and should guide the interpretation of loss-of-function variants in drug development. Analysis of predicted loss-of-function variants from 125,748 human exomes and 15,708 whole genomes in the Genome Aggregation Database (gnomAD) provides a roadmap for human ‘knockout’ studies and a guide for future research into disease biology and drug-target selection.

Variants of uncertain significance fundamentally limit the clinical utility of genetic information. The challenge they pose is epitomized by BRCA1 , a tumour suppressor gene in which germline loss-of-function variants predispose women to breast and ovarian cancer. Although BRCA1 has been sequenced in millions of women, the risk associated with most newly observed variants cannot be definitively assigned. Here we use saturation genome editing to assay 96.5% of all possible single-nucleotide variants (SNVs) in 13 exons that encode functionally critical domains of BRCA1. Functional effects for nearly 4,000 SNVs are bimodally distributed and almost perfectly concordant with established assessments of pathogenicity. Over 400 non-functional missense SNVs are identified, as well as around 300 SNVs that disrupt expression. We predict that these results will be immediately useful for the clinical interpretation of BRCA1 variants, and that this approach can be extended to overcome the challenge of variants of uncertain significance in additional clinically actionable genes. Germline BRCA1 loss-of-function variants are associated with predisposition to early-onset breast and ovarian cancer; here the authors use CRISPR/Cas9 genome editing to functionally assess thousands of BRCA1 variants in order to facilitate the clinical interpretation of these variants.

Cryptophyte and chlorarachniophyte algae are transitional forms in the widespread secondary endosymbiotic acquisition of photosynthesis by engulfment of eukaryotic algae. Unlike most secondary plastid-bearing algae, miniaturized versions of the endosymbiont nuclei (nucleomorphs) persist in cryptophytes and chlorarachniophytes. To determine why, and to address other fundamental questions about eukaryote–eukaryote endosymbiosis, we sequenced the nuclear genomes of the cryptophyte Guillardia theta and the chlorarachniophyte Bigelowiella natans . Both genomes have >21,000 protein genes and are intron rich, and B. natans exhibits unprecedented alternative splicing for a single-celled organism. Phylogenomic analyses and subcellular targeting predictions reveal extensive genetic and biochemical mosaicism, with both host- and endosymbiont-derived genes servicing the mitochondrion, the host cell cytosol, the plastid and the remnant endosymbiont cytosol of both algae. Mitochondrion-to-nucleus gene transfer still occurs in both organisms but plastid-to-nucleus and nucleomorph-to-nucleus transfers do not, which explains why a small residue of essential genes remains locked in each nucleomorph. Sequencing the nuclear genomes of Guillardia theta and Bigelowiella natans , transitional forms in the endosymbiotic acquisition of photosynthesis by engulfment of certain eukaryotic algae, reveals unprecedented alternative splicing for a single-celled organism ( B. natans ) and extensive genetic and biochemical mosaicism, shedding light on why nucleomorphs persist in these species but not other algae. Evolutionarily complex algal genomes revealed This paper presents the sequences of the nuclear genomes of two eukaryotic microbes of remarkable genetic and cellular complexity, Guillardia and Bigelowiella . These algae are transitional forms in the endosymbiotic acquisition of photosynthesis by engulfment of eukaryotic algae, and possess four genomes: mitochondrial and plastid (chloroplast) genomes, a nuclear genome of host origin and a miniaturized 'nucleomorph' genome of endosymbiotic origin. Analyses reveal unprecedented alternative splicing for a single-celled organism, and extensive genetic and biochemical mosaicism. Whereas the mitochondrion-to-nucleus gene transfer continues in both organisms, plastid-to-nucleus and nucleomorph-to-nucleus transfers have ceased, explaining nucleomorph persistence.

The core promoters of developmental and housekeeping genes are shown to have distinct specificities for different enhancer sequences in Drosophila , and this specificity separates developmental and housekeeping gene regulatory programs across the genome. Separate control of development and housekeeping genes The transcription of genes during development involves interactions between RNA polymerase II bound to core promoters, and factors bound to distal enhancers that confer cell-type-specific activation. Alexander Stark and colleagues show in Drosophila that the enhancers of developmental and housekeeping genes have distinct specificities for different core promoter sequences, and that this sequence-encoded enhancer–core-promoter specificity separates developmental and housekeeping gene regulatory programs across the genome. Gene transcription in animals involves the assembly of RNA polymerase II at core promoters and its cell-type-specific activation by enhancers that can be located more distally 1 . However, how ubiquitous expression of housekeeping genes is achieved has been less clear. In particular, it is unknown whether ubiquitously active enhancers exist and how developmental and housekeeping gene regulation is separated. An attractive hypothesis is that different core promoters might exhibit an intrinsic specificity to certain enhancers 2 , 3 , 4 , 5 , 6 . This is conceivable, as various core promoter sequence elements are differentially distributed between genes of different functions 7 , including elements that are predominantly found at either developmentally regulated or at housekeeping genes 8 , 9 , 10 . Here we show that thousands of enhancers in Drosophila melanogaster S2 and ovarian somatic cells (OSCs) exhibit a marked specificity to one of two core promoters—one derived from a ubiquitously expressed ribosomal protein gene and another from a developmentally regulated transcription factor—and confirm the existence of these two classes for five additional core promoters from genes with diverse functions. Housekeeping enhancers are active across the two cell types, while developmental enhancers exhibit strong cell-type specificity. Both enhancer classes differ in their genomic distribution, the functions of neighbouring genes, and the core promoter elements of these neighbouring genes. In addition, we identify two transcription factors—Dref and Trl—that bind and activate housekeeping versus developmental enhancers, respectively. Our results provide evidence for a sequence-encoded enhancer–core-promoter specificity that separates developmental and housekeeping gene regulatory programs for thousands of enhancers and their target genes across the entire genome.

Nature uses 64 codons to encode the synthesis of proteins from the genome, and chooses 1 sense codon—out of up to 6 synonyms—to encode each amino acid. Synonymous codon choice has diverse and important roles, and many synonymous substitutions are detrimental. Here we demonstrate that the number of codons used to encode the canonical amino acids can be reduced, through the genome-wide substitution of target codons by defined synonyms. We create a variant of Escherichia coli with a four-megabase synthetic genome through a high-fidelity convergent total synthesis. Our synthetic genome implements a defined recoding and refactoring scheme—with simple corrections at just seven positions—to replace every known occurrence of two sense codons and a stop codon in the genome. Thus, we recode 18,214 codons to create an organism with a 61-codon genome; this organism uses 59 codons to encode the 20 amino acids, and enables the deletion of a previously essential transfer RNA. High-fidelity convergent total synthesis is used to produce Escherichia coli with a 61-codon synthetic genome that uses 59 codons to encode all of the canonical amino acids.