Cytoplasmic male sterility (CMS) is important for large-scale hybrid seed production. Rearrangements in the mitochondrial DNA (mtDNA) for the cotton (Gossypium hirsutum L.) CMS line J4A were responsible for pollen abortion. However, the expression patterns of nuclear genes associated with pollen abortion and the molecular basis of CMS for J4A are unknown, and were the objectives of this study by comparing J4A with the J4B maintainer line. Cytological evaluation of J4A anthers showed that microspore abortion occurs during meiosis preventing pollen development. Changes in enzyme activity of mitochondrial respiratory chain complex IV and mitochondrial respiratory chain complex V and the content of ribosomal protein and ATP during anther abortion were observed for J4A suggesting insufficient synthesis of ATP hindered pollen production. Additionally, levels of sucrose, starch, soluble sugar, and fructose were significantly altered in J4A during the meiosis stage, suggesting reduced sugar metabolism contributed to sterility. Transcriptome and miRNAomics analyses identified 4461 differentially expressed mRNAs (DEGs) and 26 differentially expressed microRNAs (DEMIs). Pathway enrichment analysis indicated that the DEMIs were associated with starch and sugar metabolism. Six deduced target gene regulatory pairs that may participate in CMS were identified, ghi-MIR7484-10/mitogen-activated protein kinase kinase 6 (MAPKK6), ghi-undef-156/agamous-like MADS-box protein AGL19 (AGL19), ghi-MIR171-1-22/SNF1-related protein kinase regulatory subunit gamma-1 and protein trichome birefringence-like 38, and ghi-MIR156-(8/36)/WRKY transcription factor 28 (WRKY28). Overall, a putative CMS mechanism involving mitochondrial dysfunction, the ghi-MIR7484-10/MAPKK6 network, and reduced glucose metabolism was suggested, and ghi-MIR7484-10/MAPKK6 may be related to abnormal microspore meiosis and induction of excessive sucrose accumulation in anthers.

DNA methylation is defined as a conserved epigenetic modification mechanism that plays a key role in maintaining normal gene expression without altering the DNA sequence. Several studies have reported that altered methylation patterns were associated with male sterility in some plants such as rice and wheat, but global methylation profiles and their possible roles in cytoplasmic male sterility (CMS), especially in cotton near-isogenic lines, remain unclear. In this study, bisulfite sequencing technology and RNA-Seq were used to investigate CMS line 07-113A and its near-isogenic line 07-113B. Using integrated methylome and transcriptome analyses, we found that the number of hypermethylated genes in the differentially methylated regions, whether in the promoter region or in the gene region, was more in 07-113A than the number in 07-113B. The data indicated that 07-113A was more susceptible to methylation. In order to further analyze the regulatory network of male sterility, transcriptome sequencing and DNA methylation group data were used to compare the characteristics of near-isogenic lines 07-113A and 07-113B in cotton during the abortion stage. Combined methylation and transcriptome analysis showed that differentially expressed methylated genes were mainly concentrated in vital metabolic pathways including the starch and sucrose metabolism pathways and galactose metabolism. And there was a negative correlation between gene methylation and gene expression. In addition, five key genes that may be associated with CMS in cotton were identified. These data will support further understanding of the effect of DNA methylation on gene expression and their potential roles in cotton CMS.

Cytoplasmic male sterility (CMS) is an important feature for achieving heterosis in the development of hybrid crops. Mitochondria contribute to CMS, especially via mitochondrial DNA (mtDNA) rearrangements and chimeric genes. However, the mechanisms of CMS have not been fully elucidated, and the isonuclear alloplasmic lines used in previous studies have limited utility in cotton CMS research. In this study, three CMS lines (J4A-1, J4A-2 and J4A-3) and their isoplasmic maintainer line (J4B) were analyzed for mtDNA structural differences via high-throughput sequencing. The results showed that mtDNA was conserved (with similarities higher than 99%) among the three CMS lines and their isoplasmic maintainer line. All lines harbored 36 known protein-coding genes, 3 rRNAs, and 15 tRNAs. The protein-coding genes with non-synonymous mutations mainly encoded two types of proteins: ATPase and ribosomal proteins. Four new open reading frames (ORFs) (orf116b, orf186a-1, orf186a-2 and orf305a) were identified as candidate ORFs responsible for CMS. Two of the ORFs (orf186a-1 and orf186a-2) were identified as orf4 and orf4-2 of the upland cotton CMS line 2074A (a line with Gossypium harknessii Brandegee CMS-D2-2 cytoplasm), respectively. These findings provide a reference for CMS research in cotton or other crops.