Mapping-by-sequencing in complex polyploid genomes using genic sequence capture: a case study to map yellow rust resistance in hexaploid wheat

Transcriptome Analysis Reveals Potential Mechanism in Storage Protein Trafficking within Developing Grains of Common Wheat

Gluten proteins are the major storage protein fraction in the mature wheat grain. They are restricted to the starchy endosperm, which defines the viscoelastic properties of wheat dough. The synthesis of these storage proteins is controlled by the endoplasmic reticulum (ER) and is directed into the vacuole via the Golgi apparatus. In the present study, transcriptome analysis was used to explore the potential mechanism within critical stages of grain development of wheat cultivar "Shaannong 33" and its sister line used as the control (CK). Samples were collected at 10 DPA (days after anthesis), 14 DPA, 20 DPA, and 30 DPA for transcriptomic analysis. The comparative transcriptome analysis identified that a total of 18,875 genes were differentially expressed genes (DEGs) between grains of four groups "T10 vs. CK10, T14 vs. CK14, T20 vs. CK20, and T30 vs. CK30", including 2824 up-regulated and 5423 down-regulated genes in T30 vs. CK30. Further, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment highlighted the maximum number of genes regulating protein processing in the endoplasmic reticulum (ER) during grain enlargement stages (10-20 DPA). In addition, KEGG database analysis reported 1362 and 788 DEGs involved in translation, ribosomal structure, biogenesis, flavonoid biosynthesis pathway and intracellular trafficking, secretion, and vesicular transport through protein processing within ER pathway (ko04141). Notably, consistent with the higher expression of intercellular storage protein trafficking genes at the initial 10 DPA, there was relatively low expression at later stages. Expression levels of nine randomly selected genes were verified by qRT-PCR, which were consistent with the transcriptome data. These data suggested that the initial stages of "cell division" played a significant role in protein quality control within the ER, thus maintaining the protein quality characteristics at grain maturity. Furthermore, our data suggested that the protein synthesis, folding, and trafficking pathways directed by a different number of genes during the grain enlargement stage contributed to the observed high-quality characteristics of gluten protein in Shaannong 33 (Triticum aestivum L.).

Genetic Analysis of Adult Plant Resistance to Stripe Rust in Common Wheat Cultivar "Pascal"

Wheat stripe rust is an important foliar disease that affects the wheat yield globally. Breeding for resistant wheat varieties is one of the most economically and environmentally effective ways to control this disease. The common wheat (Triticum aestivum L.) cultivar "Pascal" exhibited susceptibility to stripe rust at the seedling stage but it showed high resistance to stripe rust at the adult plant stage over 20 years in Gansu, a hotspot of the disease in northwestern China. To understand the genetic mechanism of stripe rust resistance in this cultivar, a 55K SNP array was used to analyze the two parents and the 220 recombinant inbred lines (RILs) derived from the cross of "Huixianhong" × "Pascal." We detected three new stripe rust adult plant resistance (APR) quantitative trait locus (QTL) contributed by Pascal, viz. QYr.gaas-1AL, QYr.gaas-3DL, and QYr.gaas-5AS, using the inclusive composite interval mapping method. They were flanked by SNP markers AX-111218361-AX-110577861, AX-111460455-AX-108798599, and AX-111523523-AX-110028503, respectively, and explained the phenotypic variation ranging from 11.0 to 23.1%. Bulked segregant exome capture sequencing (BSE-Seq) was used for fine mapping of QYr.gaas-1AL and selection of candidate genes, TraesCS1A02G313700, TraesCS1A02G313800, and TraesCS1A02G314900 for QYr.gaas-1AL. KASP markers BSE-1A-12 and HXPA-3D for QYr.gaas-1AL and QYr.gaas-3DL were developed for breeders to develop durable stripe rust-resistant wheat varieties.

Whole-exome sequencing of selected bread wheat recombinant inbred lines as a useful resource for allele mining and bulked segregant analysis

Although wheat (Triticum aestivum L.) is the main staple crop in the world and a major source of carbohydrates and proteins, functional genomics and allele mining are still big challenges. Given the advances in next-generation sequencing (NGS) technologies, the identification of causal variants associated with a target phenotype has become feasible. For these reasons, here, by combining sequence capture and target-enrichment methods with high-throughput NGS re-sequencing, we were able to scan at exome-wide level 46 randomly selected bread wheat individuals from a recombinant inbred line population and to identify and classify a large number of single nucleotide polymorphisms (SNPs). For technical validation of results, eight randomly selected SNPs were converted into Kompetitive Allele-Specific PCR (KASP) markers. This resource was established as an accessible and reusable molecular toolkit for allele data mining. The dataset we are making available could be exploited for novel studies on bread wheat genetics and as a foundation for starting breeding programs aimed at improving different key agronomic traits.

Delaying or delivering: identification of novel NAM-1 alleles that delay senescence to extend wheat grain fill duration

Senescence is a complex trait under genetic and environmental control, in which resources are remobilized from vegetative tissue into grain. Delayed senescence, or 'staygreen' traits, can confer stress tolerance, with extended photosynthetic activity hypothetically sustaining grain filling. The genetics of senescence regulation are largely unknown, with senescence variation often correlated with phenological traits. Here, we confirm staygreen phenotypes of two Triticum aestivum cv. Paragon ethyl methane sulfonate mutants previously identified during a forward genetic screen and selected for their agronomic performance, similar phenology, and differential senescence phenotypes. Grain filling experiments confirmed a positive relationship between onset of senescence and grain fill duration, reporting an associated ~14% increase in final dry grain weight for one mutant (P<0.05). Recombinant inbred line (RIL) populations segregating for the timing of senescence were developed for trait mapping purposes and phenotyped over multiple years under field conditions. Quantification and comparison of senescence metrics aided RIL selection, facilitating exome capture-enabled bulk segregant analysis (BSA). Using BSA we mapped our two staygreen traits to two independent, dominant, loci of 4.8 and 16.7 Mb in size encompassing 56 and 142 genes, respectively. Combining association analysis with variant effect prediction, we identified single nucleotide polymorphisms encoding self-validating mutations located in NAM-1 homoeologues, which we propose as gene candidates.

A rare single nucleotide variant in Pm5e confers powdery mildew resistance in common wheat

Powdery mildew poses severe threats to wheat production. The most sustainable way to control this disease is through planting resistant cultivars. We report the map-based cloning of the powdery mildew resistance allele Pm5e from a Chinese wheat landrace. We applied a two-step bulked segregant RNA sequencing (BSR-Seq) approach in developing tightly linked or co-segregating markers to Pm5e. The first BSR-Seq used phenotypically contrasting bulks of recombinant inbred lines (RILs) to identify Pm5e-linked markers. The second BSR-Seq utilized bulks of genetic recombinants screened from a fine-mapping population to precisely quantify the associated genomic variation in the mapping interval, and identified the Pm5e candidate genes. The function of Pm5e was validated by transgenic assay, loss-of-function mutants and haplotype association analysis. Pm5e encodes a nucleotide-binding domain leucine-rich-repeat-containing (NLR) protein. A rare nonsynonymous single nucleotide variant (SNV) within the C-terminal leucine rich repeat (LRR) domain is responsible for the gain of powdery mildew resistance function of Pm5e, an allele endemic to wheat landraces of Shaanxi province of China. Results from this study demonstrate the value of landraces in discovering useful genes for modern wheat breeding. The key SNV associated with powdery mildew resistance will be useful for marker-assisted selection of Pm5e in wheat breeding programs.

A framework for gene mapping in wheat demonstrated using the Yr7 yellow rust resistance gene

We used three approaches to map the yellow rust resistance gene Yr7 and identify associated SNPs in wheat. First, we used a traditional QTL mapping approach using a double haploid (DH) population and mapped Yr7 to a low-recombination region of chromosome 2B. To fine map the QTL, we then used an association mapping panel. Both populations were SNP array genotyped allowing alignment of QTL and genome-wide association scans based on common segregating SNPs. Analysis of the association panel spanning the QTL interval, narrowed the interval down to a single haplotype block. Finally, we used mapping-by-sequencing of resistant and susceptible DH bulks to identify a candidate gene in the interval showing high homology to a previously suggested Yr7 candidate and to populate the Yr7 interval with a higher density of polymorphisms. We highlight the power of combining mapping-by-sequencing, delivering a complete list of gene-based segregating polymorphisms in the interval with the high recombination, low LD precision of the association mapping panel. Our mapping-by-sequencing methodology is applicable to any trait and our results validate the approach in wheat, where with a near complete reference genome sequence, we are able to define a small interval containing the causative gene.

Understanding DNA Methylation Patterns in Wheat

The bread wheat genome is large (17 Gb), allohexaploid, and highly repetitive (80-90% of the genome), which makes genomic and epigenomic analyses expensive to conduct and a challenge to analyze. Here we provide an overview of recent bioinformatic and experimental methods that have been developed to understand DNA methylation patterns in the complex polyploid genome of wheat.

Combining a New Exome Capture Panel With an Effective varBScore Algorithm Accelerates BSA-Based Gene Cloning in Wheat

The discovery of functional genes underlying agronomic traits is of great importance for wheat improvement. Here we designed a new wheat exome capture probe panel based on IWGSC RefSeq v1.0 genome sequence information and developed an effective algorithm, varBScore, that can sufficiently reduce the background noise in gene mapping and identification. An effective method, termed bulked segregant exome capture sequencing (BSE-Seq) for identifying causal mutations or candidate genes was established by combining the use of a newly designed wheat exome capture panel, sequencing of bulked segregant pools from segregating populations, and the robust algorithm varBScore. We evaluated the effectiveness of varBScore on SNP calling using the published dataset for mapping and cloning the yellow rust resistance gene Yr7 in wheat. Furthermore, using BSE-Seq, we rapidly identified a wheat yellow leaf mutant gene, ygl1, in an ethyl methanesulfonate (EMS) mutant population and found that a single mutation of G to A at 921 position in the wild type YGL1 gene encoding magnesium-chelatase subunit chlI caused the leaf yellowing phenotype. We further showed that mutation of YGL1 through CRISPR/Cas9 gene editing led to a yellow phenotype on the leaves of transgenic wheat, indicating that ygl1 is the correct causal gene responsible for the mutant phenotype. In summary, our approach is highly efficient for discovering causal mutations and gene cloning in wheat.

Integrating genomic resources to present full gene and putative promoter capture probe sets for bread wheat

Background

Whole-genome shotgun resequencing of wheat is expensive because of its large, repetitive genome. Moreover, sequence data can fail to map uniquely to the reference genome, making it difficult to unambiguously assign variation. Resequencing using target capture enables sequencing of large numbers of individuals at high coverage to reliably identify variants associated with important agronomic traits. Previous studies have implemented complementary DNA/exon or gene-based probe sets in which the promoter and intron sequence is largely missing alongside newly characterized genes from the recent improved reference sequences.

Results

We present and validate 2 gold standard capture probe sets for hexaploid bread wheat, a gene and a putative promoter capture, which are designed using recently developed genome sequence and annotation resources. The captures can be combined or used independently. We demonstrate that the capture probe sets effectively enrich the high-confidence genes and putative promoter regions that were identified in the genome alongside a large proportion of the low-confidence genes and associated promoters. Finally, we demonstrate successful sample multiplexing that allows generation of adequate sequence coverage for single-nucleotide polymorphism calling while significantly reducing cost per sample for gene and putative promoter capture.

Conclusions

We show that a capture design employing an “island strategy” can enable analysis of the large gene/putative promoter space of wheat with only 2 × 160 Mbp probe sets. Furthermore, these assays extend the regions of the wheat genome that are amenable to analyses beyond its exome, providing tools for detailed characterization of these regulatory regions in large populations.

Integrating genomic resources to present full gene and putative promoter capture probe sets for bread wheat

BACKGROUND

Whole-genome shotgun resequencing of wheat is expensive because of its large, repetitive genome. Moreover, sequence data can fail to map uniquely to the reference genome, making it difficult to unambiguously assign variation. Resequencing using target capture enables sequencing of large numbers of individuals at high coverage to reliably identify variants associated with important agronomic traits. Previous studies have implemented complementary DNA/exon or gene-based probe sets in which the promoter and intron sequence is largely missing alongside newly characterized genes from the recent improved reference sequences.

RESULTS

We present and validate 2 gold standard capture probe sets for hexaploid bread wheat, a gene and a putative promoter capture, which are designed using recently developed genome sequence and annotation resources. The captures can be combined or used independently. We demonstrate that the capture probe sets effectively enrich the high-confidence genes and putative promoter regions that were identified in the genome alongside a large proportion of the low-confidence genes and associated promoters. Finally, we demonstrate successful sample multiplexing that allows generation of adequate sequence coverage for single-nucleotide polymorphism calling while significantly reducing cost per sample for gene and putative promoter capture.

CONCLUSIONS

We show that a capture design employing an "island strategy" can enable analysis of the large gene/putative promoter space of wheat with only 2 × 160 Mbp probe sets. Furthermore, these assays extend the regions of the wheat genome that are amenable to analyses beyond its exome, providing tools for detailed characterization of these regulatory regions in large populations.

Integrating genomic resources to present full gene and putative promoter capture probe sets for bread wheat

BACKGROUND

Whole-genome shotgun resequencing of wheat is expensive because of its large, repetitive genome. Moreover, sequence data can fail to map uniquely to the reference genome, making it difficult to unambiguously assign variation. Resequencing using target capture enables sequencing of large numbers of individuals at high coverage to reliably identify variants associated with important agronomic traits. Previous studies have implemented complementary DNA/exon or gene-based probe sets in which the promoter and intron sequence is largely missing alongside newly characterized genes from the recent improved reference sequences.

RESULTS

We present and validate 2 gold standard capture probe sets for hexaploid bread wheat, a gene and a putative promoter capture, which are designed using recently developed genome sequence and annotation resources. The captures can be combined or used independently. We demonstrate that the capture probe sets effectively enrich the high-confidence genes and putative promoter regions that were identified in the genome alongside a large proportion of the low-confidence genes and associated promoters. Finally, we demonstrate successful sample multiplexing that allows generation of adequate sequence coverage for single-nucleotide polymorphism calling while significantly reducing cost per sample for gene and putative promoter capture.

CONCLUSIONS

We show that a capture design employing an "island strategy" can enable analysis of the large gene/putative promoter space of wheat with only 2 × 160 Mbp probe sets. Furthermore, these assays extend the regions of the wheat genome that are amenable to analyses beyond its exome, providing tools for detailed characterization of these regulatory regions in large populations.

Applying the latest advances in genomics and phenomics for trait discovery in polyploid wheat

Improving traits in wheat has historically been challenging due to its large and polyploid genome, limited genetic diversity and in-field phenotyping constraints. However, within recent years many of these barriers have been lowered. The availability of a chromosome-level assembly of the wheat genome now facilitates a step-change in wheat genetics and provides a common platform for resources, including variation data, gene expression data and genetic markers. The development of sequenced mutant populations and gene-editing techniques now enables the rapid assessment of gene function in wheat directly. The ability to alter gene function in a targeted manner will unmask the effects of homoeolog redundancy and allow the hidden potential of this polyploid genome to be discovered. New techniques to identify and exploit the genetic diversity within wheat wild relatives now enable wheat breeders to take advantage of these additional sources of variation to address challenges facing food production. Finally, advances in phenomics have unlocked rapid screening of populations for many traits of interest both in greenhouses and in the field. Looking forwards, integrating diverse data types, including genomic, epigenetic and phenomics data, will take advantage of big data approaches including machine learning to understand trait biology in wheat in unprecedented detail.

Multiple Facets of Marine Invertebrate Conservation Genomics

Conservation genomics aims to preserve the viability of populations and the biodiversity of living organisms. Invertebrate organisms represent 95% of animal biodiversity; however, few genomic resources currently exist for the group. The subset of marine invertebrates includes the most ancient metazoan lineages and possesses codes for unique gene products and possible keys to adaptation. The benefits of supporting invertebrate conservation genomics research (e.g., likely discovery of novel genes, protein regulatory mechanisms, genomic innovations, and transposable elements) outweigh the various hurdles (rare, small, or polymorphic starting materials). Here we review best conservation genomics practices in the laboratory and in silico when applied to marine invertebrates and also showcase unique features in several case studies of acroporid corals, crown-of-thorns starfish, apple snails, and abalone. Marine conservation genomics should also address how diversity can lead to unique marine innovations, the impact of deleterious variation, and how genomic monitoring and profiling could positively affect broader conservation goals (e.g., value of baseline data for in situ/ex situ genomic stocks).

Hidden variation in polyploid wheat drives local adaptation

Wheat has been domesticated into a large number of agricultural environments and has the ability to adapt to diverse environments. To understand this process, we survey genotype, repeat content, and DNA methylation across a bread wheat landrace collection representing global genetic diversity. We identify independent variation in methylation, genotype, and transposon copy number. We show that these, so far unexploited, sources of variation have had a significant impact on the wheat genome and that ancestral methylation states become preferentially "hard coded" as single nucleotide polymorphisms (SNPs) via 5-methylcytosine deamination. These mechanisms also drive local adaption, impacting important traits such as heading date and salt tolerance. Methylation and transposon diversity could therefore be used alongside SNP-based markers for breeding.

Mapping Late Leaf Spot Resistance in Peanut (Arachis hypogaea) Using QTL-seq Reveals Markers for Marker-Assisted Selection

Late leaf spot (LLS; Cercosporidium personatum) is a major fungal disease of cultivated peanut (Arachis hypogaea). A recombinant inbred line population segregating for quantitative field resistance was used to identify quantitative trait loci (QTL) using QTL-seq. High rates of false positive SNP calls using established methods in this allotetraploid crop obscured significant QTLs. To resolve this problem, robust parental SNPs were first identified using polyploid-specific SNP identification pipelines, leading to discovery of significant QTLs for LLS resistance. These QTLs were confirmed over 4 years of field data. Selection with markers linked to these QTLs resulted in a significant increase in resistance, showing that these markers can be immediately applied in breeding programs. This study demonstrates that QTL-seq can be used to rapidly identify QTLs controlling highly quantitative traits in polyploid crops with complex genomes. Markers identified can then be deployed in breeding programs, increasing the efficiency of selection using molecular tools. Key Message: Field resistance to late leaf spot is a quantitative trait controlled by many QTLs. Using polyploid-specific methods, QTL-seq is faster and more cost effective than QTL mapping.

Mapping causal mutations by exome sequencing in a wheat TILLING population: a tall mutant case study

Forward genetic screens of induced mutant plant populations are powerful tools to identify genes underlying phenotypes of interest. Using traditional techniques, mapping causative mutations from forward screens is a lengthy, multi-step process, requiring the identification of a broad genetic region followed by candidate gene sequencing to characterize the causal variant. Mapping by whole genome sequencing accelerates the identification of causal mutations by simultaneously defining a mapping region and providing information on the induced genetic variants. In wheat, although the availability of a high-quality draft genome assembly facilitates mapping and mutation calling, whole genome resequencing remains prohibitively expensive due to its large genome. In the current study, we used exome sequencing as a complexity reduction strategy to detect mutations associated with a target phenotype. In a segregating wheat EMS population, we identified a clear peak region on chromosome arm 4BS associated with increased plant height. Although none of the significant SNPs seemed causative for the mutant phenotype, they were sufficient to identify a linked ~ 1.9 Mb deletion encompassing nine genes. These genes included Rht-B1, which is known to have a strong effect on plant height and is a strong candidate for the observed phenotype. We performed simulation experiments to determine the impacts of sequencing depth and bulk size and discuss the importance of considering each factor when designing mapping-by-sequencing experiments in wheat. This approach can accelerate the identification of candidate causal point mutations or linked deletions underlying important phenotypes.

Mapping causal mutations by exome sequencing in a wheat TILLING population: a tall mutant case study

Forward genetic screens of induced mutant plant populations are powerful tools to identify genes underlying phenotypes of interest. Using traditional techniques, mapping causative mutations from forward screens is a lengthy, multi-step process, requiring the identification of a broad genetic region followed by candidate gene sequencing to characterize the causal variant. Mapping by whole genome sequencing accelerates the identification of causal mutations by simultaneously defining a mapping region and providing information on the induced genetic variants. In wheat, although the availability of a high-quality draft genome assembly facilitates mapping and mutation calling, whole genome resequencing remains prohibitively expensive due to its large genome. In the current study, we used exome sequencing as a complexity reduction strategy to detect mutations associated with a target phenotype. In a segregating wheat EMS population, we identified a clear peak region on chromosome arm 4BS associated with increased plant height. Although none of the significant SNPs seemed causative for the mutant phenotype, they were sufficient to identify a linked ~ 1.9 Mb deletion encompassing nine genes. These genes included Rht-B1, which is known to have a strong effect on plant height and is a strong candidate for the observed phenotype. We performed simulation experiments to determine the impacts of sequencing depth and bulk size and discuss the importance of considering each factor when designing mapping-by-sequencing experiments in wheat. This approach can accelerate the identification of candidate causal point mutations or linked deletions underlying important phenotypes.

Rapid cloning of genes in hexaploid wheat using cultivar-specific long-range chromosome assembly

Cereal crops such as wheat and maize have large repeat-rich genomes that make cloning of individual genes challenging. Moreover, gene order and gene sequences often differ substantially between cultivars of the same crop species. A major bottleneck for gene cloning in cereals is the generation of high-quality sequence information from a cultivar of interest. In order to accelerate gene cloning from any cropping line, we report 'targeted chromosome-based cloning via long-range assembly' (TACCA). TACCA combines lossless genome-complexity reduction via chromosome flow sorting with Chicago long-range linkage to assemble complex genomes. We applied TACCA to produce a high-quality (N50 of 9.76 Mb) de novo chromosome assembly of the wheat line CH Campala Lr22a in only 4 months. Using this assembly we cloned the broad-spectrum Lr22a leaf-rust resistance gene, using molecular marker information and ethyl methanesulfonate (EMS) mutants, and found that Lr22a encodes an intracellular immune receptor homologous to the Arabidopsis thaliana RPM1 protein.

Wheat genomics comes of age

Advances in wheat genomics have lagged behind other major cereals (e.g., rice and maize) due to its highly repetitive and large polyploid genome. Recent technological developments in sequencing and assembly methods, however, have largely overcome these barriers. The community now moves to an era centred on functional characterisation of the genome. This includes understanding sequence and structural variation as well as how information is integrated across multiple homoeologous genomes. This understanding promises to uncover variation previously hidden from natural and human selection due to the often observed functional redundancy between homoeologs. Key functional genomic resources will enable this, including sequenced mutant populations and gene editing technologies which are now available in wheat. Training the next-generation of genomics-enabled researchers will be essential to ensure these advances are quickly translated into farmers' fields.