High-resolution temporal dynamic transcriptome landscape reveals a GhCAL-mediated flowering regulatory pathway in cotton (Gossypium hirsutum L.)

GhGASA14 regulates the flowering time of upland cotton in response to GA3

KEY MESSAGE

The silencing of GhGASA14 and the identification of superior allelic variation in its coding region indicate that GhGASA14 may positively regulate flowering and the response to GA3. Gibberellic acid-stimulated Arabidopsis (GASA), a member of the gibberellin-regulated short amino acid family, has been extensively investigated in several plant species and found to be critical for plant growth and development. However, research on this topic in cotton has been limited. In this study, we identified 38 GhGASAs that were dispersed across 18 chromosomes in upland cotton, and all of these genes had a GASA core domain. Transcriptome expression patterns and qRT-PCR results revealed that GhGASA9 and GhGASA14 exhibited upregulated expression not only in the floral organs but also in the leaves of early-maturing cultivars. The two genes were functionally characterized by virus-induced gene silencing (VIGS), and the budding and flowering times after silencing the target genes were later than those of the control (TRV:00). Compared with that in the water-treated group (MOCK), the flowering period of the different fruiting branches in the GA3-treated group was more concentrated. Interestingly, allelic variation was detected in the coding sequence of GhGASA14 between early-maturing and late-maturing accessions, and the frequency of this favorable allele was greater in high-latitude cotton cultivars than in low-latitude ones. Additionally, a significant linear relationship was observed between the expression level of GhGASA14 and flowering time among the 12 upland cotton accessions. Taken together, these results indicated that GhGASA14 may positively regulate flowering time and respond to GA3. These findings could lead to the use of valuable genetic resources for breeding early-maturing cotton cultivars in the future.

A cotton mitochondrial alternative electron transporter, GhD2HGDH, induces early flowering by modulating GA and photoperiodic pathways

D-2-hydroxyglutarate dehydrogenase (D2HGDH) is a mitochondrial enzyme containing flavin adenine dinucleotide FAD, existing as a dimer, and it facilitates the specific oxidation of D-2HG to 2-oxoglutarate (2-OG), which is a key intermediate in the tricarboxylic acid (TCA) cycle. A Genome-wide expression analysis (GWEA) has indicated an association between GhD2HGDH and flowering time. To further explore the role of GhD2HGDH, we performed a comprehensive investigation encompassing phenotyping, physiology, metabolomics, and transcriptomics in Arabidopsis thaliana plants overexpressing GhD2HGDH. Transcriptomic and qRT-PCR data exhibited heightened expression of GhD2HGDH in upland cotton flowers. Additionally, early-maturing cotton exhibited higher expression of GhD2HGDH across all tissues than delayed-maturing cotton. Subcellular localization confirmed its presence in the mitochondria. Overexpression of GhD2HGDH in Arabidopsis resulted in early flowering. Using virus-induced gene silencing (VIGS), we investigated the impact of GhD2HGDH on flowering in both early- and delayed-maturing cotton plants. Manipulation of GhD2HGDH expression levels led to changes in photosynthetic pigment and gas exchange attributes. GhD2HGDH responded to gibberellin (GA3) hormone treatment, influencing the expression of GA biosynthesis genes and repressing DELLA genes. Protein interaction studies, including yeast two-hybrid, luciferase complementation (LUC), and GST pull-down assays, confirmed the interaction between GhD2HGDH and GhSOX (Sulfite oxidase). The metabolomics analysis demonstrated GhD2HGDH's modulation of the TCA cycle through alterations in various metabolite levels. Transcriptome data revealed that GhD2HGDH overexpression triggers early flowering by modulating the GA3 and photoperiodic pathways of the flowering core factor genes. Taken together, GhD2HGDH positively regulates the network of genes associated with early flowering pathways.

Multistage temporal transcriptomic atlas unveils major contributor to reproductive phase in upland cotton

Flowering is a major developmental transition in plants, but asynchronous flowering hinders the utilization of wild cotton relatives in breeding programs. We performed comparative transcriptomic profiling of early- and late-flowering Gossypium hirsutum genotypes to elucidate genetic factors influencing reproductive timing. Shoot apices were sampled from the photoperiod-sensitive landrace G. hirsutum purpurascens (GhP) and early-maturing variety ZhongMianSuo (ZMS) at five time points following the emergence of sympodial nodes. RNA-sequencing revealed extensive transcriptional differences during floral transition. Numerous flowering-associated genes exhibited genotype-specific expression, including FLOWERING LOCUS T (FT) homologs upregulated in ZMS. FT-interacting factors like SOC1 and CO-like also showed higher expression in ZMS, implicating florigen pathways in early flowering. Additionally, circadian clock and light signalling components were misregulated between varieties, suggesting altered photoperiod responses in GhP. Weighted co-expression network analysis specifically linked a module enriched for circadian-related genes to GhP's late flowering. Through an integrated transcriptome analysis, we defined a regulatory landscape of reproductive phase change in cotton. Differentially expressed genes related to photoperiod, circadian clock, and light signalling likely contribute to delayed flowering in wild cottons. Characterization of upstream flowering regulators will enable modifying photoperiod sensitivity and expand germplasm use for cotton improvement. This study provides candidate targets for elucidating interactive mechanisms that control cotton flowering time across diverse genotypes.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

GhAP1-D3 positively regulates flowering time and early maturity with no yield and fiber quality penalties in upland cotton

Flowering time (FTi) is a major factor determining how quickly cotton plants reach maturity. Early maturity greatly affects lint yield and fiber quality and is crucial for mechanical harvesting of cotton in northwestern China. Yet, few quantitative trait loci (QTLs) or genes regulating early maturity have been reported in cotton, and the underlying regulatory mechanisms are largely unknown. In this study, we characterized 152, 68, and 101 loci that were significantly associated with the three key early maturity traits-FTi, flower and boll period (FBP) and whole growth period (WGP), respectively, via four genome-wide association study methods in upland cotton (Gossypium hirsutum). We focused on one major early maturity-related genomic region containing three single nucleotide polymorphisms on chromosome D03, and determined that GhAP1-D3, a gene homologous to Arabidopsis thaliana APETALA1 (AP1), is the causal locus in this region. Transgenic plants overexpressing GhAP1-D3 showed significantly early flowering and early maturity without penalties for yield and fiber quality compared to wild-type (WT) plants. By contrast, the mutant lines of GhAP1-D3 generated by genome editing displayed markedly later flowering than the WT. GhAP1-D3 interacted with GhSOC1 (SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1), a pivotal regulator of FTi, both in vitro and in vivo. Changes in GhAP1-D3 transcript levels clearly affected the expression of multiple key flowering regulatory genes. Additionally, DNA hypomethylation and high levels of H3K9ac affected strong expression of GhAP1-D3 in early-maturing cotton cultivars. We propose that epigenetic modifications modulate GhAP1-D3 expression to positively regulate FTi in cotton through interaction of the encoded GhAP1 with GhSOC1 and affecting the transcription of multiple flowering-related genes. These findings may also lay a foundation for breeding early-maturing cotton varieties in the future.

Recent advances and future perspectives in early-maturing cotton research

Cotton's fundamental requirements for long periods of growth and specific seasonal temperatures limit the global arable areas that can be utilized to cultivate cotton. This constraint can be alleviated by breeding for early-maturing varieties. By delaying the sowing dates without impacting the boll-opening time, early-maturing varieties not only mitigate the yield losses brought on by unfavorable weathers in early spring and late autumn but also help reducing the competition between cotton and other crops for arable land, thereby optimizing the cropping system. This review presents studies and breeding efforts for early-maturing cotton, which efficiently pyramid early maturity, high-quality, multiresistance traits, and suitable plant architecture by leveraging pleiotropic genes. Attempts are also made to summarize our current understanding of the molecular mechanisms underlying early maturation, which involves many pathways such as epigenetic, circadian clock, and hormone signaling pathways. Moreover, new avenues and effective measures are proposed for fine-scale breeding of early-maturing crops to ensure the healthy development of the agricultural industry.

QTL mapping and BSA-seq map a major QTL for the node of the first fruiting branch in cotton

Understanding the genetic basis of the node of the first fruiting branch (NFFB) improves early-maturity cotton breeding. Here we report QTL mapping on 200 F₂ plants and derivative F_2:3 and F_2:4 populations by genotyping by sequencing (GBS). BC₁F₂ population was constructed by backcrossing one F_2:4 line with the maternal parent JF914 and used for BSA-seq for further QTL mapping. A total of 1,305,642 SNPs were developed between the parents by GBS, and 2,907,790 SNPs were detected by BSA-seq. A high-density genetic map was constructed containing 11,488 SNPs and spanning 4,202.12 cM in length. A total of 13 QTL were mapped in the 3 tested populations. JF914 conferred favorable alleles for 11 QTL, and JF173 conferred favorable alleles for the other 2 QTL. Two stable QTL were repeatedly mapped in F_2:3 and F_2:4, including qNFFB-D3-1 and qNFFB-D6-1. Only qNFFB-D3-1 contributed more than 10% of the phenotypic variation. This QTL covered about 24.7 Mb (17,130,008-41,839,226 bp) on chromosome D3. Two regions on D3 (41,779,195-41,836,120 bp, 41,836,768-41,872,287 bp) were found by BSA-seq and covered about 92.4 Kb. This 92.4 Kb region overlapped with the stable QTL qNFFB-D3-1 and contained 8 annotated genes. By qRT-PCR, Ghir_D03G012430 showed a lower expression level from the 1- to 2-leaf stage and a higher expression level from the 3- to 6-leaf stage in the buds of JF173 than that of JF914. Ghir_D03G012390 reached the highest level at the 3- and 5-leaf stages in the buds of JF173 and JF914, respectively. As JF173 has lower NFFB and more early maturity than JF914, these two genes might be important in cell division and differentiation during NFFB formation in the seedling stage. The results of this study will facilitate a better understanding of the genetic basis of NFFB and benefit cotton molecular breeding for improving earliness traits.

Identification of Differentially Expressed Genes Related to Floral Bud Differentiation and Flowering Time in Three Populations of Lycoris radiata

The transition from vegetative to reproductive growth is important for controlling the flowering of Lycoris radiata. However, the genetic control of this complex developmental process remains unclear. In this study, 18 shoot apical meristem (SAM) samples were collected from early-, mid- and late-flowering populations during floral bud differentiation. The histological analysis of paraffin sections showed that the floral bud differentiation could be divided into six stages; the differentiation time of the early group was earlier than that of the middle and late groups, and the late group was the latest. In different populations, some important differential genes affecting the flowering time were identified by transcriptome profiles of floral bud differentiation samples. Weighted gene co-expression network analysis (WGCNA) was performed to enrich the gene co-expression modules of diverse flowering time populations (FT) and floral bud differentiation stages (ST). In the MEyellow module, five core hub genes were identified, including CO14, GI, SPL8, SPL9, and SPL15. The correlation network of hub genes showed that they interact with SPLs, AP2, hormone response factors (auxin, gibberellin, ethylene, and abscisic acid), and several transcription factors (MADS-box transcription factor, bHLH, MYB, and NAC3). It suggests the important role of these genes and the complex molecular mechanism of floral bud differentiation and flowering time in L. radiata. These results can preliminarily explain the molecular mechanism of floral bud differentiation and provide new candidate genes for the flowering regulation of Lycoris.

Functional divergence of GhAP1.1 and GhFUL2 associated with flowering regulation in upland cotton (Gossypium hirsutum L.)

The AP1/FUL transcription factors are important for floral development, but the underlying molecular mechanisms remain unclear. In this study, we cloned and identified two AP1/FUL-like genes, GhAP1.1 and GhFUL2, in upland cotton, which is a commonly cultivated economically valuable crop. Sequence alignment and phylogenetic analysis indicated that GhAP1.1 and GhFUL2, which are encoded by genes in the AP1/FUL clade, have conserved N-terminal regions but diverse C-terminal domains. Quantitative real-time PCR analysis revealed that GhAP1.1 and GhFUL2 were expressed in the flower and root, and showed opposite expression patterns during shoot apical meristem development. The upregulated expression of GhAP1.1 in Arabidopsis did not result in significant changes to the flowering time or floral organ development, and the transcript levels of the florigen FT increased and those of LFY decreased. Overexpression of GhFUL2 in Arabidopsis delayed flowering and promoted bolting by decreasing FT and LFY transcript levels. Silencing GhFUL2 in cotton dramatically increased the expression of GhFT and GhAP1.3 and promoted flowering. Yeast two-hybrid and bimolecular fluorescence complementation assays indicated that GhAP1.1 could interact with the SVP homolog GhSVP2.2, whereas GhFUL2 formed heterodimers with GhSEP3/GhSEP4 homologs and GhSVP2.2. The present results demonstrated that the functional divergence of GhAP1.1 and GhFUL2, which involved changes in sequences and expression patterns, influenced the regulation of cotton flower development.

Genomewide Identification and Characterization of the Genes Involved in the Flowering of Cotton

Flowering is a prerequisite for flowering plants to complete reproduction, and flowering time has an important effect on the high and stable yields of crops. However, there are limited reports on flowering-related genes at the genomic level in cotton. In this study, genomewide analysis of the evolutionary relationship of flowering-related genes in different cotton species shows that the numbers of flowering-related genes in the genomes of tetraploid cotton species Gossypium hirsutum and Gossypium barbadense were similar, and that these numbers were approximately twice as much as the number in diploid cotton species Gossypium arboretum. The classification of flowering-related genes shows that most of them belong to the photoperiod and circadian clock flowering pathway. The distribution of flowering-related genes on the chromosomes of the At and Dt subgenomes was similar, with no subgenomic preference detected. In addition, most of the flowering-related core genes in Arabidopsis thaliana had homologs in the cotton genome, but the copy numbers and expression patterns were disparate; moreover, flowering-related genes underwent purifying selection throughout the evolutionary and selection processes. Although the differentiation and reorganization of many key genes of the cotton flowering regulatory network occurred throughout the evolutionary and selection processes, most of them, especially those involved in the important flowering regulatory networks, have been relatively conserved and preferentially selected.

Integrative Analyses of Transcriptomes and Metabolomes Reveal Associated Genes and Metabolites with Flowering Regulation in Common Vetch (Vicia sativa L.)

As an important source of protein for livestock and human consumption, Vicia sativa is cultivated worldwide, but its seed production is hampered at high altitudes because of the short frost-free period. Flowering represents the transition from a vegetative to a reproductive period, and early flowering benefits plant seed production at high altitudes. However, the molecular mechanisms of flowering regulation in V. sativa remain elusive. In the present study, two V. sativa accessions with different flowering characteristics were used: Lan3 (early-flowering) was cultivated by our laboratory, and 503 (late-flowering) was selected from 222 V. sativa accessions after three years of field experiments. The shoot samples (shoot tip length = 10 cm) of these two accessions were collected 63, 70, and 77 days after sowing, and the molecular regulatory mechanism of the flowering process was identified by integrative analyses of the transcriptomes and metabolomes. Kyoto Encyclopedia of Genes and Genomes enrichment showed that the synthesis and signal transduction of plant hormone pathways were the most enriched pathways in 4274 differentially expressed genes (DEGs) and in 259 differential metabolites between Lan3 and 503. Moreover, the contents of three metabolites related to salicylic acid biosynthesis and the transcription levels of two DEGs related to salicylic acid signal transduction in Lan3 were higher than those in 503. Further verification in various accessions indicated that salicylic acid metabolism may be involved in the flowering regulation process of V. sativa. These findings provide valuable information for understanding the flowering mechanism and for promoting breeding research in V. sativa.

Transcriptome analysis and identification of genes associated with oil accumulation in upland cotton

Cotton is not only the most important fiber crop but also the fifth most important oilseed crop in the world because of its oil-rich seeds as a byproduct of fiber production. By comparative transcriptome analysis between two germplasms with diverse oil accumulation, we reveal pieces of the gene expression network involved in the process of oil synthesis in cottonseeds. Approximately, 197.16 Gb of raw data from 30 RNA sequencing samples with 3 biological replicates were generated. Comparison of the high-oil and low-oil transcriptomes enabled the identification of 7682 differentially expressed genes (DEGs). Based on gene expression profiles relevant to triacylglycerol (TAG) biosynthesis, we proposed that the Kennedy pathway (diacylglycerol acyltransferase-catalyzed diacylglycerol to TAG) is the main pathway for oil production, rather than the phospholipid diacylglycerol acyltransferase-mediated pathway. Using weighted gene co-expression network analysis, 5312 DEGs were obtained and classified into 14 co-expression modules, including the MEblack module containing 10 genes involved in lipid metabolism. Among the DEGs in the MEblack module, GhCYSD1 was identified as a potential key player in oil biosynthesis. The overexpression of GhCYSD1 in yeast resulted in increased oil content and altered fatty acid composition. This study may not only shed more light on the underlying molecular mechanism of oil accumulation in cottonseed oil, but also provide a set of new gene for potential enhancement of oil content in cottonseeds.

Comparative transcriptomic analysis on compatible/incompatible grafts in citrus

Grafting is a useful cultivation technology to resist abiotic and biotic stresses and is an integral part of citrus production. However, some widely utilized rootstocks may still exhibit graft incompatibility in the orchard. "Hongmian miyou" (Citrus maxima (Burm.) Merrill) is mutated from "Guanxi miyou", but these two scions showed different compatibility with available Poncirus trifoliata rootstock. Foliage etiolation is an observed symptom of graft incompatibility, but its mechanism remains poorly understood. This study is the first to investigate the morphological, physiological, and anatomical differences between the compatible/incompatible grafts, and perform transcriptome profiling at crucial stages of the foliage etiolation process. Based on the comprehensive analyses, hormonal balance was disordered, and two rate-limiting genes, NCED3 (9-cis-epoxycarotenoid dioxygenases 3) and NCED5, being responsible for ABA (abscisic acid) accumulation, were highlighted. Further correlation analysis indicated that IAA (indole-3-acetic acid) and ABA were the most likely inducers of the expression of stresses-related genes. In addition, excessive starch accumulation was observed in lamina and midribs of incompatible grafts leaves. These results provided a new insight into the role of the hormonal balance and abscisic acid biosynthesis genes in regulation and contribution to the graft incompatibility, and will further define and deploy candidate genes to explore the mechanisms underlying citrus rootstock- scion interactions.

Insights Into MicroRNA-Mediated Regulation of Flowering Time in Cotton Through Small RNA Sequencing

The timing of flowering is a key determinant for plant reproductive. It has been demonstrated that microRNAs (miRNAs) play an important role in transition from the vegetative to reproductive stage in cotton; however, knowledge remains limited about the regulatory role of miRNAs involved in flowering time regulation in cotton. To elucidate the molecular basis of miRNAs in response to flowering time in cotton, we performed high-throughput small RNA sequencing at the fifth true leaf stage. We identified 56 and 43 miRNAs that were significantly up- and downregulated in two elite early flowering cultivars (EFC) compared with two late flowering cultivars (LFC), respectively. The miRNA targets by RNA sequencing analysis showed that GhSPL4 in SBP transcription factor family targeted by GhmiR156 was significantly upregulated in EFCs. Co-expression regulatory network analysis (WGCNA) revealed that GhSOC1, GhAP1, GhFD, GhCOL3, and GhAGL16 act as node genes in the auxin- and gibberellin-mediated flowering time regulatory networks in cotton. Therefore, elucidation of miRNA-mediated flowering time regulatory network will contribute to our understanding of molecular mechanisms underlying flowering time in cotton.

QTL and candidate gene identification of the node of the first fruiting branch (NFFB) by QTL-seq in upland cotton (Gossypium hirsutum L.)

Background

The node of the first fruiting branch (NFFB) is an important precocious trait in cotton. Many studies have been conducted on the localization of quantitative trait loci (QTLs) and genes related to fiber quality and yield, but there has been little attention to traits related to early maturity, especially the NFFB, in cotton.

Results

To identify the QTL associated with the NFFB in cotton, a BC₄F₂ population comprising 278 individual plants was constructed. The parents and two DNA bulks for high and low NFFB were whole genome sequenced, and 243.8 Gb of clean nucleotide data were generated. A total of 449,302 polymorphic SNPs and 135,353 Indels between two bulks were identified for QTL-seq. Seventeen QTLs were detected and localized on 11 chromosomes in the cotton genome, among which two QTLs (qNFFB-Dt2–1 and qNFFB-Dt3–3) were located in hotspots. Two candidate genes (GhAPL and GhHDA5) related to the NFFB were identified using quantitative real-time PCR (qRT-PCR) and virus-induced gene silencing (VIGS) experiments in this study. Both genes exhibited higher expression levels in the early-maturing cotton material RIL182 during flower bud differentiation, and the silencing of GhAPL and GhHDA5 delayed the flowering time and increased the NFFB compared to those of VA plants in cotton.

Conclusions

Our study preliminarily found that GhAPL and GhHDA5 are related to the early maturity in cotton. The findings provide a basis for the further functional verification of candidate genes related to the NFFB and contribute to the study of early maturity in cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08164-2.

Transcriptomic and Metabolomic Analyses Provide Insights Into an Aberrant Tissue of Tea Plant (Camellia sinensis)

Tea plant (Camellia sinensis (L.) O. Kuntze) is one of the most important economic crops with multiple mutants. Recently, we found a special tea germplasm that has an aberrant tissue on its branches. To figure out whether this aberrant tissue is associated with floral bud (FB) or dormant bud (DB), we performed tissue section, transcriptome sequencing, and metabolomic analysis of these tissues. Longitudinal sections indicated the aberrant tissue internal structure was more like a special bud (SB), but was similar to that of DB. Transcriptome data analysis showed that the number of heterozygous and homozygous SNPs was significantly different in the aberrant tissue compared with FB and DB. Further, by aligning the unmapped sequences of the aberrant tissue to the Non-Redundant Protein Sequences (NR) database, we observed that 36.13% of unmapped sequences were insect sequences, which suggested that the aberrant tissue might be a variation of dormant bud tissue influenced by the interaction of tea plants and insects or pathogens. Metabolomic analysis showed that the differentially expressed metabolites (DEMs) between the aberrant tissue and DB were significantly enriched in the metabolic pathways of biosynthesis of plant hormones and biosynthesis of phenylpropanoids. Subsequently, we analyzed the differentially expressed genes (DEGs) in the above mentioned two tissues, and the results indicated that photosynthetic capacity in the aberrant tissue was reduced, whereas the ethylene, salicylic acid and jasmonic acid signaling pathways were activated. We speculated that exogenous infection induced programmed cell death (PCD) and increased the lignin content in dormant buds of tea plants, leading to the formation of this aberrant tissue. This study advanced our understanding of the interaction between plants and insects or pathogens, providing important clues about biotic stress factors and key genes that lead to mutations and formation of the aberrant tissue.

Genome-Wide Identification and Functional Investigation of 1-Aminocyclopropane-1-carboxylic Acid Oxidase (ACO) Genes in Cotton

ACO is one of the rate-limiting enzymes in the biosynthesis of ethylene, and it plays a critical role in the regulation of plant growth and development. However, the function of ACO genes in cotton is not well studied. In this study, a total of 332 GhACOs, 187 GaACOs, and 181 GrACOs were identified in G. hirsutum, G. arboretum, and G. raimondii, respectively. Gene duplication analysis showed that whole-genome duplication (WGD) and tandem duplication were the major forces driving the generation of cotton ACO genes. In the promoters of GhACOs, there were cis-acting elements responding to stress, phytohormones, light, and circadian factors, indicating the possible involvement of GhACOs in these processes. Expression and co-expression analyses illustrated that most GhACOs were not only widely expressed in various tissues but also coexpressed with other genes in response to salt and drought stress. GhACO106_At overexpression in Arabidopsis promoted flowering and increased salt tolerance. These results provide a comprehensive overview of the ACO genes of cotton and lay the foundation for subsequent functional studies of these genes.

GhLUX1 and GhELF3 Are Two Components of the Circadian Clock That Regulate Flowering Time of Gossypium hirsutum

Photoperiod is an important external factor that regulates flowering time, the core mechanism of which lies in the circadian clock-controlled expression of FLOWERING LOCUS T (FT) and its upstream regulators. However, the roles of the circadian clock in regulating cotton flowering time are largely unknown. In this study, we cloned two circadian clock genes in cotton, GhLUX1 and GhELF3. The physicochemical and structural properties of their putative proteins could satisfy the prerequisites for the interaction between them, which was proved by yeast two-hybrid (Y2H) and Bimolecular Fluorescent Complimentary (BiFC) assays. Phylogenetic analysis of LUXs and ELF3s indicated that the origin of LUXs was earlier than that of ELF3s, but ELF3s were more divergent and might perform more diverse functions. GhLUX1, GhELF3, GhCOL1, and GhFT exhibited rhythmic expression and were differentially expressed in the early flowering and late-flowering cotton varieties under different photoperiod conditions. Both overexpression of GhLUX1 and overexpression of GhELF3 in Arabidopsis delayed flowering probably by changing the oscillation phases and amplitudes of the key genes in the photoperiodic flowering pathway. Both silencing of GhLUX1 and silencing of GhELF3 in cotton increased the expression of GhCOL1 and GhFT and resulted in early flowering. In summary, the circadian clock genes were involved in regulating cotton flowering time and could be the candidate targets for breeding early maturing cotton varieties.

Comparative transcriptome analysis reveals regulatory network and regulators associated with proanthocyanidin accumulation in persimmon

BACKGROUND

Proanthocyanidins (PAs) are important plant secondary metabolites that confer flavor, nutritional value, and resistance to pathogens. Persimmon is one of the PA richest crops. Mature fruits can be inedible because of the astringency caused by high PA levels and need to go through a de-astringency treatment before consumption. The molecular basis for PA accumulation is poorly known, particularly transcriptional regulators. We characterised three genotypes ('Luotiantianshi' (LT), 'Mopanshi' (MP), and 'Youhou' (YH)) with different PA accumulation patterns using an approach that combined PacBio full-length sequencing and Illumina-based RNA sequencing to build high-quality full-length transcriptomes. Additionally, we analysed transcriptome dynamics of the three genotypes (LT, MP, and YH) at four key fruit developmental stages.

RESULTS

A total of 96,463 transcripts were obtained. We identified 80,075 protein-coding sequences (CDSs), 71,137 simple sequence repeats (SSRs), and 27,845 long noncoding RNAs (lncRNAs). Pearson correlation coefficient (PCC), principal component analysis (PCA), and differentially expressed transcripts (DETs) analyses indicated that the four different developmental stages within a genotype exhibited similar transcriptome activities. A total of 2,164 transcripts specific to each fruit developmental stage were detected. The transcripts specific to early stages were attributed to phenylpropanoid and flavonoid biosynthesis. Co-expression network analyses revealed MEbrown and MEblue modules were strongly associated to PA accumulation. From these two modules, 20 hub TFs are potential regulators for PA accumulation. Among them, Cluster_78388 (SBP protein), Cluster_63454 (bZIP protein), and Cluster_66595 (MYB protein) appear to involve in the PA biosynthesis in Chinese genotypes.

CONCLUSIONS

This is the first high-quality reference transcriptome for commercial persimmon. Our work provides insights into the molecular pathways underlying PA accumulation and enhances our global understanding of transcriptome dynamics throughout fruit development.

Comprehensive identification and expression analysis of B-Box genes in cotton

BACKGROUND

B-BOX (BBX) proteins are zinc-finger transcription factors with one or two BBX domains and sometimes a CCT domain. These proteins play an essential role in regulating plant growth and development, as well as in resisting abiotic stress. So far, the BBX gene family has been widely studied in other crops. However, no one has systematically studied the BBX gene in cotton.

RESULTS

In the present study, 17, 18, 37 and 33 BBX genes were detected in Gossypium arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively, via genome-wide identification. Phylogenetic analysis showed that all BBX genes were divided into 5 main categories. The protein motifs and exon/intron structures showed that each group of BBX genes was highly conserved. Collinearity analysis revealed that the amplification of BBX gene family in Gossypium spp. was mainly through segmental replication. Nonsynonymous (Ka)/ synonymous (Ks) substitution ratios indicated that the BBX gene family had undergone purification selection throughout the long-term natural selection process. Moreover, transcriptomic data showed that some GhBBX genes were highly expressed in floral organs. The qRT-PCR results showed that there were significant differences in GhBBX genes in leaves and shoot apexes between early-maturing materials and late-maturing materials at most periods. Yeast two-hybrid results showed that GhBBX5/GhBBX23 and GhBBX8/GhBBX26 might interact with GhFT. Transcriptome data analysis and qRT-PCR verification showed that different GhBBX genes had different biological functions in abiotic stress and phytohormone response.

CONCLUSIONS

Our comprehensive analysis of BBX in G. hirsutum provided a basis for further study on the molecular role of GhBBXs in regulating flowering and cotton resistance to abiotic stress.

Phylogenetic Analysis of the Membrane Attack Complex/Perforin Domain-Containing Proteins in Gossypium and the Role of GhMACPF26 in Cotton Under Cold Stress

The membrane attack complex/perforin (MACPF) domain-containing proteins are involved in the various developmental processes and in responding to diverse abiotic stress. The function and regulatory network of the MACPF genes are rarely reported in Gossypium spp. We study the detailed identification and partial functional verification of the members of the MACPF family. Totally, 100 putative MACPF proteins containing complete MACPF domain were identified from the four cotton species. They were classified into three phylogenetic groups and underwent multifold pressure indicating that selection produced new functional differentiation. Cotton MACPF gene family members expanded mainly through the whole-genome duplication (WGD)/segmental followed by the dispersed. Expression and cis-acting elements analysis revealed that MACPFs play a role in resistance to abiotic stresses, and some selected GhMACPFs were able to respond to the PEG and cold stresses. Co-expression analysis showed that GhMACPFs might interact with valine-glutamine (VQ), WRKY, and Apetala 2 (AP2)/ethylene responsive factor (ERF) domain-containing genes under cold stress. In addition, silencing endogenous GhMACPF26 in cotton by the virus-induced gene silencing (VIGS) method indicated that GhMACPF26 negatively regulates cold tolerance. Our data provided a comprehensive phylogenetic evolutionary view of Gossypium MACPFs. The MACPFs may work together with multiple transcriptional factors and play roles in acclimation to abiotic stress, especially cold stress in cotton.

Evolutionary Relationships and Divergence of KNOTTED1-Like Family Genes Involved in Salt Tolerance and Development in Cotton (Gossypium hirsutum L.)

The KNOX (KNOTTED1-like homeobox) transcription factors play an important role in leaf, shoot apical meristem and seed development and respond to biotic and abiotic stresses. In this study, we analyzed the diversity and evolutionary history of the KNOX gene family in the genome of tetraploid cotton (Gossypium hirsutum). Forty-four putative KNOX genes were identified. All KNOX genes from seven higher plant species were classified into KNOXI, KNOXII, and KNATM clades based on a phylogenetic analysis. Chromosomal localization and collinearity analysis suggested that whole-genome duplication and a polyploidization event contributed to the expansion of the cotton KNOX gene family. Analyses of expression profiles revealed that the GhKNOX genes likely responded to diverse stresses and were involved in cotton growth developmental processes. Silencing of GhKNOX2 enhanced the salt tolerance of cotton seedlings, whereas silencing of GhKNOX10 and GhKNOX14 reduced seedling tolerance to salt stress. Silencing of GhSTM3 influenced the cotton flowering time and plant development. These findings clarify the evolution of the cotton KNOX gene family and provide a foundation for future functional studies of KNOX proteins in cotton growth and development and response to abiotic stresses.

Genome-wide identification and expression patterns analysis of the RPD3/HDA1 gene family in cotton

BACKGROUND

Histone deacetylases (HDACs) catalyze histone deacetylation and suppress gene transcription during various cellular processes. Within the superfamily of HDACs, RPD3/HDA1-type HDACs are the most studied, and it is reported that RPD3 genes play crucial roles in plant growth and physiological processes. However, there is a lack of systematic research on the RPD3/HDA1 gene family in cotton.

RESULTS

In this study, genome-wide analysis identified 9, 9, 18, and 18 RPD3 genes in Gossypium raimondii, G. arboreum, G. hirsutum, and G. barbadense, respectively. This gene family was divided into 4 subfamilies through phylogenetic analysis. The exon-intron structure and conserved motif analysis revealed high conservation in each branch of the cotton RPD3 genes. Collinearity analysis indicated that segmental duplication was the primary driving force during the expansion of the RPD3 gene family in cotton. There was at least one presumed cis-element related to plant hormones in the promoter regions of all GhRPD3 genes, especially MeJA- and ABA-responsive elements, which have more members than other hormone-relevant elements. The expression patterns showed that most GhRPD3 genes had relatively high expression levels in floral organs and performed higher expression in early-maturity cotton compared with late-maturity cotton during flower bud differentiation. In addition, the expression of GhRPD3 genes could be significantly induced by one or more abiotic stresses as well as exogenous application of MeJA or ABA.

CONCLUSIONS

Our findings reveal that GhRPD3 genes may be involved in flower bud differentiation and resistance to abiotic stresses, which provides a basis for further functional verification of GhRPD3 genes in cotton development and a foundation for breeding better early-maturity cotton cultivars in the future.