Search

Resolving the genetic basis of fertility restoration for cytoplasmic male sterility (CMS) can improve the efficiency of three-line hybrid breeding. However, the genetic determinants of male fertility restoration in cotton are still largely unknown. This study comprehensively compared the full-length transcripts of CMS-D2 and CMS-D8 systems to identify potential genes linked with fertility restorer genes Rf1 or Rf2. Target comparative analysis revealed a higher percentage of differential genes in each restorer line as compared to their corresponding sterile and maintainer lines. An array of genes with specific expression in the restorer line of CMS-D2 had functional annotations related to floral development and pathway enrichments in various secondary metabolites, while specifically expressed genes in the CMS-D8 restorer line showed functional annotations related to anther development and pathway enrichment in the biosynthesis of secondary metabolites. Further analysis identified potentially key genes located in the target region of fertility restorer genes Rf1 or Rf2. In particular, Ghir_D05G032450 can be the candidate gene related to restorer gene Rf1, and Ghir_D05G035690 can be the candidate gene associated with restorer gene Rf2. Further gene expression validation with qRT-PCR confirmed the accuracy of our results. Our findings provide useful insights into decoding the potential regulatory network that retrieves pollen fertility in cotton and will help to further reveal the differences in the genetic basis of fertility restoration for two CMS systems.

CMS-D2 is an economical and effective system for producing hybrid cotton seeds than artificial and chemical emasculation methods. However, the unstable restoring ability of restorer lines is a main barrier in the large-scale application of "three-line" hybrid cotton. Our phenotypic investigation determined that the homozygous Rf1Rf1 allelic genotype had a stronger ability to generate fertile pollen than the heterozygous Rf1rf1 allelic genotype. To decipher the genetic mechanisms that control the differential levels of pollen fertility, an integrated metabolomic and transcriptomic analysis was performed on pollen grains at two environments using four cotton genotypes differing in Rf1 alleles or cytoplasm. Totally 5,391 differential metabolite features were detected, and 369 specific differential metabolites (DMs) were identified between homozygous and heterozygous Rf1 allelic genotypes with CMS-D2 cytoplasm. Additionally, transcriptome analysis identified 2,490 differentially expressed genes (DEGs) and 96 unique hub DEGs with dynamic regulation in this comparative combination. Further integrated analyses revealed that several key DEGs and DMs involved in indole biosynthesis, flavonoid biosynthesis, and sugar metabolism had strong network linkage with fertility restoration. In vitro application of auxin analogue NAA and inhibitor Auxinole confirmed that over-activated auxin signaling might inhibit pollen development whereas suppressing auxin signaling partially promoted pollen development in CMS-D2 cotton. Our results provide new insight into how the dosage effects of the Rf1 gene regulate pollen fertility of CMS-D2 cotton.

Hybridization is useful to enhance the yield potential of agronomic crops in the world. Cotton has genome doubling due to the allotetraploid process and hybridization in coordination with duplicated genome can produce more yield and adaptability. Therefore, the expression of homoeologous gene pairs between hybrids and inbred parents is vital to characterize the genetic source of heterosis in cotton. Investigation results of homoeolog gene pairs between two contrasting hybrids and their respective inbred parents identified 36853 homoeolog genes in hybrids. It was observed both high and low hybrids had similar trends in homoeolog gene expression patterns in each tissue under study. An average of 96% of homoeolog genes had no biased expression and their expressions were derived from the equal contribution of both parents. Besides, very few homoeolog genes (an average of 1%) showed no biased or novel expression in both hybrids. The functional analysis described secondary metabolic pathways had a majority of novel biased homoeolog genes in hybrids. These results contribute preliminary knowledge about how hybridization affects expression patterns of homoeolog gene pairs in upland cotton hybrids. Our study also highlights the functional genomics of metabolic genes to explore the genetic mechanism of heterosis in cotton.

Heterosis refers to the superior phenotypes observed in hybrids. Cytoplasmic male sterility (CMS) system plays an important role in cotton heterosis utilization. However, the global gene expression patterns of CMS-D2 and its interaction with the restorer gene Rf1 remain unclear. Here, the full-length transcript sequencing was performed in anthers of the CMS-D2 restorer line using PacBio single-molecule real-time sequencing technology. Combining PacBio SMRT long-read isoforms and Illumina RNA-seq data, 107,066 isoforms from 44,338 loci were obtained, including 10,086 novel isoforms of novel genes and 66,419 new isoforms of known genes. Totally 56,572 alternative splicing (AS) events, 1146 lncRNAs, 61 fusion transcripts and 10,466 genes exhibited alternative polyadenylation (APA), and 60,995 novel isoforms with predicted open reading frames (ORFs) were further identified. Furthermore, the specifically expressed genes in restorer line were selected and confirmed by qRT-PCR. These findings provide a basis for upland cotton genome annotation and transcriptome research, and will help to reveal the molecular mechanism of interaction between Rf1 and CMS-D2 cytoplasm.

Deleterious effects on anther development and main economy traits caused by sterile genes or cytoplasms are one of the important genetic characteristics of cytoplasmic male sterility (CMS) systems in cotton, which severely hinder the large-scale application of "three-line" hybrids in production. Therefore, distinct characterization of each cytoplasmic type is mandatory to improve the breeding efficiency of cotton hybrids. In this study, four isonuclear-alloplasmic cotton male sterile lines with G. hirsutum (CMS-(AD)1), G. barbadense (CMS-(AD)2), G. harknessii (CMS-D2), and G. trilobum (CMS-D8) cytoplasms were first created by multiple backcrosses with common genotype Shikang126. Then, 64 pairs of mitochondrial simple sequence repeat (mtSSR) markers were designed to explore the mitochondrial DNA diversities among four isonuclear-alloplasmic cotton male sterile lines, and a total of nine pairs of polymorphic mtSSR molecular markers were successfully developed. Polymorphism analysis indicated that mtSSR59 marker correlated to the atp1 gene could effectively divide the CMS-D2, CMS-(AD)1, and CMS-(AD)2 in one category while the CMS-D8 in another category. Further cytological observation and determination of ATP contents also confirmed the accurate classification of CMS-D2 and CMS-D8 lines. Moreover, the mtSSR59 marker was successfully applied in the marker-assisted selection (MAS) for breeding new male sterile lines and precise differentiation or purity identification of different CMS-based "three-line" and conventional cotton hybrids. This study provides new technical measures for classifying various cytoplasmic sterile lines, and our results will significantly improve the efficiency of there-line hybrid breeding in cotton.

Cellulose synthases (CesAs) are multi-subunit enzymes found on the plasma membrane of plant cells and play a pivotal role in cellulose production. The cotton fiber is mainly composed of cellulose, and the genetic relationships between CesA genes and cotton fiber yield and quality are not fully understood. Through a phylogenetic analysis, the CesA gene family in diploid Gossypium arboreum and Gossypium raimondii, as well as tetraploid Gossypium hirsutum ('TM-1') and Gossypium barbadense ('Hai-7124' and '3-79'), was divided into 6 groups and 15 sub-groups, with each group containing two to five homologous genes. Most CesA genes in the four species are highly collinear. Among the five cotton genomes, 440 and 1929 single nucleotide polymorphisms (SNPs) in the CesA gene family were identified in exons and introns, respectively, including 174 SNPs resulting in amino acid changes. In total, 484 homeologous SNPs between the A and D genomes were identified in diploids, while 142 SNPs were detected between the two tetraploids, with 32 and 82 SNPs existing within G. hirsutum and G. barbadense, respectively. Additionally, 74 quantitative trait loci near 18 GhCesA genes were associated with fiber quality. One to four GhCesA genes were differentially expressed (DE) in ovules at 0 and 3 days post anthesis (DPA) between two backcross inbred lines having different fiber lengths, but no DE genes were identified between these lines in developing fibers at 10 DPA. Twenty-seven SNPs in above DE CesA genes were detected among seven cotton lines, including one SNP in Ghi_A08G03061 that was detected in four G. hirsutum genotypes. This study provides the first comprehensive characterization of the cotton CesA gene family, which may play important roles in determining cotton fiber quality.

Background: Heterosis has greatly improved productivity of cotton crop. However, major constraint for large-scale promotion of F1 hybrid cotton is artificial emasculation and pollination. This study through comparative evaluation of hybrid generations proposed potential utilization of F2 hybrids to improve upland cotton production.Aims：The current study aimed to compare F1 and F2 hybrids and combine with the breeding practice of strong hybrid cotton to select the best combination of heterosis, provide reference for the feasibility of parental selection and utilization of F2 heterosis in China.Methods:We analyzed eight upland cotton varieties and made crosses according to NCⅡ incomplete diallel cross breeding design in two cotton belts of China.Results: The results showed that Variance analysis revealed significant differences for agronomic, yield and fiber quality in both generations and environments. The broad-sense heritability of agronomic and yield traits were relatively higher than quality traits. Furthermore, narrow sense heritability of some traits was higher in F2 than F1 generation under both cotton belts. Overall, parental lines Zhong901, ZB, L28, and Z98 were observed with maximum combining ability while combinations with strong special combining ability were ZB ×DT, SJ48 ×Z98, L28 ×Z98, and ZB ×851. The yield traits heterosis was predominant in both generations. However, the level of heterosis was altered with trait, hybrid combination, generation, and environment. Interestingly, L28×Z98 performed outstandingly in Anyang. Its lint yield was 24.8% higher in F1 and 11.6% in F2 than that of the control Ruiza816. The performance of SJ48 × Z98 was excellent in Alar which showed 36.5% higher LY in F1 and 10.9% in F2 than control Zhongmiansuo 49. Further results revealed most hybrid combinations had shown low level of heterosis for agronomic and fiber quality traits in both generations. Comparatively, ZB ×DT and L28×Z98 showed hybrid vigor for multiple traits in both generations and cotton belts. Moreover, lint yield can be improved independently in both cotton belts in these combinations.Conclusions: Through comparative analysis of variance, combining ability, and heterosis in F1 and F2 hybrids under different cotton belts, this study proposed potential utilization of F2 hybrids to improve upland cotton productivity in China.

Cotton is the most important natural fiber used in textiles. Breeding for "three-lines", i.e., cytoplasmic male sterility (CMS)-based sterile (A), maintainer (B), and restorer (R) line, is a promising approach to harness hybrid vigor in cotton. Pentatricopeptide repeat (PPR) protein-encoding genes play an important role in plant growth and development including restoration of CMS plants to male fertility. However, PPRs, especially those contributing to CMS and fiber development, remain largely unknown in cotton. In this study, a genome-wide identification and characterization of PPR gene family in four Gossypium species with genome sequences (G. arboreum, G. raimondii, G. hirsutum, and G. barbadense) were performed, and expressed PPR genes in developing floral buds, ovules, and fibers were compared to identify possible PPRs related to CMS restoration and fiber development. A total of 539, 558, 1032, and 1055 PPRs were predicted in the above four species, respectively, which were further mapped to chromosomes for a synteny analysis. Through an RNA-seq analysis, 86% (882) PPRs were expressed in flowering buds of upland cotton (G. hirsutum); however, only 11 and 6 were differentially expressed (DE) between restorer R and its near-isogenic (NI) B and between R and its NI A line, respectively. Another RNA-seq analysis identified the expression of only 54% (556) PPRs in 0 and 3 day(s) post-anthesis (DPA) ovules and 24% (247) PPRs in 10 DPA fibers; however, only 59, 6, and 27 PPRs were DE in 0 and 3 DPA ovules, and 10 DPA fibers between two backcross inbred lines (BILs) with differing fiber length, respectively. Only 2 PPRs were DE between Xuzhou 142 and its fiberless and fuzzless mutant. Quantitative RT-PCR analysis confirmed the validity of the RNA-seq results for the gene expression pattern. Therefore, only a very small number of PPRs may be associated with fertility restoration of CMS and genetic differences in fiber initiation and elongation. These results lay a foundation for understanding the roles of PPR genes in cotton, and will be useful in the prioritization of candidate PPR gene functional validation for cotton CMS restoration and fiber development.

The objective of this study was to examine the relationship between heterosis and genetic distance in upland cotton (Gossypium hirsutum L.). Heterosis in the different environments was evaluated and the relationship between heterosis and genetic distance (GD) was determined based on SSR markers for yield quantitative trait loci (QTL). Yields of seed cotton and lint showed a linear relationship with mid‐parent heterosis (MPH) and better‐parent heterosis (BPH). The variation in heterosis for other traits and their correlation with GD may be due to environmental factors or the effort of QTLs tested in the present study may vary in different environments. The present findings provide a foundation for heterotic grouping of parental lines and breeding of new cotton hybrids with improved seed cotton yield. This study calls for more research with stable QTLs as well as advance molecular marker techniques may be used in predicting yield heterosis in a more precise and reliable manner.

Background: Hybridization is useful to enhance yield potential of agronomic crops in the world. Cotton has genome doubling due to alloteraploidy process and hybridization process in coordinate with duplicated genome can produce more yield and adaptability. Therefore, expression of homoeolog gene pairs between hybrids and inbred parents are vital to characterize genetic source of heterosis in cotton.Results: Investigation results of homoeolog gene pairs between two contrasting hybrids and their respective inbred parents identified 36853 homoeolog genes in hybrids. It was observed both high and low hybrids had similar trends in homoeolog gene expression patterns in each tissue under study. An average of 96% of homoeolog genes had no biased expression and their expressions were derived from the equal contribution of both parents. Besides, very few homoeolog genes (An average of 1%) showed no biased or novel expression in both hybrids. The functional analysis described secondary metabolic pathways had a majority of novel biased homoeolog genes in hybrids. Conclusions: These results contribute preliminary knowledge about how hybridization affects expression patterns of homoeolog gene pairs in upland cotton hybrids. Our study also highlights the functional genomics of metabolic genes to explore the genetic mechanism of heterosis in cotton.

Background: Cytoplasmic male sterile (CMS) with cytoplasm from Gossypium Trilobum (D8) fails to produce functional pollen. It is useful for commercial hybrid cotton seed production. The restore line of CMS-D8 containing Rf2 gene can restore the fertility of the corresponding sterile line. This study combined the whole genome resequencing bulked segregant analysis (BSA) with high-throughput SNP genotyping to accelerate the physical mapping of Rf2 locus in CMS-D8 cotton. Methods: The fertility of backcross population ((sterile line×restorer line)×maintainer line) comprising of 1623 individuals was investigated in the field. The fertile pool (100 plants with fertile phenotypes, F-pool) and the sterile pool (100 plants with sterile phenotypes, S-pool) were constructed for BSA resequencing. The selection of 24 single nucleotide polymorphisms (SNP) through high-throughput genotyping and the development insertion and deletion (InDel) markers were conducted to narrow down the candidate interval. The pentapeptide repeat (PPR) family genes and upregulated genes in restore line in the candidate interval were analysed by qRT-PCR. Results: The fertility investigation results showed that fertile and sterile separation ratio was consistent with 1:1. BSA resequencing technology, high-throughput SNP genotyping, and InDel markers were used to identify Rf2 locus on candidate interval of 1.48 Mb on Chromosome D05. Furthermore, it was quantified in this experiment that InDel markers co-segregated with Rf2 enhanced the selection of the restorer line. The qRT-PCR analysis revealed PPR family gene Gh_D05G3391 located in candidate interval had significantly lower expression than sterile and maintainer lines. In addition, utilization of anther RNA-Seq data of CMS-D8 identified that the expression level of Gh_D05G3374 encoding NB-ARC domain-containing disease resistance protein in restorer lines was significantly higher than that in sterile and maintainer lines. Conclusions: This study not only enabled us to precisely locate the restore gene Rf2 but also evaluated the utilization of InDel markers for marker assisted selection in the CMS-D8 Rf2 cotton breeding line. The results of this study provide an important foundation for further studies on the mapping and cloning of restorer genes.

Additional file 1: Table S1 The primers of qRT-PCR. Table S2 Genotype of SNP site. Table S3 The primers of InDel markers. Table S4 The FPKM of 67 genes in CMS-D8 system.

Background Utilization of heterosis has greatly improved the productivity of many crops worldwide. Understanding the potential molecular mechanism about how hybridization produces superior yield in upland cotton is critical for efficient breeding programs. Results In this study, high, medium, and low hybrids varying in the level of yield heterosis were screened based on field experimentation of different years and locations. Phenotypically, high hybrid produced a mean of 14% more seed cotton yield than its better parent. Whole-genome RNA sequencing of these hybrids and their four inbred parents was performed using different tissues of the squaring stage. Comparative transcriptomic differences in each hybrid parent triad revealed a higher percentage of differentially expressed genes (DEGs) in each tissue. Expression level dominance analysis identified majority of hybrids DEGs were biased towards parent like expressions. An array of DEGs involved in ATP and protein binding, membrane, cell wall, mitochondrion, and protein phosphorylation had more functional annotations in hybrids. Sugar metabolic and plant hormone signal transduction pathways were most enriched in each hybrid. Further, these two pathways had most mapped DEGs on known seed cotton yield QTLs. Integration of transcriptome, QTLs, and gene co-expression network analysis discovered genes Gh_A03G1024, Gh_D08G1440, Gh_A08G2210, Gh_A12G2183, Gh_D07G1312, Gh_D08G1467, Gh_A03G0889, Gh_A08G2199, and Gh_D05G0202 displayed a complex regulatory network of many interconnected genes. qRT-PCR of these DEGs was performed to ensure the accuracy of RNA-Seq data. Conclusions Through genome-wide comparative transcriptome analysis, the current study identified nine key genes and pathways associated with biological process of yield heterosis in upland cotton. Our results and data resources provide novel insights and will be useful for dissecting the molecular mechanism of yield heterosis in cotton.

Additional file 19: Table S6. Number of co-expressed modules and genes for hybrids.

Additional file 21: Table S8. List of gene primers used for qRT-PCR.

Key message This research based on RNA-seq, biochemical, and cytological analyses sheds that ROS may serve as important signaling molecules of cytoplasmic male sterility in CMS-D8 cotton. To understand the mechanism of cytoplasmic male sterility in cotton (Gossypium hirsutum), transcriptomic, cytological, and biochemical analysis were performed between the cytoplasmic male sterility CMS-D8 line, Zhong41A, and its maintainer line Zhong41B. A total of 2335 differentially expressed genes (DEGs) were identified in the CMS line at three different stages of anther development. Bioinformatics analysis of these DEGs indicated their relationship to reactive oxygen species (ROS) homeostasis, including reduction-oxidation reactions and the metabolism of glutathione and ascorbate. At the same time, DEGs associated with tapetum development, especially the transition to secretory tapetum, were down-regulated in the CMS line. Biochemical analysis indicated that the ability of the CMS line to eliminate ROS was decreased, which led to the rapid release of H2O2. Cytological analysis revealed that the most crucial defect in the CMS line was the abnormal tapetum. All these results are consistent with the RNA sequencing data. On the basis of our findings, we propose that ROS act as signal molecules, which are released from mitochondria and transferred to the nucleus, triggering the formation of abnormal tapetum.

Anther development in flowering plants is highly sensitive to high-temperature (HT) stress. Understanding the potential epigenetic mechanism of anther infertility induced by HT stress in cotton (Gossypium hirsutum L.) is crucial for the effective use of genetic resources to guide plant breeding. Using the whole-genome bisulfite sequencing, we map cytosine methylation at single-base resolution across the whole genome of cotton anthers, and changes in the methylome of the cytoplasmic male sterility system associated with HT stress were analysed in two cotton lines with contrasting HT stress tolerance. The cotton anther genome was found to display approximately 31.6%, 68.7%, 61.8%, and 21.8% methylation across all sequenced C sites and in the CG, CHG, and CHH sequence contexts, respectively. In an integrated global methylome and transcriptome analysis, only promoter-unmethylated genes showed higher expression levels than promoter-methylated genes, whereas gene body methylation presented an obvious positive correlation with gene expression. The methylation profiles of transposable elements in cotton anthers were characterized, and more differentially methylated transposable elements were demethylated under HT stress. HT-induced promoter methylation changes led to the up-regulation of the mitochondrial respiratory chain enzyme-associated genes GhNDUS7, GhCOX6A, GhCX5B2, and GhATPBM, ultimately promoting a series of redox processes to form ATP for normal anther development under HT stress. In vitro application of the common DNA methylation inhibitor 5-azacytidine and accelerator methyl trifluoromethanesulfonate demonstrated that DNA demethylation promoted anther development, while increased methylation only partially inhibited anther development under HT stress.

Anther development in flowering plants is highly sensitive to high-temperature (HT) stress. Understanding the potential epigenetic mechanism of anther infertility induced by HT stress in cotton (Gossypium hirsutum L.) is crucial for the effective use of genetic resources to guide plant breeding. Using the whole-genome bisulfite sequencing, we map cytosine methylation at single-base resolution across the whole genome of cotton anthers, and changes in the methylome of the cytoplasmic male sterility system associated with HT stress were analysed in two cotton lines with contrasting HT stress tolerance. The cotton anther genome was found to display approximately 31.6%, 68.7%, 61.8%, and 21.8% methylation across all sequenced C sites and in the CG, CHG, and CHH sequence contexts, respectively. In an integrated global methylome and transcriptome analysis, only promoter-unmethylated genes showed higher expression levels than promoter-methylated genes, whereas gene body methylation presented an obvious positive correlation with gene expression. The methylation profiles of transposable elements in cotton anthers were characterized, and more differentially methylated transposable elements were demethylated under HT stress. HT-induced promoter methylation changes led to the up-regulation of the mitochondrial respiratory chain enzyme-associated genes GhNDUS7, GhCOX6A, GhCX5B2, and GhATPBM, ultimately promoting a series of redox processes to form ATP for normal anther development under HT stress. In vitro application of the common DNA methylation inhibitor 5-azacytidine and accelerator methyl trifluoromethanesulfonate demonstrated that DNA demethylation promoted anther development, while increased methylation only partially inhibited anther development under HT stress.

Proportions of variance components for yield, yield components, and fiber quality traits based on 30 F1 hybrids together with their 11 parental lines of upland cotton tested at three locations in China during 2016â 2017. (DOCX 17 kb)

Predicted dominance-environment interaction effects for yield and yield component traits in three locations based on 2 years filed test. (DOCX 32 kb)

Climate statistics and planting arrangement in three locations of this study. (DOCX 15 kb)

Predicted dominance-environment interaction effects for fiber quality traits in three locations based on 2 years filed test. (DOCX 32 kb)

Background Cytoplasmic male sterility in flowering plants is a convenient way to use heterosis via hybrid breeding and may be restored by nuclear restorer-of-fertility (Rf) genes. In most cases, Rf genes encoded pentatricopeptide repeat (PPR) proteins and several Rf genes are present in clusters of similar Rf-PPR-like (RFL) genes. However, the Rf genes in cotton were not fully characterized until now. Results In total, 35 RFL genes were identified in G. hirsutum, 16 in G. arboreum, and 24 in G. raimondii. Additionally, four RFL-rich regions were identified; the RFL-rich region in Gh_D05 is the probable location of Rf-PPR genes in cotton and will be studied further in the future. Furthermore, an insertion sequence was identified in the promoter sequence of Gh_D05G3392 gene in the restorer line, as compared with the CMS-D2 line and maintainer lines. An InDel-R marker was then developed and could be used to distinguish the restorer line carrying Rf1 from other genotypes without the Rf1 allele. Conclusion In this study, genome-wide identification and analysis of RFL genes have identified the candidate Rf-PPR genes for CMS in Gossypium. The identification and analysis of RFL genes and sequence variation analysis will be useful for cloning Rf genes in the future and also for three-line hybrid breeding in cotton.

Background Cytoplasmic male sterility (CMS) in flowering plants is usually caused by incompatibility between mitochondrial and nuclear genomes, and can be restored by nuclear genes known as restorer-of-fertility (Rf). Although the CMS/Rf system is useful and convenient for economic production of commercial hybrid seed, the molecular mechanisms of CMS occurrence and fertility restoration in cotton are unclear. Results Here, a combined small RNA and transcriptome sequencing analysis was performed on floral buds at the meiosis stage in three-line hybrid cotton system, and differentially expressed microRNAs (DEMs) and their target genes were identified and further analyzed for a possible involvement in CMS and fertility restoration. Totally 10 and 30 differentially expressed miRNA-target gene pairs were identified in A-B and A-R comparison group, respectively. A putative regulatory network of CMS occurrence and fertility restoration-related miRNA-target pairs during anther development were then constructed. The RLM-RACE analysis showed that gra-miR7505b regulates a PPR gene (Gh_D05G3392) by cleaving precisely at the 643 nt and 748 nt sites. The further analysis indicated that the sequence variation in the binding regions of Gh_D05G3392 and Gh_D05G3356 may cause a lower cleavage efficiency of the PPR genes by miR7505b and miR7505 in R line, respectively, leading to the up-regulation of the PPR genes and fertility restoration. These results have established their genetic involvement in fertility restoration in the CMS-D2 system. Conclusion Our combined miRNA and transcriptome analysis in three-line hybrid cotton system provides new insights into the molecular mechanisms of CMS occurrence and fertility restoration, which will contribute to further hybrid breeding in cotton. Electronic supplementary material The online version of this article (10.1186/s12870-018-1446-7) contains supplementary material, which is available to authorized users.

DNA methylation is an important epigenetic modification involved in multiple biological processes. Altered methylation patterns have been reported to be associated with male sterility in some plants, but their role in cotton cytoplasmic male sterility (CMS) remains unclear. Here, integrated methylome and transcriptome analyses were conducted between the CMS-D2 line ZBA and its near-isogenic maintainer line ZB in upland cotton. More methylated cytosine sites (mCs) and higher methylation levels (MLs) were found among the three sequence contexts in ZB compared to ZBA. A total of 4568 differentially methylated regions (DMRs) and 2096 differentially methylated genes (DMGs) were identified. Among the differentially expressed genes (DEGs) associated with DMRs (DMEGs), 396 genes were upregulated and 281 genes were downregulated. A bioinformatics analysis of these DMEGs showed that hyper-DEGs were significantly enriched in the &ldquo, oxidative phosphorylation&rdquo, pathway. Further qRT-PCR validation indicated that these hypermethylated genes (encoding the subunits of mitochondrial electron transport chain (ETC) complexes I and V) were all significantly upregulated in ZB. Our biochemical data revealed a higher extent of H2O2 production but a lower level of adenosine triphosphate (ATP) synthesis in CMS-D2 line ZBA. On the basis of the above results, we propose that disrupted DNA methylation in ZBA may disrupt the homeostasis of reactive oxygen species (ROS) production and ATP synthesis in mitochondria, triggering a burst of ROS that is transferred to the nucleus to initiate programmed cell death (PCD) prematurely, ultimately leading to microspore abortion. This study illustrates the important role of DNA methylation in cotton CMS.

The RNA editing occurring in plant organellar genomes mainly involves the change of cytidine to uridine. This process involves a deamination reaction, with cytidine deaminase as the catalyst. Pentatricopeptide repeat (PPR) proteins with a C-terminal DYW domain are reportedly associated with cytidine deamination, similar to members of the deaminase superfamily. PPR genes are involved in many cellular functions and biological processes including fertility restoration to cytoplasmic male sterility (CMS) in plants. In this study, we identified 227 and 211 DYW deaminase-coding PPR genes for the cultivated tetraploid cotton species G. hirsutum and G. barbadense (2n = 4x = 52), respectively, as well as 126 and 97 DYW deaminase-coding PPR genes in the ancestral diploid species G. raimondii and G. arboreum (2n = 26), respectively. The 227 G. hirsutum PPR genes were predicted to encode 52-2016 amino acids, 203 of which were mapped onto 26 chromosomes. Most DYW deaminase genes lacked introns, and their proteins were predicted to target the mitochondria or chloroplasts. Additionally, the DYW domain differed from the complete DYW deaminase domain, which contained part of the E domain and the entire E+ domain. The types and number of DYW tripeptides may have been influenced by evolutionary processes, with some tripeptides being lost. Furthermore, a gene ontology analysis revealed that DYW deaminase functions were mainly related to binding as well as hydrolase and transferase activities. The G. hirsutum DYW deaminase expression profiles varied among different cotton tissues and developmental stages, and no differentially expressed DYW deaminase-coding PPRs were directly associated with the male sterility and restoration in the CMS-D2 system. Our current study provides an important piece of information regarding the structural and evolutionary characteristics of Gossypium DYW-containing PPR genes coding for deaminases and will be useful for characterizing the DYW deaminase gene family in cotton biology and breeding.

The cytoplasmic male sterility (CMS) system is convenient and efficient for hybrid seed production in Upland cotton (Gossypium hirsutum L.). However, it has not been widely used because of limited restorer lines carrying the restorer gene Rf1 in the CMS-D2 system. In this study, the fertility segregation in a backcross (BC8F1) population of 409 individuals and an F2 population of 695 plants confirmed that the fertility restoration was determined by one dominant restorer gene (Rf1). A sequence alignment showed that 13 Rf1-linked simple sequence repeat marker sequences were distributed on nine scaffolds of chromosome 9 in the sequenced D5 genome of G. raimondii Ulbrich. Ten pentotricopeptide repeat (PPR)-like genes were identified on two scaffolds, including Scaffold 333 where nine PPR-like genes were clustered in a region of about 160 kb. Among them, PPR-like gene Cotton_D_gene_10013437 was identified as the candidate for the Rf1 gene through a comparative sequence analysis of the homologous gene among sterile (A), maintainer (B) and restorer (R) lines, and co-segregation analysis. Compared with the non-restoring lines, the restorer had a 9-nucleotide (nt) insertion and a single nucleotide polymorphism (SNP) 8 nt upstream of the insertion at the 3′ untranslated regions (3′ UTRs) in this gene. A cleaved amplified polymorphic sequence (CAPS) marker named CAPS-R was developed from the SNP site using the restriction enzyme DraI, and was further used to track the restorer gene and its homozygous or heterozygous status in molecular breeding for restorer lines. A marker-assisted selection system using the Rf1-specific CAPS-R marker and a CMS-D2 cytoplasm-specific SCAR marker was established to distinguish the three-line hybrids from other genotypes.

Cytoplasmic male sterility (CMS) is a maternally inherited trait that fails to produce functional pollen grains. The CMS system is widely employed to facilitate the utilization of heterosis in major crops. However, little is known about the CMS associated genes in Upland cotton (Gossypium hirsutum). The objective of this study was to compare CMS cotton (CMS-D2) with the cytoplasm from G. harknessii and its isogenic maintainer line with the normal fertile Upland cotton cytoplasm to identify CMS-D2 specific gene(s) and to develop CMS-specific sequence characterized amplified region (SCAR) markers. Based on Southern blot analysis using 10 mitochondrial gene-specific probes (cob, cox2, atp6, atp9, nad3, cox3, atpA, cox1, nad6 and nad9), three probes (cox3, atpA, and nad6) revealed restriction fragment length polymorphisms (RFLP) between the CMS-D2 and its isogenic maintainer line. RT-PCR confirmed that the three genes were differentially expressed between the two lines. These results indicated that there existed structural and expression variations in the three genes when the mitochondrial D2 genome was transferred into Upland cotton. Genome walking and rapid amplification of cDNA ends (RACE) were further performed to analyze the sequences of these genes and their flanking regions. For cox3 and nad6, there was only one different nucleotide each in the gene regions between the two lines. Also some nucleotides upstream of the ATG codon were different. For atpA, the sequences downstream the atpA were significantly different between the two cytoplasmic lines. Furthermore, two nucleotides at the -4 and -5 position from ATG codon were also changed between the two cytoplasms (i.e., CG→AA), and this mutation also exists in RNA sequences. Interestingly, nine nucleotides (ATGCAACTA) were also inserted at the location of 899 bp upstream of ATG codon in the CMS line. The results suggest that the abnormal sequence and expression of atpA gene is associated with CMS expression in Upland cotton. According to the significant different sequences downstream the atpA gene, a CMS-D2 specific SCAR marker was developed. The CMS-specific PCR bands were verified for 10 cultivars containing either normal- or CMS-D2cytoplasm. This will allow quick and reliable identification of the cytoplasmic types of individual plants at the seedling stage, and assessment of the purity of F1 seed lots.