MicroSyn: a user friendly tool for detection of microsynteny in a gene family

Analysis of Carica papaya Informs Lineage-Specific Evolution of the Aquaporin (AQP) Family in Brassicales

Aquaporins (AQPs), a type of intrinsic membrane proteins that transport water and small solutes across biological membranes, play crucial roles in plant growth and development. This study presents a first genome-wide identification and comparative analysis of the AQP gene family in papaya (Carica papaya L.), an economically and nutritionally important fruit tree of tropical and subtropical regions. A total of 29 CpAQP genes were identified, which represent five subfamilies, i.e., nine plasma intrinsic membrane proteins (PIPs), eight tonoplast intrinsic proteins (TIPs), seven NOD26-like intrinsic proteins (NIPs), two X intrinsic proteins (XIPs), and three small basic intrinsic proteins (SIPs). Although the family is smaller than the 35 members reported in Arabidopsis, it is highly diverse, and the presence of CpXIP genes as well as orthologs in Moringa oleifera and Bretschneidera sinensis implies that the complete loss of the XIP subfamily in Arabidopsis is lineage-specific, sometime after its split with papaya but before Brassicaceae-Cleomaceae divergence. Reciprocal best hit-based sequence comparison of 530 AQPs and synteny analyses revealed that CpAQP genes belong to 29 out of 61 identified orthogroups, and lineage-specific evolution was frequently observed in Brassicales. Significantly, the well-characterized NIP3 group was completely lost; lineage-specific loss of the NIP8 group in Brassicaceae occurred sometime before the divergence with Cleomaceae, and lineage-specific loss of NIP2 and SIP3 groups in Brassicaceae occurred sometime after the split with Cleomaceae. In contrast to a predominant role of recent whole-genome duplications (WGDs) on the family expansion in B. sinensis, Tarenaya hassleriana, and Brassicaceae plants, no recent AQP repeats were identified in papaya, and ancient WGD repeats are mainly confined to the PIP subfamily. Subfamily even group-specific evolution was uncovered via comparing exon-intron structures, conserved motifs, the aromatic/arginine selectivity filter, and gene expression profiles. Moreover, down-regulation during fruit ripening and expression divergence of duplicated CpAQP genes were frequently observed in papaya. These findings will not only improve our knowledge on lineage-specific family evolution in Brassicales, but also provide valuable information for further studies of AQP genes in papaya and species beyond.

Genome Context Viewer (GCV) version 2: enhanced visual exploration of multiple annotated genomes

The Genome Context Viewer is a web application for identifying, aligning, and visualizing genomic regions based on their micro and macrosyntenic structures. By using functional elements such as gene annotations as the unit of search and comparison, the Genome Context Viewer can compute and display relationships between regions across many assemblies from federated data sources in real-time, enabling users to rapidly explore multiple annotated genomes and identify divergence and structural events that can help provide insight into evolutionary mechanisms associated with functional consequences. In this work, we introduce version 2 of the Genome Context Viewer and highlight new features that enhance usability, performance, and ease of deployment.

Comprehensive Genomic Survey, Structural Classification, and Expression Analysis of WRKY Transcription Factor Family in Rhododendron simsii

(1) Rhododendron is one of the top ten traditional flowers in China, with both high ornamental and economic values. However, with the change of the environment, Rhododendron suffers from various biological stresses. The WRKY transcription factor is a member of the most crucial transcription factor families, which plays an essential regulatory role in a variety of physiological processes and developmental stresses. (2) In this study, 57 RsWRKYs were identified using genome data and found to be randomly distributed on 13 chromosomes. Based on gene structure and phylogenetic relationships, 57 proteins were divided into three groups: I, II, and III. Multiple alignments of RsWRKYs with Arabidopsis thaliana homologous genes revealed that WRKY domains in different groups had different conserved sites. RsWRKYs have a highly conserved domain, WRKYGQK, with three variants, WRKYGKK, WRKYGEK, and WRKYGRK. Furthermore, cis-acting elements analysis revealed that all of the RsWRKYs had stress and plant hormone cis-elements, with figures varying by group. Finally, the expression patterns of nine WRKY genes treated with gibberellin acid (GA), methyl jasmonate (MeJA), heat, and drought in Rhododendron were also measured using quantitative real-time PCR (qRT-PCR). The results showed that the expression levels of the majority of RsWRKY genes changed in response to multiple phytohormones and abiotic stressors. (3) This current study establishes a theoretical basis for future studies on the response of RsWRKY transcription factors to various hormone and abiotic stresses as well as a significant foundation for the breeding of new stress-tolerant Rhododendron varieties.

Evolution and function of ubiquitin-specific proteases (UBPs): Insight into seed development roles in tung tree (Vernicia fordii)

The tung oil produced by the tung tree (Vernicia fordii) provides resources for the manufacture of biodiesel. Ubiquitin-specific proteases (UBPs) are the largest group of deubiquitinases and play key roles in regulating development and stress responses. Here, 21 UBPs were identified in V. fordii, roughly one-half the number found in Manihot esculenta and Hevea brasiliensis. Most UBP duplications are produced from whole-genome duplication (WGD), and significant differences in gene retention existed among Euphorbiaceae. The great majority of UBP-containing blocks in V. fordii, V. montana, Ricinus communis, and Jatropha curcas exhibited extensive conservation with the duplicated regions of M. esculenta and H. brasiliensis. These blocks formed 14 orthologous groups, indicating they shared WGD with UBPs in M. esculenta and H. brasiliensis, but most of these UBPs copies were lost. The UBP orthologs contained significant functional divergence which explained the susceptibility of V. fordii to Fusarium wilt and the resistance of V. montana to Fusarium wilt. The expression patterns and experiments suggested that Vf03G1417 could affect the seed-related traits and positively regulate the seed oil accumulation. This study provided important insights into the evolution of UBPs in Euphorbiaceae and identified important candidate VfUBPs for marker-assisted breeding in V. fordii.

SWEET Transporters and the Potential Functions of These Sequences in Tea (Camellia sinensis)

Tea (Camellia sinensis) is an important economic beverage crop. Its flowers and leaves could be used as healthcare tea for its medicinal value. SWEET proteins were recently identified in plants as sugar transporters, which participate in diverse physiological processes, including pathogen nutrition, seed filling, nectar secretion, and phloem loading. Although SWEET genes have been characterized and identified in model plants, such as Arabidopsis thaliana and Oryza sativa, there is very little knowledge of these genes in C. sinensis. In this study, 28 CsSWEETs were identified in C. sinensis and further phylogenetically divided into four subfamilies with A. thaliana. These identified CsSWEETs contained seven transmembrane helixes (TMHs) which were generated by an ancestral three-TMH unit with an internal duplication experience. Microsynteny analysis revealed that the large-scale duplication events were the main driving forces for members from CsSWEET family expansion in C. sinensis. The expression profiles of the 28 CsSWEETs revealed that some genes were highly expressed in reproductive tissues. Among them, CsSWEET1a might play crucial roles in the efflux of sucrose, and CsSWEET17b could control fructose content as a hexose transporter in C. sinensis. Remarkably, CsSWEET12 and CsSWEET17c were specifically expressed in flowers, indicating that these two genes might be involved in sugar transport during flower development. The expression patterns of all CsSWEETs were differentially regulated under cold and drought treatments. This work provided a systematic understanding of the members from the CsSWEET gene family, which would be helpful for further functional studies of CsSWEETs in C. sinensis.

Systematic analysis of the OFP genes in six Rosaceae genomes and their roles in stress response in Chinese pear (Pyrus bretschneideri)

OVATE family proteins (OFPs) are the plant-specific transcription factors, and have significant functions in regulating plant growth, development and resistance. The OFP genes have been investigated in several plants, but they still lack a systematic analysis of OFP genes in Chinese pear and some other five Rosaceae genomes. Here, 28 PbrOFPs were identified within Chinese pear and compared them with those of other five Rosaceae genomes. Evolutionary tree revealed that all OFP genes from six Rosaceae genomes were divided into eight groups. Seventeen conserved microsynteny regions were detected in Chinese pear genome, suggested that these PbrOFP genes might be considered to have originated from the large-scale duplication events., indicating these PbrOFP genes might contain specialized regulatory mechanisms in these tissues, such as flower, ovary and fruit. Remarkably, two PbrOFP genes (Pbr010426.1 and Pbr010427.1) were up-regulated under Venturia nashicola treatment, and five PbrOFP genes were up-regulated under PEG treatment, suggesting that these genes might play crucial roles in defence to environmental stresses. Our data presented a systematic analysis and might aid in the selection of appropriate PbrOFPs for further functional studies in Chinese pear, especially in response to the mechanism of biotic and abiotic stresses.

Genes encoding light-harvesting chlorophyll a/b-binding proteins in papaya (Carica papaya L.) and insight into lineage-specific evolution in Brassicaceae

Light-harvesting chlorophyll a/b-binding (Lhc) proteins comprise a plant-specific superfamily involved in photosynthesis and stress responses. Despite their importance, little is known in papaya (Carica papaya), an economically important tree fruit crop as well as a species close to the model plant arabidopsis (Arabidopsis thaliana). This study reports a first genome-wide analysis of Lhc superfamily genes in papaya, and a total of 28 members that represent four defined families or 26 orthologous groups were identified from the papaya genome. The superfamily number is comparable to 28 or 27 reported in castor (Ricinus communis) and jatropha (Jatropha curcas), respectively, two Euphorbiaceous plants also without any recent whole-genome duplication (WGD), but relatively less than 35, 34, 32, 32, 37, 30 or 32 present in cassava (Manihot esculenta), arabidopsis, A. lyrata, A. halleri, Capsella rubella, C. grandiflora, and Eutrema salsugineum, respectively, representative species having experienced one or two recent WGDs. Local duplication was shown to play a predominant role in gene expansion in papaya, castor, and jatropha, which is only confined to the Lhcb1 group. By contrast, WGD plays a relatively more important role in cassava, arabidopsis, and other Brassicaceous plants. Further comparison of Brassicaceous plants revealed that loss of the SEP6 group in arabidopsis is lineage-specific, occurring sometime after papaya-arabidopsis divergence but before the radiation of Brassicaceous plants. Transcriptional profiling revealed a leaf-preferential expression pattern of most CpLhc superfamily genes and their transcript levels were markedly regulated by three abiotic stresses, i.e., mimicking drought, cold, and high salt. These findings not only facilitate further functional studies in papaya, but also improve our knowledge on lineage-specific evolution of this special gene superfamily in Brassicaceae.

Genome-wide analysis of lectin receptor-like kinases family from potato (Solanum tuberosum L.)

Lectin receptor-like kinases (LecRLKs) are involved in responses to diverse environmental stresses and pathogenic microbes. A comprehensive acknowledgment of the family members in potato (Solanum tuberosum) genome is largely limited until now. In total, 113 potato LecRLKs (StLecRLKs) were first identified, including 85 G-type, 26 L-type and 2 C-type members. Based on phylogenetic analysis, StLecRLKs were sub-grouped into seven clades, including C-type, L-type, G-I, G-II, G-III G-IV and G-V. Chromosomal distribution and gene duplication analysis revealed the expansion of StLecRLKs occurred majorly through tandem duplication although the whole-genome duplication (WGD)/segmental duplication events were found. Cis-elements in the StLecRLKs promoter region responded mainly to signals of defense and stress, phytohormone, biotic or abiotic stress. Moreover, expressional investigations indicated that the family members of the clades L-type, G-I, G-IV and G-V were responsive to both bacterial and fungal infection. Based on qRT-PCR analysis, the expressions of PGSC0003DMP400055136 and PGSC0003DMP400067047 were strongly induced in all treatments by both Fusarium sulphureum (Fs) and Phytophthora infestans (Pi) inoculation. The present study provides valuable information for LecRLKs gene family in potato genome, and establishes a foundation for further research into the functional analysis.

Light-harvesting chlorophyll a/b-binding protein-coding genes in jatropha and the comparison with castor, cassava and arabidopsis

The Lhc (light-harvesting chlorophyll a/b-binding protein) superfamily represents a class of antennae proteins that play indispensable roles in capture of solar energy as well as photoprotection under stress conditions. Despite their importance, little information has been available beyond model plants. In this study, we presents a first genome-wide analysis of Lhc superfamily genes in jatropha (Jatropha curcas L., Euphorbiaceae), an oil-bearing plant for biodiesel purpose. A total of 27 members were identified from the jatropha genome, which were shown to distribute over nine out of the 11 chromosomes. The superfamily number is comparable to 28 present in castor (Ricinus communis, Euphorbiaceae), but relatively less than 35 in cassava (Manihot esculenta, Euphorbiaceae) and 34 in arabidopsis (Arabidopsis thaliana) that experienced one or two recent whole-genome duplications (WGDs), respectively. In contrast to a high number of paralogs present in cassava and arabidopsis, few duplicates were found in jatropha as observed in castor, corresponding to no recent WGD occurred in these two species. Nevertheless, 26 orthologous groups representing four defined families were found in jatropha, and nearly one-to-one orthologous relationship was observed between jatropha and castor. By contrast, a novel group named SEP6 was shown to have been lost in arabidopsis. Global transcriptome profiling revealed a predominant expression pattern of most JcLhc superfamily genes in green tissues, reflecting their key roles in photosynthesis. Moreover, their expression profiles upon hormones, drought, and salt stresses were also investigated. These findings not only improve our knowledge on species-specific evolution of the Lhc supergene family, but also provide valuable information for further studies in jatropha.

A Chinese White Pear (Pyrus bretschneideri) BZR Gene PbBZR1 Act as a Transcriptional Repressor of Lignin Biosynthetic Genes in Fruits

BZR transcription factors play essential roles in plant growth and environmental stimuli, and they are also the positive regulators of Brassinosteroid (BR) signal transduction in diverse plants. In addition, BZR TFs, as crucial regulators of BR synthesis, may have multiple stress-resistance functions and their related regulatory mechanisms have been well illustrated in model plants. Here, we carried out a genome-wide identification of BZR members in Chinese pear (Pyrus bretschneideri) and identified 13 members. By comparative analysis in five Rosaceae genomes, BZR members in the pear genome may have undergone large-scale duplication events during evolution. Purifying selection played an important role in almost all of the orthologous and paralogous gene pairs. According to the expression analysis of the PbBZRs during fruit development, three PbBZRs were selected for detailed analysis. Transcriptional activation assays presented that PbBZR1 repressed the promoters of P. bretschneideri lignin biosynthetic genes, such as PbCES9, PbCOMT3, and PbHCT6. Our study traces the evolution of BZR gene family members in Rosaceae genomes and illustrates that the rates of gene loss and gain are far from equilibrium in different species. At the same time, our results suggest that PbBZR1 may be involved in the negative regulation of lignin biosynthesis.

Integrative Analysis of the Core Fruit Lignification Toolbox in Pear Reveals Targets for Fruit Quality Bioengineering

Stone cell content is an important factor affecting pear fruit flavor. Lignin, a major component of pear stone cells, hinders the quality and value of commercial fruit. The completion of the Chinese white pear (Pyrus bretschneideri) genome sequence provides an opportunity to perform integrative analysis of the genes encoding the eleven protein families (i.e., PAL, C4H, 4CL, HCT, C3H, CSE, CCoAOMT, CCR, F5H, COMT, and CAD) in the phenylpropanoid pathway. Here, a systematic study based on expression patterns and phylogenetic analyses was performed to identify the members of each gene family potentially involved in the lignification in the Chinese white pear. The phylogenetic analysis suggested that 35 P. bretschneideri genes belong to bona fide lignification clade members. Compared to other plants, some multigene families are expanded by tandem gene duplication, such as HCT, C3H, COMT, and CCR. RNA sequencing was used to study the expression patterns of the genes in different tissues, including leaf, petal, bud, sepal, ovary, stem, and fruit. Eighteen genes presented a high expression in fruit, indicating that these genes may be involved in the biosynthesis of lignin in pear fruit. Similarly to what has been observed for Populus trichocarpa, a bimolecular fluorescence complementation (BiFC) experiment indicated that P. bretschneideri C3H and C4H might also interact with each other to regulate monolignol biosynthesis in P. bretschneideri, ultimately affecting the stone cell content in pear fruits. The identification of the major genes involved in lignin biosynthesis in pear fruits provides the basis for the development of strategies to improve fruit quality.

Metacaspase gene family in Rosaceae genomes: Comparative genomic analysis and their expression during pear pollen tube and fruit development

Metacaspase (MC), which is discovered gene family with distant caspase homologs in plants, fungi, and protozoa, may be involved in programmed cell death (PCD) processes during plant development and respond abiotic and biotic stresses. To reveal the evolutionary relationship of MC gene family in Rosaceae genomes, we identified 8, 7, 8, 12, 12, and 23 MC genes in the genomes of Fragaria vesca, Prunus mume, Prunus persica, Pyrus communis, Pyrus bretschneideri and Malus domestica, respectively. Phylogenetic analysis suggested that the MC genes could be grouped into three clades: Type I*, Type I and Type II, which was supported by gene structure and conserved motif analysis. Microsynteny analysis revealed that MC genes present in the corresponding syntenic blocks of P. communis, P. bretschneideri and M. domestica, and further suggested that large-scale duplication events play an important role in the expansion of MC gene family members in these three genomes than other Rosaceae plants (F. vesca, P. mume and P. persica). RNA-seq data showed the specific expression patterns of PbMC genes in response to drought stress. The expression analysis of MC genes demonstrated that PbMC01 and PbMC03 were able to be detected in all four pear pollen tubes and seven fruit development stages. The current study highlighted the evolutionary relationship and duplication of the MC gene family in these six Rosaceae genomes and provided appropriate candidate genes for further studies in P. bretschneideri.

Genome-wide identification and comparative evolutionary analysis of the Dof transcription factor family in physic nut and castor bean

DNA-binding with one finger (Dof) proteins comprise a plant-specific transcription factor family involved in plant growth, development and stress responses. This study presents a genome-wide comparison of Dof family genes in physic nut (Jatropha curcas) and castor bean (Ricinus communis), two Euphorbiaceae plants that have not experienced any recent whole-genome duplication. A total of 25 or 24 Dof genes were identified from physic nut and castor genomes, respectively, where JcDof genes are distributed across nine out of 11 chromosomes. Phylogenetic analysis assigned these genes into nine groups representing four subfamilies, and 24 orthologous groups were also proposed based on comparison of physic nut, castor, Arabidopsis and rice Dofs. Conserved microsynteny was observed between physic nut and castor Dof-coding scaffolds, which allowed anchoring of 23 RcDof genes to nine physic nut chromosomes. In contrast to how no recent duplicate was present in castor, two tandem duplications and one gene loss were found in the Dof gene family of physic nut. Global transcriptome profiling revealed diverse patterns of Jc/RcDof genes over various tissues, and key Dof genes involved in flower development and stress response were also identified in physic nut. These findings provide valuable information for further studies of Dof genes in physic nut and castor.

Genome-wide identification and comparative analysis of phosphate starvation-responsive transcription factors in maize and three other gramineous plants

The present study identified several important candidate Pi regulation genes of maize and provides a better understanding on the generation of PHR genes in gramineous plants. Plants have evolved adaptive responses to cope with low phosphate (Pi) soils. The previous studies have indicated that phosphate starvation response (PHR) genes play central roles in regulating plant Pi starvation responses. However, the investigation of PHR family in gramineous plants is limited. In this study, we identified 64 PHR genes in four gramineous plants, including maize, rice, sorghum, and brachypodium, and conducted systematical analyses on phylogenetic, structure, collinearity, and expression pattern of these PHR genes. Genome synteny analysis revealed that a number of PHR genes were present in the corresponding syntenic blocks of maize, rice, sorghum, and brachypodium, indicating that large-scale duplication events contributed significantly to the expansion and evolution of PHR genes in these gramineous plants. Gene expression analysis showed that many PHR genes were expressed in various tissues, suggesting that these genes are involved in Pi redistribution and allocation. In addition, the expression levels of PHR genes from maize and rice under low Pi stress conditions revealed that some PHRs may play an important role in Pi starvation response. Our results provided a better understanding on the generation of PHR genes in gramineous plants and identified several important candidate Pi regulation genes of maize.

Characterization and expression of the ABC family (G group) in 'Dangshansuli' pear (Pyrus bretschneideri Rehd.) and its russet mutant

The plant genes encoding ABCGs that have been identified to date play a role in suberin formation in response to abiotic and biotic stress. In the present study, 80 ABCG genes were identified in ‘Dangshansuli’ Chinese white pear and designated as PbABCGs. Based on the structural characteristics and phylogenetic analysis, the PbABCG family genes could be classified into seven main groups: classes A-G. Segmental and dispersed duplications were the primary forces underlying the PbABCG gene family expansion in ‘Dangshansuli’ pear. Most of the PbABCG duplicated gene pairs date to the recent whole-genome duplication that occurred 30~45 million years ago. Purifying selection has also played a critical role in the evolution of the ABCG genes. Ten PbABCG genes screened in the transcriptome of ‘Dangshansuli’ pear and its russet mutant ‘Xiusu’ were validated, and the expression levels of the PbABCG genes exhibited significant differences at different stages. The results presented here will undoubtedly be useful for better understanding of the complexity of the PbABCG gene family and will facilitate the functional characterization of suberin formation in the russet mutant.

Retention and Molecular Evolution of Lipoxygenase Genes in Modern Rosid Plants

Whole-genome duplication events have occurred more than once in the genomes of some rosids and played a significant role over evolutionary time. Lipoxygenases (LOXs) are involved in many developmental and resistance processes in plants. Our study concerns the subject of the LOX gene family; we tracked the evolutionary process of ancestral LOX genes in four modern rosids. Here we show that some members of the LOX gene family in the Arabidopsis genome are likely to be lost during evolution, leading to a smaller size than that in Populus, Vitis, and Carica. Strong purifying selection acted as a critical role in almost all of the paralogous and orthologous genes. The structure of LOX genes in Carica and Populus are relatively stable, whereas Vitis and Arabidopsis have a difference. By searching conserved motifs of LOX genes, we found that each sub-family shared similar components. Research on intraspecies gene collinearity show that recent duplication holds an important position in Populus and Arabidopsis. Gene collinearity analysis within and between these four rosid plants revealed that all LOX genes in each modern rosid were the offspring from different ancestral genes. This study traces the evolution of LOX genes which have been differentially retained and expanded in rosid plants. Our results presented here may aid in the selection of special genes retained in the rosid plants for further analysis of biological function.

Evolutionary and functional analysis of mulberry type III polyketide synthases

Background

Type III polyketide synthases are important for the biosynthesis of flavonoids and various plant polyphenols. Mulberry plants have abundant polyphenols, but very little is known about the mulberry type III polyketide synthase genes. An analysis of these genes may provide new targets for genetic improvement to increase relevant secondary metabolites and enhance the plant tolerance to biotic and abiotic stresses.

Results

Eighteen genes encoding type III polyketide synthases were identified, including six chalcone synthases (CHS), ten stilbene synthases (STS), and two polyketide synthases (PKS). Functional characterization of four genes representing most of the MnCHS and MnSTS genes by coexpression with 4-Coumaroyl-CoA ligase in Escherichia coli indicated that their products were able to catalyze p-coumaroyl-CoA and malonyl-CoA to generate naringenin and resveratrol, respectively. Microsynteny analysis within mulberry indicated that segmental and tandem duplication events contributed to the expansion of the MnCHS family, while tandem duplications were mainly responsible for the generation of the MnSTS genes. Combining the evolution and expression analysis results of the mulberry type III PKS genes indicated that MnCHS and MnSTS genes evolved mainly under purifying selection to maintain their original functions, but transcriptional subfunctionalization occurred during long-term species evolution. Moreover, mulberry leaves can rapidly accumulated oxyresveratrol after UV-C irradiation, suggesting that resveratrol was converted to oxyresveratrol.

Conclusions

Characterizing the functions and evolution of mulberry type III PKS genes is crucial for advancing our understanding of these genes and providing the basis for further studies on the biosynthesis of relevant secondary metabolites in mulberry plants.

Electronic supplementary material

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

Molecular evolution and expression divergence of the Populus euphratica Hsf genes provide insight into the stress acclimation of desert poplar

Heat shock transcription factor (Hsf) family is one of the most important regulators in the plant kingdom. Hsf has been demonstrated to be involved in various processes associated with plant growth, development as well as in response to hormone and abiotic stresses. In this study, we carried out a comprehensive analysis of Hsf family in desert poplar, Populus euphratica. Total of 32 genes encoding Hsf were identified and they were classified into three main classes (A, B, and C). Gene structure and conserved motif analyses indicated that the members in each class were relatively conserved. Total of 10 paralogous pairs were identified in PeuHsf family, in which nine pairs were generated by whole genome duplication events. Ka/Ks analysis showed that PeuHsfs underwent purifying selection pressure. In addition, various cis-acting elements involved in hormone and stress responses located in the promoter regions of PeuHsfs. Gene expression analysis indicated that several PeuHsfs were tissue-specific expression. Compared to Arabidopsis, more PeuHsf genes were significantly induced by heat, drought, and salt stresses (21, 19, and 22 PeuHsfs, respectively). Our findings are helpful in understanding the distinguished adaptability of P. euphratica to extreme environment and providing a basis for functional analysis of PeuHsfs in the future.

Characterization and Comparison of the CPK Gene Family in the Apple (Malus × domestica) and Other Rosaceae Species and Its Response to Alternaria alternata Infection

As one of the Ca²⁺ sensors, calcium-dependent protein kinase (CPK) plays vital roles in immune and stress signaling, growth and development, and hormone responses, etc. Recently, the whole genome of apple (Malus × domestica), pear (Pyrus communis), peach (Prunus persica), plum (Prunus mume) and strawberry (Fragaria vesca) in Rosaceae family has been fully sequenced. However, little is known about the CPK gene family in these Rosaceae species. In this study, 123 CPK genes were identified from five Rosaceae species, including 37 apple CPKs, 37 pear CPKs, 17 peach CPKs, 16 strawberry CPKs, and 16 plum CPKs. Based on the phylogenetic tree topology and structural characteristics, we divided the CPK gene family into 4 distinct subfamilies: Group I, II, III, and IV. Whole-genome duplication (WGD) or segmental duplication played vital roles in the expansion of the CPK in these Rosaceae species. Most of segmental duplication pairs in peach and plum may have arisen from the γ triplication (~140 million years ago [MYA]), while in apple genome, many duplicated genes may have been derived from a recent WGD (30~45 MYA). Purifying selection also played a critical role in the function evolution of CPK family genes. Expression of apple CPK genes in response to apple pathotype of Alternaria alternata was verified by analysis of quantitative real-time RT-PCR (qPCR). Expression data demonstrated that CPK genes in apple might have evolved independently in different biological contexts. The analysis of evolution history and expression profile laid a foundation for further examining the function and complexity of the CPK gene family in Rosaceae.

VGSC: A Web-Based Vector Graph Toolkit of Genome Synteny and Collinearity

BACKGROUND

In order to understand the colocalization of genetic loci amongst species, synteny and collinearity analysis is a frequent task in comparative genomics research. However many analysis software packages are not effective in visualizing results. Problems include lack of graphic visualization, simple representation, or inextensible format of outputs. Moreover, higher throughput sequencing technology requires higher resolution image output.

IMPLEMENTATION

To fill this gap, this paper publishes VGSC, the Vector Graph toolkit of genome Synteny and Collinearity, and its online service, to visualize the synteny and collinearity in the common graphical format, including both raster (JPEG, Bitmap, and PNG) and vector graphic (SVG, EPS, and PDF).

RESULT

Users can upload sequence alignments from blast and collinearity relationship from the synteny analysis tools. The website can generate the vector or raster graphical results automatically. We also provide a java-based bytecode binary to enable the command-line execution.

Genome-wide analysis of WRKY transcription factors in white pear (Pyrus bretschneideri) reveals evolution and patterns under drought stress

Background

WRKY transcription factors (TFs) constitute one of the largest protein families in higher plants, and its members contain one or two conserved WRKY domains, about 60 amino acid residues with the WRKYGQK sequence followed by a C₂H₂ or C₂HC zinc finger motif. WRKY proteins play significant roles in plant development, and in responses to biotic and abiotic stresses. Pear (Pyrus bretschneideri) is one of the most important fruit crops in the world and is frequently threatened by abiotic stress, such as drought, affecting growth, development and productivity. Although the pear genome sequence has been released, little is known about the WRKY TFs in pear, especially in respond to drought stress at the genome-wide level.

Results

We identified a total of 103 WRKY TFs in the pear genome. Based on the structural features of WRKY proteins and topology of the phylogenetic tree, the pear WRKY (PbWRKY) family was classified into seven groups (Groups 1, 2a–e, and 3). The microsyteny analysis indicated that 33 (32 %) PbWRKY genes were tandemly duplicated and 57 genes (55.3 %) were segmentally duplicated. RNA-seq experiment data and quantitative real-time reverse transcription PCR revealed that PbWRKY genes in different groups were induced by drought stress, and Group 2a and 3 were mainly involved in the biological pathways in response to drought stress. Furthermore, adaptive evolution analysis detected a significant positive selection for Pbr001425 in Group 3, and its expression pattern differed from that of other members in this group. The present study provides a solid foundation for further functional dissection and molecular evolution of WRKY TFs in pear, especially for improving the water-deficient resistance of pear through manipulation of the PbWRKYs.

Electronic supplementary material

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

SynFind: Compiling Syntenic Regions across Any Set of Genomes on Demand

The identification of conserved syntenic regions enables discovery of predicted locations for orthologous and homeologous genes, even when no such gene is present. This capability means that synteny-based methods are far more effective than sequence similarity-based methods in identifying true-negatives, a necessity for studying gene loss and gene transposition. However, the identification of syntenic regions requires complex analyses which must be repeated for pairwise comparisons between any two species. Therefore, as the number of published genomes increases, there is a growing demand for scalable, simple-to-use applications to perform comparative genomic analyses that cater to both gene family studies and genome-scale studies. We implemented SynFind, a web-based tool that addresses this need. Given one query genome, SynFind is capable of identifying conserved syntenic regions in any set of target genomes. SynFind is capable of reporting per-gene information, useful for researchers studying specific gene families, as well as genome-wide data sets of syntenic gene and predicted gene locations, critical for researchers focused on large-scale genomic analyses. Inference of syntenic homologs provides the basis for correlation of functional changes around genes of interests between related organisms. Deployed on the CoGe online platform, SynFind is connected to the genomic data from over 15,000 organisms from all domains of life as well as supporting multiple releases of the same organism. SynFind makes use of a powerful job execution framework that promises scalability and reproducibility. SynFind can be accessed at http://genomevolution.org/CoGe/SynFind.pl. A video tutorial of SynFind using Phytophthrora as an example is available at http://www.youtube.com/watch?v=2Agczny9Nyc.

Comparative genomic analysis of the WRKY III gene family in populus, grape, arabidopsis and rice

Background

WRKY III genes have significant functions in regulating plant development and resistance. In plant, WRKY gene family has been studied in many species, however, there still lack a comprehensive analysis of WRKY III genes in the woody plant species poplar, three representative lineages of flowering plant species are incorporated in most analyses: Arabidopsis (a model plant for annual herbaceous dicots), grape (one model plant for perennial dicots) and Oryza sativa (a model plant for monocots).

Results

In this study, we identified 10, 6, 13 and 28 WRKY III genes in the genomes of Populus trichocarpa, grape (Vitis vinifera), Arabidopsis thaliana and rice (Oryza sativa), respectively. Phylogenetic analysis revealed that the WRKY III proteins could be divided into four clades. By microsynteny analysis, we found that the duplicated regions were more conserved between poplar and grape than Arabidopsis or rice. We dated their duplications by Ks analysis of Populus WRKY III genes and demonstrated that all the blocks were formed after the divergence of monocots and dicots. Strong purifying selection has played a key role in the maintenance of WRKY III genes in Populus. Tissue expression analysis of the WRKY III genes in Populus revealed that five were most highly expressed in the xylem. We also performed quantitative real-time reverse transcription PCR analysis of WRKY III genes in Populus treated with salicylic acid, abscisic acid and polyethylene glycol to explore their stress-related expression patterns.

Conclusions

This study highlighted the duplication and diversification of the WRKY III gene family in Populus and provided a comprehensive analysis of this gene family in the Populus genome. Our results indicated that the majority of WRKY III genes of Populus was expanded by large-scale gene duplication. The expression pattern of PtrWRKYIII gene identified that these genes play important roles in the xylem during poplar growth and development, and may play crucial role in defense to drought stress. Our results presented here may aid in the selection of appropriate candidate genes for further characterization of their biological functions in poplar.

Reviewers

This article was reviewed by Prof Dandekar and Dr Andrade-Navarro.

Electronic supplementary material

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

Comparative genomic analysis of the WRKY III gene family in populus, grape, arabidopsis and rice

BACKGROUND

WRKY III genes have significant functions in regulating plant development and resistance. In plant, WRKY gene family has been studied in many species, however, there still lack a comprehensive analysis of WRKY III genes in the woody plant species poplar, three representative lineages of flowering plant species are incorporated in most analyses: Arabidopsis (a model plant for annual herbaceous dicots), grape (one model plant for perennial dicots) and Oryza sativa (a model plant for monocots).

RESULTS

In this study, we identified 10, 6, 13 and 28 WRKY III genes in the genomes of Populus trichocarpa, grape (Vitis vinifera), Arabidopsis thaliana and rice (Oryza sativa), respectively. Phylogenetic analysis revealed that the WRKY III proteins could be divided into four clades. By microsynteny analysis, we found that the duplicated regions were more conserved between poplar and grape than Arabidopsis or rice. We dated their duplications by Ks analysis of Populus WRKY III genes and demonstrated that all the blocks were formed after the divergence of monocots and dicots. Strong purifying selection has played a key role in the maintenance of WRKY III genes in Populus. Tissue expression analysis of the WRKY III genes in Populus revealed that five were most highly expressed in the xylem. We also performed quantitative real-time reverse transcription PCR analysis of WRKY III genes in Populus treated with salicylic acid, abscisic acid and polyethylene glycol to explore their stress-related expression patterns.

CONCLUSIONS

This study highlighted the duplication and diversification of the WRKY III gene family in Populus and provided a comprehensive analysis of this gene family in the Populus genome. Our results indicated that the majority of WRKY III genes of Populus was expanded by large-scale gene duplication. The expression pattern of PtrWRKYIII gene identified that these genes play important roles in the xylem during poplar growth and development, and may play crucial role in defense to drought stress. Our results presented here may aid in the selection of appropriate candidate genes for further characterization of their biological functions in poplar.

Comparative genomic analysis of the WRKY III gene family in populus, grape, arabidopsis and rice

BACKGROUND

WRKY III genes have significant functions in regulating plant development and resistance. In plant, WRKY gene family has been studied in many species, however, there still lack a comprehensive analysis of WRKY III genes in the woody plant species poplar, three representative lineages of flowering plant species are incorporated in most analyses: Arabidopsis (a model plant for annual herbaceous dicots), grape (one model plant for perennial dicots) and Oryza sativa (a model plant for monocots).

RESULTS

In this study, we identified 10, 6, 13 and 28 WRKY III genes in the genomes of Populus trichocarpa, grape (Vitis vinifera), Arabidopsis thaliana and rice (Oryza sativa), respectively. Phylogenetic analysis revealed that the WRKY III proteins could be divided into four clades. By microsynteny analysis, we found that the duplicated regions were more conserved between poplar and grape than Arabidopsis or rice. We dated their duplications by Ks analysis of Populus WRKY III genes and demonstrated that all the blocks were formed after the divergence of monocots and dicots. Strong purifying selection has played a key role in the maintenance of WRKY III genes in Populus. Tissue expression analysis of the WRKY III genes in Populus revealed that five were most highly expressed in the xylem. We also performed quantitative real-time reverse transcription PCR analysis of WRKY III genes in Populus treated with salicylic acid, abscisic acid and polyethylene glycol to explore their stress-related expression patterns.

CONCLUSIONS

This study highlighted the duplication and diversification of the WRKY III gene family in Populus and provided a comprehensive analysis of this gene family in the Populus genome. Our results indicated that the majority of WRKY III genes of Populus was expanded by large-scale gene duplication. The expression pattern of PtrWRKYIII gene identified that these genes play important roles in the xylem during poplar growth and development, and may play crucial role in defense to drought stress. Our results presented here may aid in the selection of appropriate candidate genes for further characterization of their biological functions in poplar.

Evolutionary Comparison of Two Combinatorial Regulators of SBP-Box Genes, MiR156 and MiR529, in Plants

A complete picture of the evolution of miRNA combinatorial regulation requires the synthesis of information on all miRNAs and their targets. MiR156 and miR529 are two combinatorial regulators of squamosa promoter binding protein-like (SBP-box) genes. Previous studies have clarified the evolutionary dynamics of their targets; however, there have been no reports on the evolutionary patterns of two miRNA regulators themselves to date. In this study, we investigated the evolutionary differences between these two miRNA families in extant land plants. Our work found that miR529 precursor, especially of its mature miRNA sequence, has a higher evolutionary rate. Such accelerating evolution of miR529 has significantly effects on its structural stability, and sequence conservation against existence of itself. By contrast, miR156 evolves more rapidly in loop region of the stable secondary structure, which may contribute to its functional diversity. Moreover, miR156 and miR529 genes have distinct rates of loss after identical duplication events. MiR529 genes have a higher average loss rate and asymmetric loss rate in duplicated gene pairs, indicating preferred miR529 gene losses become another predominant mode of inactivation, that are implicated in the contraction of this family. On the contrary, duplicated miR156 genes have a low loss rate, and could serve as another new source for functional diversity. Taken together, these results provide better insight into understanding the evolutionary divergence of miR156 and miR529 family in miRNA combinational regulation network.

Genome-wide analysis of Dof family transcription factors and their responses to abiotic stresses in Chinese cabbage

BACKGROUND

Chinese cabbage is an important leaf vegetable that experienced long-term cultivation and artificial selection. Dof (DNA-binding One Zinc Finger) transcription factors, with a highly conserved Dof domain, are members of a major plant-specific transcription factor family that play important roles in many plant biological processes. The Dof family transcription factors, one of the most important families of transcriptional regulators in higher plants, are involved in massive aspects of plant growth, development, and response to abiotic stresses. Our study will supply resources for understanding how Dof transcription factors respond to abiotic stress and the interaction network of these genes in tolerance mechanism.

RESULTS

In this study, we performed a comprehensive analysis of Dof family factors in Chinese cabbage. In total, 76 genes encoding BraDof family transcription factor were identified from Chinese cabbage, and those BraDof factors were divided into nine classes. Fifteen motifs were found based on Dof amino acid sequence alignments. Chromosome locations and gene duplications of BraDof family genes were also analyzed. Ten duplicate events of BraDof genes were discovered in Chinese cabbage chromosomes. The uneven distribution of BraDof genes in Brassica chromosomes may cause the expansion of BraDof genes. In the Dof family, 37 and 7 orthologous genes were identified between Chinese cabbage and Arabidopsis and between Chinese cabbage and Oryza sativa, respectively. The interaction networks of Dof factors in Chinese cabbage were also constructed. Expression profiles of nine selected genes from different nine classes subjected to four abiotic stresses (cold, heat, salt and drought) were further investigated by quantitative real-time PCR to obtain a better understanding of the functions and regulation mechanisms of BraDof family transcription factors in two Chinese cabbage varieties, 'Lubaisanhao' and 'Qingdao 87-114'.

CONCLUSIONS

Dof-family transcription factors were analyzed in genome of Chinese cabbage. Chromosomal locations showed that duplication might result in expansion. Response to abiotic stresses was elucidated in Chinese cabbage varieties. The results provide novel insights into the stress responses of BraDof genes and promote a better understanding of the construction and function of Dofs in Chinese cabbage.

MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity

MCScan is an algorithm able to scan multiple genomes or subgenomes in order to identify putative homologous chromosomal regions, and align these regions using genes as anchors. The MCScanX toolkit implements an adjusted MCScan algorithm for detection of synteny and collinearity that extends the original software by incorporating 14 utility programs for visualization of results and additional downstream analyses. Applications of MCScanX to several sequenced plant genomes and gene families are shown as examples. MCScanX can be used to effectively analyze chromosome structural changes, and reveal the history of gene family expansions that might contribute to the adaptation of lineages and taxa. An integrated view of various modes of gene duplication can supplement the traditional gene tree analysis in specific families. The source code and documentation of MCScanX are freely available at http://chibba.pgml.uga.edu/mcscan2/.