Integrated Methylome and Transcriptome Analysis between the CMS-D2 Line ZBA and Its Maintainer Line ZB in Upland Cotton

Crosstalk between 6-methyladenine and 4-methylcytosine in Geobacter sulfurreducens exposed to extremely low-frequency electromagnetic field

4-Methylcytosine (4mC) and 6-methyladenine (6mA) are the most prevalent types of DNA modifications in prokaryotes. However, whether there is crosstalk between 4mC and 6mA remain unknown. Here, methylomes and transcriptomes of Geobacter sulfurreducens exposed to different intensities of extremely low frequency electromagnetic fields (ELF-EMF) were investigated. Results showed that the second adenine of all the 5'-GTACAG-3' motif was modified to 6mA (M-6mA). For the other 6mA (O-6mA), the variation in their distance from the neighboring M-6mA increased with the intensity of ELF-EMF. Moreover, cytosine adjacent to O-6mA has a much higher probability of being modified to 4mC than cytosine adjacent to M-6mA, and the closer an unmodified cytosine is to 4mC, the higher the probability that the cytosine will be modified to 4mC. Furthermore, there was no significant correlation between DNA methylation and gene expression regulation. These results suggest a reference signal that goes from M-6mA to O-6mA to 4mC.

The Effects of DNA Methylation on Cytoplasmic Male Sterility in Sugar Beet

DNA methylation is widely found in higher plants and can control gene expression by regulation without changing the DNA sequence. In this study, the whole-genome methylation map of sugar beet was constructed by WGBS (whole-genome bisulfite sequencing) technology, and the results of WGBS were verified by bisulfite transformation, indicating that the results of WGBS technology were reliable. In addition, 12 differential methylation genes (DMGs) were identified, which were related to carbohydrate and energy metabolism, pollen wall development, and endogenous hormone regulation. Quantitative real-time PCR (qRT-PCR) showed that 75% of DMG expression levels showed negative feedback with methylation level, indicating that DNA methylation can affect gene expression to a certain extent. In addition, we found hypermethylation inhibited gene expression, which laid a foundation for further study on the molecular mechanism of DNA methylation at the epigenetic level in sugar beet male sterility.

Identification and characterization of circular RNAs involved in the fertility stability of cotton CMS-D2 restorer line under heat stress

BACKGROUND

As a vital type of noncoding RNAs, circular RNAs (circRNAs) play important roles in plant growth and development and stress response. However, little is known about the biological roles of circRNAs in regulating the stability of male fertility restoration for cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) cotton under high-temperature (HT) stress.

RESULTS

In this study, RNA-sequencing and bioinformatics analysis were performed on pollen grains of isonuclear alloplasmic near-isogenic restorer lines NH [N(Rf1rf1)] and SH [S(Rf1rf1)] with obvious differences in fertility stability under HT stress at two environments. A total of 967 circRNAs were identified, with 250 differentially expressed under HT stress. We confirmed the back-splicing sites of eight selected circRNAs using divergent primers and Sanger sequencing. Tissue-specific expression patterns of five differentially expressed circRNAs (DECs) were also verified by RT-PCR and qRT-PCR. Functional enrichment and metabolic pathway analysis revealed that the parental genes of DECs were significantly enriched in fertility-related biological processes such as pollen tube guidance and cell wall organization, as well as the Pentose and glucuronate interconversions, Steroid biosynthesis, and N-Glycan biosynthesis pathways. Moreover, we also constructed a putative circRNA-mediated competing endogenous RNA (ceRNA) network consisting of 21 DECs, eight predicted circRNA-binding miRNAs, and their corresponding 22 mRNA targets, especially the two ceRNA modules circRNA346-miR159a-MYB33 and circRNA484-miR319e-MYB33, which might play important biological roles in regulating pollen fertility stability of cotton CMS-D2 restorer line under HT stress.

CONCLUSIONS

Through systematic analysis of the abundance, characteristics and expression patterns of circRNAs, as well as the potential functions of their parent genes, our findings suggested that circRNAs and their mediated ceRNA networks acted vital biological roles in cotton pollen development, and might be also essential regulators for fertility stability of CMS-D2 restorer line under heat stress. This study will open a new door for further unlocking complex regulatory mechanisms underpinning the fertility restoration stability for CMS-D2 in cotton.

Heat-responsive microRNAs participate in regulating the pollen fertility stability of CMS-D2 restorer line under high-temperature stress

Anther development and pollen fertility of cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) restorer lines are susceptible to continuous high-temperature (HT) stress in summer, which seriously hinders the large-scale application of "three-line" hybrids in production. Here, integrated small RNA, transcriptome, degradome, and hormone profiling was performed to explore the roles of microRNAs (miRNAs) in regulating fertility stability in mature pollens of isonuclear alloplasmic near-isogenic restorer lines NH and SH under HT stress at two environments. A total of 211 known and 248 novel miRNAs were identified, of which 159 were differentially expressed miRNAs (DEMs). Additionally, 45 DEMs in 39 miRNA clusters (PmCs) were also identified, and most highly expressed miRNAs were significantly induced in SH under extreme HT, especially four MIR482 and six MIR6300 family miRNAs. PmC28 was located in the fine-mapped interval of the Rf1 gene and contained two DEMs, gra-miR482_L-2R + 2 and gma-miR2118a-3p_R + 1_1ss18TG. Transcriptome sequencing identified 6281 differentially expressed genes, of which heat shock protein (HSP)-related genes, such as HSP70, HSP22, HSP18.5-C, HSP18.2 and HSP17.3-B, presented significantly reduced expression levels in SH under HT stress. Through integrating multi-omics data, we constructed a comprehensive molecular network of miRNA-mRNA-gene-KEGG containing 35 pairs of miRNA/target genes involved in regulating the pollen development in response to HT, among which the mtr-miR167a_R + 1, tcc-miR167c and ghr-miR390a, tcc-miR396c_L-1 and ghr-MIR169b-p3_1ss6AG regulated the pollen fertility by influencing ARF8 responsible for the auxin signal transduction, ascorbate and aldarate metabolism, and the sugar and lipid metabolism and transport pathways, respectively. Further combination with hormone analysis revealed that HT-induced jasmonic acid signaling could activate the expression of downstream auxin synthesis-related genes and cause excessive auxin accumulation, followed by a cascade of auxin signal transduction, ultimately resulting in pollen abortion. The results provide a new understanding of how heat-responsive miRNAs regulate the stability of fertility restoration for CMS-D2 cotton under heat stress.

Dose effects of restorer gene modulate pollen fertility in cotton CMS-D2 restorer lines via auxin signaling and flavonoid biosynthesis

KEY MESSAGE

Dose effects of Rf1 gene regulated retrieval mechanism of pollen fertility for CMS-D2 cotton. Cytoplasmic male sterility conditioned by Gossypium harknessii cytoplasm (CMS-D2) is an economical pollination control system for producing hybrid cotton seeds compared to 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 in China. 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 at two environments using pollen grains of 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. In addition, 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.

ROS accumulation-induced tapetal PCD timing changes leads to microspore abortion in cotton CMS lines

BACKGROUND

Cytoplasmic male sterility (CMS) is the basis of heterosis exploitation. CMS has been used to hybrid production in cotton, but its molecular mechanism remains unclear. CMS is associated with advanced or delayed tapetal programmed cell death (PCD), and reactive oxygen species (ROS) may mediate this process. In this study, we obtained Jin A and Yamian A, two CMS lines with different cytoplasmic sources.

RESULTS

Compared with maintainer Jin B, Jin A anthers showed advanced tapetal PCD with DNA fragmentation, producing excessive ROS which accumulated around the cell membrane, intercellular space and mitochondrial membrane. The activities of peroxidase (POD) and catalase (CAT) enzymes which can scavenge ROS were significantly decreased. However, Yamian A tapetal PCD was delayed with lower ROS content, and the activities of superoxide dismutase (SOD) and POD were higher than its maintainer. These differences in ROS scavenging enzyme activities may be caused by isoenzyme gene expressions. In addition, we found the excess ROS generated in Jin A mitochondria and ROS overflow from complex III might be the source in parallel with the reduction of ATP content.

CONCLUSION

ROS accumulation or abrogation were mainly caused by the joint action of ROS generation and scavenging enzyme activities transformation, which led to the abnormal progression of tapetal PCD, affected the development of microspores, and eventually contributed to male sterility. In Jin A, tapetal PCD in advance might be caused by mitochondrial ROS overproduction, accompanied by energy deficiency. The above studies will provide new insights into the cotton CMS and guide the follow-up research ideas.

Integrated analysis of metabolome and transcriptome reveals the cytoplasmic effects of CMS-D2 on pollen fertility resulting from disrupted lipid metabolism

Using cytoplasmic male sterility of Gossypium harknesii (CMS-D2) is an economical and effective method to produce cotton hybrids. However, the detrimental cytoplasmic effects of CMS-D2 on pollen fertility and fiber yields greatly limit the further development of three-line hybrid cotton in China. In this study, an integrated non-targeted metabolomics and transcriptome analysis was performed on mature pollens of maintainer line NB, isonuclear alloplasmic near-isogenic restorer lines NH and SH under two environments. A total of 820 metabolites were obtained, of which lipids and lipid-like molecules were the most, followed by organic acids derivatives, phenylpropanoids, and polyketides. Transcriptome analysis revealed significantly more differentially expressed genes (DEGs) in SH versus NH both in Anyang and Jiujiang, and most of the DEGs were significantly upregulated. Further KEGG analysis showed that most DEGs were enriched in the biosynthesis of unsaturated fatty acids, phenylalanine metabolism, and phagosome. Based on the weighted gene co-expression network analysis, totally 74 hub genes were also identified, of which three transcription factors, i.e., WRKY22, WRKY53, and ARF18 were significantly upregulated in SH and may play a negative regulatory role in pollen development by directly or indirectly regulating the jasmonic acid synthesis and signal transduction. Moreover, we found that the negative effects of CMS-D2 cytoplasm on pollen fertility were mainly due to disturbed lipid metabolism, especially the metabolic balance of unsaturated fatty acids, ultimately resulting in the decline of pollen fertility. Meanwhile, in the presence of CMS-D2 sterile cytoplasm, the cytoplasmic-nucleus interaction effects generated a substantial quantity of flavonoids involved in the fertility restoration process. This study preliminarily clarified some of the reasons for the negative effects of CMS-D2 cytoplasm on pollen fertility, and our results will provide an important theoretical reference for further breeding and improvement of three-line hybrid cotton in the future.

Single-molecule real-time transcript sequencing of developing cotton anthers facilitates genome annotation and fertility restoration candidate gene discovery

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.

Identification and functional analysis of a pollen fertility-associated gene GhGLP4 of Gossypium hirsutum L

KEY MESSAGE

Cotton male fertility-associated gene GhGLP4, encoding a germin-like protein, is essential for anthers development by keeping ROS homeostasis through reducing H2O2 level. Utilization of heterosis is an important way to increase cotton yield and improve fiber quality in hybrid cotton development programs. Male sterility is used in the development of cotton hybrids to reduce the cost of hybrid seed production by eliminating the process of emasculation. From the transcriptome analysis of genic male sterile mutant (ms1) and its background C312 of G. hirsutum, a gene encoding germin-like protein (GhGLP4) was found significantly down-regulated in different developmental stages of ms1 anthers. To explore the gene function in cotton fertility, GhGLP4 was further studied and interfered by virus-induced gene silencing. In the GhGLP4 interfered cotton lines, the expression level of GhGLP4 was significantly decreased in the stamens, and the down-regulation of GhGLP4 resulted in pollen sac closure, stigma exertion, filament shortening, decrease in the number of anthers and complete male sterility. The expression levels of respiratory burst oxidase homologs (Rboh, NADPH oxidase) were significantly altered. Further investigation showed that the SOD activity decreased while the H2O2 content increased in the atypical stamens. These results indicated that GhGLP4 gene affected the cotton anther development through maintenance of ROS homeostasis by H2O2 reduction.

Physical mapping and InDel marker development for the restorer gene Rf2 in cytoplasmic male sterile CMS-D8 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 Rf₂ 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 Rf₂ 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 Rf₂ locus on candidate interval of 1.48 Mb on chromosome D05. Furthermore, it was quantified in this experiment that InDel markers co-segregated with Rf₂ 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 Rf₂ but also evaluated the utilization of InDel markers for marker assisted selection in the CMS-D8 Rf₂ cotton breeding line. The results of this study provide an important foundation for further studies on the mapping and cloning of restorer genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07342-y.