Expansion and diversification of the Populus R2R3-MYB family of transcription factors

Genome-wide analysis of citrus R2R3MYB genes and their spatiotemporal expression under stresses and hormone treatments

The R2R3MYB proteins represent one of the largest families of transcription factors, which play important roles in plant growth and development. Although genome-wide analysis of this family has been conducted in many species, little is known about R2R3MYB genes in citrus, In this study, 101 R2R3MYB genes has been identified in the citrus (Citrus sinesis and Citrus clementina) genomes, which are almost equal to the number of rice. Phylogenetic analysis revealed that they could be subdivided into 21 subgroups. The evolutionary relationships and the intro-exon organizations were also analyzed, revealing strong gene conservation but also the expansions of particular functional genes during the plant evolution. Tissue-specific expression profiles showed that 95 citrus R2R3MYB genes were expressed in at least one tissue and the other 6 genes showed very low expression in all tissues tested, suggesting that citrus R2R3MYB genes play important roles in the development of all citrus organs. The transcript abundance level analysis during abiotic conditions (NaCl, abscisic acid, jasmonic acid, drought and low temperature) identified a group of R2R3MYB genes that responded to one or multiple treatments, which showed a promising for improving citrus adaptation to stresses. Our results provided an essential foundation for the future selection of the citrus R2R3MYB genes for cloning and functional dissection with an aim of uncovering their roles in citrus growth and development.

Genome-wide identification and expression analysis of the IQD gene family in Populus trichocarpa

IQD proteins are downstream targets of calcium sensors, which play important roles in development and responses to environmental cues in plants. Comprehensive analyses of IQD genes have been conducted in Arabidopsis, rice, tomato, and Brachypodium distachyon, but have not been reported from poplar. The availability of the Populus trichocarpa genome sequence allowed us to conduct phylogenetic, gene structure, chromosomal location, and microarray analyses of the predicted IQD genes in P. trichocarpa. We identified 40 IQD genes (PtIQD1-40) classified in four subfamilies (I-IV). Gene structure and protein motif analyses showed that these genes are relatively conserved within each subfamily. The 40 PtIQD genes are distributed on 18 of the 19 chromosomes, with 16 gene pairs involved in segmental duplication events. The Ka/Ks ratios of the 16 segmentally-duplicated gene pairs show that the duplicated pairs underwent purifying selection with restrictive functional divergence after the duplication events. Analyses of microarray data for 38 PtIQD genes showed tissue/organ-specific expression patterns. We also performed quantitative real-time RT-PCR (qRT-PCR) analyses of twelve selected PtIQD genes in plants treated with MeJA and PEG in order to explore their stress-related expression patterns. Our results will be valuable for further analysis of poplar IQD genes to characterize their important biological functions.

The role of water channel proteins in facilitating recovery of leaf hydraulic conductance from water stress in Populus trichocarpa

Gas exchange is constrained by the whole-plant hydraulic conductance (Kplant). Leaves account for an important fraction of Kplant and may therefore represent a major determinant of plant productivity. Leaf hydraulic conductance (Kleaf) decreases with increasing water stress, which is due to xylem embolism in leaf veins and/or the properties of the extra-xylary pathway. Water flow through living tissues is facilitated and regulated by water channel proteins called aquaporins (AQPs). Here we assessed changes in the hydraulic conductance of Populus trichocarpa leaves during a dehydration-rewatering episode. While leaves were highly sensitive to drought, Kleaf recovered only 2 hours after plants were rewatered. Recovery of Kleaf was absent when excised leaves were bench-dried and subsequently xylem-perfused with a solution containing AQP inhibitors. We examined the expression patterns of 12 highly expressed AQP genes during a dehydration-rehydration episode to identify isoforms that may be involved in leaf hydraulic adjustments. Among the AQPs tested, several genes encoding tonoplast intrinsic proteins (TIPs) showed large increases in expression in rehydrated leaves, suggesting that TIPs contribute to reversing drought-induced reductions in Kleaf. TIPs were localized in xylem parenchyma, consistent with a role in facilitating water exchange between xylem vessels and adjacent living cells. Dye uptake experiments suggested that reversible embolism formation in minor leaf veins contributed to the observed changes in Kleaf.

Genome-wide classification and evolutionary and expression analyses of citrus MYB transcription factor families in sweet orange

MYB family genes are widely distributed in plants and comprise one of the largest transcription factors involved in various developmental processes and defense responses of plants. To date, few MYB genes and little expression profiling have been reported for citrus. Here, we describe and classify 177 members of the sweet orange MYB gene (CsMYB) family in terms of their genomic gene structures and similarity to their putative Arabidopsis orthologs. According to these analyses, these CsMYBs were categorized into four groups (4R-MYB, 3R-MYB, 2R-MYB and 1R-MYB). Gene structure analysis revealed that 1R-MYB genes possess relatively more introns as compared with 2R-MYB genes. Investigation of their chromosomal localizations revealed that these CsMYBs are distributed across nine chromosomes. Sweet orange includes a relatively small number of MYB genes compared with the 198 members in Arabidopsis, presumably due to a paralog reduction related to repetitive sequence insertion into promoter and non-coding transcribed region of the genes. Comparative studies of CsMYBs and Arabidopsis showed that CsMYBs had fewer gene duplication events. Expression analysis revealed that the MYB gene family has a wide expression profile in sweet orange development and plays important roles in development and stress responses. In addition, 337 new putative microsatellites with flanking sequences sufficient for primer design were also identified from the 177 CsMYBs. These results provide a useful reference for the selection of candidate MYB genes for cloning and further functional analysis forcitrus.

Populus tremula (European aspen) shows no evidence of sexual dimorphism

BACKGROUND

Evolutionary theory suggests that males and females may evolve sexually dimorphic phenotypic and biochemical traits concordant with each sex having different optimal strategies of resource investment to maximise reproductive success and fitness. Such sexual dimorphism would result in sex biased gene expression patterns in non-floral organs for autosomal genes associated with the control and development of such phenotypic traits.

RESULTS

We examined morphological, biochemical and herbivory traits to test for sexually dimorphic resource allocation strategies within collections of sexually mature and immature Populus tremula (European aspen) trees. In addition we profiled gene expression in mature leaves of sexually mature wild trees using whole-genome oligonucleotide microarrays and RNA-Sequencing.

CONCLUSIONS

We found no evidence of sexual dimorphism or differential resource investment strategies between males and females in either sexually immature or mature trees. Similarly, single-gene differential expression and machine learning approaches revealed no evidence of large-scale sex biased gene expression. However, two significantly differentially expressed genes were identified from the RNA-Seq data, one of which is a robust diagnostic marker of sex in P. tremula.

Identification and molecular characterization of MYB Transcription Factor Superfamily in C4 model plant foxtail millet (Setaria italica L.)

MYB proteins represent one of the largest transcription factor families in plants, playing important roles in diverse developmental and stress-responsive processes. Considering its significance, several genome-wide analyses have been conducted in almost all land plants except foxtail millet. Foxtail millet (Setaria italica L.) is a model crop for investigating systems biology of millets and bioenergy grasses. Further, the crop is also known for its potential abiotic stress-tolerance. In this context, a comprehensive genome-wide survey was conducted and 209 MYB protein-encoding genes were identified in foxtail millet. All 209 S. italica MYB (SiMYB) genes were physically mapped onto nine chromosomes of foxtail millet. Gene duplication study showed that segmental- and tandem-duplication have occurred in genome resulting in expansion of this gene family. The protein domain investigation classified SiMYB proteins into three classes according to number of MYB repeats present. The phylogenetic analysis categorized SiMYBs into ten groups (I - X). SiMYB-based comparative mapping revealed a maximum orthology between foxtail millet and sorghum, followed by maize, rice and Brachypodium. Heat map analysis showed tissue-specific expression pattern of predominant SiMYB genes. Expression profiling of candidate MYB genes against abiotic stresses and hormone treatments using qRT-PCR revealed specific and/or overlapping expression patterns of SiMYBs. Taken together, the present study provides a foundation for evolutionary and functional characterization of MYB TFs in foxtail millet to dissect their functions in response to environmental stimuli.

Genome-wide identification and characterisation of R2R3-MYB genes in sugar beet (Beta vulgaris)

BACKGROUND

The R2R3-MYB genes comprise one of the largest transcription factor gene families in plants, playing regulatory roles in plant-specific developmental processes, metabolite accumulation and defense responses. Although genome-wide analysis of this gene family has been carried out in some species, the R2R3-MYB genes in Beta vulgaris ssp. vulgaris (sugar beet) as the first sequenced member of the order Caryophyllales, have not been analysed heretofore.

RESULTS

We present a comprehensive, genome-wide analysis of the MYB genes from Beta vulgaris ssp. vulgaris (sugar beet) which is the first species of the order Caryophyllales with a sequenced genome. A total of 70 R2R3-MYB genes as well as genes encoding three other classes of MYB proteins containing multiple MYB repeats were identified and characterised with respect to structure and chromosomal organisation. Also, organ specific expression patterns were determined from RNA-seq data. The R2R3-MYB genes were functionally categorised which led to the identification of a sugar beet-specific clade with an atypical amino acid composition in the R3 domain, putatively encoding betalain regulators. The functional classification was verified by experimental confirmation of the prediction that the R2R3-MYB gene Bv_iogq encodes a flavonol regulator.

CONCLUSIONS

This study provides the first step towards cloning and functional dissection of the role of MYB transcription factor genes in the nutritionally and evolutionarily interesting species B. vulgaris. In addition, it describes the flavonol regulator BvMYB12, being the first sugar beet R2R3-MYB with an experimentally proven function.

Conservation and functional influence of alternative splicing in wood formation of Populus and Eucalyptus

BACKGROUND

Wood formation in tree species is regulated by multiple factors at various layers. Alternative splicing (AS) occurs within a large number of genes in wood formation. However, the functional implications and conservation of the AS occurrence are not well understood.

RESULTS

In this study, we profiled AS events in wood-forming tissues of Populus and Eucalyptus, and analyzed their functional implications as well as inter-species conservation. 28.3% and 20.7% of highly expressed transcripts in the developing xylem of Populus and Eucalyptus respectively were affected by AS events. Around 42% of the AS events resulted in changes to the original reading frame. 25.0% (in Populus) and 26.8% (in Eucalyptus) of the AS events may cause protein domain modification. In the process of wood formation, about 28% of AS-occurring genes were putative orthologs and 71 conserved AS events were identified in the two species.

CONCLUSION

Through analysis of AS events in developing xylem of two tree species, this study reveals an array of new information regarding AS occurrence and function in tree development.

Characterization of an apple TT2-type R2R3 MYB transcription factor functionally similar to the poplar proanthocyanidin regulator PtMYB134

The apple MdMYB9 gene encodes a positive regulator of proanthocyanidin synthesis that activates anthocyanidin reductase promoters from apple and poplar via interaction with basic helix-loop-helix proteins. The regulation of proanthocyanidins (PAs, condensed tannins) is of great importance in food plants due to the many benefits of PAs in the human diet. Two candidate flavonoid MYB regulators, MdMYB9 and MdMYB11, were cloned from apple (Malus × domestica) based on their similarity to known MYB PA regulators. Transcript accumulation of both MdMYB9 and MdMYB11 was induced by high light and wounding, similar to the poplar (Populus spp) PA regulator PtMYB134. In transient activation assays with various basic helix-loop-helix (bHLH) co-regulators, MdMYB9 activated apple and poplar anthocyanidin reductase (ANR) promoters, while MdMYB11 showed no activity. Potential transcription factor binding elements were found within several ANR promoters, and the importance of the bHLH binding site (E-box) on ANR promoter activation was demonstrated via mutational analysis. The ability of MdMYB9 and PtMYB134 to reciprocally activate ANR promoters from both apple and poplar and to partner with heterologous bHLH co-factors from these plants confirms the high degree of conservation of PA regulatory complexes across species. The similarity in apple and poplar PA regulation suggests that regulatory genes from poplar could be effectively employed for metabolic engineering of the PA pathway in apple.

SpMYB overexpression in tobacco plants leads to altered abiotic and biotic stress responses

The MYB transcription factors are involved in various plant biochemistry and physiology processes and play a central role in plant defense response. In the present study, a full-length cDNA sequence of a MYB gene, designated as SpMYB, was isolated from tomato. SpMYB encodes the R2R3-type protein consisting of 328 amino acids. The expression level of SpMYB was strongly induced by fungal pathogens. Transgenic tobacco plants overexpressing SpMYB had an enhanced salt and drought stress tolerance compared with wild-type plants, and showed significantly improved resistance to Alternaria alternate. Further analysis revealed that transgenic tobaccos exhibited less accumulation of malondialdehyde (MDA) and more accumulation of superoxide dismutase (SOD), peroxidase (POD) and phenylalanine ammonia-lyase (PAL) after inoculation with A. alternate. Meanwhile, changes in some photosynthetic parameters, such as photosynthetic rate (Pn), transpiration rate (Tr) and intercellular CO2 concentration (Ci) were also found in the transgenic tobaccos. Furthermore, transgenic tobaccos constitutively accumulated higher levels of pathogenesis-related (PR) gene transcripts, such as PR1 and PR2. The results suggested that the tomato SpMYB transcription factor plays an important role in responses to abiotic and biotic stress.

The transcriptional response of apple alcohol acyltransferase (MdAAT2) to salicylic acid and ethylene is mediated through two apple MYB TFs in transgenic tobacco

Volatile esters are major factors affecting the aroma of apple fruits, and alcohol acyltransferases (AATs) are key enzymes involved in the last steps of ester biosynthesis. The expression of apple AAT (MdAAT2) is known to be induced by salicylic acid (SA) or ethylene in apple fruits, although the mechanism of its transcriptional regulation remains elusive. In this study, we reveal that two apple transcription factors (TFs), MdMYB1 and MdMYB6, are involved in MdAAT2 promoter response to SA and ethylene in transgenic tobacco. According to electrophoretic mobility shift assays, MdMYB1 or MdMYB6 can directly bind in vitro to MYB binding sites in the MdAAT2 promoter. In vivo, overexpression of the two MYB TFs can greatly enhance MdAAT2 promoter activity, as demonstrated by dual luciferase reporter assays in transgenic tobacco. In contrast to the promoter of MdMYB1 or MdMYB6, the MdAAT2 promoter cannot be induced by SA or ethephon (ETH) in transgenic tobacco, even in stigmas in which the MdAAT2 promoter can be highly induced under normal conditions. However, the induced MYB TFs can dramatically enhance MdAAT2 promoter activity under SA or ETH treatment. We conclude that MdMYB1 and MdMYB6 function in MdAAT2 responses to SA and ethylene in transgenic tobacco, suggesting that a similar regulation mechanism may exist in apple.

R2R3-MYB gene pairs in Populus: evolution and contribution to secondary wall formation and flowering time

In plants, the R2R3-MYB gene family contains many pairs of paralogous genes, which play the diverse roles in developmental processes and environmental responses. The paper reports the characterization of 81 pairs of Populus R2R3-MYB genes. Chromosome placement, phylogenetic, and motif structure analyses showed that these gene pairs resulted from multiple types of gene duplications and had five different gene fates. Tissue expression patterns revealed that most duplicated genes were specifically expressed in the tissues examined. qRT-PCR confirmed that nine pairs were highly expressed in xylem, of which three pairs (PdMYB10/128, PdMYB90/167, and PdMYB92/125) were further functionally characterized. The six PdMYBs were localized to the nucleus and had transcriptional activities in yeast. The heterologous expression of PdMYB10 and 128 in Arabidopsis increased stem fibre cell-wall thickness and delayed flowering. In contrast, overexpression of PdMYB90, 167, 92, and 125 in Arabidopsis decreased stem fibre and vessel cell-wall thickness and promoted flowering. Cellulose, xylose, and lignin contents were changed in overexpression plants. The expression levels of several genes involved in secondary wall formation and flowering were affected by the overexpression of the six PdMYBs in Arabidopsis. This study addresses the diversity of gene duplications in Populus R2R3-MYBs and the roles of these six genes in secondary wall formation and flowering control.

Genome-wide identification and characterization of R2R3MYB family in Solanum lycopersicum

The R2R3MYB proteins comprise one of the largest families of transcription factors and play regulatory roles in developmental processes and defense responses in plants. However, there has been relatively little effort to systematically carry out comprehensive genomic and functional analyses of these genes in tomato (Solanum lycopersicum L.), a reference species for Solanaceae plants, and the model plant for fruit development. In this study, a total of 121 R2R3MYB genes were identified in the tomato genome released recently and further classified into 29 subgroups based on the phylogenetic analysis of the complete protein sequences. Phylogenetic comparison of the members of this superfamily among tomato, Arabidopsis, grape, rice, poplar, soybean, cucumber and apple revealed that the putative functions of some tomato R2R3MYB proteins were clustered into the Arabidopsis functional clades. The chromosome distribution pattern revealed that tomato R2R3MYB genes were enriched on several chromosomes and 52 % of the family members were tandemly duplicated genes. Tissue specificity or different expression levels of SlR2R3MYBs in different tissues suggested differential regulation of tissue development as well as metabolic regulation. The transcript abundance level analysis during abiotic conditions identified a group of R2R3MYB genes that responded to one or more treatments suggesting that the SlR2R3MYBs played major roles in the plant response to abiotic conditions and involved in signal transduction pathways. This study not only provides a solid foundation for further functional dissection of tomato R2R3MYB family genes, but may also be profitable for, in the future, the improvement of tomato stress tolerance and fruit quality.

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.

Regulation of secondary cell wall biosynthesis by poplar R2R3 MYB transcription factor PtrMYB152 in Arabidopsis

Poplar has 192 annotated R2R3 MYB genes, of which only three have been shown to play a role in the regulation of secondary cell wall formation. Here we report the characterization of PtrMYB152, a poplar homolog of the Arabidopsis R2R3 MYB transcription factor AtMYB43, in the regulation of secondary cell wall biosynthesis. The expression of PtrMYB152 in secondary xylem is about 18 times of that in phloem. When expressed in Arabidopsis under the control of either 35S or PtrCesA8 promoters, PtrMYB152 increased secondary cell wall thickness, which is likely caused by increased lignification. Accordingly, elevated expression of genes encoding sets of enzymes in secondary wall biosynthesis were observed in transgenic plants expressing PtrMYB152. Arabidopsis protoplast transfection assays suggested that PtrMYB152 functions as a transcriptional activator. Taken together, our results suggest that PtrMYB152 may be part of a regulatory network activating expression of discrete sets of secondary cell wall biosynthesis genes.

Comparative genomic analysis of the R2R3 MYB secondary cell wall regulators of Arabidopsis, poplar, rice, maize, and switchgrass

BACKGROUND

R2R3 MYB proteins constitute one of the largest plant transcription factor clades and regulate diverse plant-specific processes. Several R2R3 MYB proteins act as regulators of secondary cell wall (SCW) biosynthesis in Arabidopsis thaliana (At), a dicotyledenous plant. Relatively few studies have examined SCW R2R3 MYB function in grasses, which may have diverged from dicots in terms of SCW regulatory mechanisms, as they have in cell wall composition and patterning. Understanding cell wall regulation is especially important for improving lignocellulosic bioenergy crops, such as switchgrass.

RESULTS

Here, we describe the results of applying phylogenic, OrthoMCL, and sequence identity analyses to classify the R2R3 MYB family proteins from the annotated proteomes of Arabidposis, poplar, rice, maize and the initial genome (v0.0) and translated transcriptome of switchgrass (Panicum virgatum). We find that the R2R3 MYB proteins of the five species fall into 48 subgroups, including three dicot-specific, six grass-specific, and two panicoid grass-expanded subgroups. We observe four classes of phylogenetic relationships within the subgroups of known SCW-regulating MYB proteins between Arabidopsis and rice, ranging from likely one-to-one orthology (for AtMYB26, AtMYB103, AtMYB69) to no homologs identifiable (for AtMYB75). Microarray data for putative switchgrass SCW MYBs indicate that many maintain similar expression patterns with the Arabidopsis SCW regulators. However, some of the switchgrass-expanded candidate SCW MYBs exhibit differences in gene expression patterns among paralogs, consistent with subfunctionalization. Furthermore, some switchgrass representatives of grass-expanded clades have gene expression patterns consistent with regulating SCW development.

CONCLUSIONS

Our analysis suggests that no single comparative genomics tool is able to provide a complete picture of the R2R3 MYB protein family without leaving ambiguities, and establishing likely false-negative and -positive relationships, but that used together a relatively clear view emerges. Generally, we find that most R2R3 MYBs that regulate SCW in Arabidopsis are likely conserved in the grasses. This comparative analysis of the R2R3 MYB family will facilitate transfer of understanding of regulatory mechanisms among species and enable control of SCW biosynthesis in switchgrass toward improving its biomass quality.

Genome-wide characterization and comparative analysis of R2R3-MYB transcription factors shows the complexity of MYB-associated regulatory networks in Salvia miltiorrhiza

Background

MYB is the largest plant transcription factor gene family playing vital roles in plant growth and development. However, it has not been systematically studied in Salvia miltiorrhiza, an economically important medicinal plant.

Results

Here we report the genome-wide identification and characterization of 110 R2R3-MYBs, the largest subfamily of MYBs in S. miltiorrhiza. The MYB domain and other motifs of SmMYBs are largely conserved with Arabidopsis AtMYBs, whereas the divergence of SmMYBs and AtMYBs also exists, suggesting the conservation and diversity of plant MYBs. SmMYBs and AtMYBs may be classified into 37 subgroups, of which 31 include proteins from S. miltiorrhiza and Arabidopsis, whereas 6 are specific to a species, indicating that the majority of MYBs play conserved roles, while others may exhibit species-specialized functions. SmMYBs are differentially expressed in various tissues of S. miltiorrhiza. The expression profiles are largely consistent with known functions of their Arabidopsis counterparts. The expression of a subset of SmMYBs is regulated by microRNAs, such as miR159, miR319, miR828 and miR858. Based on functional conservation of MYBs in a subgroup, SmMYBs potentially involved in the biosynthesis of bioactive compounds were identified.

Conclusions

A total of 110 R2R3-MYBs were identified and analyzed. The results suggest the complexity of MYB-mediated regulatory networks in S. miltiorrhiza and provide a foundation for understanding the regulatory mechanism of SmMYBs.

A unique mutation in a MYB gene cosegregates with the nectarine phenotype in peach

Nectarines play a key role in peach industry; the fuzzless skin has implications for consumer acceptance. The peach/nectarine (G/g) trait was described as monogenic and previously mapped on chromosome 5. Here, the position of the G locus was delimited within a 1.1 cM interval (635 kb) based on linkage analysis of an F₂ progeny from the cross ‘Contender’ (C, peach) x ‘Ambra’ (A, nectarine). Careful inspection of the genes annotated in the corresponding genomic sequence (Peach v1.0), coupled with variant discovery, led to the identification of MYB gene PpeMYB25 as a candidate for trichome formation on fruit skin. Analysis of genomic re-sequencing data from five peach/nectarine accessions pointed to the insertion of a LTR retroelement in exon 3 of the PpeMYB25 gene as the cause of the recessive glabrous phenotype. A functional marker (indelG) developed on the LTR insertion cosegregated with the trait in the CxA F₂ progeny and was validated on a broad panel of genotypes, including all known putative donors of the nectarine trait. This marker was shown to efficiently discriminate between peach and nectarine plants, indicating that a unique mutational event gave rise to the nectarine trait and providing a useful diagnostic tool for early seedling selection in peach breeding programs.

A conserved network of transcriptional activators and repressors regulates anthocyanin pigmentation in eudicots

Plants require sophisticated regulatory mechanisms to ensure the degree of anthocyanin pigmentation is appropriate to myriad developmental and environmental signals. Central to this process are the activity of MYB-bHLH-WD repeat (MBW) complexes that regulate the transcription of anthocyanin genes. In this study, the gene regulatory network that regulates anthocyanin synthesis in petunia (Petunia hybrida) has been characterized. Genetic and molecular evidence show that the R2R3-MYB, MYB27, is an anthocyanin repressor that functions as part of the MBW complex and represses transcription through its C-terminal EAR motif. MYB27 targets both the anthocyanin pathway genes and basic-helix-loop-helix (bHLH) ANTHOCYANIN1 (AN1), itself an essential component of the MBW activation complex for pigmentation. Other features of the regulatory network identified include inhibition of AN1 activity by the competitive R3-MYB repressor MYBx and the activation of AN1, MYB27, and MYBx by the MBW activation complex, providing for both reinforcement and feedback regulation. We also demonstrate the intercellular movement of the WDR protein (AN11) and R3-repressor (MYBx), which may facilitate anthocyanin pigment pattern formation. The fundamental features of this regulatory network in the Asterid model of petunia are similar to those in the Rosid model of Arabidopsis thaliana and are thus likely to be widespread in the Eudicots.

In silico identification and characterization of N-Terminal acetyltransferase genes of poplar (Populus trichocarpa)

N-terminal acetyltransferase (Nats) complex is responsible for protein N-terminal acetylation (Nα-acetylation), which is one of the most common covalent modifications of eukaryotic proteins. Although genome-wide investigation and characterization of Nat catalytic subunits (CS) and auxiliary subunits (AS) have been conducted in yeast and humans they remain unexplored in plants. Here we report on the identification of eleven genes encoding eleven putative Nat CS polypeptides, and five genes encoding five putative Nat AS polypeptides in Populus. We document that the expansion of Nat CS genes occurs as duplicated blocks distributed across 10 of the 19 poplar chromosomes, likely only as a result of segmental duplication events. Based on phylogenetic analysis, poplar Nat CS were assigned to six subgroups, which corresponded well to the Nat CS types (CS of Nat A–F), being consistent with previous reports in humans and yeast. In silico analysis of microarray data showed that in the process of normal development of the poplar, their Nat CS and AS genes are commonly expressed at one relatively low level but share distinct tissue-specific expression patterns. This exhaustive survey of Nat genes in poplar provides important information to assist future studies on their functional role in poplar.

Evolutionary classification of ammonium, nitrate, and peptide transporters in land plants

Background

Nitrogen uptake, reallocation within the plant, and between subcellular compartments involves ammonium, nitrate and peptide transporters. Ammonium transporters are separated into two distinct families (AMT1 and AMT2), each comprised of five members on average in angiosperms. Nitrate transporters also form two discrete families (NRT1 and NRT2), with angiosperms having four NRT2s, on average. NRT1s share an evolutionary history with peptide transporters (PTRs). The NRT1/PTR family in land plants usually has more than 50 members and contains also members with distinct activities, such as glucosinolate and abscisic acid transport.

Results

Phylogenetic reconstructions of each family across 20 land plant species with available genome sequences were supplemented with subcellular localization and transmembrane topology predictions. This revealed that both AMT families diverged prior to the separation of bryophytes and vascular plants forming two distinct clans, designated as supergroups, each. Ten supergroups were identified for the NRT1/PTR family. It is apparent that nitrate and peptide transport within the NRT1/PTR family is polyphyletic, that is, nitrate and/or peptide transport likely evolved multiple times within land plants. The NRT2 family separated into two distinct clans early in vascular plant evolution. Subsequent duplications occurring prior to the eudicot/monocot separation led to the existence of two AMT1, six AMT2, 31 NRT1/PTR, and two NRT2 clans, designated as groups.

Conclusion

Phylogenetic separation of groups suggests functional divergence within the angiosperms for each family. Distinct groups within the NRT1/PTR family appear to separate peptide and nitrate transport activities as well as other activities contained within the family, for example nitrite transport. Conversely, distinct activities, such as abscisic acid and glucosinolate transport, appear to have recently evolved from nitrate transporters.

DRM1 and DRM2 expression regulation: potential role of splice variants in response to stress and environmental factors in Arabidopsis

The DRM1/ARP gene family is increasingly exhibiting associations with stress conditions at the transcript level. Traditionally correlated with dormancy, increases in transcript levels in response to various treatments have also been reported in various species. As alternative transcript splicing is common in stress conditions, the splice variants of AtDRM1 and AtDRM2 were assessed further in this study. A previously undescribed splice variant of AtDRM1 (AtDRM1.6) is introduced in this work. In silico analyses of predicted protein sequence of all splice variants showed that all variants retain the predicted intrinsically disordered nature. Transcriptional studies of AtDRM1 and AtDRM2 in response to a wide range of abiotic, physical and hormonal treatments showed that AtDRM1.6 is differentially regulated at the transcriptional level compared with other splice variants. Promoter analyses demonstrated AtDRM1 light regulation via the upstream promoter sequence.

Plant 4/1 protein: potential player in intracellular, cell-to-cell and long-distance signaling

Originally isolated as a result of its ability to interact with the movement protein of Tomato spotted wilt virus in a yeast two-hybrid system, the 4/1 protein is proving to be an excellent tool for studying intracellular protein trafficking and intercellular communication. Expression of 4/1 in vivo is tightly regulated, first appearing in the veins of the cotyledon and later in the vasculature of the leaf and stem in association with the xylem parenchyma and phloem parenchyma. Structural studies indicate that 4/1 proteins contain as many as five coiled-coil (CC) domains; indeed, the highest level of sequence identity among 4/1 proteins involves their C-terminal CC domains, suggesting that protein-protein interaction is important for biological function. Recent data predict that the tertiary structure of this C-terminal CC domain is strikingly similar to that of yeast protein She2p; furthermore, like She2p, 4/1 protein exhibits RNA-binding activity, and mutational analysis has shown that the C-terminal CC domain is responsible for RNA binding. The 4/1 protein contains a nuclear export signal. Additional microscopy studies involving leptomycin and computer prediction suggest the presence of a nuclear localization signal as well.

Transcriptome sequencing and phylogenetic analysis of floral and leaf MIKC(C) MADS-box and R2R3 MYB transcription factors from the monocot Iris fulva

The Louisiana Irises serve as an important system for the study of the evolutionary processes of speciation, including reproductive isolation, hybridization, and adaptation. Sequencing methods today allow for the generation of resources key to elucidating the genetic basis of these phenomena. Here we describe the transcriptomes of floral and young leaf tissue from Iris fulva generated by massively parallel pyrosequencing. In order to identify potential candidates for the study of reproductive isolation and adaptation in the Louisiana Irises we phylogenetically analyzed the type II MIKC(C) MADS-box and R2R3 MYB transcription factors expressed in these tissues. A total of 25 Iris MIKC(C) MADS-box genes in 9 clades and 42 Iris R2R3 MYB genes in 19 clades were identified. Through the identification of eudicot and monocot specific clades, these analyses contribute to our understanding of the evolution of these prominent transcription factor families in the angiosperms.

Insight into the genetic components of community genetics: QTL mapping of insect association in a fast-growing forest tree

Identifying genetic sequences underlying insect associations on forest trees will improve the understanding of community genetics on a broad scale. We tested for genomic regions associated with insects in hybrid poplar using quantitative trait loci (QTL) analyses conducted on data from a common garden experiment. The F2 offspring of a hybrid poplar (Populus trichocarpa x P. deltoides) cross were assessed for seven categories of insect leaf damage at two time points, June and August. Positive and negative correlations were detected among damage categories and between sampling times. For example, sap suckers on leaves in June were positively correlated with sap suckers on leaves (P<0.001) but negatively correlated with skeletonizer damage (P<0.01) in August. The seven forms of leaf damage were used as a proxy for seven functional groups of insect species. Significant variation in insect association occurred among the hybrid offspring, including transgressive segregation of susceptibility to damage. NMDS analyses revealed significant variation and modest broad-sense heritability in insect community structure among genets. QTL analyses identified 14 genomic regions across 9 linkage groups that correlated with insect association. We used three genomics tools to test for putative mechanisms underlying the QTL. First, shikimate-phenylpropanoid pathway genes co-located to 9 of the 13 QTL tested, consistent with the role of phenolic glycosides as defensive compounds. Second, two insect association QTL corresponded to genomic hotspots for leaf trait QTL as identified in previous studies, indicating that, in addition to biochemical attributes, leaf morphology may influence insect preference. Third, network analyses identified categories of gene models over-represented in QTL for certain damage types, providing direction for future functional studies. These results provide insight into the genetic components involved in insect community structure in a fast-growing forest tree.

Gene structures, classification, and expression models of the DREB transcription factor subfamily in Populus trichocarpa

We identified 75 dehydration-responsive element-binding (DREB) protein genes in Populus trichocarpa. We analyzed gene structures, phylogenies, domain duplications, genome localizations, and expression profiles. The phylogenic construction suggests that the PtrDREB gene subfamily can be classified broadly into six subtypes (DREB A-1 to A-6) in Populus. The chromosomal localizations of the PtrDREB genes indicated 18 segmental duplication events involving 36 genes and six redundant PtrDREB genes were involved in tandem duplication events. There were fewer introns in the PtrDREB subfamily. The motif composition of PtrDREB was highly conserved in the same subtype. We investigated expression profiles of this gene subfamily from different tissues and/or developmental stages. Sixteen genes present in the digital expression analysis had high levels of transcript accumulation. The microarray results suggest that 18 genes were upregulated. We further examined the stress responsiveness of 15 genes by qRT-PCR. A digital northern analysis showed that the PtrDREB17, 18, and 32 genes were highly induced in leaves under cold stress, and the same expression trends were shown by qRT-PCR. Taken together, these observations may lay the foundation for future functional analyses to unravel the biological roles of Populus' DREB genes.

Roles of JnRAP2.6-like from the transition zone of black walnut in hormone signaling

An EST sequence, designated JnRAP2-like, was isolated from tissue at the heartwood/sapwood transition zone (TZ) in black walnut (Juglans nigra L). The deduced amino acid sequence of JnRAP2-like protein consists of a single AP2-containing domain with significant similarity to conserved AP2/ERF DNA-binding domains in other species. Based on multiple sequence alignment, JnRAP2-like appears to be an ortholog of RAP2.6L (At5g13330), which encodes an ethylene response element binding protein in Arabidopsis thaliana. Real-time PCR revealed that the JnRAP2-like was expressed most abundantly in TZ of trees harvested in fall when compared with other xylem tissues harvested in the fall or summer. Independent transgenic lines over-expressing JnRAP2-like in Arabidopsis developed dramatic ethylene-related phenotypes when treated with 50 µM methyl jasmonate (MeJA). Taken together, these results indicated that JnRAP2-like may participate in the integration of ethylene and jasmonate signals in the xylem and other tissues. Given the role of ethylene in heartwood formation, it is possible JnRAP2-like expression in the transition zone is part of the signal transduction pathway leading to heartwood formation in black walnut.

Functional characterization of the poplar R2R3-MYB transcription factor PtoMYB216 involved in the regulation of lignin biosynthesis during wood formation

Because of the importance of wood in many industrial applications, tremendous studies have been performed on wood formation, especially in lignin biosynthesis. MYB transcription factors (TFs), which consist of a large family of plant TFs, have been reported to directly regulate lignin biosynthetic genes in a number of plants. In this study, we describe the cloning and functional characterization of PtoMYB216, a cDNA isolated from Chinese white poplar (Populus tomentosa Carr.). PtoMYB216 encodes a protein belonging to the R2R3-MYB family and displays significant similarity with other MYB factors shown to regulate lignin synthesis in Arabidopsis. Gene expression profiling studies showed that PtoMYB216 mRNA is specifically expressed during secondary wall formation in wood. The 1.8-kb promoter sequence of PtoMYB216 was fused to the GUS coding sequence and introduced into wild-type A. thaliana. GUS expression was shown to be restricted to tissues undergoing secondary cell wall formation. Overexpression of PtoMYB216 specifically activated the expression of the upstream genes in the lignin biosynthetic pathway and resulted in ectopic deposition of lignin in cells that are normally unligninified. These results suggest that PtoMYB216 is specific transcriptional activators of lignin biosynthesis and involved in the regulation of wood formation in poplar.

Genome-wide analysis of the CCCH zinc finger gene family in Medicago truncatula

In this study, we identified 34 CCCH Znf genes in Medicago truncatula and the results of semi-quantitative RT-PCR revealed that the expression patterns of subfamily VI members were diverse. CCCH-type zinc finger (Znf) proteins are specific transcriptional factors with a typical motif consisting of three cysteine residues and one histidine residue. Increasing evidences have revealed that CCCH Znf proteins participated in the regulation of plant growth, developmental processes and environmental responses. Survey and characterization of CCCH Znf genes in leguminous species would facilitate a better understanding of the evolutionary processes and functions of this gene family. In this study, we performed a comprehensive analysis of CCCH Znf genes in M. truncatula by describing the phylogenetic relationships, chromosomal location and gene structure of each family member. A total of 34 CCCH Znf genes were identified in the latest M. truncatula genome sequence. The 34 predicted members were clustered into nine subfamilies based on their phylogenetic analysis and structure features. In addition, the 34 Medicago CCCH Znf genes were found to be unevenly distributed on eight chromosomes. Furthermore, the expression profiles of subfamily VI were investigated under different stress conditions (PEG-6000, NaCl and ABA) by using semi-quantitative RT-PCR. The data showed that these genes displayed different expression levels in response to various stress conditions. The results presented in this study provide basic information about Medicago CCCH Znf genes and form a fundamental clue for cloning genes with specific functions in further studies and applications.

Phylogenetic, expression, and bioinformatic analysis of the ABC1 gene family in Populus trichocarpa

We studied 17 ABC1 genes in Populus trichocarpa, all of which contained an ABC1 domain consisting of about 120 amino acid residues. Most of the ABC1 gene products were located in the mitochondria or chloroplasts. All had a conserved VAVK-like motif and a DFG motif. Phylogenetic analysis grouped the genes into three subgroups. In addition, the chromosomal locations of the genes on the 19 Populus chromosomes were determined. Gene structure was studied through exon/intron organization and the MEME motif finder, while heatmap was used to study the expression diversity using EST libraries. According to the heatmap, PtrABC1P14 was highlighted because of the high expression in tension wood which related to secondary cell wall formation and cellulose synthesis, thus making a contribution to follow-up experiment in wood formation. Promoter cis-element analysis indicated that almost all of the ABC1 genes contained one or two cis-elements related to ABA signal transduction pathway and drought stress. Quantitative real-time PCR was carried out to evaluate the expression of all of the genes under abiotic stress conditions (ABA, CdCl2, high temperature, high salinity, and drought); the results showed that some of the genes were affected by these stresses and confirmed the results of promoter cis-element analysis.

The poplar MYB master switches bind to the SMRE site and activate the secondary wall biosynthetic program during wood formation

Wood is mainly composed of secondary walls, which constitute the most abundant stored carbon produced by vascular plants. Understanding the molecular mechanisms controlling secondary wall deposition during wood formation is not only an important issue in plant biology but also critical for providing molecular tools to custom-design wood composition suited for diverse end uses. Past molecular and genetic studies have revealed a transcriptional network encompassing a group of wood-associated NAC and MYB transcription factors that are involved in the regulation of the secondary wall biosynthetic program during wood formation in poplar trees. Here, we report the functional characterization of poplar orthologs of MYB46 and MYB83 that are known to be master switches of secondary wall biosynthesis in Arabidopsis. In addition to the two previously-described PtrMYB3 and PtrMYB20, two other MYBs, PtrMYB2 and PtrMYB21, were shown to be MYB46/MYB83 orthologs by complementation and overexpression studies in Arabidopsis. The functional roles of these PtrMYBs in regulating secondary wall biosynthesis were further demonstrated in transgenic poplar plants showing an ectopic deposition of secondary walls in PtrMYB overexpressors and a reduction of secondary wall thickening in their dominant repressors. Furthermore, PtrMYB2/3/20/21 together with two other tree MYBs, the Eucalyptus EgMYB2 and the pine PtMYB4, were shown to differentially bind to and activate the eight variants of the 7-bp SMRE consensus sequence, composed of ACC(A/T)A(A/C)(T/C). Together, our results indicate that the tree MYBs, PtrMYB2/3/20/21, EgMYB2 and PtMYB4, are master transcriptional switches that activate the SMRE sites in the promoters of target genes and thereby regulate secondary wall biosynthesis during wood formation.

Genome wide analysis of the apple MYB transcription factor family allows the identification of MdoMYB121 gene confering abiotic stress tolerance in plants

The MYB proteins comprise one of the largest families of transcription factors (TFs) in plants. Although several MYB genes have been characterized to play roles in secondary metabolism, the MYB family has not yet been identified in apple. In this study, 229 apple MYB genes were identified through a genome-wide analysis and divided into 45 subgroups. A computational analysis was conducted using the apple genomic database to yield a complete overview of the MYB family, including the intron-exon organizations, the sequence features of the MYB DNA-binding domains, the carboxy-terminal motifs, and the chromosomal locations. Subsequently, the expression of 18 MYB genes, including 12 were chosen from stress-related subgroups, while another 6 ones from other subgroups, in response to various abiotic stresses was examined. It was found that several of these MYB genes, particularly MdoMYB121, were induced by multiple stresses. The MdoMYB121 was then further functionally characterized. Its predicted protein was found to be localized in the nucleus. A transgenic analysis indicated that the overexpression of the MdoMYB121 gene remarkably enhanced the tolerance to high salinity, drought, and cold stresses in transgenic tomato and apple plants. Our results indicate that the MYB genes are highly conserved in plant species and that MdoMYB121 can be used as a target gene in genetic engineering approaches to improve the tolerance of plants to multiple abiotic stresses.

MYB transcription factor genes as regulators for plant responses: an overview

Regulation of gene expression at the level of transcription controls many crucial biological processes. Transcription factors (TFs) play a great role in controlling cellular processes and MYB TF family is large and involved in controlling various processes like responses to biotic and abiotic stresses, development, differentiation, metabolism, defense etc. Here, we review MYB TFs with particular emphasis on their role in controlling different biological processes. This will provide valuable insights in understanding regulatory networks and associated functions to develop strategies for crop improvement.

Genome-wide identification, characterization and expression analysis of populus leucine-rich repeat receptor-like protein kinase genes

Background

Leucine-rich repeat receptor-like kinases (LRR-RLKs) comprise the largest group within the receptor-like kinase (RLK) superfamily in plants. This gene family plays critical and diverse roles in plant growth, development and stress response. Although the LRR-RLK families in Arabidopsis and rice have been previously analyzed, no comprehensive studies have been performed on this gene family in tree species.

Results

In this work, 379 LRR-RLK genes were retrieved from the Populus trichocarpa genome and further grouped into 14 subfamilies based on their structural and sequence similarities. Approximately 82% (312 out of 379) of the PtLRR-RLK genes are located in segmental duplication blocks indicating the role of duplication process in the expansion of this gene family. The conservation and variation in motif composition and intron/exon arrangement among PtLRR-RLK subfamilies were analyzed to provide additional support for their phylogenetic relationship and more importantly to indicate the potential divergence in their functions. Expression profiling of PtLRR-RLKs showed that they were differentially expressed in different organs and tissues and some PtLRR-RLKs were specifically expressed in meristem tissues, which indicated their potential involvement in tissue development and differentiation. For most AtLRR-RLKs with defined functions, Populus homologues exhibiting similar expression patterns could be identified, which might indicate the functional conservation during evolution. Among 12 types of environmental cues analyzed by the genome-wide microarray data, PtLRR-RLKs showed specific responses to shoot organogenesis, wounding, low ammonium feeding, hypoxia and seasonal dormancy, but not to drought, re-watering after drought, flooding, AlCl₃ treatment and bacteria or fungi treatments.

Conclusions

This study provides the first comprehensive genomic analysis of the Populus LRR-RLK gene family. Segmental duplication contributes significantly to the expansion of this gene family. Populus and Arabidopsis LRR-RLK homologues not only share similar genetic structures but also exhibit comparable expression patterns which point to the possible functional conservation of these LRR-RLKs in two model systems. Transcriptome profiling provides the first insight into the functional divergence among PtLRR-RLK gene subfamilies and suggests that they might take important roles in growth and adaptation of tree species.

Transpirational demand affects aquaporin expression in poplar roots

Isohydric plants tend to maintain a water potential homeostasis primarily by controlling water loss via stomatal conductance. However, there is accumulating evidence that plants can also modulate water uptake in a dynamic manner. The dynamics of water uptake are influenced by aquaporin-mediated changes in root hydraulics. Most studies in this area have been conducted on herbaceous plants, and less is known about responses of woody plants. Here a study was conducted to determine how roots of hybrid poplar plants (Populus trichocarpa×deltoides) respond to a step change in transpirational demand. The main objective was to measure the expression of selected aquaporin genes and to assess how transcriptional responses correspond to changes in root water flow (Q R) and other parameters of water relations. A subset of plants was grown in shade and was subsequently exposed to a 5-fold increase in light level. Another group of plants was grown at ~95% relative humidity (RH) and was then subjected to lower RH while the light level remained unchanged. Both plant groups experienced a transient drop in stem water potentials. At 28h after the increase in transpirational demand, water potentials recovered. This recovery was associated with changes in the expression of PIP1 and PIP2 subfamily genes and an increase in Q R. Stomata of plants growing at high RH were larger and showed incomplete closure after application of abscisic acid. Since stomatal conductance remained high and unchanged in these plants, it is suggested that the recovery in water potential in these plants was largely driven by the increase in Q R.

Breeding with rare defective alleles (BRDA): a natural Populus nigra HCT mutant with modified lignin as a case study

Next-generation (NG) sequencing in a natural population of Populus nigra revealed a mutant with a premature stop codon in the gene encoding hydroxycinnamoyl-CoA : shikimate hydroxycinnamoyl transferase1 (HCT1), an essential enzyme in lignin biosynthesis. The lignin composition of P. nigra trees homozygous for the defective allele was compared with that of heterozygous trees and trees without the defective allele. The lignin was characterized by phenolic profiling, lignin oligomer sequencing, thioacidolysis and NMR. In addition, HCT1 was heterologously expressed for activity assays and crosses were made to introduce the mutation in different genetic backgrounds. HCT1 converts p-coumaroyl-CoA into p-coumaroyl shikimate. The mutant allele, PnHCT1-Δ73, encodes a truncated protein, and trees homozygous for this recessive allele have a modified lignin composition characterized by a 17-fold increase in p-hydroxyphenyl units. Using the lignin pathway as proof of concept, we illustrated that the capture of rare defective alleles is a straightforward approach to initiate reverse genetics and accelerate tree breeding. The proposed breeding strategy, called 'breeding with rare defective alleles' (BRDA), should be widely applicable, independent of the target gene or the species.

Prevalent role of gene features in determining evolutionary fates of whole-genome duplication duplicated genes in flowering plants

The evolution of genes and genomes after polyploidization has been the subject of extensive studies in evolutionary biology and plant sciences. While a significant number of duplicated genes are rapidly removed during a process called fractionation, which operates after the whole-genome duplication (WGD), another considerable number of genes are retained preferentially, leading to the phenomenon of biased gene retention. However, the evolutionary mechanisms underlying gene retention after WGD remain largely unknown. Through genome-wide analyses of sequence and functional data, we comprehensively investigated the relationships between gene features and the retention probability of duplicated genes after WGDs in six plant genomes, Arabidopsis (Arabidopsis thaliana), poplar (Populus trichocarpa), soybean (Glycine max), rice (Oryza sativa), sorghum (Sorghum bicolor), and maize (Zea mays). The results showed that multiple gene features were correlated with the probability of gene retention. Using a logistic regression model based on principal component analysis, we resolved evolutionary rate, structural complexity, and GC3 content as the three major contributors to gene retention. Cluster analysis of these features further classified retained genes into three distinct groups in terms of gene features and evolutionary behaviors. Type I genes are more prone to be selected by dosage balance; type II genes are possibly subject to subfunctionalization; and type III genes may serve as potential targets for neofunctionalization. This study highlights that gene features are able to act jointly as primary forces when determining the retention and evolution of WGD-derived duplicated genes in flowering plants. These findings thus may help to provide a resolution to the debate on different evolutionary models of gene fates after WGDs.

Transcriptome-wide identification of R2R3-MYB transcription factors in barley with their boron responsive expression analysis

MYB family of transcription factors (TF) comprises one of the largest transcription factors in plants and is represented in all eukaryotes. They include highly conserved MYB repeats (1R, R2R3, 3R, and 4R) in the N-terminus. In addition to this, they have diverse C-terminal sequences which help the protein gain wide distinct functions, such as controlling development, secondary metabolism, hormonal regulation and response to biotic and abiotic stress. Stress-responsive roles of the MYB TFs were reported for drought, salt, wounding, cold, freezing, dehydration and osmotic stresses. This study describes the identification of barley R2R3-MYB TFs including their expression analysis in tissues under control and Boron (B) toxic conditions. Conserved motifs for MYB proteins were searched into barley full-transcriptome RNA-seq data and a total of 320 protein sequences were filtered as MYB TFs in which 51 of them corresponded to R2R3 MYB TFs. Using various bioinformatics tools, their conserved domain structures, chromosomal distributions, gene duplications, comparative functional analysis, as well as phylogenetic relations with Arabidopsis thaliana, were conducted. Beside the RNA-seq data-based expression pattern analysis of 51 R2R3 MYB TFs, quantitative analysis of selected R2R3 MYB TF genes was assessed in control and B-stressed root and leaf tissues. Critical B-induced R2R3 MYB TFs were identified. It was concluded that the results would be useful for functional characterizations of R2R3-type MYB transcription factors that are possibly involved in both B stress and divergent regulation mechanisms in plants.

Effect of nitric oxide on gene transcription - S-nitrosylation of nuclear proteins

Nitric oxide (NO) plays an important role in many different physiological processes in plants. It mainly acts by post-translationally modifying proteins. Modification of cysteine residues termed as S-nitrosylation is believed to be the most important mechanism for transduction of bioactivity of NO. The first proteins found to be nitrosylated were mainly of cytoplasmic origin or isolated from mitochondria and peroxisomes. Interestingly, it was shown that redox-sensitive transcription factors are also nitrosylated and that NO influences the redox-dependent nuclear transport of some proteins. This implies that NO plays a role in regulating transcription and/or general nuclear metabolism which is a fascinating new aspect of NO signaling in plants. In this review, we will discuss the impact of S-nitrosylation on nuclear plant proteins with a focus on transcriptional regulation, describe the function of this modification and draw also comparisons to the animal system in which S-nitrosylation of nuclear proteins is a well characterized concept.

Isolation and molecular characterization of thirteen R2R3-MYB transcription factors from Epimedium sagittatum

Epimedium sagittatum (Sieb. et Zucc.) Maxim, a popular traditional Chinese medicinal plant, has been widely used for treating sexual dysfunction and osteoporosis in China. The main bioactive components in herba epimedii are prenylated flavonol glycosides, which are end products of a branch of the flavonoid biosynthetic pathway. The MYB transcription factors (TF) act as activators or repressors to regulate the flavonoid pathway. In this study, 13 full-length cDNA clones of R2R3-MYB TFs from E. sagittatum (designated as EsMYB1 to EsMYB13) were isolated and characterized. Sequence similarity and phylogenetic analysis placed nine R2R3-MYB members of epimedii into five subgroups of the Arabidopsis R2R3-MYB family, while four members were not clustered into a defined subgroup. The number and length of introns from epimedii R2R3-MYB genes varied significantly, but intron positions and phases were well conserved. Expression patterns of epimedii R2R3-MYB genes in various tissues showed diverse. Finally, it is suggested that five epimedii R2R3-MYB genes may be involved in regulating the flavonoid pathway and could be used as valuable candidate genes for metabolic engineering studies in future. Sequence information of 13 R2R3-MYB genes discovered here will also provide an entry point into the overview of whole R2R3-MYB family in epimedii.

Genome-wide identification, classification, and expression analysis of CDPK and its closely related gene families in poplar (Populus trichocarpa)

Calcium-dependent protein kinases (CDPKs) are Ca(2+)-binding proteins known to play crucial roles in Ca(2+) signal transduction pathways which have been identified throughout plant kingdom and in certain types of protists. Genome-wide analysis of CDPKs have been carried out in Arabidopsis, rice and wheat, and quite a few of CDPKs were proved to play crucial roles in plant stress responsive signature pathways. In this study, a comprehensive analysis of Populus CDPK and its closely related gene families was performed, including phylogeny, chromosome locations, gene structures, and expression profiles. Thirty Populus CDPK genes and twenty closely related kinase genes were identified, which were phylogenetically clustered into eight distinct subfamilies and predominately distributed across fifteen linkage groups (LG). Genomic organization analyses indicated that purifying selection has played a pivotal role in the retention and maintenance of Populus CDPK gene family. Furthermore, microarray analysis showed that a number of Populus CDPK and its closely related genes differentially expressed across disparate tissues and under various stresses. The expression profiles of paralogous pairs were also investigated to reveal their evolution fates. In addition, quantitative real-time RT-PCR was performed on nine selected CDPK genes to confirm their responses to drought stress treatment. These observations may lay the foundation for future functional analysis of Populus CDPK and its closely related gene families to unravel their biological roles.

The R2R3-MYB-like regulatory factor EOBI, acting downstream of EOBII, regulates scent production by activating ODO1 and structural scent-related genes in petunia

Flower scent is a highly dynamic trait, under developmental, spatial, and diurnal regulation. The mechanism governing scent production is only beginning to be unraveled. In petunia (Petunia hybrida), EMISSION OF BENZENOIDS II (EOBII) controls transcription of both the shikimate pathway-regulating MYB factor ODORANT1 (ODO1) and phenylpropanoid scent-related structural genes. A promoter-activation screen identified an R2R3-MYB-like regulatory factor of phenylpropanoid volatile biosynthesis acting downstream of EOBII, designated EOBI. EOBI silencing led to downregulation of ODO1 and numerous structural scent-related genes from both the shikimate and phenylpropanoid pathways. The ability of EOBI to directly activate ODO1, as revealed by electrophoretic mobility shift assay and yeast one-hybrid analysis, place EOBI upstream of ODO1 in regulating substrate availability for volatile biosynthesis. Interestingly, ODO1-silenced transgenic petunia flowers accumulated higher EOBI transcript levels than controls, suggesting a complex feedback loop between these regulatory factors. The accumulation pattern of EOBI transcript relative to EOBII and ODO1, and the effect of up/downregulation of EOBII on transcript levels of EOBI and ODO1, further support these factors' hierarchical relationships. The dependence of scent production on EOBI expression and its direct interaction with both regulatory and structural genes provide evidence for EOBI's wide-ranging involvement in the production of floral volatiles.

Evolutionary analysis of the MIXTA gene family highlights potential targets for the study of cellular differentiation

Differentiated epidermal cells such as trichomes and conical cells perform numerous essential functions in plant biology and are important for our understanding of developmental patterning and cell shape regulation. Many are also commercially significant, such as cotton fibers and trichomes that secrete pharmaceutically useful or herbivore-deterring compounds. Here, we focus on the phylogeny and evolution of the subgroup 9 R2R3 MYB gene transcription factors, which include the MIXTA gene, and that are important for the specification and regulation of plant cellular differentiation. We have sequenced 49 subgroup 9 R2R3 MYB genes from key experimental taxa and combined these sequences with those identified by an exhaustive bioinformatic search, to compile a data set of 223 subgroup 9 R2R3 MYB genes. Our phylogenetic analyses demonstrate, for the first time, the complex evolutionary history of the subgroup 9 R2R3 MYB genes. A duplication event is inferred before the origin of seed plants giving rise to two major gene lineages, here termed SBG9-A and SBG9-B. The evolutionary conservation of the SBG9-B gene lineage has not been previously recognized and its role in cellular differentiation is unknown, thus an entire clade of potential candidate genes for epidermal cell regulation remains to be explored. Using a heterologous transformation bioassay, we provide functional data that implicate members of the SBG9-B lineage in the specification of epidermal projections. Furthermore, we reveal numerous putative duplication events in both SBG9-A and SBG9-B lineages, resolving uncertainty about orthology and paralogy among the subgroup 9 R2R3 MYB genes. Finally, we provide a robust framework over which to interpret existing functional data and to direct ongoing comparative genetic research into the evolution of plant cellular diversity.

Genome-wide identification and characterization of R2R3MYB family in Cucumis sativus

BACKGROUND

The R2R3MYB proteins comprise one of the largest families of transcription factors in plants. Although genome-wide analysis of this family has been carried out in some species, little is known about R2R3MYB genes in cucumber (Cucumis sativus L.).

PRINCIPAL FINDINGS

This study has identified 55 R2R3MYB genes in the latest cucumber genome and the CsR2R3MYB family contained the smallest number of identified genes compared to other species that have been studied due to the absence of recent gene duplication events. These results were also supported by genome distribution and gene duplication analysis. Phylogenetic analysis showed that they could be classified into 11 subgroups. The evolutionary relationships and the intron-exon organizations that showed similarities with Arabidopsis, Vitis and Glycine R2R3MYB proteins were also analyzed and suggested strong gene conservation but also the expansions of particular functional genes during the evolution of the plant species. In addition, we found that 8 out of 55 (∼14.54%) cucumber R2R3MYB genes underwent alternative splicing events, producing a variety of transcripts from a single gene, which illustrated the extremely high complexity of transcriptome regulation. Tissue-specific expression profiles showed that 50 cucumber R2R3MYB genes were expressed in at least one of the tissues and the other 5 genes showed very low expression in all tissues tested, which suggested that cucumber R2R3MYB genes took part in many cellular processes. The transcript abundance level analysis during abiotic conditions (NaCl, ABA and low temperature treatments) identified a group of R2R3MYB genes that responded to one or more treatments.

CONCLUSIONS

This study has produced a comparative genomics analysis of the cucumber R2R3MYB gene family and has provided the first steps towards the selection of CsR2R3MYB genes for cloning and functional dissection that can be used in further studies to uncover their roles in cucumber growth and development.

From landing lights to mimicry: the molecular regulation of flower colouration and mechanisms for pigmentation patterning

Flower colour is a key component for plant signaling to pollinators and a staggering variety of colour variations are found in nature. Patterning of flower colour, such as pigment spots or stripes, is common and is important in promoting pollination success. Developmentally programmed pigmentation patterns are of interest with respect to the evolution of specialised plant-pollinator associations and as models for dissecting regulatory signaling in plants. This article reviews the occurrence and function of flower colour patterns, as well as the molecular genetics of anthocyanin pigmentation regulation. The transcription factors controlling anthocyanin biosynthesis have been characterised for many species and an 'MBW' regulatory complex of R2R3MYB, bHLH and WD-Repeat proteins is of central importance. In particular, R2R3MYBs are key determinants of pigmentation intensity and patterning in plants. Progress is now being made on how environmental or developmental signal pathways may in turn control the production of the MBW components. Furthermore, additional regulatory proteins that interact with the MBW activation complex are being identified, including a range of proteins that repress complex formation or action, either directly or indirectly. This review discusses some of the recent data on the regulatory factors and presents models of how patterns may be determined.