Sunish K. Sehgal

Analysis of the bread wheat genome using whole-genome shotgun sequencing

Bread wheat (Triticum aestivum) is a globally important crop, accounting for 20 per cent of the calories consumed by humans. Major efforts are underway worldwide to increase wheat production by extending genetic diversity and analysing key traits, and genomic resources can accelerate progress. But so far the very large size and polyploid complexity of the bread wheat genome have been substantial barriers to genome analysis. Here we report the sequencing of its large, 17-gigabase-pair, hexaploid genome using 454 pyrosequencing, and comparison of this with the sequences of diploid ancestral and progenitor genomes. We identified between 94,000 and 96,000 genes, and assigned two-thirds to the three component genomes (A, B and D) of hexaploid wheat. High-resolution synteny maps identified many small disruptions to conserved gene order. We show that the hexaploid genome is highly dynamic, with significant loss of gene family members on polyploidization and domestication, and an abundance of gene fragments. Several classes of genes involved in energy harvesting, metabolism and growth are among expanded gene families that could be associated with crop productivity. Our analyses, coupled with the identification of extensive genetic variation, provide a resource for accelerating gene discovery and improving this major crop.

A 4-gigabase physical map unlocks the structure and evolution of the complex genome of Aegilops tauschii, the wheat D-genome progenitor

The current limitations in genome sequencing technology require the construction of physical maps for high-quality draft sequences of large plant genomes, such as that of Aegilops tauschii, the wheat D-genome progenitor. To construct a physical map of the Ae. tauschii genome, we fingerprinted 461,706 bacterial artificial chromosome clones, assembled contigs, designed a 10K Ae. tauschii Infinium SNP array, constructed a 7,185-marker genetic map, and anchored on the map contigs totaling 4.03 Gb. Using whole genome shotgun reads, we extended the SNP marker sequences and found 17,093 genes and gene fragments. We showed that collinearity of the Ae. tauschii genes with Brachypodium distachyon, rice, and sorghum decreased with phylogenetic distance and that structural genome evolution rates have been high across all investigated lineages in subfamily Pooideae, including that of Brachypodieae. We obtained additional information about the evolution of the seven Triticeae chromosomes from 12 ancestral chromosomes and uncovered a pattern of centromere inactivation accompanying nested chromosome insertions in grasses. We showed that the density of noncollinear genes along the Ae. tauschii chromosomes positively correlates with recombination rates, suggested a cause, and showed that new genes, exemplified by disease resistance genes, are preferentially located in high-recombination chromosome regions.

A whole-genome shotgun approach for assembling and anchoring the hexaploid bread wheat genome

Polyploid species have long been thought to be recalcitrant to whole-genome assembly. By combining high-throughput sequencing, recent developments in parallel computing, and genetic mapping, we derive, de novo, a sequence assembly representing 9.1 Gbp of the highly repetitive 16 Gbp genome of hexaploid wheat, Triticum aestivum, and assign 7.1 Gb of this assembly to chromosomal locations. The genome representation and accuracy of our assembly is comparable or even exceeds that of a chromosome-by-chromosome shotgun assembly. Our assembly and mapping strategy uses only short read sequencing technology and is applicable to any species where it is possible to construct a mapping population.

Comparative Analysis of Syntenic Genes in Grass Genomes Reveals Accelerated Rates of Gene Structure and Coding Sequence Evolution in Polyploid Wheat

Cycles of whole-genome duplication (WGD) and diploidization are hallmarks of eukaryotic genome evolution and speciation. Polyploid wheat (Triticum aestivum) has had a massive increase in genome size largely due to recent WGDs. How these processes may impact the dynamics of gene evolution was studied by comparing the patterns of gene structure changes, alternative splicing (AS), and codon substitution rates among wheat and model grass genomes. In orthologous gene sets, significantly more acquired and lost exonic sequences were detected in wheat than in model grasses. In wheat, 35% of these gene structure rearrangements resulted in frame-shift mutations and premature termination codons. An increased codon mutation rate in the wheat lineage compared with Brachypodium distachyon was found for 17% of orthologs. The discovery of premature termination codons in 38% of expressed genes was consistent with ongoing pseudogenization of the wheat genome. The rates of AS within the individual wheat subgenomes (21%–25%) were similar to diploid plants. However, we uncovered a high level of AS pattern divergence between the duplicated homeologous copies of genes. Our results are consistent with the accelerated accumulation of AS isoforms, nonsynonymous mutations, and gene structure rearrangements in the wheat lineage, likely due to genetic redundancy created by WGDs. Whereas these processes mostly contribute to the degeneration of a duplicated genome and its diploidization, they have the potential to facilitate the origin of new functional variations, which, upon selection in the evolutionary lineage, may play an important role in the origin of novel traits.

Tandem Amplification of a Chromosomal Segment Harboring EPSPS Locus Confers Glyphosate Resistance in Kochia scoparia

Genes encoding enolpyruvylshikimate phosphate synthase are tandemly arranged on chromosomes of field-evolved glyphosate-resistant Kochia scoparia. Recent rapid evolution and spread of resistance to the most extensively used herbicide, glyphosate, is a major threat to global crop production. Genetic mechanisms by which weeds evolve resistance to herbicides largely determine the level of resistance and the rate of evolution of resistance. In a previous study, we determined that glyphosate resistance in Kochia scoparia is due to the amplification of the 5-Enolpyruvylshikimate-3-Phosphate Synthase (EPSPS) gene, the enzyme target of glyphosate. Here, we investigated the genomic organization of the amplified EPSPS copies using fluorescence in situ hybridization (FISH) and extended DNA fiber (Fiber FISH) on K. scoparia chromosomes. In both glyphosate-resistant K. scoparia populations tested (GR1 and GR2), FISH results displayed a single and prominent hybridization site of the EPSPS gene localized on the distal end of one pair of homologous metaphase chromosomes compared with a faint hybridization site in glyphosate-susceptible samples (GS1 and GS2). Fiber FISH displayed 10 copies of the EPSPS gene (approximately 5 kb) arranged in tandem configuration approximately 40 to 70 kb apart, with one copy in an inverted orientation in GR2. In agreement with FISH results, segregation of EPSPS copies followed single-locus inheritance in GR1 population. This is the first report of tandem target gene amplification conferring field-evolved herbicide resistance in weed populations.

A diploid wheat TILLING resource for wheat functional genomics

BackgroundTriticum monococcum L., an A genome diploid einkorn wheat, was the first domesticated crop. As a diploid, it is attractive genetic model for the study of gene structure and function of wheat-specific traits. Diploid wheat is currently not amenable to reverse genetics approaches such as insertion mutagenesis and post-transcriptional gene silencing strategies. However, TILLING offers a powerful functional genetics approach for wheat gene analysis.ResultsWe developed a TILLING population of 1,532 M2 families using EMS as a mutagen. A total of 67 mutants were obtained for the four genes studied. Waxy gene mutation frequencies are known to be 1/17.6 - 34.4 kb DNA in polyploid wheat TILLING populations. The T. monococcum diploid wheat TILLING population had a mutation frequency of 1/90 kb for the same gene. Lignin biosynthesis pathway genes- COMT1, HCT2, and 4CL1 had mutation frequencies of 1/86 kb, 1/92 kb and 1/100 kb, respectively. The overall mutation frequency of the diploid wheat TILLING population was 1/92 kb.ConclusionThe mutation frequency of a diploid wheat TILLING population was found to be higher than that reported for other diploid grasses. The rate, however, is lower than tetraploid and hexaploid wheat TILLING populations because of the higher tolerance of polyploids to mutations. Unlike polyploid wheat, most mutants in diploid wheat have a phenotype amenable to forward and reverse genetic analysis and establish diploid wheat as an attractive model to study gene function in wheat. We estimate that a TILLING population of 5, 520 will be needed to get a non-sense mutation for every wheat gene of interest with 95% probability.

Tandem Amplification of a Chromosomal Segment Harboring 5-Enolpyruvylshikimate-3-Phosphate Synthase Locus Confers Glyphosate Resistance in Kochia scoparia

Genes encoding enolpyruvylshikimate phosphate synthase are tandemly arranged on chromosomes of field-evolved glyphosate-resistant Kochia scoparia. Recent rapid evolution and spread of resistance to the most extensively used herbicide, glyphosate, is a major threat to global crop production. Genetic mechanisms by which weeds evolve resistance to herbicides largely determine the level of resistance and the rate of evolution of resistance. In a previous study, we determined that glyphosate resistance in Kochia scoparia is due to the amplification of the 5-Enolpyruvylshikimate-3-Phosphate Synthase (EPSPS) gene, the enzyme target of glyphosate. Here, we investigated the genomic organization of the amplified EPSPS copies using fluorescence in situ hybridization (FISH) and extended DNA fiber (Fiber FISH) on K. scoparia chromosomes. In both glyphosate-resistant K. scoparia populations tested (GR1 and GR2), FISH results displayed a single and prominent hybridization site of the EPSPS gene localized on the distal end of one pair of homologous metaphase chromosomes compared with a faint hybridization site in glyphosate-susceptible samples (GS1 and GS2). Fiber FISH displayed 10 copies of the EPSPS gene (approximately 5 kb) arranged in tandem configuration approximately 40 to 70 kb apart, with one copy in an inverted orientation in GR2. In agreement with FISH results, segregation of EPSPS copies followed single-locus inheritance in GR1 population. This is the first report of tandem target gene amplification conferring field-evolved herbicide resistance in weed populations.

SNP Discovery for mapping alien introgressions in wheat.

BackgroundMonitoring alien introgressions in crop plants is difficult due to the lack of genetic and molecular mapping information on the wild crop relatives. The tertiary gene pool of wheat is a very important source of genetic variability for wheat improvement against biotic and abiotic stresses. By exploring the 5Mg short arm (5MgS) of Aegilops geniculata, we can apply chromosome genomics for the discovery of SNP markers and their use for monitoring alien introgressions in wheat (Triticum aestivum L).ResultsThe short arm of chromosome 5Mg of Ae. geniculata Roth (syn. Ae. ovata L.; 2n = 4x = 28, UgUgMgMg) was flow-sorted from a wheat line in which it is maintained as a telocentric chromosome. DNA of the sorted arm was amplified and sequenced using an Illumina Hiseq 2000 with ~45x coverage. The sequence data was used for SNP discovery against wheat homoeologous group-5 assemblies. A total of 2,178 unique, 5MgS-specific SNPs were discovered. Randomly selected samples of 59 5MgS-specific SNPs were tested (44 by KASPar assay and 15 by Sanger sequencing) and 84% were validated. Of the selected SNPs, 97% mapped to a chromosome 5Mg addition to wheat (the source of t5MgS), and 94% to 5Mg introgressed from a different accession of Ae. geniculata substituting for chromosome 5D of wheat. The validated SNPs also identified chromosome segments of 5MgS origin in a set of T5D-5Mg translocation lines; eight SNPs (25%) mapped to TA5601 [T5DL · 5DS-5MgS(0.75)] and three (8%) to TA5602 [T5DL · 5DS-5MgS (0.95)]. SNPs (gsnp_5ms83 and gsnp_5ms94), tagging chromosome T5DL · 5DS-5MgS(0.95) with the smallest introgression carrying resistance to leaf rust (Lr57) and stripe rust (Yr40), were validated in two released germplasm lines with Lr57 and Yr40 genes.ConclusionThis approach should be widely applicable for the identification of species/genome-specific SNPs. The development of a large number of SNP markers will facilitate the precise introgression and monitoring of alien segments in crop breeding programs and further enable mapping and cloning novel genes from the wild relatives of crop plants.

A distorted circadian clock causes early flowering and temperature-dependent variation in spike development in the Eps-3Am mutant of einkorn wheat.

Viable circadian clocks help organisms to synchronize their development with daily and seasonal changes, thereby providing both evolutionary fitness and advantage from an agricultural perspective. A high-resolution mapping approach combined with mutant analysis revealed a cereal ortholog of Arabidopsis thaliana LUX ARRHYTHMO/PHYTOCLOCK 1 (LUX/PCL1) as a promising candidate for the earliness per se 3 (Eps-3Am) locus in einkorn wheat (Triticum monococcum L.). Using delayed fluorescence measurements it was shown that Eps-3Am containing einkorn wheat accession KT3-5 had a distorted circadian clock. The hypothesis was subsequently confirmed by performing a time course study on central and output circadian clock genes, which showed arrhythmic transcript patterns in KT3-5 under constant ambient conditions, i.e., constant light and temperature. It was also demonstrated that variation in spikelet number between wild-type and mutants is sensitive to temperature, becoming negligible at 25°. These observations lead us to propose that the distorted clock is causative for both early flowering and variation in spike size and spikelet number, and that having a dysfunctional LUX could have neutral, or even positive, effects in warmer climates. To test the latter hypothesis we ascertained sequence variation of LUX in a range of wheat germplasm. We observed a higher variation in the LUX sequence among accessions coming from the warmer climate and a unique in-frame mutation in early-flowering Chinese T. turgidum cultivar ‘Tsing Hua no. 559.’ Our results emphasize the importance of the circadian clock in temperate cereals as a promising target for adaptation to new environments.

Chromosome arm-specific BAC end sequences permit comparative analysis of homoeologous chromosomes and genomes of polyploid wheat

BackgroundBread wheat, one of the world’s staple food crops, has the largest, highly repetitive and polyploid genome among the cereal crops. The wheat genome holds the key to crop genetic improvement against challenges such as climate change, environmental degradation, and water scarcity. To unravel the complex wheat genome, the International Wheat Genome Sequencing Consortium (IWGSC) is pursuing a chromosome- and chromosome arm-based approach to physical mapping and sequencing. Here we report on the use of a BAC library made from flow-sorted telosomic chromosome 3A short arm (t3AS) for marker development and analysis of sequence composition and comparative evolution of homoeologous genomes of hexaploid wheat.ResultsThe end-sequencing of 9,984 random BACs from a chromosome arm 3AS-specific library (TaaCsp3AShA) generated 11,014,359 bp of high quality sequence from 17,591 BAC-ends with an average length of 626 bp. The sequence represents 3.2% of t3AS with an average DNA sequence read every 19 kb. Overall, 79% of the sequence consisted of repetitive elements, 1.38% as coding regions (estimated 2,850 genes) and another 19% of unknown origin. Comparative sequence analysis suggested that 70-77% of the genes present in both 3A and 3B were syntenic with model species. Among the transposable elements, gypsy/sabrina (12.4%) was the most abundant repeat and was significantly more frequent in 3A compared to homoeologous chromosome 3B. Twenty novel repetitive sequences were also identified using de novo repeat identification. BESs were screened to identify simple sequence repeats (SSR) and transposable element junctions. A total of 1,057 SSRs were identified with a density of one per 10.4 kb, and 7,928 junctions between transposable elements (TE) and other sequences were identified with a density of one per 1.39 kb. With the objective of enhancing the marker density of chromosome 3AS, oligonucleotide primers were successfully designed from 758 SSRs and 695 Insertion Site Based Polymorphisms (ISBPs). Of the 96 ISBP primer pairs tested, 28 (29%) were 3A-specific and compared to 17 (18%) for 96 SSRs.ConclusionThis work reports on the use of wheat chromosome arm 3AS-specific BAC library for the targeted generation of sequence data from a particular region of the huge genome of wheat. A large quantity of sequences were generated from the A genome of hexaploid wheat for comparative genome analysis with homoeologous B and D genomes and other model grass genomes. Hundreds of molecular markers were developed from the 3AS arm-specific sequences; these and other sequences will be useful in gene discovery and physical mapping.