An improved Nicotiana benthamiana bioproduction chassis provides novel insights into nicotine biosynthesis.

The model plant Nicotiana benthamiana is an increasingly attractive organism for the production of high-value, biologically active molecules. However, N. benthamiana accumulates high levels of pyridine alkaloids, in particular nicotine, which complicates the downstream purification processes. Here, we report a new assembly of the N. benthamiana genome as well as the generation of low-nicotine lines by CRISPR/Cas9-based inactivation of berberine bridge enzyme-like proteins (BBLs). Triple as well as quintuple mutants accumulated three to four times less nicotine than the respective control lines. The availability of lines without functional BBLs allowed us to probe their catalytic role in nicotine biosynthesis, which has remained obscure. Notably, chiral analysis revealed that the enantiomeric purity of nicotine was fully lost in the quintuple mutants. In addition, precursor feeding experiments showed that these mutants cannot facilitate the specific loss of C6 hydrogen that characterizes natural nicotine biosynthesis. Our work delivers an improved N. benthamiana chassis for bioproduction and uncovers the crucial role of BBLs in the stereoselectivity of nicotine biosynthesis.

Genomic erosion in a demographically recovered bird species during conservation rescue.

The pink pigeon (Nesoenas mayeri) is an endemic species of Mauritius that has made a remarkable recovery after a severe population bottleneck in the 1970s to early 1990s. Prior to this bottleneck, an ex situ population was established from which captive-bred individuals were released into free-living subpopulations to increase population size and genetic variation. This conservation rescue led to rapid population recovery to 400-480 individuals, and the species was twice downlisted on the International Union for the Conservation of Nature (IUCN) Red List. We analyzed the impacts of the bottleneck and genetic rescue on neutral genetic variation during and after population recovery (1993-2008) with restriction site-associated sequencing, microsatellite analyses, and quantitative genetic analysis of studbook data of 1112 birds from zoos in Europe and the United States. We used computer simulations to study the predicted changes in genetic variation and population viability from the past into the future. Genetic variation declined rapidly, despite the population rebound, and the effective population size was approximately an order of magnitude smaller than census size. The species carried a high genetic load of circa 15 lethal equivalents for longevity. Our computer simulations predicted continued inbreeding will likely result in increased expression of deleterious mutations (i.e., a high realized load) and severe inbreeding depression. Without continued conservation actions, it is likely that the pink pigeon will go extinct in the wild within 100 years. Conservation rescue of the pink pigeon has been instrumental in the recovery of the free-living population. However, further genetic rescue with captive-bred birds from zoos is required to recover lost variation, reduce expression of harmful deleterious variation, and prevent extinction. The use of genomics and modeling data can inform IUCN assessments of the viability and extinction risk of species, and it helps in assessments of the conservation dependency of populations.

La paloma rosada (Nesoenas mayeri) es una especie endémica de Mauricio que se ha recuperado impresionantemente después de un grave cuello de botella poblacional a principios de la década de 1970 que duró hasta inicios de la década de 1990. Antes de este cuello de botella se había establecido una población ex situ de la cual se liberaban individuos reproducidos en cautiverio a las subpoblaciones en libertad para incrementar la variación genética y el tamaño poblacional. Este rescate de conservación derivó en una recuperación rápida de la población (400-480 individuos) y la especie cambió positivamente de categoría dos veces en la Lista Roja de la Unión Internacional para la Conservación de la Naturaleza (UICN). Analizamos los impactos del cuello de botella y el rescate genético sobre la variación genética neutral durante y después de la recuperación poblacional (de 1993 a 2008) mediante secuenciación RAD, análisis de microsatélites y análisis genéticos cuantitativos de los datos del libro genealógico de 1112 aves ubicadas en zoológicos de Europa y los Estados Unidos. Usamos simulaciones por computadora para estudiar los cambios pronosticados en la variación genética y en la viabilidad poblacional del pasado hacia el futuro. La variación genética declinó rápidamente, a pesar de la recuperación poblacional, y el tamaño efectivo de la población fue aproximadamente un orden de magnitud más pequeño que el tamaño del censo. La especie contó con una carga genética elevada de casi 15 equivalentes letales para la longevidad. Nuestras simulaciones pronostican que la endogamia continua probablemente resultará en un incremento en la expresión de mutaciones deletéreas (es decir, una carga realizada elevada) y en una depresión endogámica severa. Sin acciones continuas para la conservación, es probable que la paloma rosada esté extinta en vida libre dentro de cien años. El rescate de conservación de la paloma rosada ha sido fundamental en la recuperación de la población silvestre; sin embargo, se requiere de un rescate genético adicional con las aves de reproducción en cautiverio de los zoológicos para recuperar la variación perdida, reducir la expresión de la variación deletérea dañina y prevenir la extinción. El uso de la genómica y los datos modelados puede orientar las valoraciones de la UICN sobre la viabilidad y el riesgo de extinción de las especies, además de que ayuda en la evaluación de la dependencia que tienen las poblaciones de la conservación.

Nanopore adaptive sampling: a tool for enrichment of low abundance species in metagenomic samples.

Adaptive sampling is a method of software-controlled enrichment unique to nanopore sequencing platforms. To test its potential for enrichment of rarer species within metagenomic samples, we create a synthetic mock community and construct sequencing libraries with a range of mean read lengths. Enrichment is up to 13.87-fold for the least abundant species in the longest read length library; factoring in reduced yields from rejecting molecules the calculated efficiency raises this to 4.93-fold. Finally, we introduce a mathematical model of enrichment based on molecule length and relative abundance, whose predictions correlate strongly with mock and complex real-world microbial communities.

An accessible, efficient and global approach for the large-scale sequencing of bacterial genomes.

We have developed an efficient and inexpensive pipeline for streamlining large-scale collection and genome sequencing of bacterial isolates. Evaluation of this method involved a worldwide research collaboration focused on the model organism Salmonella enterica, the 10KSG consortium. Following the optimization of a logistics pipeline that involved shipping isolates as thermolysates in ambient conditions, the project assembled a diverse collection of 10,419 isolates from low- and middle-income countries. The genomes were sequenced using the LITE pipeline for library construction, with a total reagent cost of less than USD$10 per genome. Our method can be applied to other large bacterial collections to underpin global collaborations.

Characteristics of Salmonella Recovered From Stools of Children Enrolled in the Global Enteric Multicenter Study.

BACKGROUND: The Global Enteric Multicenter Study (GEMS) determined the etiologic agents of moderate-to-severe diarrhea (MSD) in children under 5 years old in Africa and Asia. Here, we describe the prevalence and antimicrobial susceptibility of nontyphoidal Salmonella (NTS) serovars in GEMS and examine the phylogenetics of Salmonella Typhimurium ST313 isolates. METHODS: Salmonella isolated from children with MSD or diarrhea-free controls were identified by classical clinical microbiology and serotyped using antisera and/or whole-genome sequence data. We evaluated antimicrobial susceptibility using the Kirby-Bauer disk-diffusion method. Salmonella Typhimurium sequence types were determined using multi-locus sequence typing, and whole-genome sequencing was performed to assess the phylogeny of ST313. RESULTS: Of 370 Salmonella-positive individuals, 190 (51.4%) were MSD cases and 180 (48.6%) were diarrhea-free controls. The most frequent Salmonella serovars identified were Salmonella Typhimurium, serogroup O:8 (C2-C3), serogroup O:6,7 (C1), Salmonella Paratyphi B Java, and serogroup O:4 (B). The prevalence of NTS was low but similar across sites, regardless of age, and was similar among both cases and controls except in Kenya, where Salmonella Typhimurium was more commonly associated with cases than controls. Phylogenetic analysis showed that these Salmonella Typhimurium isolates, all ST313, were highly genetically related to isolates from controls. Generally, Salmonella isolates from Asia were resistant to ciprofloxacin and ceftriaxone, but African isolates were susceptible to these antibiotics. CONCLUSIONS: Our data confirm that NTS is prevalent, albeit at low levels, in Africa and South Asia. Our findings provide further evidence that multidrug-resistant Salmonella Typhimurium ST313 can be carried asymptomatically by humans in sub-Saharan Africa.

Multiple wheat genomes reveal global variation in modern breeding.

Advances in genomics have expedited the improvement of several agriculturally important crops but similar efforts in wheat (Triticum spp.) have been more challenging. This is largely owing to the size and complexity of the wheat genome1, and the lack of genome-assembly data for multiple wheat lines2,3. Here we generated ten chromosome pseudomolecule and five scaffold assemblies of hexaploid wheat to explore the genomic diversity among wheat lines from global breeding programs. Comparative analysis revealed extensive structural rearrangements, introgressions from wild relatives and differences in gene content resulting from complex breeding histories aimed at improving adaptation to diverse environments, grain yield and quality, and resistance to stresses4,5. We provide examples outlining the utility of these genomes, including a detailed multi-genome-derived nucleotide-binding leucine-rich repeat protein repertoire involved in disease resistance and the characterization of Sm16, a gene associated with insect resistance. These genome assemblies will provide a basis for functional gene discovery and breeding to deliver the next generation of modern wheat cultivars.

A low-cost pipeline for soil microbiome profiling.

Common bottlenecks in environmental and crop microbiome studies are the consumable and personnel costs necessary for genomic DNA extraction and sequencing library construction. This is harder for challenging environmental samples such as soil, which is rich in Polymerase Chain Reaction (PCR) inhibitors. To address this, we have established a low-cost genomic DNA extraction method for soil samples. We also present an Illumina-compatible 16S and ITS rRNA gene amplicon library preparation workflow that uses common laboratory equipment. We evaluated the performance of our genomic DNA extraction method against two leading commercial soil genomic DNA kits (MoBio PowerSoil® and MP Biomedicals™ FastDNA™ SPIN) and a recently published non-commercial extraction method by Zou et al. (PLoS Biology, 15, e2003916, 2017). Our benchmarking experiment used four different soil types (coniferous, broad-leafed, and mixed forest plus a standardized cereal crop compost mix) assessing the quality and quantity of the extracted genomic DNA by analyzing sequence variants of 16S V4 and ITS rRNA amplicons. We found that our genomic DNA extraction method compares well to both commercially available genomic DNA extraction kits in DNA quality and quantity. The MoBio PowerSoil® kit, which relies on silica column-based DNA extraction with extensive washing, delivered the cleanest genomic DNA, for example, best A260:A280 and A260:A230 absorbance ratios. The MP Biomedicals™ FastDNA™ SPIN kit, which uses a large amount of binding material, yielded the most genomic DNA. Our method fits between the two commercial kits, producing both good yields and clean genomic DNA with fragment sizes of approximately 10 kb. Comparative analysis of detected amplicon sequence variants shows that our method correlates well with the two commercial kits. Here, we present a low-cost genomic DNA extraction method for soil samples that can be coupled to an Illumina-compatible simple two-step amplicon library construction workflow for 16S V4 and ITS marker genes. Our method delivers high-quality genomic DNA at a fraction of the cost of commercial kits and enables cost-effective, large-scale amplicon sequencing projects. Notably, our extracted gDNA molecules are long enough to be suitable for downstream techniques such as full gene sequencing or even metagenomics shotgun approaches using long reads (PacBio or Nanopore), 10x Genomics linked reads, and Dovetail genomics.

A novel allele of ASY3 is associated with greater meiotic stability in autotetraploid Arabidopsis lyrata.

In this study we performed a genotype-phenotype association analysis of meiotic stability in 10 autotetraploid Arabidopsis lyrata and A. lyrata/A. arenosa hybrid populations collected from the Wachau region and East Austrian Forealps. The aim was to determine the effect of eight meiosis genes under extreme selection upon adaptation to whole genome duplication. Individual plants were genotyped by high-throughput sequencing of the eight meiosis genes (ASY1, ASY3, PDS5b, PRD3, REC8, SMC3, ZYP1a/b) implicated in synaptonemal complex formation and phenotyped by assessing meiotic metaphase I chromosome configurations. Our results reveal that meiotic stability varied greatly (20-100%) between individual tetraploid plants and associated with segregation of a novel ASYNAPSIS3 (ASY3) allele derived from A. lyrata. The ASY3 allele that associates with meiotic stability possesses a putative in-frame tandem duplication (TD) of a serine-rich region upstream of the coiled-coil domain that appears to have arisen at sites of DNA microhomology. The frequency of multivalents observed in plants homozygous for the ASY3 TD haplotype was significantly lower than in plants heterozygous for ASY3 TD/ND (non-duplicated) haplotypes. The chiasma distribution was significantly altered in the stable plants compared to the unstable plants with a shift from proximal and interstitial to predominantly distal locations. The number of HEI10 foci at pachytene that mark class I crossovers was significantly reduced in a plant homozygous for ASY3 TD compared to a plant heterozygous for ASY3 ND/TD. Fifty-eight alleles of the 8 meiosis genes were identified from the 10 populations analysed, demonstrating dynamic population variability at these loci. Widespread chimerism between alleles originating from A. lyrata/A. arenosa and diploid/tetraploids indicates that this group of rapidly evolving genes may provide precise adaptive control over meiotic recombination in the tetraploids, the very process that gave rise to them.

Metagenomic analysis of planktonic riverine microbial consortia using nanopore sequencing reveals insight into river microbe taxonomy and function.

BACKGROUND: Riverine ecosystems are biogeochemical powerhouses driven largely by microbial communities that inhabit water columns and sediments. Because rivers are used extensively for anthropogenic purposes (drinking water, recreation, agriculture, and industry), it is essential to understand how these activities affect the composition of river microbial consortia. Recent studies have shown that river metagenomes vary considerably, suggesting that microbial community data should be included in broad-scale river ecosystem models. But such ecogenomic studies have not been applied on a broad "aquascape" scale, and few if any have applied the newest nanopore technology. RESULTS: We investigated the metagenomes of 11 rivers across 3 continents using MinION nanopore sequencing, a portable platform that could be useful for future global river monitoring. Up to 10 Gb of data per run were generated with average read lengths of 3.4 kb. Diversity and diagnosis of river function potential was accomplished with 0.5-1.0 ⋅ 106 long reads. Our observations for 7 of the 11 rivers conformed to other river-omic findings, and we exposed previously unrecognized microbial biodiversity in the other 4 rivers. CONCLUSIONS: Deeper understanding that emerged is that river microbial consortia and the ecological functions they fulfil did not align with geographic location but instead implicated ecological responses of microbes to urban and other anthropogenic effects, and that changes in taxa manifested over a very short geographic space.

Sequencing smart: De novo sequencing and assembly approaches for a non-model mammal.

BACKGROUND: Whilst much sequencing effort has focused on key mammalian model organisms such as mouse and human, little is known about the relationship between genome sequencing techniques for non-model mammals and genome assembly quality. This is especially relevant to non-model mammals, where the samples to be sequenced are often degraded and of low quality. A key aspect when planning a genome project is the choice of sequencing data to generate. This decision is driven by several factors, including the biological questions being asked, the quality of DNA available, and the availability of funds. Cutting-edge sequencing technologies now make it possible to achieve highly contiguous, chromosome-level genome assemblies, but rely on high-quality high molecular weight DNA. However, funding is often insufficient for many independent research groups to use these techniques. Here we use a range of different genomic technologies generated from a roadkill European polecat (Mustela putorius) to assess various assembly techniques on this low-quality sample. We evaluated different approaches for de novo assemblies and discuss their value in relation to biological analyses. RESULTS: Generally, assemblies containing more data types achieved better scores in our ranking system. However, when accounting for misassemblies, this was not always the case for Bionano and low-coverage 10x Genomics (for scaffolding only). We also find that the extra cost associated with combining multiple data types is not necessarily associated with better genome assemblies. CONCLUSIONS: The high degree of variability between each de novo assembly method (assessed from the 7 key metrics) highlights the importance of carefully devising the sequencing strategy to be able to carry out the desired analysis. Adding more data to genome assemblies does not always result in better assemblies, so it is important to understand the nuances of genomic data integration explained here, in order to obtain cost-effective value for money when sequencing genomes.

A Genome Assembly of the Barley 'Transformation Reference' Cultivar Golden Promise.

Barley (Hordeum vulgare) is one of the most important crops worldwide and is also considered a research model for the large-genome small grain temperate cereals. Despite genomic resources improving all the time, they are limited for the cv Golden Promise, the most efficient genotype for genetic transformation. We have developed a barley cv Golden Promise reference assembly integrating Illumina paired-end reads, long mate-pair reads, Dovetail Chicago in vitro proximity ligation libraries and chromosome conformation capture sequencing (Hi-C) libraries into a contiguous reference assembly. The assembled genome of 7 chromosomes and 4.13Gb in size, has a super-scaffold N50 after Chicago libraries of 4.14Mb and contains only 2.2% gaps. Using BUSCO (benchmarking universal single copy orthologous genes) as evaluation the genome assembly contains 95.2% of complete and single copy genes from the plant database. A high-quality Golden Promise reference assembly will be useful and utilized by the whole barley research community but will prove particularly useful for CRISPR-Cas9 experiments.

Rapid MinION profiling of preterm microbiota and antimicrobial-resistant pathogens.

The MinION sequencing platform offers near real-time analysis of DNA sequence; this makes the tool attractive for deployment in fieldwork or clinical settings. We used the MinION platform coupled to the NanoOK RT software package to perform shotgun metagenomic sequencing and profile mock communities and faecal samples from healthy and ill preterm infants. Using Nanopore data, we reliably classified a 20-species mock community and captured the diversity of the immature gut microbiota over time and in response to interventions such as probiotic supplementation, antibiotic treatment or episodes of suspected sepsis. We also performed rapid real-time runs to assess gut-associated microbial communities in critically ill and healthy infants, facilitated by NanoOK RT software package, which analysed sequences as they were generated. Our pipeline reliably identified pathogenic bacteria (that is, Klebsiella pneumoniae and Enterobacter cloacae) and their corresponding antimicrobial resistance gene profiles within as little as 1 h of sequencing. Results were confirmed using pathogen isolation, whole-genome sequencing and antibiotic susceptibility testing, as well as mock communities and clinical samples with known antimicrobial resistance genes. Our results demonstrate that MinION (including cost-effective Flongle flow cells) with NanoOK RT can process metagenomic samples to a rich dataset in < 5 h, which creates a platform for future studies aimed at developing these tools and approaches in clinical settings with a focus on providing tailored patient antimicrobial treatment options.

Spatially resolved transcriptomics reveals plant host responses to pathogens.

BACKGROUND: Thorough understanding of complex model systems requires the characterisation of processes in different cell types of an organism. This can be achieved with high-throughput spatial transcriptomics at a large scale. However, for plant model systems this is still challenging as suitable transcriptomics methods are sparsely available. Here we present GaST-seq (Grid-assisted, Spatial Transcriptome sequencing), an easy to adopt, micro-scale spatial-transcriptomics workflow that allows to study expression profiles across small areas of plant tissue at a fraction of the cost of existing sequencing-based methods. RESULTS: We compare the GaST-seq method with widely used library preparation methods (Illumina TruSeq). In spatial experiments we show that the GaST-seq method is sensitive enough to identify expression differences across a plant organ. We further assess the spatial transcriptome response of Arabidopsis thaliana leaves exposed to the bacterial molecule flagellin-22, and show that with eukaryotic (Albugo laibachii) infection both host and pathogen spatial transcriptomes are obtained. CONCLUSION: We show that our method can be used to identify known, rapidly flagellin-22 elicited genes, plant immune response pathways to bacterial attack and spatial expression patterns of genes associated with these pathways.

An Ultra High-Density Arabidopsis thaliana Crossover Map That Refines the Influences of Structural Variation and Epigenetic Features.

Many environmental, genetic, and epigenetic factors are known to affect the frequency and positioning of meiotic crossovers (COs). Suppression of COs by large, cytologically visible inversions and translocations has long been recognized, but relatively little is known about how smaller structural variants (SVs) affect COs. To examine fine-scale determinants of the CO landscape, including SVs, we used a rapid, cost-effective method for high-throughput sequencing to generate a precise map of >17,000 COs between the Col-0 and Ler-0 accessions of Arabidopsis thaliana COs were generally suppressed in regions with SVs, but this effect did not depend on the size of the variant region, and was only marginally affected by the variant type. CO suppression did not extend far beyond the SV borders and CO rates were slightly elevated in the flanking regions. Disease resistance gene clusters, which often exist as SVs, exhibited high CO rates at some loci, but there was a tendency toward depressed CO rates at loci where large structural differences exist between the two parents. Our high-density map also revealed in fine detail how CO positioning relates to genetic (DNA motifs) and epigenetic (chromatin structure) features of the genome. We conclude that suppression of COs occurs over a narrow region spanning large- and small-scale SVs, representing an influence on the CO landscape in addition to sequence and epigenetic variation along chromosomes.

A critical comparison of technologies for a plant genome sequencing project.

BACKGROUND: A high-quality genome sequence of any model organism is an essential starting point for genetic and other studies. Older clone-based methods are slow and expensive, whereas faster, cheaper short-read-only assemblies can be incomplete and highly fragmented, which minimizes their usefulness. The last few years have seen the introduction of many new technologies for genome assembly. These new technologies and associated new algorithms are typically benchmarked on microbial genomes or, if they scale appropriately, on larger (e.g., human) genomes. However, plant genomes can be much more repetitive and larger than the human genome, and plant biochemistry often makes obtaining high-quality DNA that is free from contaminants difficult. Reflecting their challenging nature, we observe that plant genome assembly statistics are typically poorer than for vertebrates. RESULTS: Here, we compare Illumina short read, Pacific Biosciences long read, 10x Genomics linked reads, Dovetail Hi-C, and BioNano Genomics optical maps, singly and combined, in producing high-quality long-range genome assemblies of the potato species Solanum verrucosum. We benchmark the assemblies for completeness and accuracy, as well as DNA compute requirements and sequencing costs. CONCLUSIONS: The field of genome sequencing and assembly is reaching maturity, and the differences we observe between assemblies are surprisingly small. We expect that our results will be helpful to other genome projects, and that these datasets will be used in benchmarking by assembly algorithm developers.

Arabidopsis thaliana and Pseudomonas Pathogens Exhibit Stable Associations over Evolutionary Timescales.

Crop disease outbreaks are often associated with clonal expansions of single pathogenic lineages. To determine whether similar boom-and-bust scenarios hold for wild pathosystems, we carried out a multi-year, multi-site survey of Pseudomonas in its natural host Arabidopsis thaliana. The most common Pseudomonas lineage corresponded to a ubiquitous pathogenic clade. Sequencing of 1,524 genomes revealed this lineage to have diversified approximately 300,000 years ago, containing dozens of genetically identifiable pathogenic sublineages. There is differentiation at the level of both gene content and disease phenotype, although the differentiation may not provide fitness advantages to specific sublineages. The coexistence of sublineages indicates that in contrast to crop systems, no single strain has been able to overtake the studied A. thaliana populations in the recent past. Our results suggest that selective pressures acting on a plant pathogen in wild hosts are likely to be much more complex than those in agricultural systems.

Independent assessment and improvement of wheat genome sequence assemblies using Fosill jumping libraries.

Background: The accurate sequencing and assembly of very large, often polyploid, genomes remains a challenging task, limiting long-range sequence information and phased sequence variation for applications such as plant breeding. The 15-Gb hexaploid bread wheat (Triticum aestivum) genome has been particularly challenging to sequence, and several different approaches have recently generated long-range assemblies. Mapping and understanding the types of assembly errors are important for optimising future sequencing and assembly approaches and for comparative genomics. Results: Here we use a Fosill 38-kb jumping library to assess medium and longer-range order of different publicly available wheat genome assemblies. Modifications to the Fosill protocol generated longer Illumina sequences and enabled comprehensive genome coverage. Analyses of two independent Bacterial Artificial Chromosome (BAC)-based chromosome-scale assemblies, two independent Illumina whole genome shotgun assemblies, and a hybrid Single Molecule Real Time (SMRT-PacBio) and short read (Illumina) assembly were carried out. We revealed a surprising scale and variety of discrepancies using Fosill mate-pair mapping and validated several of each class. In addition, Fosill mate-pairs were used to scaffold a whole genome Illumina assembly, leading to a 3-fold increase in N50 values. Conclusions: Our analyses, using an independent means to validate different wheat genome assemblies, show that whole genome shotgun assemblies based solely on Illumina sequences are significantly more accurate by all measures compared to BAC-based chromosome-scale assemblies and hybrid SMRT-Illumina approaches. Although current whole genome assemblies are reasonably accurate and useful, additional improvements will be needed to generate complete assemblies of wheat genomes using open-source, computationally efficient, and cost-effective methods.

An improved assembly and annotation of the allohexaploid wheat genome identifies complete families of agronomic genes and provides genomic evidence for chromosomal translocations.

Advances in genome sequencing and assembly technologies are generating many high-quality genome sequences, but assemblies of large, repeat-rich polyploid genomes, such as that of bread wheat, remain fragmented and incomplete. We have generated a new wheat whole-genome shotgun sequence assembly using a combination of optimized data types and an assembly algorithm designed to deal with large and complex genomes. The new assembly represents >78% of the genome with a scaffold N50 of 88.8 kb that has a high fidelity to the input data. Our new annotation combines strand-specific Illumina RNA-seq and Pacific Biosciences (PacBio) full-length cDNAs to identify 104,091 high-confidence protein-coding genes and 10,156 noncoding RNA genes. We confirmed three known and identified one novel genome rearrangements. Our approach enables the rapid and scalable assembly of wheat genomes, the identification of structural variants, and the definition of complete gene models, all powerful resources for trait analysis and breeding of this key global crop.

A chromosome conformation capture ordered sequence of the barley genome.

Cereal grasses of the Triticeae tribe have been the major food source in temperate regions since the dawn of agriculture. Their large genomes are characterized by a high content of repetitive elements and large pericentromeric regions that are virtually devoid of meiotic recombination. Here we present a high-quality reference genome assembly for barley (Hordeum vulgare L.). We use chromosome conformation capture mapping to derive the linear order of sequences across the pericentromeric space and to investigate the spatial organization of chromatin in the nucleus at megabase resolution. The composition of genes and repetitive elements differs between distal and proximal regions. Gene family analyses reveal lineage-specific duplications of genes involved in the transport of nutrients to developing seeds and the mobilization of carbohydrates in grains. We demonstrate the importance of the barley reference sequence for breeding by inspecting the genomic partitioning of sequence variation in modern elite germplasm, highlighting regions vulnerable to genetic erosion.