Genetic diversity and population structure of modern wheat (Triticum aestivum L.) cultivars in Henan Province of China based on SNP markers.

BACKGROUND: Henan is the province with the greatest wheat production in China. Although more than 100 cultivars are used for production, many cultivars are still insufficient in quality, disease resistance, adaptability and yield potential. To overcome these limitations, it is necessary to constantly breed new cultivars to maintain the continuous and stable growth of wheat yield and quality. To improve breeding efficiency, it is important to evaluate the genetic diversity and population genetic structure of its cultivars. However, there are no such reports from Henan Province. Therefore, in this study, single nucleotide polymorphism (SNP) markers were used to study the population genetic structure and genetic diversity of 243 wheat cultivars included in a comparative test of wheat varieties in Henan Province, aiming to provide a reference for the utilization of backbone parents and the selection of hybrid combinations in the genetic improvement of wheat cultivars. RESULTS: In this study, 243 wheat cultivars from Henan Province of China were genotyped by the Affymetrix Axiom Wheat660K SNP chip, and 21 characteristics were investigated. The cultivars were divided into ten subgroups; each subgroup had distinct characteristics and unique utilization value. Furthermore, based on principal component analysis, Zhoumai cultivars were the main hybrid parents, followed by Aikang 58, high-quality cultivars, and Shandong cultivars. Genetic diversity analysis showed that 61.3% of SNPs had a high degree of genetic differentiation, whereas 33.4% showed a moderate degree. The nucleotide diversity of subgenome B was relatively high, with an average π value of 3.91E-5; the nucleotide diversity of subgenome D was the lowest, with an average π value of 2.44E-5. CONCLUSION: The parents used in wheat cross-breeding in Henan Province are similar, with a relatively homogeneous genetic background and low genetic diversity. These results will not only contribute to the objective evaluation and utilization of the tested cultivars but also provide insights into the current conditions and existing challenges of wheat cultivar breeding in Henan Province, thereby facilitating the scientific formulation of breeding objectives and strategies to improve breeding efficiency.

A papain-like cysteine protease-released small signal peptide confers wheat resistance to wheat yellow mosaic virus.

Wheat yellow mosaic virus (WYMV), a soil-borne pathogen, poses a serious threat to global wheat production. Here, we identify a WYMV resistance gene, TaRD21A, that belongs to the papain-like cysteine protease family. Through genetic manipulation of TaRD21A expression, we establish its positive role in the regulation of wheat to WYMV resistance. Furthermore, our investigation shows that the TaRD21A-mediated plant antiviral response relies on the release of a small peptide catalyzed by TaRD21A protease activity. To counteract wheat resistance, WYMV-encoded nuclear inclusion protease-a (NIa) suppress TaRD21A activity to promote virus infection. In resistant cultivars, a natural variant of TaRD21A features a glycine-to-threonine substitution and this substitution enables the phosphorylation of threonine, thereby weakening the interaction between NIa and TaRD21A, reinforcing wheat resistance against WYMV. Our study not only unveils a WYMV resistance gene but also offers insights into the intricate mechanisms underpinning resistance against WYMV.

Mapping QTL conferring flag leaf senescence in durum wheat cultivars.

Flag leaf senescence is a critical factor affecting the yield and quality of wheat. The aim of this study was to identify QTLs associated with flag leaf senescence in an F10 recombinant inbred line population derived from durum wheats UC1113 and Kofa. Bulked segregant analysis using the wheat 660K SNP array identified 3225 SNPs between extreme-phenotype bulks, and the differential SNPs were mainly clustered on chromosomes 1A, 1B, 3B, 5A, 5B, and 7A. BSR-Seq indicated that the significant SNPs were mainly located in two intervals of 354.0-389.0 Mb and 8.0-15.0 Mb on 1B and 3B, respectively. Based on the distribution of significant SNPs on chromosomes 1B and 3B, a total of 109 insertion/deletion (InDel) markers were developed, and 8 of them were finally used to map QTL in UC1113/Kofa population for flag leaf senescence. Inclusive composite interval mapping identified two major QTL in marker intervals Mar2005-Mar2116 and Mar207-Mar289, explaining 14.2-15.4% and 31.4-68.6% of the phenotypic variances across environments, respectively. Using BSR-Seq, gene expression and sequence analysis, the TraesCS1B02G211600 and TraesCS3B02G023000 were identified as candidate senescence-associated genes. This study has potential to be used in cloning key genes for flag leaf senescence and provides available molecular markers for genotyping and marker-assisted selection breeding. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-023-01410-3.

Mapping of quantitative trait loci for leaf rust resistance in the wheat population 'Xinmai 26/Zhoumai 22'.

Leaf rust, caused by the fungal pathogen Puccinia triticina (Pt), is one of the major and dangerous diseases of wheat, and has caused serious yield loss of wheat worldwide. Here, we investigated adult-plant resistance (APR) to leaf rust in a recombinant inbred line (RIL) population derived from 'Xinmai 26' and 'Zhoumai 22' over 3 years. Linkage mapping for APR to leaf rust revealed four quantitative trait loci (QTL) in this RIL population. Two QTL, QLr.hnau-2BS and QLr.hnau-3BS were contributed by 'Zhoumai22', whereas QLr.hnau-2DS and QLr.hnau-5AL were contributed by 'Xinmai 26'. The QLr.hnau-2BS covering a race-specific resistance gene Lr13 showed the most stable APR to leaf rust. Overexpression of Lr13 significantly increased APR to leaf rust. Interestingly, we found that a CNL(coiled coil-nucleotide-binding site-leucine-rich repeat)-like gene, TaCN, in QLr.hnau-2BS completely co-segregated with leaf rust resistance. The resistant haplotype TaCN-R possessed half the sequence of the coiled-coil domain of TaCN protein. Lr13 strongly interacted with TaCN-R, but did not interact with the full-length TaCN (TaCN-S). In addition, TaCN-R was significantly induced after Pt inoculation and changed the sub-cellular localization of Lr13 after interaction. Therefore, we hypothesized that TaCN-R mediated leaf rust resistance possibly by interacting with Lr13. This study provides important QTL for APR to leaf rust, and new insights into understanding how a CNL gene modulates disease resistance in common wheat.

HSP90.2 modulates 2Q2-mediated wheat resistance against powdery mildew.

Wheat (Triticum aestivum L.) is a critical food crop feeding the world, but pathogens threaten its production. Wheat Heat Shock Protein 90.2 (HSP90.2) is a pathogen-inducible molecular chaperone folding nascent preproteins. Here, we used wheat HSP90.2 to isolate clients regulated at the posttranslational level. Tetraploid wheat hsp90.2 knockout mutant was susceptible to powdery mildew, while the HSP90.2 overexpression line was resistant, suggesting that HSP90.2 was essential for wheat resistance against powdery mildew. We next isolated 1500 clients of HSP90.2, which contained a wide variety of clients with different biological classifications. We utilized 2Q2, a nucleotide-binding leucine repeat-rich protein, as a model to investigate the potential of HSP90.2 interactome in fungal resistance. The transgenic line co-suppressing 2Q2 was more susceptible to powdery mildew, suggesting 2Q2 as a novel Pm-resistant gene. The 2Q2 protein resided in chloroplasts, and HSP90.2 played a critical role in the accumulation of 2Q2 in thylakoids. Our data provided over 1500 HSP90.2 clients with a potential regulation at the protein folding process and contributed a nontypical approach to isolate pathogenesis-related proteins.

Genetic and molecular basis of carotenoid metabolism in cereals.

Carotenoids are vital pigments for higher plants and play a crucial function in photosynthesis and photoprotection. Carotenoids are precursors of vitamin A synthesis and contribute to human nutrition and health. However, cereal grain endosperm contains a minor carotenoid measure and a scarce supply of provitamin A content. Therefore, improving the carotenoids in cereal grain is of major importance. Carotenoid content is governed by multiple candidate genes with their additive effects. Studies on genes related to carotenoid metabolism in cereals would increase the knowledge of potential metabolic steps of carotenoids and enhance the quality of crop plants. Recognizing the metabolism and carotenoid accumulation in various staple cereal crops over the last few decades has broadened our perspective on the interdisciplinary regulation of carotenogenesis. Meanwhile, the amelioration in metabolic engineering approaches has been exploited to step up the level of carotenoid and valuable industrial metabolites in many crops, but wheat is still considerable in this matter. In this study, we present a comprehensive overview of the consequences of biosynthetic and catabolic genes on carotenoid biosynthesis, current improvements in regulatory disciplines of carotenogenesis, and metabolic engineering of carotenoids. A panoptic and deeper understanding of the regulatory mechanisms of carotenoid metabolism and genetic manipulation (genome selection and gene editing) will be useful in improving the carotenoid content of cereals.

HSP90.2 promotes CO2 assimilation rate, grain weight and yield in wheat.

Wheat fixes CO2 by photosynthesis into kernels to nourish humankind. Improving the photosynthesis rate is a major driving force in assimilating atmospheric CO2 and guaranteeing food supply for human beings. Strategies for achieving the above goal need to be improved. Here, we report the cloning and mechanism of CO2 ASSIMILATION RATE AND KERNEL-ENHANCED 1 (CAKE1) from durum wheat (Triticum turgidum L. var. durum). The cake1 mutant displayed a lower photosynthesis rate with smaller grains. Genetic studies identified CAKE1 as HSP90.2-B, encoding cytosolic molecular chaperone folding nascent preproteins. The disturbance of HSP90.2 decreased leaf photosynthesis rate, kernel weight (KW) and yield. Nevertheless, HSP90.2 over-expression increased KW. HSP90.2 recruited and was essential for the chloroplast localization of nuclear-encoded photosynthesis units, for example PsbO. Actin microfilaments docked on the chloroplast surface interacted with HSP90.2 as a subcellular track towards chloroplasts. A natural variation in the hexaploid wheat HSP90.2-B promoter increased its transcription activity, enhanced photosynthesis rate and improved KW and yield. Our study illustrated an HSP90.2-Actin complex sorting client preproteins towards chloroplasts to promote CO2 assimilation and crop production. The beneficial haplotype of Hsp90.2 is rare in modern varieties and could be an excellent molecular switch promoting photosynthesis rate to increase yield in future elite wheat varieties.

Identification of TaGL1-B1 gene controlling grain length through regulation of jasmonic acid in common wheat.

Grain length is one of the most important factors in determining wheat yield. Here, a stable QTL for grain length was mapped on chromosome 1B in a F10 recombinant inbred lines (RIL) population, and the gene TaGL1-B1 encoding carotenoid isomerase was identified in a secondary large population through multiple strategies. The genome-wide association study (GWAS) in 243 wheat accessions revealed that the marker for TaGL1-B1 was the most significant among all chromosomes. EMS mutants of TaGL1 possessed significantly reduced grain length, whereas TaGL1-B1-overexpressed lines possessed significantly increased grain length. Moreover, TaGL1-B1 strongly interacted with TaPAP6. TaPAP6-overexpressed lines had significantly increased grain length. Transcriptome analysis suggested that TaPAP6 was possibly involved in the accumulation of JA (jasmonic acid). Consistently, JA content was significantly increased in the TaGL1-B1 and TaPAP6 overexpression lines. Additionally, the role of TaGL1-B1 in regulating carotenoids was verified through QTL mapping, GWAS, EMS mutants and overexpression lines. Notably, overexpression of TaGL1-B1 significantly increased wheat yield in multiple locations. Taken together, overexpression of TaGL1-B1 enhanced grain length, probably through interaction with TaPAP6 to cause the accumulation of JA that improved carotenoid content and photosynthesis, thereby resulted in increased wheat yield. This study provided valuable genes controlling grain length to improve yield and a potential insight into the molecular mechanism of modulating JA-mediated grain size in wheat.

Reduced expression of lipoxygenase genes improves flour processing quality in soft wheat.

Lipoxygenases (Loxs) are dioxygenases that play an important role in plant growth and defense. Loxs affect flour processing quality in common wheat (Triticum aestivum). We conducted a genome-wide association study (GWAS) that identified 306 significant single-nucleotide polymorphisms (SNPs) related to Lox activity in Chinese wheat accessions. Among them, a novel lipoxygenase-encoding (Lpx) gene, TaLpx-B4, was detected on chromosome 3B in a biparental population. Analysis of mutant wheat lines induced using ethyl methanesulfonate confirmed the role of TaLpx-B4 in modulating Lox activity. A phylogenetic tree of various plant Lpx genes indicated the predominance of the 9-Lpx type in common wheat. Further analysis revealed conserved intron number, exon length, and motif number in the TaLpx gene family. GWAS, linkage mapping, and gene annotation collectively showed that 14 out of 29 annotated TaLpx genes played a critical role in regulating Lox activity in the Chinese wheat accessions. Transgenic wheat grains with knockdown of Lpx family genes by RNAi showed significantly lower Lox activity than the wild type. One TaLpx-RNAi line had significantly reduced starch content and dough stability, and thus possessed relatively superior biscuit quality in soft wheat. Further analysis of the transcriptome, lipid components, and other metabolites revealed that knockdown of TaLpx genes significantly increased biscuit quality via changes in unsaturated fatty acid content as well as in starch, sucrose, and galactose metabolism. Our results provide new insights into the role of the TaLpx gene family that will be beneficial in improving soft wheat flour quality.

Identification of Genetic Loci of Black Point in Chinese Common Wheat by Genome-Wide Association Study and Linkage Mapping.

Black point is a common disease in wheat all over the world. The disease could downgrade wheat quality and cause human health problems. In this study, 406 wheat cultivars were used to investigate black point resistance. In the field tests, 20, 65.5, and 14.5% of the tested cultivars were resistant, moderately resistant, and susceptible, respectively, suggesting that improving black point resistance is necessary in Chinese wheat breeding. A genome-wide association study (GWAS) identified 386 single-nucleotide polymorphisms (SNPs) significantly related to black point resistance in the tested wheat cultivars, and they were located on all chromosomes. Linkage mapping in a biparental population identified three quantitative trait loci (QTL) for black point resistance-QBP.hau-3A, QBP.hau-6D, and QBP.hau-7D-with 6.76, 7.79, and 8.84% phenotypic variation explained, respectively. Based on both the GWAS and linkage analyses, QBP.hau-6D covered six significant SNPs from the GWAS, and the position of these SNPs indicated that this QTL is a new locus for black point resistance. This study provides valuable germplasm for breeding wheat cultivars with resistance to black point and information for further understanding of molecular and genetic basis of black point resistance.

The Wheat 660K SNP array demonstrates great potential for marker-assisted selection in polyploid wheat.

The rapid development and application of molecular marker assays have facilitated genomic selection and genome-wide linkage and association studies in wheat breeding. Although PCR-based markers (e.g. simple sequence repeats and functional markers) and genotyping by sequencing have contributed greatly to gene discovery and marker-assisted selection, the release of a more accurate and complete bread wheat reference genome has resulted in the design of single-nucleotide polymorphism (SNP) arrays based on different densities or application targets. Here, we evaluated seven types of wheat SNP arrays in terms of their SNP number, distribution, density, associated genes, heterozygosity and application. The results suggested that the Wheat 660K SNP array contained the highest percentage (99.05%) of genome-specific SNPs with reliable physical positions. SNP density analysis indicated that the SNPs were almost evenly distributed across the whole genome. In addition, 229 266 SNPs in the Wheat 660K SNP array were located in 66 834 annotated gene or promoter intervals. The annotated genes revealed by the Wheat 660K SNP array almost covered all genes revealed by the Wheat 35K (97.44%), 55K (99.73%), 90K (86.9%) and 820K (85.3%) SNP arrays. Therefore, the Wheat 660K SNP array could act as a substitute for other 6 arrays and shows promise for a wide range of possible applications. In summary, the Wheat 660K SNP array is reliable and cost-effective and may be the best choice for targeted genotyping and marker-assisted selection in wheat genetic improvement.

Genome-wide association study of six quality traits reveals the association of the TaRPP13L1 gene with flour colour in Chinese bread wheat.

Flour colour, kernel hardness, grain protein content and wet gluten content are important quality properties that determine end use in bread wheat. Here, a wheat 90K genotyping assay was used for a genome-wide association study (GWAS) of the six quality-related traits in Chinese wheat cultivars in eight environments over four years. A total of 846 significant single nucleotide polymorphisms (SNPs) were identified, explaining approximately 30% of the phenotypic variation on average, and 103 multienvironment-significant SNPs were detected in more than four environments. Quantitative trait loci (QTL) mapping in the biparent population confirmed some important SNP loci. Moreover, it was determined that some important genes were associated with the six quality traits, including some known functional genes and annotated unknown functional genes. Of the annotated unknown functional genes, it was verified that TaRPP13L1 was associated with flour colour. Wheat cultivars or lines with TaRPP13L1-B1a showed extremely significantly higher flour redness and lower yellowness than those with TaRPP13L1-B1b in the Chinese wheat natural population and the doubled haploid (DH) population. Two tetraploid wheat lines with premature stop codons of the TaRPP13L1 gene mutagenized by ethyl methanesulfonate (EMS) showed extremely significantly higher flour redness and lower yellowness than wild type. Our data suggest that the TaRPP13L1 gene plays an important role in modulating wheat flour colour. This study provides useful information for further dissection of the genetic basis of flour colour and also provides valuable genes or genetic loci for marker-assisted selection to improve the process of breeding quality wheat in China.

Correction: Genetic dissection of main and epistatic effects of QTL based on augmented triple test cross design.

[This corrects the article DOI: 10.1371/journal.pone.0189054.].

Genetic dissection of main and epistatic effects of QTL based on augmented triple test cross design.

The use of heterosis has considerably increased the productivity of many crops; however, the biological mechanism underpinning the technique remains elusive. The North Carolina design III (NCIII) and the triple test cross (TTC) are powerful and popular genetic mating design that can be used to decipher the genetic basis of heterosis. However, when using the NCIII design with the present quantitative trait locus (QTL) mapping method, if epistasis exists, the estimated additive or dominant effects are confounded with epistatic effects. Here, we propose a two-step approach to dissect all genetic effects of QTL and digenic interactions on a whole genome without sacrificing statistical power based on an augmented TTC (aTTC) design. Because the aTTC design has more transformation combinations than do the NCIII and TTC designs, it greatly enriches the QTL mapping for studying heterosis. When the basic population comprises recombinant inbred lines (RIL), we can use the same materials in the NCIII design for aTTC-design QTL mapping with transformation combination Z1, Z2, and Z4 to obtain genetic effect of QTL and digenic interactions. Compared with RIL-based TTC design, RIL-based aTTC design saves time, money, and labor for basic population crossed with F1. Several Monte Carlo simulation studies were carried out to confirm the proposed approach; the present genetic parameters could be identified with high statistical power, precision, and calculation speed, even at small sample size or low heritability. Additionally, two elite rice hybrid datasets for nine agronomic traits were estimated for real data analysis. We dissected the genetic effects and calculated the dominance degree of each QTL and digenic interaction. Real mapping results suggested that the dominance degree in Z2 that mainly characterize heterosis showed overdominance and dominance for QTL and digenic interactions. Dominance and overdominance were the major genetic foundations of heterosis in rice.

Genome-wide association study for 13 agronomic traits reveals distribution of superior alleles in bread wheat from the Yellow and Huai Valley of China.

Bread wheat is a leading cereal crop worldwide. Limited amount of superior allele loci restricted the progress of molecular improvement in wheat breeding. Here, we revealed new allelic variation distribution for 13 yield-related traits in series of genome-wide association studies (GWAS) using the wheat 90K genotyping assay, characterized in 163 bread wheat cultivars. Agronomic traits were investigated in 14 environments at three locations over 3 years. After filtering SNP data sets, GWAS using 20 689 high-quality SNPs associated 1769 significant loci that explained, on average, ~20% of the phenotypic variation, both detected already reported loci and new promising genomic regions. Of these, repetitive and pleiotropic SNPs on chromosomes 6AS, 6AL, 6BS, 5BL and 7AS were significantly linked to thousand kernel weight, for example BS00021705_51 on 6BS and wsnp_Ex_c32624_41252144 on 6AS, with phenotypic variation explained (PVE) of ~24%, consistently identified in 12 and 13 of the 14 environments, respectively. Kernel length-related SNPs were mainly identified on chromosomes 7BS, 6AS, 5AL and 5BL. Plant height-related SNPs on chromosomes 4DS, 6DL, 2DS and 1BL were, respectively, identified in more than 11 environments, with averaged PVE of ~55%. Four SNPs were confirmed to be important genetic loci in two RIL populations. Based on repetivity and PVE, a total of 41 SNP loci possibly played the key role in modulating yield-related traits of the cultivars surveyed. Distribution of superior alleles at the 41 SNP loci indicated that superior alleles were getting popular with time and modern cultivars had integrated many superior alleles, especially for peduncle length- and plant height-related superior alleles. However, there were still 19 SNP loci showing less than percentages of 50% in modern cultivars, suggesting they should be paid more attention to improve yield-related traits of cultivars in the Yellow and Huai wheat region. This study could provide useful information for dissection of yield-related traits and valuable genetic loci for marker-assisted selection in Chinese wheat breeding programme.

Genomewide mapping reveals a combination of different genetic effects causing the genetic basis of heterosis in two elite rice hybrids.

North Carolina design III (NCIII) is one of the most powerful and widely used mating designs for understanding the genetic basis of heterosis. However, the quantitative trait mapping (QTL) conducted in previous studies with this design was mainly based on analysis of variance (ANOVA), composite interval or multiple interval mapping methods. These methodologies could not investigate all kinds of genetic effects, especially epistatic effects, simultaneously on the whole genome. In this study, with a statistical method for mapping epistatic QTL associated with heterosis using the recombinant inbred line (RIL)-based NCIII design, we conducted QTL mapping for nine agronomic traits of two elite hybrids to characterize the mode of gene action contributing to heterosis on a whole genomewide scale. In total, 23 main-effect QTL (M-QTL) and 23 digenic interactions in IJ (indica x japonica) hybrids, 11 M-QTL and 82 digenic interactions in II (indica x indica) hybrid QTLs were identified in the present study. The variation explained by individual M-QTL or interactions ranged from 2.3 to 11.0%. The number of digenic interactions and the total variation explained by interactions of each trait were larger than those of M-QTL. The augmented genetic effect ratio of most M-QTL and digenic interactions in (L1 - L2) data of two backcross populations (L1 and L2) showed complete dominance or overdominance, and in (L1 + L2) data showed an additive effect. Our results indicated that the dominance, overdominance and epistatic effect were important in conditioning the genetic basis of heterosis of the two elite hybrids. The relative contributions of the genetic components varied with traits and the genetic basis of the two hybrids was different.

Binary matrix shuffling filter for feature selection in neuronal morphology classification.

A prerequisite to understand neuronal function and characteristic is to classify neuron correctly. The existing classification techniques are usually based on structural characteristic and employ principal component analysis to reduce feature dimension. In this work, we dedicate to classify neurons based on neuronal morphology. A new feature selection method named binary matrix shuffling filter was used in neuronal morphology classification. This method, coupled with support vector machine for implementation, usually selects a small amount of features for easy interpretation. The reserved features are used to build classification models with support vector classification and another two commonly used classifiers. Compared with referred feature selection methods, the binary matrix shuffling filter showed optimal performance and exhibited broad generalization ability in five random replications of neuron datasets. Besides, the binary matrix shuffling filter was able to distinguish each neuron type from other types correctly; for each neuron type, private features were also obtained.

Informative gene selection and direct classification of tumor based on Chi-square test of pairwise gene interactions.

In efforts to discover disease mechanisms and improve clinical diagnosis of tumors, it is useful to mine profiles for informative genes with definite biological meanings and to build robust classifiers with high precision. In this study, we developed a new method for tumor-gene selection, the Chi-square test-based integrated rank gene and direct classifier (χ(2)-IRG-DC). First, we obtained the weighted integrated rank of gene importance from chi-square tests of single and pairwise gene interactions. Then, we sequentially introduced the ranked genes and removed redundant genes by using leave-one-out cross-validation of the chi-square test-based Direct Classifier (χ(2)-DC) within the training set to obtain informative genes. Finally, we determined the accuracy of independent test data by utilizing the genes obtained above with χ(2)-DC. Furthermore, we analyzed the robustness of χ(2)-IRG-DC by comparing the generalization performance of different models, the efficiency of different feature-selection methods, and the accuracy of different classifiers. An independent test of ten multiclass tumor gene-expression datasets showed that χ(2)-IRG-DC could efficiently control overfitting and had higher generalization performance. The informative genes selected by χ(2)-IRG-DC could dramatically improve the independent test precision of other classifiers; meanwhile, the informative genes selected by other feature selection methods also had good performance in χ(2)-DC.