Molecular population genetics and gene expression analysis of duplicated CBF genes of Arabidopsis thaliana

Genome-wide identification of CBF genes and their responses to cold acclimation in Taraxacum kok-saghyz

C-repeat binding factors (CBFs) are transcription factors that are known to play important roles in plant cold acclimation. They are highly conserved in most higher plants. Taraxacum kok-saghyz (TKS) is an herb native to China and Kazakhstan and is well-known for its production of rubber silk with industrial and economic value. To understand cold acclimation mechanisms, we conducted a genome-wide discovery of the CBF family genes in TKS and revealed ten CBF genes. A bioinformatic analysis of the CBF genes was carried out to analyze the phylogenetic relationship, protein conservative motifs, protein physicochemical properties, gene structure, promoter cis-acting elements, and the gene expression patterns under cold acclimation and control conditions. It was found that most of these genes were highly responsive at the late stage of cold acclimation, indicating that they play important roles in the cold acclimation processes of TKS. This study provides a theoretical basis for the study of the molecular functions of the CBF gene family in TKS, and a useful guidance for the genetic improvement of the cold tolerance traits of TKS and other plants, including crops.

Genome-wide identification and expression analysis of the CBF/DREB1 gene family in lettuce

The C-repeat binding factor (CBF)/dehydration-responsive element binding (DREB1) proteins play a prominent role in freezing tolerance and are highly conserved in higher plants. Here we performed a genome-wide search of the CBF/DREB1 gene family in lettuce (Lactuca sativa L.) and identified 14 members of the family with one member gene containing a non-sense mutation within the AP2 DNA-binding domain. A comprehensive phylogenetic analysis of the CBF/DREB1 family members in 20 plant species from the Asterid or Rosid clade provided evidence that tandem duplication played an important role in the expansion of the CBF/DREB1 family. Expression analysis showed that twelve of the lettuce CBF genes were responsive to low temperature (4 °C), and that three and six of them could also be responsive to salt and heat stresses, respectively. Unlike Arabidopsis thaliana whose members of the CBF/DREB1 family respond only to a particular stress, lettuce CBFs provide wider protection from combinations of abiotic stresses. A global transcriptome analysis revealed distinctive temporal expression patterns among the cold-regulated genes in lettuce plants exposed to low temperature. Genes induced throughout the cold treatment are enriched in functions associated with protection from UV and high-light intensity and the genes suppressed after 7 days of cold exposure are enriched in photosynthesis-associated functions. These results provide insight into the molecular evolutionary properties of the CBF/DREB1 gene family in lettuce and a reference for genetic improvement of the lettuce response to cold acclimation.

Transcriptome sequencing analysis of alfalfa reveals CBF genes potentially playing important roles in response to freezing stress

Alfalfa (Medicago sativa L.) is an important perennial forage, with high nutritional value, which is widely grown in the world. Because of low freezing tolerance, its distribution and production are threatened and limited by winter weather. To understand the complex regulation mechanisms of freezing tolerance in alfalfa, we performed transcriptome sequencing analysis under cold (4 °C) and freezing (-8 °C) stresses. More than 66 million reads were generated, and we identified 5767 transcripts differentially expressed in response to cold and/or freezing stresses. These results showed that these genes were mainly classified as response to stress, transcription regulation, hormone signaling pathway, antioxidant, nodule morphogenesis, etc., implying their important roles in response to cold and freezing stresses. Furthermore, nine CBF transcripts differentially expressed were homologous to CBF genes of Mt-FTQTL6 site, conferring freezing tolerance in M. truncatula, which indicated that a genetic mechanism controlling freezing tolerance was conservative between M. truncatula and M. sativa. In summary, this transcriptome dataset highlighted the gene regulation response to cold and/or freezing stresses in alfalfa, which provides a valuable resource for future identification and functional analysis of candidate genes in determining freezing tolerance.

The precise regulation of different COR genes by individual CBF transcription factors in Arabidopsis thaliana

The transcription factors CBF1/2/3 are reported to play a dominant role in the cold responsive network of Arabidopsis by directly regulating the expression levels of cold responsive (COR) genes. In this study, we obtained CRISPR/Cas9-mediated loss-of-function mutants of cbf1∼3. Over 3,000 COR genes identified by RNA-seq analysis showed a slight but significant change in their expression levels in the mutants compared to the wild-type plants after being treated at 4 °C for 12 h. The C-repeat (CRT) motif (5'-CCGAC-3') was enriched in promoters of genes that were up-regulated by CBF2 and CBF3 but not in promoters of genes up-regulated by CBF1. These data suggest that CBF2 and CBF3 play a more important role in directing the cold response by regulating different sets of downstream COR genes. More than 2/3 of COR genes were co-regulated by two or three CBFs and were involved mainly in cellular signal transduction and metabolic processes; less than 1/3 of the genes were regulated by one CBF, and those genes up-regulated were enriched in cold-related abiotic stress responses. Our results indicate that CBFs play an important role in the trade-off between cold tolerance and plant growth through the precise regulation of COR genes in the complicated transcriptional network.

Specific peroxidases differentiate Brachypodium distachyon accessions and are associated with drought tolerance traits

BACKGROUND AND AIMS

Brachypodium distachyon (Brachypodium) is a model system for studying cereal, bioenergy, forage and turf grasses. The genetic and evolutionary basis of the adaptation of this wild grass species to drought stress is largely unknown. Peroxidase (POD) may play a role in plant drought tolerance, but whether the allelic variations of genes encoding the specific POD isoenzymes are associated with plant response to drought stress is not well understood. The objectives of this study were to examine natural variation of POD isoenzyme patterns, to identify nucleotide diversity of POD genes and to relate the allelic variation of genes to drought tolerance traits of diverse Brachypodium accessions.

METHODS

Whole-plant drought tolerance and POD activity were examined in contrasting ecotypes. Non-denaturing PAGE and liquid chromatography-mass spectrometry were performed to detect distinct isozymes of POD in 34 accessions. Single nucleotide polymorphisms (SNPs) were identified by comparing DNA sequences of these accessions. Associations of POD genes encoding specific POD isoenzymes with drought tolerance traits were analysed using TASSEL software.

KEY RESULTS

Variations of POD isoenzymes were found among accessions with contrasting drought tolerance, while the most tolerant and susceptible accessions each had their own unique POD isoenzyme band. Eight POD genes were identified and a total of 90 SNPs were found among these genes across 34 accessions. After controlling population structure, significant associations of Bradi3g41340.1 and Bradi1g26870.1 with leaf water content or leaf wilting were identified.

CONCLUSIONS

Brachypodium ecotypes have distinct specific POD isozymes. This may contribute to natural variations of drought tolerance of this species. The role of specific POD genes in differentiating Brachypodium accessions with contrasting drought tolerance could be associated with the general fitness of Brachypodium during evolution.

Population genetics of freeze tolerance among natural populations of Populus balsamifera across the growing season

Protection against freeze damage during the growing season influences the northern range limits of plants. Freeze tolerance and freeze avoidance are the two major freeze resistance strategies. Winter survival strategies have been extensively studied in perennials, but few have addressed them and their genetic basis during the growing season. We examined intraspecific phenotypic variation in freeze resistance of Populus balsamifera across latitude and the growing season. To investigate the molecular basis of this variation, we surveyed nucleotide diversity and examined patterns of gene expression in the poplar C-repeat binding factor (CBF) gene family. Foliar freeze tolerance exhibited latitudinal and seasonal variation indicative of natural genotypic variation. CBF6 showed signatures of recent selective sweep. Of the 46 SNPs surveyed across the six CBF homologs, only CBF2_619 exhibited latitudinal differences consistent with increased freeze tolerance in the north. All six CBF genes were cold inducible, but showed varying patterns of expression across the growing season. Some Poplar CBF homologs exhibited patterns consistent with historical selection and clinal variation in freeze tolerance documented here. However, the CBF genes accounted for only a small amount of the variation, indicating that other genes in this and other molecular pathways likely play significant roles in nature.

Molecular Evidence for Functional Divergence and Decay of a Transcription Factor Derived from Whole-Genome Duplication in Arabidopsis thaliana

Functional divergence between duplicate transcription factors (TFs) has been linked to critical events in the evolution of land plants and can result from changes in patterns of expression, binding site divergence, and/or interactions with other proteins. Although plant TFs tend to be retained post polyploidization, many are lost within tens to hundreds of million years. Thus, it can be hypothesized that some TFs in plant genomes are in the process of becoming pseudogenes. Here, we use a pair of salt tolerance-conferring transcription factors, DWARF AND DELAYED FLOWERING1 (DDF1) and DDF2, that duplicated through paleopolyploidy 50 to 65 million years ago, as examples to illustrate potential mechanisms leading to duplicate retention and loss. We found that the expression patterns of Arabidopsis thaliana (At)DDF1 and AtDDF2 have diverged in a highly asymmetric manner, and AtDDF2 has lost most inferred ancestral stress responses. Consistent with promoter disablement, the AtDDF2 promoter has fewer predicted cis-elements and a methylated repetitive element. Through comparisons of AtDDF1, AtDDF2, and their Arabidopsis lyrata orthologs, we identified significant differences in binding affinities and binding site preference. In particular, an AtDDF2-specific substitution within the DNA-binding domain significantly reduces binding affinity. Cross-species analyses indicate that both AtDDF1 and AtDDF2 are under selective constraint, but among A. thaliana accessions, AtDDF2 has a higher level of nonsynonymous nucleotide diversity compared with AtDDF1. This may be the result of selection in different environments or may point toward the possibility of ongoing functional decay despite retention for millions of years after gene duplication.

Time-dependent deacclimation after cold acclimation in Arabidopsis thaliana accessions

During low temperature exposure, Arabidopsis thaliana and many other plants from temperate climates increase in freezing tolerance in a process termed cold acclimation. However, the correct timing and rate of deacclimation, resulting in loss of freezing tolerance and initiation of growth is equally important for plant fitness and survival. While the molecular basis of cold acclimation has been investigated in detail, much less information is available about deacclimation. We have characterized the responses of 10 natural accessions of Arabidopsis thaliana that vary widely in their freezing tolerance, to deacclimation conditions. Sugar, proline and transcript levels declined sharply over three days in all accessions after transfer of cold acclimated plants to ambient temperatures, while freezing tolerance only declined in tolerant accessions. Correlations between freezing tolerance and the expression levels of COR genes and the content of glucose, fructose and sucrose, as well as many correlations among transcript and solute levels, that were highly significant in cold acclimated plants, were lost during deacclimation. Other correlations persisted, indicating that after three days of deacclimation, plant metabolism had not completely reverted back to the non-acclimated state. These data provide the basis for further molecular and genetic studies to unravel the regulation of deacclimation.

Expression and molecular evolution of two DREB1 genes in black poplar (Populus nigra)

Environmental stresses such as low temperature, drought, and high salinity significantly affect plant growth and yield. As selective forces, these adverse factors play essential roles in shaping phenotypic variation in plant populations. Black poplar (Populus nigra) is an economically and ecologically important forest tree species with widely distributed populations and is thus suitable for experiments detecting evolutionary footprints left by stress. Here, we performed expression and evolutionary analysis of two duplicated DREB A1-subgroup (DREB1) genes, PnDREB68 and PnDREB69, encoding transcription factors that are involved in stress responses. The two genes showed partially overlapping but distinct expression patterns in response to stresses. These genes were strongly and rapidly induced by cold stress in leaves, stems, and roots. In leaf tissue, dehydration stress induced the expression of PnDREB68 but not PnDREB69. PnDREB69 displayed more rapid responses and longer expression durations than PnDREB68 under salt and ABA stress, respectively. Based on single nucleotide polymorphism (SNP) analysis, we found significant population genetic differentiation, with a greater F_ST value (0.09189) for PnDREB69 than for PnDREB68 (0.07743). Nucleotide diversity analysis revealed a two-fold higher π_T for PnDREB68 than for PnDREB69 (0.00563 vs. 0.00243), reflecting strong purifying selection acting on the former. The results suggest that positive selection acted on PnDREB69, as evidenced by neutral testing using Tajima’s D statistic. The distinct selective forces to which each of the genes was subjected may be associated with expression divergence. Linkage disequilibrium (LD) was low for the sequenced region, with a higher level for PnDREB68 than for PnDREB69. Additionally, analysis of the relationship among carbon isotope ratios, SNP classes and gene expression, together with motif and domain analysis, suggested that 14 polymorphisms within the two genes may be candidates for an association study of important traits such as water use efficiency/drought tolerance in black poplar.

Natural variation of C-repeat-binding factor (CBFs) genes is a major cause of divergence in freezing tolerance among a group of Arabidopsis thaliana populations along the Yangtze River in China

We used a monophyletic group of four natural populations of Arabidopsis thaliana expanded from a single ancestor along the Yangtze River c. 90,000 yr ago to study the molecular mechanism of the divergence in their freezing tolerance, in order to gain an insight into the genetic basis of their local adaption to low temperatures. Freezing tolerance assays, measurements of metabolites in the raffinose biosynthesis pathway and transactivation-activity assays of variation in forms of cold-responsive transcription factors were conducted on the four populations. Quantitative trait locus mapping was adopted with F₂ populations of the most- and least freezing-tolerant populations. The degree of freezing tolerance among the four populations was negatively correlated with the lowest monthly average temperature of January in their native habitats, and positively correlated to the expression level of some cold-regulated genes. We identified a major locus harboring three cold-responsive transcription factor genes CBF1-3, and found a nucleotide insertion in CBF2 in all populations except SXcgx, which generated a dysfunctional CBF2 protein. The CBF2 in SXcgx experienced a stronger natural selection in the cooler environment after CBF3 lost its response to low temperature, which possibly reflects a local adaptation of these populations during the expansion from a common ancestor.

Nucleotide diversity and linkage disequilibrium in antioxidant genes of Brachypodium distachyon

Brachypodium distachyon (Brachypodium) is a powerful model system for studying cereal, bioenergy, forage, and turf grasses. Nucleotide diversity (π) and linkage disequilibrium (LD) in candidate genes involved in the antioxidative pathways in this species are not known. The average π for CAT encoding catalase, GPX encoding glutathione peroxidase, DHAR encoding dehydroascorbate reductase, MDHAR encoding monodehydroascorbate reductase, and APX ecoding ascorbate peroxidase was 0.0027 among 19 accessions contrasting for drought tolerance. The highest value of π was found in APX (0.0046) and the lowest π was in MDHAR (0.0006). The average single nucleotide polymorphism (SNP) frequency across these five genes was one SNP per 131 bp between two randomly sampled sequences for the five genes in the sequence length ranging from 1,447 bp to 1,701 bp. The LD decay was slow and extended to a distance of more than 1.2kb for all genes. The neighbor-joining tree analyses of DHAR, MDHAR, and CAT generally separated accessions differing in drought tolerance. The results indicate a putative role of these candidate genes in increasing general fitness of Brachypodium.

Colinearity and similar expression pattern of rice DREB1s reveal their functional conservation in the cold-responsive pathway

The clustered genes C-repeat (CRT) binding factor (CBF)1/dehydration-responsive element binding protein (DREB)1B, CBF2/DREB1C, and CBF3/DREB1A play a central role in cold acclimation and facilitate plant resistance to freezing in Arabidopsis thaliana. Rice (Oryza sativa L.) is very sensitive to low temperatures; enhancing the cold stress tolerance of rice is a key challenge to increasing its yield. In this study, we demonstrate chilling acclimation, a phenomenon similar to Arabidopsis cold acclimation, in rice. To determine whether rice CBF/DREB1 genes participate in this cold-responsive pathway, all putative homologs of Arabidopsis DREB1 genes were filtered from the complete rice genome through a BLASTP search, followed by phylogenetic, colinearity and expression analysis. We thereby identified 10 rice genes as putative DREB1 homologs: nine of these were located in rice genomic regions with some colinearity to the Arabidopsis CBF1–CBF4 region. Expression profiling revealed that six of these genes (Os01g73770, Os02g45450, Os04g48350, Os06g03670, Os09g35010, and Os09g35030) were similarly expressed in response to chilling acclimation and cold stress and were co-expressed with genes involved in cold signalling, suggesting that these DREB1 homologs may be involved in the cold response in rice. The results presented here serve as a prelude towards understanding the function of rice homologs of DREB1 genes in cold-sensitive crops.

Clinal variation in the non-acclimated and cold-acclimated freezing tolerance of Arabidopsis thaliana accessions

Arabidopsis thaliana is a geographically widely spread species consisting of local accessions differing both genetically and phenotypically. These differences may constitute environmental adaptations and a latitudinal cline in freezing tolerance has been shown previously. Many plants, including Arabidopsis, exhibit increased freezing tolerance after cold exposure (cold acclimation). Here we present evidence for geographical clines (both latitudinal and longitudinal) in acclimated (ACC) and non-acclimated (NA) freezing tolerance, estimated from electrolyte leakage measurements on 54 accessions. Leaf Pro contents were not correlated with freezing tolerance, while sugar contents (Glc, Fru, Suc, Raf) were in the ACC, but not the NA state. Expression levels of 14 cold-induced genes were investigated before and after 2 weeks of cold acclimation by quantitative RT-PCR. Expression of the CBF1, 2 and 3 genes was not correlated with freezing tolerance. The expression of some CBF-regulated (COR) genes, however, was correlated specifically with ACC freezing tolerance. A tight correlation between CBF and COR gene expression was only observed under non-acclimating conditions, where CBF and COR expression were also correlated with the expression of PRR5, a component of the circadian clock. Collectively, this study sheds new light on the molecular determinants of plant-freezing tolerance and cold acclimation and their geographical dependence.

Inside the CBF locus in Poaceae

Several molecular evidences have been gathered in Poaceae that point out a central role of the CBF/DREB1 transcription factors in the signal transduction pathways leading to low-temperature tolerance, although to a quite different extent between crops originating from either temperate or tropical climates. A common feature of the CBF/DREB1 genes in Poaceae is their structural organization at the genome level in clusters of tandemly duplicated genes. In temperate cereals such as barley and wheat, expansion of specific multigene phylogenetic clades of CBFs that map at the Frost Resistance-2 locus has been exclusively observed. In addition, copy number variants of CBF genes between frost resistant and frost sensitive genotypes raise the question if multiple copies of the CBF/DREB1s are required to ensure freezing tolerance. On the other hand, in crops of tropical origin such as rice and maize, a smaller or less-responsive CBF regulon may have evolved, and different mechanisms might determine chilling tolerance. In this review, recent advances on the organization and diversity at the CBF cluster locus in the grasses are provided and discussed.