Multiplex sequencing of bacterial artificial chromosomes for assembling complex plant genomes

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

Natural variations of HvSRN1 modulate the spike rachis node number in barley

Grain number, one of the major determinants of yield in Triticeae crops, is largely determined by spikelet number and spike rachis node number (SRN). Here, we identified three quantitative trait loci (QTLs) for SRN using 145 recombinant inbred lines derived from a barley R90/1815D cross. qSRN1, the major-effect QTL, was mapped to chromosome 2H and explained up to 38.77% of SRN variation. Map-based cloning revealed that qSRN1 encodes the RAWUL domain-containing protein HvSRN1. Further analysis revealed that two key SNPs in the HvSRN1 promoter region (∼2 kb upstream of the transcription start site) affect the transcript level of HvSRN1 and contribute to variation in SRN. Similar to its orthologous proteins OsLAX2 and ZmBA2, HvSRN1 showed protein-protein interactions with HvLAX1, suggesting that the LAX2-LAX1 model for spike morphology regulation may be conserved in Poaceae crops. CRISPR-Cas9-induced HvSRN1 mutants showed reduced SRN but increased grain size and weight, demonstrating a trade-off effect. Our results shed light on the role of HvSRN1 variation in regulating the balance between grain number and weight in barley.

Genome Assembly of the Canadian two-row Malting Barley cultivar AAC Synergy

Barley (Hordeum vulgare L.) is one of the most important global crops. The six-row barley cultivar Morex reference genome has been used by the barley research community worldwide. However, this reference genome can have limitations when used for genomic and genetic diversity analysis studies, gene discovery, and marker development when working in two-row germplasm that is more common to Canadian barley. Here we assembled, for the first time, the genome sequence of a Canadian two-row malting barley, cultivar AAC Synergy. We applied deep Illumina paired-end reads, long mate-pair reads, PacBio sequences, 10X chromium linked read libraries, and chromosome conformation capture sequencing (Hi-C) to generate a contiguous assembly. The genome assembled from super-scaffolds had a size of 4.85 Gb, N50 of 2.32 Mb, and an estimated 93.9% of complete genes from a plant database (BUSCO, benchmarking universal single-copy orthologous genes). After removal of small scaffolds (< 300 Kb), the assembly was arranged into pseudomolecules of 4.14 Gb in size with seven chromosomes plus unanchored scaffolds. The completeness and annotation of the assembly were assessed by comparing it with the updated version of six-row Morex and recently released two-row Golden Promise genome assemblies.

TRITEX: chromosome-scale sequence assembly of Triticeae genomes with open-source tools

Chromosome-scale genome sequence assemblies underpin pan-genomic studies. Recent genome assembly efforts in the large-genome Triticeae crops wheat and barley have relied on the commercial closed-source assembly algorithm DeNovoMagic. We present TRITEX, an open-source computational workflow that combines paired-end, mate-pair, 10X Genomics linked-read with chromosome conformation capture sequencing data to construct sequence scaffolds with megabase-scale contiguity ordered into chromosomal pseudomolecules. We evaluate the performance of TRITEX on publicly available sequence data of tetraploid wild emmer and hexaploid bread wheat, and construct an improved annotated reference genome sequence assembly of the barley cultivar Morex as a community resource.

TRITEX: chromosome-scale sequence assembly of Triticeae genomes with open-source tools

Chromosome-scale genome sequence assemblies underpin pan-genomic studies. Recent genome assembly efforts in the large-genome Triticeae crops wheat and barley have relied on the commercial closed-source assembly algorithm DeNovoMagic. We present TRITEX, an open-source computational workflow that combines paired-end, mate-pair, 10X Genomics linked-read with chromosome conformation capture sequencing data to construct sequence scaffolds with megabase-scale contiguity ordered into chromosomal pseudomolecules. We evaluate the performance of TRITEX on publicly available sequence data of tetraploid wild emmer and hexaploid bread wheat, and construct an improved annotated reference genome sequence assembly of the barley cultivar Morex as a community resource.

Leaf Variegation and Impaired Chloroplast Development Caused by a Truncated CCT Domain Gene in albostrians Barley

Chloroplasts fuel plant development and growth by converting solar energy into chemical energy. They mature from proplastids through the concerted action of genes in both the organellar and the nuclear genome. Defects in such genes impair chloroplast development and may lead to pigment-deficient seedlings or seedlings with variegated leaves. Such mutants are instrumental as tools for dissecting genetic factors underlying the mechanisms involved in chloroplast biogenesis. Characterization of the green-white variegated albostrians mutant of barley (Hordeum vulgare) has greatly broadened the field of chloroplast biology, including the discovery of retrograde signaling. Here, we report identification of the ALBOSTRIANS gene HvAST (also known as Hordeum vulgare CCT Motif Family gene 7, HvCMF7) by positional cloning as well as its functional validation based on independently induced mutants by Targeting Induced Local Lesions in Genomes (TILLING) and RNA-guided clustered regularly interspaced short palindromic repeats-associated protein 9 endonuclease-mediated gene editing. The phenotypes of the independent HvAST mutants imply residual activity of HvCMF7 in the original albostrians allele conferring an imperfect penetrance of the variegated phenotype even at homozygous state of the mutation. HvCMF7 is a homolog of the Arabidopsis (Arabidopsis thaliana) CONSTANS, CO-like, and TOC1 (CCT) Motif transcription factor gene CHLOROPLAST IMPORT APPARATUS2, which was reported to be involved in the expression of nuclear genes essential for chloroplast biogenesis. Notably, in barley we localized HvCMF7 to the chloroplast, without any clear evidence for nuclear localization.

Prospects of pan-genomics in barley

The concept of a pan-genome refers to intraspecific diversity in genome content and structure, encompassing both genes and intergenic space. Pan-genomic studies employ a combination of de novo sequence assembly and reference-based alignment to discover and genotype structural variants. The large size and complex structure of Triticeae genomes were for a long time an obstacle for genomic research in barley and its relatives. Now that a reference genome is available, computational pipelines for high-quality sequence assembly are in place, and sequence costs continue to drop, investigations into the structural diversity of the barley genome seem within reach. Here, we review the recent progress on pan-genomics in the model grass Brachypodium distachyon, and the cereal crops rice and maize, and devise a multi-tiered strategy for a pan-genome project in barley. Our design involves: (1) the construction of high-quality de novo sequence assemblies for a small core set of representative genotypes, (2) short-read sequencing of a large diversity panel of genebank accessions to medium coverage and (3) the use of complementary methods such as chromosome-conformation capture sequencing and k-mer-based association genetics. The in silico representation of the barley pan-genome may inform about the mechanisms of structural genome evolution in the Triticeae and supplement quantitative genetics models of crop performance for better accuracy and predictive ability.

Genomic approaches for studying crop evolution

Understanding how crop plants evolved from their wild relatives and spread around the world can inform about the origins of agriculture. Here, we review how the rapid development of genomic resources and tools has made it possible to conduct genetic mapping and population genetic studies to unravel the molecular underpinnings of domestication and crop evolution in diverse crop species. We propose three future avenues for the study of crop evolution: establishment of high-quality reference genomes for crops and their wild relatives; genomic characterization of germplasm collections; and the adoption of novel methodologies such as archaeogenetics, epigenomics, and genome editing.

Construction and characterization of a bacterial artificial chromosome library for Gossypium mustelinum

A bacterial artificial chromosome (BAC) library for G. mustelinum Miers ex G. Watt (AD₄) was constructed. Intact nuclei from G. mustelinum (AD₄) were used to isolate high molecular weight DNA, which was partially cleaved with Hind III and cloned into pSMART BAC (Hind III) vectors. The BAC library consisted of 208,182 clones arrayed in 542 384-microtiter plates, with an average insert size of 121.72 kb ranging from 100 to 150 kb. About 2% of the clones did not contain inserts. Based on an estimated genome size of 2372 Mb for G. mustelinum, the BAC library was estimated to have a total coverage of 10.50 × genome equivalents. The high capacity library of G. mustelinum will serve as a giant gene resource for map-based cloning of quantitative trait loci or genes associated with important agronomic traits or resistance to Verticillium wilt, physical mapping and comparative genome analysis.

Construction of a map-based reference genome sequence for barley, Hordeum vulgare L

Barley (Hordeum vulgare L.) is a cereal grass mainly used as animal fodder and raw material for the malting industry. The map-based reference genome sequence of barley cv. ‘Morex’ was constructed by the International Barley Genome Sequencing Consortium (IBSC) using hierarchical shotgun sequencing. Here, we report the experimental and computational procedures to (i) sequence and assemble more than 80,000 bacterial artificial chromosome (BAC) clones along the minimum tiling path of a genome-wide physical map, (ii) find and validate overlaps between adjacent BACs, (iii) construct 4,265 non-redundant sequence scaffolds representing clusters of overlapping BACs, and (iv) order and orient these BAC clusters along the seven barley chromosomes using positional information provided by dense genetic maps, an optical map and chromosome conformation capture sequencing (Hi-C). Integrative access to these sequence and mapping resources is provided by the barley genome explorer (BARLEX).

Towards a whole-genome sequence for rye (Secale cereale L.)

We report on a whole-genome draft sequence of rye (Secale cereale L.). Rye is a diploid Triticeae species closely related to wheat and barley, and an important crop for food and feed in Central and Eastern Europe. Through whole-genome shotgun sequencing of the 7.9-Gbp genome of the winter rye inbred line Lo7 we obtained a de novo assembly represented by 1.29 million scaffolds covering a total length of 2.8 Gbp. Our reference sequence represents nearly the entire low-copy portion of the rye genome. This genome assembly was used to predict 27 784 rye gene models based on homology to sequenced grass genomes. Through resequencing of 10 rye inbred lines and one accession of the wild relative S. vavilovii, we discovered more than 90 million single nucleotide variants and short insertions/deletions in the rye genome. From these variants, we developed the high-density Rye600k genotyping array with 600 843 markers, which enabled anchoring the sequence contigs along a high-density genetic map and establishing a synteny-based virtual gene order. Genotyping data were used to characterize the diversity of rye breeding pools and genetic resources, and to obtain a genome-wide map of selection signals differentiating the divergent gene pools. This rye whole-genome sequence closes a gap in Triticeae genome research, and will be highly valuable for comparative genomics, functional studies and genome-based breeding in rye.

Rye B chromosomes encode a functional Argonaute-like protein with in vitro slicer activities similar to its A chromosome paralog

B chromosomes (Bs) are supernumerary, dispensable parts of the nuclear genome, which appear in many different species of eukaryote. So far, Bs have been considered to be genetically inert elements without any functional genes. Our comparative transcriptome analysis and the detection of active RNA polymerase II (RNAPII) in the proximity of B chromatin demonstrate that the Bs of rye (Secale cereale) contribute to the transcriptome. In total, 1954 and 1218 B-derived transcripts with an open reading frame were expressed in generative and vegetative tissues, respectively. In addition to B-derived transposable element transcripts, a high percentage of short transcripts without detectable similarity to known proteins and gene fragments from A chromosomes (As) were found, suggesting an ongoing gene erosion process. In vitro analysis of the A- and B-encoded AGO4B protein variants demonstrated that both possess RNA slicer activity. These data demonstrate unambiguously the presence of a functional AGO4B gene on Bs and that these Bs carry both functional protein coding genes and pseudogene copies. Thus, B-encoded genes may provide an additional level of gene control and complexity in combination with their related A-located genes. Hence, physiological effects, associated with the presence of Bs, may partly be explained by the activity of B-located (pseudo)genes.

A Homolog of Blade-On-Petiole 1 and 2 (BOP1/2) Controls Internode Length and Homeotic Changes of the Barley Inflorescence

Inflorescence architecture in small-grain cereals has a direct effect on yield and is an important selection target in breeding for yield improvement. We analyzed the recessive mutation laxatum-a (lax-a) in barley (Hordeum vulgare), which causes pleiotropic changes in spike development, resulting in (1) extended rachis internodes conferring a more relaxed inflorescence, (2) broadened base of the lemma awns, (3) thinner grains that are largely exposed due to reduced marginal growth of the palea and lemma, and (4) and homeotic conversion of lodicules into two stamenoid structures. Map-based cloning enforced by mapping-by-sequencing of the mutant lax-a locus enabled the identification of a homolog of BLADE-ON-PETIOLE1 (BOP1) and BOP2 as the causal gene. Interestingly, the recently identified barley uniculme4 gene also is a BOP1/2 homolog and has been shown to regulate tillering and leaf sheath development. While the Arabidopsis (Arabidopsis thaliana) BOP1 and BOP2 genes act redundantly, the barley genes contribute independent effects in specifying the developmental growth of vegetative and reproductive organs, respectively. Analysis of natural genetic diversity revealed strikingly different haplotype diversity for the two paralogous barley genes, likely affected by the respective genomic environments, since no indication for an active selection process was detected.