Arabidopsis pumila, an annual ephemeral plant, plays important roles in preventing wind and sand erosion, water and soil conservation, and microhabitat improvement in the North of Xinjiang, China. Studies of adaptive mechanisms in harsh desert environments at the genetic and genomic levels can be used to more effectively develop and protect this species. The quantitative real-time polymerase chain reaction (qRT-PCR) method is one of the essential means to achieve these goals, and the selection of an appropriate reference gene is the prerequisite for qRT-PCR. In this study, 10 candidate reference genes were identified from the full-length transcriptome data of A. pumila, and their expression stabilities under four abiotic stresses (drought, heat, cold and salt) and in seven different tissues (roots, hypocotyl, cotyledon, leaves, stems, flowers and siliques) were evaluated with four programmes geNorm, NormFinder, Bestkeeper and RefFinder. Although the most stable reference genes were variable under different treatments using different software, comprehensive ranking revealed that UEP and HAF1 under drought stress, UBQ9 and GAPDH under heat stress, UBC35 and GAPDH under cold stress, GAPDH and ACT1 under salt stress, and ACT1 and GAPDH in different tissues were the most stable reference genes. Moreover, GAPDH and UBQ9 were the most suitable reference gene combinations for all samples. The expression pattern of the K⁺ uptake permease gene KUP9 further validated that the selected reference genes were suitable for normalization of gene expression. The identification of reliable reference genes guarantees more accurate qRT-PCR quantification for A. pumila and facilitates functional genomics studies of ephemeral plants.

Arabidopsis pumila is native to the desert region of northwest China and it is extraordinarily well adapted to the local semi-desert saline soil, thus providing a candidate plant system for environmental adaptation and salt-tolerance gene mining. However, understanding of the salt-adaptation mechanism of this species is limited because of genomic sequences scarcity. In the present study, the transcriptome profiles of A. pumila leaf tissues treated with 250 mM NaCl for 0, 0.5, 3, 6, 12, 24 and 48 h were analyzed using a combination of second-generation sequencing (SGS) and third-generation single-molecule real-time (SMRT) sequencing. Correction of SMRT long reads by SGS short reads resulted in 59,328 transcripts. We found 8075 differentially expressed genes (DEGs) between salt-stressed tissues and controls, of which 483 were transcription factors and 1157 were transport proteins. Most DEGs were activated within 6 h of salt stress and their expression stabilized after 48 h; the number of DEGs was greatest within 12 h of salt stress. Gene annotation and functional analyses revealed that expression of genes associated with the osmotic and ionic phases rapidly and coordinately changed during the continuous salt stress in this species, and salt stress-related categories were highly enriched among these DEGs, including oxidation–reduction, transmembrane transport, transcription factor activity and ion channel activity. Orphan, MYB, HB, bHLH, C3H, PHD, bZIP, ARF and NAC TFs were most enriched in DEGs; ABCB1, CLC-A, CPK30, KEA2, KUP9, NHX1, SOS1, VHA-A and VP1 TPs were extensively up-regulated in salt-stressed samples, suggesting that they play important roles in slat tolerance. Importantly, further experimental studies identified a mitogen-activated protein kinase (MAPK) gene MAPKKK18 as continuously up-regulated throughout salt stress, suggesting its crucial role in salt tolerance. The expression patterns of the salt-responsive 24 genes resulted from quantitative real-time PCR were basically consistent with their transcript abundance changes identified by RNA-Seq. The full-length transcripts generated in this study provide a more accurate depiction of gene transcription of A. pumila. We identified potential genes involved in salt tolerance of A. pumila. These data present a genetic resource and facilitate better understanding of salt-adaptation mechanism for ephemeral plants.

FLOWERING LOCUS T (FT) encodes a mobile signal protein, recognized as major component of florigen, which has a central position in regulating flowering, and also plays important roles in various physiological aspects. A mode is recently emerging for the balance of indeterminate and determinate growth, which is controlled by the ratio of FT-like and TERMINAL FLOWER 1 (TFL1)-like gene activities, and has a strong influence on the floral transition and plant architecture. Orthologs of GhFT1 was previously isolated and characterized from Gossypium hirsutum. We demonstrated that ectopic overexpression of GhFT1 in tobacco, other than promoting flowering, promoted lateral shoot outgrowth at the base, induced more axillary bud at the axillae of rosette leaves, altered leaf morphology, increased chlorophyll content, had higher photosynthesis and caused flowers abscission. Analysis of gene expression suggested that flower identity genes were significantly upregulated in transgenic plants. Further analysis of tobacco FT paralogs indicated that NtFT4, acting as flower inducer, was upregulated, whereas NtFT2 and NtFT3 as flower inhibitors were upregulated in transgenic plants under long-day conditions, but downregulated under short-day conditions. Our data suggested that sufficient level of foreign FT might disturb the balance of the endogenous FT paralogs of inducers and repressors and resulted in altered phenotype in transgenic tobacco, emphasizing the expanding roles of FT in regulating shoot architecture by advancing determine growth. Manipulating the ratio for indeterminate and determinate growth factors throughout FT-like and TFL1-like gene activity holds promise to improve plant architecture and enhance crop yield.

Plant phosphatidylethanolamine-binding protein (PEBP) is comprised of three clades: FLOWERING LOCUS T (FT), TERMINAL FLOWER1 (TFL1) and MOTHER OF FT AND TFL1 (MFT). FT/TFL1-like clades regulate identities of the determinate and indeterminate meristems, and ultimately affect flowering time and plant architecture. MFT is generally considered to be the ancestor of FT/TFL1, but its function is not well understood. Here, two MFT homoeologous gene pairs in Gossypium hirsutum, GhMFT1-A/D and GhMFT2-A/D, were identified by genome-wide identification of MFT-like genes. Detailed expression analysis revealed that GhMFT1 and GhMFT2 homoeologous genes were predominately expressed in ovules, and their expression increased remarkably during ovule development but decreased quickly during seed germination. Expressions of GhMFT1 and GhMFT2 homoeologous genes in germinating seeds were upregulated in response to abscisic acid (ABA), and their expressions also responded to gibberellin (GA). In addition, ectopic overexpression of GhMFT1 and GhMFT2 in Arabidopsis inhibited seed germination at the early stage. Gene transcription analysis showed that ABA metabolism genes ABA-INSENSITIVE3 (ABI3) and ABI5, GA signal transduction pathway genes REPRESSOR OF ga1-3 (RGA) and RGA-LIKE2 (RGL2) were all upregulated in the 35S:GhMFT1 and 35S:GhMFT2 transgenic Arabidopsis seeds. GhMFT1 and GhMFT2 localize in the cytoplasm and nucleus, and both interact with a cotton bZIP transcription factor GhFD, suggesting that both of GhMFT1, 2 have similar intracellular regulation mechanisms. Taken together, the results suggest that GhMFT1 and GhMFT2 may act redundantly and differentially in the regulation of seed germination.

The CONSTANS/FLOWERING LOCUS T (CO/FT) regulon plays a central role in the control of flowering time in photoperiod-sensitive plants. Flowering time in wild cotton (Gossypium spp.) has strict photoperiod sensitivity, but domesticated cotton is day-neutral. Information on the molecular characterization of the CO and CO-like (COL) genes in cotton is very limited. In this study, we identified 42 COL homologs (GhCOLs) in the G. hirsutum genome, and many of them were previously unreported. We studied their chromosome distribution, phylogenetic relationships, and structures of genes and proteins. Our results showed that GhCOLs were classified into three groups, and 14 COLs in group I showed conserved structure when compared with other plants. Two homoeologous pairs, GhCOL1-A and GhCOL1-D in Group I, showed the highest sequence similarity to Arabidopsis thaliana CO and rice CO homologous gene Heading date1 (Hd1). Tissue-specific expression showed that 42 GhCOL genes may function as tissue-specific regulators in different cells or organs. We cloned and sequenced the 14 GhCOL genes in Group I related to flowering induction to study their diurnal expression pattern, and found that their expression showed distinct circadian regulation. Most of them peaked at dawn and decreased rapidly to their minima at dusk, then started to accumulate until following dawn under long- or short-day conditions. Transgenic study in the Arabidopsis co-2 mutant demonstrated that GhCOL1-A and GhCOL1-D fully rescued the late-flowering phenotype, whereas GhCOL3-A, GhCOL3-D, GhCOL7-A, and GhCOL7-D partially rescued the late-flowering phenotype, and the other five homoeologous pairs in Group I did not promote flowering. These results indicate that GhCOL1-A and GhCOL1-D were potential flowering inducers, and are candidate genes for research in flowering regulation in cotton.