A pair of readers of bivalent chromatin mediate formation of Polycomb-based "memory of cold" in plants.

The Polycomb-group chromatin modifiers play important roles to repress or switch off gene expression in plants and animals. How the active chromatin state is switched to a Polycomb-repressed state is unclear. In Arabidopsis, prolonged cold induces the switching of the highly active chromatin state at the potent floral repressor FLC to a Polycomb-repressed state, which is epigenetically maintained when temperature rises to confer "cold memory," enabling plants to flower in spring. We report that the cis-acting cold memory element (CME) region at FLC bears bivalent marks of active histone H3K4me3 and repressive H3K27me3 that are read and interpreted by an assembly of bivalent chromatin readers to drive cold-induced switching of the FLC chromatin state. In response to cold, the 47-bp CME and its associated bivalent chromatin feature drive the switching of active chromatin state at a recombinant gene to a Polycomb-repressed domain, conferring cold memory. We reveal a paradigm for environment-induced chromatin-state switching at bivalent loci in plants.

In vivo single-molecule RNA structural profiling.

Cellular RNAs exhibit substantial heterogeneity in structure and function. Recently, Yang et al. developed an in vivo single-molecule RNA structure profiling methodology and revealed that individual isoforms of noncoding transcripts adopt multiple diverse and functionally relevant structural conformations, which change in abundance and structure in response to temperature conditions.

Genome evolution of the ancient hexaploid Platanus × acerifolia (London planetree).

Whole-genome duplication (WGD; i.e., polyploidy) and chromosomal rearrangement (i.e., genome shuffling) significantly influence genome structure and organization. Many polyploids show extensive genome shuffling relative to their pre-WGD ancestors. No reference genome is currently available for Platanaceae (Proteales), one of the sister groups to the core eudicots. Moreover, Platanus × acerifolia (London planetree; Platanaceae) is a widely used street tree. Given the pivotal phylogenetic position of Platanus and its 2-y flowering transition, understanding its flowering-time regulatory mechanism has significant evolutionary implications; however, the impact of Platanus genome evolution on flowering-time genes remains unknown. Here, we assembled a high-quality, chromosome-level reference genome for P. × acerifolia using a phylogeny-based subgenome phasing method. Comparative genomic analyses revealed that P. × acerifolia (2n = 42) is an ancient hexaploid with three subgenomes resulting from two sequential WGD events; Platanus does not seem to share any WGD with other Proteales or with core eudicots. Each P. × acerifolia subgenome is highly similar in structure and content to the reconstructed pre-WGD ancestral eudicot genome without chromosomal rearrangements. The P. × acerifolia genome exhibits karyotypic stasis and gene sub-/neo-functionalization and lacks subgenome dominance. The copy number of flowering-time genes in P. × acerifolia has undergone an expansion compared to other noncore eudicots, mainly via the WGD events. Sub-/neo-functionalization of duplicated genes provided the genetic basis underlying the unique flowering-time regulation in P. × acerifolia. The P. × acerifolia reference genome will greatly expand understanding of the evolution of genome organization, genetic diversity, and flowering-time regulation in angiosperms.

Sex-dependent phenological responses to climate vary across species' ranges.

Anthropogenic climate change has significantly altered the flowering times (i.e., phenology) of plants worldwide, affecting their reproduction, survival, and interactions. Recent studies utilizing herbarium specimens have uncovered significant intra- and inter-specific variation in flowering phenology and its response to changes in climate but have mostly been limited to animal-pollinated species. Thus, despite their economic and ecological importance, variation in phenological responses to climate remain largely unexplored among and within wind-pollinated dioecious species and across their sexes. Using both herbarium specimens and volunteer observations of cottonwood (Populus) species, we examined how phenological sensitivity to climate varies across species, their ranges, sexes, and phenophases. The timing of flowering varied significantly across and within species, as did their sensitivity to spring temperature. In particular, male flowering generally happened earlier in the season and was more sensitive to warming than female flowering. Further, the onset of flowering was more sensitive to changes in temperature than leaf out. Increased temporal gaps between male and female flowering time and between the first open flower date and leaf out date were predicted for the future under two climate change scenarios. These shifts will impact the efficacy of sexual reproduction and gene flow among species. Our study demonstrates significant inter- and intra-specific variation in phenology and its responses to environmental cues, across species' ranges, phenophases, and sex, in wind-pollinated species. These variations need to be considered to predict accurately the effects of climate change and assess their ecological and evolutionary consequences.

Population-level gene expression can repeatedly link genes to functions in maize.

Transcriptome-wide association studies (TWAS) can provide single gene resolution for candidate genes in plants, complementing genome-wide association studies (GWAS) but efforts in plants have been met with, at best, mixed success. We generated expression data from 693 maize genotypes, measured in a common field experiment, sampled over a 2-h period to minimize diurnal and environmental effects, using full-length RNA-seq to maximize the accurate estimation of transcript abundance. TWAS could identify roughly 10 times as many genes likely to play a role in flowering time regulation as GWAS conducted data from the same experiment. TWAS using mature leaf tissue identified known true-positive flowering time genes known to act in the shoot apical meristem, and trait data from a new environment enabled the identification of additional flowering time genes without the need for new expression data. eQTL analysis of TWAS-tagged genes identified at least one additional known maize flowering time gene through trans-eQTL interactions. Collectively these results suggest the gene expression resource described here can link genes to functions across different plant phenotypes expressed in a range of tissues and scored in different experiments.

ATX1 and HUB1/2 promote recruitment of the transcription elongation factor VIP2 to modulate the floral transition in Arabidopsis.

Histone 2B ubiquitination (H2Bub) and trimethylation of H3 at lysine 4 (H3K4me3) are associated with transcription activation. However, the function of these modifications in transcription in plants remains largely unknown. Here, we report that coordination of H2Bub and H3K4me3 deposition with the binding of the RNA polymerase-associated factor VERNALIZATION INDEPENDENCE2 (VIP2) to FLOWERING LOCUS C (FLC) modulates flowering time in Arabidopsis. We found that RING domain protein HISTONE MONOUBIQUITINATION1 (HUB1) and HUB2 (we refer as HUB1/2), which are responsible for H2Bub, interact with ARABIDOPSIS TRITHORAX1 (ATX1), which is required for H3K4me3 deposition, to promote the transcription of FLC and repress the flowering time. The atx1-2 hub1-10 hub2-2 triple mutant in FRIGIDIA (FRI) background displayed early flowering like FRI hub1-10 hub2-2 and overexpression of ATX1 failed to rescue the early flowering phenotype of hub1-10 hub2-2. Mutations in HUB1 and HUB2 reduced the ATX1 enrichment at FLC, indicating that HUB1 and HUB2 are required for ATX1 recruitment and H3K4me3 deposition at FLC. We also found that the VIP2 directly binds to HUB1, HUB2, and ATX1 and that loss of VIP2 in FRI hub1-10 hub2-2 and FRI atx1-2 plants resulted in early flowering like that observed in FRI vip2-10. Loss of function of HUB2 and ATX1 impaired VIP2 enrichment at FLC, and reduced the transcription initiation and elongation of FLC. In addition, mutations in VIP2 reduced HUB1 and ATX1 enrichment and H2Bub and H3K4me3 levels at FLC. Together, our findings revealed that HUB1/2, ATX1, and VIP2 coordinately modulate H2Bub and H3K4me3 deposition, FLC transcription, and flowering time.

Loss of daylength sensitivity by splice site mutation in Cannabis pseudo-response regulator.

Photoperiod insensitivity (auto-flowering) in drug-type Cannabis sativa circumvents the need for short day (SD) flowering requirements making outdoor cultivation in high latitudes possible. However, the benefits of photoperiod insensitivity are counterbalanced by low cannabinoid content and poor flower quality in auto-flowering genotypes. Despite recent studies in cannabis flowering, a mechanistic understanding of photoperiod insensitivity is still lacking. We used a combination of genome-wide association study and genetic fine-mapping to identify the genetic cause of auto-flowering in cannabis. We then used gene expression analyses and transient transformation assays to characterize flowering time control. Herein, we identify a splice site mutation within circadian clock gene PSEUDO-RESPONSE REGULATOR 37 (CsPRR37) in auto-flowering cannabis. We show that CsPRR37 represses FT expression and its circadian oscillations transition to a less repressive state during SD as compared to long days (LD). We identify several key circadian clock genes whose expression is altered in auto-flowering cannabis, particularly under non-inductive LD. Research into the pervasiveness of this mutation and others affecting flowering time will help elucidate cannabis domestication history and advance cannabis breeding toward a more sustainable outdoor cultivation system.

A FLOWERING LOCUS T ortholog is associated with photoperiod-insensitive flowering in hemp (Cannabis sativa L.).

Hemp (Cannabis sativa L.) is an extraordinarily versatile crop, with applications ranging from medicinal compounds to seed oil and fibre products. Cannabis sativa is a short-day plant, and its flowering is highly controlled by photoperiod. However, substantial genetic variation exists for photoperiod sensitivity in C. sativa, and photoperiod-insensitive ("autoflower") cultivars are available. Using a bi-parental mapping population and bulked segregant analysis, we identified Autoflower2, a 0.5 Mbp locus significantly associated with photoperiod-insensitive flowering in hemp. Autoflower2 contains an ortholog of the central flowering time regulator FLOWERING LOCUS T (FT) from Arabidopsis thaliana which we termed CsFT1. We identified extensive sequence divergence between alleles of CsFT1 from photoperiod-sensitive and insensitive cultivars of C. sativa, including a duplication of CsFT1 and sequence differences, especially in introns. Furthermore, we observed higher expression of one of the CsFT1 copies found in the photoperiod-insensitive cultivar. Genotyping of several mapping populations and a diversity panel confirmed a correlation between CsFT1 alleles and photoperiod response, affirming that at least two independent loci involved in the photoperiodic control of flowering, Autoflower1 and Autoflower2, exist in the C. sativa gene pool. This study reveals the multiple independent origins of photoperiod insensitivity in C. sativa, supporting the likelihood of a complex domestication history in this species. By integrating the genetic relaxation of photoperiod sensitivity into novel C. sativa cultivars, expansion to higher latitudes will be permitted, thus allowing the full potential of this versatile crop to be reached.

Dynamic Phytomeric Growth Contributes to Local Adaptation in Barley.

Vascular plants have segmented body axes with iterative nodes and internodes. Appropriate node initiation and internode elongation are fundamental to plant fitness and crop yield; however, how these events are spatiotemporally coordinated remains elusive. We show that in barley (Hordeum vulgare L.), selections during domestication have extended the apical meristematic phase to promote node initiation, but constrained subsequent internode elongation. In both vegetative and reproductive phases, internode elongation displays a dynamic proximal-distal gradient, and among subpopulations of domesticated barleys worldwide, node initiation and proximal internode elongation are associated with latitudinal and longitudinal gradients, respectively. Genetic and functional analyses suggest that, in addition to their converging roles in node initiation, flowering-time genes have been repurposed to specify the timing and duration of internode elongation. Our study provides an integrated view of barley node initiation and internode elongation and suggests that plant architecture should be recognized as a collection of dynamic phytomeric units in the context of crop adaptive evolution.

OsFLZ2 interacts with OsMADS51 to fine-tune rice flowering time.

Flowering time is an important agronomic trait affecting crop yield. FCS-LIKE ZINC FINGER (FLZ) proteins are plant-specific regulatory proteins that are involved in multiple biological processes. However, their roles in plant flowering time control have not been clarified. Here, we report that OsFLZ2 is a negative regulator of rice flowering time. OsFLZ2 delays flowering by repressing the expression of key floral integrator genes. Biochemical assays showed OsFLZ2 physically interacts with OsMADS51, a flowering activator under short-day (SD) conditions. Both OsFLZ2 and OsMADS51 are highly expressed in rice leaves before floral transition under natural SD conditions, and their proteins are colocalized in the nucleus. Co-expression of OsFLZ2 can destabilize OsMADS51 and weaken its transcriptional activation of the downstream target gene Early heading date 1 (Ehd1). Taken together, these results indicate that OsFLZ2 can interfere with the function of OsMADS51 to fine-tune rice flowering time.

An adaptive teosinte mexicana introgression modulates phosphatidylcholine levels and is associated with maize flowering time.

Native Americans domesticated maize (Zea mays ssp. mays) from lowland teosinte parviglumis (Zea mays ssp. parviglumis) in the warm Mexican southwest and brought it to the highlands of Mexico and South America where it was exposed to lower temperatures that imposed strong selection on flowering time. Phospholipids are important metabolites in plant responses to low-temperature and phosphorus availability and have been suggested to influence flowering time. Here, we combined linkage mapping with genome scans to identify High PhosphatidylCholine 1 (HPC1), a gene that encodes a phospholipase A1 enzyme, as a major driver of phospholipid variation in highland maize. Common garden experiments demonstrated strong genotype-by-environment interactions associated with variation at HPC1, with the highland HPC1 allele leading to higher fitness in highlands, possibly by hastening flowering. The highland maize HPC1 variant resulted in impaired function of the encoded protein due to a polymorphism in a highly conserved sequence. A meta-analysis across HPC1 orthologs indicated a strong association between the identity of the amino acid at this position and optimal growth in prokaryotes. Mutagenesis of HPC1 via genome editing validated its role in regulating phospholipid metabolism. Finally, we showed that the highland HPC1 allele entered cultivated maize by introgression from the wild highland teosinte Zea mays ssp. mexicana and has been maintained in maize breeding lines from the Northern United States, Canada, and Europe. Thus, HPC1 introgressed from teosinte mexicana underlies a large metabolic QTL that modulates phosphatidylcholine levels and has an adaptive effect at least in part via induction of early flowering time.

The SOC1 gene plays an important role in regulating litchi flowering time.

Litchi (Litchi chinensis Sonn.) is a valuable subtropical fruit tree with high-quality fruit. However, its economic benefits and sustainable development are restrained by a number of challenges. One major challenge is the lack of extremely early and late maturing high-quality varieties due to limited availability of varieties suitable for commercial cultivation and outdated breeding methods, resulting in an imbalanced supply and low price of litchi. Flowering time is a crucial genetic factor influencing the maturation period of litchi. Our previous research has highlighted the pivotal role of the LcFT1 gene in regulating the flowering time of litchi and identified a gene associated with LcFT1 (named as LcSOC1) based on RNA-Seq and weight gene co-expression network (WGCNA) analysis. This study further investigated the function of LcSOC1. Subcellular localization analysis revealed that LcSOC1 is primarily localized in the nucleus, where it acts as a transcription factor. LcSOC1 overexpression in Nicotiana tabacum and Arabidopsis thaliana resulted in significant early flowering. Furthermore, LcSOC1 was found to be expressed in various tissues, with the highest expression in mature leaves. Analysis of spatial and temporal expression patterns of LcSOC1 in litchi varieties with different flowering time under low temperature treatment and across an annual cycle demonstrated that LcSOC1 is responsive to low temperature induction. Interestingly, early maturing varieties exhibited higher sensitivity to low temperature, with significantly premature induction of LcSOC1 expression relative to late maturing varieties. Activation of LcSOC1 triggered the transition of litchi into the flowering phase. These findings demonstrate that LcSOC1 plays a pivotal role in regulating the flowering process and determining the flowering time in litchi. Overall, this study provides theoretical guidance and important target genes for molecular breeding to regulate litchi production period.

The Arabidopsis RRM domain proteins EDM3 and IBM2 coordinate the floral transition and basal immune responses.

The transition from vegetative to reproductive development (floral transition) is a costly process in annual plants requiring increased investments in metabolic resources. The Arabidopsis thaliana (Arabidopsis) PHD finger protein EDM2 and RRM domain proteins EDM3 and IBM2 are known to form chromatin-associated complexes controlling transcript processing. We are reporting that distinct splice isoforms of EDM3 and IBM2 cooperate in the coordination of the floral transition with basal immune responses. These cooperating splice isoforms, termed EDM3L and IBM2L, control the intensity of basal immunity and, via a separate pathway, the timing of the floral transition. During the developmental phase prior to the floral transition expression of EDM3L and IBM2L strongly and gradually increases, while these isoforms simultaneously down-regulate expression of the floral suppressor gene FLC and promote the transition to reproductive growth. At the same time these accumulating EDM3 and IBM2 splice isoforms gradually suppress basal immunity against the virulent Noco2 isolate of the pathogenic oomycete Hyaloperonospora arabidopsidis and down-regulate expression of a set of defense-associated genes and immune receptor genes. We are providing clear evidence for a functional link between the floral transition and basal immunity in the annual plant Arabidopsis. Coordination of these two biological processes, which compete for metabolic resources, is likely critical for plant survival and reproductive success.

Functional divergence of FTL9 and FTL10 in flowering control in rice.

BACKGROUND: Floral transition in cereals is a critical phenomenon influenced by exogenous and endogenous signals, determining crop yield and reproduction. Flowering Locus T-like (FT-like) genes encode a mobile florigen, the main signaling molecule for flowering. RESULTS: In this study, we characterized two FT-like genes, FTL9 and FTL10, to study their functional diversity in flowering control in rice. We compared independent mutant lines of ftl10 with WT and observed negligible differences in the flowering phenotype, or agronomic traits implying potentially redundant roles of FTL10 loss-of-function in flowering control in rice. Nevertheless, we found that overexpression of FTL10, but not FTL9, substantially accelerated flowering, indicating the flowering-promoting role of FTL10 and the divergent functions between FTL9 and FTL10 in flowering. Besides flowering, additive agronomic roles were observed for FTL10-OE regulating the number of effective panicles per plant, the number of primary branches per panicle, and spikelets per panicle without regulating seed size. Mechanistically, our Y2H and BiFC analyses demonstrate that FTL10, in contrast to FTL9, can interact with FD1 and GF14c, forming a flowering activation complex and thereby regulating flowering. CONCLUSION: Altogether, our results elucidate the regulatory roles of FTL9 and FTL10 in flowering control, unveiling the molecular basis of functional divergence between FTL10 and FTL9, which provides mechanistic insights into shaping the dynamics of flowering time regulation in rice.

Identification of vernalization-related genes and cold memory element (CME) required for vernalization response in radish (Raphanus sativus L.).

Floral transition is accelerated by exposure to long-term cold like winter in plants, which is called as vernalization. Acceleration of floral transition by vernalization is observed in a diversity of biennial and perennial plants including Brassicaceae family plants. Scientific efforts to understand molecular mechanism underlying vernalization-mediated floral transition have been intensively focused in model plant Arabidopsis thaliana. To get a better understanding on floral transition by vernalization in radish (Raphanus sativus L.), we investigated transcriptomic changes taking place during vernalization in radish. Thousands of genes were differentially regulated along time course of vernalization compared to non-vernalization (NV) sample. Twelve major clusters of DEGs were identified based on distinctive expression profiles during vernalization. Radish FLC homologs were shown to exert an inhibition of floral transition when transformed into Arabidopsis plants. In addition, DNA region containing RY motifs located within a Raphanus sativus FLC homolog, RsFLC1 was found to be required for repression of RsFLC1 by vernalization. Transgenic plants harboring disrupted RY motifs were impaired in the enrichment of H3K27me3 on RsFLC1 chromatin, thus resulting in the delayed flowering in Arabidopsis. Taken together, we report transcriptomic profiles of radish during vernalization and demonstrate the requirement of RY motif for vernalization-mediated repression of RsFLC homologs in radish (Raphanus sativus L.).

Genetic Architecture of Flowering Time Differs Between Populations With Contrasting Demographic and Selective Histories.

Understanding the evolutionary factors that impact the genetic architecture of traits is a central goal of evolutionary genetics. Here, we investigate how quantitative trait variation accumulated over time in populations that colonized a novel environment. We compare the genetic architecture of flowering time in Arabidopsis populations from the drought-prone Cape Verde Islands and their closest outgroup population from North Africa. We find that trait polygenicity is severely reduced in the island populations compared to the continental North African population. Further, trait architectures and reconstructed allelic histories best fit a model of strong directional selection in the islands in accord with a Fisher-Orr adaptive walk. Consistent with this, we find that large-effect variants that disrupt major flowering time genes (FRI and FLC) arose first, followed by smaller effect variants, including ATX2 L125F, which is associated with a 4-day reduction in flowering time. The most recently arising flowering time-associated loci are not known to be directly involved in flowering time, consistent with an omnigenic signature developing as the population approaches its trait optimum. Surprisingly, we find no effect in the natural population of EDI-Cvi-0 (CRY2 V367M), an allele for which an effect was previously validated by introgression into a Eurasian line. Instead, our results suggest the previously observed effect of the EDI-Cvi-0 allele on flowering time likely depends on genetic background, due to an epistatic interaction. Altogether, our results provide an empirical example of the effects demographic history and selection has on trait architecture.

Cantil: a previously unreported organ in wild-type Arabidopsis regulated by FT, ERECTA and heterotrimeric G proteins.

We describe a previously unreported macroscopic Arabidopsis organ, the cantil, named for its 'cantilever' function of holding the pedicel at a distance from the stem. Cantil development is strongest at the first nodes after the vegetative to reproductive inflorescence transition; cantil magnitude and frequency decrease acropetally. Cantils develop in wild-type Arabidopsis accessions (e.g. Col-0, Ws and Di-G) as a consequence of delayed flowering in short days; cantil formation is observed in long days when flowering is delayed by null mutation of the floral regulator FLOWERING LOCUS T. The receptor-like kinase ERECTA is a global positive regulator of cantil formation; therefore, cantils never form in the Arabidopsis strain Ler. ERECTA functions genetically upstream of heterotrimeric G proteins. Cantil expressivity is repressed by the specific heterotrimeric complex subunits GPA1, AGB1 and AGG3, which also play independent roles: GPA1 suppresses distal spurs at cantil termini, while AGB1 and AGG3 suppress ectopic epidermal rippling. These G protein mutant traits are recapitulated in long-day flowering gpa1-3 ft-10 plants, demonstrating that cantils, spurs and ectopic rippling occur as a function of delayed phase transition, rather than as a function of photoperiod per se.

Natural variation of domestication-related genes contributed to latitudinal expansion and adaptation in soybean.

Soybean is a major source of protein and edible oil worldwide. Originating from the Huang-Huai-Hai region, which has a temperate climate, soybean has adapted to a wide latitudinal gradient across China. However, the genetic mechanisms responsible for the widespread latitudinal adaptation in soybean, as well as the genetic basis, adaptive differentiation, and evolutionary implications of theses natural alleles, are currently lacking in comprehensive understanding. In this study, we examined the genetic variations of fourteen major gene loci controlling flowering and maturity in 103 wild species, 1048 landraces, and 1747 cultivated species. We found that E1, E3, FT2a, J, Tof11, Tof16, and Tof18 were favoured during soybean improvement and selection, which explained 75.5% of the flowering time phenotypic variation. These genetic variation was significantly associated with differences in latitude via the LFMM algorithm. Haplotype network and geographic distribution analysis suggested that gene combinations were associated with flowering time diversity contributed to the expansion of soybean, with more HapA clustering together when soybean moved to latitudes beyond 35°N. The geographical evolution model was developed to accurately predict the suitable planting zone for soybean varieties. Collectively, by integrating knowledge from genomics and haplotype classification, it was revealed that distinct gene combinations improve the adaptation of cultivated soybeans to different latitudes. This study provides insight into the genetic basis underlying the environmental adaptation of soybean accessions, which could contribute to a better understanding of the domestication history of soybean and facilitate soybean climate-smart molecular breeding for various environments.

WRKY transcription factors modulate flowering time in four Arachis species: a bioinformatics analysis.

BACKGROUND: WRKY proteins are important transcription factors (TFs) in plants, involved in growth and development and responses to environmental changes. Although WRKY TFs have been studied at the genome level in Arachis genus, including oil crop and turfgrass, their regulatory networks in controlling flowering time remain unclear. The aim of this study was to predict the molecular mechanisms of WRKY TFs regulation flowering time in Arachis genus at the genome level using bioinformatics approaches. RESULTS: The flowering-time genes of Arachis genus were retrieved from the flowering-time gene database. The regulatory networks between WRKY TFs and downstream genes in Arachis genus were predicted using bioinformatics tools. The results showed that WRKY TFs were involved in aging, autonomous, circadian clock, hormone, photoperiod, sugar, temperature, and vernalization pathways to modulate flowering time in Arachis duranensis, Arachis ipaensis, Arachis monticola, and Arachis hypogaea cv. Tifrunner. The WRKY TF binding sites in homologous flowering-time genes exhibited asymmetric evolutionary pattern, indicating that the WRKY TFs interact with other transcription factors to modulate flowering time in the four Arachis species. Protein interaction network analysis showed that WRKY TFs interacted with FRUITFULL and APETALA2 to modulate flowering time in the four Arachis species. WRKY TFs implicated in regulating flowering time had low expression levels, whereas their interaction proteins had varying expression patterns in 22 tissues of A. hypogaea cv. Tifrunner. These results indicate that WRKY TFs exhibit antagonistic or synergistic interactions with the associated proteins. CONCLUSIONS: This study reveals complex regulatory networks through which WRKY TFs modulate flowering time in the four Arachis species using bioinformatics approaches.

H3K36 methyltransferase GhKMT3;1a and GhKMT3;2a promote flowering in upland cotton.

BACKGROUND: The SET domain group (SDG) genes encode histone lysine methyltransferases, which regulate gene transcription by altering chromatin structure and play pivotal roles in plant flowering determination. However, few studies have investigated their role in the regulation of flowering in upland cotton. RESULTS: A total of 86 SDG genes were identified through genome-wide analysis in upland cotton (Gossypium hirsutum). These genes were unevenly distributed across 25 chromosomes. Cluster analysis revealed that the 86 GhSDGs were divided into seven main branches. RNA-seq data and qRT‒PCR analysis revealed that lysine methyltransferase 3 (KMT3) genes were expressed at high levels in stamens, pistils and other floral organs. Using virus-induced gene silencing (VIGS), functional characterization of GhKMT3;1a and GhKMT3;2a revealed that, compared with those of the controls, the GhKMT3;1a- and GhKMT3;2a-silenced plants exhibited later budding and flowering and lower plant heightwere shorter. In addition, the expression of flowering-related genes (GhAP1, GhSOC1 and GhFT) significantly decreased and the expression level of GhSVP significantly increased in the GhKMT3;1a- and GhKMT3;2a-silenced plants compared with the control plants. CONCLUSION: A total of 86 SDG genes were identified in upland cotton, among which GhKMT3;1a and GhKMT3;2a might regulate flowering by affecting the expression of GhAP1, GhSOC1, GhFT and GhSVP. These findings will provide genetic resources for advanced molecular breeding in the future.