A chromosome bin map of 16,000 expressed sequence tag loci and distribution of genes among the three genomes of polyploid wheat

Structural and Functional Characterization at the Molecular Level of the MATE Gene Family in Wheat in Silico

A series of multidrug extransporters known as the multidrug and potentially toxic extrusion (MATE) genes are found in all living things and are crucial for the removal of heavy metal ions, metalloids, exogenous xenobiotics, endogenous secondary metabolites, and other toxic substances from the cells. However, there has only been a small amount of them in silico analysis of the MATE family of genes in plant species. In the current study, the MATE gene family was characterized in silico where two families and seven subfamilies based on their evolutionary relationships were proposed. Plant breeders may use TraesCS1D02G030400, TraesCS4B02G244400, and TraesCS1A02G029900 genes for marker-assisted or transgenic breeding to develop novel cultivars since these genes have been hypothesized from protein-protein interaction study to play a critical role in the transport of toxic chemicals across cells. The exon number varies from 01 to 14. One exon has TraesCS1A02G188100, TraesCS5B02G562500, TraesCS6A02G256400, and TraesCS6D02G384300 genes, while 14 exons have only two genes that are TraesCS6A02G418800 and TraesCS6D02G407900. Biological stress (infestations of disease) affects the expression of most of the MATE genes, with the gene TraesCS5D02G355500 having the highest expression level in the wheat expression browser tool. Using the Grain interpretation search engine tool, it is found that the vast bulk of MATE genes are voiced throughout biotic environmental stresses caused by disease pests, with the genotype TraesCS5B02G326600.1 from family 1 exhibiting the greatest level of expression throughout Fusarium head blight infection by Fusarium graminearum after 4 days of infection. The researchers constructed 39 ternary plots, each with a distinct degree of expression under biotic and abiotic stress settings, and observed that 44% of the triplets have imbalanced outputs (extreme values) due to their higher tissue specificity and increased intensity.

Fine mapping of QPm.caas-3BS, a stable QTL for adult-plant resistance to powdery mildew in wheat (Triticum aestivum L.)

A stable QTL QPm.caas-3BS for adult-plant resistance to powdery mildew was mapped in an interval of 431 kb, and candidate genes were predicted based on gene sequences and expression profiles. Powdery mildew is a devastating foliar disease occurring in most wheat-growing areas. Characterization and fine mapping of genes for powdery mildew resistance can benefit marker-assisted breeding. We previously identified a stable quantitative trait locus (QTL) QPm.caas-3BS for adult-plant resistance to powdery mildew in a recombinant inbred line population of Zhou8425B/Chinese Spring by phenotyping across four environments. Using 11 heterozygous recombinants and high-density molecular markers, QPm.caas-3BS was delimited in a physical interval of approximately 3.91 Mb. Based on re-sequenced data and expression profiles, three genes TraesCS3B02G014800, TraesCS3B02G016800 and TraesCS3B02G019900 were associated with the powdery mildew resistance locus. Three gene-specific kompetitive allele-specific PCR (KASP) markers were developed from these genes and validated in the Zhou8425B derivatives and Zhou8425B/Chinese Spring population in which the resistance gene was mapped to a 0.3 cM interval flanked by KASP14800 and snp_50465, corresponding to a 431 kb region at the distal end of chromosome 3BS. Within the interval, TraesCS3B02G014800 was the most likely candidate gene for QPm.caas-3BS, but TraesCS3B02G016300 and TraesCS3B02G016400 were less likely candidates based on gene annotations and sequence variation between the parents. These results not only offer high-throughput KASP markers for improvement of powdery mildew resistance but also pave the way to map-based cloning of the resistance gene.

Differences between diploid donors are the main contributing factor for subgenome asymmetry measured in either gene ratio or relative diversity in allopolyploids

Subgenome asymmetry (SA) has routinely been attributed to different responses between the subgenomes of a polyploid to various stimuli during evolution. Here, we compared subgenome differences in gene ratio and relative diversity between artificial and natural genotypes of several allopolyploid species. Surprisingly, consistent differences were not detected between these two types of polyploid genotypes, although they differ in times exposed to evolutionary selection. The estimated ratio of shared genes between a subgenome and its diploid donor was invariably higher for the artificial allopolyploid genotypes than those for the natural genotypes, which is expected as it is now well-known that many genes in a species are not shared among all individuals. As the exact diploid parent for a given subgenome is unknown, the estimated ratios of shared genes for the natural genotypes would also include difference among individual genotypes of the diploid donor species. Further, we detected the presence of SA in genotypes before the completion of the polyploidization events as well as in those which were not formed via polyploidization. These results indicate that SA may, to a large degree, reflect differences between its diploid donors or that changes occurred during polyploid evolution are defined by their donor genomes.

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.

Mapping of Major Fusarium Head Blight Resistance from Canadian Wheat cv. AAC Tenacious

Fusarium head blight (FHB) is one of the most devastating wheat disease due to its direct detrimental effects on grain-yield, quality and marketability. Resistant cultivars offer the most effective approach to manage FHB; however, the lack of different resistance resources is still a major bottleneck for wheat breeding programs. To identify and dissect FHB resistance, a doubled haploid wheat population produced from the Canadian spring wheat cvs AAC Innova and AAC Tenacious was phenotyped for FHB response variables incidence and severity, visual rating index (VRI), deoxynivalenol (DON) content, and agronomic traits days to anthesis (DTA) and plant height (PHT), followed by single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) marker genotyping. A high-density map was constructed consisting of 10,328 markers, mapped on all 21 chromosomes with a map density of 0.35 cM/marker. Together, two major quantitative trait loci for FHB resistance were identified on chromosome 2D from AAC Tenacious; one of these loci on 2DS also colocated with loci for DTA and PHT. Another major locus for PHT, which cosegregates with locus for low DON, was also identified along with many minor and epistatic loci. QTL identified from AAC Tenacious may be useful to pyramid FHB resistance.

Development of genome wide transposable elements based repeat junction markers in Jatropha (Jatropha curcas L.)

Jatropha curcas is a potential biodiesel crop and a highly adaptable species to various agro-climatic conditions. In this study, we have utilized transposable elements' (TE) repeat junctions (RJs) which are an important constituent of the genome, used to form a genome-wide molecular marker platform owing to its use in genomic studies of plants. We screened our previously generated Jatropha hybrid genome assembly of size 265 Mbp using RJPrimers pipeline software and identified a total of 1274 TE junctions. For the predicted RJs, we designed 2868 polymerase chain reaction (PCR) based RJ markers (RJMs) flanking the junction regions. In addition to marker design, the identified RJs were utilized to detect 225,517 TEs across the genome. The different types of transposable repeat elements mainly were scattered into Retro, LTR, Copia and Gypsy categories. The efficacy of the designed markers was tested by utilizing a subset of RJMs selected randomly. We have validated 96 randomly selected RJ primers in a group of 32 J. curcas genotypes and more than 90% of the markers effectively intensified as amplicons. Of these, 10 primers were shown to be polymorphic in estimating genetic diversity among the 32 Jatropha lines. UPGMA cluster analysis revealed the formation of two clusters such as A and B exhibiting 85.5% and 87% similarity coefficient respectively. The various RJMs identified in this study could be utilized as a significant asset in Jatropha functional genomics including genome determination, mapping and marker-assisted selection.

Partitioning and physical mapping of wheat chromosome 3B and its homoeologue 3E in Thinopyrum elongatum by inducing homoeologous recombination

We performed homoeologous recombination-based partitioning and physical mapping of wheat chromosome 3B and Th. elongatum chromosome 3E, providing a unique physical framework of this homoeologous pair for genome studies. The wheat (Triticum aestivum, 2n = 6x = 42, AABBDD) and Thinopyrum elongatum (2n = 2x = 14, EE) genomes can be differentiated from each other by fluorescent genomic in situ hybridization (FGISH) as well as molecular markers. This has facilitated homoeologous recombination-based partitioning and engineering of their genomes for physical mapping and alien introgression. Here, we constructed a special wheat genotype, which was double monosomic for wheat chromosome 3B and Th. elongatum chromosome 3E and homozygous for the ph1b mutant, to induce 3B-3E homoeologous recombination. Totally, 81 3B-3E recombinants were recovered and detected in the primary, secondary, and tertiary homoeologous recombination cycles by FGISH. Comparing to the primary recombination, the secondary and tertiary recombination shifted toward the proximal regions due to the increase in homology between the pairing partners. The 3B-3E recombinants were genotyped by high-throughput wheat 90-K single nucleotide polymorphism (SNP) arrays and their recombination breakpoints physically mapped based on the FGISH patterns and SNP results. The 3B-3E recombination physically partitioned chromosome 3B into 38 bins, and 429 SNPs were assigned to the distinct bins. Integrative analysis of FGISH and SNP results led to the construction of a composite bin map for chromosome 3B. Additionally, we developed 22 SNP-derived semi-thermal asymmetric reverse PCR markers specific for chromosome 3E and constructed a comparative map of homoeologous chromosomes 3E, 3B, 3A, and 3D. In summary, this work provides a unique physical framework for further studies of the 3B-3E homoeologous pair and diversifies the wheat genome for wheat improvement.

Genomic analysis of Spanish wheat landraces reveals their variability and potential for breeding

Background

One of the main goals of the plant breeding in the twenty-first century is the development of crop cultivars that can maintain current yields in unfavorable environments. Landraces that have been grown under varying local conditions include genetic diversity that will be essential to achieve this objective. The Center of Plant Genetic Resources of the Spanish Institute for Agriculture Research maintains a broad collection of wheat landraces. These accessions, which are locally adapted to diverse eco-climatic conditions, represent highly valuable materials for breeding. However, their efficient use requires an exhaustive genetic characterization. The overall aim of this study was to assess the diversity and population structure of a selected set of 380 Spanish landraces and 52 reference varieties of bread and durum wheat by high-throughput genotyping.

Results

The DArTseq GBS approach generated 10 K SNPs and 40 K high-quality DArT markers, which were located against the currently available bread and durum wheat reference genomes. The markers with known locations were distributed across all chromosomes with relatively well-balanced genome-wide coverage. The genetic analysis showed that the Spanish wheat landraces were clustered in different groups, thus representing genetic pools providing a range of allelic variation. The subspecies had a major impact on the population structure of the durum wheat landraces, with three distinct clusters that corresponded to subsp. durum, turgidum and dicoccon being identified. The population structure of bread wheat landraces was mainly biased by geographic origin.

Conclusions

The results showed broader genetic diversity in the landraces compared to a reference set that included commercial varieties, and higher divergence between the landraces and the reference set in durum wheat than in bread wheat. The analyses revealed genomic regions whose patterns of variation were markedly different in the landraces and reference varieties, indicating loci that have been under selection during crop improvement, which could help to target breeding efforts. The results obtained from this work will provide a basis for future genome-wide association studies.

Transgenic Wheat Harboring an RNAi Element Confers Dual Resistance Against Synergistically Interacting Wheat Streak Mosaic Virus and Triticum Mosaic Virus

Wheat streak mosaic virus (WSMV) and triticum mosaic virus (TriMV) are economically important viruses of wheat (Triticum aestivum L.), causing significant yield losses in the Great Plains region of the United States. These two viruses are transmitted by wheat curl mites, which often leads to mixed infections with synergistic interaction in grower fields that exacerbates yield losses. Development of dual-resistant wheat lines would provide effective control of these two viruses. In this study, a genetic resistance strategy employing an RNA interference (RNAi) approach was implemented by assembling a hairpin element composed of a 202-bp (404-bp in total) stem sequence of the NIb (replicase) gene from each of WSMV and TriMV in tandem and of an intron sequence in the loop. The derived RNAi element was cloned into a binary vector and was used to transform spring wheat genotype CB037. Phenotyping of T₁ lineages across eight independent transgenic events for resistance revealed that i) two of the transgenic events provided resistance to WSMV and TriMV, ii) four events provided resistance to either WSMV or TriMV, and iii) no resistance was found in two other events. T₂ populations derived from the two events classified as dual-resistant were subsequently monitored for stability of the resistance phenotype through the T₄ generation. The resistance phenotype in these events was temperature-dependent, with a complete dual resistance at temperatures ≥25°C and an increasingly susceptible response at temperatures below 25°C. Northern blot hybridization of total RNA from transgenic wheat revealed that virus-specific small RNAs (vsRNAs) accumulated progressively with an increase in temperature, with no detectable levels of vsRNA accumulation at 20°C. Thus, the resistance phenotype of wheat harboring an RNAi element was correlated with accumulation of vsRNAs, and the generation of vsRNAs can be used as a molecular marker for the prediction of resistant phenotypes of transgenic plants at a specific temperature.

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.

Production of a complete set of wheat-barley group-7 chromosome recombinants with increased grain β-glucan content

Wheat-barley group-7 recombinant chromosomes were selected using molecular cytogenetics and SNP markers; increased grain β-glucan content was observed in wheat plants with two and four copies of HvCslF6. The soluble dietary fiber (1-3)(1-4) mixed linked β-D-glucan from cereal grains is a valuable component of a healthy diet, which reduces risks of coronary disease and diabetes. Although wheat is an important cereal crop providing a substantial portion of daily calories and protein intake in the human diet, it has a low level of β-glucan. Owing to the plasticity of the polyploid wheat genome, agronomically important traits absent in the wheat primary gene pool can be introgressed from distant relatives. Barley (Hordeum vulgare L.) has a high grain β-glucan content. Earlier, we introgressed this trait into wheat in the form of whole arm compensating Robertsonian translocations (RobT) involving group-7 chromosomes of barley and all three sub-genomes of hexaploid wheat (Triticum aestivum L). In the presented research, we shortened the barley 7HL arms in these RobTs to small pericentromeric segments, using induced wheat-barley homoeologous recombination. The recombinants were selected using SNP markers and molecular cytogenetics. Plants, comprising barley cellulose synthase-like F6 gene (HvCslF6), responsible for β-glucan synthesis, had a higher grain β-glucan content than the wheat control. Three wheat-barley group-7 recombinant chromosomes involving the A, B and D sub-genomes laid the basis for a multiple-copy gene introgression to hexaploid wheat. It is hypothesized that further increases in the β-glucan content in wheat grain can be obtained by increasing the number of HvCslF6 copies through combining several recombinant chromosomes in one line. The wheat lines with four copies of HvCslF6 exceeded the β-glucan content of the lines with two copies.

Production of a complete set of wheat-barley group-7 chromosome recombinants with increased grain β-glucan content

Wheat-barley group-7 recombinant chromosomes were selected using molecular cytogenetics and SNP markers; increased grain β-glucan content was observed in wheat plants with two and four copies of HvCslF6. The soluble dietary fiber (1-3)(1-4) mixed linked β-D-glucan from cereal grains is a valuable component of a healthy diet, which reduces risks of coronary disease and diabetes. Although wheat is an important cereal crop providing a substantial portion of daily calories and protein intake in the human diet, it has a low level of β-glucan. Owing to the plasticity of the polyploid wheat genome, agronomically important traits absent in the wheat primary gene pool can be introgressed from distant relatives. Barley (Hordeum vulgare L.) has a high grain β-glucan content. Earlier, we introgressed this trait into wheat in the form of whole arm compensating Robertsonian translocations (RobT) involving group-7 chromosomes of barley and all three sub-genomes of hexaploid wheat (Triticum aestivum L). In the presented research, we shortened the barley 7HL arms in these RobTs to small pericentromeric segments, using induced wheat-barley homoeologous recombination. The recombinants were selected using SNP markers and molecular cytogenetics. Plants, comprising barley cellulose synthase-like F6 gene (HvCslF6), responsible for β-glucan synthesis, had a higher grain β-glucan content than the wheat control. Three wheat-barley group-7 recombinant chromosomes involving the A, B and D sub-genomes laid the basis for a multiple-copy gene introgression to hexaploid wheat. It is hypothesized that further increases in the β-glucan content in wheat grain can be obtained by increasing the number of HvCslF6 copies through combining several recombinant chromosomes in one line. The wheat lines with four copies of HvCslF6 exceeded the β-glucan content of the lines with two copies.

Genome-Wide Linkage Mapping Reveals Stripe Rust Resistance in Common Wheat (Triticum aestivum) Xinong1376

Stripe rust, also known as yellow rust, is a significant threat to wheat yield worldwide. Adult plant resistance (APR) is the preferred way to obtain durable protection. Chinese winter wheat cultivar Xinong1376 has maintained acceptable APR to stripe rust in field environments. To characterize APR in this cultivar, 190 F₁₀ recombinant inbred lines (RILs) developed from Xiaoyan81 × Xinong1376 were evaluated for infection type and disease severity in fields either artificially or naturally inoculated. The population along with parents were genotyped using the Illumina 90K single-nucleotide polymorphism arrays. Six quantitative trait loci (QTL) were detected using the inclusive composite interval mapping method. QYr.nwafu-4AL and QYr.nwafu-6BL.3 conferred stable resistance in all environments, and likely corresponded to a gene-rich region on the long arm of chromosomes 4A and 6B. QYr.nwafu-5AL, QYr.nwafu-5BL, QYr.nwafu-3BL.1, and QYr.nwafu-3BL.2 were detected only in some environments but enhanced the level of resistance conferred by QYr.nwafu-4AL and QYr.nwafu-6BL.3. Kompetitive allele-specific PCR (KASP) markers developed for QYr.nwafu-4AL and QYr.nwafu-6BL.3 were confirmed in a subset of RILs and 133 wheat genotypes. The QTL on 4AL and 6BL with their linked KASP markers would be useful for marker-assisted selection to improve stripe rust resistance in breeding programs.

Physical information of 2705 PCR-based molecular markers and the evaluation of their potential use in wheat

Genetic information of polymerase chain reaction (PCR)-based markers, one of the main tools of genetics and genomics research in wheat, have been well documented in wheat. However, the physical position in relation to these markers has not yet been systematically characterized. Aim of this study was to characterize the physical information of thousands of widely used molecular markers.We first assigned 2705 molecular markers to wheat physical map, of which 86.1% and 84.7% were the best hits to chromosome survey sequencing (CSS) project (CSS-contigs) and International Wheat Genome Sequencing Consortium Reference Sequence v1.0 (IWGSC RefSeq v1.0), respectively. Physical position of 96.2% markers were predicated based on BLAST analysis, were in accordance with that of the previous nullisomic/aneuploidy/linkage analysis. A suggestive high-density physical map with 4643 loci was constructed, spanning 14.01 Gb (82.4%) of the wheat genome, with 3.02 Mb between adjacent markers. Both forward and reverse primer sequences of 1166 markers had consistent best hits to IWGSC RefSeq v1.0 based on BLAST analysis, and the corresponding allele sizes were characterized. A detailed physical map with 1532 loci was released, spanning 13.93 Gb (81.9%) of the wheat genome, with 9.09 Mb between adjacent markers. Characteristic of recombination rates in different chromosomal regions was discussed. In addition, markers with multiple sites were aligned to homoeologous sites with a consistent order, confirming that a collinearity existed among A, B and D subgenomes. This study facilitates the integration of physical and genetical information of molecular markers, which could be of value for use in genetics and genomics research such as gene/QTL map-based cloning and marker-assisted selection.

A major QTL co-localized on chromosome 6BL and its epistatic interaction for enhanced wheat stripe rust resistance

Co-localization of a major QTL for wheat stripe rust resistance to a 3.9-cM interval on chromosome 6BL across both populations and another QTL on chromosome 2B with epistatic interaction. Cultivars with diverse resistance are the optimal strategy to minimize yield losses caused by wheat stripe rust (Puccinia striiformis f. sp. tritici). Two wheat populations involving resistant wheat lines P10078 and Snb"S" from CIMMYT were evaluated for stripe rust response in multiple environments. Pool analysis by Wheat660K SNP array showed that the overlapping interval on chromosome 6B likely harbored a major QTL between two populations. Then, linkage maps were constructed using KASP markers, and a co-localized locus with large effect on chromosome 6BL was detected using QTL analysis in both populations. The coincident QTL, named QYr.nwafu-6BL.2, explained 59.7% of the phenotypic maximum variation in the Mingxian 169 × P10078 and 52.5% in the Zhengmai 9023 × Snb"S" populations, respectively. This co-localization interval spanning 3.9 cM corresponds to ~ 30.5-Mb genomic region of the newest common wheat reference genome (IWGSC RefSeq v.1.0). In addition, another QTL was also detected on chromosome 2B in Zhengmai 9023 × Snb"S" population and it can accelerate expression of QYr.nwafu-6BL.2 to enhance resistance with epistatic interaction. Allowing for Pst response, marker genotypes, pedigree analysis and relative genetic distance, QYr.nwafu-6BL.2 is likely to be a distinct adult plant resistance QTL. Haplotype analysis of QYr.nwafu-6BL.2 revealed specific SNPs or alleles in the target region from a diversity panel of 176 unrelated wheat accessions. This QTL region provides opportunity for further map-based cloning, and haplotypes analysis enables pyramiding favorable alleles into commercial cultivars by marker-assisted selection.

Zinc finger nuclease-mediated precision genome editing of an endogenous gene in hexaploid bread wheat (Triticum aestivum) using a DNA repair template

Sequence-specific nucleases have been used to engineer targeted genome modifications in various plants. While targeted gene knockouts resulting in loss of function have been reported with relatively high rates of success, targeted gene editing using an exogenously supplied DNA repair template and site-specific transgene integration has been more challenging. Here, we report the first application of zinc finger nuclease (ZFN)-mediated, nonhomologous end-joining (NHEJ)-directed editing of a native gene in allohexaploid bread wheat to introduce, via a supplied DNA repair template, a specific single amino acid change into the coding sequence of acetohydroxyacid synthase (AHAS) to confer resistance to imidazolinone herbicides. We recovered edited wheat plants having the targeted amino acid modification in one or more AHAS homoalleles via direct selection for resistance to imazamox, an AHAS-inhibiting imidazolinone herbicide. Using a cotransformation strategy based on chemical selection for an exogenous marker, we achieved a 1.2% recovery rate of edited plants having the desired amino acid change and a 2.9% recovery of plants with targeted mutations at the AHAS locus resulting in a loss-of-function gene knockout. The latter results demonstrate a broadly applicable approach to introduce targeted modifications into native genes for nonselectable traits. All ZFN-mediated changes were faithfully transmitted to the next generation.

Meiotic homoeologous recombination-based mapping of wheat chromosome 2B and its homoeologues in Aegilops speltoides and Thinopyrum elongatum

We physically dissected and mapped wheat chromosome 2B and its homoeologues in Aegilops speltoides and Thinopyrum elongatum based on meiotic homoeologous recombination, providing a unique physical framework for genome studies. Common wheat has a large and complex genome with narrow genetic diversity and various degrees of recombination between the A, B, and D subgenomes. This has limited the homologous recombination-based genome studies in wheat. Here, we exploited meiotic homoeologous recombination for molecular mapping of wheat chromosome 2B and its homoeologue 2S from Aegilops speltoides and 2E from Thinopyrum elongatum. The 2B-2S and 2B-2E recombination was induced by the ph1b mutant, and recovered using molecular markers and fluorescent genomic in situ hybridization (FGISH). A total of 112 2B-2S and 87 2B-2E recombinants involving different chromosome regions were developed and physically delineated by FGISH. The 2B-2S and 2B-2E recombination hotspots mapped to the subterminal regions on both arms. Recombination hotspots with the highest recombination rates mapped to the short arms. Eighty-three 2B-2S and 67 2B-2E recombinants were genotyped using the wheat 90 K SNP arrays. Based on the genotyping results and FGISH patterns of the recombinants, chromosomes 2B, 2S, and 2E were partitioned into 93, 66, and 46 bins, respectively. In total, 1037 SNPs physically mapped onto distinct bins of these three homoeologous chromosomes. A homoeologous recombination-based bin map was constructed for chromosome 2B, providing a unique physical framework for genome studies in wheat and its relatives. Meiotic homoeologous recombination also facilitates gene introgression to diversify the wheat genome for germplasm development. Therefore, homoeologous recombination-based studies enhance understanding of the wheat genome and its homoeologous counterparts from wild grasses, and expand the genetic variability of the wheat genome.

Pm61: a recessive gene for resistance to powdery mildew in wheat landrace Xuxusanyuehuang identified by comparative genomics analysis

A single recessive powdery mildew resistance gene Pm61 from wheat landrace Xuxusanyuehuang was mapped within a 0.46-cM genetic interval spanning a 1.3-Mb interval of the genomic region of chromosome arm 4AL. Epidemics of powdery mildew incited by the biotrophic fungus Blumeria graminis f. sp. tritici (Bgt) have caused significant yield reductions in many wheat (Triticum aestivum)-producing regions. Identification of powdery mildew resistance genes is required for sustainable improvement of wheat for disease resistance. Chinese wheat landrace Xuxusanyuehuang was resistant to several Bgt isolates at the seedling stage. Genetic analysis based on the inoculation of Bgt isolate E09 on the F₁, F₂, and F_2:3 populations produced by crossing Xuxusanyuehuang to susceptible cultivar Mingxian 169 revealed that the resistance of Xuxusanyuehuang was controlled by a single recessive gene. Bulked segregant analysis and simple sequence repeat (SSR) mapping placed the gene on chromosome bin 4AL-4-0.80-1.00. Comparative genomics analysis was performed to detect the collinear genomic regions of Brachypodium distachyon, rice, sorghum, Aegilops tauschii, T. urartu, and T. turgidum ssp. dicoccoides. Based on the use of 454 contig sequences and the International Wheat Genome Sequence Consortium survey sequence of Chinese Spring wheat, four EST-SSR and seven SSR markers were linked to the gene. An F₅ recombinant inbred line population derived from Xuxusanyuehuang × Mingxian 169 cross was used to develop the genetic linkage map. The gene was localized in a 0.46-cM genetic interval between Xgwm160 and Xicsx79 corresponding to 1.3-Mb interval of the genomic region in wheat genome. This is a new locus for powdery mildew resistance on chromosome arm 4AL and is designated Pm61.

A PSTOL-like gene, TaPSTOL, controls a number of agronomically important traits in wheat

BACKGROUND

Phosphorus (P) is an essential macronutrient for plant growth, and is required in large quantities by elite varieties of crops to maintain yields. Approximately 70% of global cultivated land suffers from P deficiency, and it has recently been estimated that worldwide P resources will be exhausted by the end of this century, increasing the demand for crops more efficient in their P usage. A greater understanding of how plants are able to maintain yield with lower P inputs is, therefore, highly desirable to both breeders and farmers. Here, we clone the wheat (Triticum aestivum L.) homologue of the rice PSTOL gene (OsPSTOL), and characterize its role in phosphate nutrition plus other agronomically important traits.

RESULTS

TaPSTOL is a single copy gene located on the short arm of chromosome 5A, encoding a putative kinase protein, and shares a high level of sequence similarity to OsPSTOL. We re-sequenced TaPSTOL from 24 different wheat accessions and (3) three T. durum varieties. No sequence differences were detected in 26 of the accessions, whereas two indels were identified in the promoter region of one of the durum wheats. We characterised the expression of TaPSTOL under different P concentrations and demonstrated that the promoter was induced in root tips and hairs under P limiting conditions. Overexpression and RNAi silencing of TaPSTOL in transgenic wheat lines showed that there was a significant effect upon root biomass, flowering time independent of P treatment, tiller number and seed yield, correlating with the expression of TaPSTOL. However this did not increase PUE as elevated P concentration in the grain did not correspond to increased yields.

CONCLUSIONS

Manipulation of TaPSTOL expression in wheat shows it is responsible for many of the previously described phenotypic advantages as OsPSTOL except yield. Furthermore, we show TaPSTOL contributes to additional agronomically important traits including flowering time and grain size. Analysis of TaPSTOL sequences from a broad selection of wheat varieties, encompassing 91% of the genetic diversity in UK bread wheat, showed that there is very little genetic variation in this gene, which would suggest that this locus may have been under high selection pressure.

Genetic Analysis of Fusarium Head Blight Resistance in CIMMYT Bread Wheat Line C615 Using Traditional and Conditional QTL Mapping

Fusarium head blight (FHB) is a destructive wheat disease present throughout the world, and host resistance is an effective and economical strategy used to control FHB. Lack of adequate resistance resource is still a main bottleneck for FHB genetics and wheat breeding research. The synthetic-derived bread wheat line C615, which does not carry the Fhb1 gene, is a promising source of FHB resistance for breeding. A population of 198 recombinant inbred lines (RILs) produced by crossing C615 with the susceptible cultivar Yangmai 13 was evaluated for FHB response using point and spray inoculations. As the disease phenotype is frequently complicated by other agronomic traits, we used both traditional and multivariate conditional QTL mapping approaches to investigate the genetic relationships (at the individual QTL level) between FHB resistance and plant height (PH), spike compactness (SC), and days to flowering (FD). A linkage map was constructed from 3,901 polymorphic SNP markers, which covered 2,549.2 cM. Traditional and conditional QTL mapping analyses found 13 and 22 QTL for FHB, respectively; 10 were identified by both methods. Among these 10, three QTL from C615 were detected in multiple years; these QTL were located on chromosomes 2AL, 2DS, and 2DL. Conditional QTL mapping analysis indicated that, at the QTL level, SC strongly influenced FHB in point inoculation; whereas PH and SC contributed more to FHB than did FD in spray inoculation. The three stable QTL (QFhbs-jaas.2AL, QFhbp-jaas.2DS, and QFhbp-jaas.2DL) for FHB were partly affected by or were independent of the three agronomic traits. The QTL detected in this study improve our understanding of the genetic relationships between FHB response and related traits at the QTL level and provide useful information for marker-assisted selection for the improvement of FHB resistance in breeding.

An efficient approach for the development of genome-specific markers in allohexaploid wheat (Triticum aestivum L.) and its application in the construction of high-density linkage maps of the D genome

In common wheat, the development of genotyping platforms has been hampered by the large size of the genome, its highly repetitive elements and its allohexaploid nature. However, recent advances in sequencing technology provide opportunities to resolve these difficulties. Using next-generation sequencing and gene-targeting sequence capture, 12,551 nucleotide polymorphisms were detected in the common wheat varieties 'Hatsumochi' and 'Kitahonami' and were assigned to chromosome arms using International Wheat Genome Sequencing Consortium survey sequences. Because the number of markers for D genome chromosomes in commercially available wheat single nucleotide polymorphism arrays is insufficient, we developed markers using a genome-specific amplicon sequencing strategy. Approximately 80% of the designed primers successfully amplified D genome-specific products, suggesting that by concentrating on a specific subgenome, we were able to design successful markers as efficiently as could be done in a diploid species. The newly developed markers were uniformly distributed across the D genome and greatly extended the total coverage. Polymorphisms were surveyed in six varieties, and 31,542 polymorphic sites and 5,986 potential marker sites were detected in the D genome. The marker development and genotyping strategies are cost effective, robust and flexible and may enhance multi-sample studies in the post-genomic era in wheat.

Mining and Development of Novel SSR Markers Using Next Generation Sequencing (NGS) Data in Plants

Microsatellites, or simple sequence repeats (SSRs), are one of the most informative and multi-purpose genetic markers exploited in plant functional genomics. However, the discovery of SSRs and development using traditional methods are laborious, time-consuming, and costly. Recently, the availability of high-throughput sequencing technologies has enabled researchers to identify a substantial number of microsatellites at less cost and effort than traditional approaches. Illumina is a noteworthy transcriptome sequencing technology that is currently used in SSR marker development. Although 454 pyrosequencing datasets can be used for SSR development, this type of sequencing is no longer supported. This review aims to present an overview of the next generation sequencing, with a focus on the efficient use of de novo transcriptome sequencing (RNA-Seq) and related tools for mining and development of microsatellites in plants.

Genome-wide identification, classification, evolutionary analysis and gene expression patterns of the protein kinase gene family in wheat and Aegilops tauschii

In this study we systematically identified and classified PKs in Triticum aestivum, Triticum urartu and Aegilops tauschii. Domain distribution and exon-intron structure analyses of PKs were performed, and we found conserved exon-intron structures within the exon phases in the kinase domain. Collinearity events were determined, and we identified various T. aestivum PKs from polyploidizations and tandem duplication events. Global expression pattern analysis of T. aestivum PKs revealed that some PKs might participate in the signaling pathways of stress response and developmental processes. QRT-PCR of 15 selected PKs were performed under drought treatment and with infection of Fusarium graminearum to validate the prediction of microarray. The protein kinase (PK) gene superfamily is one of the largest families in plants and participates in various plant processes, including growth, development, and stress response. To better understand wheat PKs, we conducted genome-wide identification, classification, evolutionary analysis and expression profiles of wheat and Ae. tauschii PKs. We identified 3269, 1213 and 1448 typical PK genes in T. aestivum, T. urartu and Ae. tauschii, respectively, and classified them into major groups and subfamilies. Domain distributions and gene structures were analyzed and visualized. Some conserved intron-exon structures within the conserved kinase domain were found in T. aestivum, T. urartu and Ae. tauschii, as well as the primitive land plants Selaginella moellendorffii and Physcomitrella patens, revealing the important roles and conserved evolutionary history of these PKs. We analyzed the collinearity events of T. aestivum PKs and identified PKs from polyploidizations and tandem duplication events. Global expression pattern analysis of T. aestivum PKs revealed tissue-specific and stress-specific expression profiles, hinting that some wheat PKs may regulate abiotic and biotic stress response signaling pathways. QRT-PCR of 15 selected PKs were performed under drought treatment and with infection of F. graminearum to validate the prediction of microarray. Our results will provide the foundational information for further studies on the molecular functions of wheat PKs.

Gametocidal genes of Aegilops: segregation distorters in wheat-Aegilops wide hybridization

Aegilops is a genus belonging to the family Poaceace, which have played an indispensible role in the evolution of bread wheat and continues to do so by transferring genes by wide hybridization. Being the secondary gene pool of wheat, gene transfer from Aegilops poses difficulties and segregation distortion is common. Gametocidal genes are the most well characterized class of segregation distorters reported in interspecific crosses of wheat with Aegilops. These "selfish" genetic elements ensure their preferential transmission to progeny at the cost of gametes lacking them without providing any phenotypic benefits to the plant, thereby causing a proportional reduction in fertility. Gametocidal genes (Gc) have been reported in different species of Aegilops belonging to the sections Aegilops (Ae. geniculata and Ae. triuncialis), Cylindropyrum (Ae. caudata and Ae. cylindrica), and Sitopsis (Ae. longissima, Ae. sharonensis, and Ae. speltoides). Gametocidal activity is mostly confined to 2, 3, and 4 homeologous groups of C, S, S¹, S^sh, and M^g genomes. Removal of such genes is necessary for successful alien gene introgression and can be achieved by mutagenesis or allosyndetic pairing. However, there are some instances where Gc genes are constructively utilized for development of deletion stocks in wheat, improving genetic variability and chromosome engineering.

Evolution of Gene Expression Balance Among Homeologs of Natural Polyploids

Polyploidy is a major evolutionary process in eukaryotes, yet the expression balance of homeologs in natural polyploids is largely unknown. To study this expression balance, the expression patterns of 2180 structurally well-characterized genes of wheat were studied, of which 813 had the expected three copies and 375 had less than three. Copy numbers of the remaining 992 ranged from 4 to 14, including homeologs, orthologs, and paralogs. Of the genes with three structural copies corresponding to homeologs, 55% expressed from all three, 38% from two, and the remaining 7% expressed from only one of the three copies. Homeologs of 76-87% of the genes showed differential expression patterns in different tissues, thus have evolved different gene expression controls, possibly resulting in novel functions. Homeologs of 55% of the genes showed tissue-specific expression, with the largest percentage (14%) in the anthers and the smallest (7%) in the pistils. The highest number (1.72/3) of homeologs/gene expression was in the roots and the lowest (1.03/3) in the anthers. As the expression of homeologs changed with changes in structural copy number, about 30% of the genes showed dosage dependence. Chromosomal location also impacted expression pattern as a significantly higher proportion of genes in the proximal regions showed expression from all three copies compared to that present in the distal regions.

Wheat paleohistory created asymmetrical genomic evolution

Following the triplication reported in Brassiceae ∼10million years ago, and at the basis of rosids ∼100million years ago, bias in organization and regulation, known as subgenome dominance, has been reported between the three post-polyploidy compartments referenced to as less fractionated (LF), medium fractionated (MF1) and more fractionated (MF2), that have been proposed to derive from an hexaploidization event involving ancestors of 7-14-21 chromosomes. Modern bread wheat experienced similar paleohistory during the last half million year of evolution opening a new hypothesis where the wheat genome is at the earliest stages on the road of diploidization through subgenome dominance driving asymmetry in gene content, gene expression abundance, transposable element content as dynamics and epigenetic control between the A, B and D subgenomes.

Homoeologous recombination-based transfer and molecular cytogenetic mapping of a wheat streak mosaic virus and Triticum mosaic virus resistance gene Wsm3 from Thinopyrum intermedium to wheat

Here, we report the production of a wheat- Thinopyrum intermedium recombinant stock conferring resistance to wheat streak mosaic virus and Triticum mosaic virus. Wheat streak mosaic caused by the wheat streak mosaic virus (WSMV) is an important disease of bread wheat (Triticum aestivum) worldwide. To date, only three genes conferring resistance to WSMV have been named and two, Wsm1 and Wsm3, were derived from the distantly related wild relative Thinopyrum intermedium. Wsm3 is only available in the form of a compensating wheat-Th. intermedium whole-arm Robertsonian translocation T7BS·7S#3L. Whole-arm alien transfers usually suffer from linkage drag, which prevents their use in cultivar improvement. Here, we report ph1b-induced homoeologous recombination to shorten the Th. intermedium segment and recover a recombinant chromosome consisting of the short arm of wheat chromosome 7B, part of the long arm of 7B, and the distal 43% of the long arm derived from the Th. intermedium chromosome arm 7S#3L. The recombinant chromosome T7BS·7BL-7S#3L confers resistance to WSMV at 18 and 24 °C and also confers resistance to Triticum mosaic virus, but only at 18 °C. Wsm3 is the only gene conferring resistance to WSMV at a high temperature level of 24 °C. We also developed a user-friendly molecular marker that will allow to monitor the transfer of Wsm3 in breeding programs. Wsm3 is presently being transferred to adapted hard red winter wheat cultivars and can be used directly in wheat improvement.

A High Resolution Radiation Hybrid Map of Wheat Chromosome 4A

Bread wheat has a large and complex allohexaploid genome with low recombination level at chromosome centromeric and peri-centromeric regions. This significantly hampers ordering of markers, contigs of physical maps and sequence scaffolds and impedes obtaining of high-quality reference genome sequence. Here we report on the construction of high-density and high-resolution radiation hybrid (RH) map of chromosome 4A supported by high-density chromosome deletion map. A total of 119 endosperm-based RH lines of two RH panels and 15 chromosome deletion bin lines were genotyped with 90K iSelect single nucleotide polymorphism (SNP) array. A total of 2316 and 2695 markers were successfully mapped to the 4A RH and deletion maps, respectively. The chromosome deletion map was ordered in 19 bins and allowed precise identification of centromeric region and verification of the RH panel reliability. The 4A-specific RH map comprises 1080 mapping bins and spans 6550.9 cR with a resolution of 0.13 Mb/cR. Significantly higher mapping resolution in the centromeric region was observed as compared to recombination maps. Relatively even distribution of deletion frequency along the chromosome in the RH panel was observed and putative functional centromere was delimited within a region characterized by two SNP markers.

A High-Density Consensus Map of Common Wheat Integrating Four Mapping Populations Scanned by the 90K SNP Array

A high-density consensus map is a powerful tool for gene mapping, cloning and molecular marker-assisted selection in wheat breeding. The objective of this study was to construct a high-density, single nucleotide polymorphism (SNP)-based consensus map of common wheat (Triticum aestivum L.) by integrating genetic maps from four recombinant inbred line populations. The populations were each genotyped using the wheat 90K Infinium iSelect SNP assay. A total of 29,692 SNP markers were mapped on 21 linkage groups corresponding to 21 hexaploid wheat chromosomes, covering 2,906.86 cM, with an overall marker density of 10.21 markers/cM. Compared with the previous maps based on the wheat 90K SNP chip detected 22,736 (76.6%) of the SNPs with consistent chromosomal locations, whereas 1,974 (6.7%) showed different chromosomal locations, and 4,982 (16.8%) were newly mapped. Alignment of the present consensus map and the wheat expressed sequence tags (ESTs) Chromosome Bin Map enabled assignment of 1,221 SNP markers to specific chromosome bins and 819 ESTs were integrated into the consensus map. The marker orders of the consensus map were validated based on physical positions on the wheat genome with Spearman rank correlation coefficients ranging from 0.69 (4D) to 0.97 (1A, 4B, 5B, and 6A), and were also confirmed by comparison with genetic position on the previously 40K SNP consensus map with Spearman rank correlation coefficients ranging from 0.84 (6D) to 0.99 (6A). Chromosomal rearrangements reported previously were confirmed in the present consensus map and new putative rearrangements were identified. In addition, an integrated consensus map was developed through the combination of five published maps with ours, containing 52,607 molecular markers. The consensus map described here provided a high-density SNP marker map and a reliable order of SNPs, representing a step forward in mapping and validation of chromosomal locations of SNPs on the wheat 90K array. Moreover, it can be used as a reference for quantitative trait loci (QTL) mapping to facilitate exploitation of genes and QTL in wheat breeding.

Genome-wide linkage mapping of QTL for black point reaction in bread wheat (Triticum aestivum L.)

Nine QTL for black point resistance in wheat were identified using a RIL population derived from a Linmai 2/Zhong 892 cross and 90K SNP assay. Black point, discoloration of the embryo end of the grain, downgrades wheat grain quality leading to significant economic losses to the wheat industry. The availability of molecular markers will accelerate improvement of black point resistance in wheat breeding. The aims of this study were to identify quantitative trait loci (QTL) for black point resistance and tightly linked molecular markers, and to search for candidate genes using a high-density genetic linkage map of wheat. A recombinant inbred line (RIL) population derived from the cross Linmai 2/Zhong 892 was evaluated for black point reaction during the 2011-2012, 2012-2013 and 2013-2014 cropping seasons, providing data for seven environments. A high-density linkage map was constructed by genotyping the RILs with the wheat 90K single nucleotide polymorphism (SNP) chip. Composite interval mapping detected nine QTL on chromosomes 2AL, 2BL, 3AL, 3BL, 5AS, 6A, 7AL (2) and 7BS, designated as QBp.caas-2AL, QBp.caas-2BL, QBp.caas-3AL, QBp.caas-3BL, QBp.caas-5AS, QBp.caas-6A, QBp.caas-7AL.1, QBp.caas-7AL.2 and QBp.caas-7BS, respectively. All resistance alleles, except for QBp.caas-7AL.1 from Linmai 2, were contributed by Zhong 892. QBp.caas-3BL, QBp.caas-5AS, QBp.caas-7AL.1, QBp.caas-7AL.2 and QBp.caas-7BS probably represent new loci for black point resistance. Sequences of tightly linked SNPs were used to survey wheat and related cereal genomes identifying three candidate genes for black point resistance. The tightly linked SNP markers can be used in marker-assisted breeding in combination with the kompetitive allele specific PCR technique to improve black point resistance.

Collinearity Analysis and High-Density Genetic Mapping of the Wheat Powdery Mildew Resistance Gene Pm40 in PI 672538

The wheat powdery mildew resistance gene Pm40, which is located on chromosomal arm 7BS, is effective against nearly all prevalent races of Blumeria graminis f. sp tritici (Bgt) in China and is carried by the common wheat germplasm PI 672538. A set of the F1, F2 and F2:3 populations from the cross of the resistant PI 672538 with the susceptible line L1034 were used to conduct genetic analysis of powdery mildew resistance and construct a high-density linkage map of the Pm40 gene. We constructed a high-density linkage genetic map with a total length of 6.18 cM and average spacing between markers of 0.48 cM.Pm40 is flanked by Xwmc335 and BF291338 at genetic distances of 0.58 cM and 0.26 cM, respectively, in deletion bin C-7BS-1-0.27. Comparative genomic analysis based on EST-STS markers established a high level of collinearity of the Pm40 genomic region with a 1.09-Mbp genomic region on Brachypodium chromosome 3, a 1.16-Mbp genomic region on rice chromosome 8, and a 1.62-Mbp genomic region on sorghum chromosome 7. We further anchored the Pm40 target intervals to the wheat genome sequence. A putative linear index of 85 wheat contigs containing 97 genes on 7BS was constructed. In total, 9 genes could be considered as candidates for the resistances to powdery mildew in the target genomic regions, which encoded proteins that were involved in the plant defense and response to pathogen attack. These results will facilitate the development of new markers for map-based cloning and marker-assisted selection of Pm40 in wheat breeding programs.

Similarities and differences in the nuclear genome organization within Pooideae species revealed by comparative genomic in situ hybridization (GISH)

In this paper, we highlight the affinity between the genomes of key representatives of the Pooideae subfamily, revealed at the chromosomal level by genomic in situ hybridization (GISH). The analyses were conducted using labeled probes from each species to hybridize with chromosomes of every species used in this study based on a “round robin” rule. As a result, the whole chromosomes or chromosome regions were distinguished or variable types of signals were visualized to prove the different levels of the relationships between genomes used in this study. We observed the unexpected lack of signals in secondary constrictions of rye (RR) chromosomes probed by triticale (AABBRR) genomic DNA. We have also identified unlabeled chromosome regions, which point to species-specific sequences connected with disparate pathways of chromosome differentiation. Our results revealed a conservative character of coding sequence of 35S rDNA among selected species of the genera Aegilops, Brachypodium, Festuca, Hordeum, Lolium, Secale, and Triticum. In summary, we showed strong relationships in genomic DNA sequences between species which have been previously reported to be phylogenetically distant.

Molecular Characterization of Three GIBBERELLIN-INSENSITIVE DWARF2 Homologous Genes in Common Wheat

F-box protein is a core component of the ubiquitin E3 ligase SCF complex and is involved in the gibberellin (GA) signaling pathway. To elucidate the molecular mechanism of GA signaling in wheat, three homologous GIBBERELLIN-INSENSITIVE DWARF2 genes, TaGID2s, were isolated from the Chinese Spring wheat variety. A subcellular localization assay in onion epidermal cells and Arabidopsis mesophyll protoplasts showed that TaGID2s are localized in the nuclei. The expression profiles using quantitative real-time polymerase chain reaction showed that TaGID2s were downregulated by GA3. The interaction between TaGID2s and TSK1 (homologous to ASK1) in yeast indicated that TaGID2s might function as a component of an E3 ubiquitin-ligase SCF complex. Yeast two-hybrid assays showed that a GA-independent interaction occurred between three TaGID2s and RHT-A1a, RHT-B1a, and RHT-D1a. Furthermore, TaGID2s interact with most RHT-1s, such as RHT-B1h, RHT-B1i, RHT-D1e, RHT-D1f, etc., but cannot interact with RHT-B1b or RHT-B1e, which have a stop codon in the DELLA motif, resulting in a lack of a GRAS domain. In addition, RHT-B1k has a frame-shift mutation in the VHIID motif leading to loss of the LHRII motif in the GRAS domain and RHT-D1h has a missense mutation in the LHRII motif. These results indicate that TaGID2s, novel positive regulators of the GA response, recognize RHT-1s in the LHRII motif resulting in poly-ubiquitination and degradation of the DELLA protein.

Physical mapping of Agropyron cristatum chromosome 6P using deletion lines in common wheat background

Genetically stable deletion lines of Agropyron cristatum chromosome 6P in common wheat background were generated, which allowed for physical mapping of 255 6P-specific STS markers and leaf rust resistance gene(s). Chromosomal deletion lines are valuable tools for gene discovery and localization. The chromosome 6P of Agropyron cristatum (2n = 4x = 28, PPPP) confers many desirable agronomic traits to common wheat, such as higher grain number per spike, multiple fertile tiller number, and enhanced resistance to certain diseases. Although many elite genes from A. cristatum have been identified, their chromosomal locations were largely undetermined due to the lack of A. cristatum 6P deletion lines. In this study, various A. cristatum 6P deletion lines were developed using a wheat-A. cristatum 6P disomic addition line 4844-12 subjected to (60)Co-γ irradiation as well as an Aegilops cylindrica gametocidal chromosome. Twenty-six genetically stable A. cristatum 6P deletion lines in the genetic background of common wheat were obtained, and their genetic constitutions were elucidated by genomic in situ hybridization (GISH) and sequence-tagged site (STS) markers specific to A. cristatum chromosome 6P. Moreover, 255 novel chromosome 6P-specific STS markers were physically mapped to 14 regions of chromosome 6P. Field evaluation of leaf rust resistance of various deletion lines and BC1F2 populations indicated that the A.cristatum chromosome 6P-originated leaf rust resistance gene(s) was located in the region 6PS-0.81-1.00. This study will provide not only useful tools for characterization and utilization of wheat materials with alien chromosomal segments, but also novel wheat germplasms potentially valuable in wheat breeding and improvement.

A whole-genome, radiation hybrid mapping resource of hexaploid wheat

Generating a contiguous, ordered reference sequence of a complex genome such as hexaploid wheat (2n = 6x = 42; approximately 17 GB) is a challenging task due to its large, highly repetitive, and allopolyploid genome. In wheat, ordering of whole-genome or hierarchical shotgun sequencing contigs is primarily based on recombination and comparative genomics-based approaches. However, comparative genomics approaches are limited to syntenic inference and recombination is suppressed within the pericentromeric regions of wheat chromosomes, thus, precise ordering of physical maps and sequenced contigs across the whole-genome using these approaches is nearly impossible. We developed a whole-genome radiation hybrid (WGRH) resource and tested it by genotyping a set of 115 randomly selected lines on a high-density single nucleotide polymorphism (SNP) array. At the whole-genome level, 26 299 SNP markers were mapped on the RH panel and provided an average mapping resolution of approximately 248 Kb/cR1500 with a total map length of 6866 cR1500 . The 7296 unique mapping bins provided a five- to eight-fold higher resolution than genetic maps used in similar studies. Most strikingly, the RH map had uniform bin resolution across the entire chromosome(s), including pericentromeric regions. Our research provides a valuable and low-cost resource for anchoring and ordering sequenced BAC and next generation sequencing (NGS) contigs. The WGRH developed for reference wheat line Chinese Spring (CS-WGRH), will be useful for anchoring and ordering sequenced BAC and NGS based contigs for assembling a high-quality, reference sequence of hexaploid wheat. Additionally, this study provides an excellent model for developing similar resources for other polyploid species.

Cloning and Characterization of TaTGW-7A Gene Associated with Grain Weight in Wheat via SLAF-seq-BSA

Thousand-grain weight (TGW) of wheat (Triticum aestivum L.) contributes significantly to grain yield. In the present study, a candidate gene associated with TGW was identified through specific-locus amplified fragment sequencing (SLAF-seq) of DNA bulks of recombinant inbred lines (RIL) derived from the cross between Jing 411 and Hongmangchun 21. The gene was located on chromosome 7A, designated as TaTGW-7A with a complete genome sequence and an open reading frame (ORF). A single nucleotide polymorphism (SNP) was present in the first exon between two alleles at TaTGW-7A locus, resulting in a Val to Ala substitution, corresponding to a change from higher to lower TGW. Cleaved amplified polymorphic sequence (CAPS) (TGW7A) and InDel (TG9) markers were developed to discriminate the two alleles TaTGW-7Aa and TaTGW-7Ab for higher and lower TGW, respectively. A major QTL co-segregating with TaTGW-7A explained 21.7-27.1% of phenotypic variance for TGW in the RIL population across five environments. The association of TaTGW-7A with TGW was further validated in a natural population and Chinese mini-core collections. Quantitative real-time PCR revealed higher transcript levels of TaTGW-7Aa than those of TaTGW-7Ab during grain development. High frequencies of the superior allele TaTGW-7Aa for higher TGW in Chinese mini-core collections (65.0%) and 501 wheat varieties (86.0%) indicated a strong and positive selection of this allele in wheat breeding. The molecular markers TGW7A and TG9 can be used for improvement of TGW in breeding programs.

Sequencing-based high throughput mutation detection in bread wheat

Background

Forward genetic approaches have limited use for agronomic traits that can’t be reliably scored on a single plant basis. Thus, mutants in wheat and other crops are more useful for gene function studies by reverse genetic approach. With a long-term goal to develop a sequence-based mutation detection resource in hexaploid wheat, we conducted a feasibility study to accurately differentiate induced mutations from the homoeologs’ sequence variations present among the three wheat genomes.

Results

A reduced representation ApeKI library consisting of 21 Ethylmethane Sulfonate (EMS) induced mutants and two wild type cv. Indian plants was developed using individual barcode adapters and sequenced. A novel bioinformatics pipeline was developed to identify sequence variants using 178,464 wheat unigenes as a reference wheat transcriptome. In total, 14,130 mutational changes [Single Nucleotide Polymorphisms (SNPs) and Insertions/Deletions (INDELs)] and 150,511 homoeologous sequence changes were detected. On an average, 662 SNPs (ranging from 46 to 1,330) and 10 small INDELs (ranging from 0 to 23) were identified for each of the mutants. A mutation frequency of one per 5 Kb was observed with 70 % being transitions and 30 % transversions. The pipeline was tested using the known sequence changes in the three wheat genes. Genes present in the distal regions of the chromosomes were found to be more prone to EMS compared to genes present in the proximal regions. Redefined parameters identified a total of 28,348 mutational changes (1,349/plant).

Conclusions

We conclude that sequencing based mutation detection is a valuable method to identify induced mutations at large.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2112-1) contains supplementary material, which is available to authorized users.

Physical Mapping of Bread Wheat Chromosome 5A: An Integrated Approach

The huge size, redundancy, and highly repetitive nature of the bread wheat [Triticum aestivum (L.)] genome, makes it among the most difficult species to be sequenced. To overcome these limitations, a strategy based on the separation of individual chromosomes or chromosome arms and the subsequent production of physical maps was established within the frame of the International Wheat Genome Sequence Consortium (IWGSC). A total of 95,812 bacterial artificial chromosome (BAC) clones of short-arm chromosome 5A (5AS) and long-arm chromosome 5A (5AL) arm-specific BAC libraries were fingerprinted and assembled into contigs by complementary analytical approaches based on the FingerPrinted Contig (FPC) and Linear Topological Contig (LTC) tools. Combined anchoring approaches based on polymerase chain reaction (PCR) marker screening, microarray, and sequence homology searches applied to several genomic tools (i.e., genetic maps, deletion bin map, neighbor maps, BAC end sequences (BESs), genome zipper, and chromosome survey sequences) allowed the development of a high-quality physical map with an anchored physical coverage of 75% for 5AS and 53% for 5AL with high portions (64 and 48%, respectively) of contigs ordered along the chromosome. In the genome of grasses, Brachypodium [Brachypodium distachyon (L.) Beauv.], rice (Oryza sativa L.), and sorghum [Sorghum bicolor (L.) Moench] homologs of genes on wheat chromosome 5A were separated into syntenic blocks on different chromosomes as a result of translocations and inversions during evolution. The physical map presented represents an essential resource for fine genetic mapping and map-based cloning of agronomically relevant traits and a reference for the 5A sequencing projects.

Structure and Stability of Telocentric Chromosomes in Wheat

In most eukaryotes, centromeres assemble at a single location per chromosome. Naturally occurring telocentric chromosomes (telosomes) with a terminal centromere are rare but do exist. Telosomes arise through misdivision of centromeres in normal chromosomes, and their cytological stability depends on the structure of their kinetochores. The instability of telosomes may be attributed to the relative centromere size and the degree of completeness of their kinetochore. Here we test this hypothesis by analyzing the cytogenetic structure of wheat telosomes. We used a population of 80 telosomes arising from the misdivision of the 21 chromosomes of wheat that have shown stable inheritance over many generations. We analyzed centromere size by probing with the centromere-specific histone H3 variant, CENH3. Comparing the signal intensity for CENH3 between the intact chromosome and derived telosomes showed that the telosomes had approximately half the signal intensity compared to that of normal chromosomes. Immunofluorescence of CENH3 in a wheat stock with 28 telosomes revealed that none of the telosomes received a complete CENH3 domain. Some of the telosomes lacked centromere specific retrotransposons of wheat in the CENH3 domain, indicating that the stability of telosomes depends on the presence of CENH3 chromatin and not on the presence of CRW repeats. In addition to providing evidence for centromere shift, we also observed chromosomal aberrations including inversions and deletions in the short arm telosomes of double ditelosomic 1D and 6D stocks. The role of centromere-flanking, pericentromeric heterochromatin in mitosis is discussed with respect to genome/chromosome integrity.

Development of a D genome specific marker resource for diploid and hexaploid wheat

BACKGROUND

Mapping and map-based cloning of genes that control agriculturally and economically important traits remain great challenges for plants with complex highly repetitive genomes such as those within the grass tribe, Triticeae. Mapping limitations in the Triticeae are primarily due to low frequencies of polymorphic gene markers and poor genetic recombination in certain genetic regions. Although the abundance of repetitive sequence may pose common problems in genome analysis and sequence assembly of large and complex genomes, they provide repeat junction markers with random and unbiased distribution throughout chromosomes. Hence, development of a high-throughput mapping technology that combine both gene-based and repeat junction-based markers is needed to generate maps that have better coverage of the entire genome.

RESULTS

In this study, the available genomics resource of the diploid Aegilop tauschii, the D genome donor of bread wheat, were used to develop genome specific markers that can be applied for mapping in modern hexaploid wheat. A NimbleGen array containing both gene-based and repeat junction probe sequences derived from Ae. tauschii was developed and used to map the Chinese Spring nullisomic-tetrasomic lines and deletion bin lines of the D genome chromosomes. Based on these mapping data, we have now anchored 5,171 repeat junction probes and 10,892 gene probes, corresponding to 5,070 gene markers, to the delineated deletion bins of the D genome. The order of the gene-based markers within the deletion bins of the Chinese Spring can be inferred based on their positions on the Ae. tauschii genetic map. Analysis of the probe sequences against the Chinese Spring chromosome sequence assembly database facilitated mapping of the NimbleGen probes to the sequence contigs and allowed assignment or ordering of these sequence contigs within the deletion bins. The accumulated length of anchored sequence contigs is about 155 Mb, representing ~ 3.2 % of the D genome. A specific database was developed to allow user to search or BLAST against the probe sequence information and to directly download PCR primers for mapping specific genetic loci.

CONCLUSIONS

In bread wheat, aneuploid stocks have been extensively used to assign markers linked with genes/traits to chromosomes, chromosome arms, and their specific bins. Through this study, we added thousands of markers to the existing wheat chromosome bin map, representing a significant step forward in providing a resource to navigate the wheat genome. The database website ( http://probes.pw.usda.gov/ATRJM/ ) provides easy access and efficient utilization of the data. The resources developed herein can aid map-based cloning of traits of interest and the sequencing of the D genome of hexaploid wheat.

Marker development, saturation mapping, and high-resolution mapping of the Septoria nodorum blotch susceptibility gene Snn3-B1 in wheat

Septoria nodorum blotch (SNB), caused by Parastagonospora nodorum, is a severe foliar and glume disease on durum and common wheat. Pathogen-produced necrotrophic effectors (NEs) are the major determinants for SNB on leaves. One such NE is SnTox3, which evokes programmed cell death and leads to disease when recognized by the wheat Snn3-B1 gene. Here, we developed saturated genetic linkage maps of the Snn3-B1 region using two F2 populations derived from the SnTox3-sensitive line Sumai 3 crossed with different SnTox3-insensitive lines. Markers were identified and/or developed from various resources including previously mapped simple sequence repeats, bin-mapped expressed sequence tags, single nucleotide polymorphisms, and whole genome survey sequences. Subsequent high-resolution mapping of the Snn3-B1 locus in 5600 gametes delineated the gene to a 1.5 cM interval. Analysis of micro-colinearity of the Snn3-B1 region indicated that it was highly disrupted compared to rice and Brachypodium distachyon. The screening of a collection of durum and common wheat cultivars with tightly linked markers indicated they are not diagnostic for the presence of Snn3-B1, but can be useful for marker-assisted selection if the SnTox3 reactions of lines are first determined. Finally, we developed an ethyl methanesulfonate-induced mutant population of Sumai 3 where the screening of 408 M2 families led to the identification of 17 SnTox3-insensitive mutants. These mutants along with the markers and high-resolution map developed in this research provide a strong foundation for the map-based cloning of Snn3-B1, which will broaden our understanding of the wheat-P. nodorum system and plant-necrotrophic pathogen interactions in general.

A high-density, SNP-based consensus map of tetraploid wheat as a bridge to integrate durum and bread wheat genomics and breeding

Consensus linkage maps are important tools in crop genomics. We have assembled a high-density tetraploid wheat consensus map by integrating 13 data sets from independent biparental populations involving durum wheat cultivars (Triticum turgidum ssp. durum), cultivated emmer (T. turgidum ssp. dicoccum) and their ancestor (wild emmer, T. turgidum ssp. dicoccoides). The consensus map harboured 30 144 markers (including 26 626 SNPs and 791 SSRs) half of which were present in at least two component maps. The final map spanned 2631 cM of all 14 durum wheat chromosomes and, differently from the individual component maps, all markers fell within the 14 linkage groups. Marker density per genetic distance unit peaked at centromeric regions, likely due to a combination of low recombination rate in the centromeric regions and even gene distribution along the chromosomes. Comparisons with bread wheat indicated fewer regions with recombination suppression, making this consensus map valuable for mapping in the A and B genomes of both durum and bread wheat. Sequence similarity analysis allowed us to relate mapped gene-derived SNPs to chromosome-specific transcripts. Dense patterns of homeologous relationships have been established between the A- and B-genome maps and between nonsyntenic homeologous chromosome regions as well, the latter tracing to ancient translocation events. The gene-based homeologous relationships are valuable to infer the map location of homeologs of target loci/QTLs. Because most SNP and SSR markers were previously mapped in bread wheat, this consensus map will facilitate a more effective integration and exploitation of genes and QTL for wheat breeding purposes.

Chromosomal Localization and Contribution of Three Homoeologous Genes to Biosynthesis of Cytosolic Aspartate Aminotransferase in Common Wheat

Chromosomal localization of the three homoeologous genes encoding cytosolic aspartate aminotransferase in common wheat (Triticum aestivum cv. Chinese Spring, 2n = 6x = 42, AABBDD) was specified to: 3AL (0.42÷0.61), 3BL (0.38÷0.41) and 3DL (0.23÷0.81) by a comparative zymographic analysis of the enzymatic activities in deletion lines. It was also attempted to precisely explain the nature of the relationship between a number of genes encoding α and β subunits and a distribution of staining intensity of cytosolic aspartate aminotransferase allozyme activity bands using aneuploid lines of common wheat with modified third pair of homoeologous chromosomes from genomes A, B and D, on which the genes encoding subunit α (genome A) and β (genome B and D) are localized. The highest consistency between the experimental results and the theoretical distributions was achieved by substituting values of α = 0.57 and β = 0.43 in a theoretical model. These results demonstrate that the individual participation of the diploid genome A in the biosynthesis of the cytosolic aspartate aminotransferase allozymes subunits is greater than the individual participation of the diploid genomes B and D.

Identification and Characterization of the SnTox6-Snn6 Interaction in the Parastagonospora nodorum-Wheat Pathosystem

Parastagonospora nodorum is a necrotrophic fungal pathogen that causes Septoria nodorum blotch (SNB) (formerly Stagonospora nodorum blotch) on wheat. P. nodorum produces necrotrophic effectors (NE) that are recognized by dominant host sensitivity gene products resulting in disease development. The NE-host interaction is critical to inducing NE-triggered susceptibility (NETS). To date, seven NE-host sensitivity gene interactions, following an inverse gene-for-gene model, have been identified in the P. nodorum-wheat pathosystem. Here, we used a wheat mapping population that segregated for sensitivity to two previously characterized interactions (SnTox1-Snn1 and SnTox3-Snn3-B1) to identify and characterize a new interaction involving the NE designated SnTox6 and the host sensitivity gene designated Snn6. SnTox6 is a small secreted protein that induces necrosis on wheat lines harboring Snn6. Sensitivity to SnTox6, conferred by Snn6, was light-dependent and was shown to underlie a major disease susceptibility quantitative trait locus (QTL). No other QTL were identified, even though the P. nodorum isolate used in this study harbored both the SnTox1 and SnTox3 genes. Reverse transcription-polymerase chain reaction showed that the expression of SnTox1 was not detectable, whereas SnTox3 was expressed and, yet, did not play a significant role in disease development. This work expands our knowledge of the wheat-P. nodorum interaction and further establishes this system as a model for necrotrophic specialist pathosystems.

High-throughput physical map anchoring via BAC-pool sequencing

BACKGROUND

Physical maps created from large insert DNA libraries, typically cloned in BAC vector, are valuable resources for map-based cloning and de novo genome sequencing. The maps are most useful if contigs of overlapping DNA clones are anchored to chromosome(s), and ordered along them using molecular markers. Here we present a novel approach for anchoring physical maps, based on sequencing three-dimensional pools of BAC clones from minimum tilling path.

RESULTS

We used physical map of wheat chromosome arm 3DS to validate the method with two different DNA sequence datasets. The first comprised 567 genes ordered along the chromosome arm based on syntenic relationship of wheat with the sequenced genomes of Brachypodium, rice and sorghum. The second dataset consisted of 7,136 SNP-containing sequences, which were mapped genetically in Aegilops tauschii, the donor of the wheat D genome. Mapping of sequence reads from individual BAC pools to the first and the second datasets enabled unambiguous anchoring 447 and 311 3DS-specific sequences, respectively, or 758 in total.

CONCLUSIONS

We demonstrate the utility of the novel approach for BAC contig anchoring based on mass parallel sequencing of three-dimensional pools prepared from minimum tilling path of physical map. The existing genetic markers as well as any other DNA sequence could be mapped to BAC clones in a single in silico experiment. The approach reduces significantly the cost and time needed for anchoring and is applicable to any genomic project involving the construction of anchored physical map.

Fine mapping of shattering locus Br2 reveals a putative chromosomal inversion polymorphism between the two lineages of Aegilops tauschii

This work laid the foundation for cloning of shattering gene Br2 and provided first line of evidence that two major Aegilops tauschii lineages are differentiated by an inversion polymorphism. Chromosome inversions often accompany population differentiation and capture local adaptation during speciation. Aegilops tauschii, the D-genome donor species of hexaploid wheat, consists of two genetically isolated lineages, L1 and L2, but little is known about the genetic mechanisms underlying the population differentiation in this diploid species. During fine mapping of the shattering gene Br2 using a large F2 population derived from a cross between TA1604 (an L1 accession) and AL8/78 (an L2 accession), we found contrasting patterns of crossover distribution in the Br2 interval and neighboring regions despite the high local gene synteny with Brachypodium distachyon and rice. Br2 was localized in a 0.08-cM interval, and 13 marker loci formed a block, where single-crossovers were completely suppressed, but double-crossovers were enriched with a recombination rate of ~11 cM/Mb. In contrast, in a neighboring region no double-crossover was recovered, but single-crossover rate reached 24 cM/Mb, which is much higher than the genome-wide average. This result suggests a putative inversion polymorphism between the parental lines in the Br2 region. Genotyping using the markers from the Br2 region divided a collection of 55 randomly sampled A. tauschii accessions into two major groups, and they are largely genetically isolated. The two groups correspond to the L1 and L2 lineages based on their geographic distribution patterns. This provides first evidence that inversions may underlie the evolution of A. tauschii lineages. The presence of inter-lineage inversions may complicate map-based cloning in A. tauschii and transfer of useful traits to wheat.

Identification and genetic mapping of the putative Thinopyrum intermedium-derived dominant powdery mildew resistance gene PmL962 on wheat chromosome arm 2BS

Powdery resistance putatively derived from Thinopyrum intermedium in the wheat line L962 is controlled by a single dominant gene designated PmL962 and mapped to chromosome arm 2BS. Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is a destructive disease affecting the production of wheat (Triticum aestivum). Powdery mildew resistance was putatively transferred from Thinopyrum intermedium to the common wheat line L962, which conferred resistance to multiple Chinese Bgt isolates. Genetic analysis of the powdery mildew response was conducted by crossing the resistant line L962 with the susceptible line L983. Disease assessments of the F1, F2, and F2:3 populations from the cross L983/L962 indicated that resistance was controlled by a single dominant gene. A total of 373 F2:3 lines and 781 pairs of genomic simple sequence repeat (SSR) primers were employed to determine the chromosomal location of the resistance gene. The gene was linked to four publicly available and recently developed wheat genomic SSR markers and seven EST-STS markers. The resistance gene was mapped to chromosome arm 2BS based on the locations of the linked markers. Pedigree, molecular marker and resistance response data indicated that the powdery mildew resistance gene in L962 is novel. It was temporarily designated PmL962. It is flanked by Xwmc314 and BE443737at genetic distances of 2.09 and 3.74 cM, respectively, and located in a 20.77 cM interval that is co-linear with a 269.4 kb genomic region on chromosome 5 in Brachypodium distachyon and a 223.5 kb genomic region on rice (Oryza sativa) chromosome 4. The markers that are closely linked to this gene have potential applications in marker-assisted breeding.

A haplotype map of allohexaploid wheat reveals distinct patterns of selection on homoeologous genomes

Background

Bread wheat is an allopolyploid species with a large, highly repetitive genome. To investigate the impact of selection on variants distributed among homoeologous wheat genomes and to build a foundation for understanding genotype-phenotype relationships, we performed population-scale re-sequencing of a diverse panel of wheat lines.

Results

A sample of 62 diverse lines was re-sequenced using the whole exome capture and genotyping-by-sequencing approaches. We describe the allele frequency, functional significance, and chromosomal distribution of 1.57 million single nucleotide polymorphisms and 161,719 small indels. Our results suggest that duplicated homoeologous genes are under purifying selection. We find contrasting patterns of variation and inter-variant associations among wheat genomes; this, in addition to demographic factors, could be explained by differences in the effect of directional selection on duplicated homoeologs. Only a small fraction of the homoeologous regions harboring selected variants overlapped among the wheat genomes in any given wheat line. These selected regions are enriched for loci associated with agronomic traits detected in genome-wide association studies.

Conclusions

Evidence suggests that directional selection in allopolyploids rarely acted on multiple parallel advantageous mutations across homoeologous regions, likely indicating that a fitness benefit could be obtained by a mutation at any one of the homoeologs. Additional advantageous variants in other homoelogs probably either contributed little benefit, or were unavailable in populations subjected to directional selection. We hypothesize that allopolyploidy may have increased the likelihood of beneficial allele recovery by broadening the set of possible selection targets.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0606-4) contains supplementary material, which is available to authorized users.