Embryonic development is driven by tightly regulated patterns of gene expression, despite extensive genetic variation among individuals. Studies of expression quantitative trait loci (eQTL) indicate that genetic variation frequently alters gene expression in cell-culture models and differentiated tissues. However, the extent and types of genetic variation impacting embryonic gene expression, and their interactions with developmental programs, remain largely unknown. Here we assessed the effect of genetic variation on transcriptional (expression levels) and post-transcriptional (3' RNA processing) regulation across multiple stages of metazoan development, using 80 inbred Drosophila wild isolates, identifying thousands of developmental-stage-specific and shared QTL. Given the small blocks of linkage disequilibrium in Drosophila, we obtain near base-pair resolution, resolving causal mutations in developmental enhancers, validated transcription-factor-binding sites and RNA motifs. This fine-grain mapping uncovered extensive allelic interactions within enhancers that have opposite effects, thereby buffering their impact on enhancer activity. QTL affecting 3' RNA processing identify new functional motifs leading to transcript isoform diversity and changes in the lengths of 3' untranslated regions. These results highlight how developmental stage influences the effects of genetic variation and uncover multiple mechanisms that regulate and buffer expression variation during embryogenesis.

Lymphatic vessels are lined by lymphatic endothelial cells (LECs), and are critical for health. However, the role of metabolism in lymphatic development has not yet been elucidated. Here we report that in transgenic mouse models, LEC-specific loss of CPT1A, a rate-controlling enzyme in fatty acid β-oxidation, impairs lymphatic development. LECs use fatty acid β-oxidation to proliferate and for epigenetic regulation of lymphatic marker expression during LEC differentiation. Mechanistically, the transcription factor PROX1 upregulates CPT1A expression, which increases acetyl coenzyme A production dependent on fatty acid β-oxidation. Acetyl coenzyme A is used by the histone acetyltransferase p300 to acetylate histones at lymphangiogenic genes. PROX1-p300 interaction facilitates preferential histone acetylation at PROX1-target genes. Through this metabolism-dependent mechanism, PROX1 mediates epigenetic changes that promote lymphangiogenesis. Notably, blockade of CPT1 enzymes inhibits injury-induced lymphangiogenesis, and replenishing acetyl coenzyme A by supplementing acetate rescues this process in vivo.

Ash trees (genus Fraxinus, family Oleaceae) are widespread throughout the Northern Hemisphere, but are being devastated in Europe by the fungus Hymenoscyphus fraxineus, causing ash dieback, and in North America by the herbivorous beetle Agrilus planipennis. Here we sequence the genome of a low-heterozygosity Fraxinus excelsior tree from Gloucestershire, UK, annotating 38,852 protein-coding genes of which 25% appear ash specific when compared with the genomes of ten other plant species. Analyses of paralogous genes suggest a whole-genome duplication shared with olive (Olea europaea, Oleaceae). We also re-sequence 37 F. excelsior trees from Europe, finding evidence for apparent long-term decline in effective population size. Using our reference sequence, we re-analyse association transcriptomic data, yielding improved markers for reduced susceptibility to ash dieback. Surveys of these markers in British populations suggest that reduced susceptibility to ash dieback may be more widespread in Great Britain than in Denmark. We also present evidence that susceptibility of trees to H. fraxineus is associated with their iridoid glycoside levels. This rapid, integrated, multidisciplinary research response to an emerging health threat in a non-model organism opens the way for mitigation of the epidemic.

Chromatin structure at the length scale encompassing local nucleosome-nucleosome interactions is thought to play a crucial role in regulating transcription and access to DNA. However, this secondary structure of chromatin remains poorly understood compared with the primary structure of single nucleosomes or the tertiary structure of long-range looping interactions. Here we report the first genome-wide map of chromatin conformation in human cells at the 1-3 nucleosome (50-500 bp) scale, obtained using ionizing radiation-induced spatially correlated cleavage of DNA with sequencing (RICC-seq) to identify DNA-DNA contacts that are spatially proximal. Unbiased analysis of RICC-seq signal reveals regional enrichment of DNA fragments characteristic of alternating rather than adjacent nucleosome interactions in tri-nucleosome units, particularly in H3K9me3-marked heterochromatin. We infer differences in the likelihood of nucleosome-nucleosome contacts among open chromatin, H3K27me3-marked, and H3K9me3-marked repressed chromatin regions. After calibrating RICC-seq signal to three-dimensional distances, we show that compact two-start helical fibre structures with stacked alternating nucleosomes are consistent with RICC-seq fragmentation patterns from H3K9me3-marked chromatin, while non-compact structures and solenoid structures are consistent with open chromatin. Our data support a model of chromatin architecture in intact interphase nuclei consistent with variable longitudinal compaction of two-start helical fibres.

Although long non-coding RNAs (lncRNAs) are non-protein-coding transcripts by definition, recent studies have shown that a fraction of putative small open reading frames within lncRNAs are translated. However, the biological significance of these hidden polypeptides is still unclear. Here we identify and functionally characterize a novel polypeptide encoded by the lncRNA LINC00961. This polypeptide is conserved between human and mouse, is localized to the late endosome/lysosome and interacts with the lysosomal v-ATPase to negatively regulate mTORC1 activation. This regulation of mTORC1 is specific to activation of mTORC1 by amino acid stimulation, rather than by growth factors. Hence, we termed this polypeptide 'small regulatory polypeptide of amino acid response' (SPAR). We show that the SPAR-encoding lncRNA is highly expressed in a subset of tissues and use CRISPR/Cas9 engineering to develop a SPAR-polypeptide-specific knockout mouse while maintaining expression of the host lncRNA. We find that the SPAR-encoding lncRNA is downregulated in skeletal muscle upon acute injury, and using this in vivo model we establish that SPAR downregulation enables efficient activation of mTORC1 and promotes muscle regeneration. Our data provide a mechanism by which mTORC1 activation may be finely regulated in a tissue-specific manner in response to injury, and a paradigm by which lncRNAs encoding small polypeptides can modulate general biological pathways and processes to facilitate tissue-specific requirements, consistent with their restricted and highly regulated expression profile.

The spectrum of the hydrogen atom has played a central part in fundamental physics over the past 200 years. Historical examples of its importance include the wavelength measurements of absorption lines in the solar spectrum by Fraunhofer, the identification of transition lines by Balmer, Lyman and others, the empirical description of allowed wavelengths by Rydberg, the quantum model of Bohr, the capability of quantum electrodynamics to precisely predict transition frequencies, and modern measurements of the 1S-2S transition by Hänsch to a precision of a few parts in 1015. Recent technological advances have allowed us to focus on antihydrogen-the antimatter equivalent of hydrogen. The Standard Model predicts that there should have been equal amounts of matter and antimatter in the primordial Universe after the Big Bang, but today's Universe is observed to consist almost entirely of ordinary matter. This motivates the study of antimatter, to see if there is a small asymmetry in the laws of physics that govern the two types of matter. In particular, the CPT (charge conjugation, parity reversal and time reversal) theorem, a cornerstone of the Standard Model, requires that hydrogen and antihydrogen have the same spectrum. Here we report the observation of the 1S-2S transition in magnetically trapped atoms of antihydrogen. We determine that the frequency of the transition, which is driven by two photons from a laser at 243 nanometres, is consistent with that expected for hydrogen in the same environment. This laser excitation of a quantum state of an atom of antimatter represents the most precise measurement performed on an anti-atom. Our result is consistent with CPT invariance at a relative precision of about 2 × 10-10.

CRISPR-Cas systems provide microbes with adaptive immunity by employing short DNA sequences, termed spacers, that guide Cas proteins to cleave foreign DNA. Class 2 CRISPR-Cas systems are streamlined versions, in which a single RNA-bound Cas protein recognizes and cleaves target sequences. The programmable nature of these minimal systems has enabled researchers to repurpose them into a versatile technology that is broadly revolutionizing biological and clinical research. However, current CRISPR-Cas technologies are based solely on systems from isolated bacteria, leaving the vast majority of enzymes from organisms that have not been cultured untapped. Metagenomics, the sequencing of DNA extracted directly from natural microbial communities, provides access to the genetic material of a huge array of uncultivated organisms. Here, using genome-resolved metagenomics, we identify a number of CRISPR-Cas systems, including the first reported Cas9 in the archaeal domain of life, to our knowledge. This divergent Cas9 protein was found in little-studied nanoarchaea as part of an active CRISPR-Cas system. In bacteria, we discovered two previously unknown systems, CRISPR-CasX and CRISPR-CasY, which are among the most compact systems yet discovered. Notably, all required functional components were identified by metagenomics, enabling validation of robust in vivo RNA-guided DNA interference activity in Escherichia coli. Interrogation of environmental microbial communities combined with in vivo experiments allows us to access an unprecedented diversity of genomes, the content of which will expand the repertoire of microbe-based biotechnologies.

Respiratory dysfunction is a notorious cause of perinatal mortality in infants and sleep apnoea in adults, but the mechanisms of respiratory control are not clearly understood. Mechanical signals transduced by airway-innervating sensory neurons control respiration; however, the physiological significance and molecular mechanisms of these signals remain obscured. Here we show that global and sensory neuron-specific ablation of the mechanically activated ion channel Piezo2 causes respiratory distress and death in newborn mice. Optogenetic activation of Piezo2+ vagal sensory neurons causes apnoea in adult mice. Moreover, induced ablation of Piezo2 in sensory neurons of adult mice causes decreased neuronal responses to lung inflation, an impaired Hering-Breuer mechanoreflex, and increased tidal volume under normal conditions. These phenotypes are reproduced in mice lacking Piezo2 in the nodose ganglion. Our data suggest that Piezo2 is an airway stretch sensor and that Piezo2-mediated mechanotransduction within various airway-innervating sensory neurons is critical for establishing efficient respiration at birth and maintaining normal breathing in adults.

CLC proteins transport chloride (Cl-) ions across cellular membranes to regulate muscle excitability, electrolyte movement across epithelia, and acidification of intracellular organelles. Some CLC proteins are channels that conduct Cl- ions passively, whereas others are secondary active transporters that exchange two Cl- ions for one H+. The structural basis underlying these distinctive transport mechanisms is puzzling because CLC channels and transporters are expected to share the same architecture on the basis of sequence homology. Here we determined the structure of a bovine CLC channel (CLC-K) using cryo-electron microscopy. A conserved loop in the Cl- transport pathway shows a structure markedly different from that of CLC transporters. Consequently, the cytosolic constriction for Cl- passage is widened in CLC-K such that the kinetic barrier previously postulated for Cl-/H+ transporter function would be reduced. Thus, reduction of a kinetic barrier in CLC channels enables fast flow of Cl- down its electrochemical gradient.

Engineered spin-orbit coupling (SOC) in cold-atom systems can enable the study of new synthetic materials and complex condensed matter phenomena. However, spontaneous emission in alkali-atom spin-orbit-coupled systems is hindered by heating, limiting the observation of many-body effects and motivating research into potential alternatives. Here we demonstrate that spin-orbit-coupled fermions can be engineered to occur naturally in a one-dimensional optical lattice clock. In contrast to previous SOC experiments, here the SOC is both generated and probed using a direct ultra-narrow optical clock transition between two electronic orbital states in 87Sr atoms. We use clock spectroscopy to prepare lattice band populations, internal electronic states and quasi-momenta, and to produce spin-orbit-coupled dynamics. The exceptionally long lifetime of the excited clock state (160 seconds) eliminates decoherence and atom loss from spontaneous emission at all relevant experimental timescales, allowing subsequent momentum- and spin-resolved in situ probing of the SOC band structure and eigenstates. We use these capabilities to study Bloch oscillations, spin-momentum locking and Van Hove singularities in the transition density of states. Our results lay the groundwork for using fermionic optical lattice clocks to probe new phases of matter.

Phosphorus is an important nutrient for crop productivity. More than 60% of the total phosphorus in cereal crops is finally allocated into the grains and is therefore removed at harvest. This removal accounts for 85% of the phosphorus fertilizers applied to the field each year. However, because humans and non-ruminants such as poultry, swine and fish cannot digest phytate, the major form of phosphorus in the grains, the excreted phosphorus causes eutrophication of waterways. A reduction in phosphorus accumulation in the grain would contribute to sustainable and environmentally friendly agriculture. Here we describe a rice transporter, SULTR-like phosphorus distribution transporter (SPDT), that controls the allocation of phosphorus to the grain. SPDT is expressed in the xylem region of both enlarged- and diffuse-vascular bundles of the nodes, and encodes a plasma-membrane-localized transporter for phosphorus. Knockout of this gene in rice (Oryza sativa) altered the distribution of phosphorus, with decreased phosphorus in the grains but increased levels in the leaves. Total phosphorus and phytate in the brown de-husked rice were 20-30% lower in the knockout lines, whereas yield, seed germination and seedling vigour were not affected. These results indicate that SPDT functions in the rice node as a switch to allocate phosphorus preferentially to the grains. This finding provides a potential strategy to reduce the removal of phosphorus from the field and lower the risk of eutrophication of waterways.

Monocytes and macrophages comprise a variety of subsets with diverse functions. It is thought that these cells play a crucial role in homeostasis of peripheral organs, key immunological processes and development of various diseases. Among these diseases, fibrosis is a life-threatening disease of unknown aetiology. Its pathogenesis is poorly understood, and there are few effective therapies. The development of fibrosis is associated with activation of monocytes and macrophages. However, the specific subtypes of monocytes and macrophages that are involved in fibrosis have not yet been identified. Here we show that Ceacam1+Msr1+Ly6C-F4/80-Mac1+ monocytes, which we term segregated-nucleus-containing atypical monocytes (SatM), share granulocyte characteristics, are regulated by CCAAT/enhancer binding protein β (C/EBPβ), and are critical for fibrosis. Cebpb deficiency results in a complete lack of SatM. Furthermore, the development of bleomycin-induced fibrosis, but not inflammation, was prevented in chimaeric mice with Cebpb-/- haematopoietic cells. Adoptive transfer of SatM into Cebpb-/- mice resulted in fibrosis. Notably, SatM are derived from Ly6C-FcεRI+ granulocyte/macrophage progenitors, and a newly identified SatM progenitor downstream of Ly6C-FcεRI+ granulocyte/macrophage progenitors, but not from macrophage/dendritic-cell progenitors. Our results show that SatM are critical for fibrosis and that C/EBPβ licenses differentiation of SatM from their committed progenitor.

Approximately 1.5 billion people worldwide are overweight or affected by obesity, and are at risk of developing type 2 diabetes, cardiovascular disease and related metabolic and inflammatory disturbances. Although the mechanisms linking adiposity to associated clinical conditions are poorly understood, recent studies suggest that adiposity may influence DNA methylation, a key regulator of gene expression and molecular phenotype. Here we use epigenome-wide association to show that body mass index (BMI; a key measure of adiposity) is associated with widespread changes in DNA methylation (187 genetic loci with P < 1 × 10-7, range P = 9.2 × 10-8 to 6.0 × 10-46; n = 10,261 samples). Genetic association analyses demonstrate that the alterations in DNA methylation are predominantly the consequence of adiposity, rather than the cause. We find that methylation loci are enriched for functional genomic features in multiple tissues (P < 0.05), and show that sentinel methylation markers identify gene expression signatures at 38 loci (P < 9.0 × 10-6, range P = 5.5 × 10-6 to 6.1 × 10-35, n = 1,785 samples). The methylation loci identify genes involved in lipid and lipoprotein metabolism, substrate transport and inflammatory pathways. Finally, we show that the disturbances in DNA methylation predict future development of type 2 diabetes (relative risk per 1 standard deviation increase in methylation risk score: 2.3 (2.07-2.56); P = 1.1 × 10-54). Our results provide new insights into the biologic pathways influenced by adiposity, and may enable development of new strategies for prediction and prevention of type 2 diabetes and other adverse clinical consequences of obesity.

XRCC1 is a molecular scaffold protein that assembles multi-protein complexes involved in DNA single-strand break repair. Here we show that biallelic mutations in the human XRCC1 gene are associated with ocular motor apraxia, axonal neuropathy, and progressive cerebellar ataxia. Cells from a patient with mutations in XRCC1 exhibited not only reduced rates of single-strand break repair but also elevated levels of protein ADP-ribosylation. This latter phenotype is recapitulated in a related syndrome caused by mutations in the XRCC1 partner protein PNKP and implicates hyperactivation of poly(ADP-ribose) polymerase/s as a cause of cerebellar ataxia. Indeed, remarkably, genetic deletion of Parp1 rescued normal cerebellar ADP-ribose levels and reduced the loss of cerebellar neurons and ataxia in Xrcc1-defective mice, identifying a molecular mechanism by which endogenous single-strand breaks trigger neuropathology. Collectively, these data establish the importance of XRCC1 protein complexes for normal neurological function and identify PARP1 as a therapeutic target in DNA strand break repair-defective disease.

The heterotrimeric influenza polymerase (FluPol), comprising subunits PA, PB1 and PB2, binds to the conserved 5' and 3' termini (the 'promoter') of each of the eight single-stranded viral RNA (vRNA) genome segments and performs both transcription and replication of vRNA in the infected cell nucleus. To transcribe viral mRNAs, FluPol associates with cellular RNA polymerase II (Pol II), which enables it to take 5'-capped primers from nascent Pol II transcripts. Here we present a co-crystal structure of bat influenza A polymerase bound to a Pol II C-terminal domain (CTD) peptide mimic, which shows two distinct phosphoserine-5 (SeP5)-binding sites in the polymerase PA subunit, accommodating four CTD heptad repeats overall. Mutagenesis of the SeP5-contacting basic residues (PA K289, R454, K635 and R638) weakens CTD repeat binding in vitro without affecting the intrinsic cap-primed (transcription) or unprimed (replication) RNA synthesis activity of recombinant polymerase, whereas in cell-based minigenome assays the same mutations substantially reduce overall polymerase activity. Only recombinant viruses with a single mutation in one of the SeP5-binding sites can be rescued, but these viruses are severely attenuated and genetically unstable. Several previously described mutants that modulate virulence can be rationalized by our results, including a second site mutation (PA(C453R)) that enables the highly attenuated mutant virus (PA(R638A)) to revert to near wild-type infectivity. We conclude that direct binding of FluPol to the SeP5 Pol II CTD is fine-tuned to allow efficient viral transcription and propose that the CTD-binding site on FluPol could be targeted for antiviral drug development.

Internal bases in mRNA can be subjected to modifications that influence the fate of mRNA in cells. One of the most prevalent modified bases is found at the 5' end of mRNA, at the first encoded nucleotide adjacent to the 7-methylguanosine cap. Here we show that this nucleotide, N6,2'-O-dimethyladenosine (m6Am), is a reversible modification that influences cellular mRNA fate. Using a transcriptome-wide map of m6Am we find that m6Am-initiated transcripts are markedly more stable than mRNAs that begin with other nucleotides. We show that the enhanced stability of m6Am-initiated transcripts is due to resistance to the mRNA-decapping enzyme DCP2. Moreover, we find that m6Am is selectively demethylated by fat mass and obesity-associated protein (FTO). FTO preferentially demethylates m6Am rather than N6-methyladenosine (m6A), and reduces the stability of m6Am mRNAs. Together, these findings show that the methylation status of m6Am in the 5' cap is a dynamic and reversible epitranscriptomic modification that determines mRNA stability.

The macronutrient phosphorus is thought to limit primary productivity in the oceans on geological timescales. Although there has been a sustained effort to reconstruct the dynamics of the phosphorus cycle over the past 3.5 billion years, it remains uncertain whether phosphorus limitation persisted throughout Earth's history and therefore whether the phosphorus cycle has consistently modulated biospheric productivity and ocean-atmosphere oxygen levels over time. Here we present a compilation of phosphorus abundances in marine sedimentary rocks spanning the past 3.5 billion years. We find evidence for relatively low authigenic phosphorus burial in shallow marine environments until about 800 to 700 million years ago. Our interpretation of the database leads us to propose that limited marginal phosphorus burial before that time was linked to phosphorus biolimitation, resulting in elemental stoichiometries in primary producers that diverged strongly from the Redfield ratio (the atomic ratio of carbon, nitrogen and phosphorus found in phytoplankton). We place our phosphorus record in a quantitative biogeochemical model framework and find that a combination of enhanced phosphorus scavenging in anoxic, iron-rich oceans and a nutrient-based bistability in atmospheric oxygen levels could have resulted in a stable low-oxygen world. The combination of these factors may explain the protracted oxygenation of Earth's surface over the last 3.5 billion years of Earth history. However, our analysis also suggests that a fundamental shift in the phosphorus cycle may have occurred during the late Proterozoic eon (between 800 and 635 million years ago), coincident with a previously inferred shift in marine redox states, severe perturbations to Earth's climate system, and the emergence of animals.

Rational control of the self-assembly of large structures is one of the key challenges in chemistry, and is believed to become increasingly difficult and ultimately impossible as the number of components involved increases. So far, it has not been possible to design a self-assembled discrete molecule made up of more than 100 components. Such molecules-for example, spherical virus capsids-are prevalent in nature, which suggests that the difficulty in designing these very large self-assembled molecules is due to a lack of understanding of the underlying design principles. For example, the targeted assembly of a series of large spherical structures containing up to 30 palladium ions coordinated by up to 60 bent organic ligands was achieved by considering their topologies. Here we report the self-assembly of a spherical structure that also contains 30 palladium ions and 60 bent ligands, but belongs to a shape family that has not previously been observed experimentally. The new structure consists of a combination of 8 triangles and 24 squares, and has the symmetry of a tetravalent Goldberg polyhedron. Platonic and Archimedean solids have previously been prepared through self-assembly, as have trivalent Goldberg polyhedra, which occur naturally in the form of virus capsids and fullerenes. But tetravalent Goldberg polyhedra have not previously been reported at the molecular level, although their topologies have been predicted using graph theory. We use graph theory to predict the self-assembly of even larger tetravalent Goldberg polyhedra, which should be more stable, enabling another member of this polyhedron family to be assembled from 144 components: 48 palladium ions and 96 bent ligands.

Aromatic and antiaromatic molecules-which have delocalized circuits of [4n + 2] or [4n] electrons, respectively-exhibit ring currents around their perimeters. The direction of the ring current in an aromatic molecule is such as to generate a magnetic field that opposes the external field inside the ring (a 'diatropic' current), while the ring current in an antiaromatic molecule flows in the reverse direction ('paratropic'). Similar persistent currents occur in metal or semiconductor rings, when the phase coherence of the electronic wavefunction is preserved around the ring. Persistent currents in non-molecular rings switch direction as a function of the magnetic flux passing through the ring, so that they can be changed from diatropic ('aromatic') to paratropic ('antiaromatic') simply by changing the external magnetic field. As in molecular systems, the direction of the persistent current also depends on the number of electrons. The relationship between ring currents in molecular and non-molecular rings is poorly understood, partly because they are studied in different size regimes: the largest aromatic molecules have diameters of about one nanometre, whereas persistent currents are observed in microfabricated rings with diameters of 20-1,000 nanometres. Understanding the connection between aromaticity and quantum-coherence effects in mesoscopic rings provides a motivation for investigating ring currents in molecules of an intermediate size. Here we show, using nuclear magnetic resonance spectroscopy and density functional theory, that a six-porphyrin nanoring template complex, with a diameter of 2.4 nanometres, is antiaromatic in its 4+ oxidation state (80 π electrons) and aromatic in its 6+ oxidation state (78 π electrons). The antiaromatic state has a huge paramagnetic susceptibility, despite having no unpaired electrons. This work demonstrates that a global ring current can be promoted in a macrocycle by adjusting its oxidation state to suppress the local ring currents of its components.The discovery of ring currents around a molecule with a circumference of 7.5 nanometres, at room temperature, shows that quantum coherence can persist in surprisingly large molecular frameworks.

Packaging of the genome into a protein capsid and its subsequent delivery into a host cell are two fundamental processes in the life cycle of a virus. Unlike double-stranded DNA viruses, which pump their genome into a preformed capsid, single-stranded RNA (ssRNA) viruses, such as bacteriophage MS2, co-assemble their capsid with the genome; however, the structural basis of this co-assembly is poorly understood. MS2 infects Escherichia coli via the host 'sex pilus' (F-pilus); it was the first fully sequenced organism and is a model system for studies of translational gene regulation, RNA-protein interactions, and RNA virus assembly. Its positive-sense ssRNA genome of 3,569 bases is enclosed in a capsid with one maturation protein monomer and 89 coat protein dimers arranged in a T = 3 icosahedral lattice. The maturation protein is responsible for attaching the virus to an F-pilus and delivering the viral genome into the host during infection, but how the genome is organized and delivered is not known. Here we describe the MS2 structure at 3.6 Å resolution, determined by electron-counting cryo-electron microscopy (cryoEM) and asymmetric reconstruction. We traced approximately 80% of the backbone of the viral genome, built atomic models for 16 RNA stem-loops, and identified three conserved motifs of RNA-coat protein interactions among 15 of these stem-loops with diverse sequences. The stem-loop at the 3' end of the genome interacts extensively with the maturation protein, which, with just a six-helix bundle and a six-stranded β-sheet, forms a genome-delivery apparatus and joins 89 coat protein dimers to form a capsid. This atomic description of genome-capsid interactions in a spherical ssRNA virus provides insight into genome delivery via the host sex pilus and mechanisms underlying ssRNA-capsid co-assembly, and inspires speculation about the links between nucleoprotein complexes and the origins of viruses.