Plant spliceosomal introns: not only cut and paste

Characterization of the Metallothionein Gene Family in Avena sativa L. and the Gene Expression during Seed Germination and Heavy Metal Stress

Metallothioneins (MTs) are a family of small proteins rich in cysteine residues. The sulfhydryl group of metallothioneins can bind to metal ions, maintaining metal homeostasis and protecting the cells from damage caused by toxic heavy metals. Moreover, MTs can function as reactive oxygen species scavengers since cysteine thiols undergo reversible and irreversible oxidation. Here, we identified 21 metallothionein genes (AsMTs) in the oat (Avena sativa L.) genome, which were divided into four types depending on the amino acid sequences of putative proteins encoded by identified genes. Analysis of promoter sequences showed that MTs might respond to a variety of stimuli, including biotic and abiotic stresses and phytohormones. The results of qRT-PCR showed that all four types of AsMTs are differentially expressed during the first 48 hours of seed germination. Moreover, stress induced by the application of zinc, cadmium, and a mixture of zinc and cadmium affects the expression of oat MTs variously depending on the MT type, indicating that AsMT1-4 fulfil different roles in plant cells.

Melatonin receptor, GhCAND2-D5 motivated responding to NaCl signaling in cotton

As a receptor for plant melatonin, CAND2/PMTR plays an important role in melatonin signaling. Most of the CANDs are membrane proteins and play indispensable roles in signal transduction. In this study, the CANDs from four cotton species were characterized, and the phylogenetic relationships, expression patterns, stress responses of cotton CANDs were analyzed by bioinformatics. Through the analysis of phylogenetic and protein structure, it was found that the CANDs in clade Ⅱ might function as cotton melatonin receptors, and most of the GhCANDs in clade Ⅱ were induced by melatonin. A putative cotton melatonin receptor, GhCAND2-D5, was functionally probed by gene silencing. The plants with silenced expression of this gene exhibited decreased salt tolerance. Protein interaction prediction identified that GhCAND2-D5 interacted with several membrane proteins and played an important role in melatonin signaling. This study provided a theoretical reference for further investigation of melatonin signaling in cotton.

In Silico RNAseq and Biochemical Analyses of Glucose-6-Phosphate Dehydrogenase (G6PDH) from Sweet Pepper Fruits: Involvement of Nitric Oxide (NO) in Ripening and Modulation

Pepper (Capsicum annuum L.) fruit is a horticultural product consumed worldwide which has great nutritional and economic relevance. Besides the phenotypical changes that pepper fruit undergo during ripening, there are many associated modifications at transcriptomic, proteomic, biochemical, and metabolic levels. Nitric oxide (NO) is a recognized signal molecule that can exert regulatory functions in diverse plant processes including fruit ripening, but the relevance of NADPH as a fingerprinting of the crop physiology including ripening has also been proposed. Glucose-6-phosphate dehydrogenase (G6PDH) is the first and rate-limiting enzyme of the oxidative phase of the pentose phosphate pathway (oxiPPP) with the capacity to generate NADPH. Thus far, the available information on G6PDH and other NADPH-generating enzymatic systems in pepper plants, and their expression during the ripening of sweet pepper fruit, is very scarce. Therefore, an analysis at the transcriptomic, molecular and functional levels of the G6PDH system has been accomplished in this work for the first time. Based on a data-mining approach to the pepper genome and fruit transcriptome (RNA-seq), four G6PDH genes were identified in pepper plants and designated CaG6PDH1 to CaG6PDH4, with all of them also being expressed in fruits. While CaG6PDH1 encodes a cytosolic isozyme, the other genes code for plastid isozymes. The time-course expression analysis of these CaG6PDH genes during different fruit ripening stages, including green immature (G), breaking point (BP), and red ripe (R), showed that they were differentially modulated. Thus, while CaG6PDH2 and CaG6PDH4 were upregulated at ripening, CaG6PDH1 was downregulated, and CaG6PDH3 was slightly affected. Exogenous treatment of fruits with NO gas triggered the downregulation of CaG6PDH2, whereas the other genes were positively regulated. In-gel analysis using non-denaturing PAGE of a 50-75% ammonium-sulfate-enriched protein fraction from pepper fruits allowed for identifying two isozymes designated CaG6PDH I and CaG6PDH II, according to their electrophoretic mobility. In order to test the potential modulation of such pepper G6PDH isozymes, in vitro analyses of green pepper fruit samples in the presence of different compounds including NO donors (S-nitrosoglutathione and nitrosocysteine), peroxynitrite (ONOO-), a hydrogen sulfide (H2S) donor (NaHS, sodium hydrosulfide), and reducing agents such as reduced glutathione (GSH) and L-cysteine (L-Cys) were assayed. While peroxynitrite and the reducing compounds provoked a partial inhibition of one or both isoenzymes, NaHS exerted 100% inhibition of the two CaG6PDHs. Taken together these data provide the first data on the modulation of CaG6PDHs at gene and activity levels which occur in pepper fruit during ripening and after NO post-harvest treatment. As a consequence, this phenomenon may influence the NADPH availability for the redox homeostasis of the fruit and balance its active nitro-oxidative metabolism throughout the ripening process.

Identification and analysis of isoflavone reductase gene family in Gossypium hirsutum L

Isoflavone reductase (IFR) is a key enzyme controlling isoflavone synthesis and widely involved in response to various stresses. In this study, the IFR genes in four Gossypium species and other 7 species were identified and analyzed in the whole genome, and the physicochemical properties, gene structures, cis-acting elements, chromosomal locations, collinearity relationships and expression patterns of IFR genes were systematically analyzed. 28, 28, 14 and 15 IFR genes were identified in Gossypium hirsutum, Gossypium barbadense, Gossypium arboreum and Gossypium raimondii, respectively, which were divided into five clades according to the evolutionary tree and gene structure. Collinear analysis showed that segmental duplication and whole genome duplication were the main driving forces in the process of evolution, and most genes underwent pure selection. Gene structure analysis showed that IFR gene family was relatively conserved. Cis-element analysis of promoter showed that most GhIFR genes contain cis-elements related to abiotic stresses and plant hormones. Analysis of GhIFR gene expression under different stresses showed that GhIFR genes were involved in the response to drought, salt, heat and cold stresses through corresponding network mechanisms, especially GhIFR9A. Phenotypic analysis after silencing GhIFR9A gene by VIGS was shown that GhIFR9A gene was involved in the response to salt stress. This study laid a foundation for the subsequent functional study of cotton IFR genes.

Candidate gene based SSR and SNP markers for gynoecy in bitter gourd (Momordica charantia L.)

BACKGROUND

Even though the bitter gourd hybrids are shown to have significant heterosis for many of the economic traits, processes such as manual bagging and hand pollination make the hybrid seed production labour-intensive. Use of gynoecious line as female parent makes hybrid seed production more economical. This work was performed with the objective to identify the candidate gene based molecular markers for gynoecy in bitter gourd.

METHODS AND RESULTS

Seven putative genes for flowering and sex expression, isolated from the monoecious (MC-136) and gynoecious (KAU-MCGy-101) bitter gourd accessions, were sequence characterized. MADS-box transcription factor genes AG6 and McAG2 had nucleotide polymorphisms at five sites each and were potential candidates for marker development. An In/Del polymorphism of 48 bp ([TC]₂₄) in AG6 gene was used to develop an SSR marker and a transition mutation of [A/G] in this gene was used to develop a set of SNP markers. These markers have developed distinct polymorphism between the monoecious and gynoecious genotypes and were found suited for the marker assisted selection.

CONCLUSIONS

MADS box transcription factor genes AG6 and McAG2 are identified as candidates for sex expression in bitter gourd. Based on the InDels and transition in the intronic region of AG6, SSR marker BGAG6 and an SNP marker set segregating with the sex forms were developed. The markers have been validated using four other monoecious lines and are routinely used in our bitter gourd hybrid seed production programmes.

Genome-wide analysis of the SBT gene family involved in drought tolerance in cotton

The subtilisin-like proteases (SBTs) are a large family of serine peptidases that are unique to plants. Previous studies have shown that SBTs are associated with developmental processes and environmental responses. However, comprehensive identification and systematic analysis of the SBT family have not been conducted in cotton. We used bioinformatics methods to analyze the structural characteristics, phylogenetic relationships, gene structures, expression modes, evolutionary relationships, selection pressures and stress responses of SBT gene family members in upland cotton. In this study, we identified 120 and 112 SBTs in the tetraploid cotton species G. hirsutum and G. barbadense, while 67 and 69 SBTs were identified in the diploid species G. arboreum and G. raimondii, respectively; these SBTs were divided into five distinct subfamilies. We identified the SBT gene GhSBT27A, and explore its function through virus-induced gene silencing and transmission electron microscopy. These results suggested that the GhSBT27A gene was involved in the response to drought stress. These results lay a foundation for further study on the drought stress mechanism of cotton.

GhAAO2 was observed responding to NaHCO3 stress in cotton compared to AAO family genes

BACKGROUND

Abscisic acid (ABA) is an important stress hormone, the changes of abscisic acid content can alter plant tolerance to stress, abscisic acid is crucial for studying plant responses to abiotic stress. The abscisic acid aldehyde oxidase (AAO) plays a vital role in the final step in the synthesis of abscisic acid, therefore, understanding the function of AAO gene family is of great significance for plants to response to abiotic stresses.

RESULT

In this study, 6, 8, 4 and 4 AAO genes were identified in four cotton species. According to the structural characteristics of genes and the traits of phylogenetic tree, we divided the AAO gene family into 4 clades. Gene structure analysis showed that the AAO gene family was relatively conservative. The analysis of cis-elements showed that most AAO genes contained cis-elements related to light response and plant hormones. Tissue specificity analysis under NaHCO₃ stress showed that GhAAO2 gene was differentially expressed in both roots and leaves. After GhAAO2 gene silencing, the degree of wilting of seedlings was lighter than that of the control group, indicating that GhAAO2 could respond to NaHCO₃ stress.

CONCLUSIONS

In this study, the AAO gene family was analyzed by bioinformatics, the response of GhAAO gene to various abiotic stresses was preliminarily verified, and the function of the specifically expressed gene GhAAO2 was further verified. These findings provide valuable information for the study of potential candidate genes related to plant growth and stress.

Compensatory Effect of the ScGrf3-2R Gene in Semi-Dwarf Spring Triticale (x Triticosecale Wittmack)

The dwarfness in many triticale cultivars is provided by the dominant Ddw1 (Dominant dwarf 1) allele found in rye. However, along with conferring semi-dwarf phenotype to improve resistance to lodging, this gene also reduces grain size and weight and delays heading and flowering. Grf (Growth-regulating factors) genes are plant-specific transcription factors that regulate plant growth, including stem growth, in terms of length and thickness, and leaf and fruit size. In this work, we partially sequenced the rye gene ScGrf3 on chromosome 2R homologous to the wheat Grf3 gene, and found multiple polymorphisms in intron 3 and exon 4 complying with two alternative alleles (haplotypes ScGrf3-2Ra and ScGrf3-2Rb). For the identification of these, we developed a codominant PCR marker. Using a new marker, we studied the effect of ScGrf3-2R alleles in combination with the Ddw1 dwarf gene on economically valuable traits in F₄ and F₅ recombinant lines of spring triticale from the hybrid combination Valentin 90 x Dublet, grown in the Non-Chernozem zone for 2 years. Allele ScGrf3-2Ra was associated with greater thousand-grain weight, higher spike productivity, and earlier heading and flowering, which makes ScGrf3-2R a perspective compensator for negative effects of Ddw1 on these traits and increases prospects for its involvement in breeding semi-dwarf cultivars of triticale.

Alternative splicing as a key player in the fine-tuning of the immunity response in Arabidopsis

Plants, like animals, are constantly exposed to abiotic and biotic stresses, which often inhibit plant growth and development, and cause tissue damage, disease, and even plant death. Efficient and timely response to stress requires appropriate co- and posttranscriptional reprogramming of gene expression. Alternative pre-mRNA splicing provides an important layer of this regulation by controlling the level of factors involved in stress response and generating additional protein isoforms with specific features. Recent high-throughput studies have revealed that several defence genes undergo alternative splicing that is often affected by pathogen infection. Despite extensive work, the exact mechanisms underlying these relationships are still unclear, but the contribution of alternative protein isoforms to the defence response and the role of regulatory factors, including components of the splicing machinery, have been established. Modulation of gene expression in response to stress includes alternative splicing, chromatin remodelling, histone modifications, and nucleosome occupancy. How these processes affect plant immunity is mostly unknown, but these facets open new regulatory possibilities. Here we provide an overview of the current state of knowledge and recent findings regarding the growing importance of alternative splicing in plant response to biotic stress.

Identification of SNAT Family Genes Suggests GhSNAT3D Functional Reponse to Melatonin Synthesis Under Salinity Stress in Cotton

Serotonin N-acetyltransferase (SNAT) is a key enzyme in the biosynthesis of melatonin, and plays an important role in the regulation of melatonin synthesis. The study of SNAT is of great significance to understand the function of melatonin. In this study, we analyzed the structural characteristics, phylogenetic relationship, gene structure, expression pattern, evolutionary relationship and stress response of the members of the SNAT gene family in upland cotton through bioinformatics. A putative Serotonin n-acetyltransferase gene GhSNAT3D was identified, and preliminarily function of GhSNAT3D was verified by virus-induced gene silencing. We identified a total of 52 SNAT genes in the whole genome of G. hirsutum, and part of the GhSNATs were regulated by exogenous melatonin. The content of melatonin, antioxidant enzyme activity and Ca²⁺ content of GhSNAT3D gene silenced plants decreased, and the salt tolerance of GhSNAT3D gene silenced plants was reduced. Exogenous melatonin supplementation restored the salt tolerance of GhSNAT3D gene silenced plants. GhSNAT3D may interact with GhSNAT25D and ASMT to regulate melatonin synthesis. This study provided an important basis for further study on the regulation of melatonin in cotton against abiotic stress.

Genome-wide expression analysis of phospholipase A1 (PLA1) gene family suggests phospholipase A1-32 gene responding to abiotic stresses in cotton

Cotton is the most important crop for the production of natural fibres used in the textile industry. High salinity, drought, cold and high temperature represent serious abiotic stresses, which seriously threaten cotton production. Phospholipase A_S has an irreplaceable role in lipid signal transmission, growth and development and stress events. Phospholipase A can be divided into three families: PLA1, PLA2 and pPLA. Among them, the PLA1 family is rarely studied in plants. In order to study the potential functions of the PLA1 family in cotton, the bioinformatics analysis of the PLA1 family was correlated with cotton adversity, and tissue-specific analysis was performed. Explore the structure-function relationship of PLA1 members. It is found that the expression of GbPLA1-32 gene is affected by a variety of environmental stimuli, indicating that it plays a very important role in stress and hormone response, and closely associates the cotton adversity with this family. Through further functional verification, we found that virus-induced GbPLA1-32 gene silencing (VIGS) caused Gossypium barbadense to be sensitive to salt stress. This research provides an important basis for further research on the molecular mechanism of cotton resistance to abiotic stress.

In Silico Study of the RSH (RelA/SpoT Homologs) Gene Family and Expression Analysis in Response to PGPR Bacteria and Salinity in Brassica napus

Among several mechanisms involved in the plant stress response, synthesis of guanosine tetra and pentaphosphates (alarmones), homologous to the bacterial stringent response, is of crucial importance. Plant alarmones affect, among others, photosynthetic activity, metabolite accumulation, and nutrient remobilization, and thus regulate plant growth and development. The plant RSH (RelA/SpoT homolog) genes, that encode synthetases and/or hydrolases of alarmones, have been characterized in a limited number of plant species, e.g., Arabidopsis thaliana, Oryza sativa, and Ipomoea nil. Here, we used dry-to-wet laboratory research approaches to characterize RSH family genes in the polyploid plant Brassica napus. There are 12 RSH genes in the genome of rapeseed that belong to four types of RSH genes: 6 RSH1, 2 RSH2, 3 RSH3, and 1 CRSH. BnRSH genes contain 13-24 introns in RSH1, 2-6 introns in RSH2, 1-6 introns in RSH3, and 2-3 introns in the CRSH genes. In the promoter regions of the RSH genes, we showed the presence of regulatory elements of the response to light, plant hormones, plant development, and abiotic and biotic stresses. The wet-lab analysis showed that expression of BnRSH genes is generally not significantly affected by salt stress, but that the presence of PGPR bacteria, mostly of Serratia sp., increased the expression of BnRSH significantly. The obtained results show that BnRSH genes are differently affected by biotic and abiotic factors, which indicates their different functions in plants.

Plant Platforms for Efficient Heterologous Protein Production

Production of recombinant proteins is primarily established in cultures of mammalian, insect and bacterial cells. Concurrently, concept of using plants to produce high-value pharmaceuticals such as vaccines, antibodies, and dietary proteins have received worldwide attention. Newer technologies for plant transformation such as plastid engineering, agroinfiltration, magnifection, and deconstructed viral vectors have been used to enhance the protein production in plants along with the inherent advantage of speed, scale, and cost of production in plant systems. Production of therapeutic proteins in plants has now a more pragmatic approach when several plant-produced vaccines and antibodies successfully completed Phase I clinical trials in humans and were further scheduled for regulatory approvals to manufacture clinical grade products on a large scale which are safe, efficacious, and meet the quality standards. The main thrust of this review is to summarize the data accumulated over the last two decades and recent development and achievements of the plant derived therapeutics. It also attempts to discuss different strategies employed to increase the production so as to make plants more competitive with the established production systems in this industry.

Systematic Analysis and Functional Validation of Citrus Pectin Acetylesterases (CsPAEs) Reveals that CsPAE2 Negatively Regulates Citrus Bacterial Canker Development

The present study was designed to serve as a comprehensive analysis of Citrus sinensis (C. sinensis) pectin acetylesterases (CsPAEs), and to assess the roles of these PAEs involved in the development of citrus bacterial canker (CBC) caused by Xanthomonas citri subsp. citri (Xcc) infection. A total of six CsPAEs were identified in the genome of C. sinensis, with these genes being unevenly distributed across chromosomes 3, 6, and 9, and the unassembled scaffolds. A subset of CsPAEs were found to be involved in responses to Xcc infection. In particular, CsPAE2 was identified to be associated with such infections, as it was upregulated in CBC-susceptible variety Wanjincheng and inversely in CBC-resistant variety Calamondin. Transgenic citrus plants overexpressing CsPAE2 were found to be more susceptible to CBC, whereas the silencing of this gene was sufficient to confer CBC resistance. Together, these findings provide evolutionary insights into and functional information about the CsPAE family. This study also suggests that CsPAE2 is a potential candidate gene that negatively contributes to bacterial canker disease and can be used to breed CBC-resistant citrus plants.

Mild drought in the vegetative stage induces phenotypic, gene expression, and DNA methylation plasticity in Arabidopsis but no transgenerational effects

There is renewed interest in whether environmentally induced changes in phenotypes can be heritable. In plants, heritable trait variation can occur without DNA sequence mutations through epigenetic mechanisms involving DNA methylation. However, it remains unknown whether this alternative system of inheritance responds to environmental changes and if it can provide a rapid way for plants to generate adaptive heritable phenotypic variation. To assess potential transgenerational effects induced by the environment, we subjected four natural accessions of Arabidopsis thaliana together with the reference accession Col-0 to mild drought in a multi-generational experiment. As expected, plastic responses to drought were observed in each accession, as well as a number of intergenerational effects of the parental environments. However, after an intervening generation without stress, except for a very few trait-based parental effects, descendants of stressed and non-stressed plants were phenotypically indistinguishable irrespective of whether they were grown in control conditions or under water deficit. In addition, genome-wide analysis of DNA methylation and gene expression in Col-0 demonstrated that, while mild drought induced changes in the DNA methylome of exposed plants, these variants were not inherited. We conclude that mild drought stress does not induce transgenerational epigenetic effects.

Mild drought in the vegetative stage induces phenotypic, gene expression, and DNA methylation plasticity in Arabidopsis but no transgenerational effects

There is renewed interest in whether environmentally induced changes in phenotypes can be heritable. In plants, heritable trait variation can occur without DNA sequence mutations through epigenetic mechanisms involving DNA methylation. However, it remains unknown whether this alternative system of inheritance responds to environmental changes and if it can provide a rapid way for plants to generate adaptive heritable phenotypic variation. To assess potential transgenerational effects induced by the environment, we subjected four natural accessions of Arabidopsis thaliana together with the reference accession Col-0 to mild drought in a multi-generational experiment. As expected, plastic responses to drought were observed in each accession, as well as a number of intergenerational effects of the parental environments. However, after an intervening generation without stress, except for a very few trait-based parental effects, descendants of stressed and non-stressed plants were phenotypically indistinguishable irrespective of whether they were grown in control conditions or under water deficit. In addition, genome-wide analysis of DNA methylation and gene expression in Col-0 demonstrated that, while mild drought induced changes in the DNA methylome of exposed plants, these variants were not inherited. We conclude that mild drought stress does not induce transgenerational epigenetic effects.

Genome-wide analysis of WRKY transcription factors in Aquilaria sinensis (Lour.) Gilg

The WRKY proteins are a superfamily of transcription factor that regulate diverse developmental and physiological processes in plants. Completion of the whole-genome sequencing of Aquilaria sinensis allowed us to perform a genome-wide investigation for WRKY proteins. Here, we predicted 70 WRKY genes from the A. sinensis genome and undertaken a comprehensive bioinformatic analysis. Due to their diverse structural features, the 70 AsWRKY genes are classified into three main groups (group I-III), with five subgroups (IIa-IIe) in group II, except two belong to none of them. Distinct expression profiles of AsWRKYs with RNA sequencing data revealed their diverse expression patterns among different tissues and in the process of whole-tree-inducing agarwood formation. Based on the expression characteristics, we predict some AsWRKYs are pseudogenes, and some may be involved in the biosynthesis of agarwood sesquiterpenes as activators or repressors. Among the tested genes treated with MeJA and H₂O₂, most of them are induced by H₂O₂, but downregulated by MeJA, implying the complexity of their involvement in signal transduction regulation. Our results not only provide a basic platform for functional identification of WRKYs in A. sinensis but important clues for further analysis their regulation role in agarwood formation.

Identification and characterization of the Ipomoea nil RelA/SpoT Homologs (InRSHs) and potential directions of their transcriptional regulation

Although the stringent response has been known for more than half a century and has been well studied in bacteria, only the research of the past 19 years revealed that the homologous mechanism is conserved in plants. The plant RelA/SpoT Homolog (RSH) genes have been identified and characterized in a limited number of plant species, whereas products of their catalytic activities, (p)ppGpp (alarmones), have been shown to accumulate mainly in chloroplasts. Here, we identified full-length sequences of the Ipomoea nil RSH genes (InRSH1, InRSH2 and InCRSH), determined their copy number in the I. nil genome as well as the structural conservancy between InRSHs and their Arabidopsis and rice orthologs. We showed that InRSHs are differentially expressed in I. nil organ tissues and that only InRSH2 is upregulated in response to salt, osmotic and drought stress. Our results of the E. coli relA/spoT mutant complementation test suggest that InRSH1 is likely a (p)ppGpp hydrolase, InCRSH - synthetase and InRSH2 shows both activities. Finally, we referred our results to the recently published I. nil genomic and proteomic data and uncovered the complexity of the I. nil RSH family as well as potential ways of the InRSH transcriptional regulation.

Genome-wide characterization, molecular evolution and expression profiling of the metacaspases in potato (Solanum tuberosum L.)

Metacaspases are distant relatives of animal caspases found in plants, protozoa and fungi. Some recent studies have demonstrated that metacaspases are involved in regulating the developmental and environmentally induced programmed cell death in plants. In this study, we identified metacaspase gene family in potato (Solanum tuberosum L.) and analyzed their expression pattern in various developmental tissues and stress responses of plants. There were eight metacaspase genes identified in the Peptidase (Cysteine protease) C14 family and based upon sequence alignment and phylogenetic analysis, a systematic nomenclature of potato metacaspases (SotubMCs) has been proposed. Three of the eight candidate genes showing homology with Arabidopsis thaliana type I metacaspase, AtMC1 were given name SotubMC1, SotubMC2 and SotubMC3 as per the degree of relatedness. Similarly, the next three being homologous to A. thaliana type I metacaspase, AtMC3 were named SotubMC4, SotubMC5, and SotubMC6. The remaining two were named SotubMC7 and SotubMC8, showing significant similarity with type II metacaspases of A. thaliana, AtMC4 and AtMC9, respectively. Evolutionary divergence analysis of SotubMCs from its orthologs in seven other members of Solanaceae family as well as with A. thaliana, Vitis vinifera and Oryza sativa was also carried out. The d_N/d_S ratios of the orthologous pairs suggested the SotubMCs were under purifying (negative) selection in course of plant evolution. Splicing patterns of potato metacaspases were also analyzed. Amongst all SotubMCs, SotubMC2, SotubMC4, SotubMC6 and SotubMC7 genes appeared to produce multiple alternative spliced variants of different lengths. Furthermore using protein modeling tools, we have predicted the protein structure of identified metacaspases. The cis-regulatory elements analysis was also performed exhibiting the presence of development, stress and hormones related cis-elements in the promoter regions of the SotubMCs. This indicates that potato metacaspases might be playing important roles in the development, stress and hormone responsive pathways. Moreover, relative expression analysis of identified genes was carried out using qRT-PCR in various developmental tissues that also include stolons and tubers. The eight metacaspases showed differential expression in different tissues. Some of the tissues such as leaf undergoing senescence among different leaf developmental stages (immature, mature and senescent) displayed higher relative expression of some of the metacaspases, implying their involvement in leaf senescence. The expression pattern of SotubMCs under various abiotic, biotic and hormonal stresses was also analysed. The results showed that many members of the potato metacaspase gene family displayed differential expression patterns under various stress conditions. Taken together, the study could provide crucial resources for further investigations to understand the functional roles of the identified metacaspases in potato.

The Rapid Methylation of T-DNAs Upon Agrobacterium Inoculation in Plant Leaves

Agrobacterium tumefaciens has been foundational in the development of transgenic plants for both agricultural biotechnology and plant molecular research. However, the transformation efficiency and level of transgene expression obtained for any given construct can be highly variable. These inefficiencies often require screening of many lines to find one with consistent and heritable transgene expression. Transcriptional gene silencing is known to affect transgene expression, and is associated with DNA methylation, especially of cytosines in symmetric CG and CHG contexts. While the specificity, heritability and silencing-associated effects of DNA methylation of transgene sequences have been analyzed in many stably transformed plants, the methylation status of transgene sequences in the T-DNA during the transformation process has not been well-studied. Here we used agro-infiltration of the eGFP reporter gene in Nicotiana benthamiana leaves driven by either an AtEF1α-A4 or a CaMV-35S promoter to study early T-DNA methylation patterns of these promoter sequences. The T-DNA was examined by amplicon sequencing following sodium bisulfite treatment using three different sequencing platforms: Sanger sequencing, Ion Torrent PGM, and the Illumina MiSeq. Rapid DNA methylation was detectable in each promoter region just 2-3 days post-infiltration and the levels continued to rapidly accumulate over the first week, then steadily up to 21 days later. Cytosines in an asymmetric context (CHH) were the most heavily and rapidly methylated. This suggests that early T-DNA methylation may be important in determining the epigenetic and transcriptional fate of integrated transgenes. The Illumina MiSeq platform was the most sensitive and robust way of detecting and following the methylation profiles of the T-DNA promoters. The utility of the methods was then used to show a subtle but significant difference in promoter methylation during intron-mediated enhancement. In addition, the method was able to detect an increase in promoter methylation when the eGFP reporter gene was targeted by siRNAs generated by co-infiltration of a hairpin RNAi construct.

Sequence-guided approach to genotyping plant clones and species using polymorphic NB-ARC-related genes

Leveraging the heightened levels of polymorphism in NB-ARC-related protein encoding genes in higher plants, a bioinformatic pipeline was created to identify regions in this gene family from sequenced plant genomes that exhibit fragment length or single nucleotide differences in different accessions of the same species. Testing this approach with the aquatic plant Spirodela polyrhiza demonstrated its superior performance in comparison with currently available genotyping technologies based on PCR amplification. Rapid and economical genotyping tools that can reliably distinguish species and intraspecific variations in plants can be powerful tools for biogeographical and ecological studies. Clones of the cosmopolitan duckweed species, Spirodela polyrhiza, are difficult to distinguish morphologically due to their highly abbreviated architecture and inherently low levels of sequence variation. The use of plastidic markers and generic Amplification Fragment Length Polymorphism approaches have met with limited success in resolving clones of S. polyrhiza from diverse geographical locales. Using whole genome sequencing data from nine S. polyrhiza clones as a training set, we created an informatic pipeline to identify and rank polymorphic regions from nuclear-encoded NB-ARC-related genes to design markers for PCR, Sanger sequencing (barcoding), and fragment length analysis. With seven primer sets, we found 21 unique fingerprints from a set of 23 S. polyrhiza clones. However, three of these clones share the same fingerprint and are indistinguishable by these markers. These primer sets can also be used as interspecific barcoding tools to rapidly resolve S. polyrhiza from the closely related S. intermedia species without the need for DNA sequencing. Our work demonstrates a general approach of using hyper-polymorphic loci within genomes as a resource to produce facile tools that can have high resolving power for genotyping applications.

Identification and characterization of evolutionarily conserved alternative splicing events in a mangrove genus Sonneratia

Alternative splicing (AS), which produces multiple mRNA transcripts from a single gene, plays crucial roles in plant growth, development and environmental stress responses. Functional significances of conserved AS events among congeneric species have not been well characterized. In this study, we performed transcriptome sequencing to characterize AS events in four common species of Sonneratia, a mangrove genus excellently adaptive to intertidal zones. 7,248 to 12,623 AS events were identified in approximately 25% to 35% expressed genes in the roots of the four species. The frequency of AS events in Sonneratia was associated with genomic features, including gene expression level and intron/exon number and length. Among the four species, 1,355 evolutionarily conserved AS (ECAS) events were identified from 1,170 genes. Compared with non-ECAS events, ECAS events are of shorter length and less possibility to introduce premature stop codons (PTCs) and frameshifts. Functional annotations of the genes containing ECAS events showed that four of the 26 enriched Gene Ontology (GO) terms are involved in proton transport, signal transduction and carbon metabolism, and 60 genes from another three GO terms are implicated in responses to osmotic, oxidative and heat stresses, which may contribute to the adaptation of Sonneratia species to harsh intertidal environments.

Arabidopsis thaliana MLO genes are expressed in discrete domains during reproductive development

MLOs in Plant Reproduction. The MILDEW RESISTANCE LOCUS-O (MLO) protein family, comprised of 15 members, plays roles in diverse cell-cell communication processes such as powdery mildew susceptibility, root thigmomorphogenesis, and pollen tube reception. The NORTIA (NTA, AtMLO7) gene is expressed in the synergid cells of the female gametophyte where it functions in intercellular communication with the pollen tube. Discrepancies between previously published promoter::GUS and promoter::gene-GUS constructs expression patterns led us to explore the regulation of NTA expression. Here we found via NTA_pro::gNTA-GUS truncations that sequences within the NTA gene negatively regulate its expression in the stomata and carpel walls. This led to the hypothesis that other MLO family members may also have additional regulatory sequences within the gene. MLO_pro::gMLO-GUS constructs were examined for each family member focusing specifically on flowers in order to determine whether other MLOs could play a role in reproductive cell-cell communication. Notably, several MLOs were expressed in the pollen, in the stigma, in the pollinated style, and in the synergids and central cell. These findings indicate that other MLOs in addition to NTA could play a role in reproduction. Previous studies on the MLO family showed that phylogenetically related MLOs had redundant functions in powdery mildew infection and root thigmomorphogenesis; however, MLO expression in reproductive tissues did not strictly follow phylogenetic relationships, indicating that MLOs from different evolutionary origins may have been recruited for function in sexual reproduction.

Intron length polymorphism of β-tubulin genes of Aegilops biuncialis Vis

Intron-specific DNA polymorphism is present among plant β-tubulin gene family members and is considered to be one of the molecular markers based on the difference of tubulin introns length assayed both separately (TBP: 1st intron) or in combination (h-TBP: 1st and 2nd introns). These two approaches are possibly useful for wheat breeding programs, since TBP and h-TBP help to differentiate between the accessions of Aegilops biuncialis Vis., a wild relative of wheat. PCR-derived polymorphic fragments were resolved by PAGE electrophoresis. The length of amplicons varied significantly (395-3900 bp for TBP and 466-3440 bp for h-TBP), while the numbers of polymorphic bands were 21 for TBP and 23 for h-TBP, respectively. PIC mean value was circa 0.3. Dendrograms constructed on the basis of the Nei and Li coefficient with the high bootstrap support reveal a similar order of hierarchy for the samples analyzed using both methods. Thus, both techniques uncover DNA polymorphism level sufficiently high to distinguish different accessions of Ae. biuncialis Vis.

Rapid transporter regulation prevents substrate flow traffic jams in boron transport

Nutrient uptake by roots often involves substrate-dependent regulated nutrient transporters. For robust uptake, the system requires a regulatory circuit within cells and a collective, coordinated behaviour across the tissue. A paradigm for such systems is boron uptake, known for its directional transport and homeostasis, as boron is essential for plant growth but toxic at high concentrations. In Arabidopsis thaliana, boron uptake occurs via diffusion facilitators (NIPs) and exporters (BORs), each presenting distinct polarity. Intriguingly, although boron soil concentrations are homogenous and stable, both transporters manifest strikingly swift boron-dependent regulation. Through mathematical modelling, we demonstrate that slower regulation of these transporters leads to physiologically detrimental oscillatory behaviour. Cells become periodically exposed to potentially cytotoxic boron levels, and nutrient throughput to the xylem becomes hampered. We conclude that, while maintaining homeostasis, swift transporter regulation within a polarised tissue context is critical to prevent intrinsic traffic-jam like behaviour of nutrient flow.

Plant pectin acetylesterase structure and function: new insights from bioinformatic analysis

Background

Pectins are plant cell wall polysaccharides that can be acetylated on C2 and/or C3 of galacturonic acid residues. The degree of acetylation of pectin can be modulated by pectin acetylesterase (EC 3.1.1.6, PAE). The function and structure of plant PAEs remain poorly understood and the role of the fine-tuning of pectin acetylation on cell wall properties has not yet been elucidated.

Results

In the present study, a bioinformatic approach was used on 72 plant PAEs from 16 species among 611 plant PAEs available in plant genomic databases. An overview of plant PAE proteins, particularly Arabidopsis thaliana PAEs, based on phylogeny analysis, protein motif identification and modeled 3D structure is presented. A phylogenetic tree analysis using protein sequences clustered the plant PAEs into five clades. AtPAEs clustered in four clades in the plant kingdom PAE tree while they formed three clades when a phylogenetic tree was performed only on Arabidopsis proteins, due to isoform AtPAE9. Primitive plants that display a smaller number of PAEs clustered into two clades, while in higher plants, the presence of multiple members of PAE genes indicated a diversification of AtPAEs. 3D homology modeling of AtPAE8 from clade 2 with a human Notum protein showed an α/β hydrolase structure with the hallmark Ser-His-Asp of the active site. A 3D model of AtPAE4 from clade 1 and AtPAE10 from clade 3 showed a similar shape suggesting that the diversification of AtPAEs is unlikely to arise from the shape of the protein. Primary structure prediction analysis of AtPAEs showed a specific motif characteristic of each clade and identified one major group of AtPAEs with a signal peptide and one group without a signal peptide. A multiple sequence alignment of the putative plant PAEs revealed consensus sequences with important putative catalytic residues: Ser, Asp, His and a pectin binding site. Data mining of gene expression profiles of AtPAE revealed that genes from clade 2 including AtPAE7, AtPAE8 and AtPAE11, which are duplicated genes, are highly expressed during plant growth and development while AtPAEs without a signal peptide, including AtPAE2 and AtPAE4, are more regulated in response to plant environmental conditions.

Conclusion

Bioinformatic analysis of plant, and particularly Arabidopsis, AtPAEs provides novel insights, including new motifs that could play a role in pectin binding and catalytic sites. The diversification of AtPAEs is likely to be related to neofunctionalization of some AtPAE genes.

Electronic supplementary material

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

Cyclopiazonic Acid Is a Pathogenicity Factor for Aspergillus flavus and a Promising Target for Screening Germplasm for Ear Rot Resistance

Aspergillus flavus, an opportunistic pathogen, contaminates maize and other key crops with carcinogenic aflatoxins (AFs). Besides AFs, A. flavus makes many more secondary metabolites (SMs) whose toxicity in insects or vertebrates has been studied. However, the role of SMs in the invasion of plant hosts by A. flavus remains to be investigated. Cyclopiazonic acid (CPA), a neurotoxic SM made by A. flavus, is a nanomolar inhibitor of endoplasmic reticulum calcium ATPases (ECAs) and a potent inducer of cell death in plants. We hypothesized that CPA, by virtue of its cytotoxicity, may serve as a key pathogenicity factor that kills plant cells and supports the saprophytic life style of the fungus while compromising the host defense response. This proposal was tested by two complementary approaches. A comparison of CPA levels among A. flavus isolates indicated that CPA may be a determinant of niche adaptation, i.e., isolates that colonize maize make more CPA than those restricted only to the soil. Further, mutants in the CPA biosynthetic pathway are less virulent in causing ear rot than their wild-type parent in field inoculation assays. Additionally, genes encoding ECAs are expressed in developing maize seeds and are induced by A. flavus infection. Building on these results, we developed a seedling assay in which maize roots were exposed to CPA, and cell death was measured as Evans Blue uptake. Among >40 maize inbreds screened for CPA tolerance, inbreds with proven susceptibility to ear rot were also highly CPA sensitive. The publicly available data on resistance to silk colonization or AF contamination for many of the lines was also broadly correlated with their CPA sensitivity. In summary, our studies show that i) CPA serves as a key pathogenicity factor that enables the saprophytic life style of A. flavus and ii) maize inbreds are diverse in their tolerance to CPA. Taking advantage of this natural variation, we are currently pursuing both genome-wide and candidate gene approaches to identify novel components of maize resistance to Aspergillus ear rot.

A Leader Intron of a Soybean Elongation Factor 1A (eEF1A) Gene Interacts with Proximal Promoter Elements to Regulate Gene Expression in Synthetic Promoters

Introns, especially the first intron in the 5' untranslated region (5'UTR), can significantly impact gene expression via intron-mediated enhancement (IME). In this study, we demonstrate the leader intron of a soybean elongation factor 1A (eEF1A) gene (GmScreamM8) was essential for the high activity of the native promoter. Furthermore, the interaction of the GmScreamM8 leader intron with regulatory element sequences from several soybean eEF1A promoters was studied using synthetic promoters, which consisted of element tetramers upstream of a core promoter used to regulate a green fluorescent protein (gfp) reporter gene. Element tetramers, placed upstream of a GmScreamM8 core promoter, showed very high activity using both transient expression in lima bean cotyledons and stable expression in soybean hairy roots, only if the native leader intron was included, suggesting an interaction between intronic sequences and promoter elements. Partial deletions of the leader intron showed that a 222 bp intronic sequence significantly contributed to very high levels of GFP expression. Generation of synthetic intron variants with a monomeric or trimeric repeat of the 222 bp intronic sequence, yielded almost two-fold higher expression compared to the original intron, while partial deletion of the 222 bp intronic repeated sequence significantly decreased gene expression, indicating that this intronic sequence was essential for the intron-element interaction enhancement.

Genomic Dissection and Expression Profiling Revealed Functional Divergence in Triticum aestivum Leucine Rich Repeat Receptor Like Kinases (TaLRRKs)

The leucine rich repeat receptor like kinases (LRRK) constitute the largest subfamily of receptor like kinases (RLK), which play critical roles in plant development and stress responses. Herein, we identified 531 TaLRRK genes in Triticum aestivum (bread wheat), which were distributed throughout the A, B, and D sub-genomes and chromosomes. These were clustered into 233 homologous groups, which were mostly located on either homeologous chromosomes from various sub-genomes or in proximity on the same chromosome. A total of 255 paralogous genes were predicted which depicted the role of duplication events in expansion of this gene family. Majority of TaLRRKs consisted of trans-membrane region and localized on plasma-membrane. The TaLRRKs were further categorized into eight phylogenetic groups with numerous subgroups on the basis of sequence homology. The gene and protein structure in terms of exon/intron ratio, domains, and motifs organization were found to be variably conserved across the different phylogenetic groups/subgroups, which indicated a potential divergence and neofunctionalization during evolution. High-throughput transcriptome data and quantitative real time PCR analyses in various developmental stages, and biotic and abiotic (heat, drought, and salt) stresses provided insight into modus operandi of TaLRRKs during these conditions. Distinct expression of majority of stress responsive TaLRRKs homologous genes suggested their specified role in a particular condition. These results provided a comprehensive analysis of various characteristic features including functional divergence, which may provide the way for future functional characterization of this important gene family in bread wheat.

Genomic and Post-Translational Modification Analysis of Leucine-Rich-Repeat Receptor-Like Kinases in Brassica rapa

Among several receptor-like kinases (RLKs), leucine-rich-repeat receptor-like kinases (LRR-RLKs) are a major group of genes that play crucial roles in growth, development and stress responses in plant systems. Given that they have several functional roles, it is important to investigate their roles in Brassica rapa. In the present study, 303 LRR-RLKs were identified in the genome of B. rapa and comparative phylogenetic analysis of 1213 combined LRR-RLKs of B. rapa, Arabidopsis thaliana, Oryza sativa and Populus trichocarpa helped us to categorize the gene family into 15 subfamilies based on their sequence and structural similarities. The chromosome localizations of 293 genes allowed the prediction of duplicates, and motif conservation and intron/exon patterns showed differences among the B. rapa LRR-RLK (BrLRR-RLK) genes. Additionally, computational function annotation and expression analysis was used to predict their possible functional roles in the plant system. Biochemical results for 11 selected genes showed variations in phosphorylation activity. Interestingly, BrBAK1 showed strong auto-phosphorylation and trans-phosphorylation on its tyrosine and threonine residues compared with AtBAK1 in previous studies. The AtBAK1 receptor kinase is involved in plant growth and development, plant innate immunity, and programmed cell death, and our results suggest that BrBAK1 might also be involved in the same functions. Another interesting result was that BrBAK1, BrBRI1, BrPEPR1 and BrPEPR2 showed activity with both anti-phosphotyrosine and anti-phosphothreonine antibodies, indicating that they might have dual-specificity kinase activity. This study provides comprehensive results for the BrLRR-RLKs, revealing expansion of the gene family through gene duplications, structural similarities and variations among the genes, and potential functional roles according to gene ontology, transcriptome profiling and biochemical analysis.

Genome-wide cataloging and analysis of alternatively spliced genes in cereal crops

Background

Protein functional diversity at the post-transcriptional level is regulated through spliceosome mediated pre-mRNA alternative splicing (AS) events and that has been widely demonstrated to be a key player in regulating the functional diversity in plants. Identification and analysis of AS genes in cereal crop plants are critical for crop improvement and understanding regulatory mechanisms.

Results

We carried out the comparative analyses of the functional landscapes of the AS using the consensus assembly of expressed sequence tags and available mRNA sequences in four cereal plants. We identified a total of 8,734 in Oryza sativa subspecies (ssp) japonica, 2,657 in O. sativa ssp indica, 3,971 in Sorghum bicolor, and 10,687 in Zea mays AS genes. Among the identified AS events, intron retention remains to be the dominant type accounting for 23.5 % in S. bicolor, and up to 55.8 % in O. sativa ssp indica. We identified a total of 887 AS genes that were conserved among Z. mays, S. bicolor, and O. sativa ssp japonica; and 248 AS genes were found to be conserved among all four studied species or ssp. Furthermore, we identified 53 AS genes conserved with Brachypodium distachyon. Gene Ontology classification of AS genes revealed functional assignment of these genes in many biological processes with diverse molecular functions.

Conclusions

AS is common in cereal plants. The AS genes identified in four cereal crops in this work provide the foundation for further studying the roles of AS in regulation of cereal plant growth and development. The data can be accessed at Plant Alternative Splicing Database (http://proteomics.ysu.edu/altsplice/).

Electronic supplementary material

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

Intronic promoter-mediated feedback loop regulates bean PvSR2 gene expression

Intronic promoter usage is a widespread regulatory mechanism for increasing the diversity of gene expression regulation. However, its trans-acting mechanisms that modulate host gene expression are poorly characterized. We previously showed that bean (Phaseolus vulgaris) stress-related gene 2 (PvSR2) produces a metal-responsive element (MRE)-binding transcription factor (PvMTF-1) by using an intronic promoter. Interestingly, a functional MRE resides within the upstream PvSR2 promoter, implying that PvMTF-1 is likely to regulate its host gene expression via this MRE. Through DNA-protein interaction studies combined with an analysis of expression of PvMTF-1 in bean leaves, we now report that PvMTF-1 binds to this MRE and thereby up-regulates PvSR2 expression, thus constituting a positive regulatory feedback loop. Our findings provide evidence for a trans-acting mechanism of an intronic promoter in host gene expression, hence, adding a new layer of complexity in intronic promoter-mediated gene regulation. This study may also contribute to the understanding of the intron-mediated regulation of gene expression in trans.

Insights into the evolution and diversification of the AT-hook Motif Nuclear Localized gene family in land plants

BACKGROUND

Members of the ancient land-plant-specific transcription factor AT-Hook Motif Nuclear Localized (AHL) gene family regulate various biological processes. However, the relationships among the AHL genes, as well as their evolutionary history, still remain unexplored.

RESULTS

We analyzed over 500 AHL genes from 19 land plant species, ranging from the early diverging Physcomitrella patens and Selaginella to a variety of monocot and dicot flowering plants. We classified the AHL proteins into three types (Type-I/-II/-III) based on the number and composition of their functional domains, the AT-hook motif(s) and PPC domain. We further inferred their phylogenies via Bayesian inference analysis and predicted gene gain/loss events throughout their diversification. Our analyses suggested that the AHL gene family emerged in embryophytes and further evolved into two distinct clades, with Type-I AHLs forming one clade (Clade-A), and the other two types together diversifying in another (Clade-B). The two AHL clades likely diverged before the separation of Physcomitrella patens from the vascular plant lineage. In angiosperms, Clade-A AHLs expanded into 5 subfamilies; while, the ones in Clade-B expanded into 4 subfamilies. Examination of their expression patterns suggests that the AHLs within each clade share similar expression patterns with each other; however, AHLs in one monophyletic clade exhibit distinct expression patterns from the ones in the other clade. Over-expression of a Glycine max AHL PPC domain in Arabidopsis thaliana recapitulates the phenotype observed when over-expressing its Arabidopsis thaliana counterpart. This result suggests that the AHL genes from different land plant species may share conserved functions in regulating plant growth and development. Our study further suggests that such functional conservation may be due to conserved physical interactions among the PPC domains of AHL proteins.

CONCLUSIONS

Our analyses reveal a possible evolutionary scenario for the AHL gene family in land plants, which will facilitate the design of new studies probing their biological functions. Manipulating the AHL genes has been suggested to have tremendous effects in agriculture through increased seedling establishment, enhanced plant biomass and improved plant immunity. The information gleaned from this study, in turn, has the potential to be utilized to further improve crop production.

The most deeply conserved noncoding sequences in plants serve similar functions to those in vertebrates despite large differences in evolutionary rates

In vertebrates, conserved noncoding elements (CNEs) are functionally constrained sequences that can show striking conservation over >400 million years of evolutionary distance and frequently are located megabases away from target developmental genes. Conserved noncoding sequences (CNSs) in plants are much shorter, and it has been difficult to detect conservation among distantly related genomes. In this article, we show not only that CNS sequences can be detected throughout the eudicot clade of flowering plants, but also that a subset of 37 CNSs can be found in all flowering plants (diverging ∼170 million years ago). These CNSs are functionally similar to vertebrate CNEs, being highly associated with transcription factor and development genes and enriched in transcription factor binding sites. Some of the most highly conserved sequences occur in genes encoding RNA binding proteins, particularly the RNA splicing-associated SR genes. Differences in sequence conservation between plants and animals are likely to reflect differences in the biology of the organisms, with plants being much more able to tolerate genomic deletions and whole-genome duplication events due, in part, to their far greater fecundity compared with vertebrates.

Genome-wide analysis and expression profiling of the Solanum tuberosum aquaporins

Aquaporins belongs to the major intrinsic proteins involved in the transcellular membrane transport of water and other small solutes. A comprehensive genome-wide search for the homologues of Solanum tuberosum major intrinsic protein (MIP) revealed 41 full-length potato aquaporin genes. All potato aquaporins are grouped into five subfamilies; plasma membrane intrinsic proteins (PIPs), tonoplast intrinsic proteins (TIPs), NOD26-like intrinsic proteins (NIPs), small basic intrinsic proteins (SIPs) and x-intrinsic proteins (XIPs). Functional predictions based on the aromatic/arginine (ar/R) selectivity filters and Froger's positions showed a remarkable difference in substrate transport specificity among subfamilies. The expression pattern of potato aquaporins, examined by qPCR analysis, showed distinct expression profiles in various organs and tuber developmental stages. Furthermore, qPCR analysis of potato plantlets, subjected to various abiotic stresses revealed the marked effect of stresses on expression levels of aquaporins. Taken together, the expression profiles of aquaporins imply that aquaporins play important roles in plant growth and development, in addition to maintaining water homeostasis in response to environmental stresses.

Impact of an intragenic retrotransposon on the structural integrity and evolution of a major isoprenoid biosynthesis pathway gene in Hevea brasiliensis

Isoprenoids belong to a large family of structurally and functionally different natural compounds found universally from prokaryotes to higher animals and plants. In Hevea brasiliensis, the commercially important cis-polyisoprene (rubber) is synthesised as part of its defence mechanism in addition to other common isoprenoids like phytosterols, growth hormones etc. Farnesyl diphosphate synthase (FDPS) is a key enzyme in this process which catalyses the conversion of isoprene units into polyisoprene. Although prior sequence information is available, the structural variants of the FDPS gene presently existing in Hevea population are largely unknown. Since gene structure has a major role in gene regulation, extensive sequence analysis of this gene from different genotypes was carried out to identify the prevailing structural variants. We identified several SNPs and large indels which were associated with a partial transposable element (TE). Modification of key regulatory motifs and splice sites induced by the retroelement was also identified in the first intron. Screening of popular rubber clones, wild germplasm accessions and Hevea species revealed that the retroelement is responsible for the generation of new alleles with varying degrees of sequence homology. Segregation analysis of a progeny population confirmed that the alleles are not paralogs and are inherited in a Mendelian mode. Our findings suggest that the first intron of the FDPS gene has been subjected to various chromosomal rearrangements due to the interaction of a retrotransposon, resulting in novel alleles which may substantially contribute towards the evolution of this major gene in rubber. Moreover, the results indicate the possible existence of a retrotransposon-mediated epigenetic gene regulatory mechanism in Hevea.

Yin Yang 1 intronic binding sequences and splicing elicit intron-mediated enhancement of ubiquitin C gene expression

In a number of organisms, introns affect expression of the gene in which they are contained. Our previous studies revealed that the 5′-UTR intron of human ubiquitin C (UbC) gene is responsible for the boost of reporter gene expression and is able to bind, in vitro, Yin Yang 1 (YY1) trans-acting factor. In this work, we demonstrate that intact YY1 binding sequences are required for maximal promoter activity and YY1 silencing causes downregulation of luciferase mRNA levels. However, YY1 motifs fail to enhance gene expression when the intron is moved upstream of the proximal promoter, excluding the typical enhancer hypothesis and supporting a context-dependent action, like intron-mediated enhancement (IME). Yet, almost no expression is seen in the construct containing an unspliceable version of UbC intron, indicating that splicing is essential for promoter activity. Moreover, mutagenesis of YY1 binding sites and YY1 knockdown negatively affect UbC intron removal from both endogenous and reporter transcripts. Modulation of splicing efficiency by YY1 cis-elements and protein factor may thus be part of the mechanism(s) by which YY1 controls UbC promoter activity. Our data highlight the first evidence of the involvement of a sequence-specific DNA binding factor in IME.

Genome-wide survey of cold stress regulated alternative splicing in Arabidopsis thaliana with tiling microarray

Alternative splicing plays a major role in expanding the potential informational content of eukaryotic genomes. It is an important post-transcriptional regulatory mechanism that can increase protein diversity and affect mRNA stability. Alternative splicing is often regulated in a tissue-specific and stress-responsive manner. Cold stress, which adversely affects plant growth and development, regulates the transcription and splicing of plant splicing factors. This can affect the pre-mRNA processing of many genes. To identify cold regulated alternative splicing we applied Affymetrix Arabidopsis tiling arrays to survey the transcriptome under cold treatment conditions. A novel algorithm was used for detection of statistically relevant changes in intron expression within a transcript between control and cold growth conditions. A reverse transcription polymerase chain reaction (RT-PCR) analysis of a number of randomly selected genes confirmed the changes in splicing patterns under cold stress predicted by tiling array. Our analysis revealed new types of cold responsive genes. While their expression level remains relatively unchanged under cold stress their splicing pattern shows detectable changes in the relative abundance of isoforms. The majority of cold regulated alternative splicing introduced a premature termination codon (PTC) into the transcripts creating potential targets for degradation by the nonsense mediated mRNA decay (NMD) process. A number of these genes were analyzed in NMD-defective mutants by RT-PCR and shown to evade NMD. This may result in new and truncated proteins with altered functions or dominant negative effects. The results indicate that cold affects both quantitative and qualitative aspects of gene expression.

The production of male-sterile wheat plants through split barnase expression is promoted by the insertion of introns and flexible peptide linkers

The successful use of transgenic plants depends on the strong and stable expression of the heterologous genes. In this study, three introns (PSK7-i1 and PSK7-i3 from Petunia and UBQ10-i1 from Arabidopsis) were tested for their ability to enhance the tapetum-specific expression of a split barnase transgene. We also analyzed the effects of introducing multiple copies of flexible peptide linkers that bridged the fusion domains of the assembled protein. The barnase fragments were assembled into a functional cytotoxin via intein-mediated trans-splicing, thus leading to male sterility through pollen ablation. A total of 14 constructs carrying different combinations of introns and peptide linkers were transformed into wheat plants. The resulting populations (between 41 and 301 independent plants for each construct) were assayed for trait formation. Depending on which construct was used, there was an increase of up to fivefold in the proportion of plants exhibiting male sterility compared to the populations harboring unmodified constructs. Furthermore, the average barnase copy number in the plants displaying male sterility could be reduced. The metabolic profiles of male-sterile transgenic plants and non-transgenic plants were compared using gas chromatography-mass spectrometry. The profiles generated from leaf tissues displayed no differences, thus corroborating the anther specificity of barnase expression. The technical advances achieved in this study may be a valuable contribution for future improvement of transgenic crop systems.

Enhanced microRNA accumulation through stemloop-adjacent introns

MicroRNAs (miRNAs) originate from stemloop-forming precursor RNAs found in longer primary transcripts that often contain introns. We show that in plants, those introns, when located 3' of the stemloop, can promote mature miRNA accumulation, through a mechanism that likely operates at the level of miRNA processing or stability. Reversely, when miRNA production is reduced such as in dicer-like 1 mutants, splicing of introns that promote miRNA processing is considerably increased, pointing to a tight physical and temporal coordination of intron splicing and miRNA processing in plants. Our findings further suggest that miRNA transcripts without introns generated through alternative polyA-site usage might contribute to the differential adjustment of miRNA levels, possibly at a tissue-specific level.

Genome-wide landscape of alternative splicing events in Brachypodium distachyon

Recently, Brachypodium distachyon has emerged as a model plant for studying monocot grasses and cereal crops. Using assembled expressed transcript sequences and subsequent mapping to the corresponding genome, we identified 1219 alternative splicing (AS) events spanning across 2021 putatively assembled transcripts generated from 941 genes. Approximately, 6.3% of expressed genes are alternatively spliced in B. distachyon. We observed that a majority of the identified AS events were related to retained introns (55.5%), followed by alternative acceptor sites (16.7%). We also observed a low percentage of exon skipping (5.0%) and alternative donor site events (8.8%). The 'complex event' that consists of a combination of two or more basic splicing events accounted for ∼14.0%. Comparative AS transcript analysis revealed 163 and 39 homologous pairs between B. distachyon and Oryza sativa and between B. distachyon and Arabidopsis thaliana, respectively. In all, we found 16 AS transcripts to be conserved in all 3 species. AS events and related putative assembled transcripts annotation can be systematically browsed at Plant Alternative Splicing Database (http://proteomics.ysu.edu/altsplice/plant/).

Analysis of intron sequence features associated with transcriptional regulation in human genes

Although some preliminary work has revealed the potential transcriptional regulatory function of the introns in eukaryotes, additional evidences are needed to support this conjecture. In this study, we perform systemic analyses of the sequence characteristics of human introns. The results show that the first introns are generally longer and C, G and their dinucleotide compositions are over-represented relative to other introns, which are consistent with the previous findings. In addition, some new phenomena concerned with transcriptional regulation are found: i) the first introns are enriched in CpG islands; and ii) the percentages of the first introns containing TATA, CAAT and GC boxes are relatively higher than other position introns. The similar features of introns are observed in tissue-specific genes. The results further support that the first introns of human genes are likely to be involved in transcriptional regulation, and give an insight into the transcriptional regulatory regions of genes.

A candidate-gene association study for berry colour and anthocyanin content in Vitis vinifera L

Anthocyanin content is a trait of major interest in Vitis vinifera L. These compounds affect grape and wine quality, and have beneficial effects on human health. A candidate-gene approach was used to identify genetic variants associated with anthocyanin content in grape berries. A total of 445 polymorphisms were identified in 5 genes encoding transcription factors and 10 genes involved in either the biosynthetic pathway or transport of anthocyanins. A total of 124 SNPs were selected to examine association with a wide range of phenotypes based on RP-HPLC analysis and visual characterization. The phenotypes were total skin anthocyanin (TSA) concentration but also specific types of anthocyanins and relative abundance. The visual assessment was based on OIV (Organisation Internationale de la Vigne et du Vin) descriptors for berry and skin colour. The genes encoding the transcription factors MYB11, MYBCC and MYC(B) were significantly associated with TSA concentration. UFGT and MRP were associated with several different types of anthocyanins. Skin and pulp colour were associated with nine genes (MYB11, MYBCC, MYC(B), UFGT, MRP, DFR, LDOX, CHI and GST). Pulp colour was associated with a similar group of 11 genes (MYB11, MYBCC, MYC(B), MYC(A), UFGT, MRP, GST, DFR, LDOX, CHI and CHS(A)). Statistical interactions were observed between SNPs within the transcription factors MYB11, MYBCC and MYC(B). SNPs within LDOX interacted with MYB11 and MYC(B), while SNPs within CHI interacted with MYB11 only. Together, these findings suggest the involvement of these genes in anthocyanin content and on the regulation of anthocyanin biosynthesis. This work forms a benchmark for replication and functional studies.

T-DNA insertion mutants reveal complex expression patterns of the aldehyde dehydrogenase 3H1 locus in Arabidopsis thaliana

The Arabidopsis thaliana aldehyde dehydrogenase 3H1 gene (ALDH3H1; AT1G44170) belongs to family 3 of the plant aldehyde dehydrogenase superfamily. The full-length transcript of the corresponding gene comprises an open reading frame of 1583 bp and encodes a protein of 484 amino acid residues. Gene expression studies have shown that this transcript accumulates mainly in the roots of 4-week-old plants following abscisic acid, dehydration, and NaCl treatments. The current study provided experimental data that the ALDH3H1 locus generates at least five alternative transcript variants in addition to the previously described ALDH3H1 mRNA. The alternative transcripts accumulated in wild-type plants at a low level but were upregulated in a mutant that carried a T-DNA insertion in the first exon of the gene. Expression of the transcript isoforms involved alternative gene splicing combined with an alternative promoter. The transcript isoforms were differentially expressed in the roots and shoots and showed developmental stage- and tissue-specific expression patterns. These data support the hypothesis that alternative isoforms produced by gene splicing or alternative promoters regulate the abundance of the constitutively spliced and functional variants.

Intron retention regulates the expression of pectin methyl esterase inhibitor (Pmei) genes during wheat growth and development

Pectin is an important component of the plant cell wall and its remodelling occurs during normal plant growth or following stress responses. Pectin is secreted into the cell wall in a highly methyl-esterified form and subsequently de-methyl-esterified by pectin methyl esterase (PME), whose activity is controlled by the pectin methyl esterase inhibitor protein (PMEI). Cereal cell wall contains a low amount of pectin; nonetheless the level and pattern of pectin methyl esterification play a primary role during development or pathogen infection. Since few data are available on the role of PMEI in plant development and defence of cereal species, we isolated and characterised three Pmei genes (Tdpmei2.1, Tdpmei2.2 and Tdpmei3) and their encoded products in wheat. Sequence comparisons showed a low level of intra- and inter-specific sequence conservation of PMEIs. Tdpmei2.1 and Tdpmei2.2 share 94% identity at protein level, but only 20% identity with the product of Tdpmei3. All three Tdpmei genes code for functional inhibitors of plant PMEs and do not inhibit microbial PMEs or a plant invertase. RT-PCR analyses demonstrated, for the first time to our knowledge, that Pmei genes are regulated by intron retention. Processed and unprocessed transcripts of Tdpmei2.1 and Tdpmei2.2 accumulated in several organs, but anthers contained only mature transcripts. Tdpmei3 lacks introns and its transcript accumulated mainly in stem internodes. These findings suggest that products encoded by these Tdpmei genes control organ- or tissue-specific activity of specific PME isoforms in wheat.

Functional analysis of a rice late pollen-abundant UDP-glucose pyrophosphorylase (OsUgp2) promoter

OsUgp2, a rice UDP-glucose pyrophosphorylase gene, has previously been shown to preferentially express in maturing pollens and plays an important role in pollen starch accumulation. Here, a 1943 bp promoter fragment (P1943) of OsUgp2 was characterized by 5' deletion and gain-of-function experiments. P1943 and its 5' deletion derivatives (P1495, P1005, P665 and P159) were fused to GUS reporter gene and stably introduced into rice plants. Histochemical analyses of different tissues and pollens at different developmental stages of the transgenic plants showed that P1943 could only direct GUS expression in binucleate pollens. P1495 and P1005 could still drive GUS expression in binucleate pollens but at a lower level. On the other hand, neither P665 nor P159 transformant exhibited any GUS activity in pollens. Gain-of-function analyses showed that the region (-1005 to -665 relative to translation start site) combined with a minimal CaMV 35S promoter could direct GUS expression in pollens. Further analysis of 5' deletion truncated at -952, -847 and -740 delimited a 53 bp region (-1005 to -952) essential for pollen-specific expression. The 53 bp sequence contains two motifs of TTTCT and TTTC, which were known to be pollen-specific cis-elements. In addition, the same P1943-GUS fusion construct was introduced into tobacco to analyze its specificity in dicotyledon. Interestingly, the GUS expression pattern in transgenic tobacco was quite different from that in rice. High level of GUS expression was detected in mature pollens as well as leaves, roots, sepals and stigmas. These findings suggested a complicated transcriptional regulation of OsUgp2.