A new role for plant R2R3-MYB transcription factors in cell cycle regulation

Genome-wide identification and characterization of flowering genes in Citrus sinensis (L.) Osbeck: a comparison among C. Medica L., C. Reticulata Blanco, C. Grandis (L.) Osbeck and C. Clementina

BACKGROUND

Flowering plays an important role in completing the reproductive cycle of plants and obtaining next generation of plants. In case of citrus, it may take more than a year to achieve progeny. Therefore, in order to fasten the breeding processes, the juvenility period needs to be reduced. The juvenility in plants is regulated by set of various flowering genes. The citrus fruit and leaves possess various medicinal properties and are subjected to intensive breeding programs to produce hybrids with improved quality traits. In order to break juvenility in Citrus, it is important to study the role of flowering genes. The present study involved identification of genes regulating flowering in Citrus sinensis L. Osbeck via homology based approach. The structural and functional characterization of these genes would help in targeting genome editing techniques to induce mutations in these genes for producing desirable results.

RESULTS

A total of 43 genes were identified which were located on all the 9 chromosomes of citrus. The in-silico analysis was performed to determine the genetic structure, conserved motifs, cis-regulatory elements (CREs) and phylogenetic relationship of the genes. A total of 10 CREs responsible for flowering were detected in 33 genes and 8 conserved motifs were identified in all the genes. The protein structure, protein-protein interaction network and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was performed to study the functioning of these genes which revealed the involvement of flowering proteins in circadian rhythm pathways. The gene ontology (GO) and gene function analysis was performed to functionally annotate the genes. The structure of the genes and proteins were also compared among other Citrus species to study the evolutionary relationship among them. The expression study revealed the expression of flowering genes in floral buds and ovaries. The qRT-PCR analysis revealed that the flowering genes were highly expressed in bud stage, fully grown flower and early stage of fruit development.

CONCLUSIONS

The findings suggested that the flowering genes were highly conserved in citrus species. The qRT-PCR analysis revealed the tissue specific expression of flowering genes (CsFT, CsCO, CsSOC, CsAP, CsSEP and CsLFY) which would help in easy detection and targeting of genes through various forward and reverse genetic approaches.

Genome-Wide Identification and Analysis of MYB Transcription Factors in Pyropia yezoensis

MYB transcription factors are one of the largest transcription factor families in plants, and they regulate numerous biological processes. Red algae are an important taxonomic group and have important roles in economics and research. However, no comprehensive analysis of the MYB gene family in any red algae, including Pyropia yezoensis, has been conducted. To identify the MYB gene members of Py. yezoensis, and to investigate their family structural features and expression profile characteristics, a study was conducted. In this study, 3 R2R3-MYBs and 13 MYB-related members were identified in Py. yezoensis. Phylogenetic analysis indicated that most red algae MYB genes could be clustered with green plants or Glaucophyta MYB genes, inferring their ancient origins. Synteny analysis indicated that 13 and 5 PyMYB genes were orthologous to Pyropia haitanensis and Porphyra umbilicalis, respectively. Most Bangiaceae MYB genes contain several Gly-rich motifs, which may be the result of an adaptation to carbon limitations and maintenance of important regulatory functions. An expression profile analysis showed that PyMYB genes exhibited diverse expression profiles. However, the expression patterns of different members appeared to be diverse, and PyMYB5 was upregulated in response to dehydration, low temperature, and Pythium porphyrae infection. This is the first comprehensive study of the MYB gene family in Py. Yezoensis and it provides vital insights into the functional divergence of MYB genes.

MIR159 regulates multiple aspects of stamen and carpel development and requires dissection and delimitation of differential downstream regulatory network for manipulating fertility traits

Unravelling genetic networks regulating developmental programs are key to devising and implementing genomics assisted trait modification strategies. It is crucial to understand the role of small RNAs, and the basis of their ability to modify traits. MIR159 has been previously reported to cause defects in anther development in Arabidopsis; however, the complete spectrum and basis of the defects remained unclear. The present study was therefore undertaken to comprehensively investigate the role of miR159 from Brassica juncea in modulating vegetative and reproductive traits. Owing to the polyploid nature of Brassica, paralogous and homeologous copies of MIR159A, MIR159B, and, MIR159C were identified and analysis of the precursor uncovered extensive structural and sequence variation. The MIR159 locus with mature miR159 with perfect target complimentarily with MYB65, was cloned from Brassica juncea var. Varuna for functional characterization by generating constitutively over-expressing lines in Arabidopsis thaliana Col-0. Apart from statistically significant difference in multiple vegetative traits, drastic differences were observed in stamen and pistil. Over-expression of miR159a led to shortening of filament length and loss of tetradynamous condition. Anthers were apiculate, with improper lobe formation, and unsynchronized cellular growth between connective tissue and another lobe development. Analysis revealed arrested meiosis/cytokinesis in microspores, and altered lignin deposition pattern in endothecial walls thus affecting anther dehiscence. In the gynoecium, flaccid, dry stigmatic papillae, and large embryo sac in the female gametophyte was observed. Over-expression of miR159a thus severely affected pollination and seed-set. Analysis of the transcriptome data revealed components of regulatory networks of anther and carpel developmental pathway, and lignin metabolism that are affected. Expression analysis allowed us to position the miR159a-MYB65 module in the genetic network of stamen development, involved in pollen-grain maturation; in GA-mediated regulation of stamen development, and in lignin metabolism. The study, on one hand indicates role of miR159a-MYB65 in regulating multiple aspects of reproductive organ development that can be manipulated for trait modification, but also raises several unaddressed questions such as relationship between miR159a and male-meiosis, miR159a and filament elongation for future investigations. Accession numbers: KC204951-KC204960. Project number PRJNA1035268.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01377-7.

The Over-Expression of Two R2R3-MYB Genes, PdMYB2R089 and PdMYB2R151, Increases the Drought-Resistant Capacity of Transgenic Arabidopsis

The R2R3-MYB genes in plants play an essential role in the drought-responsive signaling pathway. Plenty of R2R3-MYB S21 and S22 subgroup genes in Arabidopsis have been implicated in dehydration conditions, yet few have been covered in terms of the role of the S21 and S22 subgroup genes in poplar under drought. PdMYB2R089 and PdMYB2R151 genes, respectively belonging to the S21 and S22 subgroups of NL895 (Populus deltoides × P. euramericana cv. 'Nanlin895'), were selected based on the previous expression analysis of poplar R2R3-MYB genes that are responsive to dehydration. The regulatory functions of two target genes in plant responses to drought stress were studied and speculated through the genetic transformation of Arabidopsis thaliana. PdMYB2R089 and PdMYB2R151 could promote the closure of stomata in leaves, lessen the production of malondialdehyde (MDA), enhance the activity of the peroxidase (POD) enzyme, and shorten the life cycle of transgenic plants, in part owing to their similar conserved domains. Moreover, PdMYB2R089 could strengthen root length and lateral root growth. These results suggest that PdMYB2R089 and PdMYB2R151 genes might have the potential to improve drought adaptability in plants. In addition, PdMYB2R151 could significantly improve the seed germination rate of transgenic Arabidopsis, but PdMYB2R089 could not. This finding provides a clue for the subsequent functional dissection of S21 and S22 subgroup genes in poplar that is responsive to drought.

Regulation of the regulators: Transcription factors controlling biosynthesis of plant secondary metabolites during biotic stresses and their regulation by miRNAs

Biotic stresses threaten to destabilize global food security and cause major losses to crop yield worldwide. In response to pest and pathogen attacks, plants trigger many adaptive cellular, morphological, physiological, and metabolic changes. One of the crucial stress-induced adaptive responses is the synthesis and accumulation of plant secondary metabolites (PSMs). PSMs mitigate the adverse effects of stress by maintaining the normal physiological and metabolic functioning of the plants, thereby providing stress tolerance. This differential production of PSMs is tightly orchestrated by master regulatory elements, Transcription factors (TFs) express differentially or undergo transcriptional and translational modifications during stress conditions and influence the production of PSMs. Amongst others, microRNAs, a class of small, non-coding RNA molecules that regulate gene expression post-transcriptionally, also play a vital role in controlling the expression of many such TFs. The present review summarizes the role of stress-inducible TFs in synthesizing and accumulating secondary metabolites and also highlights how miRNAs fine-tune the differential expression of various stress-responsive transcription factors during biotic stress.

Genome-wide identification and comprehensive analysis of tubby-like protein gene family in multiple crops

INTRODUCTION

The highly conserved tubby-like proteins (TLPs) play key roles in animal neuronal development and plant growth. The abiotic stress tolerance function of TLPs has been widely explored in plants, however, little is known about comparative studies of TLPs within crops.

METHODS

Bioinformatic identification, phylogenetic analysis, Cis-element analysis, expression analysis, Cis-element analysis, expression analysis and so on were explored to analysis the TLP gene family of multiple crops.

RESULTS

In this study, a comprehensive analysis of TLP genes were carried out in seven crops to explore whether similar function of TLPs in rice could be achieved in other crops. We identified 20, 9, 14, 11, 12, 35, 14 and 13 TLP genes in Glycine max, Hordeum vulgare, Sorghum bicolor, Arabidopsis thaliana, Oryza sativa Japonica, Triticum aestivum, Setaria italic and Zea mays, respectively. All of them were divided into two groups and ten orthogroups (Ors) based on amino acids. A majority of TLP genes had two domains, tubby-like domain and F-box domain, while members of Or5 only had tubby-like domain. In addition, Or5 had more exons and shorter DNA sequences, showing that characteristics of different Ors reflected the differentiated function and feature of TLP genes in evolutionary process, and Or5 was the most different from the other Ors. Besides, we recognized 25 cis-elements in the promoter of TLP genes and explored multiple new regulation pathway of TLPs including light and hormone response. The bioinformatic and transcriptomic analysis implied the stresses induced expression and possible functional redundancy of TLP genes. We detected the expression level of 6 OsTLP genes at 1 to 6 days after seed germination in rice, and the most obvious changes in these days were appeared in OsTLP10 and OsTLP12.

DISCUSSION

Combined yeast two-hybrid system and pull down assay, we suggested that the TLP genes of Or1 may have similar function during seed germination in different species. In general, the results of comprehensive analysis of TLP gene family in multiple species provide valuable evolutionary and functional information of TLP gene family which are useful for further application and study of TLP genes.

Global Integrated Genomic and Transcriptomic Analyses of MYB Transcription Factor Superfamily in C3 Model Plant Oryza sativa (L.) Unravel Potential Candidates Involved in Abiotic Stress Signaling

Plant transcription factors (TFs) are significant players in transcriptional regulations, signal transduction, and constitute an integral part of signaling networks. MYB TFs are major TF superfamilies that play pivotal roles in regulation of transcriptional reprogramming, physiological processes, and abiotic stress (AbS) responses. To explore the understanding of MYB TFs, genome and transcriptome-wide identification was performed in the C3 model plant, Oryza sativa (OsMYB). This study retrieved 114 OsMYB TFs that were computationally analyzed for their expression profiling, gene organization, cis-acting elements, and physicochemical properties. Based on the microarray datasets, six OsMYB genes which were sorted out and identified by a differential expression pattern were noted in various tissues. Systematic expression profiling of OsMYB TFs showed their meta-differential expression of different AbS treatments, spatio-temporal gene expression of various tissues and their growth in the field, and gene expression profiling in responses to phytohormones. In addition, the circular ideogram of OsMYB genes in related C4 grass plants conferred the gene synteny. Protein–protein interactions of these genes revealed the molecular crosstalk of OsMYB TFs. Transcriptional analysis (qPCR) of six OsMYB players in response to drought and salinity stress suggested the involvement in individual and combined AbS responses. To decipher how these OsMYB play functional roles in AbS dynamics, further research is a prerequisite.

Comprehensive analysis of the MYB transcription factor gene family in Morus alba

BACKGROUND

The V-myb myeloblastosis viral oncogene homolog (MYB) family of proteins is large, containing functionally diverse transcription factors. However, MYBs in Morus are still poorly annotated and a comprehensive functional analysis of these transcription factors is lacking.

RESULTS

In the present study, a genome-wide identification of MYBs in Morus alba was performed. In total 166 MaMYBs were identified, including 103 R2R3-MYBs and four 3R-MaMYBs. Comprehensive analyses, including the phylogenetic analysis with putative functional annotation, motif and structure analysis, gene structure organization, promoter analysis, chromosomal localization, and syntenic relationships of R2R3-MaMYBs and 3R-MaMYBs, provided primary characterization for these MaMYBs. R2R3-MaMYBs covered the subgroups reported for R2R3-MYBs in Arabidopsis and Populus, and had two Morus-specific subgroups, indicating the high retention of MYBs in Morus. Motif analysis revealed high conservative residues at the start and end of each helix and residues consisting of the third helix in R2 and R3 repeats. Thirteen intron/exon patterns (a-m) were summarized, and the intron/exon pattern of two introns with phase numbers of 0 and 2 was the prevalent pattern for R2R3-MaMYBs. Various cis-elements in promoter regions were identified, and were mainly related to light response, development, phytohormone response, and abiotic and biotic stress response and secondary metabolite production. Expression patterns of R2R3-MaMYBs in different organs showed that MaMYBs involved in secondary cell wall components and stress responsiveness were preferentially expressed in roots or stems. R2R3-MaMYBs involved in flavonoid biosynthesis and anthocyanin accumulation were identified and characterized based on functional annotation and correlation of their expression levels with anthocyanin contents.

CONCLUSION

Based on a comprehensive analysis, this work provided functional annotation for R2R3-MYBs and an informative reference for further functional dissection of MYBs in Morus.

MYB transcription factor family in sweet cherry (Prunus avium L.): genome-wide investigation, evolution, structure, characterization and expression patterns

BACK GROUND

MYB Transcription factors (TFs) are most imperative and largest gene family in plants, which participate in development, metabolism, defense, differentiation and stress response. The MYB TFs has been studied in various plant species. However, comprehensive studies of MYB gene family in the sweet cherry (Prunus avium L.) are still unknown.

RESULTS

In the current study, a total of 69 MYB genes were investigated from sweet cherry genome and classified into 28 subfamilies (C1-C28 based on phylogenetic and structural analysis). Microcollinearity analysis revealed that dispersed duplication (DSD) events might play an important role in the MYB genes family expansion. Chromosomal localization, the synonymous (Ks) and nonsynonymous (Ka) analysis, molecular characteristics (pI, weight and length of amino acids) and subcellular localization were accomplished using several bioinformatics tools. Furthermore, the members of distinct subfamilies have diverse cis-acting regions, conserved motifs, and intron-exon architectures, indicating functional heterogeneity in the MYB family. Moreover, the transcriptomic data exposed that MYB genes might play vital role in bud dormancy. The quantitative real-time qRT-PCR was carried out and the expression pattern indicated that MYB genes significantly expressed in floral bud as compared to flower and fruit.

CONCLUSION

Our comprehensive findings provide supportive insights into the evolutions, expansion complexity and functionality of PavMYB genes. These PavMYB genes should be further investigated as they seem to be brilliant candidates for dormancy manipulation in sweet cherry.

R2R3-MYBs in Durum Wheat: Genome-Wide Identification, Poaceae-Specific Clusters, Expression, and Regulatory Dynamics Under Abiotic Stresses

MYB transcription factors (TFs) represent one of the biggest TF families in plants, being involved in various specific plant processes, such as responses to biotic and abiotic stresses. The implication of MYB TFs in the tolerance mechanisms to abiotic stress is particularly interesting for crop breeding, since environmental conditions can negatively affect growth and productivity. Wheat is a worldwide-cultivated cereal, and is a major source of plant-based proteins in human food. In particular, durum wheat plays an important role in global food security improvement, since its adaptation to hot and dry conditions constitutes the base for the success of wheat breeding programs in future. In the present study, a genome-wide identification of R2R3-MYB TFs in durum wheat was performed. MYB profile search and phylogenetic analyses based on homology with Arabidopsis and rice MYB TFs led to the identification of 233 R2R3-TdMYB (Triticum durum MYB). Three Poaceae-specific MYB clusters were detected, one of which had never been described before. The expression of eight selected genes under different abiotic stress conditions, revealed that most of them responded especially to salt and drought stress. Finally, gene regulatory network analyses led to the identification of 41 gene targets for three TdR2R3-MYBs that represent novel candidates for functional analyses. This study provides a detailed description of durum wheat R2R3-MYB genes and contributes to a deeper understanding of the molecular response of durum wheat to unfavorable climate conditions.

Impact of CuO nanoparticles on maize: Comparison with CuO bulk particles with special reference to oxidative stress damages and antioxidant defense status

The present study aimed to systematically investigate the particle size effects of copper (II) oxide [CuO nanoparticles (<50 nm) and CuO bulk particles (<10 μm)] on maize (Zea mays L.). Bioaccumulation of Cu, in vivo ROS generation, membrane damage, transcriptional modulation of antioxidant genes, cellular redox status of glutathione and ascorbate pool, expression patterns of COPPER TRANSPORTER 4 and stress responsive miRNAs (miR398a, miR171b, miR159f-3p) with their targets were investigated for better understanding of the underlying mechanisms and the extent of CuO nanoparticles and CuO bulk particles induced oxidative stress damages. More restricted seedling growth, comparatively higher membrane injury, marked decline in the levels of chlorophylls and carotenoids and severe oxidative burst were evident in CuO bulk particles challenged leaves. Dihydroethidium and CM-H2DCFDA staining further supported elevated reactive oxygen species generation in CuO bulk particles stressed roots. CuO bulk particles exposed seedlings accumulated much higher amount of Cu in roots as compared to CuO nanoparticles stressed plants with low root-to-shoot Cu translocation. Moderately high GR expression with maintenance of a steady GSH-GSSG ratio in CuO nanoparticles challenged leaves might be accountable for their rather improved performance under stressed condition. miR171b-mediated enhanced expression of SCARECROW 6 might participate in the marked decline of chlorophyll content in CuO bulk particles exposed leaves. Ineffective recycling of AsA pool is another decisive feature of inadequate performance of CuO bulk particles stressed seedlings in combating oxidative stress damages. Taken together, our findings revealed that toxicity of CuO bulk particles was higher than CuO nanoparticles and the adverse effects of CuO bulk particles on maize seedlings might be due to higher Cu ions dissolution.

Low Light/Darkness as Stressors of Multifactor-Induced Senescence in Rice Plants

Leaf senescence, as an integral part of the final development stage for plants, primarily remobilizes nutrients from the sources to the sinks in response to different stressors. The premature senescence of leaves is a critical challenge that causes significant economic losses in terms of crop yields. Although low light causes losses of up to 50% and affects rice yield and quality, its regulatory mechanisms remain poorly elucidated. Darkness-mediated premature leaf senescence is a well-studied stressor. It initiates the expression of senescence-associated genes (SAGs), which have been implicated in chlorophyll breakdown and degradation. The molecular and biochemical regulatory mechanisms of premature leaf senescence show significant levels of redundant biomass in complex pathways. Thus, clarifying the regulatory mechanisms of low-light/dark-induced senescence may be conducive to developing strategies for rice crop improvement. This review describes the recent molecular regulatory mechanisms associated with low-light response and dark-induced senescence (DIS), and their effects on plastid signaling and photosynthesis-mediated processes, chloroplast and protein degradation, as well as hormonal and transcriptional regulation in rice.

Identification and characterization of genes related to salt stress tolerance within segregation distortion regions of genetic map in F2 population of upland cotton

Segregation distortion (SD) is a genetic mechanism commonly found in segregating or stable populations. The principle behind this puzzles many researchers. The F₂ generation developed from wild Gossypium darwinii and G. hirsutum CCRI12 species was used to investigate the possible transcription factors within the segregation distortion regions (SDRs). The 384 out of 2763 markers were distorted in 29 SDRs on 18 chromosomes. Good collinearity was observed among genetic and physical maps of G. hirsutum and G. barbadense syntenic blocks. Total 568 genes were identified from SDRs of 18 chromosomes. Out of these genes, 128 belonged to three top-ranked salt-tolerant gene families. The DUF597 contained 8 uncharacterized genes linked to Pkinase (PF00069) gene family in the phylogenetic tree, while 15 uncharacterized genes clustered with the zinc finger gene family. Two hundred thirty four miRNAs targeted numerous genes, including ghr-miR156, ghr-miR399 and ghr-miR482, while others targeted top-ranked stress-responsive transcription factors. Moreover, these genes were involved in the regulation of numerous stress-responsive cis-regulatory elements. The RNA sequence data of fifteen upregulated genes were verified through the RT-qPCR. The expression profiles of two highly upregulated genes (Gh_D01G2015 and Gh_A01G1773) in salt-tolerant G. darwinii showed antagonistic expression in G. hirsutum. The results indicated that salt-tolerant genes have been possibly transferred from the wild G. darwinii species. A detailed functional analysis of these genes can be carried out which might be helpful in the future for gene cloning, transformation, gene editing and the development of salt-resistant cotton varieties.

Transcriptomes analysis reveals novel insight into the molecular mechanisms of somatic embryogenesis in Hevea brasiliensis

BACKGROUND

Somatic embryogenesis (SE) is a promising technology for plant vegetative propagation, which has an important role in tree breeding. Though rubber tree (Hevea brasiliensis Muell. Arg.) SE has been founded, few late SE-related genes have been identified and the molecular regulation mechanisms of late SE are still not well understood.

RESULTS

In this study, the transcriptomes of embryogenic callus (EC), primary embryo (PE), cotyledonary embryo (CE), abnormal embryo (AE), mature cotyledonary embryo (MCE) and withered abnormal embryo (WAE) were analyzed. A total of 887,852,416 clean reads were generated, 85.92% of them were mapped to the rubber tree genome. The de novo assembly generated 36,937 unigenes. The differentially expressed genes (DEGs) were identified in the pairwise comparisons of CE vs. AE and MCE vs. WAE, respectively. The specific common DEGs were mainly involved in the phytohormones signaling pathway, biosynthesis of phenylpropanoid and starch and sucrose metabolism. Among them, hormone signal transduction related genes were significantly enriched, especially the auxin signaling factors (AUX-like1, GH3.1, SAUR32-like, IAA9-like, IAA14-like, IAA27-like, IAA28-like and ARF5-like). The transcription factors including WRKY40, WRKY70, MYBS3-like, MYB1R1-like, AIL6 and bHLH93-like were characterized as molecular markers for rubber tree late SE. CML13, CML36, CAM-7, SERK1 and LEAD-29-like were also related to rubber tree late SE. In addition, histone modification had crucial roles during rubber tree late SE.

CONCLUSIONS

This study provides important information to elucidate the molecular regulation during rubber tree late SE.

Systematic investigation and expression profiles of the GbR2R3-MYB transcription factor family in ginkgo (Ginkgo biloba L.)

As one of the largest families of transcription factors, the R2R3-MYB family plays a significant role in plant growth, development, and response to hormone and environmental stress. To explore its evolutionary mechanism and potential function in Ginkgo biloba, a gymnosperm of great economic and ecological value, we presented a comprehensive analysis of the R2R3-MYB genes in ginkgo. Sixty-nine GbR2R3-MYB genes were identified and these genes could be classified into 33 groups based on the characteristics of the amino acid sequence of the R2R3-MYB domain and gene structure. Syntenic analyses indicated that few tandem and segmental duplications possibly resulted in the contraction of the GbR2R3-MYB gene family. Based on the transcriptome data, expression profiles of eight different tissues and different developmental stages of leaf and kernel showed that GbR2R3-MYB genes had distinct temporal and spatial expression characteristics. Specific expression patterns of the sixteen GbR2R3-MYB genes were also identified in response to different abiotic stresses and hormonal exposures. Further investigation revealed that GbR2R3-MYB19 was located in the nucleus and possessed transcriptional activity, implying its potential roles in the regulation of multiple biological processes. Our findings provide a robust basis for future comprehensive evolutionary and functional analyses of GbR2R3-MYB genes in ginkgo.

Identification and Expression Analysis of the CsMYB Gene Family in Root Knot Nematode-Resistant and Susceptible Cucumbers

MYB (myeloblastosis) transcription factors (TFs) play important roles in controlling various physiological processes in plants, such as responses to biotic and abiotic stress, metabolism, and defense. A previous study identified a gene, Csa6G410090, encoding a plant lipid transfer protein (LTP), as a possible regulator in cucumber (Cucumis sativus L.) of the resistance response to root-knot nematode (RKN) [Meloidogyne incognita Kofoid and White (Chitwood)]. Myb-type DNA-binding TFs were presumed to regulate downstream genes expression, including LTPs, however, the regulation mechanism remained unclear. To elucidate whether and which MYB TFs may be involved in regulation of the resistance response, this study identified 112 genes as candidate members of the CsMYB gene family by combining CDD and SMART databases, using the Hidden Markov Model (HMM) and manual calibration. Within this group, ten phylogenetic subgroups were resolved according to sequence-based classification, consistent with results from comprehensive investigation of gene structure, conserved motifs, chromosome locations, and cis-element analysis. Distribution and collinearity analysis indicated that amplification of the CsMYB gene family in cucumber has occurred mainly through tandem repeat events. Spatial gene expression analysis showed that 8 CsMYB genes were highly expressed at differing levels in ten different tissues or organs. The roots of RKN-resistant and susceptible cucumbers were inoculated with M. incognita, finding that CsMYB (Csa6G538700, Csa1G021940, and Csa5G641610) genes showed up-regulation coincident with upregulation of the "hub" gene LTP (Csa6G410090) previously implicated as a major gene in the resistance response to RKN in cucumber. Results of this study suggest hypotheses regarding the elements and regulation of the resistant response as well as possible RKN resistance-enhancing strategies in cucumber and perhaps more broadly in plants.

Genomic Comparison and Population Diversity Analysis Provide Insights into the Domestication and Improvement of Flax

Flax has been cultivated for its oil and fiber for thousands of years. However, it remains unclear how the modifications of agronomic traits occurred on the genetic level during flax cultivation. In this study, we conducted genome-wide variation analyses on multiple accessions of oil-use, fiber-use, landraces, and pale flax to identify the genomic variations during flax cultivation. Our findings indicate that, during flax domestication, genes relevant to flowering, dehiscence, oil production, and plant architecture were preferentially selected. Furthermore, regardless of origins, the improvement of the modern oil-use flax preceded that of the fiber-use flax, although the dual selection on oil-use and fiber-use characteristics might have occurred in the early flax domestication. We also found that the expansion of MYB46/MYB83 genes may have contributed to the unique secondary cell wall biosynthesis in flax and the directional selections on MYB46/MYB83 may have shaped the morphological profile of the current oil-use and fiber-use flax.

•

Assemblies of genomes, including oil-use flax, fiber-use flax and pale flax

•

Comparative genomic analysis between pale flax and cultivated flax

•

Dual-selection mode on oil-use and fiber-use characteristics might be existing

•

Expansion and selection of MYB46/MYB83 may shape the morphological profile of flax

Open-Bud Duplicate Loci Are Identified as MML10s, Orthologs of MIXTA-Like Genes on Homologous Chromosomes of Allotetraploid Cotton

The open-bud (ob) mutants in cotton display abnormal flower buds with the stigma and upper anthers exposed before blooming. This characteristic is potentially useful for the efficient production of hybrid seeds. The recessive inheritance pattern of the ob phenotype in allotetraploid cotton is determined by duplicated recessive loci (ob1ob1ob2ob2). In this study, ob1, which is a MIXTA-like MYB gene on chromosome D13 (MML10_Dt), was identified by map-based cloning. In Gossypium barbadense (Gb) acc. 3-79, a single nucleotide polymorphism (SNP) (G/A) at the splice site of the first intron and an 8-bp deletion in the third exon of MML10_Dt were found, which are the causative mutations at the ob1 loci. A 1783-bp deletion that leads to the loss of the third exon and accounts for the causal variation at the ob2 loci was found in MML10_At of Gossypium hirsutum (Gh) acc. TM-1. The ob phenotype results from the combination of these two loss-of-function loci. Genotyping assays showed that the ob1 and ob2 loci appeared after the formation of allotetraploid cotton and were specific for Gb and Gh, respectively. All Gb lines and most Gh cultivars carry the single corresponding mutant alleles. Genome-wide transcriptome analysis showed that some of the MYB genes and genes related to cell wall biogenesis, trichome differentiation, cytokinin signal transduction, and cell division were repressed in the ob mutants, which may lead to suppression of petal growth. These findings should be of value for breeding superior ob lines in cotton.

Transgenic expression of rice MYB102 (OsMYB102) delays leaf senescence and decreases abiotic stress tolerance in Arabidopsis thaliana

MYB-type transcription factors (TFs) play important roles in plant growth and development, and in the rapid responses to unfavorable environmental conditions. We recently reported the isolation and characterization of a rice (Oryza sativa) MYB TF, OsMYB102, which is involved in the regulation of leaf senescence by downregulating abscisic acid (ABA) biosynthesis and the downstream signaling response. Based on the similarities of their sequences and expression patterns, OsMYB102 appears to be a homolog of the Arabidopsis thaliana AtMYB44 TF. Since AtMYB44 is a key regulator of leaf senescence and abiotic stress responses, it is important to examine whether AtMYB44 homologs in other plants also act similarly. Here, we generated transgenic Arabidopsis plants expressing OsMYB102 (OsMYB102-OX). The OsMYB102-OX plants showed a delayed senescence phenotype during dark incubation and were more susceptible to salt and drought stresses, considerably similar to Arabidopsis plants overexpressing AtMYB44. Real-time quantitative PCR (RT-qPCR) revealed that, in addition to known senescence-associated genes, genes encoding the ABA catabolic enzymes AtCYP707A3 and AtCYP707A4 were also significantly upregulated in OsMYB102-OX, leading to a significant decrease in ABA accumulation. Furthermore, protoplast transient expression and chromatin immunoprecipitation assays revealed that OsMYB102 directly activated AtCYP707A3 expression. Based on our findings, it is probable that the regulatory functions of AtMYB44 homologs in plants are highly conserved and they have vital roles in leaf senescence and the abiotic stress responses.

Transcriptome profiling revealed diverse gene expression patterns in poplar (Populus × euramericana) under different planting densities

Certain plant genotypes can achieve optimal growth under appropriate environmental conditions. Under high planting density conditions, plants undergo competition for uptake and utilization of light and nutrients. However, the relationship between whole-genome expression patterns and the planting density in perennial woody plants remains unknown. In this study, whole-genome RNA sequencing of poplar (Populus × euramericana) was carried out at three different sampling heights to determine gene expression patterns under high (HD) and low (LD) planting densities. As a result, 4,004 differentially expressed genes (DEGs) were detected between HD and LD, of which 2,300, 701, and 1,003 were detected at the three positions, upper, middle and bottom, respectively. Function annotation results further revealed that a large number of the DEGs were involved in distinct biological functions. There were significant changes in the expression of metabolism-related and stimulus-related genes in response to planting density. There were 37 DEGs that were found at all three positions and were subsequently screened. Several DEGs related to plant light responses and photosynthesis were observed at different positions. Meanwhile, numbers of genes related to auxin/indole-3-acetic acid, and carbon and nitrogen metabolism were also revealed, displaying overall trends of upregulation under HD. These findings provide a basis for identifying candidate genes related to planting density and could increase our molecular understanding of the effect of planting density on gene expression.

Rice transcription factor OsMYB102 delays leaf senescence by down-regulating abscisic acid accumulation and signaling

MYB-type transcription factors (TFs) play important roles in plant growth and development, and in the responses to several abiotic stresses. In rice (Oryza sativa), the roles of MYB-related TFs in leaf senescence are not well documented. Here, we examined rice MYB TF gene OsMYB102 and found that an OsMYB102 T-DNA activation-tagged line (termed osmyb102-D), which constitutively expresses OsMYB102 under the control of four tandem repeats of the 35S promoter, and OsMYB102-overexpressing transgenic lines (35S:OsMYB102 and 35S:GFP-OsMYB102) maintain green leaves much longer than the wild-type under natural, dark-induced, and abscisic acid (ABA)-induced senescence conditions. Moreover, an osmyb102 knockout mutant showed an accelerated senescence phenotype under dark-induced and ABA-induced leaf senescence conditions. Microarray analysis showed that a variety of senescence-associated genes (SAGs) were down-regulated in the osmyb102-D line. Further studies demonstrated that overexpression of OsMYB102 controls the expression of SAGs, including genes associated with ABA degradation and ABA signaling (OsABF4, OsNAP, and OsCYP707A6), under dark-induced senescence conditions. OsMYB102 inhibits ABA accumulation by directly activating the transcription of OsCYP707A6, which encodes the ABA catabolic enzyme ABSCISIC ACID 8'-HYDROXYLASE. OsMYB102 also indirectly represses ABA-responsive genes, such as OsABF4 and OsNAP. Collectively, these results demonstrate that OsMYB102 plays a critical role in leaf senescence by down-regulating ABA accumulation and ABA signaling responses.

Genome-wide analysis and identification of the low potassium stress responsive gene SiMYB3 in foxtail millet (Setariaitalica L.)

BACKGROUND

Potassium (K) is essential to plant growth and development. Foxtail millet (Setaria italic L.) is an important fodder grain crop in arid and semi-arid regions of Asia and Africa because of its strong tolerance to drought and barren stresses. The molecular mechanisms of physiological and biochemical responses and regulations to various abiotic stresses such as low potassium conditions in foxtail millet are not fully understood, which hinders the research and exploitation of this valuable resource.

RESULTS

In this research, we demonstrated that the millet variety Longgu 25 was the most insensitive variety to low potassium stress among other five varieties. The transcriptome analysis of Longgu 25 variety revealed a total of 26,192 and 26,849 genes from the K⁺-deficient and normal transcriptomic libraries by RNA-seq, respectively. A total of 1982 differentially expressed genes (DEGs) were identified including 866 up-regulated genes and 1116 down-regulated genes. We conducted a comparative analysis of these DEGs under low-K⁺ stress conditions and discovered 248 common DEGs for potassium deprivation among foxtail millet, rice and Arabidopsis. Further Gene Ontology (GO) enrichment analysis identified a series of candidate genes that may involve in K⁺-deficient response and in intersection of molecular functions among foxtail millet, rice and Arabidopsis. The expression profiles of randomly selected 18 candidate genes were confirmed as true DEGs with RT-qPCR. Furthermore, one of the 18 DEGs, SiMYB3, is specifically expressed only in the millet under low-K⁺ stress conditions. Overexpression of SiMYB3 promoted the main root elongation and improved K⁺ deficiency tolerance in transgenic Arabidopsis plants. The fresh weight of the transgenic plants was higher, the primary root length was longer and the root surface-area was larger than those of control plants after K⁺ deficiency treatments.

CONCLUSIONS

This study provides a global view of transcriptomic resources relevant to the K⁺-deficient tolerance in foxtail millet, and shows that SiMYB3 is a valuable genetic resource for the improvement of K⁺ deficiency tolerance in foxtail millet.

Genome-wide characterization and expression analyses of the MYB superfamily genes during developmental stages in Chinese jujube

The MYB transcription factor (TF) superfamily, one of the largest gene superfamilies, regulates a variety of physiological processes in plants. Although many MYB superfamily genes have been identified in plants, the MYB TFs in Chinese jujube (Ziziphus jujuba Mill.) have not been fully identified and characterized. Additionally, the functions of these genes remain unclear. In total, we identified 171 MYB superfamily genes in jujube and divided them into five subfamilies containing 99 genes of the R2R3-MYB subfamily, 58 genes of the MYB-related subfamily, four genes of the R1R2R3-MYB subfamily, one gene of the 4R-MYB subfamily, and nine genes of the atypical MYB subfamily. The 99 R2R3-MYB genes of jujube were divided into 35 groups, C1–C35, and the 58 MYB-related genes were divided into the following groups: the R-R-type, CCA1-like, I-box-binding-like, TBP-like, CPC-like, and Chinese jujube-specific groups. ZjMYB genes in jujube were well supported by additional highly conserved motifs and exon/intron structures. Most R1 repeats of MYB-related proteins comprised the R2 repeat and had highly conserved EED and EEE residue groups in jujube. Three tandem duplicated gene pairs were found on 12 chromosomes in jujube. According to an expression analysis of 126 ZjMYB genes, MYB-related genes played important roles in jujube development and fruit-related biological processes. The total flavonoid content of jujube fruit decreased as ripening progressed. A total of 93 expressed genes were identified in the RNA-sequencing data from jujube fruit, and 56 ZjMYB members presented significant correlations with total flavonoid contents by correlation analysis. Five pairs of paralogous MYB genes within jujube were composed of nine jujube MYB genes. A total of 14 ZjMYB genes had the same homology to the MYB genes of Arabidopsis and peach, indicating that these 14 MYB genes and their orthologs probably existed before the ancestral divergence of the MYB superfamily. We used a synteny analysis of MYB genes between jujube and Arabidopsis to predict that the functions of the ZjMYBs involve flavonoid/phenylpropanoid metabolism, the light signaling pathway, auxin signal transduction, and responses to various abiotic stresses (cold, drought, and salt stresses). Additionally, we speculate that ZjMYB108 is an important TF involved in the flavonoid metabolic pathway. This comprehensive analysis of MYB superfamily genes in jujube lay a solid foundation for future comprehensive analyses of ZjMYB gene functions.

Structural, Functional, and Evolutionary Characterization of Major Drought Transcription Factors Families in Maize

Drought is one of the major threats to the maize yield especially in subtropical production systems. Understanding the genes and regulatory mechanisms of drought tolerance is important to sustain the yield. Transcription factors (TFs) play a major role in gene regulation under drought stress. In the present study, a set of 15 major TF families comprising 1,436 genes was structurally and functionally characterized. The functional annotation indicated that the genes were involved in ABA signaling, ROS scavenging, photosynthesis, stomatal regulation, and sucrose metabolism. Duplication was identified as the primary force in divergence and expansion of TF families. Phylogenetic relationship was developed for individual TF and combined TF families. Phylogenetic analysis clustered the genes into specific and mixed groups. Gene structure analysis revealed that more number of genes were intron-rich as compared to intron-less. Drought-responsive cis-regulatory elements such as ABREA, ABREB, DRE1, and DRECRTCOREAT have been identified. Expression and interaction analyses identified leaf-specific bZIP TF, GRMZM2G140355, as a potential contributor toward drought tolerance in maize. Protein-protein interaction network of 269 drought-responsive genes belonging to different TFs has been provided. The information generated on structural and functional characteristics, expression, and interaction of the drought-related TF families will be useful to decipher the drought tolerance mechanisms and to breed drought-tolerant genotypes in maize.

Rice OsMYB5P improves plant phosphate acquisition by regulation of phosphate transporter

Myeloblastosis (MYB) transcription factors play central roles in plant developmental processes and in responses to nutrient deficiency. In this study, OsMYB5P, an R2R3-MYB transcription factor, was isolated and identified from rice (Oryza sativa L. 'Dongjin') under inorganic phosphate (Pi)-deficient conditions. OsMYB5P protein is localized to the nucleus and functions as a transcription activator in plant development. Overexpression of OsMYB5P in rice and Arabidopsis (Arabidopsis thaliana Col-0) increases tolerance to phosphate starvation, whereas OsMYB5P knock-out through RNA interference increases sensitivity to Pi depletion in rice. Furthermore, shoots and roots of transgenic rice plants overexpressing OsMYB5P were longer than those of wild plants under both normal and Pi-deficient conditions. These results indicate that OsMYB5P is associated with the regulation of shoot development and root- system architecture. Overexpression of OsMYB5P led to increased Pi accumulation in shoots and roots. Interestingly, OsMYB5P directly bound to MBS (MYB binding site) motifs on the OsPT5 promoter and induced transcription of OsPT5 in rice. In addition, overexpression of OsMYB5P in Arabidopsis triggered increased expression of AtPht1;3, an Arabidopsis Pi transporter, in shoots and roots under normal and Pi-deficient conditions. Together, these results demonstrate that overexpression of OsMYB5P increases tolerance to Pi deficiency in plants by modulating Pi transporters at the transcriptional level in monocots and dicots.

Transcriptome-Wide Identification and Characterization of MYB Transcription Factor Genes in the Laticifer Cells of Hevea brasiliensis

MYB transcription factors hold vital roles in the regulation of plant secondary metabolic pathways. Laticifers in rubber trees (Hevea brasiliensis) are of primary importance in natural rubber production because natural rubber is formed and stored within these structures. To understand the role of MYB transcription factors in the specialized cells, we identified 44 MYB genes (named HblMYB1 to HblMYB44) by using our previously obtained transcriptome database of rubber tree laticifer cells and the public rubber tree genome database. Expression profiles showed that five MYB genes were highly expressed in the laticifers. HblMYB19 and HblMYB44 were selected for further study. HblMYB19 and HblMYB44 bound the promoters of HbFDPS1, HbSRPP, and HRT1 in yeast. Furthermore, the transient overexpression of HblMYB19 and HblMYB44 in tobacco plants significantly increased the activity of the promoters of HbFDPS1, HbSRPP, and HRT1. Basing on this information, we proposed that HblMYB19 and HblMYB44 are the regulators of HbFDPS1, HbSRPP, and HRT1, which are involved in the biosynthesis pathway of natural rubber.

Evidence for a Strong Correlation Between Transcription Factor Protein Disorder and Organismic Complexity

Studies of diverse phylogenetic lineages reveal that protein disorder increases in concert with organismic complexity but that differences nevertheless exist among lineages. To gain insight into this phenomenology, we analyzed all of the transcription factor (TF) families for which sequences are known for 17 species spanning bacteria, yeast, algae, land plants, and animals and for which the number of different cell types has been reported in the primary literature. Although the fraction of disordered residues in TF sequences is often moderately or poorly correlated with organismic complexity as gauged by cell-type number (r² < 0.5), an unbiased and phylogenetically broad analysis shows that organismic complexity is positively and strongly correlated with the total number of TFs, the number of their spliced variants and their total disordered residues content (r² > 0.8). Furthermore, the correlation between the fraction of disordered residues and cell-type number becomes stronger when confined to the TF families participating in cell cycle, cell size, cell division, cell differentiation, or cell proliferation, and other important developmental processes. The data also indicate that evolutionarily simpler organisms allow for the detection of subtle differences in the conserved IDRs of TFs as well as changes in variable IDRs, which can influence the DNA recognition and multifunctionality of TFs through direct or indirect mechanisms. Although strong correlations cannot be taken as evidence for cause-and-effect relationships, we interpret our data to indicate that increasing TF disorder likely was an important factor contributing to the evolution of organismic complexity and not merely a concurrent unrelated effect of increasing organismic complexity.

Evolution of RAD- and DIV-Like Genes in Plants

Developmental genetic studies of Antirrhinum majus demonstrated that two transcription factors from the MYB gene family, RADIALIS (RAD) and DIVIRICATA (DIV), interact through antagonism to regulate floral dorsoventral asymmetry. Interestingly, similar antagonistic interaction found among proteins of FSM1 (RAD-like) and MYBI (DIV-like) in Solanum lycopersicum is involved in fruit development. Here, we report the reconstruction of the phylogeny of I-box-like and R-R-type clades, where RAD- and DIV-like genes belong, respectively. We also examined the homology of these antagonistic MYB proteins using these phylogenies. The results show that there are likely three paralogs of RAD-/I-box-like genes, RAD1, RAD2, and RAD3, which originated in the common ancestor of the core eudicots. In contrast, R-R-type sequences fall into two major clades, RR1 and RR2, the result of gene duplication in the common ancestor of both monocots and dicots. RR1 was divided into clades RR1A, RR1B, and RR1C, while RR2 was divided into clades RR2A/DIV1, RR2B/DIV2, and RR2C/DIV3. We demonstrate that among similar antagonistic interactions in An. Majus and So. lycopersicum, RAD-like genes originate from the RAD2 clade, while DIV-like genes originate from distantly related paralogs of the R-R-type lineage. The phylogenetic analyses of these two MYB clades lay the foundation for future comparative studies including testing the evolution of the antagonistic relationship of proteins.

Transcription dynamics of Saltol QTL localized genes encoding transcription factors, reveals their differential regulation in contrasting genotypes of rice

Salinity is one of the major environmental factors affecting the growth and yield of rice crop. Salinity stress response is a multigenic trait and numerous approaches have been used to dissect out the key determinants of complex salt tolerance trait and their regulation in plant. In the current study, we have investigated expression dynamics of the genes encoding transcription factors (SalTFs) localized within a major salinity tolerance related QTL-'Saltol' in the contrasting cultivars of rice. SalTFs were found to be differentially regulated between the contrasting genotypes of rice, with higher constitutive expression in the salt tolerant landrace, Pokkali than the cultivar IR64. Moreover, SalTFs were found to exhibit inducibility in the salt sensitive cultivar at late duration (after 24 h) of salinity stress. Further, the transcript abundance analysis of these SalTFs at various developmental stages of rice revealed that low expressing genes may be involved in developmental responses, while high expressing genes can be linked with the salt stress response. Grouping of these genes was well supported by in silico protein-protein interaction studies and distribution of single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) in the promoter and genic regions of these genes. Taken together, we propose that out of 14 SalTFs, eight members are strongly correlated with the salinity stress tolerance in rice and six are involved in plant growth and development.

AtDOF5.4/OBP4, a DOF Transcription Factor Gene that Negatively Regulates Cell Cycle Progression and Cell Expansion in Arabidopsis thaliana

In contrast to animals, plant development involves continuous organ formation, which requires strict regulation of cell proliferation. The core cell cycle machinery is conserved across plants and animals, but plants have developed new mechanisms that precisely regulate cell proliferation in response to internal and external stimuli. Here, we report that the DOF transcription factor OBP4 negatively regulates cell proliferation and expansion. OBP4 is a nuclear protein. Constitutive and inducible overexpression of OBP4 reduced the cell size and number, resulting in dwarf plants. Inducible overexpression of OBP4 in Arabidopsis also promoted early endocycle onset and inhibited cell expansion, while inducible overexpression of OBP4 fused to the VP16 activation domain in Arabidopsis delayed endocycle onset and promoted plant growth. Furthermore, gene expression analysis showed that cell cycle regulators and cell wall expansion factors were largely down-regulated in the OBP4 overexpression lines. Short-term inducible analysis coupled with in vivo ChIP assays indicated that OBP4 targets the CyclinB1;1, CDKB1;1 and XTH genes. These results strongly suggest that OBP4 is a negative regulator of cell cycle progression and cell growth. These findings increase our understanding of the transcriptional regulation of the cell cycle in plants.

An iNTT system for the large-scale screening of differentially expressed, nuclear-targeted proteins: cold-treatment-induced nucleoproteins in Rye (Secale cereale L.)

Background

Nuclear proteins play critical roles in regulating mRNA transcription and processing, DNA replication, and epigenetic genome modification. Therefore, the ability to monitor changes in nuclear proteins is helpful not only to identify important regulatory proteins but also to study the mechanisms of actions of nuclear proteins. However, no effective methods have been developed yet. Rye is strongly resistant to various biotic and abiotic stresses; however, few genes have been functionally characterized to date due to the complexity of its genome and a lack of genomic sequence information.

Results

We developed an integrative Nuclear Transportation Trap (iNTT) system that includes an improved nuclear transportation trap and utilizes the “after suppression subtraction” method. Oligonucleotides encoding a nuclear localization signal (NLS) or a transcription factor, GmAREB, were inserted into pLexAD or pLexAD-NES, respectively, and then transformed into yeast cells (EGY48). We showed that the pLexAD vector expressing a cDNA library in the iNTT system was more efficient for screening than the vector pLexAD-NES, which has previously been used in an NTT system. We used the iNTT system to screen a cDNA library of cold-treated rye. A total of 241 unique genes were identified, including 169 differentially expressed proteins; of these, 106 were of known and 63 were of unknown function. Moreover, 82 genes (49 %) among the 169 differentially expressed genes were predicted to contain an NLS domain. Thirty-three (31 %) of the 106 functionally known proteins have DNA-binding activity. To test the specificity of the nuclear proteins identified using the iNTT screen, four of the proteins differentially expressed in response to temperature stress, ScT1 (a heat shock protein), ScT36 (a MYB-like transcription factor), ScT133 (an ERF-like transcription factor) and ScT196 (a protein of unknown function), were studied in more depth. These proteins were shown to exclusively localize to the nucleus, and their expression levels were increased in response to low-temperature stress. To identify the function of these screened nuclear proteins, ScT1- and ScT36-transgenic Arabidopsis plants were constructed, and ScT1 or ScT36 overexpression was found to enhance tolerance to high-temperature or freezing stresses, respectively.

Conclusions

The newly developed iNTT system provides an effective method for identifying nuclear-targeted proteins and monitoring induced expression levels. ScT1 and ScT36 might be good candidate genes for improving the stress tolerance of plants by genetic transformation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2548-y) contains supplementary material, which is available to authorized users.

Comparative cell-specific transcriptomics reveals differentiation of C4 photosynthesis pathways in switchgrass and other C4 lineages

Almost all C4 plants require the co-ordination of the adjacent and fully differentiated cell types, mesophyll (M) and bundle sheath (BS). The C4 photosynthetic pathway operates through two distinct subtypes based on how malate is decarboxylated in BS cells; through NAD-malic enzyme (NAD-ME) or NADP-malic enzyme (NADP-ME). The diverse or unique cell-specific molecular features of M and BS cells from separate C4 subtypes of independent lineages remain to be determined. We here provide an M/BS cell type-specific transcriptome data set from the monocot NAD-ME subtype switchgrass (Panicum virgatum). A comparative transcriptomics approach was then applied to compare the M/BS mRNA profiles of switchgrass, monocot NADP-ME subtype C4 plants maize and Setaria viridis, and dicot NAD-ME subtype Cleome gynandra. We evaluated the convergence in the transcript abundance of core components in C4 photosynthesis and transcription factors to establish Kranz anatomy, as well as gene distribution of biological functions, in these four independent C4 lineages. We also estimated the divergence between NAD-ME and NADP-ME subtypes of C4 photosynthesis in the two cell types within C4 species, including differences in genes encoding decarboxylating enzymes, aminotransferases, and metabolite transporters, and differences in the cell-specific functional enrichment of RNA regulation and protein biogenesis/homeostasis. We suggest that C4 plants of independent lineages in both monocots and dicots underwent convergent evolution to establish C4 photosynthesis, while distinct C4 subtypes also underwent divergent processes for the optimization of M and BS cell co-ordination. The comprehensive data sets in our study provide a basis for further research on evolution of C4 species.

RNA-seq for gene identification and transcript profiling in relation to root growth of bermudagrass (Cynodon dactylon) under salinity stress

BACKGROUND

Soil salinity is one of the most significant abiotic stresses affecting plant shoots and roots growth. The adjustment of root architecture to spatio-temporal heterogeneity in salinity is particularly critical for plant growth and survival. Bermudagrass (Cynodon dactylon) is a widely used turf and forage perennial grass with a high degree of salinity tolerance. Salinity appears to stimulate the growth of roots and decrease their mortality in tolerant bermudagrass. To estimate a broad spectrum of genes related to root elongation affected by salt stress and the molecular mechanisms that control the positive response of root architecture to salinity, we analyzed the transcriptome of bermudagrass root tips in response to salinity.

RESULTS

RNA-sequencing was performed in root tips of two bermudagrass genotypes contrasting in salt tolerance. A total of 237,850,130 high quality clean reads were generated and 250,359 transcripts were assembled with an average length of 1115 bp. Totally, 103,324 unigenes obtained with 53,765 unigenes (52 %) successfully annotated in databases. Bioinformatics analysis indicated that major transcription factor (TF) families linked to stress responses and growth regulation (MYB, bHLH, WRKY) were differentially expressed in root tips of bermudagrass under salinity. In addition, genes related to cell wall loosening and stiffening (xyloglucan endotransglucosylase/hydrolases, peroxidases) were identified.

CONCLUSIONS

RNA-seq analysis identified candidate genes encoding TFs involved in the regulation of lignin synthesis, reactive oxygen species (ROS) homeostasis controlled by peroxidases, and the regulation of phytohormone signaling that promote cell wall loosening and therefore root growth under salinity.

Functional characterization of drought-responsive modules and genes in Oryza sativa: a network-based approach

Drought is one of the major environmental stress conditions affecting the yield of rice across the globe. Unraveling the functional roles of the drought-responsive genes and their underlying molecular mechanisms will provide important leads to improve the yield of rice. Co-expression relationships derived from condition-dependent gene expression data is an effective way to identify the functional associations between genes that are part of the same biological process and may be under similar transcriptional control. For this purpose, vast amount of freely available transcriptomic data may be used. In this study, we consider gene expression data for different tissues and developmental stages in response to drought stress. We analyze the network of co-expressed genes to identify drought-responsive genes modules in a tissue and stage-specific manner based on differential expression and gene enrichment analysis. Taking cues from the systems-level behavior of these modules, we propose two approaches to identify clusters of tightly co-expressed/co-regulated genes. Using graph-centrality measures and differential gene expression, we identify biologically informative genes that lack any functional annotation. We show that using orthologous information from other plant species, the conserved co-expression patterns of the uncharacterized genes can be identified. Presence of a conserved neighborhood enables us to extrapolate functional annotation. Alternatively, we show that single 'guide-gene' approach can help in understanding tissue-specific transcriptional regulation of uncharacterized genes. Finally, we confirm the predicted roles of uncharacterized genes by the analysis of conserved cis-elements and explain the possible roles of these genes toward drought tolerance.

Characterization of the cork oak transcriptome dynamics during acorn development

BACKGROUND

Cork oak (Quercus suber L.) has a natural distribution across western Mediterranean regions and is a keystone forest tree species in these ecosystems. The fruiting phase is especially critical for its regeneration but the molecular mechanisms underlying the biochemical and physiological changes during cork oak acorn development are poorly understood. In this study, the transcriptome of the cork oak acorn, including the seed, was characterized in five stages of development, from early development to acorn maturation, to identify the dominant processes in each stage and reveal transcripts with important functions in gene expression regulation and response to water.

RESULTS

A total of 80,357 expressed sequence tags (ESTs) were de novo assembled from RNA-Seq libraries representative of the several acorn developmental stages. Approximately 7.6 % of the total number of transcripts present in Q. suber transcriptome was identified as acorn specific. The analysis of expression profiles during development returned 2,285 differentially expressed (DE) transcripts, which were clustered into six groups. The stage of development corresponding to the mature acorn exhibited an expression profile markedly different from other stages. Approximately 22 % of the DE transcripts putatively code for transcription factors (TF) or transcriptional regulators, and were found almost equally distributed among the several expression profile clusters, highlighting their major roles in controlling the whole developmental process. On the other hand, carbohydrate metabolism, the biological pathway most represented during acorn development, was especially prevalent in mid to late stages as evidenced by enrichment analysis. We further show that genes related to response to water, water deprivation and transport were mostly represented during the early (S2) and the last stage (S8) of acorn development, when tolerance to water desiccation is possibly critical for acorn viability.

CONCLUSIONS

To our knowledge this work represents the first report of acorn development transcriptomics in oaks. The obtained results provide novel insights into the developmental biology of cork oak acorns, highlighting transcripts putatively involved in the regulation of the gene expression program and in specific processes likely essential for adaptation. It is expected that this knowledge can be transferred to other oak species of great ecological value.

Transcriptional Regulatory Network Analysis of MYB Transcription Factor Family Genes in Rice

MYB transcription factor (TF) is one of the largest TF families and regulates defense responses to various stresses, hormone signaling as well as many metabolic and developmental processes in plants. Understanding these regulatory hierarchies of gene expression networks in response to developmental and environmental cues is a major challenge due to the complex interactions between the genetic elements. Correlation analyses are useful to unravel co-regulated gene pairs governing biological process as well as identification of new candidate hub genes in response to these complex processes. High throughput expression profiling data are highly useful for construction of co-expression networks. In the present study, we utilized transcriptome data for comprehensive regulatory network studies of MYB TFs by "top-down" and "guide-gene" approaches. More than 50% of OsMYBs were strongly correlated under 50 experimental conditions with 51 hub genes via "top-down" approach. Further, clusters were identified using Markov Clustering (MCL). To maximize the clustering performance, parameter evaluation of the MCL inflation score (I) was performed in terms of enriched GO categories by measuring F-score. Comparison of co-expressed cluster and clads analyzed from phylogenetic analysis signifies their evolutionarily conserved co-regulatory role. We utilized compendium of known interaction and biological role with Gene Ontology enrichment analysis to hypothesize function of coexpressed OsMYBs. In the other part, the transcriptional regulatory network analysis by "guide-gene" approach revealed 40 putative targets of 26 OsMYB TF hubs with high correlation value utilizing 815 microarray data. The putative targets with MYB-binding cis-elements enrichment in their promoter region, functional co-occurrence as well as nuclear localization supports our finding. Specially, enrichment of MYB binding regions involved in drought-inducibility implying their regulatory role in drought response in rice. Thus, the co-regulatory network analysis facilitated the identification of complex OsMYB regulatory networks, and candidate target regulon genes of selected guide MYB genes. The results contribute to the candidate gene screening, and experimentally testable hypotheses for potential regulatory MYB TFs, and their targets under stress conditions.

Early transcriptional response of soybean contrasting accessions to root dehydration

Drought is a significant constraint to yield increase in soybean. The early perception of water deprivation is critical for recruitment of genes that promote plant tolerance. DeepSuperSAGE libraries, including one control and a bulk of six stress times imposed (from 25 to 150 min of root dehydration) for drought-tolerant and sensitive soybean accessions, allowed to identify new molecular targets for drought tolerance. The survey uncovered 120,770 unique transcripts expressed by the contrasting accessions. Of these, 57,610 aligned with known cDNA sequences, allowing the annotation of 32,373 unitags. A total of 1,127 unitags were up-regulated only in the tolerant accession, whereas 1,557 were up-regulated in both as compared to their controls. An expression profile concerning the most representative Gene Ontology (GO) categories for the tolerant accession revealed the expression "protein binding" as the most represented for "Molecular Function", whereas CDPK and CBL were the most up-regulated protein families in this category. Furthermore, particular genes expressed different isoforms according to the accession, showing the potential to operate in the distinction of physiological behaviors. Besides, heat maps comprising GO categories related to abiotic stress response and the unitags regulation observed in the expression contrasts covering tolerant and sensitive accessions, revealed the unitags potential for plant breeding. Candidate genes related to "hormone response" (LOX, ERF1b, XET), "water response" (PUB, BMY), "salt stress response" (WRKY, MYB) and "oxidative stress response" (PER) figured among the most promising molecular targets. Additionally, nine transcripts (HMGR, XET, WRKY20, RAP2-4, EREBP, NAC3, PER, GPX5 and BMY) validated by RT-qPCR (four different time points) confirmed their differential expression and pointed that already after 25 minutes a transcriptional reorganization started in response to the new condition, with important differences between both accessions.

High-throughput sequencing identification of novel and conserved miRNAs in the Brassica oleracea leaves

Background

Plant microRNAs are short (~21 nt) non-coding molecules that regulate gene expression by targeting the mRNA cleavage or protein translation inhibition. In this manner, they play many important roles in the cells of living organisms. One of the plant species in which the entire set of miRNAs has not been yet completely identified is Brassica oleracea var. capitata (cabbage). For this reason and for the economic and nutritional importance of this food crop, high-throughput small RNAs sequencing has been performed to discover the novel and conserved miRNAs in mature cabbage leaves.

Results

In this study, raw reads generated from three small RNA libraries were bioinformatically processed and further analyzed to select sequences homologous to known B. oleracea and other plant miRNAs. As a result of this analysis, 261 conserved miRNAs (belonging to 62 families) have been discovered. MIR169, MIR167 and MIR166 were the largest miRNA families, while the highest abundance molecules were miR167, miR166, miR168c and miR157a. Among the generated sequencing reads, miRNAs* were also found, such as the miR162c*, miR160a* and miR157a*. The unannotated tags were used in the prediction and evaluation of novel miRNAs, which resulted in the 26 potential miRNAs proposal. The expressions of 13 selected miRNAs were analyzed by northern blot hybridization. The target prediction and annotation for identified miRNAs were performed, according to which discovered molecules may target mRNAs encoding several potential proteins – e.g., transcription factors, polypeptides that regulate hormone stimuli and abiotic stress response, and molecules participating in transport and cell communication. Additionally, KEGG maps analysis suggested that the miRNAs in cabbage are involved in important processing pathways, including glycolysis, glycerolipid metabolism, flavonoid biosynthesis and oxidative phosphorylation.

Conclusions

Conclusively, for the first time, the large set of miRNAs was identified in mature cabbage leaves. Potential targets designation for these miRNAs may suggest their essential role in many plants primary biological processes. Presented study not only supplements the knowledge about B. oleracea miRNAs, but additionally it may be used in other research concerning the improvement of the cabbage cultivation.

Genome-wide classification and expression analysis of MYB transcription factor families in rice and Arabidopsis

BACKGROUND

The MYB gene family comprises one of the richest groups of transcription factors in plants. Plant MYB proteins are characterized by a highly conserved MYB DNA-binding domain. MYB proteins are classified into four major groups namely, 1R-MYB, 2R-MYB, 3R-MYB and 4R-MYB based on the number and position of MYB repeats. MYB transcription factors are involved in plant development, secondary metabolism, hormone signal transduction, disease resistance and abiotic stress tolerance. A comparative analysis of MYB family genes in rice and Arabidopsis will help reveal the evolution and function of MYB genes in plants.

RESULTS

A genome-wide analysis identified at least 155 and 197 MYB genes in rice and Arabidopsis, respectively. Gene structure analysis revealed that MYB family genes possess relatively more number of introns in the middle as compared with C- and N-terminal regions of the predicted genes. Intronless MYB-genes are highly conserved both in rice and Arabidopsis. MYB genes encoding R2R3 repeat MYB proteins retained conserved gene structure with three exons and two introns, whereas genes encoding R1R2R3 repeat containing proteins consist of six exons and five introns. The splicing pattern is similar among R1R2R3 MYB genes in Arabidopsis. In contrast, variation in splicing pattern was observed among R1R2R3 MYB members of rice. Consensus motif analysis of 1kb upstream region (5' to translation initiation codon) of MYB gene ORFs led to the identification of conserved and over-represented cis-motifs in both rice and Arabidopsis. Real-time quantitative RT-PCR analysis showed that several members of MYBs are up-regulated by various abiotic stresses both in rice and Arabidopsis.

CONCLUSION

A comprehensive genome-wide analysis of chromosomal distribution, tandem repeats and phylogenetic relationship of MYB family genes in rice and Arabidopsis suggested their evolution via duplication. Genome-wide comparative analysis of MYB genes and their expression analysis identified several MYBs with potential role in development and stress response of plants.

Genome-wide classification and expression analysis of MYB transcription factor families in rice and Arabidopsis

BACKGROUND

The MYB gene family comprises one of the richest groups of transcription factors in plants. Plant MYB proteins are characterized by a highly conserved MYB DNA-binding domain. MYB proteins are classified into four major groups namely, 1R-MYB, 2R-MYB, 3R-MYB and 4R-MYB based on the number and position of MYB repeats. MYB transcription factors are involved in plant development, secondary metabolism, hormone signal transduction, disease resistance and abiotic stress tolerance. A comparative analysis of MYB family genes in rice and Arabidopsis will help reveal the evolution and function of MYB genes in plants.

RESULTS

A genome-wide analysis identified at least 155 and 197 MYB genes in rice and Arabidopsis, respectively. Gene structure analysis revealed that MYB family genes possess relatively more number of introns in the middle as compared with C- and N-terminal regions of the predicted genes. Intronless MYB-genes are highly conserved both in rice and Arabidopsis. MYB genes encoding R2R3 repeat MYB proteins retained conserved gene structure with three exons and two introns, whereas genes encoding R1R2R3 repeat containing proteins consist of six exons and five introns. The splicing pattern is similar among R1R2R3 MYB genes in Arabidopsis. In contrast, variation in splicing pattern was observed among R1R2R3 MYB members of rice. Consensus motif analysis of 1kb upstream region (5' to translation initiation codon) of MYB gene ORFs led to the identification of conserved and over-represented cis-motifs in both rice and Arabidopsis. Real-time quantitative RT-PCR analysis showed that several members of MYBs are up-regulated by various abiotic stresses both in rice and Arabidopsis.

CONCLUSION

A comprehensive genome-wide analysis of chromosomal distribution, tandem repeats and phylogenetic relationship of MYB family genes in rice and Arabidopsis suggested their evolution via duplication. Genome-wide comparative analysis of MYB genes and their expression analysis identified several MYBs with potential role in development and stress response of plants.

Features of a unique intronless cluster of class I small heat shock protein genes in tandem with box C/D snoRNA genes on chromosome 6 in tomato (Solanum lycopersicum)

Physical clustering of genes has been shown in plants; however, little is known about gene clusters that have different functions, particularly those expressed in the tomato fruit. A class I 17.6 small heat shock protein (Sl17.6 shsp) gene was cloned and used as a probe to screen a tomato (Solanum lycopersicum) genomic library. An 8.3-kb genomic fragment was isolated and its DNA sequence determined. Analysis of the genomic fragment identified intronless open reading frames of three class I shsp genes (Sl17.6, Sl20.0, and Sl20.1), the Sl17.6 gene flanked by Sl20.1 and Sl20.0, with complete 5' and 3' UTRs. Upstream of the Sl20.0 shsp, and within the shsp gene cluster, resides a box C/D snoRNA cluster made of SlsnoR12.1 and SlU24a. Characteristic C and D, and C' and D', boxes are conserved in SlsnoR12.1 and SlU24a while the upstream flanking region of SlsnoR12.1 carries TATA box 1, homol-E and homol-D box-like cis sequences, TM6 promoter, and an uncharacterized tomato EST. Molecular phylogenetic analysis revealed that this particular arrangement of shsps is conserved in tomato genome but is distinct from other species. The intronless genomic sequence is decorated with cis elements previously shown to be responsive to cues from plant hormones, dehydration, cold, heat, and MYC/MYB and WRKY71 transcription factors. Chromosomal mapping localized the tomato genomic sequence on the short arm of chromosome 6 in the introgression line (IL) 6-3. Quantitative polymerase chain reaction analysis of gene cluster members revealed differential expression during ripening of tomato fruit, and relatively different abundances in other plant parts.

R2R3-NaMYB8 regulates the accumulation of phenylpropanoid-polyamine conjugates, which are essential for local and systemic defense against insect herbivores in Nicotiana attenuata

Although phenylpropanoid-polyamine conjugates (PPCs) occur ubiquitously in plants, their biological roles remain largely unexplored. The two major PPCs of Nicotiana attenuata plants, caffeoylputrescine (CP) and dicaffeoylspermidine, increase dramatically in local and systemic tissues after herbivore attack and simulations thereof. We identified NaMYB8, a homolog of NtMYBJS1, which in BY-2 cells regulates PPC biosynthesis, and silenced its expression by RNA interference in N. attenuata (ir-MYB8), to understand the ecological role(s) of PPCs. The regulatory role of NaMYB8 in PPC biosynthesis was validated by a microarray analysis, which revealed that transcripts of several key biosynthetic genes in shikimate and polyamine metabolism accumulated in a NaMYB8-dependent manner. Wild-type N. attenuata plants typically contain high levels of PPCs in their reproductive tissues; however, NaMYB8-silenced plants that completely lacked CP and dicaffeoylspermidine showed no changes in reproductive parameters of the plants. In contrast, a defensive role for PPCs was clear; both specialist (Manduca sexta) and generalist (Spodoptera littoralis) caterpillars feeding on systemically preinduced young stem leaves performed significantly better on ir-MYB8 plants lacking PPCs compared with wild-type plants expressing high levels of PPCs. Moreover, the growth of M. sexta caterpillars was significantly reduced when neonates were fed ir-MYB8 leaves sprayed with synthetic CP, corroborating the role of PPCs as direct plant defense. The spatiotemporal accumulation and function of PPCs in N. attenuata are consistent with the predictions of the optimal defense theory: plants preferentially protect their most fitness-enhancing and vulnerable parts, young tissues and reproductive organs, to maximize their fitness.