Divergence in function and expression of the NOD26-like intrinsic proteins in plants

Structural assessment of OsNIP2;1 highlighted critical residues defining solute specificity and functionality of NIP class aquaporins

INTRODUCTION

Nodulin-26-like intrinsic proteins (NIPs) are integral membrane proteins belonging to the aquaporin family, that facilitate the transport of neutral solutes across the bilayer. The OsNIP2;1 a member of NIP-III class of aquaporins is permeable to beneficial elements like silicon and hazardous arsenic. However, the atomistic cross-talk of these molecules traversing the OsNIP2;1 channel is not well understood.

OBJECTIVE

Due to the lack of genomic variation but the availability of high confidence crystal structure, this study aims to highlight structural determinants of metalloid permeation through OsNIP2;1.

METHODS

The molecular simulations, combined with site-directed mutagenesis were used to probe the role of specific residues in the metalloid transport activity of OsNIP2;1.

RESULTS

We drew energetic landscape of OsNIP2;1, for silicic and arsenous acid transport. Potential Mean Force (PMF) construction illuminate three prominent energetic barriers for metalloid passage through the pore. One corresponds to the extracellular molecular entry in the channel, the second located on ar/R filter, and the third size constriction in the cytoplasmic half. Comparative PMF for silicic acid and arsenous acid elucidate a higher barrier for silicic acid at the cytoplasmic constrict resulting in longer residence time for silicon. Furthermore, our simulation studies explained the importance of conserved residues in loop-C and loop-D with a direct effect on pore dynamics and metalloid transport. Next we assessed contribution of predicted key residues for arsenic uptake, by functional complementation in yeast. With the aim of reducing arsenic uptake while maintaining beneficial elements uptake, we identified novel OsNIP2;1 mutants with substantial reduction in arsenic uptake in yeast.

CONCLUSION

We provide a comprehensive assessment of pore lining residues of OsNIP2;1 with respect to metalloid uptake. The findings will expand mechanistic understanding of aquaporin's metalloid selectivity and facilitate variant interpretation to develop novel alleles with preference for beneficial metalloid species and reducing hazardous ones.

Integrative Analysis of miRNAs Involved in Fat Deposition in Different Pig Breeds

BACKGROUND

miRNAs are a set of small, noncoding RNAs that bind to partially complementary sequences on target mRNAs. This leads to the post-transcriptional regulation of gene expression. Many studies have shown that microRNAs play critical roles in adipose cell differentiation and fat metabolism. The aim of this study was to explore the regulatory functions of miRNAs in fat deposition for the prevention and therapy of lipid metabolism-related diseases.

METHODS

The significant differences in the fat deposition of Laiwu (LW) pigs and Large White (LY) pigs were studied. To investigate the genetic relationships of miRNAs that regulate fat deposition, we performed a genome-wide analysis of miRNAs derived from subcutaneous adipose tissue of LW and LY pigs using RNA-seq.

RESULTS

There were 39 known miRNAs and 56 novel miRNAs significantly differential expressed between the two breeds of pigs. In the analysis of the Gene Ontology and KEGG pathways, predicted targets of these differentially expressed miRNAs were involved in several fat-associated pathways, such as the peroxisome proliferator-activated receptor (PPAR), mitogen-activated protein kinases (MAPK) and Wnt signaling pathways. In addition, ssc-miR-133a-3p, ssc-miR-486 and ssc-miR-1 each had a great impact on the development of porcine subcutaneous fat through the PPAR signaling pathway.

CONCLUSIONS

We explored the role of differentially expressed miRNAs and studied the mechanisms of adipogenesis and fat deposition between two different pig breeds. In addition, these results also contribute to research relevant to human obesity.

Genetic Determinants of Fiber-Associated Traits in Flax Identified by Omics Data Integration

In this paper, we explore potential genetic factors in control of flax phenotypes associated with fiber by mining a collection of 306 flax accessions from the Federal Research Centre of the Bast Fiber Crops, Torzhok, Russia. In total, 11 traits were assessed in the course of 3 successive years. A genome-wide association study was performed for each phenotype independently using six different single-locus models implemented in the GAPIT3 R package. Moreover, we applied a multivariate linear mixed model implemented in the GEMMA package to account for trait correlations and potential pleiotropic effects of polymorphisms. The analyses revealed a number of genomic variants associated with different fiber traits, implying the complex and polygenic control. All stable variants demonstrate a statistically significant allelic effect across all 3 years of the experiment. We tested the validity of the predicted variants using gene expression data available for the flax fiber studies. The results shed new light on the processes and pathways associated with the complex fiber traits, while the pinpointed candidate genes may be further used for marker-assisted selection.

Three regions of the NIP5;1 promoter are required for expression in different cell types in Arabidopsis thaliana root

Arabidopsis thaliana NIP5;1, a boric acid diffusion facilitator, is involved in the acquisition of boron (B) from soil for growth under B limitation. AtNIP5;1 is expressed mainly in roots, where its expression is highest in the root cap and elongation zone. Here, we studied the role of the AtNIP5;1 promoter in the expression of this gene in roots. We fused a series of AtNIP5;1 promoter variants with deleted 5'-fragments to the GUS reporter gene and investigated the expression patterns by histochemical staining. We found that three regions of the AtNIP5;1 promoter are required for specific expression in the root cap and elongation zone (-880 to -863 bp from the translation start site), distal side of the differentiation zone (-747 to -722 bp), and basal side of the differentiation zone (-661 and -621 bp). The results suggest that at least three regions of the AtNIP5;1 promoter each confer different cell-type-specific expression.

Transcript profiles analysis of citrus aquaporins in response to fruit water loss during storage

Water loss is an essential factor that affects the maintenance of quality of citrus fruit during postharvest handling and storage. Aquaporins (AQPs) play an important role in the transport of water across membranes. However, the expression profiling of AQPs is incomplete for citrus fruits during storage. In this study, a post-harvest storage experiment was performed using sweet orange fruits to determine changes in water loss and fruit quality. Also, genome-wide expression analysis of CsAQP genes was carried out in fruit of different citrus varieties during storage. Low humidity storage conditions accelerated the postharvest water loss and texture decline and increased the TSS content in the fruit. A total of 39 non-redundant CsAQP genes were identified. A comprehensive analysis of these genes demonstrated that all AQPs had conserved filter motifs in the different citrus varieties examined. Moreover, multiple expression analysis revealed AQPs had complex expression profiles upon water loss in citrus fruit, being time-specific in tight-skin varieties (orange and pomelo varieties), tissue-specific between peel and pulp, and variety-specific between loose-skin (mandarin varieties) and tight-skin varieties (such as sweet orange and pummelo). These results indicated that the relative humidity in storage environment affected the postharvest water loss and quality of citrus fruit. Besides, the alternation in AQPs expression may partially account for the different water loss ratio in citrus varieties and the transfer of water between the peel and the pulp of citrus fruit during storage.

Aquaporins in regulation of plant protective responses to drought

Plasmolemma permeability is an integral indicator of the functional state of plant cells under stress. Aquaporins (AQPs), specialized transmembrane proteins that form water channels and play an important role in the adaptation of plants to adverse conditions and, in particular, to lack or excess of water, are involved in the formation of the response to drought. The main function of AQPs is to facilitate the movement of water across cell membranes and maintain aqueous cell homeostasis. Under stressful conditions, there is both an increase and decrease in the expression of individual aquaporin genes. Analysis of the data revealed differences in the expression of AQPs genes in stable and sensitive plant genotypes. It turned out that aquaporins in different stress-resistant varieties of the same species also respond differently to drought. The review provides brief information on the history of the discovery of aquaporins, the structure and function of these proteins, summarizes the latest information on the role of aquaporins in the regulation of metabolism and the response of plants to stressors, with particular emphasis on aquaporins in drought protection. The discovery and study of AQPs expands the possibilities of using genetic engineering methods for the selection of new plant species, in particular, more resistant to drought and salinization of the soil, as well as to increase their productivity. The use of aquaporins in biotechnology to improve drought resistance of various species has many prospects.

Grapevine aquaporins: Diversity, cellular functions, and ecophysiological perspectives

High-scored premium wines are typically produced under moderate drought stress, suggesting that the water status of grapevine is crucial for wine quality. Aquaporins greatly influence the plant water status by facilitating water diffusion across the plasma membrane in a tightly regulated manner. They adjust the hydraulic conductance of the plasma membrane rapidly and reversibly, which is essential in specific physiological events, including adaptation to soil water scarcity. The comprehension of the sophisticated plant-water relations at the molecular level are thus important to optimize agricultural practices or to assist plant breeding programs. This review explores the recent progresses in understanding the water transport in grapevine at the cellular level through aquaporins and its regulation. Important aspects, including aquaporin structure, diversity, cellular localization, transport properties, and regulation at the cellular and whole plant level are addressed. An ecophysiological perspective about the roles of grapevine aquaporins in plant response to drought stress is also provided.

Single-cell RNA sequencing of developing maize ears facilitates functional analysis and trait candidate gene discovery

Crop productivity depends on activity of meristems that produce optimized plant architectures, including that of the maize ear. A comprehensive understanding of development requires insight into the full diversity of cell types and developmental domains and the gene networks required to specify them. Until now, these were identified primarily by morphology and insights from classical genetics, which are limited by genetic redundancy and pleiotropy. Here, we investigated the transcriptional profiles of 12,525 single cells from developing maize ears. The resulting developmental atlas provides a single-cell RNA sequencing (scRNA-seq) map of an inflorescence. We validated our results by mRNA in situ hybridization and by fluorescence-activated cell sorting (FACS) RNA-seq, and we show how these data may facilitate genetic studies by predicting genetic redundancy, integrating transcriptional networks, and identifying candidate genes associated with crop yield traits.

Delineation of the Crucial Evolutionary Amino Acid Sites in Trehalose-6-Phosphate Synthase From Higher Plants

Trehalose-6-phosphate synthase (TPS) is a key enzyme in the biosynthesis of trehalose, with its direct product, trehalose-6-phosphate, playing important roles in regulating whole-plant carbohydrate allocation and utilization. Genes encoding TPS constitute a multigene family in which functional divergence appears to have occurred repeatedly. To identify the crucial evolutionary amino acid sites of TPS in higher plants, a series of bioinformatics tools were applied to investigate the phylogenetic relationships, functional divergence, positive selection, and co-evolution of TPS proteins. First, we identified 150 TPS genes from 13 higher plant species. Phylogenetic analysis placed these TPS proteins into 2 clades: clades A and B, of which clade B could be further divided into 4 subclades (B1-B4). This classification was supported by the intron-exon structures, with more introns present in clade A. Next, detection of the critical functionally divergent amino acid sites resulted in the isolation of a total of 286 sites reflecting nonredundant radical shifts in amino acid properties with a high posterior probability cutoff among subclades. In addition, positively selected sites were identified using a codon substitution model, from which 46 amino acid sites were isolated as exhibiting positive selection at a significant level. Moreover, 18 amino acid sites were highlighted both for functional divergence and positive selection; these may thus potentially represent crucial evolutionary sites in the TPS family. Further co-evolutionary analysis revealed 3 pairs of sites: 11S and 12H, 33S and 34N, and 109G and 110E as demonstrating co-evolution. Finally, the 18 crucial evolutionary amino acid sites were mapped in the 3-dimensional structure. A total of 77 sites harboring functionally and structurally important residues of TPS proteins were found by using the CLIPS-4D online tool; notably, no overlap was observed with the identified crucial evolutionary sites, providing positive evidence supporting their designation. A total of 18 sites were isolated as key amino acids by using multiple bioinformatics tools based on their concomitant functional divergence and positive selection. Almost all these key sites are located in 2 domains of this protein family where they exhibit no overlap with the structurally and functionally conserved sites. These results will provide an improved understanding of the complexity of the TPS gene family and of its function and evolution in higher plants. Moreover, this knowledge may facilitate the exploitation of these sites for protein engineering applications.

Molecular and Functional Characterization of Grapevine NIPs through Heterologous Expression in aqy-Null Saccharomyces cerevisiae

Plant Nodulin 26-like Intrinsic Proteins (NIPs) are multifunctional membrane channels of the Major Intrinsic Protein (MIP) family. Unlike other homologs, they have low intrinsic water permeability. NIPs possess diverse substrate selectivity, ranging from water to glycerol and to other small solutes, depending on the group-specific amino acid composition at aromatic/Arg (ar/R) constriction. We cloned three NIPs (NIP1;1, NIP5;1, and NIP6;1) from grapevine (cv. Touriga Nacional). Their expression in the membrane of aqy-null Saccharomyces cerevisiae enabled their functional characterization for water and glycerol transport through stopped-flow spectroscopy. VvTnNIP1;1 demonstrated high water as well as glycerol permeability, whereas VvTnNIP6;1 was impermeable to water but presented high glycerol permeability. Their transport activities were declined by cytosolic acidification, implying that internal-pH can regulate NIPs gating. Furthermore, an extension of C-terminal in VvTnNIP6;1M homolog, led to improved channel activity, suggesting that NIPs gating is putatively regulated by C-terminal. Yeast growth assays in the presence of diverse substrates suggest that the transmembrane flux of metalloids (As, B, and Se) and the heavy metal (Cd) are facilitated through grapevine NIPs. This is the first molecular and functional characterization of grapevine NIPs, providing crucial insights into understanding their role for uptake and translocation of small solutes, and extrusion of toxic compounds in grapevine.

Identification of CCCH Zinc Finger Proteins Family in Moso Bamboo (Phyllostachys edulis), and PeC3H74 Confers Drought Tolerance to Transgenic Plants

CCCH zinc finger proteins are a class of important zinc-finger transcription factors and have functions in various plant growth and stress responses, but their functions in moso bamboo (Phyllostachys edulis) are unclear. In this current study, we main investigated the structures, phylogenetic relationships, promoter elements and microsynteny of PeC3Hs. In this research, 119 CCCH zinc finger proteins (PeC3H1-119) identified genes in moso bamboo were divided into 13 subfamilies (A-M) based on phylogenetic analysis. Meanwhile, moso bamboo were treated with abscisic acid (ABA), methyl jasmonate (Me-JA) and gibberellic acid (GA) and 12 CCCH genes expression levels were assayed using qRT-PCR. In the three hormone treatments, 12 genes were up-regulated or down-regulated, respectively. In addition, PeC3H74 was localized on the cytomembrane, and it had self-activation activities. Phenotypic and physiological analysis showed that PeC3H74 (PeC3H74-OE) conferred drought tolerance of transgenic Arabidopsis, including H₂O₂ content, survival rate, electrolyte leakage as well as malondialdehyde content. Additionally, compared with wild-type plants, transgenic Arabidopsis thaliana seedling roots growth developed better under 10 μM ABA; Moreover, the stomatal of over-expressing PeC3H74 in Arabidopsis changed significantly under ABA treatment. The above results suggest that PeC3H74 was quickly screened by bioinformatics, and it may enhanced drought tolerance in plants through the ABA-dependent signaling pathway.

Comprehensive Analysis of SnRK Gene Family and their Responses to Salt Stress in Eucalyptus grandis

The sucrose non-fermentation-related protein kinase (SnRK) is a kind of Ser/Thr protein kinase, which plays a crucial role in plant stress response by phosphorylating the target protein to regulate the interconnection of various signaling pathways. However, little is known about the SnRK family in Eucalyptus grandis. Thirty-four putative SnRK sequences were identified in E. grandis and divided into three subgroups (SnRK1, SnRK2 and SnRK3) based on phylogenetic analysis and the type of domain. Chromosome localization showed that SnRK family members are unevenly distributed in the remaining 10 chromosomes, with the notable exception of chromosome 11. Gene structure analysis reveal that 10 of the 24 SnRK3 genes contained no introns. Moreover, conserved motif analyses showed that SnRK sequences belonged to the same subgroup that contained the same motif type of motif. The Ka/Ks ratio of 17 paralogues suggested that the EgrSnRK gene family underwent a purifying selection. The upstream region of EgrSnRK genes enriched with different type and numbers of cis-elements indicated that EgrSnRK genes are likely to play a role in the response to diverse stresses. Quantitative real-time PCR showed that the majority of the SnRK genes were induced by salt treatment. Genome-wide analyses and expression pattern analyses provided further understanding on the function of the SnRK family in the stress response to different environmental salt concentrations.

Integrated approaches to study the drought tolerance mechanism in hot pepper (Capsicum annuum L.)

Drought is one of the predominant abiotic stresses which have phenomenal impact on crop productivity. Alterations in aquaporin gene expressions are part of complex molecular responses by plant in response to drought. To better understand the role of aquaporins in economically important crop chilli (Capsicum annuum), drought induced gene expression of twelve aquaporins was determined in drought tolerant-KCa-4884 and drought susceptible-G-4 genotypes. Conjointly, the effect of drought on leaf water status and photosynthetic parameters were evaluated. Gene expression of all examined 12 aquaporins was up-regulated in KCa-4884 and in contrast, all the aquaporin genes were down-regulated in G-4 under drought stress. Significant variations among two chilli genotypes have been recorded in photosynthetic rate (P _n ), stomatal conductance (G _s ), and relative water content (RWC), sub-stomatal CO₂ concentration (C _i ). KCa-4884 revealed significantly high rates of P _n and RWC and decreased G _s under water deficit conditions providing evidence for superior drought adaptive strategies. Differences in physiological parameters illustrate prevention of water loss during drought. Up-regulation of aquaporins in drought tolerant genotype implicates their possible role in water relations and photosynthetic performance even under extended drought conditions.

The trihelix family of transcription factors: functional and evolutionary analysis in Moso bamboo (Phyllostachys edulis)

BACKGROUND

Trihelix transcription factors (TTFs) are photoresponsive proteins that have a representative three-helix structure (helix-loop-helix-loop-helix). Members of this gene family have been reported to play roles in many plant processes.

RESULTS

In this study, we performed a functional and evolutionary analysis of the TTFs in Moso bamboo (Phyllostachys edulis). A total of 35 genes were identified and grouped into five subfamilies (GT-1, GT-γ, GT-2, SIP1 and SH4) according to their structural properties. Gene structure analysis showed that most genes in the PeTTF family had fewer introns. A unique motif (Motif 16) to the GT-γ subfamily was identified by conserved motif analysis. Promoter analysis revealed various cis-acting elements related to plant growth and development, abiotic and biotic stresses, and phytohormone responses. Data for the 35 Moso bamboo TTF genes were used to generate heat maps, which indicated that these genes were expressed in different tissues or developmental stages. Most of the TTF genes identified here had high expression in leaves and panicles according to the expression profile analysis. The expression levels of the TTF members in young leaves were studied using quantitative real-time PCR to determine their tissue specificity and stress-related expression patterns to help functionally characterize individual members.

CONCLUSIONS

The results indicated that members of the TTF gene family may be involved in plant responses to stress conditions. Additionally, PeTTF29 was shown to be located in the nucleus by subcellular localization analysis and to have transcriptional activity in a transcriptional activity assay. Our research provides a comprehensive summary of the PeTTF gene family, including functional and evolutionary perspectives, and provides a basis for functionally characterizing these genes.

Study on aquaporins of Setaria italica suggests the involvement of SiPIP3;1 and SiSIP1;1 in abiotic stress response

Aquaporins are versatile proteins involved in several biological as well as molecular functions, and they have been extensively studied in various plant systems. Increasing evidences indicate their role in biotic and abiotic stresses, and therefore, studying these proteins in a naturally stress-tolerant crop would provide further insights into the roles of this important protein family. Given this, the present study was performed in foxtail millet (Setaria italica), a model plant for studying biofuel, stress tolerance, and C₄ photosynthetic traits. The study identified 12 plasma membrane intrinsic proteins (PIPs), 11 tonoplast intrinsic proteins (TIPs), 13 NOD26-like intrinsic proteins (NIPs), and 3 small basic intrinsic proteins (SIPs) in foxtail millet. The identified proteins and their corresponding genes were characterized using in silico approaches such as chromosomal localization, analysis of gene and protein properties, phylogenetic analysis, promoter analysis, and RNA-seq-derived expression profiling. The candidate genes identified through these analyses were studied for their expression in response to abiotic stresses (dehydration, salinity, and heat) as well as hormone treatments (abscisic acid, methyl jasmonate, and salicylic acid) in two contrasting cultivars of foxtail millet. The study showed that SiPIP3;1 and SiSIP1;1 were differentially expressed in both the cultivars in response to stress and hormone treatments. Overexpression of these genes in a heterologous yeast system also demonstrated that the transgenic cells were able to tolerate dehydration as well as salt stress which suggests the involvement of these proteins in the tolerance mechanism. Overall, the present study provides insights into structure and organization of the aquaporin gene family in foxtail millet and highlights the potential candidate genes for further functional characterizations.

Genome-wide identification and expression analysis of expansin gene family in common wheat (Triticum aestivum L.)

Background

Expansin loosens plant cell walls and involves in cell enlargement and various abiotic stresses. Plant expansin superfamily contains four subfamilies: α-expansin (EXPA), β-expansin (EXPB), expansin-like A (EXLA), and expansin-like B (EXLB). In this work, we performed a comprehensive study on the molecular characterization, phylogenetic relationship and expression profiling of common wheat (Triticum aestivum) expansin gene family using the recently released wheat genome database (IWGSC RefSeq v1.1 with a coverage rate of 94%).

Results

Genome-wide analysis identified 241 expansin genes in the wheat genome, which were grouped into three subfamilies (EXPA, EXPB and EXLA) by phylogenetic tree. Molecular structure analysis showed that wheat expansin gene family showed high evolutionary conservation although some differences were present in different subfamilies. Some key amino acid sites that contribute to functional divergence, positive selection, and coevolution were detected. Evolutionary analysis revealed that wheat expansin gene superfamily underwent strong positive selection. The transcriptome map and qRT-PCR analysis found that wheat expansin genes had tissue/organ expression specificity and preference, and generally highly expressed in the roots. The expression levels of some expansin genes were significantly induced by NaCl and polyethylene glycol stresses, which was consistent with the differential distribution of the cis-elements in the promoter region.

Conclusions

Wheat expansin gene family showed high evolutionary conservation and wide range of functional divergence. Different selection constraints may influence the evolution of the three expansin subfamilies. The different expression patterns demonstrated that expansin genes could play important roles in plant growth and abiotic stress responses. This study provides new insights into the structures, evolution and functions of wheat expansin gene family.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5455-1) contains supplementary material, which is available to authorized users.

Basic helix-loop-helix gene family: Genome wide identification, phylogeny, and expression in Moso bamboo

Studies have shown that basic helix-loop-helix (bHLH) transcription factors play important roles in plant growth and survival, and response to various biotic/abiotic stresses. We identified a total of 448 bHLH genes. These genes were classified into 21 bHLH subfamilies, and most genes in a given subfamily had similar gene structures and conserved motifs. We identified 176 homologous pairs in the three species. We calculated Ka, Ks, and Ka/Ks to analyze the replication relationships among the three species. Multiple sequence analysis revealed that the PebHLH genes had the distinct bHLH structure. The gene ontology annotation analysis showed that the PebHLH genes had many molecular functions. Promoter cis-element analysis revealed that most of the PebHLH genes contained cis-elements that can respond to various biotic/abiotic stress-related events. The tissue expression patterns of the PebHLH genes indicated that most members were expressed in leaves, roots, and stems. Quantitative real-time PCR analysis showed that 21 selected PebHLH genes were differentially regulated after abscisic acid, drought, and methyl jasmonate treatments. This study has laid the basis for studying the functions of AtbHLH, OsbHLH, and PebHLH genes, and will contribute to future studies of the functions of bHLH genes in other plant species.

Assessing the correlation of genotypic and phenotypic responses of indica rice varieties under drought stress

Drought is one of the severe abiotic stress that affects the productivity of rice, an important staple crop that is consumed all over the world. The traits responsible for enhancing or adapting drought resistance in rice plants can be selected and studied to improve their growth under stress conditions. Experiments have been conducted on indica rice varieties comprising Sahabhagidhan as drought tolerant variety and IR64, MTU1010 categorized as drought sensitive varieties. Various root related biochemical and morphological traits such as root length, relative water content (RWC), xylem number, xylem area, proline content, and malondialdehyde content have been investigated for a comparative study of the plant response to drought stress in different rice varieties. The results of differential root transcriptome analysis have revealed that there is a notable difference in gene expression of OsPIP2;5 and OsNIP2;1 in various indica varieties of rice at different time periods of stress. The present work aims at assessing the correlation between genotypic and phenotypic traits that can contribute towards the emerging field of rice phenomics.

Genome-wide identification of PHD-finger genes and expression pattern analysis under various treatments in moso bamboo (Phyllostachys edulis)

Plant homeodomain (PHD)-finger proteins are a class of important zinc-finger transcription factors responsible for regulating transcription and the chromatin state and responsive to various stresses. The family genes have been reported in many plants, but there is little information about PHD-finger genes in moso bamboo. In this study, 60 PHD-finger genes (PePHD1-60) were identified in moso bamboo and classified into 11 subfamilies (A-K) based on phylogenetic analysis. Gene structure and conserved motif analysis showed that these genes contained different numbers of introns but had similar motif organizations within each subfamily. Multiple sequence alignment revealed that the PHD-finger proteins possessed conserved structural domain sequences. In addition, the family underwent purifying selection during evolution and experienced a large-scale duplication event around 7.69-15.4 million years ago. Most importantly, the expression profiles of young leaves (YL), mature leaves (L), roots (R), stems (S), shoots (Sh) and rhizomes (Rh) displayed that they might involve in the formation of these tissues. Based on promoter analysis of 16 putative stress-related genes, quantitative real-time PCR assays were performed using moso bamboo leaves and showed that these genes were differentially regulated under abscisic acid (ABA), drought, low temperature and NaCl treatments. Therefore, the results reveal that PePHD genes play crucial roles in organ formation and response to multiple environmental stress conditions of moso bamboo, which will make for further function analysis of PHD-finger genes in plants.

Genome-wide analysis of VQ motif-containing proteins in Moso bamboo (Phyllostachys edulis)

29 Moso bamboo VQ proteins were genome-wide identified for the first time, and bioinformatics analysis was performed to investigate phylogenetic relationships and evolutionary divergence. The qRT-PCR data show that PeVQ genes response to different stress treatments. Accumulating evidence suggests that VQ motif-containing proteins in rice (Oryza sativa), Arabidopsis (Arabidopsis thaliana), and maize (Zea mays) play fundamental roles in response to various biotic and abiotic stresses. However, little is known about the functions of VQ family proteins in Moso bamboo (Phyllostachys edulis). In this study, we performed a genome-wide bioinformatic analysis and expression profiling of PeVQ genes. A total of 29 VQ genes was identified and divided into seven subgroups (I-VII) based on phylogenetic analysis. Gene structure and conserved motif analysis revealed that 25 of 29 VQ genes contained no introns. Multiple sequence alignment showed that Moso bamboo VQ motif-containing proteins contained five variations of the conserved motif. The time of duplication and divergence of Moso bamboo from rice and maize was calculated using K _s analysis. A heat map was generated using microarray data from 29 Moso bamboo VQ genes suggesting that these genes were expressed in different tissues or developmental stages. Quantitative real-time PCR (qRT-PCR) and promoter analysis indicated that PeVQ genes were differentially regulated following treatment with polyethylene glycol, abscisic acid and salicylic acid. Our results provide a solid foundation for further research of the specific functions of VQ motif-containing proteins in Moso bamboo.

Comparative transcriptome profiling of upland (VS16) and lowland (AP13) ecotypes of switchgrass

KEY MESSAGE

Transcriptomes of two switchgrass genotypes representing the upland and lowland ecotypes will be key tools in switchgrass genome annotation and biotic and abiotic stress functional genomics. Switchgrass (Panicum virgatum L.) is an important bioenergy feedstock for cellulosic ethanol production. We report genome-wide transcriptome profiling of two contrasting tetraploid switchgrass genotypes, VS16 and AP13, representing the upland and lowland ecotypes, respectively. A total of 268 million Illumina short reads (50 nt) were generated, of which, 133 million were obtained in AP13 and the rest 135 million in VS16. More than 90% of these reads were mapped to the switchgrass reference genome (V1.1). We identified 6619 and 5369 differentially expressed genes in VS16 and AP13, respectively. Gene ontology and KEGG pathway analysis identified key genes that regulate important pathways including C4 photosynthesis, photorespiration and phenylpropanoid metabolism. A series of genes (33) involved in photosynthetic pathway were up-regulated in AP13 but only two genes showed higher expression in VS16. We identified three dicarboxylate transporter homologs that were highly expressed in AP13. Additionally, genes that mediate drought, heat, and salinity tolerance were also identified. Vesicular transport proteins, syntaxin and signal recognition particles were seen to be up-regulated in VS16. Analyses of selected genes involved in biosynthesis of secondary metabolites, plant-pathogen interaction, membrane transporters, heat, drought and salinity stress responses confirmed significant variation in the relative expression reflected in RNA-Seq data between VS16 and AP13 genotypes. The phenylpropanoid pathway genes identified here are potential targets for biofuel conversion.

Phylogenetic Analysis, Lineage-Specific Expansion and Functional Divergence of seed dormancy 4-Like Genes in Plants

The rice gene seed dormancy 4 (OsSdr4) functions in seed dormancy and is a major factor associated with pre-harvest sprouting (PHS). Although previous studies of this protein family were reported for rice and other species, knowledge of the evolution of genes homologous to OsSdr4 in plants remains inadequate. Fifty four Sdr4-like (hereafter designated Sdr4L) genes were identified in nine plant lineages including 36 species. Phylogenetic analysis placed these genes in eight subfamilies (I-VIII). Genes from the same lineage clustered together, supported by analysis of conserved motifs and exon-intron patterns. Segmental duplications were present in both dicot and monocot clusters, while tandemly duplicated genes occurred only in monocot clusters indicating that both tandem and segmental duplications contributed to expansion of the grass I and II subfamilies. Estimation of the approximate ages of the duplication events indicated that ancestral Sdr4 genes evolved from a common angiosperm ancestor, about 160 million years ago (MYA). Moreover, diversification of Sdr4L genes in mono and dicot plants was mainly associated with genome-wide duplication and speciation events. Functional divergence was observed in all subfamily pairs, except IV/VIIIa. Further analysis indicated that functional constraints between subfamily pairs I/II, I/VIIIb, II/VI, II/VIIIb, II/IV, and VI/VIIIb were statistically significant. Site and branch-site model analyses of positive selection suggested that these genes were under strong adaptive selection pressure. Critical amino acids detected for both functional divergence and positive selection were mostly located in the loops, pointing to functional importance of these regions in this protein family. In addition, differential expression studies by transcriptome atlas of 11 Sdr4L genes showed that the duplicated genes may have undergone divergence in expression between plant species. Our findings showed that Sdr4L genes are functionally divergent and positively selected. These may contribute to further functional analysis and molecular evolution of Sdr4L gene families in land plants.

CsNIP2;1 is a Plasma Membrane Transporter from Cucumis sativus that Facilitates Urea Uptake When Expressed in Saccharomyces cerevisiae and Arabidopsis thaliana

Urea is an important source of nitrogen (N) for the growth and development of plants. It occurs naturally in soils, is the major N source in agricultural fertilizers and is an important N metabolite in plants. Therefore, the identification and characterization of urea transporters in higher plants is important for the fundamental understanding of urea-based N nutrition in plants and for designing novel strategies for improving the N-use efficiency of urea based-fertilizers. Progress in this area, however, is hampered due to scarce knowledge of plant urea transporters. From what is known, urea uptake from the soil into plant roots is mediated by two types of transporters: the major intrinsic proteins (MIPs) and the DUR3 orthologs, mediating low- and high-affinity urea transport, respectively. Here we characterized a MIP family member from Cucumis sativus, CsNIP2;1, with regard to its contribution to urea transport. We show that CsNIP2;1 is a plasma membrane transporter that mediates pH-dependent urea uptake when expressed in yeast. We also found that ectopic expression of CsNIP2;1 improves growth of wild-type Arabidopsis thaliana and rescues growth and development of the atdur3-3 mutant on medium with urea as the sole N source. In addition, CsNIP2;1 is transcriptionally up-regulated by N deficiency, urea and NO3 (-). These data and results from the analyses of the pattern of CsNIP2;1 expression in A. thaliana and cucumber suggest that CsNIP2;1 might be involved in multiple steps of urea-based N nutrition, including urea uptake and internal transport during N remobilization throughout seed germination and N delivery to developing tissues.

The divergence and positive selection of the plant-specific BURP-containing protein family

BURP domain-containing proteins belong to a plant-specific protein family and have diverse roles in plant development and stress responses. However, our understanding about the genetic divergence patterns and evolutionary rates of these proteins remain inadequate. In this study, 15 plant genomes were explored to elucidate the genetic origins, divergence, and functions of these proteins. One hundred and twenty-five BURP protein-encoding genes were identified from four main plant lineages, including 13 higher plant species. The absence of BURP family genes in unicellular and multicellular algae suggests that this family (1) appeared when plants shifted from relatively stable aquatic environments to land, where conditions are more variable and stressful, and (2) is critical in the adaptation of plants to adverse environments. Promoter analysis revealed that several responsive elements to plant hormones and external environment stresses are concentrated in the promoter region of BURP protein-encoding genes. This finding confirms that these genes influence plant stress responses. Several segmentally and tandem-duplicated gene pairs were identified from eight plant species. Thus, in general, BURP domain-containing genes have been subject to strong positive selection, even though these genes have conformed to different expansion models in different species. Our study also detected certain critical amino acid sites that may have contributed to functional divergence among groups or subgroups. Unexpectedly, all of the critical amino acid residues of functional divergence and positive selection were exclusively located in the C-terminal region of the BURP domain. In conclusion, our results contribute novel insights into the genetic divergence patterns and evolutionary rates of BURP proteins.

Molecular Characterization and Expression Profiling of NAC Transcription Factors in Brachypodium distachyon L

NAC (NAM, ATAF1/2, CUC2) transcription factors are involved in regulating plant developmental processes and response to environmental stresses. Brachypodium distachyon is an emerging model system for cereals, temperate grasses and biofuel crops. In this study, a comprehensive investigation of the molecular characterizations, phylogenetics and expression profiles under various abiotic stresses of the NAC gene family in Brachypodium distachyon was performed. In total, 118 BNAC genes in B. distachyon were identified, of which 22 (18.64%) were tandemly duplicated and segmentally duplicated, respectively. The Bayesian phylogenetic inference using Markov Chain Monte Carlo (MCMC) algorithms showed that they were divided into two clades and fourteen subfamilies, supported by similar motif compositions within one subfamily. Some critical amino acids detected using DIVERGE v3.0 might contribute to functional divergence among subfamilies. The different exon-intron organizations among subfamilies revealed structural differentiation. Promoter sequence predictions showed that the BNAC genes were involved in various developmental processes and diverse stress responses. Three NAC domain-encoding genes (BNAC012, BNAC078 and BNAC108), orthologous of NAC1, were targeted by five miRNA164 (Bdi-miR164a-c, e, f), suggesting that they might function in lateral organ enlargement, floral development and the responses to abiotic stress. Eleven (~9.32%) BNAC proteins containing α-helical transmembrane motifs were identified. 23 representative BNAC genes were analyzed by quantitative real-time PCR, showing different expression patterns under various abiotic stresses, of which 18, 17 and 11 genes were up-regulated significantly under drought, H2O2 and salt stresses, respectively. Only four and two genes were up-regulated under cold and cadmium stresses, respectively. Dynamic transcriptional expression analysis revealed that six genes showed constitutive expression and period-specific expression. The current results provide novel insights into the structure and function of the plant NAC gene family.

Metalloido-porins: Essentiality of Nodulin 26-like intrinsic proteins in metalloid transport

Metalloids are a group of physiologically important elements ranging from the essential to the highly toxic. Arsenic, antimony, germanium, and tellurium are highly toxic to plants themselves and to consumers of metalloid-contaminated plants. Boron, silicon, and selenium fulfill essential or beneficial functions in plants. However, when present at high concentrations, boron and selenium cause toxicity symptoms that are detrimental to plant fitness and yield. Consequently, all plants require efficient membrane transport systems to control the uptake and extrusion of metalloids into or out of the plant and their distribution within the plant body. Several Nodulin 26-like intrinsic proteins (NIPs) that belong to the aquaporin plant water channel protein family facilitate the diffusion of uncharged metalloid species. Genetic, physiological, and molecular evidence is that NIPs from primitive to higher plants not only transport all environmentally important metalloids, but that these proteins have a major role in the uptake, translocation, and extrusion of metalloids in plants. As most of the metalloid-permeable NIP aquaporins are impermeable or are poorly permeable to water, these NIP channel proteins should be considered as physiologically essential metalloido-porins.

Functional characterization of FaNIP1;1 gene, a ripening-related and receptacle-specific aquaporin in strawberry fruit

Strawberry fruit (Fragaria × ananassa) is a soft fruit with high water content at ripe stage (more than 90% of its fresh weight). Aquaporins play an important role in plant water homeostasis, through the facilitation of water transport and solutes. We report the role played by FaNIP1;1 in the receptacle ripening process. The analysis by qRT-PCR of FaNIP1;1 showed that this gene is mainly expressed in fruit receptacle and has a ripening-related expression pattern that was accompanied by an increase in both the abscisic acid and water content of the receptacle throughout fruit ripening. Moreover, FaNIP1;1 was induced in situations of water deficit. Additionally, we show that FaNIP1;1 expression was positively regulated by abscisic acid and negatively regulated by auxins. The water transport capacity of FaNIP1;1 was determined by a stopped-flow spectroscopy in yeast over-expressing FaNIP1;1. Glycerol, H2O2 and boron transport were also demonstrated in yeast. On the other hand, GFP-FaNIP1;1 fusion protein was located in plasma membrane. In conclusion, FaNIP1;1 seems to play an important role increasing the plasma membrane permeability, that allows the water accumulation in the strawberry fruit receptacle throughout the ripening process.

Evolution and functional diversity of aquaporins

In this review, we provide a brief synopsis of the evolution and functional diversity of the aquaporin gene superfamily in prokaryotic and eukaryotic organisms. Based upon the latest data, we discuss the expanding list of molecules shown to permeate the central pore of aquaporins, and the unexpected diversity of water channel genes in Archaea and Bacteria. We further provide new insight into the origin by horizontal gene transfer of plant glycerol-transporting aquaporins (NIPs), and the functional co-option and gene replacement of insect glycerol transporters. Finally, we discuss the origins of four major grades of aquaporins in Eukaryota, together with the increasing repertoires of aquaporins in vertebrates.

A comprehensive survey of the grapevine VQ gene family and its transcriptional correlation with WRKY proteins

WRKY proteins are a class of transcription factors (TFs) involved in the regulation of various physiological processes, including the plant response to biotic and abiotic stresses. Recent studies in Arabidopsis have revealed that some WRKY TFs interact with a class of proteins designed as VQ proteins because of their typical conserved motif (FxxxVQxLTG). So far, no information is available about the genomic organization and the function of VQ motif-containing protein in grapevine (Vitis vinifera L). In the current study, we analyzed the 12X V1 prediction of the nearly homozygous PN40024 genotype identifying up to 18 predicted VQ genes (VvVQ). VvVQs phylogenetic and bioinformatic analyses indicated that the intron-exon structures and motif distribution are highly divergent between different members of the grapevine VQ family. Moreover, the analysis of the V. vinifera cv. Corvina expression atlas revealed a tissue- and stage-specific expression of several members of the family which also showed a significant correlation with WRKY TFs. Grapevine VQ genes also exhibited altered expression in response to drought, powdery mildew infection, salicylic acid (SA) and ethylene (ETH) treatments. The present study represents the first characterization of VQ genes in a grapevine genotype and it is a pivotal foundation for further studies aimed at functionally characterizing this mostly unknown grapevine multigenic family.

Unusual tandem expansion and positive selection in subgroups of the plant GRAS transcription factor superfamily

BACKGROUND

GRAS proteins belong to a plant transcription factor family that is involved with multifarious roles in plants. Although previous studies of this protein family have been reported for Arabidopsis, rice, Chinese cabbage and other species, investigation of expansion patterns and evolutionary rate on the basis of comparative genomics in different species remains inadequate.

RESULTS

A total of 289 GRAS genes were identified in Arabidopsis, B. distachyon, rice, soybean, S. moellendorffii, and P. patens and were grouped into seven subfamilies, supported by the similarity of their exon-intron patterns and structural motifs. All of tandem duplicated genes were found in group II except one cluster of rice, indicating that tandem duplication greatly promoted the expansion of group II. Furthermore, segment duplications were mainly found in the soybean genome, whereas no single expansion pattern dominated in other plant species indicating that GRAS genes from these five species might be subject to a more complex evolutionary mechanism. Interestingly, branch-site model analyses of positive selection showed that a number of sites were positively selected under foreground branches I and V. These results strongly indicated that these groups were experiencing higher positive selection pressure. Meanwhile, the site-specific model revealed that the GRAS genes were under strong positive selection in P. patens. DIVERGE v2.0 was used to detect critical amino acid sites, and the results showed that the shifted evolutionary rate was mainly attributed to the functional divergence between the GRAS genes in the two groups. In addition, the results also demonstrated the expression divergence of the GRAS duplicated genes in the evolution. In short, the results above provide a solid foundation for further functional dissection of the GRAS gene superfamily.

CONCLUSIONS

In this work, differential expression, evolutionary rate, and expansion patterns of the GRAS gene family in the six species were predicted. Especially, tandem duplication events played an important role in expansion of group II. Together, these results contribute to further functional analysis and the molecular evolution of the GRAS gene superfamily.

Screening and identification of soybean seed-specific genes by using integrated bioinformatics of digital differential display, microarray, and RNA-seq data

Soybean is one of the most economically important crops in the world. Soybean seeds have abundant protein and lipid content and very high economic value. In this study, a total of 184 seed-specific genes were obtained using online microarray databases, DDD, and RNA-seq data. The reported seed-specific genes in soybean and the 184 seed-specific genes analyzed in this paper were compared. Of the screened genes, 26 were common to both previous reports and the current screening. Meanwhile, 90 of the 184 genes have homologous counterparts in Arabidopsis, among which 24 have seed-specific expression, as indicated by microarray data for Arabidopsis. Furthermore, promoter analysis showed that almost all seed-specific genes contain at least one seed specific-related element. Seed-specific element Skn-1 motif exists in most, if not all, of the seed-specific genes screened. Five genes were randomly selected from 184 soybean seed specific gene pool and their expressions were quantified using quantitative real time polymerase chain reaction (qRT-PCR) to further confirm the specificity of the screened genes. The results indicated that all five genes showed seed-specific expression. Moreover, the identification of genes with seed-specific expression screened in this study provides information valuable to the in-depth study of soybean.

Soybean (Glycine max) expansin gene superfamily origins: segmental and tandem duplication events followed by divergent selection among subfamilies

BACKGROUND

Expansins are plant cell wall loosening proteins that are involved in cell enlargement and a variety of other developmental processes. The expansin superfamily contains four subfamilies; namely, α-expansin (EXPA), β-expansin (EXPB), expansin-like A (EXLA), and expansin-like B (EXLB). Although the genome sequencing of soybeans is complete, our knowledge about the pattern of expansion and evolutionary history of soybean expansin genes remains limited.

RESULTS

A total of 75 expansin genes were identified in the soybean genome, and grouped into four subfamilies based on their phylogenetic relationships. Structural analysis revealed that the expansin genes are conserved in each subfamily, but are divergent among subfamilies. Furthermore, in soybean and Arabidopsis, the expansin gene family has been mainly expanded through tandem and segmental duplications; however, in rice, segmental duplication appears to be the dominant process that generates this superfamily. The transcriptome atlas revealed notable differential expression in either transcript abundance or expression patterns under normal growth conditions. This finding was consistent with the differential distribution of the cis-elements in the promoter region, and indicated wide functional divergence in this superfamily. Moreover, some critical amino acids that contribute to functional divergence and positive selection were detected. Finally, site model and branch-site model analysis of positive selection indicated that the soybean expansin gene superfamily is under strong positive selection, and that divergent selection constraints might have influenced the evolution of the four subfamilies.

CONCLUSION

This study demonstrated that the soybean expansin gene superfamily has expanded through tandem and segmental duplication. Differential expression indicated wide functional divergence in this superfamily. Furthermore, positive selection analysis revealed that divergent selection constraints might have influenced the evolution of the four subfamilies. In conclusion, the results of this study contribute novel detailed information about the molecular evolution of the expansin gene superfamily in soybean.

The large soybean (Glycine max) WRKY TF family expanded by segmental duplication events and subsequent divergent selection among subgroups

BACKGROUND

WRKY genes encode one of the most abundant groups of transcription factors in higher plants, and its members regulate important biological process such as growth, development, and responses to biotic and abiotic stresses. Although the soybean genome sequence has been published, functional studies on soybean genes still lag behind those of other species.

RESULTS

We identified a total of 133 WRKY members in the soybean genome. According to structural features of their encoded proteins and to the phylogenetic tree, the soybean WRKY family could be classified into three groups (groups I, II, and III). A majority of WRKY genes (76.7%; 102 of 133) were segmentally duplicated and 13.5% (18 of 133) of the genes were tandemly duplicated. This pattern was not apparent in Arabidopsis or rice. The transcriptome atlas revealed notable differential expression in either transcript abundance or in expression patterns under normal growth conditions, which indicated wide functional divergence in this family. Furthermore, some critical amino acids were detected using DIVERGE v2.0 in specific comparisons, suggesting that these sites have contributed to functional divergence among groups or subgroups. In addition, site model and branch-site model analyses of positive Darwinian selection (PDS) showed that different selection regimes could have affected the evolution of these groups. Sites with high probabilities of having been under PDS were found in groups I, II c, II e, and III. Together, these results contribute to a detailed understanding of the molecular evolution of the WRKY gene family in soybean.

CONCLUSIONS

In this work, all the WRKY genes, which were generated mainly through segmental duplication, were identified in the soybean genome. Moreover, differential expression and functional divergence of the duplicated WRKY genes were two major features of this family throughout their evolutionary history. Positive selection analysis revealed that the different groups have different evolutionary rates. Together, these results contribute to a detailed understanding of the molecular evolution of the WRKY gene family in soybean.

Herbivory of maize by southern corn rootworm induces expression of the major intrinsic protein ZmNIP1;1 and leads to the discovery of a novel aquaporin ZmPIP2;8

Aquaporins channel water and other neutral molecules through cell membranes. Aquaporin gene expression is subject to transcriptional control and can be modulated by factors affecting water balance such as salt, abscisic acid and drought. During infestation of maize by southern corn rootworm (SCR), an insect that chews into and significantly damages maize roots, three maize aquaporins were differentially expressed upon prolonged infestation. Using a brief infestation of maize roots ZmNIP1;1 transcript abundance again increased under infestation while expression of a new aquaporin, ZmPIP2;8 and ZmTIP2;2 expression did not change. Since ZmPIP2;8 has not been described previously, the deduced protein sequence was analyzed in silico and found to contain the hallmarks of plant aquaporins, with a predicted protein structure similar to other functionally characterized PIP2s. NIPs characterized to date have been implicated in facilitating the movement of a variety of small molecules, while TIPs and PIPs often have the capacity to facilitate trans-membrane movement of water. Functional assays (using heterologous expression in Xenopus laevis oocytes) of ZmTIP2;2 and ZmPIP2;8 confirmed that these aquaporins demonstrate water channel capacity.

Expression dynamics and genome distribution of osmoprotectants in soybean: identifying important components to face abiotic stress

BACKGROUND

Despite the importance of osmoprotectants, no previous in silico evaluation of high throughput data is available for higher plants. The present approach aimed at the identification and annotation of osmoprotectant-related sequences applied to short transcripts from a soybean HT-SuperSAGE (High Throughput Super Serial Analysis of Gene Expression; 26-bp tags) database, and also its comparison with other transcriptomic and genomic data available from different sources.

METHODS

A curated set of osmoprotectants related sequences was generated using text mining and selected seed sequences for identification of the respective transcripts and proteins in higher plants. To test the efficiency of the seed sequences, these were aligned against four HT-SuperSAGE contrasting libraries generated by our group using soybean tolerant and sensible plants against water deficit, considering only differentially expressed transcripts (p ≤ 0.05). Identified transcripts from soybean and their respective tags were aligned and anchored against the soybean virtual genome.

RESULTS

The workflow applied resulted in a set including 1,996 seed sequences that allowed the identification of 36 differentially expressed genes related to the biosynthesis of osmoprotectants [Proline (P5CS: 4, P5CR: 2), Trehalose (TPS1: 9, TPPB: 1), Glycine betaine (BADH: 4) and Myo-inositol (MIPS: 7, INPS1: 8)], also mapped in silico in the soybean genome (25 loci). Another approach considered matches using Arabidopsis full length sequences as seed sequences, and allowed the identification of 124 osmoprotectant-related sequences, matching ~10.500 tags anchored in the soybean virtual chromosomes. Osmoprotectant-related genes appeared clustered in all soybean chromosomes, with higher density in some subterminal regions and synteny among some chromosome pairs.

CONCLUSIONS

Soybean presents all searched osmoprotectant categories with some important members differentially expressed among the comparisons considered (drought tolerant or sensible vs. control; tolerant vs. sensible), allowing the identification of interesting candidates for biotechnological inferences. The identified tags aligned to corresponding genes that matched 19 soybean chromosomes. Osmoprotectant-related genes are not regularly distributed in the soybean genome, but clustered in some regions near the chromosome terminals, with some redundant clusters in different chromosomes indicating their involvement in previous duplication and rearrangements events. The seed sequences, transcripts and map represent the first transversal evaluation for osmoprotectant-related genes and may be easily applied to other plants of interest.

Characterization of physiological response and identification of associated genes under heat stress in rice seedlings

Global warming, which is caused by greenhouse gas emissions, makes food crops more vulnerable to heat stress. Understanding the heat stress-related mechanisms in crops and classifying heat stress-related genes can increase our knowledge in heat-resistant molecular biology and propel developments in molecular design breeding, which can help rice cope with unfavorable temperatures. In this study, we carried out a physiological analysis of rice plants after heat stress. The results show a dramatic increase in malondialdehyde contents and SOD activities. We successfully isolated 11 heat-related rice genes with known function annotation through DNSH, which is an improved SSH method for screening long cDNA fragments. The reanalysis of microarray data from public database revealed that all these genes displayed various expression patterns after heat stress, drought, cold and salt. Quantitative real-time reverse transcription PCR was also performed to validate the expression of these genes after heat stress. The expressions in 10 genes were all significantly changed except for contig 77, which is a CBL-interacting protein kinase. Several reports have been published about the members of the same gene family.

Two putative-aquaporin genes are differentially expressed during arbuscular mycorrhizal symbiosis in Lotus japonicus

BACKGROUND

Arbuscular mycorrhizas (AM) are widespread symbioses that provide great advantages to the plant, improving its nutritional status and allowing the fungus to complete its life cycle. Nevertheless, molecular mechanisms that lead to the development of AM symbiosis are not yet fully deciphered. Here, we have focused on two putative aquaporin genes, LjNIP1 and LjXIP1, which resulted to be upregulated in a transcriptomic analysis performed on mycorrhizal roots of Lotus japonicus.

RESULTS

A phylogenetic analysis has shown that the two putative aquaporins belong to different functional families: NIPs and XIPs. Transcriptomic experiments have shown the independence of their expression from their nutritional status but also a close correlation with mycorrhizal and rhizobial interaction. Further transcript quantification has revealed a good correlation between the expression of one of them, LjNIP1, and LjPT4, the phosphate transporter which is considered a marker gene for mycorrhizal functionality. By using laser microdissection, we have demonstrated that one of the two genes, LjNIP1, is expressed exclusively in arbuscule-containing cells. LjNIP1, in agreement with its putative role as an aquaporin, is capable of transferring water when expressed in yeast protoplasts. Confocal analysis have demonstrated that eGFP-LjNIP1, under its endogenous promoter, accumulates in the inner membrane system of arbusculated cells.

CONCLUSIONS

Overall, the results have shown different functionality and expression specificity of two mycorrhiza-inducible aquaporins in L. japonicus. One of them, LjNIP1 can be considered a novel molecular marker of mycorrhizal status at different developmental stages of the arbuscule. At the same time, LjXIP1 results to be the first XIP family aquaporin to be transcriptionally regulated during symbiosis.

Discovery of a multigene family of aquaporin silicon transporters in the primitive plant Equisetum arvense

Plants benefit greatly from silicon (Si) absorption provided that they contain Si transporters. The latter have recently been identified in the roots of some higher plants known to accumulate high concentrations of Si, and all share a high level of sequence identity. In this study, we searched for transporters in the primitive vascular plant Equisetum arvense (horsetail), which is a valuable but neglected model plant for the study of Si absorption, as it has one of the highest Si concentrations in the plant kingdom. Our initial attempts to identify Si transporters based on sequence homology with transporters from higher plants proved unsuccessful, suggesting a divergent structure or property in horsetail transporters. Subsequently, through sequencing of the horsetail root transcriptome and a search using amino acid sequences conserved in plant aquaporins, we were able to identify a multigene family of aquaporin Si transporters. Comparison of known functional domains and phylogenetic analysis of sequences revealed that the horsetail proteins belong to a different group than higher-plant Si transporters. In particular, the newly identified proteins contain a STAR pore as opposed to the GSGR pore common to all previously identified Si transporters. In order to determine its functionality, the proteins were heterologously expressed in both Xenopus oocytes and Arabidopsis, and the results showed that the horsetail proteins are extremely efficient a transporting Si. These findings offer new insights into the elusive properties of Si and its absorption by plants.

Mining plant genome browsers as a means for efficient connection of physical, genetic and cytogenetic mapping: An example using soybean

Physical maps are important tools to uncover general chromosome structure as well as to compare different plant lineages and species, helping to elucidate genome structure, evolution and possibilities regarding synteny and colinearity. The increasing production of sequence data has opened an opportunity to link information from mapping studies to the underlying sequences. Genome browsers are invaluable platforms that provide access to these sequences, including tools for genome analysis, allowing the integration of multivariate information, and thus aiding to explain the emergence of complex genomes. The present work presents a tutorial regarding the use of genome browsers to develop targeted physical mapping, providing also a general overview and examples about the possibilities regarding the use of Fluorescent In Situ Hybridization (FISH) using bacterial artificial chromosomes (BAC), simple sequence repeats (SSR) and rDNA probes, highlighting the potential of such studies for map integration and comparative genetics. As a case study, the available genome of soybean was accessed to show how the physical and in silico distribution of such sequences may be compared at different levels. Such evaluations may also be complemented by the identification of sequences beyond the detection level of cytological methods, here using members of the aquaporin gene family as an example. The proposed approach highlights the complementation power of the combination of molecular cytogenetics and computational approaches for the anchoring of coding or repetitive sequences in plant genomes using available genome browsers, helping in the determination of sequence location, arrangement and number of repeats, and also filling gaps found in computational pseudochromosome assemblies.

Genetic variation and association mapping of silica concentration in rice hulls using a germplasm collection

An association analysis on the genetic variability for silica concentration in rice hulls was performed using a "Mini-Core" set of 174 accessions representative of the germplasm diversity found in the USDA world collection of rice. Hull silica concentration was determined in replicated trials conducted in two southern states in the USA and was analyzed for its association with 164 genome-wide DNA markers. Among the accessions, the average silica concentration ranged from 120 to 251 mg g(-1). Ample variation was seen within each of the five sub-populations of rice, as well as the 14 geographic regions that the accessions originated from. There was also an effect due to location and accession × location (G × E) interaction demonstrating the importance of assessing silica concentration across multiple environments. Twelve markers on ten chromosomes were significantly associated with hull silica concentration. Six markers (RM5644, RM5371, RM1335, RM283, RM263, and RM178) corroborated quantitative trait locus for silica concentration identified in other mapping studies. Our results provide germplasm and genetic markers that will assist breeding efforts to develop cultivars that have either high or low hull silica concentration. High silica hulls are good raw material for silica based industrial compounds, while low silica hulls are more biodegradable.

Phylogenetic perspectives on the origins of nodulation

Recent refinements to the phylogeny of rosid angiosperms support the conclusion that nodulation has evolved several times in the so-called N(2)-fixing clade (NFC), and provide dates for these origins. The hypothesized predisposition that enabled the evolution of nodulation occurred approximately 100 million years ago (MYA), was retained in the various lineages that radiated rapidly shortly thereafter, and was functional in its non-nodulation role for at least an additional 30 million years in each nodulating lineage. Legumes radiated rapidly shortly after their origin approximately 60 MYA, and nodulation most likely evolved several times during this radiation. The major lineages of papilionoid legumes diverged close to the time of origin of nodulation, accounting for the diversity of nodule biology in the group. Nodulation symbioses exemplify the concept of "deep homology," sharing various homologous components across nonhomologous origins of nodulation, largely due to recruitment from existing functions, notably the older arbuscular mycorrhizal symbiosis. Although polyploidy may have played a role in the origin of papilionoid legume nodules, it did not do so in other legumes, nor did the prerosid whole-genome triplication lead directly to the predisposition of nodulation.

Arabidopsis thaliana NIP7;1: an anther-specific boric acid transporter of the aquaporin superfamily regulated by an unusual tyrosine in helix 2 of the transport pore

Plant nodulin-26 intrinsic proteins (NIPs) are members of the aquaporin superfamily that serve as multifunctional transporters of uncharged metabolites. In Arabidopsis thaliana, a specific NIP pore subclass, known as the NIP II proteins, is represented by AtNIP5;1 and AtNIP6;1, which encode channel proteins expressed in roots and leaf nodes, respectively, that participate in the transport of the critical cell wall nutrient boric acid. Modeling of the protein encoded by the AtNIP7;1 gene shows that it is a third member of the NIP II pore subclass in Arabidopsis. However, unlike AtNIP5;1 and AtNIP6;1 proteins, which form constitutive boric acid channels, AtNIP7;1 forms a channel with an extremely low intrinsic boric acid transport activity. Molecular modeling and molecular dynamics simulations of AtNIP7;1 suggest that a conserved tyrosine residue (Tyr81) located in transmembrane helix 2 adjacent to the aromatic arginine (ar/R) pore selectivity region stabilizes a closed pore conformation through interaction with the canonical Arg220 in ar/R region. Substitution of Tyr81 with a Cys residue, characteristic of established NIP boric acid channels, results in opening of the AtNIP7;1 pore that acquires a robust, transport activity for boric acid as well as other NIP II test solutes (glycerol and urea). Substitution of a Phe for Tyr81 also opens the channel, supporting the prediction from MD simulations that hydrogen bond interaction between the Tyr81 phenol group and the ar/R Arg may contribute to the stabilization of a closed pore state. Expression analyses show that AtNIP7;1 is selectively expressed in developing anther tissues of young floral buds of A. thaliana, principally in developing pollen grains of stage 9-11 anthers. Because boric acid is both an essential nutrient as well as a toxic compound at high concentrations, it is proposed that Tyr81 modulates transport and may provide an additional level of regulation for this transporter in male gametophyte development.

Dynamic programming procedure for searching optimal models to estimate substitution rates based on the maximum-likelihood method

The substitution rate in a gene can provide valuable information for understanding its functionality and evolution. A widely used method to estimate substitution rates is the maximum-likelihood method implemented in the CODEML program in the PAML package. A limited number of branch models, chosen based on a priori information or an interest in a particular lineage(s), are tested, whereas a large number of potential models are neglected. A complementary approach is also needed to test all or a large number of possible models to search for the globally optional model(s) of maximum likelihood. However, the computational time for this search even in a small number of sequences becomes impractically long. Thus, it is desirable to explore the most probable spaces to search for the optimal models. Using dynamic programming techniques, we developed a simple computational method for searching the most probable optimal branch-specific models in a practically feasible computational time. We propose three search methods to find the optimal models, which explored O(n) (method 1) to O(n(2)) (method 2 and method 3) models when the given phylogeny has n branches. In addition, we derived a formula to calculate the number of all possible models, revealing the complexity of finding the optimal branch-specific model. We show that in a reanalysis of over 50 previously published studies, the vast majority obtained better models with significantly higher likelihoods than the conventional hypothesis model methods.

Physiology and gene expression of the rice landrace Horkuch under salt stress

Good donors in breeding for salt tolerance are a prerequisite for food security under changing climatic conditions. Horkuch, a farmer-popular salt tolerant rice (Oryza sativa L.) variety from the south-west coast of Bangladesh was characterised up to maturity under NaCl stress, together with a modern variety (BRRI dhan41), a sensitive control (BRRI dhan29) and Pokkali, the salt-tolerant benchmark for rice. Horkuch had low reduction in shoot biomass, a low Na:K ratio in flag leaves, a low percent reduction in yield and good partitioning of Na in the older leaves, and maintained high levels of Ca and Mg in the flag leaves. In order to understand the physiology at the molecular level, the expression of salt-responsive genes was investigated using microarray analysis. Salt-stressed cDNA of Horkuch seedlings were hybridised with cDNA probes synthesised mainly from database sequences of Arabidopsis thaliana (L.) Heynh. The upregulated genes included transcription factors, signal transducers, metabolic enzymes, reactive oxygen species (ROS) scavengers, osmoprotectants and some specific salt-induced transcripts. An increase in expression of photosynthesis-related genes as well ROS scavengers suggested that this could be the reason for the better yield performance of Horkuch. The data therefore indicate Horkuch as a potential donor alternative to Pokkali in breeding programs for salt tolerance.

Functional divergence of the NIP III subgroup proteins involved altered selective constraints and positive selection

BACKGROUND

Nod26-like intrinsic proteins (NIPs) that belong to the aquaporin superfamily are unique to plants. According to homology modeling and phylogenetic analysis, the NIP subfamily can be further divided into three subgroups with distinct biological functions (NIP I, NIP II, and NIP III). In some grasses, the NIP III subgroup proteins (NIP2s) were demonstrated to be permeable to solutes with larger diameter, such as silicic acid and arsenous acids. However, to date there is no data-mining or direct experimental evidences for the permeability of such larger solutes for dicot NIP2s, although they exhibit similar three-dimensional structures as those in grasses. It is therefore intriguing to investigate the molecular mechanisms that drive the evolution of plant NIP2s.

RESULTS

The NIP III subgroup is more ancient with a divergence time that predates the monocot-dicot split. The proliferation of NIP2 genes in modern grass species is primarily attributed to whole genome and segmental chromosomal duplication events. The structure of NIP2 genes is relatively conserved, possessing five exons and four introns. All NIP2s possess an ar/R filter consisting of G, S, G, and R, except for the cucumber CsNIP2;2, where a small G in the H2 is substituted with the bulkier C residue. Our maximum likelihood analysis revealed that NIP2s, especially the loop A (LA) region, have undergone strong selective pressure for adaptive evolution. The analysis at the amino acid level provided strong statistical evidences for the functional divergence between monocot and dicot NIP III subgroup proteins. In addition, several SDPs (Specificity Determining Positions) responsible for functional specificity were predicted.

CONCLUSIONS

The present study provides the first evidences of functional divergence between dicot and monocot NIP2s, and suggests that positive selection, as well as a radical shift of evolutionary rate at some critical amino acid sites is the primary driver. These findings will expand our understanding to evolutionary mechanisms driving the functional diversification of monocot and dicot NIP III subgroup proteins.

Identification of the family of aquaporin genes and their expression in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Cotton (Gossypium spp.) is produced in over 30 countries and represents the most important natural fiber in the world. One of the primary factors affecting both the quantity and quality of cotton production is water. A major facilitator of water movement through cell membranes of cotton and other plants are the aquaporin proteins. Aquaporin proteins are present as diverse forms in plants, where they function as transport systems for water and other small molecules. The plant aquaporins belong to the large major intrinsic protein (MIP) family. In higher plants, they consist of five subfamilies including plasma membrane intrinsic proteins (PIP), tonoplast intrinsic proteins (TIP), NOD26-like intrinsic proteins (NIP), small basic intrinsic proteins (SIP), and the recently discovered X intrinsic proteins (XIP). Although a great deal is known about aquaporins in plants, very little is known in cotton.

RESULTS

From a molecular cloning effort, together with a bioinformatic homology search, 71 upland cotton (G. hirsutum) aquaporin genes were identified. The cotton aquaporins consist of 28 PIP and 23 TIP members with high sequence similarity. We also identified 12 NIP and 7 SIP members that showed more divergence. In addition, one XIP member was identified that formed a distinct 5th subfamily. To explore the physiological roles of these aquaporin genes in cotton, expression analyses were performed for a select set of aquaporin genes from each subfamily using semi-quantitative reverse transcription (RT)-PCR. Our results suggest that many cotton aquaporin genes have high sequence similarity and diverse roles as evidenced by analysis of sequences and their expression.

CONCLUSION

This study presents a comprehensive identification of 71 cotton aquaporin genes. Phylogenetic analysis of amino acid sequences divided the large and highly similar multi-gene family into the known 5 aquaporin subfamilies. Together with expression and bioinformatic analyses, our results support the idea that the genes identified in this study represent an important genetic resource providing potential targets to modify the water use properties of cotton.