RESUMEN
The erect leaf plays a crucial role in determining plant architecture, with its growth and development regulated by genetic factors. However, there has been a lack of comprehensive studies on the regulatory mechanisms governing wheat lamina joint development, thus failing to meet current breeding demands. In this study, a wheat erect leaf mutant, mths29, induced via fast neutron mutagenesis, was utilized for QTL fine mapping and investigation of lamina joint development. Genetic analysis of segregating populations derived from mths29 and Jimai22 revealed that the erect leaf trait was controlled by a dominant single gene. Using BSR sequencing and map-based cloning techniques, the QTL responsible for the erect leaf trait was mapped to a 1.03 Mb physical region on chromosome 5A. Transcriptome analysis highlighted differential expression of genes associated with cell division and proliferation, as well as several crucial transcription factors and kinases implicated in lamina joint development, particularly in the boundary cells of the preligule zone in mths29. These findings establish a solid foundation for understanding lamina joint development and hold promise for potential improvements in wheat plant architecture.
RESUMEN
As climate change continues to negatively impact our farmlands, abiotic factors like salinity and drought stress increasingly threaten global food security. The development of elite germplasms with resistance to multiple abiotic stresses is essential for breeding climate-resilient wheat cultivars. In this study, we determined that the previously reported salt-tolerant st1 mutant, obtained via spaceflight mutagenesis, may also resist to drought stress at the seedling stage. Moreover, our field trial revealed that yield-related traits including plant height, 1000-grain weight, and spike number per plant were significantly increased in st1 compared to the wild type. An F2 population of 334 individuals derived from a cross between the wild type and st1 displayed a bimodal distribution indicating that st1 plant height is controlled by a single major gene. Our Bulked Segregant Analysis and exome capture sequencing indicate that this gene is located on chromosome 4D. Further genetic linkage and gene sequence analysis suggests that a reverse mutation of Rht2 is putatively responsible for plant height variation in st1. Our genotypic and phenotypic analysis of the F2 population and F3 lines indicate that this reverse mutation significantly increases plant height and thousand grain weight but slightly decreases spike number per plant. Together, these results supply helpful information for the utilization of Rht2 in wheat breeding and provide an important material for breeding environmentally resilient, high-yield wheat varieties.
Asunto(s)
Vuelo Espacial , Triticum , Humanos , Triticum/genética , Fitomejoramiento , Fenotipo , Mutación/genética , Grano ComestibleRESUMEN
Tiller number greatly contributes to grain yield in wheat. Using ethylmethanesulfonate mutagenesis, we previously discovered the oligo-tillering mutant ot1. The tiller number was significantly lower in ot1 than in the corresponding wild type from the early tillering stage until the heading stage. Compared to the wild type, the thousand-grain weight and grain length were increased by 15.41% and 31.44%, respectively, whereas the plant height and spike length were decreased by 26.13% and 37.25%, respectively. Transcriptomic analysis was conducted at the regreening and jointing stages to identify differential expressed genes (DEGs). Functional enrichment analysis with the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) databases showed differential expression of genes associated with ADP binding, transmembrane transport, and transcriptional regulation during tiller development. Differences in tiller number in ot1 led to the upregulation of genes in the strigolactone (SL) and abscisic acid (ABA) pathways. Specifically, the SL biosynthesis genes DWARF (D27), D17, D10, and MORE AXILLARY GROWTH 1 (MAX1) were upregulated by 3.37- to 8.23-fold; the SL signal transduction genes D14 and D53 were upregulated by 1.81- and 1.32-fold, respectively; the ABA biosynthesis genes 9-CIS-EPOXICAROTENOID DIOXIGENASE 3 (NCED3) and NCED5 were upregulated by 1.66- and 3.4-fold, respectively; and SNF1-REGULATED PROTEIN KINASE2 (SnRK2) and PROTEIN PHOSPHATASE 2C (PP2C) genes were upregulated by 1.30- to 4.79-fold. This suggested that the tiller number reduction in ot1 was due to alterations in plant hormone pathways. Genes known to promote tillering growth were upregulated, whereas those known to inhibit tillering growth were downregulated. For example, PIN-FORMED 9 (PIN9), which promotes tiller development, was upregulated by 8.23-fold in ot1; Ideal Plant Architecture 1 (IPA1), which inhibits tiller development, was downregulated by 1.74-fold. There were no significant differences in the expression levels of TILLER NUMBER 1 (TN1) or TEOSINTE BRANCHED 1 (TB1), indicating that the tiller reduction in ot1 was not controlled by known genes. Our findings provide valuable data for subsequent research into the genetic bases and regulatory mechanisms of wheat tillering. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-024-01450-3.
RESUMEN
Homeodomain proteins encoded by BEL1- and KNAT1-type genes are ubiquitously distributed across plant species and play important roles in growth and development, whereby a comprehensive investigation of their molecular interactions and potential functions in wheat is of great significance. In this study, we systematically investigated the phylogenetic relationships, gene structures, conserved domains, and cis-acting elements of 34 TaBEL and 34 TaKNAT genes in the wheat genome. Our analysis revealed these genes evolved under different selective pressures and showed variable transcript levels in different wheat tissues. Subcellular localization analysis further indicated the proteins encoded by these genes were either exclusively located in the nucleus or both in the nucleus and the cytoplasm. Additionally, a comprehensive protein-protein interaction network was constructed with representative genes in which each TaBEL or TaKNAT proteins interact with at least two partners. The evaluation of wheat mutants identified key genes, including TaBEL-5B, TaBEL-4A.4, and TaKNAT6, which are involved in grain-related traits. Finally, haplotype analysis suggests TaKNAT-6B is associated with grain-related traits and is preferentially selected among a large set of wheat accessions. Our study provides important information on BEL1- and KNAT1-type gene families in wheat, and lays the foundation for functional research in the future.
RESUMEN
Hexaploid wheat (Triticum aestivum), a major staple crop, has a remarkably large genome of ~14.4 Gb (containing 106 913 high-confidence [HC] and 159 840 low-confidence [LC] genes in the Chinese Spring v2.1 reference genome), which poses a major challenge for functional genomics studies. To overcome this hurdle, we performed whole-exome sequencing to generate a nearly saturated wheat mutant database containing 18 025 209 mutations induced by ethyl methanesulfonate (EMS), carbon (C)-ion beams, or γ-ray mutagenesis. This database contains an average of 47.1 mutations per kb in each gene-coding sequence: the potential functional mutations were predicted to cover 96.7% of HC genes and 70.5% of LC genes. Comparative analysis of mutations induced by EMS, γ-rays, or C-ion beam irradiation revealed that γ-ray and C-ion beam mutagenesis induced a more diverse array of variations than EMS, including large-fragment deletions, small insertions/deletions, and various non-synonymous single nucleotide polymorphisms. As a test case, we combined mutation analysis with phenotypic screening and rapidly mapped the candidate gene responsible for the phenotype of a yellow-green leaf mutant to a 2.8-Mb chromosomal region. Furthermore, a proof-of-concept reverse genetics study revealed that mutations in gibberellic acid biosynthesis and signalling genes could be associated with negative impacts on plant height. Finally, we built a publically available database of these mutations with the corresponding germplasm (seed stock) repository to facilitate advanced functional genomics studies in wheat for the broad plant research community.
Asunto(s)
Genómica , Triticum , Triticum/genética , Secuenciación del Exoma , Mutación/genética , Mutagénesis , Metanosulfonato de Etilo/farmacología , Genoma de Planta/genéticaRESUMEN
BACKGROUND: Plant height (PH) and spike compactness (SC) are important agronomic traits that affect yield improvement in wheat crops. The identification of the loci or genes responsible for these traits is thus of great importance for marker-assisted selection in wheat breeding. RESULTS: In this study, we used a recombinant inbred line (RIL) population with 139 lines derived from a cross between the mutant Rht8-2 and the local wheat variety NongDa5181 (ND5181) to construct a high-density genetic linkage map by applying the Wheat 40 K Panel. We identified seven stable QTLs for PH (three) and SC (four) in two environments using the RIL population, and found that Rht8-B1 is the causal gene of qPH2B.1 by further genetic mapping, gene cloning and gene editing analyses. Our results also showed that two natural variants from GC to TT in the coding region of Rht8-B1 resulted in an amino acid change from G (ND5181) to V (Rht8-2) at the 175th position, reducing PH by 3.6%~6.2% in the RIL population. Moreover, gene editing analysis suggested that the height of T2 generation in Rht8-B1 edited plants was reduced by 5.6%, and that the impact of Rht8-B1 on PH was significantly lower than Rht8-D1. Additionally, analysis of the distribution of Rht8-B1 in various wheat resources suggested that the Rht8-B1b allele has not been widely utilized in modern wheat breeding. CONCLUSIONS: The combination of Rht8-B1b with other favorable Rht genes might be an alternative approach for developing lodging-resistant crops. Our study provides important information for marker-assisted selection in wheat breeding.
Asunto(s)
Fitomejoramiento , Triticum , Triticum/genética , Mapeo Cromosómico , Fenotipo , Sitios de Carácter Cuantitativo/genéticaRESUMEN
KEY MESSAGE: A minor-effect QTL, Qhd.2AS, that affects heading date in wheat was mapped to a genomic interval of 1.70-Mb on 2AS, and gene analysis indicated that the C2H2-type zinc finger protein gene TraesCS2A02G181200 is the best candidate for Qhd.2AS. Heading date (HD) is a complex quantitative trait that determines the regional adaptability of cereal crops, and identifying the underlying genetic elements with minor effects on HD is important for improving wheat production in diverse environments. In this study, a minor QTL for HD that we named Qhd.2AS was detected on the short arm of chromosome 2A by Bulked Segregant Analysis and validated in a recombinant inbred population. Using a segregating population of 4894 individuals, Qhd.2AS was further delimited to an interval of 0.41 cM, corresponding to a genomic region spanning 1.70 Mb (from 138.87 to 140.57 Mb) that contains 16 high-confidence genes based on IWGSC RefSeq v1.0. Analyses of sequence variations and gene transcription indicated that TraesCS2A02G181200, which encodes a C2H2-type zinc finger protein, is the best candidate gene for Qhd.2AS that influences HD. Screening a TILLING mutant library identified two mutants with premature stop codons in TraesCS2A02G181200, both of which exhibited a delay in HD of 2-4 days. Additionally, variations in its putative regulatory sites were widely present in natural accession, and we also identified the allele which was positively selected during wheat breeding. Epistatic analyses indicated that Qhd.2AS-mediated HD variation is independent of VRN-B1 and environmental factors. Phenotypic investigation of homozygous recombinant inbred lines (RILs) and F2:3 families showed that Qhd.2AS has no negative effect on yield-related traits. These results provide important cues for refining HD and therefore improving yield in wheat breeding programs and will deepen our understanding of the genetic regulation of HD in cereal plants.
Asunto(s)
Sitios de Carácter Cuantitativo , Triticum , Humanos , Mapeo Cromosómico/métodos , Triticum/genética , Fitomejoramiento , Fenotipo , Dedos de Zinc/genéticaRESUMEN
Plant height is one of the most important agronomic traits that affects yield in wheat, owing to that the utilization of dwarf or semi-dwarf genes is closely associated with lodging resistance. In this study, we identified a semi-dwarf mutant, jg0030, induced by γ-ray mutagenesis of the wheat variety 'Jing411' (wild type). Compared with the 'Jing411', plant height of the jg0030 mutant was reduced by 7%-18% in two years' field experiments, and the plants showed no changes in yield-related traits. Treatment with gibberellic acid (GA) suggested that jg0030 is a GA-sensitive mutant. Analysis of the frequency distribution of plant height in 297 F3 families derived from crossing jg0030 with the 'Jing411' indicated that the semi-dwarf phenotype is controlled by a major gene. Using the wheat 660K SNP array-based Bulked Segregant Analysis (BSA) and the exome capture sequencing-BSA assay, the dwarf gene was mapped on the long arm of chromosome 2B. We developed a set of KASP markers and mapped the dwarf gene to a region between marker PH1 and PH7. This region encompassed a genetic distance of 55.21 cM, corresponding to a physical distance of 98.3 Mb. The results of our study provide a new genetic resource and linked markers for wheat improvement in molecular breeding programs.
RESUMEN
Heading date (HD) is an important trait for wide adaptability and yield stability in wheat. The Vernalization 1 (VRN1) gene is a key regulatory factor controlling HD in wheat. The identification of allelic variations in VRN1 is crucial for wheat improvement as climate change becomes more of a threat to agriculture. In this study, we identified an EMS-induced late-heading wheat mutant je0155 and crossed it with wide-type (WT) Jing411 to construct an F2 population of 344 individuals. Through Bulk Segregant Analysis (BSA) of early and late-heading plants, we identified a Quantitative Trait Locus (QTL) for HD on chromosome 5A. Further genetic linkage analysis limited the QTL to a physical region of 0.8 Mb. Cloning and sequencing revealed three copies of VRN-A1 in the WT and mutant lines; one copy contained a missense mutation of C changed to T in exon 4 and another copy contained a mutation in intron 5. Genotype and phenotype analysis of the segregation population validated that the mutations in VRN-A1 contributed to the late HD phenotype in the mutant. Expression analysis of C- or T-type alleles in exon 4 of the WT and mutant lines indicated that this mutation led to lower expression of VRN-A1, which resulted in the late-heading of je0155. This study provides valuable information for the genetic regulation of HD and many important resources for HD refinement in wheat breeding programs.
Asunto(s)
Mutación Missense , Triticum , Triticum/genética , Fitomejoramiento , Mapeo Cromosómico , Sitios de Carácter Cuantitativo , AlelosRESUMEN
Novel genetic variations can be obtained by inducing mutations in the plant which help to achieve novel traits. The useful mutant can be obtained through radiation mutation in a short period which can be used as a new material to produce new varieties with high yield and good quality wheat. In this paper, the proteomic analysis of wheat treated with different doses of 12C and 7Li ion beam radiation at the seedling stage was carried out through a Tandem Mass Tag (TMT) tagging quantitative proteomic analysis platform based on high-resolution liquid chromatography-mass spectrometry, and the traditional 60Co-γ-ray radiation treatment for reference. A total of 4,764 up-regulated and 5,542 down-regulated differentially expressed proteins were identified. These proteins were mainly enriched in the KEGG pathway associated with amino acid metabolism, fatty acid metabolism, carbon metabolism, photosynthesis, signal transduction, protein synthesis, and DNA replication. Functional analysis of the differentially expressed proteins showed that the oxidative defense system in the plant defense system was fully involved in the defense response after 12C ion beam and 7Li ion beam radiation treatments. Photosynthesis and photorespiration were inhibited after 12C ion beam and 60Co-γ-ray irradiation treatments, while there was no effect on the plant with 7Li ion beam treatment. In addition, the synthesis of biomolecules such as proteins, as well as multiple signal transduction pathways also respond to radiations. Some selected differentially expressed proteins were verified by Parallel Reaction Monitoring (PRM) and qPCR, and the experimental results were consistent with the quantitative results of TMT. The present study shows that the physiological effect of 12C ion beam radiation treatment is different as compared to the 7Li ion beam, but its similar to the 60Co-γ ray depicting a significant effect on the plant by using the same dose. The results of this study will provide a theoretical basis for the application of 12C and 7Li ion beam radiation in the mutation breeding of wheat and other major crops and promote the development of heavy ion beam radiation mutation breeding technology.
RESUMEN
Ionizing radiation (IR) is an effective approach for mutation breeding. Understanding the mutagenesis and transcriptional profiles induced by different mutagens is of great significance for improving mutation breeding efficiency. Here, using RNA sequencing and methylation-sensitive amplification polymorphism (MSAP) approaches, we compared the genetic variations, epigenetics, and transcriptional responses induced by the mixed high-energy particle field (CR) and 7Li-ion beam (LR) radiation in M1 seedlings of two wheat genotypes (Yangmai 18 and Yangmai 20). The results showed that, in both wheat genotypes, CR displayed significantly a higher mutation efficiency (1.79 × 10-6/bp) than that by LR (1.56 × 10-6/bp). The induced mutations were not evenly distributed across chromosomes and varied across wheat genotypes. In Y18 M1, the highest number of mutations were detected on Chr. 6B and Chr. 6D, whilst in Y20 M1, Chr. 7A and Chr. 3A had the highest mutations. The transcript results showed that total of 4,755 CR-regulated and 1,054 LR-regulated differentially expressed genes (DEGs) were identified in the both genotypes. Gene function enrichment analysis of DEGs showed that these DEGs overlapped or diverged in the cascades of molecular networks involved in "phenylpropanoid biosynthesis" and "starch and sucrose metabolism" pathways. Moreover, IR type specific responses were observed between CR an LR irradiation, including specific TFs and response pathways. MSAP analysis showed that DNA methylation level increased in LR treatment, while decreased at CR. The proportion of hypermethylation was higher than that of hypomethylation at LR, whereas a reverse pattern was observed at CR, indicating that DNA methylation plays critical roles in response to IR irradiation. All these results support that the response to different IRs in wheat includes both common and unique pathways, which can be served as a useful resource to better understand the mechanisms of responses to different IRs in other plants.
Asunto(s)
Genes de Plantas , Triticum , Alelos , Clonación Molecular , Genes de Plantas/genética , Triticum/genéticaRESUMEN
Plant height is a key agronomic trait that is closely to the plant morphology and lodging resistance in wheat. However, at present, the few dwarf genes widely used in wheat breeding have narrowed wheat genetic diversity. In this study, we selected a semi-dwarf wheat mutant dwarf33 that exhibits decreased plant height with little serious negative impact on other agronomic traits. Genetic analysis and mutant gene mapping indicated that dwarf33 contains a new recessive semi-dwarf gene Rht-SN33d, which was mapped into ~1.3 Mb interval on the 3DL chromosome. The gibberellin metabolism-related gene TraesCS3D02G542800, which encodes gibberellin 2-beta-dioxygenase, is considered a potential candidate gene of Rht-SN33d. Rht-SN33d reduced plant height by approximately 22.4% in mutant dwarf33. Further study revealed that shorter stem cell length may be the main factor causing plant height decrease. In addition, the coleoptile length of dwarf33 was just 9.3% shorter than that of wild-type Shaannong33. These results will help to expand our understanding of new mechanisms of wheat height regulation, and obtain new germplasm for wheat improvement.
Asunto(s)
Giberelinas , Triticum , Triticum/genética , Fitomejoramiento , Mapeo Cromosómico/métodos , FenotipoRESUMEN
Stem elongation is a critical phase for yield determination and, as a major trait, is targeted for manipulation for improvement in bread wheat (Triticum aestivum L.). In a previous study, we characterized a mutant showing rapid stem elongation but with no effect on plant height at maturity. The present study aimed to finely map the underlying mutated gene, qd1, in this mutant. By analyzing an F2 segregating population consisting of 606 individuals, we found that the qd1 gene behaved in a dominant manner. Moreover, by using the bulked segregant RNA sequencing (BSR-seq)-based linkage analysis method, we initially mapped the qd1 gene to a 13.55 Mb region on chromosome 4B (from 15.41 to 28.96 Mb). This result was further confirmed in F2 and BC3F2 segregating populations. Furthermore, by using transcriptome sequencing data, we developed 14 Kompetitive Allele-Specific PCR (KASP) markers and then mapped the qd1 gene to a smaller and more precise 5.08 Mb interval from 26.80 to 31.88 Mb. To develop additional markers to finely map the qd1 gene, a total of 4,481 single-nucleotide polymorphisms (SNPs) within the 5.08 Mb interval were screened, and 25 KASP markers were developed based on 10x-depth genome resequencing data from both wild-type (WT) and mutant plants. The qd1 gene was finally mapped to a 1.33 Mb interval from 28.86 to 30.19 Mb on chromosome 4B. Four candidate genes were identified in this region. Among them, the expression pattern of only TraesCS4B02G042300 in the stems was concurrent with the stem development of the mutant and WT. The qd1 gene could be used in conjunction with molecular markers to manipulate stem development in the future.
RESUMEN
The chemical method for the determination of the resistant starch (RS) content in grains is time-consuming and labor intensive. Near-infrared (NIR) and attenuated total reflectance mid-infrared (ATR-MIR) spectroscopy are rapid and nondestructive analytical techniques for determining grain quality. This study was the first report to establish and compare these two spectroscopic techniques for determining the RS content in wheat grains. Calibration models with four preprocessing techniques based on the partial least squares (PLS) algorithm were built. In the NIR technique, the mean normalization + Savitzky-Golay smoothing (MN + SGS) preprocessing technique had a higher coefficient of determination (R c 2 = 0.672; R p 2 = 0.552) and a relative lower root mean square error value (RMSEC = 0.385; RMSEP = 0.459). In the ATR-MIR technique, the baseline preprocessing method exhibited a better performance regarding to the values of coefficient of determination (R c 2 = 0.927; R p 2 = 0.828) and mean square error value (RMSEC = 0.153; RMSEP = 0.284). The validation of the developed best NIR and ATR-MIR calibration models showed that the ATR-MIR best calibration model has a better RS prediction ability than the NIR best calibration model. Two high grain RS content wheat mutants were screened out by the ATR-MIR best calibration model from the wheat mutant library. There was no significant difference between the predicted values and chemical measured values in the two high RS content mutants. It proved that the ATR-MIR model can be a perfect substitute in RS measuring. All the results indicated that the ATR-MIR spectroscopy with improved screening efficiency can be used as a fast, rapid, and nondestructive method in high grain RS content wheat breeding.
RESUMEN
Manipulation of genes involved in starch synthesis could significantly affect wheat grain weight and yield. The starch-branching enzyme (SBE) catalyzes the formation of branch points by cleaving the α-1,4 linkage in polyglucans and reattaching the chain via an α-1,6 linkage. Three types of SBE isoforms (SBEI, SBEII, and SBEIII) exist in higher plants, with the number of SBE isoforms being species-specific. In this study, the coding sequence of the wheat TaSBEIII gene was amplified. After the multiple sequence alignment of TaSBEIII genome from 20 accessions in a wheat diversity panel, one SNP was observed in TaSBEIII-A, which formed the allelic marker allele-T. Based on this SNP at 294 bp (C/T), a KASP molecular marker was developed to distinguish allelic variation among the wheat genotypes for thousand grain weight (TGW). The results were validated using 262 accessions of mini core collection (MCC) from China, 153 from Pakistan, 53 from CIMMYT, and 17 diploid and 18 tetraploid genotypes. Association analysis between TaSBEIII-A allelic variation and agronomic traits found that TaSBEIII-A was associated with TGW in mini core collection of China (MCC). The accessions possessing Allele-T had higher TGW than those possessing Allele-C; thus, Allele-T was a favorable allelic variation. By analyzing the frequency of the favorable allelic variation Allele-T in MCC, it increased from pre-1950 (25%) to the 1960s (45%) and increased continuously from 1960 to 1990 (80%). The results suggested that the KASP markers can be utilized in grain weight improvement, which ultimately improves wheat yield by marker-assisted selection in wheat breeding. The favorable allelic variation allele-T should be valuable in enhancing grain yield by improving the source and sink simultaneously. Furthermore, the newly developed KASP marker validated in different genetic backgrounds could be integrated into a breeding kit for screening high TGW wheat.
RESUMEN
Blast is caused by the host-specific lineages of the fungus Magnaporthe oryzae and is the most important destructive disease in major crop plants, including rice and wheat. The first wheat blast outbreak that occurred in Bangladesh in 2016 and the recent epidemic in Zambia were caused by the M. oryzae Triticum (MoT) pathotype, a fungal lineage belonging to M. oryzae. Although a few reported wheat cultivars show modest resistance to MoT, the patterns of genetic variation and diversity of this pathotype make it crucial to identify additional lines of resistant wheat germplasm. Nearly 40 rice blast resistant and susceptible genes have so far been cloned. Here, we used BLAST analysis to locate two rice blast susceptible genes in the wheat reference genome, bsr-d1 and bsr-k1, and identified six identical homologous genes located on subgenomes A, B, and D. We uncovered a total of 171 single nucleotide polymorphisms (SNPs) in an ethyl methanesulfonate (EMS)-induced population, with mutation densities ranging from 1/1107.1 to 1/230.7 kb through Targeting Induced Local Lesions IN Genomes (TILLING) by sequencing. These included 81 SNPs located in exonic and promoter regions, and 13 coding alleles that are predicted to have severe effects on protein function, including two pre-mature mutants that might affect wheat blast resistance. The loss-of-function alleles identified in this study provide insights into new wheat blast resistant lines, which represent a valuable breeding resource.
RESUMEN
Agronomic traits such as heading date (HD), plant height (PH), thousand grain weight (TGW), and spike length (SL) are important factors affecting wheat yield. In this study, we constructed a high-density genetic linkage map using the Wheat55K SNP Array to map quantitative trait loci (QTLs) for these traits in 207 recombinant inbred lines (RILs). A total of 37 QTLs were identified, including 9 QTLs for HD, 7 QTLs for PH, 12 QTLs for TGW, and 9 QTLs for SL, which explained 3.0-48.8% of the phenotypic variation. Kompetitive Allele Specific PCR (KASP) markers were developed based on sequencing data and used for validation of the stably detected QTLs on chromosomes 3A, 4B and 6A using 400 RILs. A QTL cluster on chromosome 4B for PH and TGW was delimited to a 0.8 Mb physical interval explaining 12.2-22.8% of the phenotypic variation. Gene annotations and analyses of SNP effects suggested that a gene encoding protein Photosynthesis Affected Mutant 68, which is essential for photosystem II assembly, is a candidate gene affecting PH and TGW. In addition, the QTL for HD on chromosome 3A was narrowed down to a 2.5 Mb interval, and a gene encoding an R3H domain-containing protein was speculated to be the causal gene influencing HD. The linked KASP markers developed in this study will be useful for marker-assisted selection in wheat breeding, and the candidate genes provide new insight into genetic study for those traits in wheat.
RESUMEN
The stems of cereal crops provide both mechanical support for lodging resistance and a nutrient supply for reproductive organs. Elongation, which is considered a critical phase for yield determination in winter wheat (Triticum aestivum L.), begins from the first node detectable to anthesis. Previously, we characterized a heavy ion beam triggered wheat mutant qd, which exhibited an altered stem elongation pattern without affecting mature plant height. In this study, we further analyzed mutant stem developmental characteristics by using transcriptome data. More than 40.87 Mb of clean reads including at least 36.61 Mb of unique mapped reads were obtained for each biological sample in this project. We utilized our transcriptome data to identify 124,971 genes. Among these genes, 4,340 differentially expressed genes (DEG) were identified between the qd and wild-type (WT) plants. Compared to their WT counterparts, qd plants expressed 2,462 DEGs with downregulated expression levels and 1878 DEGs with upregulated expression levels. Using DEXSeq, we identified 2,391 counting bins corresponding to 1,148 genes, and 289 of them were also found in the DEG analysis, demonstrating differences between qd and WT. The 5,199 differentially expressed genes between qd and WT were employed for GO and KEGG analyses. Biological processes, including protein-DNA complex subunit organization, protein-DNA complex assembly, nucleosome organization, nucleosome assembly, and chromatin assembly, were significantly enriched by GO analysis. However, only benzoxazinoid biosynthesis pathway-associated genes were enriched by KEGG analysis. Genes encoding the benzoxazinoid biosynthesis enzymes Bx1, Bx3, Bx4, Bx5, and Bx8_9 were confirmed to be differentially expressed between qd and WT. Our results suggest that benzoxazinoids could play critical roles in regulating the stem elongation phenotype of qd.