Comparative DNA sequence analysis of wheat and rice genomes

Are cereal grasses a single genetic system?

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

Identification and characterization of RuvBL DNA helicase genes for tolerance against abiotic stresses in bread wheat (Triticum aestivum L.) and related species

Recombination UVB (sensitivity) like (RuvBL) helicase genes represent a conserved family of genes, which are known to be involved in providing tolerance against abiotic stresses like heat and drought. We identified nine wheat RuvBL genes, one each on nine different chromosomes, belonging to homoeologous groups 2, 3, and 4. The lengths of genes ranged from 1647 to 2197 bp and exhibited synteny with corresponding genes in related species including Ae. tauschii, Z. mays, O. sativa, H. vulgare, and B. distachyon. The gene sequences were associated with regulatory cis-elements and transposable elements. Two genes, namely TaRuvBL1a-4A and TaRuvBL1a-4B, also carried targets for a widely known miRNA, tae-miR164. Gene ontology revealed that these genes were closely associated with ATP-dependent formation of histone acetyltransferase complex. Analysis of the structure and function of RuvBL proteins revealed that the proteins were localized mainly in the cytoplasm. A representative gene, namely TaRuvBL1a-4A, was also shown to be involved in protein-protein interactions with ten other proteins. On the basis of phylogeny, RuvBL proteins were placed in two sub-divisions, namely RuvBL1 and RuvBL2, which were further classified into clusters and sub-clusters. In silico studies suggested that these genes were differentially expressed under heat/drought. The qRT-PCR analysis confirmed that expression of TaRuvBL genes differed among wheat cultivars, which differed in the level of thermotolerance. The present study advances our understanding of the biological role of wheat RuvBL genes and should help in planning future studies on RuvBL genes in wheat including use of RuvBL genes in breeding thermotolerant wheat cultivars.

Optimizing genomic selection of agricultural traits using K-wheat core collection

The agricultural traits that constitute basic plant breeding information are usually quantitative or complex in nature. This quantitative and complex combination of traits complicates the process of selection in breeding. This study examined the potential of genome-wide association studies (GWAS) and genomewide selection (GS) for breeding ten agricultural traits by using genome-wide SNPs. As a first step, a trait-associated candidate marker was identified by GWAS using a genetically diverse 567 Korean (K)-wheat core collection. The accessions were genotyped using an Axiom^® 35K wheat DNA chip, and ten agricultural traits were determined (awn color, awn length, culm color, culm length, ear color, ear length, days to heading, days to maturity, leaf length, and leaf width). It is essential to sustain global wheat production by utilizing accessions in wheat breeding. Among the traits associated with awn color and ear color that showed a high positive correlation, a SNP located on chr1B was significantly associated with both traits. Next, GS evaluated the prediction accuracy using six predictive models (G-BLUP, LASSO, BayseA, reproducing kernel Hilbert space, support vector machine (SVM), and random forest) and various training populations (TPs). With the exception of the SVM, all statistical models demonstrated a prediction accuracy of 0.4 or better. For the optimization of the TP, the number of TPs was randomly selected (10%, 30%, 50% and 70%) or divided into three subgroups (CC-sub 1, CC-sub 2 and CC-sub 3) based on the subpopulation structure. Based on subgroup-based TPs, better prediction accuracy was found for awn color, culm color, culm length, ear color, ear length, and leaf width. A variety of Korean wheat cultivars were used for validation to evaluate the prediction ability of populations. Seven out of ten cultivars showed phenotype-consistent results based on genomics-evaluated breeding values (GEBVs) calculated by the reproducing kernel Hilbert space (RKHS) predictive model. Our research provides a basis for improving complex traits in wheat breeding programs through genomics assisted breeding. The results of our research can be used as a basis for improving wheat breeding programs by using genomics-assisted breeding.

Breeding for Higher Yields of Wheat and Rice through Modifying Nitrogen Metabolism

Wheat and rice produce nutritious grains that provide 32% of the protein in the human diet globally. Here, we examine how genetic modifications to improve assimilation of the inorganic nitrogen forms ammonium and nitrate into protein influence grain yield of these crops. Successful breeding for modified nitrogen metabolism has focused on genes that coordinate nitrogen and carbon metabolism, including those that regulate tillering, heading date, and ammonium assimilation. Gaps in our current understanding include (1) species differences among candidate genes in nitrogen metabolism pathways, (2) the extent to which relative abundance of these nitrogen forms across natural soil environments shape crop responses, and (3) natural variation and genetic architecture of nitrogen-mediated yield improvement. Despite extensive research on the genetics of nitrogen metabolism since the rise of synthetic fertilizers, only a few projects targeting nitrogen pathways have resulted in development of cultivars with higher yields. To continue improving grain yield and quality, breeding strategies need to focus concurrently on both carbon and nitrogen assimilation and consider manipulating genes with smaller effects or that underlie regulatory networks as well as genes directly associated with nitrogen metabolism.

Unravelling consensus genomic regions associated with quality traits in wheat using meta-analysis of quantitative trait loci

Meta-analysis in wheat for three major quality traits identified 110 meta-QTL (MQTL) with reduced confidence interval (CI). Five GWAS validated MQTL (viz., 1A.1, 1B.2, 3B.4, 5B.2, and 6B.2), each involving more than 20 initial QTL and reduced CI (95%) (< 2 cM), were selected for quality breeding programmes. Functional characterization including candidate gene mining and expression analysis discovered 44 high confidence candidate genes associated with quality traits. A meta-analysis of quantitative trait loci (QTL) associated with dough rheology properties, nutritional traits, and processing quality traits was conducted in wheat. For this purpose, as many as 2458 QTL were collected from 50 interval mapping studies published during 2013-2020. Of the total QTL, 1126 QTL were projected onto the consensus map saturated with 249,603 markers which led to the identification of 110 meta-QTL (MQTL). These MQTL exhibited an 18.84-fold reduction in the average CI compared to the average CI of the initial QTL (ranging from 14.87 to 95.55 cM with an average of 40.35 cM). Of the 110, 108 MQTL were physically anchored to the wheat reference genome, including 51 MQTL verified with marker-trait associations (MTAs) reported from earlier genome-wide association studies. Candidate gene (CG) mining allowed the identification of 2533 unique gene models from the MQTL regions. In-silico expression analysis discovered 439 differentially expressed gene models with > 2 transcripts per million expressions in grains and related tissues, which also included 44 high-confidence CGs involved in the various cellular and biochemical processes related to quality traits. Nine functionally characterized wheat genes associated with grain protein content, high-molecular-weight glutenin, and starch synthase enzymes were also found to be co-localized with some of the MQTL. Synteny analysis between wheat and rice MQTL regions identified 23 wheat MQTL syntenic to 16 rice MQTL associated with quality traits. Furthermore, 64 wheat orthologues of 30 known rice genes were detected in 44 MQTL regions. Markers flanking the MQTL identified in the present study can be used for marker-assisted breeding and as fixed effects in the genomic selection models for improving the prediction accuracy during quality breeding. Wheat orthologues of rice genes and other CGs available from MQTLs can be promising targets for further functional validation and to better understand the molecular mechanism underlying the quality traits in wheat.

The Potential Role of RNA Structure in Crop Molecular Breeding

The continually growing human population creates a concomitantly increasing demand for nutritious crops with high yields. Advances in high throughput sequencing technologies have revealed the genetic architecture of major crops. This includes extensive information enabling comprehensive genetic markers for breeding selection, new gene discoveries, and novel gene regulatory strategies for crop editing. RNA structure is an important type of genetic feature, essential for post-transcriptional regulation of gene expression. Here, we summarize recent advances in genome-wide RNA structure studies in crops and review the associated RNA structure-mediated regulation of gene expression. We also discuss the functional importance of those single nucleotide variations that induce large RNA structure disparities. Lastly, we discuss the potential role of RNA structure in crop molecular breeding.

Comparative Study of Starch Phosphorylase Genes and Encoded Proteins in Various Monocots and Dicots with Emphasis on Maize

Starch phosphorylase (PHO) is a multimeric enzyme with two distinct isoforms: plastidial starch phosphorylase (PHO1) and cytosolic starch phosphorylase (PHO2). PHO1 specifically resides in the plastid, while PHO2 is found in the cytosol. Both play a critical role in the synthesis and degradation of starch. This study aimed to report the detailed structure, function, and evolution of genes encoding PHO1 and PHO2 and their protein ligand-binding sites in eight monocots and four dicots. "True" orthologs of PHO1 and PHO2 of Oryza sativa were identified, and the structure of the enzyme at the protein level was studied. The genes controlling PHO2 were found to be more conserved than those controlling PHO1; the variations were mainly due to the variable sequence and length of introns. Cis-regulatory elements in the promoter region of both genes were identified, and the expression pattern was analyzed. The real-time quantitative polymerase chain reaction indicated that PHO2 was expressed in all tissues with a uniform pattern of transcripts, and the expression pattern of PHO1 indicates that it probably contributes to the starch biosynthesis during seed development in Zea mays. Under abscisic acid (ABA) treatment, PHO1 was found to be downregulated in Arabidopsis and Hordeum vulgare. However, we found that ABA could up-regulate the expression of both PHO1 and PHO2 within 12 h in Zea mays. In all monocots and dicots, the 3D structures were highly similar, and the ligand-binding sites were common yet fluctuating in the position of aa residues.

Meta-QTLs, ortho-MQTLs and candidate genes for the traits contributing to salinity stress tolerance in common wheat (Triticum aestivum L.)

A meta-analysis of QTLs associated with the traits contributing to salinity tolerance was undertaken in wheat to detect consensus and robust meta-QTLs (MQTLs) using 844 known QTLs retrieved from 26 earlier studies. A consensus map with a total length of 4621.56 cM including 7710 markers was constructed using 21 individual linkage maps and three previously published integrated genetic maps. Out of 844 QTLs, 571 QTLs were projected on the consensus map which gave origin to 100 MQTLs. Interestingly, 49 MQTLs were co-located with marker-trait associations reported in wheat genome-wide association studies for the traits contributing to salinity stress tolerance. Five potential MQTLs associated with the major salinity-responsive traits were also identified to be utilized in the breeding programme. In the resulted MQTLs, the average confidence interval (CI, 3.58 cM) was reduced up to 4.16 folds compared to the mean CI of the initial QTLs. Furthermore, as many as 617 gene models including 81 most likely candidate genes (CGs) were identified in the high confidence MQTL regions. These most likely CGs encoded proteins mainly belonging to the following families: B-box-type zinc finger, cytochrome P450 protein, pentatricopeptide repeat, phospholipid/glycerol acyltransferase, F-box protein, small auxin-up RNA, UDP-glucosyltransferase, glutathione S-transferase protein, etc. In addition, ortho-MQTL analysis based on synteny among wheat, rice and barley was also performed which permitted the identification of six ortho-MQTLs among these three cereals. This meta-analysis defines a genome-wide landscape on the most stable and consistent loci associated with reliable molecular markers and candidate genes for salinity tolerance in wheat.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01112-0.

The MYB transcription factor TaPHR3-A1 is involved in phosphate signaling and governs yield-related traits in bread wheat

Improved inorganic phosphate (Pi) use efficiency in crops will be important for sustainable agriculture. Exploring molecular mechanisms that regulate Pi uptake could provide useful information for breeding wheat with improved Pi use efficiency. Here, a TaPHR3-A1 (Gene ID: TraesCS7A02G415800) ortholog of rice OsPHR3 that functions in transcriptional regulation of Pi signaling was cloned from wheat chromosome 7A. Ectopic expression of TaPHR3-A1 in Arabidopsis and rice produced enhanced vegetative growth and more seeds. Overexpression in transgenic rice led to increased biomass, grain number, and primary panicle branching by 61.23, 42.12, and 36.34% compared with the wild type. Transgenic wheat lines with down-regulation of TaPHR3-A1 exhibited retarded growth and root hair development at the seedling stage, and showed yield-related effects at the adult stage when grown in both low- and sufficient Pi conditions, indicating that TaPHR3-A1 positively regulated tolerance to low Pi. Introgression lines further confirmed the effect of TaPHR3-A1 in improving grain number. The Chinese wheat mini core collection and a recombinant inbred line analysis demonstrated that the favorable allele TaPHR3-A1-A associated with higher grain number was positively selected in breeding. A TaPHR3-A1-derived cleaved amplified polymorphic sequence marker effectively identified haplotype TaPHR3-A1-A. Our results suggested that TaPHR3-A1 was a functional regulatory factor for Pi uptake and provided useful information for marker-assisted selection for high yield in wheat.

Molecular genetic analysis of spring wheat core collection using genetic diversity, population structure, and linkage disequilibrium

Background

Wheat (Triticum aestivium L.) is an important crop globally which has a complex genome. To identify the parents with useful agronomic characteristics that could be used in the various breeding programs, it is very important to understand the genetic diversity among global wheat genotypes. Also, understanding the genetic diversity is useful in breeding studies such as marker-assisted selection (MAS), genome-wide association studies (GWAS), and genomic selection.

Results

To understand the genetic diversity in wheat, a set of 103 spring wheat genotypes which represented five different continents were used. These genotypes were genotyped using 36,720 genotyping-by-sequencing derived SNPs (GBS-SNPs) which were well distributed across wheat chromosomes. The tested 103-wheat genotypes contained three different subpopulations based on population structure, principle coordinate, and kinship analyses. A significant variation was found within and among the subpopulations based on the AMOVA. Subpopulation 1 was found to be the more diverse subpopulation based on the different allelic patterns (Na, Ne, I, h, and uh). No high linkage disequilibrium was found between the 36,720 SNPs. However, based on the genomic level, D genome was found to have the highest LD compared with the two other genomes A and B. The ratio between the number of significant LD/number of non-significant LD suggested that chromosomes 2D, 5A, and 7B are the highest LD chromosomes in their genomes with a value of 0.08, 0.07, and 0.05, respectively. Based on the LD decay, the D genome was found to be the lowest genome with the highest number of haplotype blocks on chromosome 2D.

Conclusion

The recent study concluded that the 103-spring wheat genotypes and their GBS-SNP markers are very appropriate for GWAS studies and QTL-mapping. The core collection comprises three different subpopulations. Genotypes in subpopulation 1 are the most diverse genotypes and could be used in future breeding programs if they have desired traits. The distribution of LD hotspots across the genome was investigated which provides useful information on the genomic regions that includes interesting genes.

Gene editing of the wheat homologs of TONNEAU1-recruiting motif encoding gene affects grain shape and weight in wheat

Grain size and weight are important components of a suite of yield-related traits in crops. Here, we showed that the CRISPR-Cas9 gene editing of TaGW7, a homolog of rice OsGW7 encoding a TONNEAU1-recruiting motif (TRM) protein, affects grain shape and weight in allohexaploid wheat. By editing the TaGW7 homoeologs in the B and D genomes, we showed that mutations in either of the two or both genomes increased the grain width and weight but reduced the grain length. The effect sizes of mutations in the TaGW7 gene homoeologs coincided with the relative levels of their expression in the B and D genomes. The effects of gene editing on grain morphology and weight traits were dosage dependent with the double-copy mutant showing larger effect than the respective single copy mutants. The TaGW7-centered gene co-expression network indicated that this gene is involved in the pathways regulating cell division and organ growth, also confirmed by the cellular co-localization of TaGW7 with α- and β-tubulin proteins, the building blocks of microtubule arrays. The analyses of exome capture data in tetraploid domesticated and wild emmer, and hexaploid wheat revealed the loss of diversity around TaGW7-associated with domestication selection, suggesting that TaGW7 is likely to play an important role in the evolution of yield component traits in wheat. Our study showed how integrating CRISPR-Cas9 system with cross-species comparison can help to uncover the function of a gene fixed in wheat for allelic variants targeted by domestication selection and select targets for engineering new gene variants for crop improvement.

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.

Alteration of flavonoid pigmentation patterns during domestication of food crops

Flavonoids are plant pigments that provide health benefits for human and animal consumers. Understanding why domesticated crops have altered pigmentation patterns and unraveling the molecular/genetic mechanisms that underlie this will facilitate the breeding of new (healthier) varieties. We present an overview of changes in flavonoid pigmentation patterns that have occurred during crop domestication and, where possible, link them to the molecular changes that brought about the new phenotypes. We consider species that lost flavonoid pigmentation in the edible part of the plant at some point during domestication (like cereals). We also consider the converse situation, for example eggplant (aubergine), which instead gained strong anthocyanin accumulation in the skin of the fruit during domestication, and some varieties of citrus and apple that acquired anthocyanins in the fruit flesh. Interestingly, the genes responsible for such changes are sometimes closely linked to, or have pleiotropic effects on, important domestication genes, suggesting accidental and perhaps inevitable changes of anthocyanin patterning during domestication. In other cases, flavonoid pigmentation patterns in domesticated crops are the result of cultural preferences, with examples being found in varieties of citrus, barley, wheat, and maize. Finally, and more recently, in some species, anthocyanins seem to have been the direct target of selection in a second wave of domestication that followed the introduction of industrial food processing.

Genotyping by sequencing provides new insights into the diversity of Napier grass (Cenchrus purpureus) and reveals variation in genome-wide LD patterns between collections

Napier grass is an important tropical forage-grass and of growing potential as an energy crop. One-hundred-five Napier grass accessions, encompassing two independent collections, were subjected to genotyping by sequencing which generated a set of high-density genome-wide markers together with short sequence reads. The reads, averaging 54 nucleotides, were mapped to the pearl millet genome and the closest genes and annotation information were used to select candidate genes linked to key forage traits. 980 highly polymorphic SNP markers, distributed across the genome, were used to assess population structure and diversity with seven-subgroups identified. A few representative accessions were selected with the objective of distributing subsets of a manageable size for further evaluation. Genome-wide linkage disequilibrium (LD) analyses revealed a fast LD-decay, on average 2.54 kbp, in the combined population with a slower LD-decay in the ILRI collection compared with the EMBRAPA collection, the significance of which is discussed. This initiative generated high-density markers with a good distribution across the genome. The diversity analysis revealed the existence of a substantial amount of variation in the ILRI collection and identified some unique materials from the EMBRAPA collection, demonstrating the potential of the overall population for further genetic and marker-trait-association studies.

Physical Mapping of Peroxidase Genes and Development of Functional Markers for TaPod-D1 on Bread Wheat Chromosome 7D

Peroxidase (POD) activity in wheat (Triticum aestivum L.) grain influences natural carotenoid pigment content and is associated with the color of flour, and processing and product quality. Here, we report the molecular characterization and physical mapping of POD genes in bread wheat. The complete genomic DNA (gDNA) sequence of two POD genes (TaPod-A2 and TaPod-D1), and the partial gDNA sequence of two additional POD genes (TaPod-A3 and TaPod-B1) from wheat were characterized using in silico cloning and validated through laboratory experiments. Using a set of 21 nullisomic-tetrasomic (NT) lines, six group-7 ditelosomic (Dt) lines, and 38 group-7 deletion (Del) lines of Chinese Spring (CS), TaPod-A2 and TaPod-D1 were found to be physically located on 0.73-0.83 and on the most distal 0.39 fraction arm length (FL) of 7AS and 7DS in cv. CS, respectively; whereas, TaPod-A3 and TaPod-B1 were assigned to the 0.40-0.49 and 0.40-0.48 FL of 7AL and 7BL, respectively. Based on single nucleotide polymorphisms (SNPs) of two alleles at the TaPod-D1 locus, two functional markers POD-7D1 and POD-7D6 were developed, amplifying 540- and 640-bp, fragments in varieties with higher and lower POD activities, respectively. A total of 224 wheat varieties were analyzed and showed a significant association between the polymorphic fragments and POD activity using POD-7D1 and POD-7D6 markers. The analysis of variance (ANOVA) indicated the average POD activities of 115 varieties with TaPod-D1a were significantly lower than 109 varieties with TaPod-D1b (P < 0.01). This study provides useful information of the POD genes in bread wheat, insight into wheat genome synteny and structure, gene-specific markers, and contributes a valuable resource for quality improvement in wheat breeding programs.

Genome-Wide Identification and Transcriptional Regulation of Aquaporin Genes in Bread Wheat (Triticum aestivum L.) under Water Stress

Aquaporins (AQPs) are transmembrane proteins essential for controlling the flow of water and other molecules required for development and stress tolerance in plants, including important crop species such as wheat (Triticum aestivum). In this study, we utilized a genomic approach for analyzing the information about AQPs available in public databases to characterize their structure and function. Furthermore, we validated the expression of a suite of AQP genes, at the transcriptional level, including accessions with contrasting responses to drought, different organs and water stress levels. We found 65 new AQP genes, from which 60% are copies expanded by polyploidization. Sequence analysis of the AQP genes showed that the purifying selection pressure acted on duplicate genes, which was related to a high conservation of the functions. This situation contrasted with the expression patterns observed for different organs, developmental stages or genotypes under water deficit conditions, which indicated functional divergence at transcription. Expression analyses on contrasting genotypes showed high gene transcription from Tonoplast Intrinsic Protein 1 (TIP1) and 2 (TIP2), and Plasma Membrane Intrinsic Protein 1 (PIP1) and 2 (PIP2) subfamilies in roots and from TIP1 and PIP1 subfamilies in leaves. Interestingly, during severe drought stress, 4 TIP genes analyzed in leaves of the tolerant accession reached up to 15-fold the level observed at the susceptible genotype, suggesting a positive relationship with drought tolerance. The obtained results extend our understanding of the structure and function of AQPs, particularly under water stress conditions.

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

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

Genetic Diversity and Population Structure of F3:6 Nebraska Winter Wheat Genotypes Using Genotyping-By-Sequencing

The availability of information on the genetic diversity and population structure in wheat (Triticum aestivum L.) breeding lines will help wheat breeders to better use their genetic resources and manage genetic variation in their breeding program. The recent advances in sequencing technology provide the opportunity to identify tens or hundreds of thousands of single nucleotide polymorphism (SNPs) in large genome species (e.g., wheat). These SNPs can be utilized for understanding genetic diversity and performing genome wide association studies (GWAS) for complex traits. In this study, the genetic diversity and population structure were investigated in a set of 230 genotypes (F_3:6) derived from various crosses as a prerequisite for GWAS and genomic selection. Genotyping-by-sequencing provided 25,566 high-quality SNPs. The polymorphism information content (PIC) across chromosomes ranged from 0.09 to 0.37 with an average of 0.23. The distribution of SNPs markers on the 21 chromosomes ranged from 319 on chromosome 3D to 2,370 on chromosome 3B. The analysis of population structure revealed three subpopulations (G1, G2, and G3). Analysis of molecular variance identified 8% variance among and 92% within subpopulations. Of the three subpopulations, G2 had the highest level of genetic diversity based on three genetic diversity indices: Shannon’s information index (I) = 0.494, diversity index (h) = 0.328 and unbiased diversity index (uh) = 0.331, while G3 had lowest level of genetic diversity (I = 0.348, h = 0.226 and uh = 0.236). This high genetic diversity identified among the subpopulations can be used to develop new wheat cultivars.

Utilization of a Wheat660K SNP array-derived high-density genetic map for high-resolution mapping of a major QTL for kernel number

In crop plants, a high-density genetic linkage map is essential for both genetic and genomic researches. The complexity and the large size of wheat genome have hampered the acquisition of a high-resolution genetic map. In this study, we report a high-density genetic map based on an individual mapping population using the Affymetrix Wheat660K single-nucleotide polymorphism (SNP) array as a probe in hexaploid wheat. The resultant genetic map consisted of 119 566 loci spanning 4424.4 cM, and 119 001 of those loci were SNP markers. This genetic map showed good collinearity with the 90 K and 820 K consensus genetic maps and was also in accordance with the recently released wheat whole genome assembly. The high-density wheat genetic map will provide a major resource for future genetic and genomic research in wheat. Moreover, a comparative genomics analysis among gramineous plant genomes was conducted based on the high-density wheat genetic map, providing an overview of the structural relationships among theses gramineous plant genomes. A major stable quantitative trait locus (QTL) for kernel number per spike was characterized, providing a solid foundation for the future high-resolution mapping and map-based cloning of the targeted QTL.

Utilization of a Wheat660K SNP array-derived high-density genetic map for high-resolution mapping of a major QTL for kernel number

In crop plants, a high-density genetic linkage map is essential for both genetic and genomic researches. The complexity and the large size of wheat genome have hampered the acquisition of a high-resolution genetic map. In this study, we report a high-density genetic map based on an individual mapping population using the Affymetrix Wheat660K single-nucleotide polymorphism (SNP) array as a probe in hexaploid wheat. The resultant genetic map consisted of 119 566 loci spanning 4424.4 cM, and 119 001 of those loci were SNP markers. This genetic map showed good collinearity with the 90 K and 820 K consensus genetic maps and was also in accordance with the recently released wheat whole genome assembly. The high-density wheat genetic map will provide a major resource for future genetic and genomic research in wheat. Moreover, a comparative genomics analysis among gramineous plant genomes was conducted based on the high-density wheat genetic map, providing an overview of the structural relationships among theses gramineous plant genomes. A major stable quantitative trait locus (QTL) for kernel number per spike was characterized, providing a solid foundation for the future high-resolution mapping and map-based cloning of the targeted QTL.

QTL mapping of wheat plant architectural characteristics and their genetic relationship with seven QTLs conferring resistance to sheath blight

Sheath blight is one of the most devastating wheat diseases worldwide. Breeding resistant cultivars is the most powerful strategy to defeat the disease. Plant resistance on “disease escape” works through modulation of morphological traits and shows sustainable resistance to disease. Plant architectural traits have been reported to play a significant role in disease response. Therefore, exploring the genetic relationship between plant architecture and disease resistance is of importance to the understanding of plant resistance via “disease escape”. Using an F₉ population of 266 RILs (Recombinant Inbred Lines) derived from the cross of Luke × AQ24788-83, we have generated a linkage map of 631 markers on 21 chromosomes. In this study, we present the QTL identification of fourteen plant architectural characteristics and heading time from two years and analyze their genetic relationships with seven previously published QTLs to sheath blight (QSBs, QSe.cau), including plant height (PH), the space between the flag leaf and penultimate leaf (fdR), heading date (Hd), and other traits. Twelve stable QTLs of the morphological traits were identified with good consistency across five replicates. For the seven previously published QSBs, we found no significant association with plant height. However, some of the QSBs displayed strong associations with plant architectural traits and heading date. Especially, Q_fdR.cau-1AS, Q_Hd.cau-2BS, Q_fdR.cau-5DL, and Q_fdR.cau-6BL were respectively mapped to the same regions as QSe.cau-1AS, QSe.cau-2BS, QSe.cau-5DL, and QSe.cau-6BL. Taken together, we have demonstrated that plant height did not exert a direct influence on the resistance to sheath blight conferred by the seven QSBs and that the plant architecture and heading date did exhibit a tight relationship with the resistance. Therefore, this study provides a novel evidence to help understand sheath blight resistance in wheat. In addition, the linked morphological characteristics and the generated flanking markers will facilitate breeding for resistance to sheath blight in wheat.

Comparative Analysis of AGPase Genes and Encoded Proteins in Eight Monocots and Three Dicots with Emphasis on Wheat

ADP-glucose pyrophosphorylase (AGPase) is a heterotetrameric enzyme with two large subunits (LS) and two small subunits (SS). It plays a critical role in starch biosynthesis. We are reporting here detailed structure, function and evolution of the genes encoding the LS and the SS among monocots and dicots. "True" orthologs of maize Sh2 (AGPase LS) and Bt2 (AGPase SS) were identified in seven other monocots and three dicots; structure of the enzyme at protein level was also studied. Novel findings of the current study include the following: (i) at the DNA level, the genes controlling the SS are more conserved than those controlling the LS; the variation in both is mainly due to intron number, intron length and intron phase distribution; (ii) at protein level, the SS genes are more conserved relative to those for LS; (iii) "QTCL" motif present in SS showed evolutionary differences in AGPase belonging to wheat 7BS, T. urartu, rice and sorghum, while "LGGG" motif in LS was present in all species except T. urartu and chickpea; SS provides thermostability to AGPase, while LS is involved in regulation of AGPase activity; (iv) heterotetrameric structure of AGPase was predicted and analyzed in real time environment through molecular dynamics simulation for all the species; (v) several cis-acting regulatory elements were identified in the AGPase promoters with their possible role in regulating spatial and temporal expression (endosperm and leaf tissue) and also the expression, in response to abiotic stresses; and (vi) expression analysis revealed downregulation of both subunits under conditions of heat and drought stress. The results of the present study have allowed better understanding of structure and evolution of the genes and the encoded proteins and provided clues for exploitation of variability in these genes for engineering thermostable AGPase.

TaXA21-A1 on chromosome 5AL is associated with resistance to multiple pests in wheat

KEY MESSAGE

The wheat ortholog of the rice gene OsXA21 against bacterial leaf blight showed resistance to multiple pests in bread wheat but different interacting proteins.

ABSTRACT

A quantitative trait locus QYr.osu-5A on the long arm of chromosome 5A in bread wheat (Triticum aestivum L., 2n = 6x = 42; AABBDD) was previously reported to confer consistent resistance in adult plants to predominant stripe rust races, but the gene causing the quantitative trait locus (QTL) is not known. Single-nucleotide polymorphism (SNP) markers were used to saturate the QTL region. Comparative and syntenic regions between wheat and rice (Oryza sativa) were applied to identify candidate genes for QYr.osu-5A. TaXA21-A1, which is referred to as a wheat ortholog of OsXA21-like gene on chromosome 9 in rice, was mapped under the peak of the QYr.osu-5A. TaXA21-A1 not only explained the phenotypic variation in reaction to different stripe rust races but also showed significant effects on resistance to powdery mildew and Hessian fly biotype BP. The natural allelic variation resulted in the alternations of four amino acids in deduced TaXA21-A1 proteins. The interacting proteins of TaXA21-A1 were different from those identified by OsXA21 on rice chromosome 11 against bacterial leaf blight. TaXA21-A1 confers unique resistance against multiple pests in wheat but might not have common protein interactors or thus overlapping functions with OsXA21 in rice. XA21 function has diverged during evolution of cereal crops. The molecular marker developed for TaXA21-A1 would accelerate its application of the candidate gene at the QYr.osu-5A locus in wheat breeding programs.

Genome-Wide QTL Mapping for Wheat Processing Quality Parameters in a Gaocheng 8901/Zhoumai 16 Recombinant Inbred Line Population

Dough rheological and starch pasting properties play an important role in determining processing quality in bread wheat (Triticum aestivum L.). In the present study, a recombinant inbred line (RIL) population derived from a Gaocheng 8901/Zhoumai 16 cross grown in three environments was used to identify quantitative trait loci (QTLs) for dough rheological and starch pasting properties evaluated by Mixograph, Rapid Visco-Analyzer (RVA), and Mixolab parameters using the wheat 90 and 660 K single nucleotide polymorphism (SNP) chip assays. A high-density linkage map constructed with 46,961 polymorphic SNP markers from the wheat 90 and 660 K SNP assays spanned a total length of 4121 cM, with an average chromosome length of 196.2 cM and marker density of 0.09 cM/marker; 6596 new SNP markers were anchored to the bread wheat linkage map, with 1046 and 5550 markers from the 90 and 660 K SNP assays, respectively. Composite interval mapping identified 119 additive QTLs on 20 chromosomes except 4D; among them, 15 accounted for more than 10% of the phenotypic variation across two or three environments. Twelve QTLs for Mixograph parameters, 17 for RVA parameters and 55 for Mixolab parameters were new. Eleven QTL clusters were identified. The closely linked SNP markers can be used in marker-assisted wheat breeding in combination with the Kompetitive Allele Specific PCR (KASP) technique for improvement of processing quality in bread wheat.

A genetic strategy generating wheat with very high amylose content

Resistant starch (RS), a type of dietary fibre, plays an important role in human health; however, the content of RS in most modern processed starchy foods is low. Cereal starch, when structurally manipulated through a modified starch biosynthetic pathway to greatly increase the amylose content, could be an important food source of RS. Transgenic studies have previously revealed the requirement of simultaneous down-regulation of two starch branching enzyme (SBE) II isoforms both located on the long arm of chromosome 2, namely SBEIIa and SBEIIb, to elevate the amylose content in wheat from ~25% to ~75%. The current study revealed close proximity of genes encoding SBEIIa and SBEIIb isoforms in wheat with a genetic distance of 0.5 cM on chromosome 2B. A series of deletion and single nucleotide polymorphism (SNP) loss of function alleles in SBEIIa, SBEIIb or both was isolated from two different wheat populations. A breeding strategy to combine deletions and SNPs generated wheat genotypes with altered expression levels of SBEIIa and SBEIIb, elevating the amylose content to an unprecedented ~85%, with a marked concomitant increase in RS content. Biochemical assays were used to confirm the complete absence in the grain of expression of SBEIIa from all three genomes in combination with the absence of SBEIIb from one of the genomes.

Small-scale gene duplications played a major role in the recent evolution of wheat chromosome 3B

Background

Bread wheat is not only an important crop, but its large (17 Gb), highly repetitive, and hexaploid genome makes it a good model to study the organization and evolution of complex genomes. Recently, we produced a high quality reference sequence of wheat chromosome 3B (774 Mb), which provides an excellent opportunity to study the evolutionary dynamics of a large and polyploid genome, specifically the impact of single gene duplications.

Results

We find that 27 % of the 3B predicted genes are non-syntenic with the orthologous chromosomes of Brachypodium distachyon, Oryza sativa, and Sorghum bicolor, whereas, by applying the same criteria, non-syntenic genes represent on average only 10 % of the predicted genes in these three model grasses. These non-syntenic genes on 3B have high sequence similarity to at least one other gene in the wheat genome, indicating that hexaploid wheat has undergone massive small-scale interchromosomal gene duplications compared to other grasses. Insertions of non-syntenic genes occurred at a similar rate along the chromosome, but these genes tend to be retained at a higher frequency in the distal, recombinogenic regions. The ratio of non-synonymous to synonymous substitution rates showed a more relaxed selection pressure for non-syntenic genes compared to syntenic genes, and gene ontology analysis indicated that non-syntenic genes may be enriched in functions involved in disease resistance.

Conclusion

Our results highlight the major impact of single gene duplications on the wheat gene complement and confirm the accelerated evolution of the Triticeae lineage among grasses.

Electronic supplementary material

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

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.

Chromosomal genomics facilitates fine mapping of a Russian wheat aphid resistance gene

Making use of wheat chromosomal resources, we developed 11 gene-associated markers for the region of interest, which allowed reducing gene interval and spanning it by four BAC clones. Positional gene cloning and targeted marker development in bread wheat are hampered by high complexity and polyploidy of its nuclear genome. Aiming to clone a Russian wheat aphid resistance gene Dn2401 located on wheat chromosome arm 7DS, we have developed a strategy overcoming problems due to polyploidy and enabling efficient development of gene-associated markers from the region of interest. We employed information gathered by GenomeZipper, a synteny-based tool combining sequence data of rice, Brachypodium, sorghum and barley, and took advantage of a high-density linkage map of Aegilops tauschii. To ensure genome- and locus-specificity of markers, we made use of survey sequence assemblies of isolated wheat chromosomes 7A, 7B and 7D. Despite the low level of polymorphism of the wheat D subgenome, our approach allowed us to add in an efficient and cost-effective manner 11 new gene-associated markers in the Dn2401 region and narrow down the target interval to 0.83 cM. Screening 7DS-specific BAC library with the flanking markers revealed a contig of four BAC clones that span the Dn2401 region in wheat cultivar 'Chinese Spring'. With the availability of sequence assemblies and GenomeZippers for each of the wheat chromosome arms, the proposed strategy can be applied for focused marker development in any region of the wheat genome.

Heterologous expression and transcript analysis of gibberellin biosynthetic genes of grasses reveals novel functionality in the GA3ox family

BACKGROUND

The gibberellin (GA) pathway plays a central role in the regulation of plant development, with the 2-oxoglutarate-dependent dioxygenases (2-ODDs: GA20ox, GA3ox, GA2ox) that catalyse the later steps in the biosynthetic pathway of particularly importance in regulating bioactive GA levels. Although GA has important impacts on crop yield and quality, our understanding of the regulation of GA biosynthesis during wheat and barley development remains limited. In this study we identified or assembled genes encoding the GA 2-ODDs of wheat, barley and Brachypodium distachyon and characterised the wheat genes by heterologous expression and transcript analysis.

RESULTS

The wheat, barley and Brachypodium genomes each contain orthologous copies of the GA20ox, GA3ox and GA2ox genes identified in rice, with the exception of OsGA3ox1 and OsGA2ox5 which are absent in these species. Some additional paralogs of 2-ODD genes were identified: notably, a novel gene in the wheat B genome related to GA3ox2 was shown to encode a GA 1-oxidase, named as TaGA1ox-B1. This enzyme is likely to be responsible for the abundant 1β-hydroxylated GAs present in developing wheat grains. We also identified a related gene in barley, located in a syntenic position to TaGA1ox-B1, that encodes a GA 3,18-dihydroxylase which similarly accounts for the accumulation of unusual GAs in barley grains. Transcript analysis showed that some paralogs of the different classes of 2-ODD were expressed mainly in a single tissue or at specific developmental stages. In particular, TaGA20ox3, TaGA1ox1, TaGA3ox3 and TaGA2ox7 were predominantly expressed in developing grain. More detailed analysis of grain-specific gene expression showed that while the transcripts of biosynthetic genes were most abundant in the endosperm, genes encoding inactivation and signalling components were more highly expressed in the seed coat and pericarp.

CONCLUSIONS

The comprehensive expression and functional characterisation of the multigene families encoding the 2-ODD enzymes of the GA pathway in wheat and barley will provide the basis for a better understanding of GA-regulated development in these species. This analysis revealed the existence of a novel, endosperm-specific GA 1-oxidase in wheat and a related GA 3,18-dihydroxylase enzyme in barley that may play important roles during grain expansion and development.

Does morphological and anatomical plasticity during the vegetative stage make wheat more tolerant of water deficit stress than rice?

Water scarcity and the increasing severity of water deficit stress are major challenges to sustaining irrigated rice (Oryza sativa) production. Despite the technologies developed to reduce the water requirement, rice growth is seriously constrained under water deficit stress compared with other dryland cereals such as wheat (Triticum aestivum). We exposed rice cultivars with contrasting responses to water deficit stress and wheat cultivars well adapted to water-limited conditions to the same moisture stress during vegetative growth to unravel the whole-plant (shoot and root morphology) and organ/tissue (root anatomy) responses. Wheat cultivars followed a water-conserving strategy by reducing specific leaf area and developing thicker roots and moderate tillering. In contrast, rice 'IR64' and 'Apo' adopted a rapid water acquisition strategy through thinner roots under water deficit stress. Root diameter, stele and xylem diameter, and xylem number were more responsive and varied with different positions along the nodal root under water deficit stress in wheat, whereas they were relatively conserved in rice cultivars. Increased metaxylem diameter and lower metaxylem number near the root tips and exactly the opposite phenomena at the root-shoot junction facilitated the efficient use of available soil moisture in wheat. Tolerant rice 'Nagina 22' had an advantage in root morphological and anatomical attributes over cultivars IR64 and Apo but lacked plasticity, unlike wheat cultivars exposed to water deficit stress. The key traits determining the adaptation of wheat to dryland conditions have been summarized and discussed.

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.

Precisely mapping a major gene conferring resistance to Hessian fly in bread wheat using genotyping-by-sequencing

Background

One of the reasons hard red winter wheat cultivar ‘Duster’ (PI 644016) is widely grown in the southern Great Plains is that it confers a consistently high level of resistance to biotype GP of Hessian fly (Hf). However, little is known about the genetic mechanism underlying Hf resistance in Duster. This study aimed to unravel complex structures of the Hf region on chromosome 1AS in wheat by using genotyping-by-sequencing (GBS) markers and single nucleotide polymorphism (SNP) markers.

Results

Doubled haploid (DH) lines generated from a cross between two winter wheat cultivars, ‘Duster’ and ‘Billings’ , were used to identify genes in Duster responsible for effective and consistent resistance to Hf. Segregation in reaction of the 282 DH lines to Hf biotype GP fit a one-gene model. The DH population was genotyped using 2,358 markers developed using the GBS approach. A major QTL, explaining 88% of the total phenotypic variation, was mapped to a chromosome region that spanned 178 cM and contained 205 GBS markers plus 1 SSR marker and 1 gene marker, with 0.86 cM per marker in genetic distance. The analyses of GBS marker sequences and further mapping of SSR and gene markers enabled location of the QTL-containing linkage group on the short arm of chromosome 1A. Comparative mapping of the common markers for the gene for QHf.osu-1A^d in Duster and the Hf-resistance gene for QHf.osu-1A⁷⁴ in cultivar ‘2174’ showed that the two Hf resistance genes are located on the same chromosome arm 1AS, only 11.2 cM apart in genetic distance. The gene at QHf.osu-1A^d in Duster has been delimited within a 2.7 cM region.

Conclusion

Two distinct resistance genes exist on the short arm of chromosome 1A as found in the two hard red winter cultivars, 2174 and Duster. Whereas the Hf resistance gene in 2174 is likely allelic to one or more of the previously mapped resistance genes (H9, H10, H11, H16, or H17) in wheat, the gene in Duster is novel and confers a more consistent phenotype than 2174 in response to biotype GP infestation in controlled-environment assays.

Electronic supplementary material

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

Isolation of candidate genes for the barley Ant1 and wheat Rc genes controlling anthocyanin pigmentation in different vegetative tissues

MYB transcription factors exist in a large copy number and control various plant phenotypes. We cloned R2R3 MYB-type transcription factors that determine the coloration of basal sheaths in barley and wheat coleoptiles. These genes are highly homologous to maize C1 and rice OsC1, regulators for anthocyanin biosynthesis, but they control seed pigmentation in maize and rice. On the basis of high homology, barley and wheat counterparts are designated HvC1 and TaC1, respectively. HvC1 gene is located on the short arm of chromosome 7H, and TaC1 genes are located on the short arms of chromosomes 7A, 7B, and 7D (TaC1-A1, B1, and D1, respectively). HvC1 is a strong candidate for Ant1 because of (1) complete co-segregation of anthocyanin pigmentation phenotype of the basal sheath with the HvC1 genotype in genetic mapping, and (2) complete deletion of the HvCl gene in two anthocyanin-decreased allelic mutants (ant1.1 and ant1.2) that were induced by irradiation. In contrast, colorless coleoptile wheat lines had lesions in all three genomes consisting of a single-nucleotide substitution or a 1-bp deletion of TaC1-A1, a 1.7-kb insertion of TaC1-B1, and a 2.0-kb insertion of TaC1-D1. At least one normal TaC1 gene appears to be sufficient to produce anthocyanin pigments in wheat coleoptiles. Previous crossing experiments localized Rc (red coleoptile) genes to homoeologous group 7 chromosomes and deduced Rc genotypes of several wheat lines. Their TaC1 gene sequence variation coincided with deduced Rc genotypes; therefore, the present molecular genetic study demonstrates that TaC1 is a strong candidate for Rc in wheat.

Identification of candidate genes, regions and markers for pre-harvest sprouting resistance in wheat (Triticum aestivum L.)

BACKGROUND

Pre-harvest sprouting (PHS) of wheat grain leads to a reduction in grain yield and quality. The availability of markers for marker-assisted selection (MAS) of PHS resistance will serve to enhance breeding selection and advancement of lines for cultivar development. The aim of this study was to identify candidate regions and develop molecular markers for PHS resistance in wheat. This was achieved via high density mapping of single nucleotide polymorphism (SNP) markers from an Illumina 90 K Infinium Custom Beadchip in a doubled haploid (DH) population derived from a RL4452/'AC Domain' cross and subsequent detection of quantitative trait loci (QTL) for PHS related traits (falling number [FN], germination index [GI] and sprouting index [SI]). SNP marker sequences flanking QTL were used to locate colinear regions in Brachypodium and rice, and identify genic markers associated with PHS resistance that can be utilized for MAS in wheat.

RESULTS

A linkage map spanning 2569.4 cM was constructed with a total of 12,201 SNP, simple sequence repeat (SSR), diversity arrays technology (DArT) and expressed sequence tag (EST) markers. QTL analyses using Multiple Interval Mapping (MIM) identified four QTL for PHS resistance traits on chromosomes 3B, 4A, 7B and 7D. Sequences of SNPs flanking these QTL were subject to a BLASTN search on the International Wheat Genome Sequencing Consortium (IWGSC) database (http://wheat-urgi.versailles.inra.fr/Seq-Repository). Best survey sequence hits were subject to a BLASTN search on Gramene (www.gramene.org) against both Brachypodium and rice databases, and candidate genes and regions for PHS resistance were identified. A total of 18 SNP flanking sequences on chromosomes 3B, 4A, 7B and 7D were converted to KASP markers and validated with matching genotype calls of Infinium SNP data.

CONCLUSIONS

Our study identified candidate genes involved in abscissic acid (ABA) and gibberellin (GA) metabolism, and flowering time in four genomic regions of Brachypodium and rice respectively, in addition to 18 KASP markers for PHS resistance in wheat. These markers can be deployed in future genetic studies of PHS resistance and might also be useful in the evaluation of PHS in germplasm and breeding material.

TEF-7A, a transcript elongation factor gene, influences yield-related traits in bread wheat (Triticum aestivum L.)

In this study, TaTEF-7A, a member of the transcript elongation factor gene family, and its flanking sequences were isolated. TaTEF-7A was located on chromosome 7A and was flanked by markers Xwmc83 and XP3156.3. Subcellular localization revealed that TaTEF-7A protein was localized in the nucleus. This gene was expressed in all organs, but the highest expression occurred in young spikes and developing seeds. Overexpression of TaTEF-7A in Arabidopsis thaliana produced pleiotropic effects on vegetative and reproductive development that enhanced grain length, silique number, and silique length. No diversity was found in the coding region of TaTEF-7A, but 16 single nucleotide polymorphisms and Indels were detected in the promoter regions of different cultivars. Markers based on sequence variations in the promoter regions (InDel-629 and InDel-604) were developed, and three haplotypes were identified based on those markers. Haplotype-trait association analysis of the Chinese wheat mini core collection revealed that TaTEF-7A was significantly associated with grain number per spike. Phenotyping of near-isogenic lines (NILs) confirmed that TaTEF-7A increases potential grain yield and yield-related traits. Frequency changes in favoured haplotypes gradually increased in cultivars released in China from the 1940s. Geographic distributions of favoured haplotypes were characterized in six major wheat production regions worldwide. The presence of Hap-7A-3, the favoured haplotype, showed a positive correlation with yield in a global set of breeding lines. These results suggest that TaTEF-7A is a functional regulatory factor for grain number per spike and provide a basis for marker-assisted selection.

TEF-7A, a transcript elongation factor gene, influences yield-related traits in bread wheat (Triticum aestivum L.)

In this study, TaTEF-7A, a member of the transcript elongation factor gene family, and its flanking sequences were isolated. TaTEF-7A was located on chromosome 7A and was flanked by markers Xwmc83 and XP3156.3. Subcellular localization revealed that TaTEF-7A protein was localized in the nucleus. This gene was expressed in all organs, but the highest expression occurred in young spikes and developing seeds. Overexpression of TaTEF-7A in Arabidopsis thaliana produced pleiotropic effects on vegetative and reproductive development that enhanced grain length, silique number, and silique length. No diversity was found in the coding region of TaTEF-7A, but 16 single nucleotide polymorphisms and Indels were detected in the promoter regions of different cultivars. Markers based on sequence variations in the promoter regions (InDel-629 and InDel-604) were developed, and three haplotypes were identified based on those markers. Haplotype-trait association analysis of the Chinese wheat mini core collection revealed that TaTEF-7A was significantly associated with grain number per spike. Phenotyping of near-isogenic lines (NILs) confirmed that TaTEF-7A increases potential grain yield and yield-related traits. Frequency changes in favoured haplotypes gradually increased in cultivars released in China from the 1940s. Geographic distributions of favoured haplotypes were characterized in six major wheat production regions worldwide. The presence of Hap-7A-3, the favoured haplotype, showed a positive correlation with yield in a global set of breeding lines. These results suggest that TaTEF-7A is a functional regulatory factor for grain number per spike and provide a basis for marker-assisted selection.

Genetic and molecular bases of yield-associated traits: a translational biology approach between rice and wheat

Transferring the knowledge bases between related species may assist in enlarging the yield potential of crop plants. Being cereals, rice and wheat share a high level of gene conservation; however, they differ at metabolic levels as a part of the environmental adaptation resulting in different yield capacities. This review focuses on the current understanding of genetic and molecular regulation of yield-associated traits in both crop species, highlights the similarities and differences and presents the putative knowledge gaps. We focus on the traits associated with phenology, photosynthesis, and assimilate partitioning and lodging resistance; the most important drivers of yield potential. Currently, there are large knowledge gaps in the genetic and molecular control of such major biological processes that can be filled in a translational biology approach in transferring genomics and genetics informations between rice and wheat.

Genetic rearrangements of six wheat-agropyron cristatum 6P addition lines revealed by molecular markers

Agropyron cristatum (L.) Gaertn. (2n = 4x = 28, PPPP) not only is cultivated as pasture fodder but also could provide many desirable genes for wheat improvement. It is critical to obtain common wheat–A. cristatum alien disomic addition lines to locate the desired genes on the P genome chromosomes. Comparative analysis of the homoeologous relationships between the P genome chromosome and wheat genome chromosomes is a key step in transferring different desirable genes into common wheat and producing the desired alien translocation line while compensating for the loss of wheat chromatin. In this study, six common wheat–A. cristatum disomic addition lines were produced and analyzed by phenotypic examination, genomic in situ hybridization (GISH), SSR markers from the ABD genomes and STS markers from the P genome. Comparative maps, six in total, were generated and demonstrated that all six addition lines belonged to homoeologous group 6. However, chromosome 6P had undergone obvious rearrangements in different addition lines compared with the wheat chromosome, indicating that to obtain a genetic compensating alien translocation line, one should recombine alien chromosomal regions with homoeologous wheat chromosomes. Indeed, these addition lines were classified into four types based on the comparative mapping: 6P_I, 6P_II, 6P_III, and 6P_IV. The different types of chromosome 6P possessed different desirable genes. For example, the 6P_I type, containing three addition lines, carried genes conferring high numbers of kernels per spike and resistance to powdery mildew, important traits for wheat improvement. These results may prove valuable for promoting the development of conventional chromosome engineering techniques toward molecular chromosome engineering.

Structural-functional dissection and characterization of yield-contributing traits originating from a group 7 chromosome of the wheatgrass species Thinopyrum ponticum after transfer into durum wheat

For the first time, using chromosome engineering of durum wheat, the underlying genetic determinants of a yield-improving segment from Thinopyrum ponticum (7AgL) were dissected. Three durum wheat-Th. ponticum near-isogenic recombinant lines (NIRLs), with distal portions of their 7AL arm (fractional lengths 0.77, 0.72, and 0.60) replaced by alien chromatin, were field-tested for two seasons under rainfed conditions. Yield traits and other agronomic characteristics of the main shoot and whole plant were measured. Loci for seed number per ear and per spikelet were detected in the proximal 7AgL segment (0.60-0.72). Loci determining considerable increases of flag leaf width and area, productive tiller number per plant, biomass per plant, and grain yield per plant were located in the distally adjacent 0.72-0.77 7AgL segment, while in the most distal portion (0.77-1.00) genetic effects on spikelet number per ear were identified. Contrary to previous reports, trials with the bread wheat T4 translocation line, carrying on 7DL a sizeable 7AgL segment of which those present in the durum wheat-Th. ponticum NIRLs represent fractions, gave no yield advantage. The hypothesis that ABA might be a factor contributing to the 7AgL effects was tested by analysing endogenous ABA contents of the NIRLs and their responses to exogenous ABA application. The 7AgL yield-related loci were shown to be ABA-independent. This study highlights the value of wheat-alien recombinant lines for dissecting the genetic and physiological basis of complex traits present in wild germplasm, and provides a basis for their targeted exploitation in wheat breeding.

Preparation of high molecular weight gDNA and bacterial artificial chromosome (BAC) libraries in plants

Bacterial artificial chromosome (BAC) libraries are extremely valuable large-insert DNA libraries for physical mapping, positional cloning, comparative genomic analysis, complete genome sequencing, and evolutionary studies. Due to their stability and relative simplicity BAC libraries are most preferred over other approaches for cloning large genomic DNA fragments for large-insert libraries. Isolation of intact high molecular weight (HMW) DNA is a critical step underlying the success of large-insert genomic DNA library construction. It requires the isolation of purified nuclei, embedding them into LMP agarose plugs, restriction digestion of the plugs, and quite often size selection using PFGE and electro-elution of insert DNA. The construction of BAC libraries is complex and challenging for most molecular laboratories. To facilitate the construction of BAC libraries, we present a step-by-step protocol for isolation of HMW DNA and construction of plant BAC libraries.

Genetic mapping and marker development for resistance of wheat against the root lesion nematode Pratylenchus neglectus

BACKGROUND

The Rlnn1 locus, which resides on chromosome 7A of bread wheat (Triticum aestivum L.) confers moderate resistance against the root lesion nematode Pratylenchus neglectus. Prior to this research, the exact linkage relationships of Rlnn1 with other loci on chromosome 7A were not clear and there were no simple codominant markers available for selection of Rlnn1 in wheat breeding. The objectives of the research reported here were to (1) develop an improved genetic map of the Rlnn1 region of chromosome 7A and (2) develop molecular markers that could be used in marker-assisted selection to improve resistance of wheat against P. neglectus.

RESULTS

A large-effect quantitative trait locus (QTL) for resistance against P. neglectus was genetically mapped using a population of Excalibur/Kukri doubled haploid lines. This QTL coincides in position with the rust resistance gene(s) Lr20/Sr15, the phytoene synthase gene Psy-A1 and 10 molecular markers, including five new markers designed using wheat-rice comparative genomics and wheat expressed sequence tags. Two of the new markers are suitable for use as molecular diagnostic tools to distinguish plants that carry Rlnn1 and Lr20/Sr15 from those that do not carry these resistance genes.

CONCLUSIONS

The genomic location of Rlnn1 was confirmed to be in the terminal region of the long arm of chromosome 7A. Molecular markers were developed that provide simple alternatives to costly phenotypic assessment of resistance against P. neglectus in wheat breeding. In Excalibur, genetic recombination seems to be completely suppressed in the Rlnn1 region.

Comparative high-resolution mapping of the wax inhibitors Iw1 and Iw2 in hexaploid wheat

The wax (glaucousness) on wheat leaves and stems is mainly controlled by two sets of genes: glaucousness loci (W1 and W2) and non-glaucousness loci (Iw1 and Iw2). The non-glaucousness (Iw) loci act as inhibitors of the glaucousness loci (W). High-resolution comparative genetic linkage maps of the wax inhibitors Iw1 originating from Triticum dicoccoides, and Iw2 from Aegilops tauschii were developed by comparative genomics analyses of Brachypodium, sorghum and rice genomic sequences corresponding to the syntenic regions of the Iw loci in wheat. Eleven Iw1 and eight Iw2 linked EST markers were developed and mapped to linkage maps on the distal regions of chromosomes 2BS and 2DS, respectively. The Iw1 locus mapped within a 0.96 cM interval flanked by the BE498358 and CA499581 EST markers that are collinear with 122 kb, 202 kb, and 466 kb genomic regions in the Brachypodium 5S chromosome, the sorghum 6S chromosome and the rice 4S chromosome, respectively. The Iw2 locus was located in a 4.1 to 5.4-cM interval in chromosome 2DS that is flanked by the CJ886319 and CJ519831 EST markers, and this region is collinear with a 2.3 cM region spanning the Iw1 locus on chromosome 2BS. Both Iw1 and Iw2 co-segregated with the BF474014 and CJ876545 EST markers, indicating they are most likely orthologs on 2BS and 2DS. These high-resolution maps can serve as a framework for chromosome landing, physical mapping and map-based cloning of the wax inhibitors in wheat.

Comparative genome analysis between Agrostis stolonifera and members of the Pooideae subfamily, including Brachypodium distachyon

Creeping bentgrass (Agrostis stolonifera, allotetraploid 2n = 4x = 28) is one of the major cool-season turfgrasses. It is widely used on golf courses due to its tolerance to low mowing and aggressive growth habit. In this study, we investigated genome relationships of creeping bentgrass relative to the Triticeae (a consensus map of Triticum aestivum, T. tauschii, Hordeum vulgare, and H. spontaneum), oat, rice, and ryegrass maps using a common set of 229 EST-RFLP markers. The genome comparisons based on the RFLP markers revealed large-scale chromosomal rearrangements on different numbers of linkage groups (LGs) of creeping bentgrass relative to the Triticeae (3 LGs), oat (4 LGs), and rice (8 LGs). However, we detected no chromosomal rearrangement between creeping bentgrass and ryegrass, suggesting that these recently domesticated species might be closely related, despite their memberships to different Pooideae tribes. In addition, the genome of creeping bentgrass was compared with the complete genome sequence of Brachypodium distachyon in Pooideae subfamily using both sequences of the above-mentioned mapped EST-RFLP markers and sequences of 8,470 publicly available A. stolonifera ESTs (AgEST). We discovered large-scale chromosomal rearrangements on six LGs of creeping bentgrass relative to B. distachyon. Also, a total of 24 syntenic blocks based on 678 orthologus loci were identified between these two grass species. The EST orthologs can be utilized in further comparative mapping of Pooideae species. These results will be useful for genetic improvement of Agrostis species and will provide a better understanding of evolution within Pooideae species.