Integrated Analysis of Microarray, Small RNA, and Degradome Datasets Uncovers the Role of MicroRNAs in Temperature-Sensitive Genic Male Sterility in Wheat.

Temperature-sensitive genic male sterile (TGMS) line Beijing Sterility 366 (BS366) has been utilized in hybrid breeding for a long time, but the molecular mechanism underlying male sterility remains unclear. Expression arrays, small RNA, and degradome sequencing were used in this study to explore the potential role of miRNA in the cold-induced male sterility of BS366. Microspore observation showed defective cell plates in dyads and tetrads and shrunken microspores at the vacuolated stage. Differential regulation of Golgi vesicle transport, phragmoplast formation, sporopollenin biosynthesis, pollen exine formation, and lipid metabolism were observed between cold and control conditions. Pollen development was significantly represented in the 352 antagonistic miRNA-target pairs in the integrated analysis of miRNA and mRNA profiles. The specific cleavage of ARF17 and TIR1 by miR160 and miR393 were found in the cold-treated BS366 degradome, respectively. Thus, the cold-mediated miRNAs impaired cell plate formation through repression of Golgi vesicle transport and phragmoplast formation. The repressed expression of ARF17 and TIR1 impaired pollen exine formation. The results of this study will contribute to our understanding of the roles of miRNAs in male sterility in wheat.

Comparative transcriptome and DNA methylation analysis in temperature-sensitive genic male sterile wheat BS366.

BACKGROUND: Known as the prerequisite component for the heterosis breeding system, the male sterile line determines the hybrid yield and seed purity. Therefore, a deep understanding of the mechanism and gene network that leads to male sterility is crucial. BS366, a temperature-sensitive genic male sterile (TGMS) line, is male sterile under cold conditions (12 °C with 12 h of daylight) but fertile under normal temperature (20 °C with 12 h of daylight). RESULTS: During meiosis, BS366 was defective in forming tetrads and dyads due to the abnormal cell plate. During pollen development, unusual vacuolated pollen that could not accumulate starch grains at the binucleate stage was also observed. Transcriptome analysis revealed that genes involved in the meiotic process, such as sister chromatid segregation and microtubule-based movement, were repressed, while genes involved in DNA and histone methylation were induced in BS366 under cold conditions. MethylRAD was used for reduced DNA methylation sequencing of BS366 spikes under both cold and control conditions. The differentially methylated sites (DMSs) located in the gene region were mainly involved in carbohydrate and fatty acid metabolism, lipid metabolism, and transport. Differentially expressed and methylated genes were mainly involved in cell division. CONCLUSIONS: These results indicated that the methylation of genes involved in carbon metabolism or fatty acid metabolism might contribute to male sterility in BS366 spikes, providing novel insight into the molecular mechanism of wheat male sterility.

Genome-wide identification and transcriptional characterization of DNA methyltransferases conferring temperature-sensitive male sterility in wheat.

BACKGROUND: DNA methyltransferase (DMT) genes contribute to plant stress responses and development by de novo establishment and subsequent maintenance of DNA methylation during replication. The photoperiod and/or temperature-sensitive genic male sterile (P/TGMS) lines play an important role in hybrid seed production of wheat. However, only a few studies have reported on the effect of DMT genes on temperature-sensitive male sterility of wheat. Although DMT genes have been investigated in some plant species, the identification and analysis of DMT genes in wheat (Triticum aestivum L.) based on genome-wide levels have not been reported. RESULTS: In this study, a detailed overview of phylogeny of 52 wheat DMT (TaDMT) genes was presented. Homoeolog retention for TaDMT genes was significantly above the average retention rate for whole-wheat genes, indicating the functional importance of many DMT homoeologs. We found that the strikingly high number of TaDMT genes resulted mainly from the significant expansion of the TaDRM subfamily. Intriguingly, all 5 paralogs belonged to the wheat DRM subfamily, and we speculated that tandem duplications might play a crucial role in the TaDRM subfamily expansion. Through the transcriptional analysis of TaDMT genes in a TGMS line BS366 and its hybrids with the other six fertile lines under sterile and fertile conditions, we concluded that TaCMT-D2, TaMET1-B1, and TaDRM-U6 might be involved in male sterility in BS366. Furthermore, a correlation analysis showed that TaMET1-B1 might negatively regulate the expression of TaRAFTIN1A, an important gene for pollen development, so we speculated regarding an epigenetic regulatory mechanism underlying the male sterility of BS366 via the interaction between TaMET1-B1 and TaRAFTIN1A. CONCLUSIONS: Our findings presented a detailed phylogenic overview of the DMT genes and could provide novel insights into the effects of DMT genes on TGMS wheat.

A high-throughput KASP assay provides insights into the evolution of multiple resistant mutations in populations of the two-spotted spider mite Tetranychus urticae across China.

BACKGROUND: The two-spotted spider mite (TSSM), Tetranychus urticae (Acari: Tetranychidae), is a cosmopolitan phytophagous pest in agriculture and horticulture. It has developed resistance to many acaricides by target-site mutations. Understanding the status and evolution of resistant mutations in the field is essential for resistance management. Here, we applied a high-throughput Kompetitive allele-specific polymerase chain reaction (KASP) method for detecting six mutations conferring resistance to four acaricides of the TSSM. We genotyped 3274 female adults of TSSM from 43 populations collected across China in 2017, 2020, and 2021. RESULTS: The KASP genotyping of 24 testing individuals showed 99% agreement with Sanger sequencing results. KASP assays showed that most populations had a high frequency of mutations conferring avermectin (G314D and G326E) and pyridaben (H92R) resistance. The frequency of mutation conferring bifenazate (A269V and G126S) and etoxazole (I1017F) resistance was relatively low. Multiple mutations were common in the TSSM, with 70.2% and 24.6% of individuals having 2-6 and 7-10 of 10 possible resistant alleles, respectively. No loci were linked in most populations among the six mutations, indicating the development of multiple resistance is mainly by independent selection. However, G314D and I1017F on the nuclear genome deviated from Hardy-Weinberg equilibrium in most populations, indicating significant selective pressure on TSSM populations by acaricides or fitness cost of the mutations in the absence of acaricide selection. CONCLUSION: Our study revealed that the high frequency of TSSMs evolved multiple resistant mutations in population and individual levels by independent selection across China, alarming for managing multiple-acaricides resistance. © 2023 Society of Chemical Industry.

KASP genotyping and semi-quantitation of G275E mutation in the α6 subunit of Thrips palmi nAChR gene conferring spinetoram resistance.

BACKGROUND: Pesticide resistance is a long-standing and growing problem in the chemical control of invertebrate pests. Molecular diagnostic methods can facilitate pesticide resistance management by accurately and efficiently detecting resistant mutations and their frequency. In this study, the kompetitive allele specific PCR (KASP) approach, a technology for high-throughput single nucleotide polymorphism (SNP) genotyping, is validated as a useful method for characterizing genotypes at a pesticide-resistance locus for the first time. We focus on the spinetoram resistance mutation of G275E in the nicotinic acetylcholine receptor alpha 6 (nAChR α6) subunit gene of Thrips palmi. RESULTS: Of the 341 individuals of Thrips palmi tested, 98.24% were successfully genotyped, with 100% concordance with Sanger sequencing results. We then quantitatively mixed genomic DNA of known genotypes to establish 21 DNA mixtures with a resistant allele frequency ranging from 0 to 100% at steps of 5%. The linear discriminant analysis (LDA) showed that 75.8% of original grouped cases were correctly classified; six groups had no overlap in membership (resistant allele frequency: 0%, 5%, 10-75%, 80-85%, 90-95%, and 100%). When we chose 11 pooled samples with 10% steps for LDA, 84.4% of original grouped cases were correctly classified; seven groups had no overlap in membership (0%, 10%, 20-30%, 40-70%, 80%, 90%, 100%). The results indicated that KASP applied to pooled samples may provide a semi-quantitative estimate of resistance. CONCLUSIONS: Our study points to the suitability of KASP for high-throughput genotyping of genotypes affecting pesticide resistance and semi-quantitative assessments of resistance allele frequencies in populations. © 2023 Society of Chemical Industry.

Characterization and evolutionary analysis of phosphate starvation response genes in wheat and other major gramineous plants.

Developing cultivars with improved Pi use efficiency is essential for the sustainability of agriculture as well as the environment. Phosphate starvation response (PHR) regulators have not yet been systematically studied in wheat. This study provides the detailed characteristics of PHRs in hexaploid wheat as well as other major gramineous plants at the genome-wide level. The identified PHR proteins were divided into six subfamilies through phylogeny analysis, and a total of 63 paralogous TaPHR pairs were designated as arising from duplication events, with strong purifying selection. The promoters of TaPHRs were identified as stations for many transcription factors. Protein-protein interaction network and gene ontology enrichment analysis indicated a core biological process of cellular response to phosphate starvation. The three-dimensional structures of core PHR proteins showed a high phylogenetic relationship, but amino acid deletions in core protein domains may cause functional differentiation between rice and wheat. TaPHR3 could interact with TaSPX1 and TaSPX5 proteins, which is regarded as a novel interaction mode. Under different Pi gradient treatments, TaPHRs showed low inducible expression patterns among all subfamilies. Our study is the first to comprehensively clarify the basic properties of TaPHR proteins and might accumulate basic data for improving grain yield and environmental homeostasis.

Determination of wheat moisture using terahertz spectroscopy combined with the tabu search algorithm.

The detection of the wheat moisture content plays a key role in grain storage and classification. Harvested wheat grains were taken as samples in the current research. A total of 240 reaped wheat samples with different moisture contents were tested by applying terahertz (THz) spectroscopy. The frequency domain spectra and absorption coefficient spectra of wheat were obtained in the band of 0.1-1.2 THz, and the spectra were pretreated by mean centering, Savitzky-Golay (S-G), Multiplicative Scatter Correction (MSC) and Stand Normal Variate (SNV), respectively. Then a special algorithm of Tabu Search (TS) was used to find out the effective variables and remove the useless variables from the terahertz spectrum of the sample. Finally, the partial least squares (PLS) of chemometrics were used for quantitative model building and prediction. The correlation coefficient of calibration (Rc) is 0.9522. The root mean square error of calibration (RMSEC) is 0.4730. The correlation coefficient of prediction (Rp) is 0.9531. The root mean square error of prediction (RMSEP) is 0.5396. The results demonstrated that an accurate quantitative analysis of moisture in wheat samples could be achieved by terahertz time-domain spectroscopy combined with the TS algorithm. In addition, the results show that the model S-G + MSC + TS + PLS can effectively predict wheat moisture, and provide a rapid quantitative detection and analysis method for the detection of wheat moisture.

Isolation and molecular characterization of high molecular weight glutenin subunit genes 1Bx13 and 1By16 from hexaploid wheat.

The high molecular weight glutenin subunit (HMW-GS) pair 1Bx13 + 1By16are recognized to positively correlate with bread-making quality; however, their molecular data remain unknown. In order to reveal the mechanism by which 1By16 and 1Bx13 creates high quality, their open reading frames (ORFs) were amplified from common wheat Atlas66 and Jimai 20 using primers that were designed based on published sequences of HMW glutenin genes. The ORF of 1By16 was 2,220 bp, deduced into 738 amino acid residues with seven cysteines including 59 hexapeptides and 22 nanopeptides motifs. The ORF of 1Bx13 was 2,385 bp, deduced into 795 amino acid residues with four cysteines including 68 hexapeptides, 25 nanopeptides and six tripeptides motifs. We found that 1By16 was the largest y-type HMW glutenin gene described to date in common wheat. The 1By16 had 36 amino acid residues inserted in the central repetitive domain compared with 1By15. Expression in bacteria and western-blot tests confirmed that the sequence cloned was the ORF of HMW-GS 1By16, and that 1Bx13 was one of the largest 1Bx genes that have been described so far in common wheat, exhibiting a hexapeptide (PGQGQQ) insertion in the end of central repetitive domain compared with 1Bx7. A phylogenetic tree based on the deduced full-length amino acid sequence alignment of the published HMW-GS genes showed that the 1By16 was clustered with Glu-1B-2, and that the 1Bx13 was clustered with Glu-1B-1 alleles.

Production of a monoclonal antibody specific for high molecular weight glutenin subunits (HMW-GS) in wheat and its antigenic determinant.

Wheat high molecular weight glutenin subunits (HMW-GS) 1Bx14 and 1By15 isolated by preparative SDS-PAGE are used as antigen to immunize BALB/c mice. Subcutaneous inoculation of the antigen is performed. The intra-peritoneal injection is completed 3 days before fusion with myeloma cell (SP2/0) via PEG-1500. The fusion cells are selected by indirect enzyme-linked immuno-sorbent assay (ELISA). Positive hybrid cells are further verified three times by limit dilution of the culture cells. A hybridoma cell line is successfully obtained. The monoclonal antibody belongs to IgG1 subclass. In immunoblotting, the antibody binds to all HMW-GS of T. aestivum cultivars, but does not bind to other storage proteins in seeds of wheat. This result is consisting with the high homology in amino acid sequences among the HMW glutenin subunits in wheat. The antibody also binds to HMW-GS storage proteins in Aegilops squarrosa and T. durum (durum wheat). Furthermore, it also binds to HMW storage proteins in Secale cereale (rye), Hordeum vulgare (barley). However, it never binds seed storage proteins in other cereals such as maize, oat, rice, foxtail millet, sorghum etc. The antigen determinant recognized by the antibody has been located within hexapeptide [PGQGQQ] or / and nonapeptide [GYYPTSPQQ] in the central repetitive region of HMW-GS.

Expression of wheat high molecular weight glutenin subunit 1Bx is affected by large insertions and deletions located in the upstream flanking sequences.

To better understand the transcriptional regulation of high molecular weight glutenin subunit (HMW-GS) expression, we isolated four Glu-1Bx promoters from six wheat cultivars exhibiting diverse protein expression levels. The activities of the diverse Glu-1Bx promoters were tested and compared with ß-glucuronidase (GUS) reporter fusions. Although all the full-length Glu-1Bx promoters showed endosperm-specific activities, the strongest GUS activity was observed with the 1Bx7OE promoter in both transient expression assays and stable transgenic rice lines. A 43 bp insertion in the 1Bx7OE promoter, which is absent in the 1Bx7 promoter, led to enhanced expression. Analysis of promoter deletion constructs confirmed that a 185 bp MITE (miniature inverted-repeat transposable element) in the 1Bx14 promoter had a weak positive effect on Glu-1Bx expression, and a 54 bp deletion in the 1Bx13 promoter reduced endosperm-specific activity. To investigate the effect of the 43 bp insertion in the 1Bx7OE promoter, a functional marker was developed to screen 505 Chinese varieties and 160 European varieties, and only 1Bx7-type varieties harboring the 43 bp insertion in their promoters showed similar overexpression patterns. Hence, the 1Bx7OE promoter should be important tool in crop genetic engineering as well as in molecular assisted breeding.