Functional and Evolutionary Characterization of the CONSTANS Gene Family in Short-Day Photoperiodic Flowering in Soybean Faqiang Wu1., Brian William Price1., Waseem Haider1, Gabriela Seufferheld1, Randall Nelson2, Yoshie Hanzawa1* 1 Department of Crop Sciences, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America, 2 USDA-Agricultural Research Service, Soybean/ Maize Germplasm, Pathology, and Genetics Research Unit, Department of Crop Sciences, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America

Abstract CONSTANS (CO) plays a central role in photoperiodic flowering control of plants. However, much remains unknown about the function of the CO gene family in soybean and the molecular mechanisms underlying short-day photoperiodic flowering of soybean. We identified 26 CO homologs (GmCOLs) in the soybean genome, many of them previously unreported. Phylogenic analysis classified GmCOLs into three clades conserved among flowering plants. Two homeologous pairs in Clade I, GmCOL1a/GmCOL1b and GmCOL2a/GmCOL2b, showed the highest sequence similarity to Arabidopsis CO. The mRNA abundance of GmCOL1a and GmCOL1b exhibited a strong diurnal rhythm under flowering-inductive short days and peaked at dawn, which coincided with the rise of GmFT5a expression. In contrast, the mRNA abundance of GmCOL2a and GmCOL2b was extremely low. Our transgenic study demonstrated that GmCOL1a, GmCOL1b, GmCOL2a and GmCOL2b fully complemented the late flowering effect of the co-1 mutant in Arabidopsis. Together, these results indicate that GmCOL1a and GmCOL1b are potential inducers of flowering in soybean. Our data also indicate rapid regulatory divergence between GmCOL1a/GmCOL1b and GmCOL2a/GmCOL2b but conservation of their protein function. Dynamic evolution of GmCOL regulatory mechanisms may underlie the evolution of photoperiodic signaling in soybean. Citation: Wu F, Price BW, Haider W, Seufferheld G, Nelson R, et al. (2014) Functional and Evolutionary Characterization of the CONSTANS Gene Family in Short-Day Photoperiodic Flowering in Soybean. PLoS ONE 9(1): e85754. doi:10.1371/journal.pone.0085754 Editor: David E. Somers, Ohio State University, United States of America Received July 19, 2013; Accepted November 29, 2013; Published January 21, 2014 This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication. Funding: This project is supported by the Agriculture and Food Research Initiative Competitive Grants Program from the USDA National Institute of Food and Agriculture (USDA-NIFA2011-00078). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist. * E-mail: [email protected] . These authors contributed equally to this work.

BINDING, KELCH REPEAT, F-BOX 1 (FKF1) in the late afternoon under long day (LD). On the CO promoter, the GIFKF1 complex degrades CYCLING DOF FACTOR 1 (CDF1), a repressor of CO mRNA expression, resulting in activation of CO transcription at the end of the day under LD. An example of the ‘‘external coincidental mechanism’’, CO acts as a point of integration of the internal circadian clock and the external daynight cycles. This integration occurs when CO mRNA oscillation coincides with specific light quality at a specific time of the day. In the early morning under LD, the red-light receptor PHYB promotes degradation of CO protein, but the far-red receptor PHYA and the blue-light receptors CRYPTOCHROME 1 (CRY1) and CRY2 stabilize CO protein toward the end of the day. FKF1 is also known to stabilize CO protein in the afternoon by interacting with CO [13]. Other factors, including SUPRESSOR OF PHYA 1 (SPA1), SPA3, SPA4 and CONSTITUTIVE MORPHOGENIC 1 (COP1), regulate CO protein stability during the night. Through the concerted action of light controlling CO mRNA oscillation and protein stability, CO protein becomes highly accumulated only at the end of the day under LD. This high CO accumulation activates the floral inducer FT, a mobile

Introduction Soybean, one of the most important sources of protein and oil in the world, flowers in response to a photoperiod change from long day to short day at the summer solstice [1]. Despite this crop’s agronomic importance, the molecular basis of the photoperiodic flowering response in soybean remains only partially characterized. Several studies have identified eight Maturity Loci, E1 to E8, that affect flowering and seed maturity time [2–6]. The causal genes for four of these maturity loci, E1– E4, were recently identified. The photoreceptor gene PHYTOCHROME A (PHYA) was isolated as the causal gene of E3 and E4 [7,8], and a homolog of the key flowering gene GIGANTEA (GI) was isolated as the causal gene of E2 [9]. A transcription factor containing a plantspecific B3 domain was also recently cloned as the causal gene of E1 [10]. In the long-day flowering model plant Arabidopsis thaliana (Arabidopsis), the nuclear protein CONSTANS (CO) plays a central role in photoperiodic flowering control [11,12]. Expression of CO mRNA is partly regulated by GIGANTEA (GI), a component of the circadian clock and a regulator of photoperiodic flowering. Under blue light, GI forms a complex with FLAVIN-

PLOS ONE | www.plosone.org

1

January 2014 | Volume 9 | Issue 1 | e85754

Soybean CONSTANS Genes in Photoperiodic Flowering

florigen signal synthesized in the leaves and transmitted to the shoot apical meristem to cause flowering transition. Multiple factors regulate FT in a CO-independent manner in photoperiodic flowering. For example, GI regulates the amount of the small non-coding RNA microRNA172 (miR172), which induces FT expression through repression of APETALA2-like (AP2-like) transcription factors independently of CO [11,12]. GI also binds directly to the FT promoter and enhances FT expression. In addition, an Arabidopsis histone deacetylate (HDAC) complex accumulated at dusk under LD is known to modulate histone deacetylation of FT chromatin, leading to FT expression [14]. Moreover, recent studies suggest that various chromatin modifiers are involved in the control of FT expression [15]. The function of CO in flowering induction is conserved among distantly related flowering plants [16–18]. In the short-day plant rice, the HEADING DATE 1 (Hd1) gene encoding a CO ortholog induces expression of FT homologs Hd3a and RFT1 and promotes flowering under short day (SD), analogous to FT induction by CO in Arabidopsis [16]. Similarly, in Pharbitis nil, a Japanese morning glory whose short-day flowering habit exemplifies a classical model to study photoperiodism [19], possible involvement of the Pharbitis CO homolog (PnCO) and FT homologs (PnFT) in short-day flowering induction has been suggested [20,21], although unknown factors other than PnCO may induce PnFT expression [21]. CO homologs in other species, including potato [22], wheat [23], ryegrass [24], grape [25], alfalfa [26] and barley [27,28], are also thought to be involved in photoperiodic flowering induction. Despite the conserved role of CO in flowering, considerable variation in the CO function in flowering has also been reported. Known as a ‘‘dual-functional’’ flowering regulator, rice Hd1 induces flowering under SD but represses flowering under LD. In poplar, ectopically expressed poplar CO1 and CO2, which show highest sequence similarity to Arabidopsis CO among poplar CO homologs, cause changes in metabolic gene expression but no changes in flowering time [29]. Whether poplar recruits other CO homologs or currently unknown genes in flowering remains unclear. Because of limited information about CO function in soybean, expression analysis has been performed to aid functional characterization of soybean CO homologs (GmCOLs) and FT homologs (GmFTs) [30–37]. The function of two of these homologs, GmFT2a and GmFT5a, in flowering control was demonstrated by transgenic approaches in Arabidopsis and soybean [35,36]. These results suggest a conserved role of GmFTs in flowering in soybean; however, further information is required to clarify the function of GmCOLs and the regulatory interaction between GmCOLs and GmFTs. Here we report comprehensive characterization of the soybean CO gene family and its divergence in gene and protein structure, mRNA expression, regulatory sequence and protein function. Expression of GmCOL1a and GmCOL1b, which clustered together with Arabidopsis CO in our phylogenetic analysis, exhibited a strong diurnal rhythm under SD and showed a peak at dawn that overlapped with expression of GmFT5a. In addition, GmCOL1a and GmCOL1b, as well as GmCOL2a and GmCOL2b, complemented the late flowering effect of the co mutant in Arabidopsis. Our results support a hypothesis that GmCOL1a and GmCOL1b are key inducers of flowering, likely through the induction of GmFT5a in the morning under SD. Our data also indicate dynamic evolution of regulatory sequences of the GmCOL family.

PLOS ONE | www.plosone.org

Materials and Methods Phylogenetic Analysis A BLAST (Basic Local Alignment Search Tool) search was carried out using the full-length amino acid sequences of CO and COL homologs from Arabidopsis, rice and Chlamidomonas reinhardtii against the soybean genome from the Phytozome database (www. phytozome.net). Amino acid sequences similar in length to Arabidopsis CO and COLs were chosen and used for BLAST search for additional CO and COL homologs in the soybean genome. BLAST search continued until no more new homologs appeared. In total, 26 soybean CO homologs (GmCOLs) were obtained. Full-length protein sequences of CO homologs from soybean, Arabidopsis, rice and Chlamidomonas reinhardtii were aligned using ClustalW in MEGA5.1 (http://www. megasoftware.net/) with default parameters (pairwise alignment gap opening = 10.0, pairwise alignment gap extension = 0.1, multiple alignment gap opening = 10.0, multiple alignment gap extension = 0.2, and minimum gap separation distance = 4). The phylogenic tree was generated using MEGA5.1 by the NeighborJoining method. The bootstrap analysis employed 2,000 replicates.

Plant Growth Condition and Sampling Seeds of the seven soybean genotypes were provided by the USDA Soybean Germplasm Collection. Two common North American cultivars, Clark (PI 548533) and Williams 82 (PI 518671); the four near isogenic lines (NILs) of E loci carrying contrasting alleles in E1, E2, E3 and E5 (PI 547431; PI 547432; PI 547610; PI 591490); and Glycine soja (PI 549046) were used (Table 1). Plants were grown in the greenhouse under SD (10 hours of light, 6:45–16:45) and LD (16 hours of light, 6:45–22:45) conditions at 25uC and were sampled every four hours at six time points, T1– T6 (6:30, 10:30, 14:30, 18:30, 22:30 and 2:30), over a 24-hour time period three weeks after germination. For a shift experiment, plants were first grown under LD for three weeks and then transferred to SD for 5 days. A whole shoot above the cotyledon, including 3–4 trifoliates, stem and shoot meristems, was harvested from each plant. Three to four biological replications were sampled for each time point and photoperiod condition.

Table 1. Seven soybean inbred lines, their allele types for Maturity Loci (E loci) and the days to flowering.

Variety

Maturity Loci (E loci)

Days to flowering Field

Greenhouse SD

LD

Clark

e1E2E3E4e5E7

33

2761

6362

Williams 82

E2

33

2761

6163

L65-3366

E1E2E3E4e5E7

61

2962

97610

L66-432

e1e2E3E4e5E7

50

2961

7464

L74-441

e1E2e3E4e5E7

56

3061

5263

L92-1195

e1E2e3E4E5E7

36

2860

7069

Glycine soja

N/A

N/A

2761

10768

The E loci genotype information is based on the record in the USDA soybean germplasm collection. The flowering time in the field is the average of flowering time data from two consecutive years based on the record in the USDA soybean germplasm collection. doi:10.1371/journal.pone.0085754.t001

2

January 2014 | Volume 9 | Issue 1 | e85754

Soybean CONSTANS Genes in Photoperiodic Flowering

RNA Preparation

Motif Logo Analysis, Dot Plot Analysis and Hierarchical Clustering Analysis

Samples stored in 280uC were first disrupted in a mortar for homogenization and then ground by the TissueLyser II (QiagenH) while frozen. The fine powder was then used for RNA preparation following the protocol of the RNeasyH Plant Mini Kit (QiagenH) with a few minor modifications. The RNA samples were measured for quality and quantity using the NanoDropTM 1000 Spectrophotometer (Thermo ScientificH), and cDNA libraries were synthesized using SuperScriptTM III First-Strand Synthesis System for RT-PCR kit (InvitrogenTM). The quality of the cDNA was tested by RT-PCR using a pair of primers specific for the housekeeping gene GmPBB2 [30] (Table S1).

For the motif logo diagrams, sequences corresponding to representative motifs were aligned and then submitted to Weblogo [41] for visualization. The dot plot was created using Nucleic Acid Dot Plots program (http://www.vivo.colostate.edu/molkit/ dnadot/index.html) with the following parameters: window size: 11; mismatch limit: 0. Coding sequence and 2.5 kb upstream noncoding region (introns were excluded) were used for generating the diagram. Hierarchical clustering analysis was performed using MeV software [42]. The hierarchical clustering method [43] was used to construct a binary tree by grouping expression of GmCOLs in RPKM values obtained by RNA sequencing.

Quantitative RT-PCR Gene Cloning, Vector Construction and Arabidopsis Transformation

The quantitative RT-PCR (qRT-PCR) was performed using a 7900 HT Fast Real-Time PCR System (Applied BiosystemsH), following the manufacturer’s manual of FastStart Universal SYBR Green Master (Rox) (Roche Applied ScienceH). Gene-specific primers used for qRT-PCR are listed in Table S1. All reactions were carried out in a 384-well Clear Optical Reaction Plate (Applied BiosystemsH), with a volume of 15 ml per well that, consisted of 7.5 ml 26FastStart Master Mix, 6.4 ml sterilized distilled water, 0.3 ml of each primer (10 M), and 0.5 ml template. The resulting data were recorded and analyzed by the 7900 HT RT-PCR System software. Transcript levels were calculated relative to those of the reference gene GmPBB2 (Table S1). Each qRT-PCR reaction was performed in triplicate and all data were presented as means 6 SEM.

DNA fragments containing the coding sequence of soybean COL1a, COL1b, COL2a and COL2b were amplified from cDNA of the soybean Williams 82 using gene-specific primers (Table S1) and then cloned into the pCR2.1 vector (Invitrogen) following manufacturer’s instructions. Several independent clones were sequenced for validation purposes. The coding sequence of the four genes was subcloned into the entry vector pCR8 (Invitrogen) and finally transferred to the overexpression binary vector pEarley100 [44] via the LR Gateway recombination reaction following manufacturer’s instructions (Invitrogen). Arabidopsis plants (Col-0 wild type and co-1 mutant introgressed into Col-0 by outcrossing four times) were infected with Agrobacterium tumefaciens strain pGV3101 transformed with the obtained pEarley100 clones using the floral dipping method [45]. For screening T1 tranformants, T1 seeds obtained from the dipped plants were grown on soil and sprayed with Basta (glufosinate ammonium, Sigma, 100 mg/L) at least three times.

Transcriptome Sequencing Data Analysis 189 RNA samples, 3 photoperiod treatments (SD, LD, and a shift from LD to SD), 3 time points (6:30, 14:30, and 22:30), 7 genotypes (Table 1) and 3 biological replications were used for RNA sequencing. Quality of RNA was tested using electrophoresis on a 1% agarose gel, and 5 ug total RNA for each sample was submitted to the W. M. Keck Center for Comparative and Functional Genomics in the Roy J. Carver Biotechnology Center at the University of Illinois in Urbana-Champaign. The cDNA libraries were prepared using Illumina’s TruSeq RNAseq Sample Prep kit, and RNA sequencing was performed using Illumina sequencing technology with 12 lanes. Samples were assigned randomly to each lane in order to reduce the effect of lanes. To increase the number of reads per sample, we allocated fewer samples to each lane (15 or 16 samples per lane) than the maximum possible number of samples. Obtained reads were 100 nt in length. Number of reads ranged from 150 million to 185 million per lane and 9.5 million to 12.3 million per sample. This will result in approximately 175–228 reads per gene if all genes are equally expressed. Data quality was ensured by FastX (http:// hannonlab.cshl.edu/fastx_toolkit/index.html) and fastQC (http:// www.bioinformatics.bbsrc.ac.uk/projects/download. html#fastqc). The Bowtie package [38] was used to map and align the reads to the soybean transcriptome (Williams 82) obtained from the Phytozome database [39]. The latest 54,175 predicted protein coding loci and their 73,320 transcripts were used. Counts of aligned read were obtained by using an in-house python script. Normalization of mapped reads was performed by RPKM (Reads Per Kilobase per Million mapped reads) [40]. The RNA sequencing data (the accession numbers GSM1234545 – GSM1234733) were deposited in the Gene Expression Omnibus (GEO) site at the National Canter for Biotechnology Information (NCBI) website (http://www.ncbi.nlm.nih.gov/geo/). PLOS ONE | www.plosone.org

Arabidopsis Growth Conditions and Measurement of Flowering Time Arabidopsis plants were grown in growth chamber under LD conditions (16 h light/8 h dark) with white fluorescence light at 22uC. For flowering time determination, the number of rosette leaves was counted at the time the first flowers bloomed [46].

Results Identification of CO Homologs in Soybean A BLAST (Basic Local Alignment Search Tool) search using the full-length amino acid sequences of CO and COL homologs from Arabidopsis, rice and Chlamydomonas reinhardtii identified 26 soybean CO homologs (GmCOLs). Phylogenetic analysis clarified the evolutionary relationships of GmCOLs (Figure 1A). Each GmCOL was then assigned a name based on its homology level to Arabidopsis CO and COLs, with a designation of a or b for homeologous pairs originating from the most recent duplication event (e.g. GmCOL1a and GmCOL1b). The 26 GmCOLs comprised 13 such homeologous pairs that were classified into three clades: 4 pairs in Clade I, three pairs in Clade II, and six pairs in Clade III. Clade I contained two sub-clades that nested two such homeologous pairs, and Clade III contained two subclades that nested three such homeologous pairs. This observation corresponds well with the multiple genome duplication events that have occurred in soybean [39]. Among the eight GmCOLs in Clade I, GmCOL1a, GmCOL1b, GmCOL2a and GmCOL2b clustered together with Arabidopsis CO and rice CO (Hd1), wellcharacterized flowering inducers. 3

January 2014 | Volume 9 | Issue 1 | e85754

Soybean CONSTANS Genes in Photoperiodic Flowering

Figure 1. CO homologs in soybean. (A) Phylogenetic analysis of full-length amino acid sequences of CO homologs of soybean (Gm), Arabidopsis (At), rice (Os) and Chlamydomonas reinhardtii (Cr). Numbers at nodes indicate the value of 1,000 bootstrap analyses. The domain structure of B-box 1 (red rectangles), B-box 2 (blue rectangles), CCT (yellow rectangles) and the introns (black triangle) of the genes are shown. Followed are the accession numbers of the proteins and their corresponding databases (in square brackets). (B) BLAST best hits of CO homologs in Arabidopsis (A. thaliana) and soybean (G. max). The BLAST best hits are shown by dotted lines with an arrow head, and the protein pairs of the reciprocal best hits are shown by lines with double-head arrows. doi:10.1371/journal.pone.0085754.g001

PLOS ONE | www.plosone.org

4

January 2014 | Volume 9 | Issue 1 | e85754

Soybean CONSTANS Genes in Photoperiodic Flowering

CO homologs, which constitute a subgroup of zinc finger proteins, contain one or more zinc-binding B-box domains at the N terminus and a CCT domain at the C terminus [47,48]. CO homologs in each clade possess specific gene structure and protein domains [48,49]. GmCOL homologs in each clade showed structural characteristics similar to each other (Figure 1A). GmCOLs in Clades I and II contained one intron, with the exception of GmCOL5b, which contained three introns. GmCOLs in Clade III contained three introns. GmCOLs in Clades I and III possessed two B-boxes and a CCT domain in the protein structure, whereas GmCOLs in Clade II possessed only one B-box domain and a CCT domain. Comparison of the amino acid sequences of the domains indicated that B-box 2 domain was less conserved compared to B-box 1 and CCT domains between Arabidopsis and soybean (Figure S1). A reciprocal BLAST search then identified potential orthologous relationships in CO homologs between Arabidopsis and soybean (Figure 1B). Supporting the observed clades being conserved among species in the phylogenetic analysis, the BLAST best hits of all Arabidopsis CO homologs were identified from soybean homologs within the same clade and vice versa. The analysis revealed seven reciprocal BLAST best hit pairs: 2 in Clade I, 1 in Clade II and 4 in Clade III. GmCOL2a and GmCOL2b were the BLAST best hits of Arabidopsis CO; however, no reciprocal BLAST best hit pairs existed between Arabidopsis CO and GmCOLs.

points based on the RNA sequencing data under SD and LD (Figure 2). In contrast, striking variation appeared among the genotypes in the shift experiment. GmCOL1a, GmCOL1b and GmCOL3b in Clade I showed elevated mRNA abundance at all three time points in Clark and Williams 82 and at T1 (6:30) in a NIL, L65-3366 (E1E2E3E4e5E7), in response to the photoperiod shift. A very similar response in mRNA abundance to the photoperiod shift was observed in GmCOL5a in Clade II. Another notable variation among genotypes was the elevated mRNA abundance of GmCOL7a, GmCOL7b and many GmCOLs in Clade III at T3 (22:30) in L74-441 (e1E2e3E4e5E7) and at all three time points in L92-1195 (e1E2e3E4E5E7) and G. soja.

Clustering of GmCOL Expression Patterns A hierarchical clustering analysis of the mRNA expression patterns of GmCOLs under LD and SD revealed rapid divergence of regulatory mechanisms among GmCOLs as well as high conservation among recently derived homeologs. Two main branches appeared: a cluster of GmCOLs from Clades I and II, and a cluster of GmCOLs from Clade III (Figure 3), in agreement with the coding sequence similarity noted in the phylogenetic analysis (Figure 1A). An exception was GmCOL7a and GmCOL7b from Clade II clustering together with GmCOL8a/GmCOL8b and GmCOL9a/GmCOL9b from Clade III, creating a sub-clade nested under the cluster of Clade III genes (Figure 3). This observation is consistent with the similarity in their mRNA abundance in response to the photoperiod shift (Figure 2). All pairs of homeologous genes derived from the last genome duplication clustered together, except the GmCOL7a/mCOL7b, GmCOL8a/ GmCOL8b and GmCOL9a/GmCOL9b pairs (Figure 3).

Photoperiods Affect Flowering Time Variation All seven genotypes used in this study flowered earlier under SD than under LD. G. soja flowered extremely late under LD (106.7 days after germination) but as early as other genotypes under SD (27.2 days), exhibiting the strongest response to photoperiods among the genotypes (Table 1). Flowering time among the genotypes showed large variation under LD (51.5 days to 106.7 days), whereas flowering time under SD showed more subtle variation (27.2 days to 30.2 days).

GmCOL1a and GmCOL1b Show a Diurnal Rhythm with a Peak before Dawn under SD GmCOL1a, GmCOL1b, GmCOL2a and GmCOL2b mRNA accumulation was further characterized by qRT-PCR at all six time points. The diurnal expression patterns of GmCOL1a and GmCOL1b observed by qRT-PCR corresponded very well with the RNA sequencing results (Figure 4). However, GmCOL2a and GmCOL2b showed less conserved patterns between RNA sequencing and qRT-PCR, likely due to their extreme low abundance. GmCOL1a and GmCOL1b demonstrated similar expression patterns to each other, with a diurnal rhythm in both SD and LD (Figure 4). Under SD, GmCOL1a and GmCOL1b expression rose in late night (T6: 2:30), peaked at dawn (T1: 6:30) and declined sharply over the next four hours (T2: 10:30). The abundance remained low during the day and after dusk until late night. Under LD, GmCOL1a and GmCOL1b expression remained low during the day and early evening, peaked in the middle of the night (T6: 2:30) and declined before dawn. In general, expression patterns of GmCOL1a and GmCOL1b demonstrated similarity among the seven genotypes under SD in general. Notable variations among the genotypes, however, appeared in LD: higher peak expression of GmCOL1a at T6 (2:30) in Clark and Williams 82 than in other genotypes, and higher expression of GmCOL1b at T6 (2:30) in a near isogenic line (NIL), L74-441 (e1E2e3E4e5E7). GmCOL2a and GmCOL2b exhibited a similar diurnal expression rhythm under SD, but their rhythmic expression patterns were less clear under LD (Figure 4). In clear contrast to GmCOL1a and GmCOL1b, expression of GmCOL2a and GmCOL2b in SD peaked after dusk (T4: 18:30) and declined during the night. Under LD, their expression appeared to peak at two time points: T4 (18:30) and T6 (2:30).

Photoperiodic Regulation of GmCOLs To survey the mRNA abundance of GmCOLs among different genotypes and photoperiods, an RNA sequencing approach was employed with three representative time points: T1 (6:30), T3 (2:30) and T5 (22:30) (Figure 2). A wide range of mRNA abundance was observed among GmCOLs. GmCOL1a, GmCOL1b, GmCOL3b and GmCOL5a showed notably high expression with strong rhythmic patterns particularly under SD with a peak at T1 (6:30), while GmCOL2a, GmCOL2b, GmCOL10a and GmCOL10b displayed faint expression under both SD and LD. Expression patterns of homeologous genes were generally well conserved. All expressed GmCOLs showed different expression patterns between SD and LD, suggesting photoperiodic regulation of their mRNA abundance. Under SD, the peak expression levels of GmCOL1a, GmCOL1b, GmCOL3b and GmCOL5a at T1 (6:30) exceeded 250 RPKM (Reads Per Kilobase of transcript per Million mapped reads), with GmCOL1b showing the highest level of approximately 700 RPKM. Under LD, the expression levels of these genes at T1 (6:30) were more than 2-fold lower than under SD, with GmCOL1b showing the strongest reduction of more than 7-fold. GmCOL2a and GmCOL2b demonstrated much lower expression (below 2.0 RPKM) than GmCOL1a and GmCOL1b under both SD and LD.

Variation in Expression of GmCOLs among Genotypes No significant variation in the mRNA abundance of GmCOLs appeared among the seven genotypes at the representative time PLOS ONE | www.plosone.org

5

January 2014 | Volume 9 | Issue 1 | e85754

Soybean CONSTANS Genes in Photoperiodic Flowering

PLOS ONE | www.plosone.org

6

January 2014 | Volume 9 | Issue 1 | e85754

Soybean CONSTANS Genes in Photoperiodic Flowering

Figure 2. mRNA abundance of GmCOLs measured by RNA sequencing. RPKM values are displayed on the left. The top panel shows GmCOLs in Clade I, the middle panel shows GmCOLs in Clade II, and the bottom panel shows GmCOLs in Clade III. SD is 10 hours light (6:45–16:45), LD is 16 hours light (6:45–22:45), and LD-SD is a shift from three weeks LD to 5 days SD. Samples are three representative time points: T1 (6:30), T3 (14:30) and T5 (22:30). doi:10.1371/journal.pone.0085754.g002

sequencing and qRT-PCR results indicated that GmFT5a expression rose in late night (T6: 2:30), peaked at dawn (T1: 6:30) and declined sharply over the next four hours (T2: 10:30) under SD, whereas no expression was detected under LD (Figure 5). Peak expression of GmFT5a occurred at T2 (10:30) in all genotypes

GmCOL1a and GmCOL1b Coincide with GmFT5a under SD Abundance of GmFT2a and GmFT5a mRNA accumulation peaks at 3–4 hours after dawn under SD [33]. Similarly, our RNA

Figure 3. Hierarchical clustering of GmCOLs based on their mRNA abundance under LD and SD. GmCOL2a, GmCOL2b, GmCOL10a and GmCOL10b are excluded due to their extreme low abundance. To demonstrate the variation in the mRNA abundance, a rainbow color scheme was used in which red indicates high expression and blue indicates low expression. doi:10.1371/journal.pone.0085754.g003

PLOS ONE | www.plosone.org

7

January 2014 | Volume 9 | Issue 1 | e85754

Soybean CONSTANS Genes in Photoperiodic Flowering

Figure 4. mRNA abundance of GmCOL1a, GmCOL1b, GmCOL2a and GmCOL2b. Circles and solid lines represent the data obtained by qRT-PCR, and gray bars represent the data obtained by RNA sequencing. Relative expression levels of qRT-PCR among the samples in each panel are shown on

PLOS ONE | www.plosone.org

8

January 2014 | Volume 9 | Issue 1 | e85754

Soybean CONSTANS Genes in Photoperiodic Flowering

the left, and RPKM values of RNA sequencing are shown on the right. SD is 10 hours light (6:45–16:45) and LD is 16 hours light (6:45–22:45). Time points T1–T6 represent 6:30, 10:30, 14:30, 18:30, 22:30 and 2:30. Note that RNA sequencing samples contain three representative time points: T1, T3 and T5. doi:10.1371/journal.pone.0085754.g004

except L66-432 (e1e2E3E4e5E7), in which GmFT5a peaked at T1 (6:30). The mRNA abundance of GmFT5a corresponds very well with that of GmCOL1a and GmCOL1b, which peaked at dawn under SD preceding the GmFT5a peak (Figures 4 and 5). Like GmFT5a, GmFT2a showed expression under SD but not under LD. However, the expression pattern of GmFT2a under SD appeared different from that of GmFT5a: it peaked at T3 (14:30), declined at T4 (18:30) and increased again at T5 (22:30). The second peak was more apparent in L65-3366 (E1E2E3E4e5E7), L74-441 (e1E2e3E4e5E7), L92-1195 (e1E2e3E4E5E7) and G. soja than in other genotypes. This second peak of GmFT2a seems to overlap with the peak expression of GmCOL2a and GmCOL2b (Figures 4 and 5).

Discussion GmCOL1a and GmCOL1b are Potential Flowering Inducers We identified 26 COL homologs from the soybean genome and characterized their photoperiod response in mRNA accumulation, sequence divergence and protein function. The CO gene family participates in a wide range of events in plant development in response to photoperiodic signaling, including flowering [50], seedling growth [51,52], dormancy [53] and tuberization [54] in vascular plants, as well as cell growth in the green algae Chlamidomonas [18]. Functional divergence of the CO gene family in Arabidopsis has been reported. For instance, despite their high sequence similarity with CO, COL1 and COL2 cause no changes in flowering time when ectopically expressed [55]. COL3 influences root growth and lateral root formation [51] and acts also as a repressor of flowering, unlike CO. Similar to COL3, COL9 represses flowering [56]. We observed that the soybean CO gene family, like its counterpart in Arabidopsis, showed highly diverse expression patterns while responding to photoperiods, strongly suggesting GmCOLs’ functional divergence in multiple aspects of soybean photoperiodic response, including flowering. We gathered evidence indicating that homeologous GmCOL1a and GmCOL1b are inducers of flowering in soybean. A key piece of evidence supporting the role of GmCOL1a and GmCOL1b in flowering is that they complemented the function of Arabidopsis CO in our transgenic study (Figure 7). Further evidence for the role of GmCOL1a and GmCOL1b in flowering stems from their high sequence similarity to Arabidopsis CO (Figure 1A). Our phylogenetic analysis showed that two homeologous pairs of GmCOLs, GmCOL1a/GmCOL1b and GmCOL2a/GmCOL2b, clustered together with Arabidopsis CO and rice Hd1. Among GmCOLs, GmCOL1a and GmCOL1b showed the highest expression, while GmCOL2a and GmCOL2b showed extremely low expression and a weak peak after dusk (Figure 4), indicating the possibility of GmCOL2a and GmCOL2b becoming pseudogenes or evolving novel function. Offering more support for their role in flowering, GmCOL1a and GmCOL1b display strong photoperiod responsiveness in two ways. First, the mRNA abundance level and patterns of GmCOL1a and GmCOL1b showed striking discrepancies between SD and LD, with a high peak under flowering-inductive SD (Figure 4). Second, GmCOL1a and GmCOL1b, as well as GmCOL3a and GmCOL5a, exhibited a strong response to the photoperiod shift (Figure 2). Taken together, our results suggest that GmCOL1a and GmCOL1b are inducers of photoperiodic flowering. Our work may also indicate that photoperiodic flowering regulation by the CO gene family in flowering plants stems from an ancestral gene of Clade I that harbors Arabidopsis CO and rice Hd1 as well as GmCOL1a and GmCOL1b.

Rapid Divergence of GmCOL Regulatory Mechanisms Comparison of the 2.5 kb upstream intergenic region of GmCOL1a, GmCOL1b, GmCOL2a and GmCOL2b showed highly conserved promoter sequences of homeologous genes derived from the last genome duplication (Figure 6). GmCOL1a and GmCOL1b shared high sequence similarity throughout the coding region and the 2.5 kb upstream intergenic region. Similarly, GmCOL2a and GmCOL2b displayed high similarity in the coding region and a restricted upstream intergenic region (0.5 kb). In contrast, comparison between homeologous genes derived from the earlier genome duplication, GmCOL1a/GmCOL1b and GmCOL2a/ GmCOL2b, showed less or no obvious similarity in the coding region and the upstream intergenic region, which is consistent with their diverse mRNA expression patterns. Other GmCOLs demonstrated similar trends in sequence divergence of the upstream intergenic region: high similarity between homeologous gene pairs derived from the last genome duplication and no similarity among homeologous genes derived from the earlier genome duplication (data not shown), due to the longer evolutionary time. Despite the striking similarity observed in the mRNA abundance of GmCOL5a and GmCOL1a/GmCOL1b in response to photoperiods (Figure 3), no sequence similarity was observed in the upstream intergenic region (data not shown). Similarly, no sequence similarity was observed in the upstream intergenic region among GmCOL7a/GmCOL7b, GmCOL8a/GmCOL8b and GmCOL9a/GmCOL9b.

GmCOL1a, GmCOL1b, GmCOL2a and GmCOL2b Share Conserved Protein Function In marked contrast to GmCOL’s rapid evolution of regulatory sequences, their functional evolution is slower. To examine the protein function of GmCOL1a, GmCOL1b, GmCOL2a and GmCOL2b derived from the past genomic duplication events, we expressed the coding sequences of GmCOL1a, GmCOL1b, GmCOL2a and GmCOL2b under the control of the CaMV 35S promoter in Arabidopsis plants carrying a mutation in the CO gene. While the co-1 mutant showed significant late flowering compared to wild-type, the transgenic plants expressing GmCOL1a, GmCOL1b, GmCOL2a and GmCOL2b fully complemented the late flowering effect of co-1 (Figure 7). This observation suggests that the protein function of GmCOL1a, GmCOL1b, GmCOL2a and GmCOL2b was highly conserved, despite the significant regulatory divergence between GmCOL1a/GmCOL1b and GmCOL2a/ GmCOL2b. PLOS ONE | www.plosone.org

GmCOL1a and GmCOL1b are Potential Activators of GmFTs Our expression study showed that the peak expression of GmCOL1a and GmCOL1b at dawn overlapped with the peak expression of GmFT5a in the morning under flowering inductive SD (Figures 4 and 5). This result may indicate that the regulatory interaction between CO and FT reported in Arabidopsis, rice and 9

January 2014 | Volume 9 | Issue 1 | e85754

Soybean CONSTANS Genes in Photoperiodic Flowering

Figure 5. mRNA abundance of GmFT2a and GmFT5a. Circles and solid lines represent the data obtained by qRT-PCR, and gray bars represent the data obtained by RNA sequencing. Relative expression levels of qRT-PCR among the samples in each panel are shown on the left, and RPKM values of RNA sequencing are shown on the right. SD is 10 hours light (6:45–16:45), LD is 16 hours light (6:45–22:45) and LD-SD is a shift from three weeks LD to 5 days SD. Time points T1–T6 represent 6:30, 10:30, 14:30, 18:30, 22:30 and 2:30. Note that RNA sequencing samples contain three representative time points: T1, T3 and T5. doi:10.1371/journal.pone.0085754.g005

other flowering plants is conserved in soybean. However, further clarification is required to confirm this hypothesis. In the SD flowering plants rice and Pharbitis, expression of FT homologs peaks at the end of the night under SD, consistent with SD flowering plants’ requirement of exposure to long nights. Our results and published findings show that GmCOL1a and GmCOL1b

PLOS ONE | www.plosone.org

expression peaks at the end of the night when GmFT5a expression rises (Figures 4 and 5), indicating that a mechanism similar to the one in rice and Pharbitis measures photoperiod to regulate flowering in soybean. In contrast, GmFT2a expression peaks several hours after dawn in the afternoon, unlike the peak expression of GmFT5a in the

10

January 2014 | Volume 9 | Issue 1 | e85754

Soybean CONSTANS Genes in Photoperiodic Flowering

Figure 6. Comparison of the 2.5 kb upstream intergenic region and the coding region of GmCOL1a, GmCOL1b, GmCOL2a and GmCOL2b. Empty arrows indicate the coding sequence (introns were excluded); lines indicate the upstream intergenic region. Conserved sequences between two sequences are evident from the detected diagonal dotted lines. The dot plot was created using the Nucleic Acid Dot Plots program, with the following parameters: window size: 11; mismatch limit: 0. doi:10.1371/journal.pone.0085754.g006

regulatory divergence, as well as protein functional divergence, is the primary evolutionary source of novel gene function [58–60]. Our results show that regulatory sequence divergence among GmCOL1a/GmCOL1b and GmCOL2a/GmCOL2b preceded their protein functional divergence (Figures 4, 6 and 7), supporting the importance of regulatory divergence that may lead to functional evolution of duplicated genes. Providing further insight into regulatory evolution of the GmCOL gene family, the mRNA accumulation of GmCOL5a, which belongs to Clade II, showed a photoperiodic response resembling that of GmCOL1a and GmCOL1b in Clade I (Figure 2). The striking similarity in the photoperiodic regulation of mRNA accumulation among distantly related GmCOL1a/GmCOL1b and GmCOL5a may indicate rewiring of genetic networks through evolution of co-regulation mechanisms among them. Similarly, remarkable resemblance observed in the mRNA abundance of GmCOL7a/GmCOL7b in Clade II and GmCOLs in Clade III may also imply that regulatory evolution leads to re-arrangement of genetic interactions and, hence, evolution of novel function. Although further study is required to clarify their regulatory interaction and function, our work demonstrates dynamic evolution of the soybean COL genes and their regulatory mechanisms that may underlie the evolution of the photoperiodic signaling in soybean.

morning. The time lag between the peak expression of GmCOLs and GmFT2a may suggest that unidentified mechanisms are necessary for GmFT2a induction. Such mechanisms may include soybean GI homologs. In Arabidopsis, GI is known to induce FT expression independently of CO [11,12]. Similarly, the GI homolog LATE1 in pea regulates the pea FT homolog FTL but has only minor effect on the expression of a pea CO homolog [57]. Supporting this hypothesis, we observed that GmGIa, the causal gene of the maturity locus E2, peaked at T3 (14:30), corresponding well with the peak expression of GmFT2a (Figure S2). In addition, GmFT2a expression was considerably low in the NIL L66-432 (e1e2E3E4e5E7) that carried the recessive e2 allele and showed lower accumulation of GmGIa at T3 (14:30) (Figures 5 and S2). Another clear difference between GmFT2a and GmFT5a is that GmFT2a exhibits a second expression peak during the night under SD (Figure 5), which coincides with the peak expression of GmCOL2a and GmCOL2b (Figures 4 and 5). The role of the GmFT2a second peak in flowering regulation and the possible involvement of GmCOL2a and GmCOL2b in GmFT2a induction are currently unclear.

Accelerated Regulatory Evolution of GmCOLs How gene function diversifies following gene duplication is a major question in evolutionary biology. It has been suggested that PLOS ONE | www.plosone.org

11

January 2014 | Volume 9 | Issue 1 | e85754

Soybean CONSTANS Genes in Photoperiodic Flowering

Figure 7. Overexpression of GmCOL1a, GmCOL1b, GmCOL2a and GmCOL2b rescued the late flowering phenotype of Arabidopsis co-1 mutant. (A) The vertical axis indicates the number of T1 transgenic plants, and the horizontal axis indicates the number of rosette leaves. (B) Col-0 (left), co-1 mutant (middle) and co-1 mutant overexpressing GmCOL1a (right). doi:10.1371/journal.pone.0085754.g007

visualization created by Weblogo (Crooks et al., 2004). CO homologs in Clade II do not contain B-box 2. (DOCX)

Supporting Information Figure S1 Comparison of B-box 1, B-box 2 and CCT domains of COL homologs. Amino acid sequences of B-box 1, B-box 2 and CCT domains of COL homologs in Clades I, II and III in Arabidopsis (At) and soybean (Gm) are compared. Conserved amino acids are shown in large characters for

PLOS ONE | www.plosone.org

Figure S2 mRNA abundance of Glyma06g23026/E1, Glyma10g36600/E2/GmGIa, Glyma19g41210/E3/ GmPHYA3 and Glyma20g22160/E4/GmPHYA2 mea-

12

January 2014 | Volume 9 | Issue 1 | e85754

Soybean CONSTANS Genes in Photoperiodic Flowering

sured by RNA sequencing. RPKM values are displayed on the left. SD is 10 hours light (6:45–16:45), LD is 16 hours light (6:45– 22:45), and LD-SD is a shift from three weeks LD to 5 days SD. Samples are three representative time points: T1 (6:30), T3 (14:30) and T5 (22:30). (DOCX)

Acknowledgments

Table S1 Primers used for this work.

Conceived and designed the experiments: YH. Performed the experiments: FW BWP GS. Analyzed the data: FB BWP WH. Contributed reagents/ materials/analysis tools: RN. Wrote the paper: FB YH.

We thank Matthew E. Hudson (University of Illinois) for helpful assistance in the sequencing data analysis, and Daniel Wickland (University of Illinois) for critical reading of the manuscript.

Author Contributions

(XLSX)

References 23. Nemoto Y, Kisaka M, Fuse T, Yano M, Ogihara Y (2003) Characterization and functional analysis of three wheat genes with homology to the CONSTANS flowering time gene in transgenic rice. The Plant Journal 36: 82–93. 24. Martin J, Storgaard M, Andersen CH, Nielsen KK (2004) Photoperiodic regulation of flowering in perennial ryegrass involving a CONSTANS-like homolog. Plant molecular biology 56: 159–169. 25. Almada R, Cabrera N, Casaretto JA, Ruiz-Lara S, Villanueva EG (2009) VvCO and VvCOL1, two CONSTANS homologous genes, are regulated during flower induction and dormancy in grapevine buds. Plant cell reports 28: 1193–1203. 26. Herrmann D, Barre P, Santoni S, Julier B (2010) Association of a CONSTANSLIKE gene to flowering and height in autotetraploid alfalfa. Theoretical and applied genetics 121: 865–876. 27. Kikuchi R, Kawahigashi H, Oshima M, Ando T, Handa H (2012) The differential expression of HvCO9, a member of the CONSTANS-like gene family, contributes to the control of flowering under short-day conditions in barley. Journal of experimental botany 63: 773–784. 28. Campoli C, Drosse B, Searle I, Coupland G, von Korff M (2012) Functional characterisation of HvCO1, the barley (Hordeum vulgare) flowering time ortholog of CONSTANS. The Plant Journal 69: 868–880. 29. Hsu C-Y, Adams JP, No K, Liang H, Meilan R, et al. (2012) Overexpression of CONSTANS homologs CO1 and CO2 fails to alter normal reproductive onset and fall bud set in woody perennial poplar. PLoS One 7: e45448. 30. Thakare D, Kumudini S, Dinkins RD (2010) Expression of flowering-time genes in soybean E1 near-isogenic lines under short and long day conditions. Planta 231: 951–963. 31. Liu L, Ma J, Han Y, Chen X, Fu YF (2011) The isolation and analysis of a soybean CO Homologue GmCOL10. Russian Journal of Plant Physiology 58: 330–336. 32. Jiang Y, Han Y, Zhang X (2011) Expression profiles of a CONSTANS homologue GmCOL11 in Glycine max. Russian Journal of Plant Physiology 58: 928–935. 33. Zhang QZ, Ma JH, Chen XJ, Fu YF (2010) Cloning and Analysis of GmCOL4 Gene in Glycine max. Acta Agronomica Sinica 36: 539–548. 34. Ma J, Zhang Q, Chen J, Fu Y (2009) Cloning and expression analysis of GmCOL8 gene in soyabean. Plant Physiology Communications 46: 17–23. 35. Kong F, Liu B, Xia Z, Sato S, Kim BM, et al. (2010) Two coordinately regulated homologs of FLOWERING LOCUS T are involved in the control of photoperiodic flowering in soybean. Plant physiology 154: 1220–1231. 36. Sun H, Jia Z, Cao D, Jiang B, Wu C, et al. (2011) GmFT2a, a soybean homolog of FLOWERING LOCUS T, is involved in flowering transition and maintenance. PLoS One 6: e29238. 37. Jung C-H, Wong CE, Singh MB, Bhalla PL (2012) Comparative genomic analysis of soybean flowering genes. PLoS One 7: e38250. 38. Langmead B, Trapnell C, Pop M, Salzberg SL (2009) Ultrafast and memoryefficient alignment of short DNA sequences to the human genome. Genome Biol 10: R25. 39. Schmutz J, Cannon SB, Schlueter J, Ma J, Mitros T, et al. (2010) Genome sequence of the palaeopolyploid soybean. Nature 463: 178–183. 40. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature methods 5: 621– 628. 41. Crooks GE, Hon G, Chandonia J-M, Brenner SE (2004) WebLogo: a sequence logo generator. Genome research 14: 1188–1190. 42. Saeed AI, Bhagabati NK, Braisted JC, Liang W, Sharov V, et al. (2006) [9] TM4 Microarray Software Suite. Methods in enzymology 411: 134–193. 43. Eisen MB, Spellman PT, Brown PO, Botstein D (1998) Cluster analysis and display of genome-wide expression patterns. Proceedings of the National Academy of Sciences 95: 14863–14868. 44. Earley KW, Haag JR, Pontes O, Opper K, Juehne T, et al. (2006) Gatewaycompatible vectors for plant functional genomics and proteomics. The Plant Journal 45: 616–629. 45. Clough SJ, Bent AF (1998) Floral dip: a simplified method for Agrobacteriummediated transformation of Arabidopsis thaliana. The Plant Journal 16: 735– 743. 46. Hanano S, Goto K (2011) Arabidopsis TERMINAL FLOWER1 is involved in the regulation of flowering time and inflorescence development through transcriptional repression. The Plant Cell Online 23: 3172–3184.

1. Borthwick H, Parker M (1939) Photoperiodic responses of several varieties of soybeans. Botanical Gazette: 341–365. 2. Cober E, Tanner J, Voldeng H (1996) Genetic control of photoperiod response in early-maturing, near-isogenic soybean lines. Crop science 36: 601–605. 3. Abe J, Xu D, Suzuki Y, Kanazawa A, Shimamoto Y (2003) Soybean germplasm pools in Asia revealed by nuclear SSRs. TAG Theoretical and Applied Genetics 106: 445–453. 4. Molnar SJ, Rai S, Charette M, Cober ER (2003) Simple sequence repeat (SSR) markers linked to E1, E3, E4, and E7 maturity genes in soybean. Genome 46: 1024–1036. 5. Zhang WK, Wang YJ, Luo GZ, Zhang JS, He CY, et al. (2004) QTL mapping of ten agronomic traits on the soybean (Glycine max L. Merr.) genetic map and their association with EST markers. TAG Theoretical and Applied Genetics 108: 1131–1139. 6. Liu B, Abe J (2010) QTL mapping for photoperiod insensitivity of a Japanese soybean landrace Sakamotowase. Journal of Heredity 101: 251–256. 7. Liu B, Kanazawa A, Matsumura H, Takahashi R, Harada K, et al. (2008) Genetic redundancy in soybean photoresponses associated with duplication of the phytochrome A gene. Genetics 180: 995–1007. 8. Watanabe S, Hideshima R, Xia Z, Tsubokura Y, Sato S, et al. (2009) Mapbased cloning of the gene associated with the soybean maturity locus E3. Genetics 182: 1251–1262. 9. Watanabe S, Xia Z, Hideshima R, Tsubokura Y, Sato S, et al. (2011) A mapbased cloning strategy employing a residual heterozygous line reveals that the GIGANTEA gene is involved in soybean maturity and flowering. Genetics 188: 395–407. 10. Xia Z, Watanabe S, Yamada T, Tsubokura Y, Nakashima H, et al. (2012) Positional cloning and characterization reveal the molecular basis for soybean maturity locus E1 that regulates photoperiodic flowering. Proceedings of the National Academy of Sciences 109: E2155–E2164. 11. Turck F, Fornara F, Coupland G (2008) Regulation and identity of florigen: FLOWERING LOCUS T moves center stage. Annu Rev Plant Biol 59: 573– 594. 12. Wu F, Hanzawa Y (2014) Photoperiodic control of flowering in plants. In: Pessarakli M, editor. Handbook of plant and crop physiology. 3rd ed. Boca Raton, Florida: CRC press. 129–139. 13. Song SH, Smith RW, To BJ, Millar AJ, Imaizumi T (2012) FKF1 Conveys Timing Information for CONSTANS Stabilization in Photoperiodic Flowering. Science 336: 1045. 14. Gu X, Wang Y, He Y (2013) Photoperiodic Regulation of Flowering Time through Periodic Histone Deacetylation of the Florigen Gene FT. PLoS biology 11: e1001649. 15. He Y (2012) Chromatin Regulation of Flowering. Trends Plant Sci 17: 556–562. 16. Yano M, Katayose Y, Ashikari M, Yamanouchi U, Monna L, et al. (2000) Hd1, a major photoperiod sensitivity quantitative trait locus in rice, is closely related to the Arabidopsis flowering time gene CONSTANS. The Plant Cell Online 12: 2473–2483. 17. Bohlenius H, Huang T, Charbonnel-Campaa L, Brunner AM, Jansson S, et al. (2006) CO/FT regulatory module controls timing of flowering and seasonal growth cessation in trees. Science Signalling 312: 1040. 18. Serrano G, Herrera-Palau R, Romero JM, Serrano A, Coupland G, et al. (2009) Chlamydomonas CONSTANS and the Evolution of Plant Photoperiodic Signaling. Current Biology 19: 359–368. 19. Imamura S, Gakkai NSS (1967) Physiology of Flowering in Pharbitis nil: Japanese Society of Plant Physiologists. 20. Liu J, Yu J, McIntosh L, Kende H, Zeevaart JA (2001) Isolation of a CONSTANS ortholog from Pharbitis nil and its role in flowering. Plant physiology 125: 1821–1830. 21. Hayama R, Agashe B, Luley E, King R, Coupland G (2007) A circadian rhythm set by dusk determines the expression of FT homologs and the short-day photoperiodic flowering response in Pharbitis. The Plant Cell Online 19: 2988– 3000. 22. Martı´nez-Garcı´a JF, Virgo´s-Soler A, Prat S (2002) Control of photoperiodregulated tuberization in potato by the Arabidopsis flowering-time gene CONSTANS. Proceedings of the National Academy of Sciences 99: 15211– 15216.

PLOS ONE | www.plosone.org

13

January 2014 | Volume 9 | Issue 1 | e85754

Soybean CONSTANS Genes in Photoperiodic Flowering

54. Gonza´lez-Schain N, Sua´rez-Lo´pez P (2008) CONSTANS delays flowering and affects tuber yield in potato. Biologia Plantarum 52: 251–258. 55. Ledger S, Strayer C, Ashton F, Kay SA, Putterill J (2001) Analysis of the function of two circadian-regulated CONSTANS-LIKE genes. The Plant Journal 26: 15–22. 56. Cheng X-F, Wang Z-Y (2005) Overexpression of COL9, a CONSTANS-LIKE gene, delays flowering by reducing expression of CO and FT in Arabidopsis thaliana. The Plant Journal 43: 758–768. 57. Hecht V, Knowles CL, Vander Schoor JK, Liew LC, Jones SE, et al. (2007) Pea LATE BLOOMER1 is a GIGANTEA ortholog with roles in photoperiodic flowering, deetiolation, and transcriptional regulation of circadian clock gene homologs. Plant physiology 144: 648–661. 58. Adams KL, Wendel JF (2005) Polyploidy and genome evolution in plants. Current opinion in plant biology 8: 135–141. 59. Duarte JM, Cui L, Wall PK, Zhang Q, Zhang X, et al. (2006) Expression pattern shifts following duplication indicative of subfunctionalization and neofunctionalization in regulatory genes of Arabidopsis. Molecular Biology and Evolution 23: 469–478. 60. Tirosh I, Barkai N (2007) Comparative analysis indicates regulatory neofunctionalization of yeast duplicates. Genome biology 8: R50.

47. Klug A, Schwabe J (1995) Protein motifs 5. Zinc fingers. The FASEB journal 9: 597–604. 48. Khanna R, Kronmiller B, Maszle DR, Coupland G, Holm M, et al. (2009) The Arabidopsis B-box zinc finger family. Plant Cell 21: 3416–3420. 49. Griffiths S, Dunford RP, Coupland G, Laurie DA (2003) The evolution of CONSTANS-like gene families in barley, rice, and Arabidopsis. Plant physiology 131: 1855–1867. 50. Valverde F (2011) CONSTANS and the evolutionary origin of photoperiodic timing of flowering. Journal of experimental botany 62: 2453–2463. 51. Datta S, Hettiarachchi G, Deng X-W, Holm M (2006) Arabidopsis CONSTANS-LIKE3 is a positive regulator of red light signaling and root growth. The Plant Cell Online 18: 70–84. 52. Datta S, Hettiarachchi C, Johansson H, Holm M (2007) SALT TOLERANCE HOMOLOG2, a B-box protein in Arabidopsis that activates transcription and positively regulates light-mediated development. The Plant Cell Online 19: 3242–3255. 53. Bohlenius H, Huang T, Charbonnel-Campaa L, Brunner AM, Jansson S, et al. (2006) CO/FT regulatory module controls timing of flowering and seasonal growth cessation in trees. Science Signaling 312: 1040.

PLOS ONE | www.plosone.org

14

January 2014 | Volume 9 | Issue 1 | e85754

Functional and evolutionary characterization of the CONSTANS gene family in short-day photoperiodic flowering in soybean.

CONSTANS (CO) plays a central role in photoperiodic flowering control of plants. However, much remains unknown about the function of the CO gene famil...
1MB Sizes 0 Downloads 0 Views