Genetic Analysis and Fine Mapping of QTL for the Erect Leaf in Mutant mths29 Induced through Fast Neutron in Wheat.

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.

Strigolactone and abscisic acid synthesis and signaling pathways are enhanced in the wheat oligo-tillering mutant ot1.

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.

Mapping and identification of a reverse mutation of Rht2 that enhances plant height and thousand grain weight in an elite wheat mutant induced by spaceflight.

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.

Liver Stiffness Measurement is Useful in Predicting Type 2 Diabetes Mellitus Among Nonalcohol Fatty Liver Disease Patients.

Background: Type 2 diabetes mellitus (T2DM) and non-alcoholic fatty liver disease (NAFLD) are closely related conditions. Aim: This study investigated a group of individuals with NAFLD to evaluate if liver fibrosis, identified by FibroScan, correlated with T2DM. Methods: 154 NAFLD patients obtained FibroScan, liver ultrasonography (US), and a thorough assessment of clinical implications and chemical biomarkers. Results: In comparison to the NAFLD without T2DM group, the hemoglobin A1c(HBA1c)(mmol/mol%), homeostasis model of assessment for insulin resistance index (HOMA-IR), gamma-glutamyl transferase (GGT), fibrosis indices, and liver stiffness measurement (LSM) values were all considerably higher in the NAFLD with T2DM group. Patients with NAFLD and T2DM had considerably lower serum uric acid(SUA) levels than those with NAFLD alone.Those with severe fibrosis (79.3%, 23/29) in the NAFLD group showed a greater frequency of T2DM than those with mild fibrosis (45.6%, 21/46) or no fibrosis (27.85%, 22/79) (P=0.000). LSM value and elements of the metabolic syndrome (MetS) were independent risk factors for incident T2DM among NAFLD patients (OR=1.466, 95% CI [1.139-1.888], P=0.003; and OR=0.273, 95% CI [0.081-0.916], P=0.036). Conclusion: FibroScan can identify significant fibrosis, which is independently linked to a higher prevalence of T2DM. As a result, it is crucial to make use of this technology to predict T2DM in NAFLD patients.

Clinical Classification of Obesity and Implications for Metabolic Dysfunction-Associated Fatty Liver Disease and Treatment.

Obesity,and metabolic dysfunction-associated fatty liver disease (MAFLD) have reached epidemic proportions globally. Obesity and MAFLD frequently coexist and act synergistically to increase the risk of adverse clinical outcomes (both hepatic and extrahepatic). Type 2 diabetes mellitus (T2DM) is the most important risk factor for rapid progression of steatohepatitis and advanced fibrosis. Conversely, the later stages of MAFLD are associated with an increased risk of T2DM incident. According to the proposed criteria, MAFLD is diagnosed in patients with liver steatosis and in at least one in three: overweight or obese, T2DM, or signs of metabolic dysregulation if they are of normal weight. However, the clinical classification and correlation between obesity and MAFLD is more complex than expected. In addition, treatment for obesity and MAFLD are associated with a reduced risk of T2DM, suggesting that liver-based treatments could reduce the risk of developing T2DM. This review describes the clinical classification of obesity and MAFLD, discusses the clinical features of various types of obesity and MAFLD, emphasizes the role of visceral obesity and insulin resistance (IR) in the development of MAFLD,and summarizes the existing treatments for obesity and MAFLD that reduce the risk of developing T2DM.

FIB-4 is closer to FibroScan screen results to detecting advanced liver fibrosis and maybe facilitates NAFLD warning.

To assess the relationship between clinical biochemical characteristics and steatosis or fibrosis by Fibroscan in non-alcoholic fatty liver disease (NAFLD) patients in order to seek the simple effective screening method closed to the results of the fibroScan measurement. A cross-sectional study was conducted on 188 patients with NAFLD who underwent FibroScan examinations. Demographic data and clinical biochemical characteristics were collected and analyzed. The result showed elevated serum uric acid (SUA) (P = .023, odds ratio [OR = 1.005, 95% CI (1.001-1.009) and metabolic syndrome (MetS) (P = .000, OR = 4.549, 95%CI (1.974-10.484) were associated with severe steatosis (controlled attenuation parameter, CAP ≥ 300 dB/m). The magnitude of liver stiffness measured using FibroScan was positively correlated with aspartate transaminase/alanine aminotransferase (AST/ALT) ratio (R = 0.419, P = .000), AST to platelet ratio index (APRI) score (R = 0.309, P = .000), and Fibrosis-4 score (FIB-4) (R = 0.507, P = .000). The areas under the receiver operating curve (ROC) of AST/ALT, APRI, and FIB-4 for mild or severe fibrosis were 0.563, 0.696, and 0.728, respectively, and those for advanced fibrosis were 0.648, 0.750, and 0.821, respectively. The FIB-4 index cutoff value was 1.65 with a sensitivity of 68.3% and specificity of 89.8% during the diagnosis of advanced fibrosis. MetS and elevated SUA are associated with severe steatosis according to the CAP value screen, whereas FIB-4, as the fibrosis score method, is closer to the liver stiffness measurement results from FibroScan, which may facilitate early warning of NAFLD in the community or in remote areas.

Genome-wide characterization of two homeobox families identifies key genes associated with grain-related traits in wheat.

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.

A large-scale whole-exome sequencing mutant resource for functional genomics in wheat.

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.

Genetic mapping and identification of Rht8-B1 that regulates plant height in wheat.

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.

Fine mapping and genetic analysis identified a C2H2-type zinc finger as a candidate gene for heading date regulation in wheat.

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.

Genetic analysis and mapping of dwarf gene without yield penalty in a γ-ray-induced wheat mutant.

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.

G-quadruplex-enhanced circular single-stranded DNA (G4-CSSD) adsorption of miRNA to inhibit colon cancer progression.

BACKGROUND: Chromosomal heterogeneity leads to the abnormal expression and mutation of tumor-specific genes. Drugs targeting oncogenes have been extensively developed. However, given the random mutation of tumor suppressor genes, the development of its targeted drugs is difficult. METHODS: Our early research revealed that artificial circular single-stranded DNA (CSSD) can restore multiple tumor suppressor genes to inhibit tumor malignant progression by adsorbing miRNA. Here, we improved CSSD to a fully closed single-stranded DNA with G quadruplex DNA secondary structure (G4-CSSD), which made G4-CSSD with higher acquisition rate and decreased degradation. The Cancer Genome Atlas (TCGA) and Human Protein Atlas database were used to predict tumour suppressor genes in colon cancer. Cellular and animal experiments were performed to validate the role of G4-CSSD in cancer cell progression. RESULTS: In colon cancer, we observed the simultaneous low expressions of chloride channel accessory 1 (CLCA1), UDP-GlcNAc:betaGal beta-1,3-N-acetylglucosaminyltransferase 6 (B3GNT6) and UDP glucuronosyltransferase family 2 member A3 (UGT2A3), which indicated an favourable prognosis. After repressing miR-590-3p with G4-CSSD590, the upregulation of CLCA1, B3GNT6 and UGT2A3 inhibited the proliferation and metastasis of colon cancer cells. CONCLUSIONS: This study may provide basis for new treatment methods for colon cancer by restoration of tumor suppressor genes.

Gene Mapping and Identification of a Missense Mutation in One Copy of VRN-A1 Affects Heading Date Variation in Wheat.

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.

Physiological and Differential Proteomic Analysis at Seedling Stage by Induction of Heavy-Ion Beam Radiation in Wheat Seeds.

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.

Fine Mapping of qd1, a Dominant Gene that Regulates Stem Elongation in Bread Wheat.

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.

Genetic Mapping by Integration of 55K SNP Array and KASP Markers Reveals Candidate Genes for Important Agronomic Traits in Hexaploid Wheat.

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.

Transcriptome Analysis Reveals a Potential Role of Benzoxazinoid in Regulating Stem Elongation in the Wheat Mutant qd.

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.

Genome-Wide and Exome-Capturing Sequencing of a Gamma-Ray-Induced Mutant Reveals Biased Variations in Common Wheat.

Induced mutagenesis is a powerful approach for the creation of novel germplasm and the improvement of agronomic traits. The evaluation of mutagenic effects and functional variations in crops is needed for breeding mutant strains. To investigate the mutagenic effects of gamma-ray irradiation in wheat, this study characterized genomic variations of wheat early heading mutant (eh1) as compared to wild-type (WT) Zhongyuan 9 (ZY9). Whole-genome resequencing of eh1 and ZY9 produced 737.7 Gb sequencing data and identified a total of 23,537,117 homozygous single nucleotide polymorphism (SNP) and 1,608,468 Indel. Analysis of SNP distribution across the chromosome suggests that mutation hotspots existed in certain chromosomal regions. Among the three subgenomes, the variation frequency in subgenome D was significantly lower than in subgenomes A and B. A total of 27.8 Gb data were obtained by exome-capturing sequencing, while 217,948 SNP and 13,554 Indel were identified. Variation annotation in the gene-coding sequences demonstrated that 5.0% of the SNP and 5.3% of the Indel were functionally important. Characterization of exomic variations in 12 additional gamma-ray-induced mutant lines further provided additional insights into the mutagenic effects of this approach. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genome (KEGG) analysis suggested that genes with functional variations were enriched in several metabolic pathways, including plant-pathogen interactions and ADP binding. Kompetitive allele-specific PCR (KASP) genotyping with selected SNP within functional genes indicated that 85.7% of the SNPs were polymorphic between the eh1 and wild type. This study provides a basic understanding of the mechanism behind gamma-ray irradiation in hexaploid wheat.

Identification of the Q Gene Playing a Role in Spike Morphology Variation in Wheat Mutants and Its Regulatory Network.

The wheat AP2 family gene Q controls domestication traits, including spike morphology and threshability, which are critical for the widespread cultivation and yield improvement of wheat. Although many studies have investigated the molecular mechanisms of the Q gene, its direct target genes, especially those controlling spike morphology, are not clear, and its regulatory pathways are not well established. In this study, we conducted gene mapping of a wheat speltoid spike mutant and found that a new allele of the Q gene with protein truncation played a role in spike morphology variation in the mutant. Dynamic expression levels of the Q gene throughout the spike development process suggested that the transcript abundances of the mutant were decreased at the W6 and W7 scales compared to those of the WT. We identified several mutation sites on the Q gene and showed that mutations in different domains resulted in distinct phenotypes. In addition, we found that the Q gene produced three transcripts via alternative splicing and that they exhibited differential expression patterns in nodes, internodes, flag leaves, and spikes. Finally, we identified several target genes directly downstream of Q, including TaGRF1-2D and TaMGD-6B, and proposed a possible regulatory network. This study uncovered the target genes of Q, and the results can help to clarify the mechanism of wheat spike morphology and thereby improve wheat grain yield.