PIECE: a database for plant gene structure comparison and evolution

Predicting the role of β-GAL genes in bean under abiotic stress and genome-wide characterization of β-GAL gene family members

Β-Gals are a subgroup of the glycoside hydrolase (GH) family of enzymes, which possess the Glyco_hydro_35 (GH35) domain. Although studies have been conducted on the β-Gal gene family in numerous plant species, no such research has been conducted on beans. The purpose of this study was to determine the gene expression levels of β-Gal genes in the leaf tissue of P. vulgaris under salt and drought stress using quantitative real-time polymerase chain reaction (qRT-PCR) and to perform a comprehensive analysis of β-Gal gene family members using bioinformatics tools. In the bean genome, 25 Pvul-βGAL proteins with amino acid numbers ranging from 291 to 1119, molecular weights from 32.94 to 126.56 kDa, and isoelectric points from 5.46 to 9.08 were identified. Both segmental and tandem duplication have occurred in β-Gal genes in the bean genome, and Pvul-BGAL genes have been subject to negative selection in the evolutionary process. For a deeper comprehension of the evolutionary proximity of Pvul-BGAL genes, a phylogenetic tree and synteny map were drawn together with Arabidopsis thaliana and Glycine max β-Gal genes. The expression profiles of β-Gal genes in different tissues of the bean were determined in silico. In addition, the expression profiles of β-Gal genes in the leaves of bean plants subjected to drought and salt stress were analyzed, and the role of β-Gal genes in salt and drought stress was estimated. In this study, the role of β-Gal gene family in abiotic stress response and the characterization of β-Gal genes in beans were determined for the first time and will provide a basis for future functional genomics studies.

Understanding saffron biology using omics- and bioinformatics tools: stepping towards a better Crocus phenome

Saffron is a unique plant in many aspects, and its cellular processes are regulated at multiple levels. The genetic makeup in the form of eight chromosome triplets (2n = 3x = 24) with a haploid genetic content (genome size) of 3.45 Gbp is decoded into different types of RNA by transcription. The RNA then translates into peptides and functional proteins, sometimes involving post-translational modifications too. The interactions of the genome, transcriptome, proteome and other regulatory molecules ultimately result in the complex set of primary and secondary metabolites of saffron metabolome. These complex interactions manifest in the form of a set of traits 'phenome' peculiar to saffron. The phenome responds to the environmental changes occurring in and around saffron and modify its response in respect of growth, development, disease response, stigma quality, apocarotenoid biosynthesis, and other processes. Understanding these complex relations between different yet interconnected biological activities is quite challenging in saffron where classical genetics has a very limited role owing to its sterility, and the absence of a whole-genome sequence. Omics-based technologies are immensely helpful in overcoming these limitations and developing a better understanding of saffron biology. In addition to creating a comprehensive picture of the molecular mechanisms involved in apocarotenoid synthesis, stigma biogenesis, corm activity, and flower development, omics-technologies will ultimately lead to the engineering of saffron plants with improved phenome.

Insights into the gene and protein structures of the CaSWEET family members in chickpea (Cicer arietinum), and their gene expression patterns in different organs under various stress and abscisic acid treatments

'Sugars Will Eventually be Exported Transporters' (SWEETs) are a group of sugar transporters that play crucial roles in various biological processes, particularly plant stress responses. However, no information is available yet for the CaSWEET family in chickpea. Here, we identified all putative CaSWEET members in chickpea, and obtained their major characteristics, including physicochemical patterns, chromosomal distribution, subcellular localization, gene organization, conserved motifs and three-dimensional protein structures. Subsequently, we explored available transcriptome data to compare spatiotemporal transcript abundance of CaSWEET genes in various major organs. Finally, we studied the changes in their transcript levels in leaves and/or roots following dehydration and exogenous abscisic acid treatments using RT-qPCR to obtain valuable information underlying their potential roles in chickpea responses to water-stress conditions. Our results provide the first insights into the characteristics of the CaSWEET family members and a foundation for further functional characterizations of selected candidate genes for genetic engineering of chickpea.

PID: An integrative and comprehensive platform of plant intron

Intron is a non-coding sequence of a broken gene and participates in important biological processes, such as transcription regulation, alternative splicing, and nuclear export. With the development of plant genomes, a comprehensive platform for intron analysis in plants must be established. Plant Intron Database (PID), a publicly available searchable database, was developed to efficiently store, query, analyze, and integrate intron resources in plants. The information of intron, exon, and gene can be searched by key words in PID. Users cannot only view intron length distribution pie chart and 5' and 3' splice site sequence feature maps in a statistical interface but can also browse intron information in a graphical visualization interface through JBrowse. ViroBlast for sequence homology searches, Intron detection and sequence interception tools were also provided. PID contains annotated genes from 118 sequenced plants, 24,782,048 introns, 30,843,049 exons, and 414 visual maps. This tool will greatly accelerate research on the distribution, length characteristics, and functions of introns in plants. PID is accessible at http://biodb.sdau.edu.cn/PID/index.php.

SilkDB 3.0: visualizing and exploring multiple levels of data for silkworm

SilkDB is an open-accessibility database and powerful platform that provides comprehensive information on the silkworm (Bombyx mori) genome. Since SilkDB 2.0 was released 10 years ago, vast quantities of data about multiple aspects of the silkworm have been generated, including genome, transcriptome, Hi-C and pangenome. To visualize data at these different biological levels, we present SilkDB 3.0 (https://silkdb.bioinfotoolkits.net), a visual analytic tool for exploring silkworm data through an interactive user interface. The database contains a high-quality chromosome-level assembly of the silkworm genome, and its coding sequences and gene sets are more accurate than those in the previous version. SilkDB 3.0 provides a view of the information for each gene at the levels of sequence, protein structure, gene family, orthology, synteny, genome organization and gives access to gene expression information, genetic variation and genome interaction map. A set of visualization tools are available to display the abundant information in the above datasets. With an improved interactive user interface for the integration of large data sets, the updated SilkDB 3.0 database will be a valuable resource for the silkworm and insect research community.

Natural variation in HsfA2 pre-mRNA splicing is associated with changes in thermotolerance during tomato domestication

Wild relatives of crops thrive in habitats where environmental conditions can be restrictive for productivity and survival of cultivated species. The genetic basis of this variability, particularly for tolerance to high temperatures, is not well understood. We examined the capacity of wild and cultivated accessions to acclimate to rapid temperature elevations that cause heat stress (HS). We investigated genotypic variation in thermotolerance of seedlings of wild and cultivated accessions. The contribution of polymorphisms associated with thermotolerance variation was examined regarding alterations in function of the identified gene. We show that tomato germplasm underwent a progressive loss of acclimation to strong temperature elevations. Sensitivity is associated with intronic polymorphisms in the HS transcription factor HsfA2 which affect the splicing efficiency of its pre-mRNA. Intron splicing in wild species results in increased synthesis of isoform HsfA2-II, implicated in the early stress response, at the expense of HsfA2-I which is involved in establishing short-term acclimation and thermotolerance. We propose that the selection for modern HsfA2 haplotypes reduced the ability of cultivated tomatoes to rapidly acclimate to temperature elevations, but enhanced their short-term acclimation capacity. Hence, we provide evidence that alternative splicing has a central role in the definition of plant fitness plasticity to stressful conditions.

Differential interaction of Or proteins with the PSY enzymes in saffron

Colored apocarotenoids accumulate at high concentrations in few plant species, where display a role in attraction of pollinators and seed dispersers. Among these apocarotenoids, crocins accumulate at high concentrations in the stigma of saffron and are responsible for the organoleptic and medicinal properties of this spice. Phytoene synthase and Orange protein are key for carotenoid biosynthesis and accumulation. We previously isolated four phytoene synthase genes from saffron with differential roles in carotenoid and apocarotenoid biosynthesis. However, the implications of Orange genes in the regulation of apocarotenoid accumulation are unknown. Here, we have identified two Orange genes from saffron, with different expression patterns. CsOr-a was mainly expressed in vegetative tissues and was induced by light and repressed by heat stress. Both CsOr-a and CsOr-b were expressed in stigmas but showed a different profile during the development of this tissue. The interactions of CsOr-a and CsOr-b were tested with all the four phytoene synthase proteins from saffron and with CsCCD2. None interactions were detected with CCD2 neither with the phytoene synthase 2, involved in apocarotenoid biosynthesis in saffron. The obtained results provide evidence of different mechanisms regulating the phytoene synthase enzymes in saffron by Orange for carotenoid and apocarotenoid accumulation in saffron.

Global Investigation of Cytochrome P450 Genes in the Chicken Genome

Cytochrome P450 (CYP) superfamily enzymes are broadly involved in a variety of physiological and toxicological processes. However, genome-wide analysis of this superfamily has never been investigated in the chicken genome. In this study, genome-wide analyses identified 45 chicken CYPs (cCYPs) from the chicken genome, and their classification and evolutionary relationships were investigated by phylogenetic, conserved protein motif, and gene structure analyses. The comprehensive evolutionary data revealed several remarkable characteristics of cCYPs, including the highly divergent and rapid evolution of the cCYPs, and the loss of cCYP2AF in the chicken genome. Furthermore, the cCYP expression profile was investigated by RNA-sequencing. The differential expression of cCYPs in developing embryos revealed the involvement of cCYPs in embryonic development. The significantly regulated cCYPs suggested its potential role in hepatic metabolism. Additionally, 11 cCYPs, including cCYP2AC1, cCYP2C23a, and cCYP2C23b, were identified as estrogen-responsive genes, which indicates that these cCYPs are involved in the estrogen-signaling pathway. Meanwhile, an expression profile analysis highlights the divergent role of different cCYPs. These data expand our view of the phylogeny and evolution of cCYPs, provide evolutionary insight, and can help elucidate the roles of cCYPs in physiological and toxicological processes in chicken.

Characterization of Cellulose Synthase A (CESA) Gene Family in Eudicots

Cellulose synthase A (CESA) is a key enzyme involved in the complex process of plant cell wall biosynthesis, and it remains a productive subject for research. We employed systems biology approaches to explore structural diversity of eudicot CESAs by exon-intron organization, mode of duplication, synteny, and splice site analyses. Using a combined phylogenetics and comparative genomics approach coupled with co-expression networks we reconciled the evolution of cellulose synthase gene family in eudicots and found that the basic forms of CESA proteins are retained in angiosperms. Duplications have played an important role in expansion of CESA gene family members in eudicots. Co-expression networks showed that primary and secondary cell wall modules are duplicated in eudicots. We also identified 230 simple sequence repeat markers in 103 eudicot CESAs. The 13 identified conserved motifs in eudicots will provide a basis for gene identification and functional characterization in other plants. Furthermore, we characterized (in silico) eudicot CESAs against senescence and found that expression levels of CESAs decreased during leaf senescence.

Genome-wide characterization and expression profiling of diacylglycerol acyltransferase genes from maize

Diacylglycerol acyltransferase (DGAT) catalyzes the only rate-limiting step in the pathway of plant oil (TAG) biosynthesis and is involved in plant development. In this study, five DGAT family members were identified from maize genome database. Phylogenetic analysis classified the ZmDGATs into type-I, II, and III clusters. Conserved functional domain analysis revealed that the proteins encoded by ZmDGAT1 contained conserved MBOAT domains, while two ZmDGAT2-encoding proteins harbored LPLAT domains. qRT-PCR analysis showed that ZmDGAT genes exhibited very high relative expression in developing seeds, especially at the early stage of seed development. Under various abiotic stress conditions, differential responses of ZmDGAT genes were observed. An overall significant induction of ZmDGAT genes under cold stress in leaves and a quick and strong response to osmotic stresses in roots were highlighted. This study provides useful information for understanding the roles of DGATs in oil accumulation and stress responses in higher plants.

Genome-wide analysis of the Glycerol-3-Phosphate Acyltransferase (GPAT) gene family reveals the evolution and diversification of plant GPATs

sn-Glycerol-3-phosphate 1-O-acyltransferase (GPAT) is an important enzyme that catalyzes the transfer of an acyl group from acyl-CoA or acyl-ACP to the sn-1 or sn-2 position of sn-glycerol-3-phosphate (G3P) to generate lysophosphatidic acids (LPAs). The functional studies of GPAT in plants demonstrated its importance in controlling storage and membrane lipid. Identifying genes encoding GPAT in a variety of plant species is crucial to understand their involvement in different metabolic pathways and physiological functions. Here, we performed genome-wide and evolutionary analyses of GPATs in plants. GPAT genes were identified in all algae and plants studied. The phylogenetic analysis showed that these genes group into three main clades. While clades I (GPAT9) and II (soluble GPAT) include GPATs from algae and plants, clade III (GPAT1-8) includes GPATs specific from plants that are involved in the biosynthesis of cutin or suberin. Gene organization and the expression pattern of GPATs in plants corroborate with clade formation in the phylogeny, suggesting that the evolutionary patterns is reflected in their functionality. Overall, our results provide important insights into the evolution of the plant GPATs and allowed us to explore the evolutionary mechanism underlying the functional diversification among these genes.

Characterization of OglDREB2A gene from African rice (Oryza glaberrima), comparative analysis and its transcriptional regulation under salinity stress

In this study, AP2 DNA-binding domain-containing transcription factor, OglDREB2A, was cloned from the African rice (Oryza glaberrima) and compared with 3000 rice genotypes. Further, the phylogenetic and various structural analysis was performed using in silico approaches. Further, to understand its allelic variation in rice, SNPs and indels were detected among the 3000 rice genotypes which indicated that while coding region is highly conserved, yet noncoding regions such as UTR and intron contained most of the variation. Phylogenetic analysis of the OglDREB2A sequence in different Oryza as well as in diverse eudicot species revealed that DREB from various Oryza species were diversed much earlier than other genes. Further, structural features and in silico analyses provided insights into different properties of OglDREB2A protein. The neutrality test on the coding region of OglDREB2A from different genotypes of O. glaberrima showed the lack of selection in this gene. Among the different developmental stages, it was upregulated at tillering and flag leaf under salinity treatment indicating its positive role in seedling and reproductive stage tolerance. Real-time PCR analysis also indicated the conserve expression pattern of this gene under salinity stress across the three different Oryza species having different degree of salinity tolerance.

PrLPAAT4, a Putative Lysophosphatidic Acid Acyltransferase from Paeonia rockii, Plays an Important Role in Seed Fatty Acid Biosynthesis

Lysophosphatidic acid acyltransferases (LPAATs) are essential for the acylation of lysophosphatidic acid (LPA) and the synthesis of phosphatidic acid (PA), a key intermediate in the synthesis of membrane phospholipids and storage lipids. Here, a putative lysophosphatidic acid acyltransferase gene, designated PrLPAAT4, was isolated from seed unsaturated fatty acid (UFA)-rich P. rockii. The complete PrLPAAT4 cDNA contained a 1116-bp open reading frame (ORF), encoding a 42.9 kDa protein with 371 amino acid residues. Bioinformatic analysis indicates that PrLPAAT4 is a plasma membrane protein belonging to acyl-CoA:1-acylglycerol-sn-3-phosphate acyltranferases (AGPAT) family. PrLPAAT4 shared high sequence similarity with its homologs from Citrus clementina, Populus trichocarpa, Manihot esculenta, and Ricinus communis. In Arabidopsis, overexpression of PrLPAAT4 resulted in a significant increase in the content of oleic acid (OA) and total fatty acids (FAs) in seeds. AtDGAT1, AtGPAT9, and AtOleosin, involved in TAG assembly, were upregulated in PrLPAAT4-overexpressing lines. These results indicated that PrLPAAT4 functions may be as a positive regulator in seed FA biosynthesis.

Functional divergence and origin of the DAG-like gene family in plants

The nuclear-encoded DAG-like (DAL) gene family plays critical roles in organelle C-to-U RNA editing in Arabidopsis thaliana. However, the origin, diversification and functional divergence of DAL genes remain unclear. Here, we analyzed the genomes of diverse plant species and found that: DAL genes are specific to spermatophytes, all DAL genes share a conserved gene structure and protein similarity with the inhibitor I9 domain of subtilisin genes found in ferns and mosses, suggesting that DAL genes likely arose from I9-containing proproteases via exon shuffling. Based on phylogenetic inference, DAL genes can be divided into five subfamilies, each composed of putatively orthologous and paralogous genes from different species, suggesting that all DAL genes originated from a common ancestor in early seed plants. Significant type I functional divergence was observed in 6 of 10 pairwise comparisons, indicating that shifting functional constraints have contributed to the evolution of DAL genes. This inference is supported by the finding that functionally divergent amino acids between subfamilies are predominantly located in the DAL domain, a critical part of the RNA editosome. Overall, these findings shed light on the origin of DAL genes in spermatophytes and outline functionally important residues involved in the complexity of the RNA editosome.

Characterization of the maize lipoxygenase gene family in relation to aflatoxin accumulation resistance

Maize (Zea mays L.) is a globally important staple food crop prone to contamination by aflatoxin, a carcinogenic secondary metabolite produced by the fungus Aspergillus flavus. An efficient approach to reduce accumulation of aflatoxin is the development of germplasm resistant to colonization and toxin production by A. flavus. Lipoxygenases (LOXs) are a group of non-heme iron containing dioxygenase enzymes that catalyze oxygenation of polyunsaturated fatty acids (PUFAs). LOX derived oxylipins play critical roles in plant defense against pathogens including A. flavus. The objectives of this study were to summarize sequence diversity and expression patterns for all LOX genes in the maize genome, and map their effect on aflatoxin accumulation via linkage and association mapping. In total, 13 LOX genes were identified, characterized, and mapped. The sequence of one gene, ZmLOX10, is reported from 5 inbred lines. Genes ZmLOX1/2, 5, 8, 9, 10 and 12 (GRMZM2G156861, or V4 numbers ZM00001D042541 and Zm00001D042540, GRMZM2G102760, GRMZM2G104843, GRMZM2G017616, GRMZM2G015419, and GRMZM2G106748, respectively) fell under previously published QTL in one or more mapping populations and are linked to a measurable reduction of aflatoxin in maize grains. Association mapping results found 28 of the 726 SNPs tested were associated with reduced aflatoxin levels at p ≤ 9.71 x 10-4 according to association statistics. These fell within or near nine of the ZmLOX genes. This work confirms the importance of some lipoxygenases for resistance to aflatoxin accumulation and may be used to direct future genetic selection in maize.

Evolutionary expansion and structural functionalism of the ancient family of profilin proteins

Structure and functional similarities of a recent protein's orthologs with its ancient counterpart are largely determined by the configuration of evolutionary preservation of amino acids. The emergence of genome sequencing databases allowed dissecting the evolutionarily important gene families at a comprehensive and genome-wide scale. The profilin multi-gene family is an ancient, universal, and functionally diverged across kingdoms, which regulates various aspects of cellular development in both prokarya and eukarya, especially cell-wall maintenance through actin sequestering, nucleation and cytokinesis. We performed a meta-analysis of the evolutionary expansion, structural conservation, evolution of function motifs, and transcriptional biases of profilin proteins across kingdoms. An exhaustive search of various genome databases of cyanobacteria, fungi, animalia and plantae kingdoms revealed 172 paralogous/orthologous profilins those were phylogenetically clustered in various groups. Orthologous gene comparisons indicated that segmental and tandem duplication events under strong purifying selection are predominantly responsible for their convoluted structural divergences. Evidently, structural divergences were more prevalent in the paralogs than orthologs, and evolutionary variations in the exon/intron architecture were accomplished by 'exon/intron-gain' and insertion/deletion during sequence-exonization. Remarkably, temporal expression evolution of profilin paralogs/homeologs during cotton fiber domestication provides evolutionary impressions of the selection of highly diverged transcript abundance notably in the fiber morpho-evolution. These results provide global insights into the profilin evolution, their structural design across taxa; and their future utilization in translational research.

Cloning, Characterization, and Expression Analysis of a Gene Encoding a Putative Lysophosphatidic Acid Acyltransferase from Seeds of Paeonia rockii

Tree peony (Paeonia section Moutan DC.) is an excellent woody oil crop, and the cloning and functional analysis of genes related to fatty acid (FA) metabolism from this organism has not been reported. Lysophosphatidic acid acyltransferase (LPAAT), which converts lysophosphatidic acid (LPA) to phosphatidic acid (PA), catalyzes the addition of fatty acyl moieties to the sn-2 position of the LPA glycerol backbone in triacylglycerol (TAG) biosynthesis. This project reports a putative lysophosphatidic acid acyltransferase gene PrLPAAT1 isolated from Paeonia rockii. Our data indicated that PrLPAAT1 has 1047 nucleotides and encodes a putative 38.8 kDa protein with 348 amino acid residues. Bioinformatic analysis demonstrated that PrLPAAT1 contains two transmembrane domains (TMDs). Subcellular localization analysis confirmed that PrLPAAT1 is a plasma membrane protein. Phylogenetic analysis revealed that PrLPAAT1 shared 74.3 and 65.5% amino acid sequence identities with the LPAAT1 sequences from columbine and grape, respectively. PrLPAAT1 belongs to AGPAT family, and may have acyltransferase activity. PrLPAAT1 was ubiquitously expressed in diverse tissues, and PrLPAAT1 expression was higher in the flower and developing seed. PrLPAAT1 is probably an important component in the FA accumulation process, especially during the early stages of seed development. PrLPAAT1 overexpression using a seed-specific promoter increased total FA content and the main FA accumulation in Arabidopsis transgenic plants.

PIECE 2.0: an update for the plant gene structure comparison and evolution database

PIECE (Plant Intron Exon Comparison and Evolution) is a web-accessible database that houses intron and exon information of plant genes. PIECE serves as a resource for biologists interested in comparing intron–exon organization and provides valuable insights into the evolution of gene structure in plant genomes. Recently, we updated PIECE to a new version, PIECE 2.0 (http://probes.pw.usda.gov/piece or http://aegilops.wheat.ucdavis.edu/piece). PIECE 2.0 contains annotated genes from 49 sequenced plant species as compared to 25 species in the previous version. In the current version, we also added several new features: (i) a new viewer was developed to show phylogenetic trees displayed along with the structure of individual genes; (ii) genes in the phylogenetic tree can now be also grouped according to KOG (The annotation of Eukaryotic Orthologous Groups) and KO (KEGG Orthology) in addition to Pfam domains; (iii) information on intronless genes are now included in the database; (iv) a statistical summary of global gene structure information for each species and its comparison with other species was added; and (v) an improved GSDraw tool was implemented in the web server to enhance the analysis and display of gene structure. The updated PIECE 2.0 database will be a valuable resource for the plant research community for the study of gene structure and evolution.

Computational Tools for Genomic Studies in Plants

In recent years, the genomic sequence of numerous plant species including the main crop species has been determined. Computational tools have been developed to deal with the issue of which plant has been sequenced and where is the sequence hosted. In this mini-review, the databases for genome projects, the databases created to host species/clade projects and the databases developed to perform plant comparative genomics are revised. Because of their importance in modern research, an in-depth analysis of the plant comparative genomics databases has been performed. This comparative analysis is focused in the common and specific computational tools developed to achieve the particular objectives of each database. Besides, emerging high-performance bioinformatics tools specific for plant research are commented. What kind of computational approaches should be implemented in next years to efficiently analyze plant genomes is discussed.

Diversity and evolution of plant diacylglycerol acyltransferase (DGATs) unveiled by phylogenetic, gene structure and expression analyses

Since the first diacylglycerol acyltransferase (DGAT) gene was characterized in plants, a number of studies have focused on understanding the role of DGAT activity in plant triacylglycerol (TAG) biosynthesis. DGAT enzyme is essential in controlling TAGs synthesis and is encoded by different genes. DGAT1 and DGAT2 are the two major types of DGATs and have been well characterized in many plants. On the other hand, the DGAT3 and WS/DGAT have received less attention. In this study, we present the first general view of the presence of putative DGAT3 and WS/DGAT in several plant species and report on the diversity and evolution of these genes and its relationships with the two main DGAT genes (DGAT1 and DGAT2). According to our analyses DGAT1, DGAT2, DGAT3 and WS/DGAT are very divergent genes and may have distinct origin in plants. They also present divergent expression patterns in different organs and tissues. The maintenance of several types of genes encoding DGAT enzymes in plants demonstrates the importance of DGAT activity for TAG biosynthesis. Evolutionary history studies of DGATs coupled with their expression patterns help us to decipher their functional role in plants, helping to drive future biotechnological studies.

Diversity and evolution of plant diacylglycerol acyltransferase (DGATs) unveiled by phylogenetic, gene structure and expression analyses

Since the first diacylglycerol acyltransferase (DGAT) gene was characterized in plants, a number of studies have focused on understanding the role of DGAT activity in plant triacylglycerol (TAG) biosynthesis. DGAT enzyme is essential in controlling TAGs synthesis and is encoded by different genes. DGAT1 and DGAT2 are the two major types of DGATs and have been well characterized in many plants. On the other hand, the DGAT3 and WS/DGAT have received less attention. In this study, we present the first general view of the presence of putative DGAT3 and WS/DGAT in several plant species and report on the diversity and evolution of these genes and its relationships with the two main DGAT genes (DGAT1 and DGAT2). According to our analyses DGAT1, DGAT2, DGAT3 and WS/DGAT are very divergent genes and may have distinct origin in plants. They also present divergent expression patterns in different organs and tissues. The maintenance of several types of genes encoding DGAT enzymes in plants demonstrates the importance of DGAT activity for TAG biosynthesis. Evolutionary history studies of DGATs coupled with their expression patterns help us to decipher their functional role in plants, helping to drive future biotechnological studies.

Evolutionary and Expression Analysis Provides Evidence for the Plant Glutamate-like Receptors Family is Involved in Woody Growth-related Function

Glutamate-like receptors (GLRs) is a highly conserved family of ligand-gated ion channels, which have been associated with various physiological and developmental processes. Here, we investigated the evolutionary pattern of GLRs in plants. We observed that tandem duplications occupied the largest proportion of the plant GLR gene family expansion. Based on a phylogenetic tree, we suggested a new subfamily, GLR4, which is widespread in angiosperm but absence on Brassicales. Meanwhile, because GLR1 and GLR2 subfamilies were potential sister clades, we combined them into the GLR1&2 subfamily. A comparative analysis of plant GLR subfamilies revealed that selective forces shaped the GLR1&2 repertoires in the stems of eudicotyledons with distinct functional preferences. Moreover, GLR1&2 formed a species-specific highwoody-expanded subfamily, with preferential expression in the cambial-enriched and shoot apical meristem fractions of the highwood species. Together, these findings lay the foundation for evolutionary analysis of plant GLRs over the entire plant timescale and identified unique targets for manipulating the woody-growth behaviours of plant GLRs.

Intron gain by tandem genomic duplication: a novel case in a potato gene encoding RNA-dependent RNA polymerase

The origin and subsequent accumulation of spliceosomal introns are prominent events in the evolution of eukaryotic gene structure. However, the mechanisms underlying intron gain remain unclear because there are few proven cases of recently gained introns. In an RNA-dependent RNA polymerase (RdRp) gene, we found that a tandem duplication occurred after the divergence of potato and its wild relatives among other Solanum plants. The duplicated sequence crosses the intron-exon boundary of the first intron and the second exon. A new intron was detected at this duplicated region, and it includes a small previously exonic segment of the upstream copy of the duplicated sequence and the intronic segment of the downstream copy of the duplicated sequence. The donor site of this new intron was directly obtained from the small previously exonic segment. Most of the splicing signals were inherited directly from the parental intron/exon structure, including a putative branch site, the polypyrimidine tract, the 3' splicing site, two putative exonic splicing enhancers, and the GC contents differed between the intron and exon. In the widely cited model of intron gain by tandem genomic duplication, the duplication of an AGGT-containing exonic segment provides the GT and AG splicing sites for the new intron. Our results illustrate that the tandem duplication model of intron gain should be diverse in terms of obtaining the proper splicing signals.

Transcript structure and domain display: a customizable transcript visualization tool

Transcript Structure and Domain Display (TSDD) is a publicly available, web-based program that provides publication quality images of transcript structures and domains. TSDD is capable of producing transcript structures from GFF/GFF3 and BED files. Alternatively, the GFF files of several model organisms have been pre-loaded so that users only needs to enter the locus IDs of the transcripts to be displayed. Visualization of transcripts provides many benefits to researchers, ranging from evolutionary analysis of DNA-binding domains to predictive function modeling.

AVAILABILITY AND IMPLEMENTATION

TSDD is freely available for non-commercial users at http://shenlab.sols.unlv.edu/shenlab/software/TSD/transcript_display.html

CONTACT

: jeffery.shen@unlv.nevada.edu.

Novel Structural and Functional Motifs in cellulose synthase (CesA) Genes of Bread Wheat (Triticum aestivum, L.)

Cellulose is the primary determinant of mechanical strength in plant tissues. Late-season lodging is inversely related to the amount of cellulose in a unit length of the stem. Wheat is the most widely grown of all the crops globally, yet information on its CesA gene family is limited. We have identified 22 CesA genes from bread wheat, which include homoeologs from each of the three genomes, and named them as TaCesAXA, TaCesAXB or TaCesAXD, where X denotes the gene number and the last suffix stands for the respective genome. Sequence analyses of the CESA proteins from wheat and their orthologs from barley, maize, rice, and several dicot species (Arabidopsis, beet, cotton, poplar, potato, rose gum and soybean) revealed motifs unique to monocots (Poales) or dicots. Novel structural motifs CQIC and SVICEXWFA were identified, which distinguished the CESAs involved in the formation of primary and secondary cell wall (PCW and SCW) in all the species. We also identified several new motifs specific to monocots or dicots. The conserved motifs identified in this study possibly play functional roles specific to PCW or SCW formation. The new insights from this study advance our knowledge about the structure, function and evolution of the CesA family in plants in general and wheat in particular. This information will be useful in improving culm strength to reduce lodging or alter wall composition to improve biofuel production.

JuncDB: an exon-exon junction database

Intron positions upon the mRNA transcript are sometimes remarkably conserved even across distantly related eukaryotic species. This has made the comparison of intron-exon architectures across orthologous transcripts a very useful tool for studying various evolutionary processes. Moreover, the wide range of functions associated with introns may confer biological meaning to evolutionary changes in gene architectures. Yet, there is currently no database that offers such comparative information. Here, we present JuncDB (http://juncdb.carmelab.huji.ac.il/), an exon-exon junction database dedicated to the comparison of architectures between orthologous transcripts. It covers nearly 40,000 sets of orthologous transcripts spanning 88 eukaryotic species. JuncDB offers a user-friendly interface, access to detailed information, instructive graphical displays of the comparative data and easy ways to download data to a local computer. In addition, JuncDB allows the analysis to be carried out either on specific genes, or at a genome-wide level for any selected group of species.

Characterization of the Maize Chitinase Genes and Their Effect on Aspergillus flavus and Aflatoxin Accumulation Resistance

Maize (Zea mays L.) is a crop of global importance, but prone to contamination by aflatoxins produced by fungi in the genus Aspergillus. The development of resistant germplasm and the identification of genes contributing to resistance would aid in the reduction of the problem with a minimal need for intervention by farmers. Chitinolytic enzymes respond to attack by potential pathogens and have been demonstrated to increase insect and fungal resistance in plants. Here, all chitinase genes in the maize genome were characterized via sequence diversity and expression patterns. Recent evolution within this gene family was noted. Markers from within each gene were developed and used to map the phenotypic effect on resistance of each gene in up to four QTL mapping populations and one association panel. Seven chitinase genes were identified that had alleles associated with increased resistance to aflatoxin accumulation and A. flavus infection in field grown maize. The chitinase in bin 1.05 identified a new and highly significant QTL, while chitinase genes in bins 2.04 and 5.03 fell directly beneath the peaks of previously published QTL. The expression patterns of these genes corroborate possible grain resistance mechanisms. Markers from within the gene sequences or very closely linked to them are presented to aid in the use of marker assisted selection to improve this trait.

Genome-wide identification and comparative analysis of the heat shock transcription factor family in Chinese white pear (Pyrus bretschneideri) and five other Rosaceae species

BACKGROUND

Heat shock transcription factors (Hsfs), which act as important transcriptional regulatory proteins in eukaryotes, play a central role in controlling the expression of heat-responsive genes. At present, the genomes of Chinese white pear ('Dangshansuli') and five other Rosaceae fruit crops have been fully sequenced. However, information about the Hsfs gene family in these Rosaceae species is limited, and the evolutionary history of the Hsfs gene family also remains unresolved.

RESULTS

In this study, 137 Hsf genes were identified from six Rosaceae species (Pyrus bretschneideri, Malus × domestica, Prunus persica, Fragaria vesca, Prunus mume, and Pyrus communis), 29 of which came from Chinese white pear, designated as PbHsf. Based on the structural characteristics and phylogenetic analysis of these sequences, the Hsf family genes could be classified into three main groups (classes A, B, and C). Segmental and dispersed duplications were the primary forces underlying Hsf gene family expansion in the Rosaceae. Most of the PbHsf duplicated gene pairs were dated back to the recent whole-genome duplication (WGD, 30-45 million years ago (MYA)). Purifying selection also played a critical role in the evolution of Hsf genes. Transcriptome data demonstrated that the expression levels of the PbHsf genes were widely different. Six PbHsf genes were upregulated in fruit under naturally increased temperature.

CONCLUSION

A comprehensive analysis of Hsf genes was performed in six Rosaceae species, and 137 full length Hsf genes were identified. The results presented here will undoubtedly be useful for better understanding the complexity of the Hsf gene family and will facilitate functional characterization in future studies.

GSDS 2.0: an upgraded gene feature visualization server

Summary: Visualizing genes’ structure and annotated features helps biologists to investigate their function and evolution intuitively. The Gene Structure Display Server (GSDS) has been widely used by more than 60 000 users since its first publication in 2007. Here, we reported the upgraded GSDS 2.0 with a newly designed interface, supports for more types of annotation features and formats, as well as an integrated visual editor for editing the generated figure. Moreover, a user-specified phylogenetic tree can be added to facilitate further evolutionary analysis. The full source code is also available for downloading.

Availability and implementation: Web server and source code are freely available at http://gsds.cbi.pku.edu.cn.

Contact: gaog@mail.cbi.pku.edu.cn or gsds@mail.cbi.pku.edu.cn

Supplementary information: Supplementary data are available at Bioinformatics online.

Characterization of the lipoxygenase (LOX) gene family in the Chinese white pear (Pyrus bretschneideri) and comparison with other members of the Rosaceae

Background

Lipoxygenases (LOXs), a type of non-haem iron-containing dioxygenase, are ubiquitous enzymes in plants and participate in the formation of fruit aroma which is a very important aspect of fruit quality. Amongst the various aroma volatiles, saturated and unsaturated alcohols and aldehydes provide the characteristic aroma of the fruit. These compounds are formed from unsaturated fatty acids through oxidation, pyrolysis and reduction steps. This biosynthetic pathway involves at least four enzymes, including LOX, the enzyme responsible for lipid oxidation. Although some studies have been conducted on the LOX gene family in several species including Arabidopsis, soybean, cucumber and apple, there is no information from pear; and the evolutionary history of this gene family in the Rosaceae is still not resolved.

Results

In this study we identified 107 LOX homologous genes from five Rosaceous species (Pyrus bretschneideri, Malus × domestica, Fragaria vesca, Prunus mume and Prunus persica); 23 of these sequences were from pear. By using structure analysis, phylogenic analysis and collinearity analysis, we identified variation in gene structure and revealed the phylogenetic evolutionary relationship of this gene family. Expression of certain pear LOX genes during fruit development was verified by analysis of transcriptome data.

Conclusions

23 LOX genes were identified in pear and these genes were found to have undergone a duplication 30–45 MYA; most of these 23 genes are functional. Specific gene duplication was found on chromosome4 in the pear genome. Useful information was provided for future research on the evolutionary history and transgenic research on LOX genes.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-444) contains supplementary material, which is available to authorized users.

Uncovering legumain genes in rice

Legumains are Asn specific cysteine proteases physiologically related to the biosynthesis of vacuolar components, degradation of storage proteins and programmed cell death. The present work identifies and characterizes the genic family of legumains in rice (Oryza sativa), which comprises five different loci. Rice legumains (OsaLegs) were ubiquitously detected in all plant tissues analyzed. However, phylogenetic analyses and gene expression studies demonstrated greater association of OsaLeg2 and OsaLeg3 to seed-related legumains, whereas OsaLeg1, 4 and 5 would act as vegetative-related proteases. Additionally, OsaLeg1 mRNA is strongly induced in senescent leaves. All rice legumain genes respond in different ways to environmental conditions such as wounding, salt and abscisic acid treatments. Mainly, wounding is capable of inducing all the four expressed genes OsaLeg1, 2, 3 and 4. Alternative splicing isoforms, with potential to generate pre-activated OsaLeg1 and OsaLeg2 nonvacuolar enzymes under different environmental situations were also observed.

Presence-absence variation in A. thaliana is primarily associated with genomic signatures consistent with relaxed selective constraints

The sequencing of multiple genomes of the same plant species has revealed polymorphic gene and exon loss. Genes associated with disease resistance are overrepresented among those showing structural variations, suggesting an adaptive role for gene and exon presence–absence variation (PAV). To shed light on the possible functional relevance of polymorphic coding region loss and the mechanisms driving this process, we characterized genes that have lost entire exons or their whole coding regions in 17 fully sequenced Arabidopsis thaliana accessions. We found that although a significant enrichment in genes associated with certain functional categories is observed, PAV events are largely restricted to genes with signatures of reduced essentiality: PAV genes tend to be newer additions to the genome, tissue specific, and lowly expressed. In addition, PAV genes are located in regions of lower gene density and higher transposable element density. Partial coding region PAV events were associated with only a marginal reduction in gene expression level in the affected accession and occurred in genes with higher levels of alternative splicing in the Col-0 accession. Together, these results suggest that although adaptive scenarios cannot be ruled out, PAV events can be explained without invoking them.