Next Article in Journal
Targeted-Deletion of a Tiny Sequence via Prime Editing to Restore SMN Expression
Previous Article in Journal
Identification and Characterization of Specific Nanobodies against Trop-2 for Tumor Targeting
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 14
10.3390/ijms23147940
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Xiao Li
1,2,
Yuanlong Wu
1,
Huabin Chi
1,
Hengling Wei
2,
Hantao Wang
2 and
Shuxun Yu
2,3,*
1
National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University, Wuhan 430000, China
2
State Key Laboratory of Cotton Biology, Institute of Cotton Research of CAAS, Anyang 455000, China
3
State Key Laboratory of Subtropical Silviculture, Zhejiang A & F University, Lin’an 311300, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(14), 7940; https://doi.org/10.3390/ijms23147940
Submission received: 20 June 2022 / Revised: 12 July 2022 / Accepted: 16 July 2022 / Published: 19 July 2022
(This article belongs to the Section Molecular Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
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.

1. Introduction

The transition from vegetative growth to reproductive development is a major physiological change in the life cycle of higher plants. Flowering regulatory networks in model plant species such as Arabidopsis thaliana [1,2,3,4,5,6,7] and rice [1,8,9,10,11] were thoroughly studied. According to environmental factors affecting flowering, flowering regulatory networks can be divided into five types: those involving the photoperiod pathway [1,7,12,13,14,15], hormone pathway [16,17,18,19,20,21], vernalization pathway [2,22,23], aging pathway, and autonomous pathway [24]. Furthermore, each flowering pathway has core flowering regulatory genes [25] that undergo cross-talk among each other [25,26,27,28,29].
In the photoperiod pathway, day length is perceived by the leaves and induces a systemic signal called florigen that moves through the phloem sieve elements to the shoot apical meristem (SAM); the FLOWERING LOCUST (FT) protein is a major component of florigen [4,30] in different species, such as Arabidopsis, rice, and temperate cereals [8]. At the SAM of various plant species, such as Arabidopsis, rice, wheat, temperate crops, and cereals [1,2,4,30], the FT-FLOWERING LOCUS D (FD) complex is involved in the network associated with the SAM, and causes changes in the expression of floral meristem-identity genes, which reprogram the stem cell differentiation fate and change the fate of the axillary primordia to form floral primordia instead of leaf primordia [4,5,31]. SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1), FRUITFULL (FUL), AGAMOUS-LIKE 24 (AGL24), LEAFY (LFY) and APETALA1 (AP1) involved in the initiation of flower primordia in the SAM are major meristem-identity genes in different species, such as Arabidopsis, rice, wheat, and various other cereals [1,2,4,30]. SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 3 (SPL3) is a direct upstream activator of LFY and AP1 [32], and it was investigated if OsSPL14 can be regulated by OsmiR156 to define the ideal plant architecture in rice via facilitating branch primordia differentiation [33]. Moreover, in the leaves, there are varieties of photoreceptors that sense light signals; the blue light photoreceptors cryptochrome (CRYs) [34], the other blue light photoreceptors ZEITLUPE/FLAVIN-BINDING, KELCH REPEAT, and F-BOX/LOV KELCH PROTEIN 2 (ZTL/FKF1/LKP2) [15,35,36,37], and the red/far-red light photoreceptor phytochromes (PHYA to PHYE) [13,38,39,40] are three major photoreceptors that regulate flowering time in Arabidopsis [41,42]. Recently, significant progress has been achieved on light receptors functioning on flowering in species such as rice [43], wheat [1,2], soybean [44], longan [45], Brassica juncea [46] and tomato [47]. Additionally, under long-day conditions, GIGANTEA (GI) can also directly facilitate FT transcription in the leaves [48,49].
The hormone pathway refers to the requirement of the pathway of hormones, especially gibberellic acid (GA) [25], auxin (IAA), and brassinosteroids (BRs) [21], for normal flowering patterns. In flowering, BR signals strongly undergo cross-talk among the GA, IAA, and light signaling pathways [21]. DELLA proteins, which are the transcription inhibitors of GA signals, inhibit both BRASSINAZOLE RESISTANT 1 (BZR1) and PHYTOCHROME-INTERACTING FACTORs (PIFs) to integrate the light, GA, and BR pathways [46]. Moreover, DELLA proteins are highly conserved among different species, including Arabidopsis, rice, maize, grape, wheat, barley, and Brassica [50,51]. BRs can repress the expression of PHYB/D, and activate the expression of CONSTITUTIVELY PHOTOMORPHOGENIC1 (COP1) and SUPPRESSOR OF PHYTOCHROME A 1 (SPA1) to regulate photomorphogenesis in Arabidopsis [21]. However, although BRs were first discovered in 1979 in extracts of Brassica napus pollen [51], the basic signaling pathway appears to be highly conserved in higher plants, the current BR network is mainly based on Arabidopsis, and the BR-signaling network in different plant species remains unclear [21].
Major floral initiation repressor FLOWERING LOCUS C (FLC), which is the core component of the vernalization pathway, can perceive changes in environmental temperature, regulate the expression of FT in leaves both in Arabidopsis and other species such as Aquilegia and Brassica napus [52,53,54,55], and respond to FT signals by directly repressing the expression of SOC1 and preventing the upregulation of FD in the SAM [22,52]. The autonomous pathway primarily facilitates flowering by downregulating FLC through RNA processing and epigenetic regulation. The crucial components involved in this pathway include FLOWERING LOCUS CA (FCA), FLOWERING LOCUS KH DOMAIN (FLK), FLOWERING LOCUS PA (FPA), FLOWERING LOCUS D (FLD), FLOWERING LOCUS VE (FVE), LUMINIDEPENDENS (LD), and FLOWERING LOCUS Y (FY) [24,25,54,56]. Most of the abovementioned genes are major components of autonomous pathways in various plant species, such as rice [57], canola [58], Doritaenopsis [59], and soybean [60].
These flowering pathways are not independent of each other; there is a complex link between them [25]. Flower repressor TEMPRANILLO (TEM) plays a pivotal role in the direct repression of FT by directly repressing the expression of the GA biosynthesis-related genes GA3-oxidase 1 and 2 (GA3OX1 and GA3OX2) to link the photoperiod and hormone pathways for flowering in A. thaliana [28,61]. Additionally, in soybean, BR can reduce the transcript level of GmRAV, the family to which TEMs belong [62]. Moreover, the balance between CONSTANS (CO) and TEM proteins hastens FT expression to trigger flowering [26]. GA biosynthesis-inhibiting gene SHORT VEGETATIVE PHASE (SVP) is involved in the control of flowering in Arabidopsis [29,63]. It was also investigated in other plant specises such as wheat [64], and temperate fruit trees including apple [65], Prunus mume (mei) [66], sweet cherries [67], and kiwifruit [68].
Cotton plants, which are members of the Malvaceae family and Gossypium genus, are considered to be a model species for the study of plant polyploidy [69]. Since most of the planted varieties are allotetraploids, the cotton genome is more complex, coupled with complex traits and a long growth period, and studies on flowering are key to finding a way to regulate flowering time to breed early-maturity cotton. Cotton flowering is positively correlated with the differentiation of flower primordia. However, cotton flower primordia differentiation is driven by axillary bud primordia, which is governed by SAM differentiation. Axillary bud primordia differentiate into leaf primordia during early growth and development; however, they differentiate into flower primordia during later growth and differentiation, marking the beginning of the cotton transition from vegetative growth to reproductive development.
Previous genomewide association study (GWAS) results showed that flowering time, the length of the period from the first flower bloom to the first boll opening, the length of the whole growth period, the yield percentage before frost, and the plant height and height of the first branch node were early-maturity traits, and some important flowering-related alleles have been identified, such as GhSPY, GhZTL, GhELF6, GhSVP, GhELF4, GhGA2OX6, and GhPHYA [69]. Furthermore, the GhCAL protein can form heterodimers with GhAP1 and GhAGL6 to regulate these two genes’ expression and facilitate cotton flowering [70]. Coincidentally, these flowering-related genes are highly homologous to those in A. thaliana, suggesting that the key flowering mechanism was relatively conserved during the evolution and domestication of plant species. However, how strong selection pressure alters the genome of cotton, particularly the genetic components underlying the adaptation to environmental changes and improved early maturity, remains unknown.
In comparison with extensive studies of the flowering mechanisms in model plant species A. thaliana, the flowering mechanisms in cotton remain unknown. In this study, to better understand the flowering mechanism of cotton, flowering-related candidate genes were identified with genomewide identification-based techniques. A detailed analysis of the flowering-related genes, including their functional classification, chromosomal locations, phylogenetic relationships, conserved motifs, and expression patterns, was performed. Moreover, the nonsynonymous substitution rate (Ka), synonymous substitution rate (Ks), and ω (Ka/Ks) ratio were calculated to determine the divergence time and selection pressure of homologous gene pairs. This systematic analysis of flowering-related genes in cotton provides a basis for further studies on the flowering mechanism of cotton.

2. Results

2.1. Comparison of Flowering-Related Genes from Arabidopsis thaliana, Gossypium hirsutum, Gossypium barbadense, and Gossypium arboretum

With the development and improvement of the A. thaliana flowering-related genomic network, 306 flowering time-related genes have been found, and information about their evolution and an interactive database, FLOR-ID, can be found at http://www.flor-id.org (accessed on 3 November 2021) [71]. On the basis of the protein sequences of the 306 flowering-related genes in A. thaliana, we investigated the flowering-related homologous genes in G. hirsutum and G. barbadense. In total, 636 and 673 flowering-related genes were identified in G. hirsutum (Table S1) and G. barbadense (Table S2), respectively, and more than half of the A. thaliana flowering-related genes had putative G. hirsutum (176 out of 306 genes) and G. barbadense (174 out of 306 genes) homologs. As in the G. barbadense genome, homologs in the G. hirsutum genome were lacking for 130 genes, and most of these genes (120 of 130 genes) have a single functional effect on flowering (Tables S3 and S4). Furthermore, 350 flowering-related genes in G. arboreum were identified (Table S5), and more than half of the flowering-related genes in A. thaliana had putative homologs in G. arboretum (173 out of 306 genes; Table S5); in addition, 133 genes have been lost (Table S6).
To explore the flowering mechanism among the cultivated cotton species (G. hirsutum, G. arboretum, and G. barbadense), the putative flowering-related genes were classified into the eight following gene sets according to the classification of A. thaliana genes: those involved in the photoperiod pathway, circadian clock, and light signaling (Ph); sugar (Su); vernalization (Ve); the autonomous pathway (Au); hormone signaling and metabolism (Ho); the aging pathway (Ag); flower development and the apical meristem response pathway (Fd); and the ambient temperature pathway (At) [71]. Interestingly, most of the lost (53) genes belonged to the photoperiod pathway in G. hirsutum. There were 320 genes involved in the photoperiod pathway, and there were too many overlapping genes between the photoperiod and vernalization (62 genes), and hormone (39) pathways. The fewest common genes were detected between the photoperiod and ambient temperature pathways (three genes) (Figure 1A), and similar results were found in G. barbadense (Figure 1B).
The number of flowering-related genes in G. hirsutum (636 genes) and G. barbadense (673 genes) was approximately twice that in G. arboretum (350 genes, Figure 1 and Figure S1). Moreover, there were no significant differences in the number of flowering-related genes between G. hirsutum (636 genes) and G. barbadense (673 genes). By classifying the putative flowering-related genes in G. hirsutum, G. barbadense, and G. arboretum into the eight gene sets according to the A. thaliana dataset, we found that most of the genes participate in the photoperiod pathway (Figure S1). Additionally, there was no significant difference in the number of G. hirsutum and G. barbadense genes associated with each gene set (Figure S1).
Additionally, the hormone signaling gene sets (16.07% were lost) or meristem response and development gene sets (20% were lost) were preferentially retained instead of the other pathway genes. However, more than 30% of the light signaling (38.78%), vernalization (30.86%), and aging (47.5%) genes were lost (Table S3).

2.2. Chromosomal Localization and Duplication Analysis of Flowering-Related Genes in G. hirsutum and G. barbadense

The 636 and 673 flowering-related genes in G. hirsutum (Figure 2 and Table S7) and G. barbadense (Figure S2 and Table S8), respectively, were mapped onto pseudomolecular chromosomes. These genes exhibited equal genomic distribution, with 324 (G. hirsutum) and 329 (G. barbadense) genes localized on the At subgenome, and 310 (G. hirsutum) and 329 (G. barbadense) genes localized on the Dt subgenome. However, the chromosomal distribution of these genes was uneven. For the At subgenome, 37 (G. hirsutum) and 34 (G. barbadense) genes were localized on chromosome A05, with most located on the top half, representing 11.42% (G. hirsutum) and 10.33% (G. barbadense) of the flowering-related genes in the At subgenome (Figure 2A and Figure S2A). Only 12 flowering-related genes were detected on chromosomes A04 (3.70%, G. hirsutum) and A02 (3.65%, G. barbadense), with most distributed at both ends (Figure 2A and Figure S2A). In addition, 34 (10.97%, G. hirsutum) and 33 (10.03%, G. barbadense) genes were located on chromosome D05, with most located on the top half. Only 11 genes (3.55%) were detected on chromosome D03 of G. hirsutum, and 12 genes (3.65%) were detected on D01 and D04 of G. barbadense, with most located on the bottom half (Figure 2B and Figure S2B). Interestingly, most of the other genes were distributed near the ends of chromosomes (Figure 2 and Figure S2).

2.3. Floral Pathway Integrators

In flowering, endogenous elements always undergo cross-talk with exogenous factors through certain integrators [4]. The functional classification of flowering-related genes showed that systemic signal FT [1,4,8,28], which is called GhFT or TWIN SISTER OF FT (GhTSF) in G. hirsutum, functions in all eight flowering pathways, namely, Ph, Ho, Ag, At, Ve, Fd, Au, and Su [1,2,4,7,8] (Figure 3). In addition, floral pathway integrator SOC1 [72,73], and floral meristem-identity genes AP1 [2,74,75], LFY [76,77,78] and SVP [79] (Figures S3 and S4) are found in G. hirsutum. GhSOC1, also named GhAGL20 and GhLFY, participates in the Ph, Ho, Ag, At, Ve, Fd and Su pathways. GhAP1, also called GhAGL8, is regulated by the Ph, Fd, Ag, At, Su and Ve flowering pathways, and GhSVP, also called GhAGL24, works in the At, Ag, Ph, Fd and Ve flowering pathways in G. hirsutum (Figure S3).
Furthermore, phylogenetic analysis revealed that floral pathway integrators GhFT, GhSOC1, GhSVP and GhAP1 have a close genetic relationship (Figure 3A), which is consistent with the findings of our exon–intron structural analysis; 3 identified GhFT homologs contained 4 exons, 4 GhSVP homologs contained 7 exons, and 9 GhAP1 homologs contained 7 to 8 exons (Figure 3B).
To systematically determine the conserved motifs of the flowering-related genes involved in the Ph pathway, the distribution of conserved motifs of these genes was estimated by using the online MEME server in conjunction with 25 putative conserved motifs (Figure S5). Results show that the length of the conserved motifs ranged from 19 to 50 amino acids, and the number of conserved motifs within each flowering-related gene involved in the Ph pathway ranged from 0 to 10 (Figure 3C,D), showing that all GhSVP, GhAP1 and GhSOC1 homologs contained three conserved motifs. Interestingly, the three GhFT homologs contained no conserved motifs with other Ph pathway flowering-related genes (Figure 3C).

2.4. Flowering-Related Genes Involved in the Photoperiod Pathway, Circadian Clock, and Light Signalling

Light is indispensable for plant life, and the perception of the light environment dictates seed germination, photomorphogenesis, phototropism, shade avoidance, and flowering [61]. In this study, 320 genes in G. hirsutum homologous to 174 genes related to the photoperiod pathway, circadian clock, and light signaling in A. thaliana were identified (Figure 1 and Table S1). To explore the evolution and genomic structure of Ph pathway-related genes, the phylogenetic analysis and intron diagrams of Ph pathway flowering-related genes were generated on the basis of their sequences (Figure 3A,B). The findings show that GhSPL3, GhGI, and photoreceptors GhFKF1, GhCRY1/2, and GhPHYB are closely genetically related (Figure 3A). The exon number of each gene widely ranged from 1 to 15 exons. Significantly, 3 GhELF4 genes were identified, and all of them contained only 1 exon, 2 GhSPL3 copies had 10 exons, and 4 GhPHYB copies had 3 to 5 exons. Four GhGIs were identified and contained 14 to 15 exons (Figure 3B). FKF1 belongs to a family of F-box proteins of which ZTL and LKP2 are members, and all of them comprise light, oxygen, or voltage (LOV) domains [35]. Like in G. barbadense, there are eight copies of each of the genes GhZTL, GhLKP2, and GhFKF1 in G. hirsutum (Table S2). Moreover, seven GhFKF1 genes contained two exons (Figure 3B). Blue light photoreceptors CRY1 and CRY2, which monitor light signals to regulate plant flowering [34,80], were identified in G. hirsutum (Table S1), and four GhCRY1 and GhCRY2 contained four exons (Figure 3B).
Furthermore, most GA biosynthesis homologs, such as GA2ox, GA3ox, and GA20ox [81], were retained (Figure 3); they were simultaneously induced in response to light signaling and contained six conserved motifs in G. hirsutum (Figure 3C and Table S1). Our conserved motif analysis showed that genes GhELF4, GhLFY, GhSPL3, and GhPHYB had no conserved motifs. However, the five GhTEM2 genes had two conserved motifs, the GhGI genes contained one conserved motif, the GhFKF1 genes had eight conserved motifs, and the GhCRYs had seven conserved motifs (Figure 3C).

2.5. Expression Analysis of Ph Pathway Genes in G. hirsutum

To investigate the divergence in the expression levels of homologous genes, and their participation in the photoperiod and circadian clock pathways, we analyzed the expression patterns of these supposed flowering-related genes in G. hirsutum, including the root, stem, leaf, petal, anther, stigma, ovule 0D (Day), ovule 1D, ovule 3D, ovule 10D, ovule 20D, fiber 5D, fiber 10D, fiber 20D and fiber 25D tissues (Figure 3D and Table S9). The results showed that most of the expressed genes were in the roots, petals, stems and leaves (Figure 3D and Table S9). Four GhCRY1/2, four GhGA20ox1, two GhFKF1, two GhPHYB, three GhTEM2 and one GhSPL3 were specifically or preferentially expressed in ovule 1D tissue. Four GhSVP, four GhFKF1, four GhGI, two GhSOC1, two GhGA2ox, two GhELF4, one GhAP1, and one GhFT were specifically or preferentially expressed in the stems and leaves. Eight GhAP1, three GhCRY1, two GhSOC1, two GhPHYB, two GhSPL3, and two GhFT genes were specifically or preferentially expressed in the roots or leaves. Only one copy of GhLFY was found in G. hirsutum, and it was preferentially expressed in the roots (Figure 3D).

2.6. Hormone-Pathway-Related Flowering-Related Genes

In A. thaliana, by controlling the spatial expression of floral regulatory genes throughout the plant independent of light signals, GA signals promote the initiation of floral primordia [27]. 85 genes likely related to the hormone pathway were identified in G. hirsutum genome. GA is perceived by its receptor, GID1, and there are 12 GhGID1s in G. hirsutum (Figure 4A–C and Table S1). CID1s can bind to bioactive GA through the conformational changes, which then facilitates the interaction between CID1 and DELLA proteins, which are major flowering repressors [82,83,84]. There are five DELLA family genes, namely, GIBBERELLIC ACID INSENSITIVE (GAI), REPRESSOR OF ga13 (RGA), RGA-LIKE1 (RGL1), RGL2 and RGL3, in the A. thaliana genome [85,86]. We found four GhGAIs, namely, GhRGA, GhRGL1, GhRGL2 and GhRGL3 in G. hirsutum (Figure 4A–C and Table S1).
To further study the structural features of these genes, the exon–intron structure was analyzed, the results of which show that the 12 GhGID1s possessed 2 exons, and the 4 GhGAI only had 1 exon (Figure 4A). To visualize the conserved residues, the Ho pathway flowering-related gene motifs were analyzed with 10 putative conserved motifs (Figure S6), the results of which show that the GhCID1 genes have eight conserved motifs, two GhGAI genes have four conserved motifs, and one GhRGA and one GhRGL have three conserved motifs (Figure 4B). Phylogenetic analysis revealed that the GhGID1s and GhGAI genes are closely related to their homologs (Figure 4C). The expression pattern analysis of 15 different cotton tissues showed that GhGAI has divergent expression patterns among their different copies coupled with sequence similarity, five GhGAI genes were highly expressed in petals, and six other genes were expressed in the leaf, stem, and ovule 0D tissues. Two GhGAI genes were highly expressed in ovule 0D; in contrast, two other genes showed nearly no expression in ovule 0D (Figure 4D and Table S9).
We detected content changes in hormones zeatin riboside (ZR), abscisic acid (ABA), GA, and IAA in the SAM and leaves of two early-maturing cultivars Zhong50 and Zhong74; two late-maturing cultivars Zhong60 and Lu28 at different true-leaf flattening stages (first true-leaf stage 1TLS, second true-leaf stage 2TLS, third true-leaf stage 3TLS, and fourth true-leaf stage 4TLS). Results show that, during the development of the SAMs in cotton, the contents of hormones ZR, ABA, and GA are maximized; conversely, hormone IAA reached its minimal value during the second to third true-leaf stages (Figure 4E–H). Moreover, the balance between various endogenous hormones peaked at 2 TLS in cotton cultivars Zhong74, and at 3 TLS in Zhong50, Zhong60 and Lu28 (Figure S7).

2.7. Selection Pressure on Flowering Pathway Gene Sets

To investigate the evolution process, the nonsynonymous/synonymous substitution ratio (Ka/Ks), which is related to the evolutionary selection patterns of the corresponding genome [87], of flowering-related gene orthologous pairs between the G. hirsutum and A. thaliana genomes as well as between the G. barbadense and A. thaliana genomes were computed to determine the selection pressure. The calculation of the selection pressure showed that the mean Ka/Ks ratios between A. thaliana and G. hirsutum of different flowering-related gene sets ranged from 0.001 to 0.200 (Figure 5, Tables S10 and S11). In particular, genes involved in flower development and the apical meristem response pathway appear to have been subjected to less negative selection pressure than the other pathway genes were. Moreover, the range of the selection pressure variation was largest in Ph pathway genes during the evolution of cotton (Figure 5A,B, Tables S10 and S11).

3. Discussion

Compared with the A. thaliana flowering mechanism, which has been comprehensively described, the exploration of cotton flowering mechanisms is still deficient. Cotton species G. hirsutum, G. barbadense, and G. arboretum were sequenced, and the release of genomic data makes it greatly convenient for studying the genes and functions of important agronomic traits [88]. Here, a bioinformatics approach was used to analyze the G. hirsutum and G. barbadense genes that potentially participate in flowering, and implied a possible regulatory model of the multiple feedbacks and inputs during the regulation of cotton flowering in the leaf and meristem (Figure 6).
In this study, the number of flowering-related genes identified in G. hirsutum was nearly equal to that in G. barbadense and twice that in G. arboretum, which is consistent with the results that G. hirsutum and G. barbadense, both of which are allotetraploid cotton species, were formed by the two closest diploids extant progenitor hybridization, evolution and artificial domestication [88,89,90]. Interestingly, many of the detected flowering-related genes were functionally redundant, suggesting that these genes may play an integrated role during evolution and domestication; moreover, the photoperiod pathway may undergo cross-talk with the vernalization and hormone pathways and play a substantial role in regulating cotton flowering. Inversely, most of the flowering-related homologous that were found in A. thaliana but not in G. hirsutum had a single function, implying that these genes may not participate in the core flowering pathways. Moreover, the genetic localization results show that the flowering mechanism of cotton might be relatively conserved among the allotetraploid species and between the subgenomes. Meanwhile, more than 30% of the light signaling (38.78%), vernalization (30.86%) and aging (47.5%) gene sets were lost. Considering that currently grown cotton varieties evolved and were domesticated from perennial shrub plants in the subtropics [69,89,90], the loss of genes in the photoperiod pathway, vernalization, and aging gene sets may be related to changes in photoperiod and environmental temperature during their northward movement. Notably, all of the floral pathway integrator SOC1 [72,73] and the floral meristem-identity genes AP1 [2,74], LFY [75,76,77,78] and SVP [79] are activated to regulate floral primordium initiation, facilitate the transition from vegetative growth to reproductive growth, and accelerate early flowering in SAM during the transition from vegetative growth to reproductive development in Arabidopsis [4,73], and all of them were identified in G. hirsutum. In addition, a previous study had investigated that GhCAL can form heterodimers with GhAP1/GhAGL6 to regulate their expression and facilitate cotton flowering [70]. GWAS results [69] indicate that, despite the loss of some of the flowering-related genes, the core mode of flowering is evolutionarily conserved in G. hirsutum and A. thaliana to a certain extent, and the basic light signaling-related flowering pathways are likely relatively conserved. These results indicated that the differentiation fate transition from a vegetative SAM to a reproductive SAM might be a pivotal cell biological characteristic affecting the early flowering of cotton. The flowering mechanisms have been relatively conserved during evolution and domestication, and all of these core factors may also play an indispensable role during floral primordium initiation in cotton SAM, which seems to explain why G. hirsutum can flower in temperate areas. Interestingly, combined with the core role of FT in the flowering of Arabidopsis, rice, and temperate crops [1,2,8], the analysis of gene conserved motifs in G. hirsutum, showed that GhFT have no conserved motifs with other flowering-related genes that identified in this study, which indicated that GhFT homologs in cotton likely also play a core part different with other genes during the evolution, and the function of which might can’t be reproduced by other genes.
Additionally, plants perceive light signaling through their leaves [7] through several photoreceptors, such as FKF1/ZTL/LKP2 [35,36,37,91,92], CRY1/2 [34,42] and PHYA/B/C/D/E [38,39,40], and induce florigen, a systemic signal that consists of primarily FT protein and moves through the phloem to SAM [5,31,93]. At the SAM, via reprograming the expression of floral identity genes, such as AP1, SOC1, LFY and SVP, FT drives SAM differentiation such that floral primordia are produced [2,4]. In this study, the homologs of GhCRY1/2, GhFKF1/ZTL/LKP2 and GhPHYA/B/C/D/E were identified in the G. hirsutum genome. Consistently, the identification of GA biosynthesis- and metabolism-related genes GID1, GA20ox1/2, GA30x1, and GA2ox1/2/3/6/8 [82,83,84] as well as the DELLA family members GAI, RGA, RGL1, RGL2 and RGL3 [85,86], coupled with the changes in hormone contents in the SAM and leaves during the transition of floral primordia in cotton. Furthermore, GA signaling can not only function in leaves by promoting the expression of flowering time integrator genes such as FT and TSF independently of CO and GI [27], but also regulates the expression of SPL genes in both the leaves and the SAM [27]. All of the findings suggest that hormones, especially GA signals and light signals may participate in the initiation of the floral primordia during the transition from vegetative growth to the reproductive stage.
Furthermore, expression pattern analysis indicates that, in cotton, during evolution and domestication, duplicated genes can undergo nonfunctionalization, neofunctionalization, or subfunctionalization. Genes preferentially expressed in stems and leaves may perceive a light signal in the leaves, transport it through the phloem to the SAM, and stimulate the initiation of floral primordia [4,5,31,93]. Remarkably, selective pressure analysis shows that all Ka/Ks ratio values were lower than 0.200, which indicates that the flowering-related genes for cotton underwent strong purifying selection and that the negative selection had acted against extreme polymorphic variants of flowering-related genes in cotton. Additionally, the sequences of flower development and apical meristem response pathway genes exhibit great sequence diversity, with higher Ka/Ks ratios than those of other flowering pathways between not only the G. hirsutum and A. thaliana genomes but also the G. barbadense and A. thaliana genomes, indicating that the cotton flowering-related genes underwent strong purification during evolution, and the same was found between A. thaliana and G. barbadense, and variation in flower development and apical meristem response pathway genes may contribute to the rapid evolutionary capacity to change the thermal requirement to flowering in G. hirsutum.

4. Materials and Methods

4.1. Data Resources

Details regarding the annotated A. thaliana genome were downloaded from the TAIR10 website (http://www.arabidopsis.org/index.jsp, accessed on 10 December 2021). G. hirsutum, G. barbadense, and G. arboretum genomic data and the fragments per kilobase of transcript per million mapped reads (FPKM) of the flowering-related genes for G. hirsutum were downloaded from the CottonFGD database (https://cottonfgd.org/about/download.html, accessed on 13 January 2022). The genomic, coding, and protein sequences, and the gene expression of flowering-related genes were downloaded from the CottonFGD database (accessed on 25 January 2022).

4.2. Identification of Flowering-Related Gene Homologs in Cotton

Homologous genes were detected by using a combination of similarity- and synteny-based approaches. In the similarity-based approach, the flowering-related genes in A. thaliana were used as queries, and the BLASTP program was used to perform searches against G. hirsutum, G. barbadense, and G. arboretum protein sequences with the following conditions: E-value < 1 × 10−20, identity > 60%, coverage > 75%, and match length > 70 amino acids.

4.3. Chromosomal Localization of Flowering-Related Genes in the Cotton Genome

The locations on the chromosomes or scaffolds indicate the distribution of flowering-related genes, and detailed information on genome localization for the predicted G. hirsutum, G. barbadense and G. arboretum flowering-related genes was obtained from the cotton genome data. A map was constructed using the software of MapChart (MapChart 2.32, which is supported by dr.ir. RE (Roeland) Voorrips in Wageningen, The Netherlands; https://www.wur.nl/en/show/Mapchart.htm, accessed on 5 February 2022) to visualize the putative flowering-related genes on pseudomolecular chromosomes.

4.4. Analysis of Multiple Alignments and Phylogenetic, Gene Structure and Motif Recognition for Flowering-Related Genes

Multiple gene protein sequences were aligned by using ClustalW 1.81 software, which is currently maintained at the Conway Institute UCD Dublin by Des Higgins, Fabian Sievers, David Dineen, and Andreas Wilm (Dublin, Germany). A phylogenetic tree was constructed using the MEGA 6.0 program with the maximum likelihood (ML) method (1000 bootstrap replicates) on the basis of full-length protein sequences.
To visualize the structure of flowering-related genes and illustrate the organization of exons and introns in G. hirsutum, Gene Structure Display Server 2.0 (GSDS 2.0, http://gads.cbi.pku.edu.cn/index.php, accessed on 16 February 2022) was used.
Multiple Em for Motif Elicitation (MEME Version 4.11.4, http://meme-suite.org/tools/meme, accessed on 17 February 2022) was used to identify the conserved motifs of the protein sequences of significant flowering-related genes in G. hirsutum with default parameters except for the maximum number of motifs, which was set to 25.

4.5. Ka/Ks Ratios of Flowering-Related Gene Pairs between A. thaliana and G. hirsutum/G. barbadense

The Ka/ Ks ratios of homologous gene pairs were calculated to investigate the molecular evolution of homologous gene pairs. In the calculations of Ka/Ks, the full-length amino acid sequences of the G. hirsutum/G. barbadense and A. thaliana flowering-related genes underwent pairwise alignments using MUSCLE. Then, the aligned amino acid sequences were translated into the corresponding nucleotide coding sequences using Perl scripts derived from ParaAT_2.0 software. Lastly, the translated nucleotide coding sequences were used as input files for computing Ka/Ks values by using the Phylogenetic Analysis by Maximum Likelihood (PAML) program (PAML4.9j, http://www.bork.embl.de/pal2nal/index.cgi, accessed on 25 February 2022) with the default parameters. All variable sites of the alignment pairs were used in the calculation of Ka/Ks.

5. Conclusions

These findings provide genomic evidence and genomewide information for the investigation of the mechanism of early cotton maturation, and the flowering mechanism of cotton might be relatively conserved among the allotetraploid species and between the subgenomes. The core mode of flowering is also evolutionarily conserved between G. hirsutum and A. thaliana, and the basic light and hormone-signaling-related flowering pathways are likely relatively conserved. Additionally, the cotton flowering-related genes underwent strong purification during evolution. On the basis of these results, we provide a possible regulatory model of the multiple feedbacks and inputs during the regulation of cotton flowering in the leaves and SAM, which are further systematic and comprehensive flowering regulatory molecular networks of allotetraploid cotton species.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23147940/s1.

Author Contributions

Conceptualization and funding acquisition, S.Y.; Investigation, X.L.; Data curation and formal analysis, X.L. and Y.W.; Formal analysis, H.C.; Funding acquisition and resources, S.Y., H.W. (Hengling Wei) and H.W. (Hantao Wang); Writing—original draft preparation, X.L.; Writing—review and editing, S.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This project was financially supported by the National Key R&D Program of China (2020YFD1001004) and the China Agriculture Research System (grant no. CARS-15-06).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would particularly like to acknowledge the cotton team of the State Key Laboratory of Huazhong Agricultural University for their support and help in this study.

Conflicts of Interest

The authors declare that there are no conflict of interest.

References

  1. Song, Y.H.; Shim, J.S.; Kinmonth-Schultz, H.A.; Imaizumi, T. Photoperiodic flowering: Time measurement mechanisms in leaves. Annu. Rev. Plant Biol. 2015, 66, 441–464. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Kim, D.H.; Doyle, M.R.; Sung, S.; Amasino, R.M. Vernalization: Winter and the timing of flowering in plants. Annu. Rev. Cell Dev. Biol. 2009, 25, 277–299. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Fletcher, J.C. Shoot and floral meristem maintenance in arabidopsis. Annu. Rev. Plant Biol. 2002, 53, 45–66. [Google Scholar] [CrossRef]
  4. Turck, F.; Fornara, F.; Coupland, G. Regulation and identity of florigen: FLOWERING LOCUS T moves center stage. Annu. Rev. Plant Biol. 2008, 59, 573–594. [Google Scholar] [CrossRef] [Green Version]
  5. Turnbull, C. Long-distance regulation of flowering time. J. Exp. Bot. 2011, 62, 4399–4413. [Google Scholar] [CrossRef] [Green Version]
  6. Demarsy, E.; Fankhauser, C. Higher plants use LOV to perceive blue light. Curr. Opin. Plant Biol. 2009, 12, 69–74. [Google Scholar] [CrossRef] [Green Version]
  7. Song, Y.H.; Ito, S.; Imaizumi, T. Flowering time regulation: Photoperiod- and temperature-sensing in leaves. Trends Plant Sci. 2013, 18, 575–583. [Google Scholar] [CrossRef] [Green Version]
  8. Shrestha, R.; Gomez-Ariza, J.; Brambilla, V.; Fornara, F. Molecular control of seasonal flowering in rice, arabidopsis and temperate cereals. Ann. Bot. 2014, 114, 1445–1458. [Google Scholar] [CrossRef] [Green Version]
  9. Itoh, H.; Nonoue, Y.; Yano, M.; Izawa, T. A pair of floral regulators sets critical day length for Hd3a florigen expression in rice. Nat. Genet. 2010, 42, 635–638. [Google Scholar] [CrossRef]
  10. Zhang, D.; Yuan, Z. Molecular control of grass inflorescence development. Annu. Rev. Plant Biol. 2014, 65, 553–578. [Google Scholar] [CrossRef]
  11. Yang, R.; Li, P.; Mei, H.; Wang, D.; Sun, J.; Yang, C.; Hao, L.; Cao, S.; Chu, C.; Hu, S.; et al. Fine-Tuning of MiR528 Accumulation Modulates Flowering Time in Rice. Mol. Plant 2019, 12, 1103–1113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Giakountis, A.; Coupland, G. Phloem transport of flowering signals. Curr. Opin. Plant Biol. 2008, 11, 687–694. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Valverde, F.; Mouradov, A.; Soppe, W.; Ravenscroft, D.; Samach, A.; Coupland, G. Photoreceptor regulation of CONSTANS protein in photoperiodic flowering. Science 2004, 303, 1003–1006. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Lee, B.D.; Cha, J.Y.; Kim, M.R.; Paek, N.C.; Kim, W.Y. Photoperiod sensing system for timing of flowering in plants. BMB Rep. 2018, 51, 163–164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Christie, J.M.; Salomon, M.; Nozue, K.; Wada, M.; Briggs, W.R. LOV (light, oxygen, or voltage) domains of the blue-light photoreceptor phototropin (nph1): Binding sites for the chromophore flavin mononucleotide. Proc. Natl. Acad. Sci. USA 1999, 96, 8779–8783. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Yamaguchi, N.; Wu, M.F.; Winter, C.M.; Berns, M.C.; Nole-Wilson, S.; Yamaguchi, A.; Coupland, G.; Krizek, B.A.; Wagner, D. A molecular framework for auxin-mediated initiation of flower primordia. Dev. Cell 2013, 24, 271–282. [Google Scholar] [CrossRef] [Green Version]
  17. Yamaguchi, N.; Winter, C.M.; Wu, M.F.; Kanno, Y.; Yamaguchi, A.; Seo, M.; Wagner, D. Gibberellin acts positively then negatively to control onset of flower formation in Arabidopsis. Science 2014, 344, 638–641. [Google Scholar] [CrossRef]
  18. Jokela, V.; Virkajarvi, P.; Tanskanen, J.; Seppanen, M.M. Vernalization, gibberellic acid and photo period are important signals of yield formation in timothy (Phleum pratense). Physiol. Plant. 2014, 152, 152–163. [Google Scholar] [CrossRef]
  19. Campos-Rivero, G.; Osorio-Montalvo, P.; Sanchez-Borges, R.; Us-Camas, R.; Duarte-Ake, F.; De-la-Pena, C. Plant hormone signaling in flowering: An epigenetic point of view. J. Plant Physiol. 2017, 214, 16–27. [Google Scholar] [CrossRef]
  20. Yu, S.; Galvao, V.C.; Zhang, Y.C.; Horrer, D.; Zhang, T.Q.; Hao, Y.H.; Feng, Y.Q.; Wang, S.; Schmid, M.; Wang, J.W. Gibberellin regulates the Arabidopsis floral transition through miR156-targeted SQUAMOSA promoter binding-like transcription factors. Plant Cell 2012, 24, 3320–3332. [Google Scholar] [CrossRef] [Green Version]
  21. Wang, Z.Y.; Bai, M.Y.; Oh, E.; Zhu, J.Y. Brassinosteroid signaling network and regulation of photomorphogenesis. Annu. Rev. Genet. 2012, 46, 701–724. [Google Scholar] [CrossRef] [PubMed]
  22. Searle, I.; He, Y.; Turck, F.; Vincent, C.; Fornara, F.; Krober, S.; Amasino, R.A.; Coupland, G. The transcription factor FLC confers a flowering response to vernalization by repressing meristem competence and systemic signaling in Arabidopsis. Genes Dev. 2006, 20, 898–912. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Sheldon, C.C.; Finnegan, E.J.; Peacock, W.J.; Dennis, E.S. Mechanisms of gene repression by vernalization in Arabidopsis. Plant J. 2009, 59, 488–498. [Google Scholar] [CrossRef]
  24. Simpson, G.G. The autonomous pathway: Epigenetic and post-transcriptional gene regulation in the control of Arabidopsis flowering time. Curr. Opin. Plant Biol. 2004, 7, 570–574. [Google Scholar] [CrossRef] [PubMed]
  25. Srikanth, A.; Schmid, M. Regulation of flowering time: All roads lead to Rome. Cell Mol. Life Sci. 2011, 68, 2013–2037. [Google Scholar] [CrossRef]
  26. Castillejo, C.; Pelaz, S. The balance between CONSTANS and TEMPRANILLO activities determines FT expression to trigger flowering. Curr. Biol. 2008, 18, 1338–1343. [Google Scholar] [CrossRef] [Green Version]
  27. Galvao, V.C.; Horrer, D.; Kuttner, F.; Schmid, M. Spatial control of flowering by DELLA proteins in Arabidopsis thaliana. Development 2012, 139, 4072–4082. [Google Scholar] [CrossRef] [Green Version]
  28. Osnato, M.; Castillejo, C.; Matias-Hernandez, L.; Pelaz, S. TEMPRANILLO genes link photoperiod and gibberellin pathways to control flowering in Arabidopsis. Nat. Commun. 2012, 3, 808. [Google Scholar] [CrossRef] [Green Version]
  29. Andres, F.; Porri, A.; Torti, S.; Mateos, J.; Romera-Branchat, M.; Garcia-Martinez, J.L.; Fornara, F.; Gregis, V.; Kater, M.M.; Coupland, G. SHORT VEGETATIVE PHASE reduces gibberellin biosynthesis at the Arabidopsis shoot apex to regulate the floral transition. Proc. Natl. Acad. Sci. USA 2014, 111, E2760–E2769. [Google Scholar] [CrossRef] [Green Version]
  30. Andres, F.; Coupland, G. The genetic basis of flowering responses to seasonal cues. Nat. Rev. Genet. 2012, 13, 627–639. [Google Scholar] [CrossRef]
  31. Notaguchi, M.; Abe, M.; Kimura, T.; Daimon, Y.; Kobayashi, T.; Yamaguchi, A.; Tomita, Y.; Dohi, K.; Mori, M.; Araki, T. Long-distance, graft-transmissible action of Arabidopsis FLOWERING LOCUS T protein to promote flowering. Plant Cell Physiol. 2008, 49, 1645–1658. [Google Scholar] [CrossRef] [Green Version]
  32. Yamaguchi, A.; Wu, M.F.; Yang, L.; Wu, G.; Poethig, R.S.; Wagner, D. The microRNA-regulated SBP-Box transcription factor SPL3 is a direct upstream activator of LEAFY, FRUITFULL, and APETALA1. Dev. Cell 2009, 17, 268–278. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Jiao, Y.; Wang, Y.; Xue, D.; Wang, J.; Yan, M.; Liu, G.; Dong, G.; Zeng, D.; Lu, Z.; Zhu, X.; et al. Regulation of OsSPL14 by OsmiR156 defines ideal plant architecture in rice. Nat. Genet. 2010, 42, 541–544. [Google Scholar] [CrossRef] [PubMed]
  34. Ahmad, M.; Cashmore, A.R. HY4 gene of A. thaliana encodes a protein with characteristics of a blue-light photoreceptor. Nature 1993, 366, 162–166. [Google Scholar] [CrossRef] [PubMed]
  35. Ito, S.; Song, Y.H.; Imaizumi, T. LOV domain-containing F-box proteins: Light-dependent protein degradation modules in Arabidopsis. Mol. Plant 2012, 5, 573–582. [Google Scholar] [CrossRef] [Green Version]
  36. Suetsugu, N.; Wada, M. Evolution of three LOV blue light receptor families in green plants and photosynthetic stramenopiles: Phototropin, ZTL/FKF1/LKP2 and aureochrome. Plant Cell Physiol. 2013, 54, 8–23. [Google Scholar] [CrossRef] [Green Version]
  37. Takase, T.; Nishiyama, Y.; Tanihigashi, H.; Ogura, Y.; Miyazaki, Y.; Yamada, Y.; Kiyosue, T. LOV KELCH PROTEIN2 and ZEITLUPE repress Arabidopsis photoperiodic flowering under non-inductive conditions, dependent on FLAVIN-BINDING KELCH REPEAT F-BOX1. Plant J. 2011, 67, 608–621. [Google Scholar] [CrossRef]
  38. Franklin, K.A.; Quail, P.H. Phytochrome functions in Arabidopsis development. J. Exp. Bot. 2010, 61, 11–24. [Google Scholar] [CrossRef] [Green Version]
  39. Franklin, K.A.; Praekelt, U.; Stoddart, W.M.; Billingham, O.E.; Halliday, K.J.; Whitelam, G.C. Phytochromes B, D, and E act redundantly to control multiple physiological responses in Arabidopsis. Plant Physiol. 2003, 131, 1340–1346. [Google Scholar] [CrossRef] [Green Version]
  40. Endo, M.; Tanigawa, Y.; Murakami, T.; Araki, T.; Nagatani, A. PHYTOCHROME-DEPENDENT LATE-FLOWERING accelerates flowering through physical interactions with phytochrome B and CONSTANS. Proc. Natl. Acad. Sci. USA 2013, 110, 18017–18022. [Google Scholar] [CrossRef] [Green Version]
  41. Moglich, A.; Yang, X.; Ayers, R.A.; Moffat, K. Structure and function of plant photoreceptors. Annu. Rev. Plant Biol. 2010, 61, 21–47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Endo, M.; Araki, T.; Nagatani, A. Tissue-specific regulation of flowering by photoreceptors. Cell Mol. Life Sci. 2016, 73, 829–839. [Google Scholar] [CrossRef] [PubMed]
  43. Han, S.H.; Yoo, S.C.; Lee, B.D.; An, G.; Paek, N.C. Rice FLAVIN-BINDING, KELCH REPEAT, F-BOX 1 (OsFKF1) promotes flowering independent of photoperiod. Plant Cell Environ. 2015, 38, 2527–2540. [Google Scholar] [CrossRef] [PubMed]
  44. Xue, Z.G.; Zhang, X.M.; Lei, C.F.; Chen, X.J.; Fu, Y.F. Molecular cloning and functional analysis of one ZEITLUPE homolog GmZTL3 in soybean. Mol. Biol. Rep. 2012, 39, 1411–1418. [Google Scholar] [CrossRef] [PubMed]
  45. Jia, T.; Wei, D.; Meng, S.; Allan, A.C.; Zeng, L. Identification of regulatory genes implicated in continuous flowering of longan (Dimocarpus longan L.). PLoS ONE 2014, 9, 1–24. [Google Scholar]
  46. Zeng, J.; Wu, A.; Yang, Y.; Zhao, L.; Chen, L.; Wang, C.; Yang, X.; Chen, Y.; Liu, S.; Luo, P.; et al. The expression and interaction proteins analysis of BjuFKF1/LKP2 in B. juncea. Mol. Genet. Genom. 2022, 297, 75–85. [Google Scholar] [CrossRef]
  47. Giliberto, L.; Perrotta, G.; Pallara, P.; Weller, J.L.; Fraser, P.D.; Bramley, P.M.; Fiore, A.; Tavazza, M.; Giuliano, G. Manipulation of the blue light photoreceptor cryptochrome 2 in tomato affects vegetative development, flowering time, and fruit antioxidant content. Plant Physiol. 2005, 137, 199–208. [Google Scholar] [CrossRef] [Green Version]
  48. Sawa, M.; Nusinow, D.A.; Kay, S.A.; Imaizumi, T. FKF1 and GIGANTEA complex formation is required for day-length measurement in Arabidopsis. Science 2007, 318, 261–265. [Google Scholar] [CrossRef] [Green Version]
  49. Sawa, M.; Kay, S.A. GIGANTEA directly activates Flowering Locus T in Arabidopsis thaliana. Proc. Natl. Acad. Sci. USA 2011, 108, 11698–11703. [Google Scholar] [CrossRef] [Green Version]
  50. Fleet, C.M.; Sun, T.P. A DELLAcate balance: The role of gibberellin in plant morphogenesis. Curr. Opin. Plant Biol. 2005, 8, 77–85. [Google Scholar] [CrossRef]
  51. Kim, E.J.; Russinova, E. Brassinosteroid signalling. Curr. Biol. 2020, 30, R294–R298. [Google Scholar] [CrossRef] [PubMed]
  52. Helliwell, C.A.; Wood, C.C.; Robertson, M.; James Peacock, W.; Dennis, E.S. The Arabidopsis FLC protein interacts directly in vivo with SOC1 and FT chromatin and is part of a high-molecular-weight protein complex. Plant J. 2006, 46, 183–192. [Google Scholar] [CrossRef] [PubMed]
  53. Ballerini, E.S.; Kramer, E.M. Environmental and molecular analysis of the floral transition in the lower eudicot Aquilegia formosa. EvoDevo 2011, 2, 1–20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Whittaker, C.; Dean, C. The FLC Locus: A Platform for Discoveries in Epigenetics and Adaptation. Annu. Rev. Cell Dev. Biol. 2017, 33, 555–575. [Google Scholar] [CrossRef]
  55. Calderwood, A.; Lloyd, A.; Hepworth, J.; Tudor, E.H.; Jones, D.M.; Woodhouse, S.; Bilham, L.; Chinoy, C.; Williams, K.; Corke, F.; et al. Total FLC transcript dynamics from divergent paralogue expression explains flowering diversity in Brassica napus. New Phytol. 2021, 229, 3534–3548. [Google Scholar] [CrossRef]
  56. Marquardt, S.; Boss, P.K.; Hadfield, J.; Dean, C. Additional targets of the Arabidopsis autonomous pathway members, FCA and FY. J. Exp. Bot. 2006, 57, 3379–3386. [Google Scholar] [CrossRef] [Green Version]
  57. Jang, Y.H.; Park, H.Y.; Kim, S.K.; Lee, J.H.; Suh, M.C.; Chung, Y.S.; Paek, K.H.; Kim, J.K. Survey of rice proteins interacting with OsFCA and OsFY proteins which are homologous to the Arabidopsis flowering time proteins, FCA and FY. Plant Cell Physiol. 2009, 50, 1479–1492. [Google Scholar] [CrossRef]
  58. Raman, H.; Raman, R.; Coombes, N.; Song, J.; Prangnell, R.; Bandaranayake, C.; Tahira, R.; Sundaramoorthi, V.; Killian, A.; Meng, J.; et al. Genome-wide association analyses reveal complex genetic architecture underlying natural variation for flowering time in canola. Plant Cell Environ. 2016, 39, 1228–1239. [Google Scholar] [CrossRef]
  59. Sun, X.; Qin, Q.; Zhang, J.; Zhang, C.; Zhou, M.; Paek, K.Y.; Cui, Y. Isolation and characterization of the FVE gene of a Doritaenopsis hybrid involved in the regulation of flowering. Plant Growth Regul. 2012, 68, 77–86. [Google Scholar] [CrossRef]
  60. Hu, Q.; Jin, Y.; Shi, H.; Yang, W. GmFLD, a soybean homolog of the autonomous pathway gene FLOWERING LOCUS D, promotes flowering in Arabidopsis thaliana. BMC Plant Biol. 2014, 14, 263. [Google Scholar] [CrossRef] [Green Version]
  61. de Wit, M.; Galvao, V.C.; Fankhauser, C. Light-Mediated Hormonal Regulation of Plant Growth and Development. Annu. Rev. Plant Biol. 2016, 67, 513–537. [Google Scholar] [CrossRef] [PubMed]
  62. Matias-Hernandez, L.; Aguilar-Jaramillo, A.E.; Marin-Gonzalez, E.; Suarez-Lopez, P.; Pelaz, S. RAV genes: Regulation of floral induction and beyond. Ann. Bot. 2014, 114, 1459–1470. [Google Scholar] [CrossRef] [PubMed]
  63. Lee, J.H.; Yoo, S.J.; Park, S.H.; Hwang, I.; Lee, J.S.; Ahn, J.H. Role of SVP in the control of flowering time by ambient temperature in Arabidopsis. Genes Dev. 2007, 21, 397–402. [Google Scholar] [CrossRef] [Green Version]
  64. Xie, L.; Zhang, Y.; Wang, K.; Luo, X.; Xu, D.; Tian, X.; Li, L.; Ye, X.; Xia, X.; Li, W.; et al. TaVrt2, an SVP-like gene, cooperates with TaVrn1 to regulate vernalization-induced flowering in wheat. New Phytol. 2021, 231, 834–848. [Google Scholar] [CrossRef] [PubMed]
  65. Wu, R.; Tomes, S.; Karunairetnam, S.; Tustin, S.D.; Hellens, R.P.; Allan, A.C.; Macknight, R.C.; Varkonyi-Gasic, E. SVP-like MADS Box Genes Control Dormancy and Budbreak in Apple. Front. Plant Sci. 2017, 8, 477. [Google Scholar] [CrossRef]
  66. Li, Y.; Zhou, Y.; Yang, W.; Cheng, T.; Wang, J.; Zhang, Q. Isolation and functional characterization of SVP-like genes in Prunus mume. Sci. Hortic. 2017, 215, 91–101. [Google Scholar] [CrossRef]
  67. Wang, J.; Jiu, S.; Xu, Y.; Sabir, I.A.; Wang, L.; Ma, C.; Xu, W.; Wang, S.; Zhang, C. SVP-like gene PavSVP potentially suppressing flowering with PavSEP, PavAP1, and PavJONITLESS in sweet cherries (Prunus avium L.). Plant Physiol. Biochem. 2021, 159, 277–284. [Google Scholar] [CrossRef]
  68. Falavigna, V.D.S.; Guitton, B.; Costes, E.; Andres, F. I Want to (Bud) Break Free: The Potential Role of DAM and SVP-Like Genes in Regulating Dormancy Cycle in Temperate Fruit Trees. Front. Plant Sci. 2018, 9, 1990. [Google Scholar] [CrossRef] [Green Version]
  69. Li, L.; Zhang, C.; Huang, J.; Liu, Q.; Wei, H.; Wang, H.; Liu, G.; Gu, L.; Yu, S. Genomic analyses reveal the genetic basis of early maturity and identification of loci and candidate genes in upland cotton (Gossypium hirsutum L.). Plant Biotechnol. J. 2021, 19, 109–123. [Google Scholar] [CrossRef]
  70. Cheng, S.; Chen, P.; Su, Z.; Ma, L.; Hao, P.; Zhang, J.; Ma, Q.; Liu, G.; Liu, J.; Wang, H.; et al. High-resolution temporal dynamic transcriptome landscape reveals a GhCAL-mediated flowering regulatory pathway in cotton (Gossypium hirsutum L.). Plant Biotechnol. J. 2021, 19, 153–166. [Google Scholar] [CrossRef]
  71. Bouche, F.; Lobet, G.; Tocquin, P.; Perilleux, C. FLOR-ID: An interactive database of flowering-time gene networks in Arabidopsis thaliana. Nucleic Acids Res. 2016, 44, D1167–D1171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Hyun, Y.; Richter, R.; Vincent, C.; Martinez-Gallegos, R.; Porri, A.; Coupland, G. Multi-layered Regulation of SPL15 and Cooperation with SOC1 Integrate Endogenous Flowering Pathways at the Arabidopsis Shoot Meristem. Dev. Cell 2016, 37, 254–266. [Google Scholar] [CrossRef] [PubMed]
  73. Tao, Z.; Shen, L.; Liu, C.; Liu, L.; Yan, Y.; Yu, H. Genome-wide identification of SOC1 and SVP targets during the floral transition in Arabidopsis. Plant J. 2012, 70, 549–561. [Google Scholar] [CrossRef]
  74. Mandel, M.A.; Yanofsky, M.F. A gene triggering flower formation in Arabidopsis. Nature 1995, 377, 522–524. [Google Scholar] [CrossRef]
  75. Han, Y.; Zhang, C.; Yang, H.; Jiao, Y. Cytokinin pathway mediates APETALA1 function in the establishment of determinate floral meristems in Arabidopsis. Proc. Natl. Acad. Sci. USA 2014, 111, 6840–6845. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Jin, R.; Klasfeld, S.; Zhu, Y.; Fernandez Garcia, M.; Xiao, J.; Han, S.K.; Konkol, A.; Wagner, D. LEAFY is a pioneer transcription factor and licenses cell reprogramming to floral fate. Nat. Commun. 2021, 12, 626. [Google Scholar] [CrossRef]
  77. Engelhorn, J.; Moreau, F.; Fletcher, J.C.; Carles, C.C. ULTRAPETALA1 and LEAFY pathways function independently in specifying identity and determinacy at the Arabidopsis floral meristem. Ann. Bot. 2014, 114, 1497–1505. [Google Scholar] [CrossRef] [Green Version]
  78. Weigel, D.; Alvarez, J.; Smyth, D.R.; Yanofsky, M.F.; Meyerowitz, E.M. LEAFY Controls Floral Meristem Identity in Arabidopsis. Cell 1992, 69, 843–859. [Google Scholar] [CrossRef] [Green Version]
  79. Gregis, V.; Andrés, F.; Sessa, A.; Guerra, R.F.; Simonini, S.; Mateos, J.L.; Torti, S.; Zambelli, F.; Prazzoli, G.M.; Bjerkan, K.N.; et al. Identification of pathways directly regulated by SHORT VEGETATIVE PHASE during vegetative and reproductive development in Arabidopsis. Genome Biol. 2013, 14, R56. [Google Scholar] [CrossRef] [Green Version]
  80. Zuo, Z.; Liu, H.; Liu, B.; Liu, X.; Lin, C. Blue light-dependent interaction of CRY2 with SPA1 regulates COP1 activity and floral initiation in Arabidopsis. Curr. Biol. 2011, 21, 841–847. [Google Scholar] [CrossRef] [Green Version]
  81. Yamaguchi, S. Gibberellin metabolism and its regulation. Annu. Rev. Plant Biol. 2008, 59, 225–251. [Google Scholar] [CrossRef] [PubMed]
  82. Griffiths, J.; Murase, K.; Rieu, I.; Zentella, R.; Zhang, Z.L.; Powers, S.J.; Gong, F.; Phillips, A.L.; Hedden, P.; Sun, T.P.; et al. Genetic characterization and functional analysis of the GID1 gibberellin receptors in Arabidopsis. Plant Cell 2006, 18, 3399–3414. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Murase, K.; Hirano, Y.; Sun, T.P.; Hakoshima, T. Gibberellin-induced DELLA recognition by the gibberellin receptor GID1. Nature 2008, 456, 459–463. [Google Scholar] [CrossRef] [PubMed]
  84. Willige, B.C.; Ghosh, S.; Nill, C.; Zourelidou, M.; Dohmann, E.M.; Maier, A.; Schwechheimer, C. The DELLA domain of GA INSENSITIVE mediates the interaction with the GA INSENSITIVE DWARF1A gibberellin receptor of Arabidopsis. Plant Cell 2007, 19, 1209–1220. [Google Scholar] [CrossRef] [Green Version]
  85. Gallego-Bartolome, J.; Minguet, E.G.; Marin, J.A.; Prat, S.; Blazquez, M.A.; Alabadi, D. Transcriptional diversification and functional conservation between DELLA proteins in Arabidopsis. Mol. Biol. Evol. 2010, 27, 1247–1256. [Google Scholar] [CrossRef] [Green Version]
  86. Lee, S.; Cheng, H.; King, K.E.; Wang, W.; He, Y.; Hussain, A.; Lo, J.; Harberd, N.P.; Peng, J. Gibberellin regulates Arabidopsis seed germination via RGL2, a GAI/RGA-like gene whose expression is up-regulated following imbibition. Genes Dev. 2002, 16, 646–658. [Google Scholar] [CrossRef] [Green Version]
  87. Wang, J.; Qiu, Y.; Cheng, F.; Chen, X.; Zhang, X.; Wang, H.; Song, J.; Duan, M.; Yang, H.; Li, X. Genome-wide identification, characterization, and evolutionary analysis of flowering genes in radish (Raphanus sativus L.). BMC Genom. 2017, 18, 981. [Google Scholar] [CrossRef] [Green Version]
  88. Wang, M.; Tu, L.; Yuan, D.; Zhu, D.; Shen, C.; Li, J.; Liu, F.; Pei, L.; Wang, P.; Zhao, G.; et al. Reference genome sequences of two cultivated allotetraploid cottons, Gossypium hirsutum and Gossypium barbadense. Nat. Genet. 2019, 51, 224–229. [Google Scholar] [CrossRef] [Green Version]
  89. Fang, L.; Wang, Q.; Hu, Y.; Jia, Y.; Chen, J.; Liu, B.; Zhang, Z.; Guan, X.; Chen, S.; Zhou, B.; et al. Genomic analyses in cotton identify signatures of selection and loci associated with fiber quality and yield traits. Nat. Genet. 2017, 49, 1089–1098. [Google Scholar] [CrossRef]
  90. Huang, G.; Wu, Z.; Percy, R.G.; Bai, M.; Li, Y.; Frelichowski, J.E.; Hu, J.; Wang, K.; Yu, J.Z.; Zhu, Y. Genome sequence of Gossypium herbaceum and genome updates of Gossypium arboreum and Gossypium hirsutum provide insights into cotton A-genome evolution. Nat. Genet. 2020, 52, 516–524. [Google Scholar] [CrossRef] [Green Version]
  91. Nelson, D.C.; Lasswell, J.; Rogg, L.E.; Cohen, M.A.; Bartel, B. FKF1, a Clock-Controlled Gene that Regulates the Transition to Flowering in Arabidopsis. Cell 2000, 101, 331–340. [Google Scholar] [CrossRef] [Green Version]
  92. Baudry, A.; Ito, S.; Song, Y.H.; Strait, A.A.; Kiba, T.; Lu, S.; Henriques, R.; Pruneda-Paz, J.L.; Chua, N.H.; Tobin, E.M.; et al. F-box proteins FKF1 and LKP2 act in concert with ZEITLUPE to control Arabidopsis clock progression. Plant Cell 2010, 22, 606–622. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Corbesier, L.; Vincent, C.; Jang, S.; Fornara, F.; Fan, C.; Searle, I.; Giakountis, A.; Farrona, S.; Gissot, L.; Turnbull, C.; et al. FT protein movement contributes to long-distance signaling in floral induction of Arabidopsis. Science 2007, 316, 1030–1033. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijms 23 07940 g001 550
Figure 1. Identification and classification of flowering-related genes in G. hirsutum, G. barbadense, and G. arboretum. (A) G. hirsutum. (B) G. barbadense. (C) G. arboretum. The histogram on the left represents the number of flowering-related genes of each flowering-regulated pathway, a single black point in the middle matrix represents the number of genes specific to a pathway, the lines between the black points represent the intersection of the different pathways, and the vertical blue histogram represents the corresponding intersecting gene numbers. Ph: photoperiodism pathway; Au: autonomous pathway; Ho: hormone pathway; Ve: vernalization; Ag: aging pathway; Su: sugar signal; Fd: flower development and apical meristem response pathway; At: ambient temperature pathway.
Figure 1. Identification and classification of flowering-related genes in G. hirsutum, G. barbadense, and G. arboretum. (A) G. hirsutum. (B) G. barbadense. (C) G. arboretum. The histogram on the left represents the number of flowering-related genes of each flowering-regulated pathway, a single black point in the middle matrix represents the number of genes specific to a pathway, the lines between the black points represent the intersection of the different pathways, and the vertical blue histogram represents the corresponding intersecting gene numbers. Ph: photoperiodism pathway; Au: autonomous pathway; Ho: hormone pathway; Ve: vernalization; Ag: aging pathway; Su: sugar signal; Fd: flower development and apical meristem response pathway; At: ambient temperature pathway.
Ijms 23 07940 g001
Ijms 23 07940 g002 550
Figure 2. Distribution of flowering-related genes on G. hirsutum chromosomes. (A) Brownish bars: pseudochromosomes of the At subgenome; green bars: pseudochromosomes of the Dt subgenome. (B) Black lines on brownish and green bars indicate the locations of flowering-related genes on the pseudochromosomes. Values corresponding to the scales on the black vertical line indicate the physical distances.
Figure 2. Distribution of flowering-related genes on G. hirsutum chromosomes. (A) Brownish bars: pseudochromosomes of the At subgenome; green bars: pseudochromosomes of the Dt subgenome. (B) Black lines on brownish and green bars indicate the locations of flowering-related genes on the pseudochromosomes. Values corresponding to the scales on the black vertical line indicate the physical distances.
Ijms 23 07940 g002
Ijms 23 07940 g003 550
Figure 3. Distribution of the conserved motifs and exon–introns of photoperiod pathway flowering-related genes in G. hirsutum. (A) Phylogenetic tree of photoperiod pathway flowering-related genes in G. hirsutum constructed via the MEGA 6.0 program with the maximum likelihood (ML) method and 1000 bootstrap replicates. (B) Exon–intron structure of G. hirsutum flowering-related genes involved in the Ph pathway, as revealed by TBtools software. (C) Distribution of the conserved motifs of the flowering-related genes; the conserved motifs are indicated by the colored boxes. (D) Expression profiles of flowering-related genes. The expression profiles of flowering-related genes were detected in 15 different tissues: root, stigma, ovule 0D, ovule 1D, ovule 3D, ovule 10D, ovule 20D, petal, anther, fiber 5D, fiber 10D, fiber 20D, fiber 25D, stem and leaf tissues. Red arrows indicate the tissues in which the genes were highly expressed. D, day.
Figure 3. Distribution of the conserved motifs and exon–introns of photoperiod pathway flowering-related genes in G. hirsutum. (A) Phylogenetic tree of photoperiod pathway flowering-related genes in G. hirsutum constructed via the MEGA 6.0 program with the maximum likelihood (ML) method and 1000 bootstrap replicates. (B) Exon–intron structure of G. hirsutum flowering-related genes involved in the Ph pathway, as revealed by TBtools software. (C) Distribution of the conserved motifs of the flowering-related genes; the conserved motifs are indicated by the colored boxes. (D) Expression profiles of flowering-related genes. The expression profiles of flowering-related genes were detected in 15 different tissues: root, stigma, ovule 0D, ovule 1D, ovule 3D, ovule 10D, ovule 20D, petal, anther, fiber 5D, fiber 10D, fiber 20D, fiber 25D, stem and leaf tissues. Red arrows indicate the tissues in which the genes were highly expressed. D, day.
Ijms 23 07940 g003
Ijms 23 07940 g004 550
Figure 4. Distribution of conserved motifs and exon–intron structure, and expression level analysis of Ho pathway flowering-related genes. (A) Exon–intron structure of G. hirsutum flowering-related genes involved in the Ho pathway, as revealed by TBtools software. (B) Distribution of conserved motifs of flowering-related genes; conserved motifs are indicated by colored boxes. (C) Unrooted phylogenetic tree of Ho pathway flowering-related genes in G. hirsutum according to the MEGA 6.0 program with the maximum likelihood (ML) method and 1000 bootstrap replicates. (D) Tissue expression patterns of key genes involved in the G. hirsutum photoperiod pathway. The abscissa represents different tissues, and the rightmost ordinate represents the gene ID and name. Red arrows indicate the tissues in which the genes were highly expressed. (E–H) Changes in hormone contents ((E) ZR; (F) ABA; (G) GA3; (H) IAA) in the SAM and leaves of cotton at different true-leaf flattening stages (first true-leaf stage 1TLS, second true-leaf stage 2TLS, third true-leaf stage 3TLS, fourth true-leaf stage 4TLS). The abscissa represents the development stages for cotton. o, cotyledon period; 1, the first true-leaf stage; 2, the second true-leaf stage; 3, the third true-leaf stage; 4, the fourth true-leaf stage. The vertical axis represents hormone content. Two early-maturing cultivars: Chinese Cotton Research Institute 50 (CCRI50, also named Zhong50) and Chinese Cotton Research Institute 74 (CCRI74, also named Zhong74); two late-maturing cultivars: Chinese Cotton Research Institute 60 (CCRI60, also named Zhong60) and Luyanmian28 (Lu28). SAM, shoot apical meristem.
Figure 4. Distribution of conserved motifs and exon–intron structure, and expression level analysis of Ho pathway flowering-related genes. (A) Exon–intron structure of G. hirsutum flowering-related genes involved in the Ho pathway, as revealed by TBtools software. (B) Distribution of conserved motifs of flowering-related genes; conserved motifs are indicated by colored boxes. (C) Unrooted phylogenetic tree of Ho pathway flowering-related genes in G. hirsutum according to the MEGA 6.0 program with the maximum likelihood (ML) method and 1000 bootstrap replicates. (D) Tissue expression patterns of key genes involved in the G. hirsutum photoperiod pathway. The abscissa represents different tissues, and the rightmost ordinate represents the gene ID and name. Red arrows indicate the tissues in which the genes were highly expressed. (E–H) Changes in hormone contents ((E) ZR; (F) ABA; (G) GA3; (H) IAA) in the SAM and leaves of cotton at different true-leaf flattening stages (first true-leaf stage 1TLS, second true-leaf stage 2TLS, third true-leaf stage 3TLS, fourth true-leaf stage 4TLS). The abscissa represents the development stages for cotton. o, cotyledon period; 1, the first true-leaf stage; 2, the second true-leaf stage; 3, the third true-leaf stage; 4, the fourth true-leaf stage. The vertical axis represents hormone content. Two early-maturing cultivars: Chinese Cotton Research Institute 50 (CCRI50, also named Zhong50) and Chinese Cotton Research Institute 74 (CCRI74, also named Zhong74); two late-maturing cultivars: Chinese Cotton Research Institute 60 (CCRI60, also named Zhong60) and Luyanmian28 (Lu28). SAM, shoot apical meristem.
Ijms 23 07940 g004
Ijms 23 07940 g005 550
Figure 5. Direction and magnitude of natural selection acting on different flowering-related gene sets. Ka/Ks values of flowering-related genes in G. hirsutum (A) and G. barbadense (B). Quantile boxplots showing the distribution of Ka/Ks values for homologous gene pairs. Horizontal bars in each box indicate median values. Upper and lower bars correspond to the upper and lower adjacent values 1.5 times past the interquartile range. The outliers are plotted as discrete dots. Ve, vernalization; Su, sugar; Ph, photoperiod pathway, circadian clock, and light signaling; Ho, hormone signaling and metabolism; Au, autonomous pathways; At, ambient temperature pathway; Ag, aging pathway; Fd, flower development and apical meristem response pathway.
Figure 5. Direction and magnitude of natural selection acting on different flowering-related gene sets. Ka/Ks values of flowering-related genes in G. hirsutum (A) and G. barbadense (B). Quantile boxplots showing the distribution of Ka/Ks values for homologous gene pairs. Horizontal bars in each box indicate median values. Upper and lower bars correspond to the upper and lower adjacent values 1.5 times past the interquartile range. The outliers are plotted as discrete dots. Ve, vernalization; Su, sugar; Ph, photoperiod pathway, circadian clock, and light signaling; Ho, hormone signaling and metabolism; Au, autonomous pathways; At, ambient temperature pathway; Ag, aging pathway; Fd, flower development and apical meristem response pathway.
Ijms 23 07940 g005
Ijms 23 07940 g006 550
Figure 6. Possible regulatory model of the multiple feedbacks and inputs during the regulation of cotton flowering in the leaves and SAM.
Figure 6. Possible regulatory model of the multiple feedbacks and inputs during the regulation of cotton flowering in the leaves and SAM.
Ijms 23 07940 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, X.; Wu, Y.; Chi, H.; Wei, H.; Wang, H.; Yu, S. Genomewide Identification and Characterization of the Genes Involved in the Flowering of Cotton. Int. J. Mol. Sci. 2022, 23, 7940. https://doi.org/10.3390/ijms23147940

AMA Style

Li X, Wu Y, Chi H, Wei H, Wang H, Yu S. Genomewide Identification and Characterization of the Genes Involved in the Flowering of Cotton. International Journal of Molecular Sciences. 2022; 23(14):7940. https://doi.org/10.3390/ijms23147940

Chicago/Turabian Style

Li, Xiao, Yuanlong Wu, Huabin Chi, Hengling Wei, Hantao Wang, and Shuxun Yu. 2022. "Genomewide Identification and Characterization of the Genes Involved in the Flowering of Cotton" International Journal of Molecular Sciences 23, no. 14: 7940. https://doi.org/10.3390/ijms23147940

APA Style

Li, X., Wu, Y., Chi, H., Wei, H., Wang, H., & Yu, S. (2022). Genomewide Identification and Characterization of the Genes Involved in the Flowering of Cotton. International Journal of Molecular Sciences, 23(14), 7940. https://doi.org/10.3390/ijms23147940

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop