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

Full-Length Transcriptome Construction and Systematic Characterization of Virulence Factor-Associated Isoforms in Vairimorpha (Nosema) Ceranae

Vairimorpha (Nosema) ceranae is a single-cellular fungus that obligately infects the midgut epithelial cells of adult honeybees, causing bee microsporidiosis and jeopardizing bee health and production. This work aims to construct the full-length transcriptome of V. ceranae and conduct a relevant investigation using PacBio single-molecule real-time (SMRT) sequencing technology. Following PacBio SMRT sequencing, 41,950 circular consensus (CCS) were generated, and 25,068 full-length non-chimeric (FLNC) reads were then detected. After polishing, 4387 high-quality, full-length transcripts were gained. There are 778, 2083, 1202, 1559, 1457, 1232, 1702, and 3896 full-length transcripts that could be annotated to COG, GO, KEGG, KOG, Pfam, Swiss-Prot, eggNOG, and Nr databases, respectively. Additionally, 11 alternative splicing (AS) events occurred in 6 genes were identified, including 1 alternative 5' splice-site and 10 intron retention. The structures of 225 annotated genes in the V. ceranae reference genome were optimized, of which 29 genes were extended at both 5' UTR and 3' UTR, while 90 and 106 genes were, respectively, extended at the 5' UTR as well as 3' UTR. Furthermore, a total of 29 high-confidence lncRNAs were obtained, including 12 sense-lncRNAs, 10 lincRNAs, and 7 antisense-lncRNAs. Taken together, the high-quality, full-length transcriptome of V. ceranae was constructed and annotated, the structures of annotated genes in the V. ceranae reference genome were improved, and abundant new genes, transcripts, and lncRNAs were discovered. Findings from this current work offer a valuable resource and a crucial foundation for molecular and omics research on V. ceranae.

Chromosome genomics uncovers plant genome organization and function

The identification of causal genomic loci and their interactions underlying various traits in plants has been greatly aided by progress in understanding the organization of the nuclear genome. This provides clues to the responses of plants to environmental stimuli at the molecular level. Apart from other uses, these insights are needed to fully explore the potential of new breeding techniques that rely on genome editing. However, genome analysis and sequencing is not straightforward in the many agricultural crops and their wild relatives that possess large and complex genomes. Chromosome genomics streamlines this task by dissecting the genome to single chromosomes whose DNA is then used instead of nuclear DNA. This results in a massive and lossless reduction in DNA sample complexity, reduces the time and cost of the experiment, and simplifies data interpretation. Flow cytometric sorting of condensed mitotic chromosomes makes it possible to purify single chromosomes in large quantities, and as the DNA remains intact this process can be coupled successfully with many techniques in molecular biology and genomics. Since the first experiments with flow cytometric sorting in the late 1980s, numerous applications have been developed, and chromosome genomics has been having a significant impact in many areas of research, including the sequencing of complex genomes of important crops and gene cloning. This review discusses these applications, describes their contribution to advancements in plant genome analysis and gene cloning, and outlines future directions.

Reduced chromatin accessibility underlies gene expression differences in homologous chromosome arms of diploid Aegilops tauschii and hexaploid wheat

Background

Polyploidy is centrally important in the evolution and domestication of plants because it leads to major genomic changes, such as altered patterns of gene expression, which are thought to underlie the emergence of new traits. Despite the common occurrence of these globally altered patterns of gene expression in polyploids, the mechanisms involved are not well understood.

Results

Using a precisely defined framework of highly conserved syntenic genes on hexaploid wheat chromosome 3DL and its progenitor 3 L chromosome arm of diploid Aegilops tauschii, we show that 70% of these gene pairs exhibited proportionately reduced gene expression, in which expression in the hexaploid context of the 3DL genes was ∼40% of the levels observed in diploid Ae tauschii. Several genes showed elevated expression during the later stages of grain development in wheat compared with Ae tauschii. Gene sequence and methylation differences probably accounted for only a few cases of differences in gene expression. In contrast, chromosome-wide patterns of reduced chromatin accessibility of genes in the hexaploid chromosome arm compared with its diploid progenitor were correlated with both reduced gene expression and the imposition of new patterns of gene expression.

Conclusions

Our pilot-scale analyses show that chromatin compaction may orchestrate reduced gene expression levels in the hexaploid chromosome arm of wheat compared to its diploid progenitor chromosome arm.

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.

Pairing and Exchanging between Daypyrum villosum Chromosomes 6V#2 and 6V#4 in the Hybrids of Two Different Wheat Alien Substitution Lines

Normal pairing and exchanging is an important basis to evaluate the genetic relationship between homologous chromosomes in a wheat background. The pairing behavior between 6V#2 and 6V#4, two chromosomes from different Dasypyrum villosum accessions, is still not clear. In this study, two wheat alien substitution lines, 6V#2 (6A) and 6V#4 (6D), were crossed to obtain the F₁ hybrids and F₂ segregating populations, and the testcross populations were obtained by using the F₁ as a parent crossed with wheat variety Wan7107. The chromosomal behavior at meiosis in pollen mother cells (PMCs) of the F₁ hybrids was observed using a genomic in situ hybridization (GISH) technique. Exchange events of two alien chromosomes were investigated in the F₂ populations using nine polymerase chain reaction (PCR) markers located on the 6V short arm. The results showed that the two alien chromosomes could pair with each other to form ring- or rod-shaped bivalent chromosomes in 79.76% of the total PMCs, and most were pulled to two poles evenly at anaphase I. Investigation of the F₂ populations showed that the segregation ratios of seven markers were consistent with the theoretical values 3:1 or 1:2:1, and recombinants among markers were detected. A genetic linkage map of nine PCR markers for 6VS was accordingly constructed based on the exchange frequencies and compared with the physical maps of wheat and barley based on homologous sequences of the markers, which showed that conservation of sequence order compared to 6V was 6H and 6B > 6A > 6D. In the testcross populations with 482 plants, seven showed susceptibility to powdery mildew (PM) and lacked amplification of alien chromosomal bands. Six other plants had amplification of specific bands of both the alien chromosomes at multiple sites, which suggested that the alien chromosomes had abnormal separation behavior in about 1.5% of the PMCs in F₁, which resulted in some gametes containing two alien chromosomes. In addition, three new types of chromosome substitution were developed. This study lays a foundation for alien allelism tests and further assessment of the genetic relationship among 6V#2, 6V#4, and their wheat homoeologous chromosomes.

High-throughput long paired-end sequencing of a Fosmid library by PacBio

BACKGROUND

Large insert paired-end sequencing technologies are important tools for assembling genomes, delineating associated breakpoints and detecting structural rearrangements. To facilitate the comprehensive detection of inter- and intra-chromosomal structural rearrangements or variants (SVs) and complex genome assembly with long repeats and segmental duplications, we developed a new method based on single-molecule real-time synthesis sequencing technology for generating long paired-end sequences of large insert DNA libraries.

RESULTS

A Fosmid vector, pHZAUFOS3, was developed with the following new features: (1) two 18-bp non-palindromic I-SceI sites flank the cloning site, and another two sites are present in the skeleton of the vector, allowing long DNA inserts (and the long paired-ends in this paper) to be recovered as single fragments and the vector (~ 8 kb) to be fragmented into 2-3 kb fragments by I-SceI digestion and therefore was effectively removed from the long paired-ends (5-10 kb); (2) the chloramphenicol (Cm) resistance gene and replicon (oriV), necessary for colony growth, are located near the two sides of the cloning site, helping to increase the proportion of the paired-end fragments to single-end fragments in the paired-end libraries. Paired-end libraries were constructed by ligating the size-selected, mechanically sheared pooled Fosmid DNA fragments to the Ampicillin (Amp) resistance gene fragment and screening the colonies with Cm and Amp. We tested this method on yeast and Setaria italica Yugu1. Fosmid-size paired-ends with an average length longer than 2 kb for each end were generated. The N50 scaffold lengths of the de novo assemblies of the yeast and S. italica Yugu1 genomes were significantly improved. Five large and five small structural rearrangements or assembly errors spanning tens of bp to tens of kb were identified in S. italica Yugu1 including deletions, inversions, duplications and translocations.

CONCLUSIONS

We developed a new method for long paired-end sequencing of large insert libraries, which can efficiently improve the quality of de novo genome assembly and identify large and small structural rearrangements or assembly errors.

Physical location of tandem repeats in the wheat genome and application for chromosome identification

A general distribution of tandem repeats (TRs) in the wheat genome was predicted and a new web page combined with fluorescence in situ hybridization experiments, and the newly developed Oligo probes will improve the resolution for wheat chromosome identification. Comprehensive sequence analysis of tandem repeats (TR) in the wheat reference genome permits discovery and application of TRs for chromosome identification. Genome-wide localization of TRs was identified in the reference sequences of Chinese Spring using Tandem Repeat Finder (TRF). A database of repeats unit size, array number, and physical coverage length of TRs in the wheat genome was built. The distribution of TRs occupied 3-5% of the wheat chromosomes, with non-random dispersal across the A, B, and D genomes. Three classes of TRs surrounding the predicted genes were compared. An optimized computer-assisted website page B2DSC was constructed for the general distribution and chromosomally enriched zones of TR sequences to be displayed graphically. The physical distribution of predicted TRs in the wheat genome by B2DSC matched well with the corresponding hybridization signals obtained with fluorescence in situ hybridization (FISH). We developed 20 oligonucleotide probes representing 20-60 bp lengths of high copy number of TRs and verified by FISH. An integrated physical map of TR-Oligo probes for wheat chromosome identification was constructed. Our results suggest that the combination of both molecular cytogenetics and genomic research will significantly benefit wheat breeding through chromosome manipulation and engineering.

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.

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.

Shifting the limits in wheat research and breeding using a fully annotated reference genome

An annotated reference sequence representing the hexaploid bread wheat genome in 21 pseudomolecules has been analyzed to identify the distribution and genomic context of coding and noncoding elements across the A, B, and D subgenomes. With an estimated coverage of 94% of the genome and containing 107,891 high-confidence gene models, this assembly enabled the discovery of tissue- and developmental stage-related coexpression networks by providing a transcriptome atlas representing major stages of wheat development. Dynamics of complex gene families involved in environmental adaptation and end-use quality were revealed at subgenome resolution and contextualized to known agronomic single-gene or quantitative trait loci. This community resource establishes the foundation for accelerating wheat research and application through improved understanding of wheat biology and genomics-assisted breeding.