Characterization of wheat DArT markers: genetic and functional features

Genome-wide analysis for root and leaf architecture traits associated with drought tolerance at the seedling stage in a highly ecologically diverse wheat population

Drought stress occurred at early growth stages in wheat affecting the following growth stages. Therefore, selecting promising drought-tolerant genotypes with highly adapted traits at the seedling stage is an important task for wheat breeders and geneticists. Few research efforts were conducted on the genetic control for drought-adaptive traits at the seedling stage in wheat. In this study, a set of 146 highly diverse spring wheat core collections representing 28 different countries was evaluated under drought stress at the seedling stage. All genotypes were exposed to drought stress for 13 days by water withholding. Leaf traits including seedling length, leaf wilting, days to wilting, leaf area, and leaf rolling were scored. Moreover, root traits such as root length, maximum width, emergence angle, tip angle, and number of roots were scored. Considerable significant genetic variation was found among all genotypes tested in these experiments. The heritability estimates ranged from 0.74 (leaf witling) to 0.99 (root tip angle). A set of nine genotypes were selected and considered drought-tolerant genotypes. Among all leaf traits, shoot length had significant correlations with all root traits under drought stress. The 146 genotypes were genotyped using the Infinium Wheat 15 K single nucleotide polymorphism (SNP) array and diversity arrays technology (DArT) marker platform. The result of genotyping revealed 12,999 SNPs and 2150 DArT markers which were used to run a genome-wide association study (GWAS). The results of GWAS revealed 169 markers associated with leaf and root traits under drought stress. Out of the 169 markers, 82 were considered major quantitative trait loci (QTL). The GWAS revealed 95 candidate genes were identified with 53 genes showing evidence for drought tolerance in wheat, while the remaining candidate genes were considered novel. No shared markers were found between leaf and root traits. The results of the study provided mapping novel markers associated with new root traits at the seedling stage. Also, the selected genotypes from different countries could be employed in future wheat breeding programs not only for improving adaptive drought-tolerant traits but also for expanding genetic diversity.

Haplotype-Based Genome-Wide Association Analysis Using Exome Capture Assay and Digital Phenotyping Identifies Genetic Loci Underlying Salt Tolerance Mechanisms in Wheat

Soil salinity can impose substantial stress on plant growth and cause significant yield losses. Crop varieties tolerant to salinity stress are needed to sustain yields in saline soils. This requires effective genotyping and phenotyping of germplasm pools to identify novel genes and QTL conferring salt tolerance that can be utilised in crop breeding schemes. We investigated a globally diverse collection of 580 wheat accessions for their growth response to salinity using automated digital phenotyping performed under controlled environmental conditions. The results show that digitally collected plant traits, including digital shoot growth rate and digital senescence rate, can be used as proxy traits for selecting salinity-tolerant accessions. A haplotype-based genome-wide association study was conducted using 58,502 linkage disequilibrium-based haplotype blocks derived from 883,300 genome-wide SNPs and identified 95 QTL for salinity tolerance component traits, of which 54 were novel and 41 overlapped with previously reported QTL. Gene ontology analysis identified a suite of candidate genes for salinity tolerance, some of which are already known to play a role in stress tolerance in other plant species. This study identified wheat accessions that utilise different tolerance mechanisms and which can be used in future studies to investigate the genetic and genic basis of salinity tolerance. Our results suggest salinity tolerance has not arisen from or been bred into accessions from specific regions or groups. Rather, they suggest salinity tolerance is widespread, with small-effect genetic variants contributing to different levels of tolerance in diverse, locally adapted germplasm.

A new approach fits multivariate genomic prediction models efficiently

Background

Fast, memory-efficient, and reliable algorithms for estimating genomic estimated breeding values (GEBV) for multiple traits and environments are needed to make timely decisions in breeding. Multivariate genomic prediction exploits genetic correlations between traits and environments to increase accuracy of GEBV compared to univariate methods. These genetic correlations are estimated simultaneously with GEBV, because they are specific to year, environment, and management. However, estimating genetic parameters is computationally demanding with restricted maximum likelihood (REML) and Bayesian samplers, and canonical transformations or orthogonalizations cannot be used for unbalanced experimental designs.

Methods

We propose a multivariate randomized Gauss–Seidel algorithm for simultaneous estimation of model effects and genetic parameters. Two previously proposed methods for estimating genetic parameters were combined with a Gauss–Seidel (GS) solver, and were called Tilde-Hat-GS (THGS) and Pseudo-Expectation-GS (PEGS). Balanced and unbalanced experimental designs were simulated to compare runtime, bias and accuracy of GEBV, and bias and standard errors of estimates of heritabilities and genetic correlations of THGS, PEGS, and REML. Models with 10 to 400 response variables, 1279 to 42,034 genetic markers, and 5990 to 1.85 million observations were fitted.

Results

Runtime of PEGS and THGS was a fraction of REML. Accuracies of GEBV were slightly lower than those from REML, but higher than those from the univariate approach, hence THGS and PEGS exploited genetic correlations. For 500 to 600 observations per response variable, biases of estimates of genetic parameters of THGS and PEGS were small, but standard errors of estimates of genetic correlations were higher than for REML. Bias and standard errors decreased as sample size increased. For balanced designs, GEBV and estimates of genetic correlations from THGS were unbiased when only an intercept and eigenvectors of genotype scores were fitted.

Conclusions

THGS and PEGS are fast and memory-efficient algorithms for multivariate genomic prediction for balanced and unbalanced experimental designs. They are scalable for increasing numbers of environments and genetic markers. Accuracy of GEBV was comparable to REML. Estimates of genetic parameters had little bias, but their standard errors were larger than for REML. More studies are needed to evaluate the proposed methods for datasets that contain selection.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00730-w.

Genetic mapping of adult-plant resistance genes to powdery mildew in triticale

Triticale is a cereal of high economic importance; however, along with the increase in the area of this cereal, it is more often infected by the fungal pathogen Blumeria graminis, which causes powdery mildew. The rapid development of molecular biology techniques, in particular methods based on molecular markers may be an important tool used in modern plant breeding. Development of genetic maps, location of the QTLs defining the region of the genome associated with resistance and selection of markers linked to particular trait can be used to select resistant genotypes as well as to pyramidize several resistance genes in one variety. In this paper, we present a new, high-density genetic map of triticale doubled haploids (DH) population “Grenado” × “Zorro” composed of DArT, silicoDArT, and SNP markers. Composite interval mapping method was used to detect eight QTL regions associated with the area under disease progress curve (AUDPC) and 15 regions with the average value of powdery mildew infection (avPM) based on observation conducted in 3-year period in three different locations across the Poland. Two regions on rye chromosome 4R, and single loci on 5R and 6R were reported for the first time as regions associated with powdery mildew resistance. Among all QTLs, 14 candidate genes were identified coded cyclin-dependent kinase, serine/threonine-protein kinase-like protein as well as AMEIOTIC 1 homolog DYAD-like protein, DETOXIFICATION 16-like protein, and putative disease resistance protein RGA3. Three of identified candidate genes were found among newly described QTL regions associated with powdery mildew resistance in triticale.

Identifying suitable tester for evaluating Striga resistant lines using DArTseq markers and agronomic traits

A desirable tester that elicits greater genetic difference in Striga resistance among test crosses in a breeding program has not been reported. Therefore, this study was conducted to characterize 30 Striga resistant yellow endosperm maize inbred lines and three testers with varying resistance levels to Striga using DArTseq SNP markers and agronomic traits to identify a suitable tester for resistance hybrid breeding. Marker-based and agronomic trait-based genetic distances were estimated for yellow endosperm maize inbred lines and testers with varying resistance levels to Striga. The Marker-based cluster analysis separated the Striga resistant lines and testers into two distinct groups. Although the susceptible tester (T3) was the most distantly related to the 30 Striga resistant inbred lines, it exhibited a narrower range in genetic distance estimates and poor agronomic performance under Striga infestation in crosses with the resistant lines. In contrast, the resistant tester (T2) showed a broader range in genetic distance estimates in pairs with the 30 resistant lines. Also, it formed many high yielding hybrids with desirable traits under parasite pressure. Furthermore, the most significant positive association between agronomic trait-based and marker-based distance estimates (r = 0.389, P = 0.01) was observed when T2 has paired with the Striga resistant maize inbred lines. It thus appears that T2 may be used as a suitable tester to determine the breeding value of lines in hybrid maize resistance breeding programs. T2 was the most suitable tester, with a tolerant tester (T1) as an alternative tester to characterize the combining ability of Striga resistant maize inbred lines. This result can also encourage other breeders to investigate testers relative discriminating ability with varying levels of resistance in hybrid breeding for resistance to diseases, pests, and other parasitic plants.

Toward the Genetic Basis and Multiple QTLs of Kernel Hardness in Wheat

Kernel hardness is one of the most important single traits of wheat seed. It classifies wheat cultivars, determines milling quality and affects many end-use qualities. Starch granule surfaces, polar lipids, storage protein matrices and Puroindolines potentially form a four-way interaction that controls wheat kernel hardness. As a genetic factor, Puroindoline polymorphism explains over 60% of the variation in kernel hardness. However, genetic factors other than Puroindolines remain to be exploited. Over the past two decades, efforts using population genetics have been increasing, and numerous kernel hardness-associated quantitative trait loci (QTLs) have been identified on almost every chromosome in wheat. Here, we summarize the state of the art for mapping kernel hardness. We emphasize that these steps in progress have benefitted from (1) the standardized methods for measuring kernel hardness, (2) the use of the appropriate germplasm and mapping population, and (3) the improvements in genotyping methods. Recently, abundant genomic resources have become available in wheat and related Triticeae species, including the high-quality reference genomes and advanced genotyping technologies. Finally, we provide perspectives on future research directions that will enhance our understanding of kernel hardness through the identification of multiple QTLs and will address challenges involved in fine-tuning kernel hardness and, consequently, food properties.

QTL Mapping Using a High-Density Genetic Map to Identify Candidate Genes Associated With Metribuzin Tolerance in Hexaploid Wheat (Triticum aestivum L.)

Tolerance to metribuzin, a broad-spectrum herbicide, is an important trait for weed control in wheat breeding. However, the genetics of metribuzin tolerance in relation to the underlying quantitative trait loci (QTL) and genes is limited. This study developed F₈ recombinant inbred lines (RILs) from a cross between a highly resistant genotype (Chuan Mai 25) and highly susceptible genotype (Ritchie), which were used for QTL mapping of metribuzin tolerance. Genotyping was done using a diversity arrays technology sequencing (DArTseq) platform, and phenotyping was done in controlled environments. Herbicide tolerance was measured using three traits, visual score (VS), reduction of chlorophyll content (RCC), and mean value of chlorophyll content for metribuzin-treated plants (MCC). A high-density genetic linkage map was constructed using 2,129 DArTseq markers. Inclusive composite interval mapping (ICIM) identified seven QTL, one each on chromosomes 2A, 2D, 3A, 3B, 4A, 5A, and 6A. Three major QTL-Qrcc.uwa.2AS, Qrcc.uwa.5AL, and Qrcc.uwa.6AL-explained 11.39%, 11.06%, and 11.45% of the phenotypic variation, respectively. The 5A QTL was further validated using kompetitive allele-specific PCR (KASP) assays in an F₃ validation population developed from Chuan Mai 25 × Dagger. Blasting the single-nucleotide polymorphisms (SNPs) flanking the QTL in the wheat reference genome RefV1.0 revealed SNP markers within or very close to annotated genes which could be candidate genes responsible for metribuzin tolerance. Most of the candidate genes were related to metabolic detoxification, especially those of P450 pathway and xenobiotic transmembrane transporter activity, which are reportedly key molecules responsible for herbicide tolerance. This study is the first to use specially developed populations to conduct QTL mapping on the metribuzin tolerance trait. The three major QTL and candidate genes identified in this study could facilitate marker-assisted metribuzin breeding in wheat. The QTL could be fine-mapped to locate the genes responsible for metribuzin tolerance, which could be introgressed into elite wheat cultivars.

Cereal Root Interactions with Soilborne Pathogens—From Trait to Gene and Back

Realizing the yield potential of crop plants in the presence of shifting pathogen populations, soil quality, rainfall, and other agro-environmental variables remains a challenge for growers and breeders worldwide. In this review, we discuss current approaches for combatting the soilborne phytopathogenic nematodes, Pratylenchus and Heterodera of wheat and barley, and Meloidogyne graminicola Golden and Birchfield, 1965 of rice. The necrotrophic fungal pathogens, Rhizoctonia solani Kühn 1858 AG-8 and Fusarium spp. of wheat and barley, also are discussed. These pathogens constitute major causes of yield loss in small-grain cereals of the Pacific Northwest, USA and throughout the world. Current topics include new sources of genetic resistance, molecular leads from whole genome sequencing and genome-wide patterns of hosts, nematode or fungal gene expression during root-pathogen interactions, host-induced gene silencing, and building a molecular toolbox of genes and regulatory sequences for deployment of resistance genes. In conclusion, improvement of wheat, barley, and rice will require multiple approaches.

Characterizing Palestinian snake melon (Cucumis melo var. flexuosus) germplasm diversity and structure using SNP and DArTseq markers

BACKGROUND

Crop landraces embody a source of beneficial genes potentially providing endurance to environmental stress and other agronomic qualities including yield. Our study included 88 snake melon accessions (Cucumis melo var. flexuosus) collected from 9 districts in the Palestinian West-Bank. These accessions represent four landraces of Palestinian snake melon: Green, and White Baladi, and Green, and White Sahouri.

RESULTS

This is the first report on successful application of genotyping by sequencing in snake melon. Nine thousand seven hundred fifty single-nucleotide polymorphism (SNP) and 7400 DArTseq genetic markers were employed to evaluate genetic biodiversity and population structure of Palestinian snake melon germplasm collection. Clustering based on neighbor-joining-analysis, principle coordinate and Bayesian model implemented in Structure showed that patterns of genetic diversity of snake melon landraces depends on their geographical source and unraveled the presence of two major local landraces (Sahouri, and Baladi) with accessions from each group clustering together. A significant correlation was observed between both types of markers in Mantel correlation test. A significant association between genetic and geographic matrices (P < 0.0001) was also detected. AMOVA indicated that majority of variation (90%) was due to the difference within accessions.

CONCLUSION

The Palestinian landraces seem to have unique genes that may allow the enhancement of the global snake melon gene pool and developments of the plant production worldwide. Our subsequent objective is to detect genotypes with promising qualities and to conduct association mapping studies concentrating on Fusarium-wilt resistance, yield, and environmental stresses.

Genetic Diversity and Molecular Characterization of Iranian Durum Wheat Landraces (Triticum turgidum durum (Desf.) Husn.) Using DArT Markers

Durum wheat is grown mainly in rain-fed regions of Iran and the Mediterranean district under stressful conditions. Different environmental conditions and agricultural practices among ancient communities have led to the development of locally adapted genotypes known as landraces. Landraces are a valued source of genetic variety and show definite adaptation to local environmental conditions according to their home of origin. This study aimed to explore linkage disequilibrium (LD) analysis and the population structure and genetic diversity of Iranian durum wheat landraces. In this study, population structure and genome-wide LD were investigated in 129 durum landrace accessions using 1500 DArT markers. Both structure and discriminant analysis of principal components obviously subdivided the sample collection into seven distinct groups centered on key ancestors and regions of origin of the germplasm. Genetic diversity among the populations was primarily within population (68 vs. 32%). Mean LD values across the entire population sample decayed below r² of 0.11 after 1 cM. LD decay of genomes A and B of Iranian durum wheat landrace is approximately 2-3 cM (r² = 0.11) and approximately 0.5 cM (r² = 0.12), respectively. Altogether, low LD decay, a high number of subpopulations, and the high existence of genetic diversity among and within populations were characteristics of the Iranian durum landrace collection. Hence, the existing genetic diversity within the population can be associated with the very long evolutionary history of plants in Iran. The populations we studied are hence presented as a valuable resource that can be used in basic and applied research in durum wheat breeding.

Regulation and Evolution of NLR Genes: A Close Interconnection for Plant Immunity

NLR (NOD-like receptor) genes belong to one of the largest gene families in plants. Their role in plants' resistance to pathogens has been clearly described for many members of this gene family, and dysregulation or overexpression of some of these genes has been shown to induce an autoimmunity state that strongly affects plant growth and yield. For this reason, these genes have to be tightly regulated in their expression and activity, and several regulatory mechanisms are described here that tune their gene expression and protein levels. This gene family is subjected to rapid evolution, and to maintain diversity at NLRs, a plethora of genetic mechanisms have been identified as sources of variation. Interestingly, regulation of gene expression and evolution of this gene family are two strictly interconnected aspects. Indeed, some examples have been reported in which mechanisms of gene expression regulation have roles in promotion of the evolution of this gene family. Moreover, co-evolution of the NLR gene family and other gene families devoted to their control has been recently demonstrated, as in the case of miRNAs.

Mapping QTL for Root and Shoot Morphological Traits in a Durum Wheat × T. dicoccum Segregating Population at Seedling Stage

A segregating population of 136 recombinant inbred lines derived from a cross between the durum wheat cv. "Simeto" and the T. dicoccum accession "Molise Colli" was grown in soil and evaluated for a number of shoot and root morphological traits. A total of 17 quantitative trait loci (QTL) were identified for shoot dry weight, number of culms, and plant height and for root dry weight, volume, length, surface area, and number of forks and tips, on chromosomes 1B, 2A, 3A, 4B, 5B, 6A, 6B, and 7B. LODs were 2.1 to 21.6, with percent of explained phenotypic variability between 0.07 and 52. Three QTL were mapped to chromosome 4B, one of which corresponds to the Rht-B1 locus and has a large impact on both shoot and root traits (LOD 21.6). Other QTL that have specific effects on root morphological traits were also identified. Moreover, meta-QTL analysis was performed to compare the QTL identified in the "Simeto" × "Molise Colli" segregating population with those described in previous studies in wheat, with three novel QTL defined. Due to the complexity of phenotyping for root traits, further studies will be helpful to validate these regions as targets for breeding programs for optimization of root function for field performance.

DArT Markers Effectively Target Gene Space in the Rye Genome

Large genome size and complexity hamper considerably the genomics research in relevant species. Rye (Secale cereale L.) has one of the largest genomes among cereal crops and repetitive sequences account for over 90% of its length. Diversity Arrays Technology is a high-throughput genotyping method, in which a preferential sampling of gene-rich regions is achieved through the use of methylation sensitive restriction enzymes. We obtained sequences of 6,177 rye DArT markers and following a redundancy analysis assembled them into 3,737 non-redundant sequences, which were then used in homology searches against five Pooideae sequence sets. In total 515 DArT sequences could be incorporated into publicly available rye genome zippers providing a starting point for the integration of DArT- and transcript-based genomics resources in rye. Using Blast2Go pipeline we attributed putative gene functions to 1101 (29.4%) of the non-redundant DArT marker sequences, including 132 sequences with putative disease resistance-related functions, which were found to be preferentially located in the 4RL and 6RL chromosomes. Comparative analysis based on the DArT sequences revealed obvious inconsistencies between two recently published high density consensus maps of rye. Furthermore we demonstrated that DArT marker sequences can be a source of SSR polymorphisms. Obtained data demonstrate that DArT markers effectively target gene space in the large, complex, and repetitive rye genome. Through the annotation of putative gene functions and the alignment of DArT sequences relative to reference genomes we obtained information, that will complement the results of the studies, where DArT genotyping was deployed, by simplifying the gene ontology and microcolinearity based identification of candidate genes.

Identification of New Resistance Loci to African Stem Rust Race TTKSK in Tetraploid Wheats Based on Linkage and Genome-Wide Association Mapping

Stem rust, caused by Puccinia graminis Pers. f. sp. tritici Eriks. and E. Henn. (Pgt), is one of the most destructive diseases of wheat. Races of the pathogen in the "Ug99 lineage" are of international concern due to their virulence for widely used stem rust resistance genes and their spread throughout Africa. Disease resistant cultivars provide one of the best means for controlling stem rust. To identify quantitative trait loci (QTL) conferring resistance to African stem rust race TTKSK at the seedling stage, we evaluated an association mapping (AM) panel consisting of 230 tetraploid wheat accessions under greenhouse conditions. A high level of phenotypic variation was observed in response to race TTKSK in the AM panel, allowing for genome-wide association mapping of resistance QTL in wild, landrace, and cultivated tetraploid wheats. Thirty-five resistance QTL were identified on all chromosomes, and seventeen are of particular interest as identified by multiple associations. Many of the identified resistance loci were coincident with previously identified rust resistance genes; however, nine on chromosomes 1AL, 2AL, 4AL, 5BL, and 7BS may be novel. To validate AM results, a biparental population of 146 recombinant inbred lines was also considered, which derived from a cross between the resistant cultivar "Cirillo" and susceptible "Neodur." The stem rust resistance of Cirillo was conferred by a single gene on the distal region of chromosome arm 6AL in an interval map coincident with the resistance gene Sr13, and confirmed one of the resistance loci identified by AM. A search for candidate resistance genes was carried out in the regions where QTL were identified, and many of them corresponded to NBS-LRR genes and protein kinases with LRR domains. The results obtained in the present study are of great interest as a high level of genetic variability for resistance to race TTKSK was described in a germplasm panel comprising most of the tetraploid wheat sub-species.

Genome wide identification of C1-2i zinc finger proteins and their response to abiotic stress in hexaploid wheat

The C1-2i wheat Q-type C2H2 zinc finger protein (ZFP) transcription factor subclass has been reported to play important roles in plant stress responses. This subclass of ZFPs has not been studied in hexaploid wheat (Triticum aestivum) and we aimed to identify all members of this subclass and evaluate their responses to different abiotic stresses causing oxidative stress. Exploiting the recently published wheat draft genome sequence, we identified 53 members (including homoeologs from A, B and D genomes) of the C1-2i wheat Q-type C2H2 ZFPs (TaZFPs) representing 21 genes. Evolution analysis revealed that 9 TaZFPs members are directly inherited from the parents Triticum urartu and Aegilops tauschii, while 15 diverged through neoploidization events. This TaZFP subclass is responsive to the oxidative stress generator H2O2 and to high light, drought stress and flooding. Most TaZFPs are responsive to H2O2 (37/53), high light (44/53), flooding (31/53) or drought (37/53); 32 TaZFPs were up-regulated by at least 3 stresses and 16 were responsive to all stresses tested. A large number of these TaZFPs were physically mapped on different wheat draft genome sequences with known markers useful for QTL mapping. Our results show that the C1-2i subclass of TaZFPs is associated with responses to different abiotic stresses and that most TaZFPs (30/53 or 57 %) are located on group 5 chromosomes known to be involved in environment adaptation. Detailed characterization of these novel wheat TaZFPs and their association to QTL or eQTL may help to design wheat cultivars with improved tolerance to abiotic stress.

Radiation hybrid maps of the D-genome of Aegilops tauschii and their application in sequence assembly of large and complex plant genomes

BACKGROUND

The large and complex genome of bread wheat (Triticum aestivum L., ~17 Gb) requires high resolution genome maps with saturated marker scaffolds to anchor and orient BAC contigs/ sequence scaffolds for whole genome assembly. Radiation hybrid (RH) mapping has proven to be an excellent tool for the development of such maps for it offers much higher and more uniform marker resolution across the length of the chromosome compared to genetic mapping and does not require marker polymorphism per se, as it is based on presence (retention) vs. absence (deletion) marker assay.

METHODS

In this study, a 178 line RH panel was genotyped with SSRs and DArT markers to develop the first high resolution RH maps of the entire D-genome of Ae. tauschii accession AL8/78. To confirm map order accuracy, the AL8/78-RH maps were compared with:1) a DArT consensus genetic map constructed using more than 100 bi-parental populations, 2) a RH map of the D-genome of reference hexaploid wheat 'Chinese Spring', and 3) two SNP-based genetic maps, one with anchored D-genome BAC contigs and another with anchored D-genome sequence scaffolds. Using marker sequences, the RH maps were also anchored with a BAC contig based physical map and draft sequence of the D-genome of Ae. tauschii.

RESULTS

A total of 609 markers were mapped to 503 unique positions on the seven D-genome chromosomes, with a total map length of 14,706.7 cR. The average distance between any two marker loci was 29.2 cR which corresponds to 2.1 cM or 9.8 Mb. The average mapping resolution across the D-genome was estimated to be 0.34 Mb (Mb/cR) or 0.07 cM (cM/cR). The RH maps showed almost perfect agreement with several published maps with regard to chromosome assignments of markers. The mean rank correlations between the position of markers on AL8/78 maps and the four published maps, ranged from 0.75 to 0.92, suggesting a good agreement in marker order. With 609 mapped markers, a total of 2481 deletions for the whole D-genome were detected with an average deletion size of 42.0 Mb. A total of 520 markers were anchored to 216 Ae. tauschii sequence scaffolds, 116 of which were not anchored earlier to the D-genome.

CONCLUSION

This study reports the development of first high resolution RH maps for the D-genome of Ae. tauschii accession AL8/78, which were then used for the anchoring of unassigned sequence scaffolds. This study demonstrates how RH mapping, which offered high and uniform resolution across the length of the chromosome, can facilitate the complete sequence assembly of the large and complex plant genomes.

Identification of genes induced by Fusarium graminearum inoculation in the resistant durum wheat line Langdon(Dic-3A)10 and the susceptible parental line Langdon

The wheat recombinant chromosome inbred line LDN(Dic-3A)10, obtained through introgression of a Triticum dicoccoides disomic chromosome 3A fragment into Triticum turgidum spp. durum var. Langdon, is resistant to fusarium head blight (FHB) caused by Fusarium graminearum. To identify genes involved in FHB resistance, we used a cDNA-AFLP approach to compare gene expression between LDN(Dic-3A)10 and the susceptible parental line LDN at different time points post-inoculation. In total, 85 out of the ∼ 500 transcript-derived fragments (TDFs) were found to be differentially expressed: 36 and 19% were upregulated in LDN(Dic-3A)10 and LDN, respectively, whereas 45% were induced in both genotypes. Several of the cloned TDFs showed similarity to proteins involved in specific recognition of plant pathogens or associated with early responses to infection. Some TDFs specific to the inoculation response did not show similarity to characterized proteins. The availability of T. aestivum genome sequences allowed the in silico mapping of 28 TDFs and the acquirement of the corresponding gene sequences and, in some cases, their regulatory regions. Analysis of promoter regions revealed the potential existence of shared transcription regulation mechanisms. For instance, three TDF-associated genes contained binding sites for WRKY transcription factors, which have been implicated in the regulation of genes associated with pathogen defense, and three for abscisic acid-responsive element (ABRE). Collectively, our results revealed specific pathogen recognition in the interactions of LDN and LDN(Dic-3A)10 with F. graminearum. Such recognition leads to changes in the expression of several transcripts, attributable to the presence of the wheat QTL Qfhs.ndsu-3AS.

Genetic analysis of root morphological traits in wheat

Traits related to root architecture are of great importance for yield performance of crop species, although they remain poorly understood. The present study is aimed at identifying the genomic regions involved in the control of root morphological traits in durum wheat (Triticum durum Desf.). A set of 123 recombinant inbred lines derived from the durum wheat cross of cvs. 'Creso' × 'Pedroso' were grown hydroponically to two growth stages, and were phenotypically evaluated for a number of root traits. In addition, meta-(M)QTL analysis was performed that considered the results of other root traits studies in wheat, to compare with the 'Creso' × 'Pedroso' cross and to increase the QTL detection power. Eight quantitative trait loci (QTL) for traits related to root morphology were identified on chromosomes 1A, 1B, 2A, 3A, 6A and 6B in the 'Creso' × 'Pedroso' segregating population. Twenty-two MQTL that comprised from two to six individual QTL that had widely varying confidence intervals were found on 14 chromosomes. The data from the present study provide a detailed analysis of the genetic basis of morphological root traits in wheat. This study of the 'Creso' × 'Pedroso' durum-wheat population has revealed some QTL that had not been previously identified.

Genome-wide association mapping for seedling and adult plant resistance to stripe rust in synthetic hexaploid wheat

Use of genetic diversity from related wild and domesticated species has made a significant contribution to improving wheat productivity. Synthetic hexaploid wheats (SHWs) exhibit natural genetic variation for resistance and/or tolerance to biotic and abiotic stresses. Stripe rust caused by (Puccinia striiformis f. sp. tritici; Pst), is an important disease of wheat worldwide. To characterise loci conferring resistance to stripe rust in SHWs, we conducted a genome-wide association study (GWAS) with a panel of 181 SHWs using the wheat 9 K SNP iSelect array. The SHWs were evaluated for their response to the prevailing races of Pst at the seedling and adult plant stages, the latter in replicated field trials at two sites in Ethiopia in 2011. About 28% of the SHWs exhibited immunity at the seedling stage while 56% and 83% were resistant to Pst at the adult plant stage at Meraro and Arsi Robe, respectively. A total of 27 SNPs in nine genomic regions (1 BS, 2 AS, 2 BL, 3 BL, 3 DL, 5A, 5 BL, 6DS and 7A) were linked with resistance to Pst at the seedling stage, while 38 SNPs on 18 genomic regions were associated with resistance at the adult plant stage. Six genomic regions were commonly detected at both locations using a mixed linear model corrected for population structure, kinship relatedness and adjusted for false discovery rate (FDR). The loci on chromosome regions 1 AS, 3 DL, 6 DS and 7 AL appeared to be novel QTL; our results confirm that resynthesized wheat involving its progenitor species is a rich source of new stripe (yellow) rust resistance that may be useful in choosing SHWs and incorporating diverse yellow rust (YR) resistance loci into locally adapted wheat cultivars.

Genome-wide association for grain morphology in synthetic hexaploid wheats using digital imaging analysis

BACKGROUND

Grain size and shape greatly influence grain weight which ultimately enhances grain yield in wheat. Digital imaging (DI) based phenomic characterization can capture the three dimensional variation in grain size and shape than has hitherto been possible. In this study, we report the results from using digital imaging of grain size and shape to understand the relationship among different components of this trait, their contribution to enhance grain weight, and to identify genomic regions (QTLs) controlling grain morphology using genome wide association mapping with high density diversity array technology (DArT) and allele-specific markers.

RESULTS

Significant positive correlations were observed between grain weight and grain size measurements such as grain length (r = 0.43), width, thickness (r = 0.64) and factor from density (FFD) (r = 0.69). A total of 231 synthetic hexaploid wheats (SHWs) were grouped into five different sub-clusters by Bayesian structure analysis using unlinked DArT markers. Linkage disequilibrium (LD) decay was observed among DArT loci > 10 cM distance and approximately 28% marker pairs were in significant LD. In total, 197 loci over 60 chromosomal regions and 79 loci over 31 chromosomal regions were associated with grain morphology by genome wide analysis using general linear model (GLM) and mixed linear model (MLM) approaches, respectively. They were mainly distributed on homoeologous group 2, 3, 6 and 7 chromosomes. Twenty eight marker-trait associations (MTAs) on the D genome chromosomes 2D, 3D and 6D may carry novel alleles with potential to enhance grain weight due to the use of untapped wild accessions of Aegilops tauschii. Statistical simulations showed that favorable alleles for thousand kernel weight (TKW), grain length, width and thickness have additive genetic effects. Allelic variations for known genes controlling grain size and weight, viz. TaCwi-2A, TaSus-2B, TaCKX6-3D and TaGw2-6A, were also associated with TKW, grain width and thickness. In silico functional analysis predicted a range of biological functions for 32 DArT loci and receptor like kinase, known to affect plant development, appeared to be common protein family encoded by several loci responsible for grain size and shape.

CONCLUSION

Conclusively, we demonstrated the application and integration of multiple approaches including high throughput phenotyping using DI, genome wide association studies (GWAS) and in silico functional analysis of candidate loci to analyze target traits, and identify candidate genomic regions underlying these traits. These approaches provided great opportunity to understand the breeding value of SHWs for improving grain weight and enhanced our deep understanding on molecular genetics of grain weight in wheat.

Genetic diversity and population structure analysis of European hexaploid bread wheat (Triticum aestivum L.) varieties

Progress in plant breeding is facilitated by accurate information about genetic structure and diversity. Here, Diversity Array Technology (DArT) was used to characterize a population of 94 bread wheat (Triticum aestivum L.) varieties of mainly European origin. In total, 1,849 of 7,000 tested markers were polymorphic and could be used for population structure analysis. Two major subgroups of wheat varieties, GrI and GrII, were identified using the program STRUCTURE, and confirmed by principal component analysis (PCA). These subgroups were largely separated according to origin; GrI comprised varieties from Southern and Eastern Europe, whereas GrII contained mostly modern varieties from Western and Northern Europe. A large proportion of the markers contributing most to the genetic separation of the subgroups were located on chromosome 2D near the Reduced height 8 (Rht8) locus, and PCR-based genotyping suggested that breeding for the Rht8 allele had a major impact on subgroup separation. Consistently, analysis of linkage disequilibrium (LD) suggested that different selective pressures had acted on chromosome 2D in the two subgroups. Our data provides an overview of the allele composition of bread wheat varieties anchored to DArT markers, which will facilitate targeted combination of alleles following DArT-based QTL studies. In addition, the genetic diversity and distance data combined with specific Rht8 genotypes can now be used by breeders to guide selection of crossing parents.

Genome-wide association mapping of yield and yield components of spring wheat under contrasting moisture regimes

A stable QTL that may be used in marker-assisted selection in wheat breeding programs was detected for yield, yield components and drought tolerance-related traits in spring wheat association mapping panel. Genome-wide association mapping has become a widespread method of quantitative trait locus (QTL) identification for many crop plants including wheat (Triticum aestivum L.). Its benefit over traditional bi-parental mapping approaches depends on the extent of linkage disequilibrium in the mapping population. The objectives of this study were to determine linkage disequilibrium decay rate and population structure in a spring wheat association mapping panel (n = 285-294) and to identify markers associated with yield and yield components, morphological, phenological, and drought tolerance-related traits. The study was conducted under fully irrigated and rain-fed conditions at Greeley, CO, USA and Melkassa, Ethiopia in 2010 and 2011 (five total environments). Genotypic data were generated using diversity array technology markers. Linkage disequilibrium decay rate extended over a longer genetic distance for the D genome (6.8 cM) than for the A and B genomes (1.7 and 2.0 cM, respectively). Seven subpopulations were identified with population structure analysis. A stable QTL was detected for grain yield on chromosome 2DS both under irrigated and rain-fed conditions. A multi-trait region significant for yield and yield components was found on chromosome 5B. Grain yield QTL on chromosome 1BS co-localized with harvest index QTL. Vegetation indices shared QTL with harvest index on chromosome 1AL and 5A. After validation in relevant genetic backgrounds and environments, QTL detected in this study for yield, yield components and drought tolerance-related traits may be used in marker-assisted selection in wheat breeding programs.

Description of durum wheat linkage map and comparative sequence analysis of wheat mapped DArT markers with rice and Brachypodium genomes

Background

The importance of wheat to the world economy, together with progresses in high-throughput next-generation DNA sequencing, have accelerated initiatives of genetic research for wheat improvement. The availability of high density linkage maps is crucial to identify genotype-phenotype associations, but also for anchoring BAC contigs to genetic maps, a strategy followed for sequencing the wheat genome.

Results

Here we report a genetic linkage map in a durum wheat segregating population and the study of mapped DArT markers. The linkage map consists of 126 gSSR, 31 EST-SSR and 351 DArT markers distributed in 24 linkage groups for a total length of 1,272 cM. Through bioinformatic approaches we have analysed 327 DArT clones to reveal their redundancy, syntenic and functional aspects. The DNA sequences of 174 DArT markers were assembled into a non-redundant set of 60 marker clusters. This explained the generation of clusters in very small chromosome regions across genomes. Of these DArT markers, 61 showed highly significant (Expectation < E-10) BLAST similarity to gene sequences in public databases of model species such as Brachypodium and rice. Based on sequence alignments, the analysis revealed a mosaic gene conservation, with 54 and 72 genes present in rice and Brachypodium species, respectively.

Conclusions

In the present manuscript we provide a detailed DArT markers characterization and the basis for future efforts in durum wheat map comparing.

Different stress responsive strategies to drought and heat in two durum wheat cultivars with contrasting water use efficiency

Background

Durum wheat often faces water scarcity and high temperatures, two events that usually occur simultaneously in the fields. Here we report on the stress responsive strategy of two durum wheat cultivars, characterized by different water use efficiency, subjected to drought, heat and a combination of both stresses.

Results

The cv Ofanto (lower water use efficiency) activated a large set of well-known drought-related genes after drought treatment, while Cappelli (higher water use efficiency) showed the constitutive expression of several genes induced by drought in Ofanto and a modulation of a limited number of genes in response to stress. At molecular level the two cvs differed for the activation of molecular messengers, genes involved in the regulation of chromatin condensation, nuclear speckles and stomatal closure. Noteworthy, the heat response in Cappelli involved also the up-regulation of genes belonging to fatty acid β-oxidation pathway, glyoxylate cycle and senescence, suggesting an early activation of senescence in this cv. A gene of unknown function having the greatest expression difference between the two cultivars was selected and used for expression QTL analysis, the corresponding QTL was mapped on chromosome 6B.

Conclusion

Ofanto and Cappelli are characterized by two opposite stress-responsive strategies. In Ofanto the combination of drought and heat stress led to an increased number of modulated genes, exceeding the simple cumulative effects of the two single stresses, whereas in Cappelli the same treatment triggered a number of differentially expressed genes lower than those altered in response to heat stress alone. This work provides clear evidences that the genetic system based on Cappelli and Ofanto represents an ideal tool for the genetic dissection of the molecular response to drought and other abiotic stresses.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-821) contains supplementary material, which is available to authorized users.

Genetic basis of qualitative and quantitative resistance to powdery mildew in wheat: from consensus regions to candidate genes

Background

Powdery mildew (Blumeria graminis f. sp. tritici) is one of the most damaging diseases of wheat. The objective of this study was to identify the wheat genomic regions that are involved in the control of powdery mildew resistance through a quantitative trait loci (QTL) meta-analysis approach. This meta-analysis allows the use of collected QTL data from different published studies to obtain consensus QTL across different genetic backgrounds, thus providing a better definition of the regions responsible for the trait, and the possibility to obtain molecular markers that will be suitable for marker-assisted selection.

Results

Five QTL for resistance to powdery mildew were identified under field conditions in the durum-wheat segregating population Creso × Pedroso. An integrated map was developed for the projection of resistance genes/ alleles and the QTL from the present study and the literature, and to investigate their distribution in the wheat genome. Molecular markers that correspond to candidate genes for plant responses to pathogens were also projected onto the map, particularly considering NBS-LRR and receptor-like protein kinases. More than 80 independent QTL and 51 resistance genes from 62 different mapping populations were projected onto the consensus map using the Biomercator statistical software. Twenty-four MQTL that comprised 2–6 initial QTL that had widely varying confidence intervals were found on 15 chromosomes. The co-location of the resistance QTL and genes was investigated. Moreover, from analysis of the sequences of DArT markers, 28 DArT clones mapped on wheat chromosomes have been shown to be associated with the NBS-LRR genes and positioned in the same regions as the MQTL for powdery mildew resistance.

Conclusions

The results from the present study provide a detailed analysis of the genetic basis of resistance to powdery mildew in wheat. The study of the Creso × Pedroso durum-wheat population has revealed some QTL that had not been previously identified. Furthermore, the analysis of the co-localization of resistance loci and functional markers provides a large list of candidate genes and opens up a new perspective for the fine mapping and isolation of resistance genes, and for the marker-assisted improvement of resistance in wheat.

Genetic Diversity and Population Structure of Tetraploid Wheats (Triticum turgidum L.) Estimated by SSR, DArT and Pedigree Data

Levels of genetic diversity and population genetic structure of a collection of 230 accessions of seven tetraploid Triticum turgidum L. subspecies were investigated using six morphological, nine seed storage protein loci, 26 SSRs and 970 DArT markers. The genetic diversity of the morphological traits and seed storage proteins was always lower in the durum wheat compared to the wild and domesticated emmer. Using Bayesian clustering (K = 2), both of the sets of molecular markers distinguished the durum wheat cultivars from the other tetraploid subspecies, and two distinct subgroups were detected within the durum wheat subspecies, which is in agreement with their origin and year of release. The genetic diversity of morphological traits and seed storage proteins was always lower in the improved durum cultivars registered after 1990, than in the intermediate and older ones. This marked effect on diversity was not observed for molecular markers, where there was only a weak reduction. At K >2, the SSR markers showed a greater degree of resolution than for DArT, with their identification of a greater number of groups within each subspecies. Analysis of DArT marker differentiation between the wheat subspecies indicated outlier loci that are potentially linked to genes controlling some important agronomic traits. Among the 211 loci identified under selection, 109 markers were recently mapped, and some of these markers were clustered into specific regions on chromosome arms 2BL, 3BS and 4AL, where several genes/quantitative trait loci (QTLs) are involved in the domestication of tetraploid wheats, such as the tenacious glumes (Tg) and brittle rachis (Br) characteristics. On the basis of these results, it can be assumed that the population structure of the tetraploid wheat collection partially reflects the evolutionary history of Triticum turgidum L. subspecies and the genetic potential of landraces and wild accessions for the detection of unexplored alleles.

Plant nucleotide binding site-leucine-rich repeat (NBS-LRR) genes: active guardians in host defense responses

The most represented group of resistance genes are those of the nucleotide binding site-leucine-rich repeat (NBS-LRR) class. These genes are very numerous in the plant genome, and they often occur in clusters at specific loci following gene duplication and amplification events. To date, hundreds of resistance genes and relatively few quantitative trait loci for plant resistance to pathogens have been mapped in different species, with some also cloned. When these NBS-LRR genes have been physically or genetically mapped, many cases have shown co-localization between resistance loci and NBS-LRR genes. This has allowed the identification of candidate genes for resistance, and the development of molecular markers linked to R genes. This review is focused on recent genomics studies that have described the abundance, distribution and evolution of NBS-LRR genes in plant genomes. Furthermore, in terms of their expression, NBS-LRR genes are under fine regulation by cis- and trans-acting elements. Recent findings have provided insights into the roles of alternative splicing, the ubiquitin/ proteasome system, and miRNAs and secondary siRNAs in the regulation of NBS-LRR gene expression at the post-transcriptional, post-translational and epigenetic levels. The possibility to use this knowledge for genetic improvement of plant resistance to pathogens is discussed.