Genome-wide chromatin accessibility landscape and dynamics of transcription factor networks during ovule and fiber development in cotton

BACKGROUND

The development of cotton fiber is regulated by the orchestrated binding of regulatory proteins to cis-regulatory elements associated with developmental genes. The cis-trans regulatory dynamics occurred throughout the course of cotton fiber development are elusive. Here we generated genome-wide high-resolution DNase I hypersensitive sites (DHSs) maps to understand the regulatory mechanisms of cotton ovule and fiber development.

RESULTS

We generated DNase I hypersensitive site (DHS) profiles from cotton ovules at 0 and 3 days post anthesis (DPA) and fibers at 8, 12, 15, and 18 DPA. We obtained a total of 1185 million reads and identified a total of 199,351 DHSs through ~ 30% unique mapping reads. It should be noted that more than half of DNase-seq reads mapped multiple genome locations and were not analyzed in order to achieve a high specificity of peak profile and to avoid bias from repetitive genomic regions. Distinct chromatin accessibilities were observed in the ovules (0 and 3 DPA) compared to the fiber elongation stages (8, 12, 15, and 18 DPA). Besides, the chromatin accessibility during ovules was particularly elevated in genomic regions enriched with transposable elements (TEs) and genes in TE-enriched regions were involved in ovule cell division. We analyzed cis-regulatory modules and revealed the influence of hormones on fiber development from the regulatory divergence of transcription factor (TF) motifs. Finally, we constructed a reliable regulatory network of TFs related to ovule and fiber development based on chromatin accessibility and gene co-expression network. From this network, we discovered a novel TF, WRKY46, which may shape fiber development by regulating the lignin content.

CONCLUSIONS

Our results not only reveal the contribution of TEs in fiber development, but also predict and validate the TFs related to fiber development, which will benefit the research of cotton fiber molecular breeding.

Genome-Wide Association Study of Lint Percentage in Gossypium hirsutum L. Races

Lint percentage is one of the most essential yield components and an important economic index for cotton planting. Improving lint percentage is an effective way to achieve high-yield in cotton breeding worldwide, especially upland cotton (Gossypium hirsutum L.). However, the genetic basis controlling lint percentage has not yet been systematically understood. Here, we performed a genome-wide association mapping for lint percentage using a natural population consisting of 189 G. hirsutum accessions (188 accessions of G. hirsutum races and one cultivar TM-1). The results showed that 274 single-nucleotide polymorphisms (SNPs) significantly associated with lint percentage were detected, and they were distributed on 24 chromosomes. Forty-five SNPs were detected at least by two models or at least in two environments, and their 5 Mb up- and downstream regions included 584 makers related to lint percentage identified in previous studies. In total, 11 out of 45 SNPs were detected at least in two environments, and their 550 Kb up- and downstream region contained 335 genes. Through RNA sequencing, gene annotation, qRT-PCR, protein-protein interaction analysis, the cis-elements of the promotor region, and related miRNA prediction, Gh_D12G0934 and Gh_A08G0526 were selected as key candidate genes for fiber initiation and elongation, respectively. These excavated SNPs and candidate genes could supplement marker and gene information for deciphering the genetic basis of lint percentage and facilitate high-yield breeding programs of G. hirsutum ultimately.

The cotton protein GhIQD21 interacts with GhCaM7 and modulates organ morphogenesis in Arabidopsis by influencing microtubule stability

KEY MESSAGE

GhIQD21, a cotton IQ67-domain protein, interacts with GhCaM7 and alters organ shape in Arabidopsis by modulating microtubule stability. Calcium ion (Ca²⁺) and the calcium sensor calmodulin play crucial roles in the growth and development of plants. GhCaM7, a calmodulin in upland cotton (Gossypium hirsutum L.), is highly expressed in cotton fiber cells during the rapid elongation period and plays an important role in fiber cell development. In this study, we screened for GhCaM7-interacting proteins and identified GhIQD21, which contains a typical IQ67-domain. GhIQD21 was preferentially expressed at the fiber rapid elongation stage, and the protein localized to microtubules (MTs). Ectopic expression of GhIQD21 in Arabidopsis resulted in shorter leaves, petals, siliques, and plant height, thicker inflorescences, and more trichomes when compared with wild type (WT). Further investigation indicated that the morphogenesis of leaf epidermal cells and silique cells was altered. There was less consistency in the orientation of cortical microtubules in cotyledon and hypocotyl epidermal cells. Furthermore, compared with WT, transgenic seedling hypocotyls were more sensitive to oryzalin, a MT depolymerization drug. These results indicated that GhIQD21 is a GhCaM7-interacting protein located in MTs and that it plays a role in plant growth and potentially cotton fiber development. This study provides a foundation for further studies of the function and regulatory mechanism of GhIQD21 in fiber cell development.

Identification of the Flavone-Inducible Counter-Defense Genes and Their cis-Elements in Helicoverpa armigera

Flavone is widely found in plants and plays an important role in plant defense against pests. Many pests, such as Helicoverpa armigera, use flavone as a cue to upregulate counter-defense genes for detoxification of flavone. Yet the spectrum of the flavone-inducible genes and their linked cis-regulatory elements remains unclear. In this study, 48 differentially expressed genes (DEGs) were found by RNA-seq. These DEGs were mainly concentrated in the retinol metabolism and drug metabolism-cytochrome P450 pathways. Further in silico analysis of the promoter regions of 24 upregulated genes predicted two motifs through MEME and five previously characterized cis-elements including CRE, TRE, EcRE, XRE-AhR and ARE. Functional analysis of the two predicted motifs and two different versions of ARE (named ARE1 and ARE2) in the promoter region of the flavone-inducible carboxylesterase gene CCE001j verified that the two motifs and ARE2 are not responsible for flavone induction of H. armigera counter-defense genes, whereas ARE1 is a new xenobiotic response element to flavone (XRE-Fla) and plays a decisive role in flavone induction of CCE001j. This study is of great significance for further understanding the antagonistic interaction between plants and herbivorous insects.

Identification and analysis of isoflavone reductase gene family in Gossypium hirsutum L

Isoflavone reductase (IFR) is a key enzyme controlling isoflavone synthesis and widely involved in response to various stresses. In this study, the IFR genes in four Gossypium species and other 7 species were identified and analyzed in the whole genome, and the physicochemical properties, gene structures, cis-acting elements, chromosomal locations, collinearity relationships and expression patterns of IFR genes were systematically analyzed. 28, 28, 14 and 15 IFR genes were identified in Gossypium hirsutum, Gossypium barbadense, Gossypium arboreum and Gossypium raimondii, respectively, which were divided into five clades according to the evolutionary tree and gene structure. Collinear analysis showed that segmental duplication and whole genome duplication were the main driving forces in the process of evolution, and most genes underwent pure selection. Gene structure analysis showed that IFR gene family was relatively conserved. Cis-element analysis of promoter showed that most GhIFR genes contain cis-elements related to abiotic stresses and plant hormones. Analysis of GhIFR gene expression under different stresses showed that GhIFR genes were involved in the response to drought, salt, heat and cold stresses through corresponding network mechanisms, especially GhIFR9A. Phenotypic analysis after silencing GhIFR9A gene by VIGS was shown that GhIFR9A gene was involved in the response to salt stress. This study laid a foundation for the subsequent functional study of cotton IFR genes.

Genome-Wide Investigation and Co-Expression Network Analysis of SBT Family Gene in Gossypium

Subtilases (SBTs), which belong to the serine peptidases, control plant development by regulating cell wall properties and the activity of extracellular signaling molecules, and affect all stages of the life cycle, such as seed development and germination, and responses to biotic and abiotic environments. In this study, 146 Gossypium hirsutum, 138 Gossypium barbadense, 89 Gossypium arboreum and 84 Gossypium raimondii SBTs were identified and divided into six subfamilies. Cotton SBTs are unevenly distributed on chromosomes. Synteny analysis showed that the members of SBT1 and SBT4 were expanded in cotton compared to Arabidopsis thaliana. Co-expression network analysis showed that six Gossypium arboreum SBT gene family members were in a network, among which five SBT1 genes and their Gossypium hirsutum and Arabidopsis thaliana direct homologues were down-regulated by salt treatment, indicating that the co-expression network might share conserved functions. Through co-expression network and annotation analysis, these SBTs may be involved in the biological processes of auxin transport, ABA signal transduction, cell wall repair and root tissue development. In summary, this study provides valuable information for the study of SBT genes in cotton and excavates SBT genes in response to salt stress, which provides ideas for cotton breeding for salinity resistance.

Sphingosine Promotes Embryo Biomass in Upland Cotton: A Biochemical and Transcriptomic Analysis

Sphingolipids are essential membrane components and signal molecules, but their regulatory role in cotton embryo growth is largely unclear. In this study, we evaluated the effects of treatment with the sphingolipid synthesis inhibitor fumonisin B1 (FB1), the serine palmityl transferase (SPT) inhibitor myriocin, the SPT sphingolipid product DHS (d18:0 dihydrosphingosine), and the post-hydroxylation DHS product PHS (t18:0 phytosphingosine) on embryo growth in culture, and performed comparative transcriptomic analysis on control and PHS-treated samples. We found that FB1 could inhibit cotton embryo development. At the five-day ovule/embryo developmental stage, PHS was the most abundant sphingolipid. An SPT enzyme inhibitor reduced the fresh weight of embryos, while PHS had the opposite effect. The transcriptomic analysis identified 2769 differentially expressed genes (1983 upregulated and 786 downregulated) in the PHS samples. A large number of transcription factors were highly upregulated, such as zinc finger, MYB, NAC, bHLH, WRKY, MADS, and GRF in PHS-treated samples compared to controls. The lipid metabolism and plant hormone (auxin, brassinosteroid, and zeatin) related genes were also altered. Our findings provide target metabolites and genes for cotton seed improvement.

Transcriptomic analyses show that 24-epibrassinolide (EBR) promotes cold tolerance in cotton seedlings

In plants, brassinosteroids (BRs) are a class of steroidal hormones that are involved in numerous physiological responses. However, the function of BRs in cold tolerance in cotton has not been explored. In this study, cotton seedlings were treated with five concentrations (0, 0.05, 0.1, 0.2, 0.5 and 1.0 mg/L) of 24-Epibrassinolide (EBR) at 4°C. We measured the electrolyte leakage, malondialdehyde (MDA) content, proline content, and net photosynthesis rate (Pn) of the seedlings, which showed that EBR treatment increased cold tolerance in cotton in a dose-dependent manner, and that 0.2 mg/L is an optimum concentration for enhancing cold tolerance. The function of EBR in cotton cotyledons was investigated in the control 0 mg/L (Cold+water) and 0.2 mg/L (Cold+EBR) treatments using RNA-Seq. A total of 4,001 differentially expressed genes (DEGs), including 2,591 up-regulated genes and 1,409 down-regulated genes were identified. Gene Ontology (GO) and biochemical pathway enrichment analyses showed that EBR is involved in the genetic information process, secondary metabolism, and also inhibits abscisic acid (ABA) and ethylene (ETH) signal transduction. In this study, physiological experiments showed that EBR can increase cold tolerance in cotton seedlings, and the comprehensive RNA-seq data shed light on the mechanisms through which EBR increases cold tolerance in cotton seedlings.

Comprehensive analysis of WOX genes uncovers that WOX13 is involved in phytohormone-mediated fiber development in cotton

BACKGROUND

The WOX (WUSCHEL-RELATED HOMEOBOX) gene family encodes a class of transcription factors that are unique to green plants, where they are involved in regulating the development of plant tissues and organs by determining cell fate. Although the importance of the WOX gene is well known, there are few studies describing their functions in cotton.

RESULTS

In this study, 32 WOX genes were found in Gossypium hirsutum. Phylogenetic analysis showed that WOX proteins of cotton can be divided into three clades: the ancient, intermediate, and WUS clades. The number of WOX proteins in the WUS clade was greater than the sum of the proteins in the other two clades. Our analysis revealed that 20 GhWOX genes are distributed on 16 cotton chromosomes and that duplication events are likely to have contributed to the expansion of the GhWOX family. All GhWOX genes have introns, and each GhWOX protein contains multiple motifs. RNA-seq data and real-time PCR showed that GhWOX13 gene subfamily is specifically expressed at a high level in cotton fibers. We also identified putative GA, NAA, and BR response elements in the promoter regions of the GhWOX13 genes and GhWOX13 transcripts were significantly induced by GA, NAA, and BR.

CONCLUSIONS

Our data provides a useful resource for future studies on the functional roles of cotton WOX genes and shows that the GhWOX13 genes may influence cotton fiber development. Our results also provide an approach for identifying and characterizing WOX protein genes in other species.

An integrated omics analysis reveals molecular mechanisms that are associated with differences in seed oil content between Glycine max and Brassica napus

BACKGROUND

Rapeseed (Brassica napus L.) and soybean (Glycine max L.) seeds are rich in both protein and oil, which are major sources of biofuels and nutrition. Although the difference in seed oil content between soybean (~ 20%) and rapeseed (~ 40%) exists, little is known about its underlying molecular mechanism.

RESULTS

An integrated omics analysis was performed in soybean, rapeseed, Arabidopsis (Arabidopsis thaliana L. Heynh), and sesame (Sesamum indicum L.), based on Arabidopsis acyl-lipid metabolism- and carbon metabolism-related genes. As a result, candidate genes and their transcription factors and microRNAs, along with phylogenetic analysis and co-expression network analysis of the PEPC gene family, were found to be largely associated with the difference between the two species. First, three soybean genes (Glyma.13G148600, Glyma.13G207900 and Glyma.12G122900) co-expressed with GmPEPC1 are specifically enriched during seed storage protein accumulation stages, while the expression of BnPEPC1 is putatively inhibited by bna-miR169, and two genes BnSTKA and BnCKII are co-expressed with BnPEPC1 and are specifically associated with plant circadian rhythm, which are related to seed oil biosynthesis. Then, in de novo fatty acid synthesis there are rapeseed-specific genes encoding subunits β-CT (BnaC05g37990D) and BCCP1 (BnaA03g06000D) of heterogeneous ACCase, which could interfere with synthesis rate, and β-CT is positively regulated by four transcription factors (BnaA01g37250D, BnaA02g26190D, BnaC01g01040D and BnaC07g21470D). In triglyceride synthesis, GmLPAAT2 is putatively inhibited by three miRNAs (gma-miR171, gma-miR1516 and gma-miR5775). Finally, in rapeseed there was evidence for the expansion of gene families, CALO, OBO and STERO, related to lipid storage, and the contraction of gene families, LOX, LAH and HSI2, related to oil degradation.

CONCLUSIONS

The molecular mechanisms associated with differences in seed oil content provide the basis for future breeding efforts to improve seed oil content.