Background: Real-time quantitative reverse transcription PCR (RT-qPCR) data needs to be normalized for its proper interpretation. Housekeeping genes are routinely employed for this purpose, but their expression level cannot be assumed to remain constant under all possible experimental conditions. Thus, a systematic validation of reference genes is required to ensure proper normalization. For soybean, only a small number of validated reference genes are available to date. Results: A systematic comparison of 14 potential reference genes for soybean is presented. These included seven commonly used (ACT2, ACT11, TUB4, TUA5, CYP, UBQ10, EF1b) and seven new candidates (SKIP16, MTP, PEPKR1, HDC, TIP41, UKN1, UKN2). Expression stability was examined by RT-qPCR across 116 biological samples, representing tissues at various developmental stages, varied photoperiodic treatments, and a range of soybean cultivars. Expression of all 14 genes was variable to some extent, but that of SKIP16, UKN1 and UKN2 was overall the most stable. A combination of ACT11, UKN1 and UKN2 would be appropriate as a reference panel for normalizing gene expression data among different tissues, whereas the combination SKIP16, UKN1 and MTP was most suitable for developmental stages. ACT11, TUA5 and TIP41 were the most stably expressed when the photoperiod was altered, and TIP41, UKN1 and UKN2 when the light quality was changed. For six different cultivars in long day (LD) and short day (SD), their expression stability did not vary significantly with ACT11, UKN2 and TUB4 being the most stable genes. The relative gene expression level of GmFTL3, an ortholog of Arabidopsis FT (FLOWERING LOCUS T) was detected to validate the reference genes selected in this study. Conclusion: None of the candidate reference genes was uniformly expressed across all experimental conditions, and the most suitable reference genes are conditional-, tissue-specific-, developmental-, and cultivar-dependent. Most of the new reference genes performed better than the conventional housekeeping genes. These results should guide the selection of reference genes for gene expression studies in soybean.

As photosynthetic organisms, microalgae can efficiently convert solar energy into biomass. Microalgae are currently used as an important source of valuable natural biologically active molecules, such as carotenoids, chlorophyll, long-chain polyunsaturated fatty acids, phycobiliproteins, carotenoids and enzymes. Significant advances have been achieved in microalgae biotechnology over the last decade, and the use of microalgae as bioreactors for expressing recombinant proteins is receiving increased interest. Compared with the bioreactor systems that are currently in use, microalgae may be an attractive alternative for the production of pharmaceuticals, recombinant proteins and other valuable products. Products synthesized via the genetic engineering of microalgae include vaccines, antibodies, enzymes, blood-clotting factors, immune regulators, growth factors, hormones, and other valuable products, such as the anticancer agent Taxol. In this paper, we briefly compare the currently used bioreactor systems, summarize the progress in genetic engineering of microalgae, and discuss the potential for microalgae as bioreactors to produce pharmaceuticals.

Photoperiodic control of flowering time is believed to affect latitudinal distribution of plants. The blue light receptor CRY2 regulates photoperiodic flowering in the experimental model plant Arabidopsis thaliana . However, it is unclear whether genetic variations affecting cryptochrome activity or expression is broadly associated with latitudinal distribution of plants. We report here an investigation of the function and expression of two cryptochromes in soybean, GmCRY1a and GmCRY2a. Soybean is a short-day (SD) crop commonly cultivated according to the photoperiodic sensitivity of cultivars. Both cultivated soybean ( Glycine max ) and its wild relative ( G. soja ) exhibit a strong latitudinal cline in photoperiodic flowering. Similar to their Arabidopsis counterparts, both GmCRY1a and GmCRY2a affected blue light inhibition of cell elongation, but only GmCRY2a underwent blue light- and 26S proteasome-dependent degradation. However, in contrast to Arabidopsis cryptochromes, soybean GmCRY1a, but not GmCRY2a, exhibited a strong activity promoting floral initiation, and the level of protein expression of GmCRY1a, but not GmCRY2a, oscillated with a circadian rhythm that has different phase characteristics in different photoperiods. Consistent with the hypothesis that GmCRY1a is a major regulator of photoperiodic flowering in soybean, the photoperiod-dependent circadian rhythmic expression of the GmCRY1a protein correlates with photoperiodic flowering and latitudinal distribution of soybean cultivars. We propose that genes affecting protein expression of the GmCRY1a protein play an important role in determining latitudinal distribution of soybeans.

Background: CO and FT orthologs, belonging to the BBX and PEBP family, respectively, have important and conserved roles in the photoperiod regulation of flowering time in plants. Soybean genome experienced at least three rounds of whole genome duplications (WGDs), which resulted in multiple copies of about 75% of genes. Subsequent subfunctionalization is the main fate for paralogous gene pairs during the evolutionary process. Results: The phylogenic relationships revealed that CO orthologs were widespread in the plant kingdom while FT orthologs were present only in angiosperms. Twenty-eight CO homologous genes and twenty-four FT homologous genes were gained in the soybean genome. Based on the collinear relationship, the soybean ancestral CO ortholog experienced three WGD events, but only two paralogous gene pairs (GmCOL1/2 and GmCOL5/13) survived in the modern soybean. The paralogous gene pairs, GmCOL1/2 or GmCOL5/13, showed similar expression patterns in pair but different between pairs, indicating that they functionally diverged. GmFTL1 to 7 were derived from the same ancestor prior to the whole genome triplication (WGT) event, and after the Legume WGD event the ancestor diverged into two branches, GmFTL3/5/7 and GmFTL1/2/4/6. GmFTL7 were truncated in the N-terminus compared to other FT-lineage genes, but ubiquitously expressed. Expressions of GmFTL1 to 6 were higher in leaves at the flowering stage than that at the seedling stage. GmFTL3 was expressed at the highest level in all tissues except roots at the seedling stage, and its circadian pattern was different from the other five ones. The transcript of GmFTL6 was highly accumulated in seedling roots. The circadian rhythms of GmCOL5/13 and GmFT1/2/4/5/6 were synchronized in a day, demonstrating the complicate relationship of CO-FT regulons in soybean leaves. Over-expression of GmCOL2 did not rescue the flowering phenotype of the Arabidopsis co mutant. However, ectopic expression of GmCOL5 did rescue the co mutant phenotype. All GmFTL1 to 6 showed flower-promoting activities in Arabidopsis. Conclusions: After three recent rounds of whole genome duplications in the soybean, the paralogous genes of CO-FT regulons showed subfunctionalization through expression divergence. Then, only GmCOL5/13 kept flowering-promoting activities, while GmFTL1 to 6 contributed to flowering control. Additionally, GmCOL5/13 and GmFT1/2/3/4/5/6 showed similar circadian expression profiles. Therefore, our results suggested that GmCOL5/13 and GmFT1/2/3/4/5/6 formed the complicate CO-FT regulons in the photoperiod regulation of flowering time in soybean.

Background: The Phosphate transporter 1 (PHT1) gene family has crucial roles in phosphate uptake, translocation, remobilization, and optimization of metabolic processes using of Pi. Gene duplications expand the size of gene families, and subfunctionalization of paralog gene pairs is a predominant tendency after gene duplications. To date, experimental evidence for the evolutionary relationships among different paralog gene pairs of a given gene family in soybean is limited. Results: All potential Phosphate transporter 1 genes in Glycine max L. (GmPHT1) were systematically analyzed using both bioinformatics and experimentation. The soybean PHT1 genes originated from four distinct ancestors prior to the Gamma WGT and formed 7 paralog gene pairs and a singleton gene. Six of the paralog gene pairs underwent subfunctionalization, and while GmPHT1;4 paralog gene experienced pseudogenization. Examination of long-term evolutionary changes, six GmPHT1 paralog gene pairs diverged at multiple levels, in aspects of spatio-temporal expression patterns and/or quanta, phosphates affinity properties, subcellular localization, and responses to phosphorus stress. Conclusions: These characterized divergences occurred in tissue- and/or development-specific modes, or conditional modes. Moreover, they have synergistically shaped the evolutionary rate of GmPHT1 family, as well as maintained phosphorus homeostasis at cells and in the whole plant.

Tissue-specific differentiation of phosphorus transporters coordinated with their regulators confers phosphorus homeostasis, regulating nodule development.