Comparative analysis of GT14/GT14-like gene family in Arabidopsis, Oryza, Populus, Sorghum and Vitis

Glycosyltransferase family14 (GT14) belongs to the glycosyltransferase (GT) superfamily that plays important roles in the biosynthesis of cell walls, the most abundant source of cellulosic biomass for bioethanol production. It has been hypothesized that DUF266 proteins are a new class of GTs related to GT14. In this study, we identified 62 GT14 and 106 DUF266 genes (named GT14-like herein) in Arabidopsis, Oryza, Populus, Sorghum and Vitis. Our phylogenetic analysis separated GT14 and GT14-like genes into two distinct clades, which were further divided into eight and five groups, respectively. Similarities in protein domain, 3D structure and gene expression were uncovered between the two phylogenetic clades, supporting the hypothesis that GT14 and GT14-like genes belong to one family. Therefore, we proposed a new family name, GT14/GT14-like family that combines both subfamilies. Variation in gene expression and protein subcellular localization within the GT14-like subfamily were greater than those within the GT14 subfamily. One-half of the Arabidopsis and Populus GT14/GT14-like genes were found to be preferentially expressed in stem/xylem, indicating that they are likely involved in cell wall biosynthesis. This study provided new insights into the evolution and functional diversification of the GT14/GT14-like family genes.     

Identification of candidate genes in Arabidopsis and Populus cell wall biosynthesis using text-mining, co-expression network analysis and comparative genomics

Populus is an important bioenergy crop for bioethanol production. A greater understanding of cell wall biosynthesis processes is critical in reducing biomass recalcitrance, a major hindrance in efficient generation of biofuels from lignocellulosic biomass. Here, we report the identification of candidate cell wall biosynthesis genes through the development and application of a novel bioinformatics pipeline. As a first step, via text-mining of PubMed publications, we obtained 121 Arabidopsis genes that had the experimental evidence supporting their involvement in cell wall biosynthesis or remodeling. The 121 genes were then used as bait genes to query an Arabidopsis co-expression database, and additional genes were identified as neighbors of the bait genes in the network, increasing the number of genes to 548. The 548 Arabidopsis genes were then used to re-query the Arabidopsis co-expression database and re-construct a network that captured additional network neighbors, expanding to a total of 694 genes. The 694 Arabidopsis genes were computationally divided into 22 clusters. Queries of the Populus genome using the Arabidopsis genes revealed 817 Populus orthologs. Functional analysis of gene ontology and tissue-specific gene expression indicated that these Arabidopsis and Populus genes are high likelihood candidates for functional characterization in relation to cell wall biosynthesis.     

A physical map for the Amborella trichopoda genome sheds light on the evolution of angiosperm genome structure

Background

Recent phylogenetic analyses have identified Amborella trichopoda, an understory tree species endemic to the forests of New Caledonia, as sister to a clade including all other known flowering plant species. The Amborella genome is a unique reference for understanding the evolution of angiosperm genomes because it can serve as an outgroup to root comparative analyses. A physical map, BAC end sequences and sample shotgun sequences provide a first view of the 870 Mbp Amborella genome.

Results

Analysis of Amborella BAC ends sequenced from each contig suggests that the density of long terminal repeat retrotransposons is negatively correlated with that of protein coding genes. Syntenic, presumably ancestral, gene blocks were identified in comparisons of the Amborella BAC contigs and the sequenced Arabidopsis thaliana, Populus trichocarpa, Vitis vinifera and Oryza sativa genomes. Parsimony mapping of the loss of synteny corroborates previous analyses suggesting that the rate of structural change has been more rapid on lineages leading to Arabidopsis and Oryza compared with lineages leading to Populus and Vitis. The gamma paleohexiploidy event identified in the Arabidopsis, Populus and Vitis genomes is shown to have occurred after the divergence of all other known angiosperms from the lineage leading to Amborella.

Conclusions

When placed in the context of a physical map, BAC end sequences representing just 5.4% of the Amborella genome have facilitated reconstruction of gene blocks that existed in the last common ancestor of all flowering plants. The Amborella genome is an invaluable reference for inferences concerning the ancestral angiosperm and subsequent genome evolution.     

Genome-wide identification and expression analysis of the mitogen-activated protein kinase gene family from banana suggest involvement of specific members in different stages of fruit ripening

Mitogen-activated protein kinases (MAPKs) are important components of the tripartite mitogen-activated protein kinase signaling cascade and play an important role in plant growth and development. Although members of the MAPK gene family have been identified in model plants, little information is available regarding this gene family in fruit crops. In this study, we carried out a computational analysis using the Musa Genome database to identify members of the MAPK gene family in banana, an economically important crop and the most popular fruit worldwide. Our analysis identified 25 members of the MAP kinase (MAPK or MPK) gene family. Phylogenetic analyses of MPKs in Arabidopsis, Oryza, and Populus have classified these MPKs into four subgroups. The presence of conserved domains in the deduced amino acid sequences, phylogeny, and genomic organization strongly support their identity as members of the MPK gene family. Expression analysis during ethylene-induced banana fruit ripening suggests the involvement of several MPKs in the ethylene signal transduction pathway that are necessary for banana fruit ripening. Analysis of the cis-regulatory elements in the promoter regions and the involvement of the identified MPKs in various cellular processes, as analyzed using Pathway Studio, suggest a role for the banana MPK gene family in diverse functions related to growth, development, and the stress response. This report is the first concerning the identification of members of a gene family and the elucidation of their role in various processes using the Musa Genome database.     

GST-PRIME: a genome-wide primer design software for the generation of gene sequence tags

The availability of sequenced genomes has generated a need for experimental approaches that allow the simultaneous analysis of large, or even complete, sets of genes. To facilitate such analyses, we have developed GST-PRIME, a software package for retrieving and assembling gene sequences, even from complex genomes, using the NCBI public database, and then designing sets of primer pairs for use in gene amplification. Primers were designed by the program for the direct amplification of gene sequence tags (GSTs) from either genomic DNA or cDNA. Test runs of GST-PRIME on 2000 randomly selected Arabidopsis and Drosophila genes demonstrate that 93 and 88% of resulting GSTs, respectively, fulfilled imposed length criteria. GST-PRIME primer pairs were tested on a set of 1900 Arabidopsis genes coding for chloroplast-targeted proteins: 95% of the primer pairs used in PCRs with genomic DNA generated the correct amplicons. GST-PRIME can thus be reliably used for large-scale or specific amplification of intron-containing genes of multicellular eukaryotes.     

The F-box gene family is expanded in herbaceous annual plants relative to woody perennial plants

F-box proteins are generally responsible for substrate recognition in the Skp1-Cullin-F-box complexes that are involved in protein degradation via the ubiquitin-26S proteasome pathway. In plants, F-box genes influence a variety of biological processes, such as leaf senescence, branching, self-incompatibility, and responses to biotic and abiotic stresses. The number of F-box genes in Populus (Populus trichocarpa; approximately 320) is less than half that found in Arabidopsis (Arabidopsis thaliana; approximately 660) or Oryza (Oryza sativa; approximately 680), even though the total number of genes in Populus is equivalent to that in Oryza and 1.5 times that in Arabidopsis. We performed comparative genomics analysis between the woody perennial plant Populus and the herbaceous annual plants Arabidopsis and Oryza in order to explicate the functional implications of this large gene family. Our analyses reveal interspecific differences in genomic distribution, orthologous relationship, intron evolution, protein domain structure, and gene expression. The set of F-box genes shared by these species appear to be involved in core biological processes essential for plant growth and development; lineage-specific differences primarily occurred because of an expansion of the F-box genes via tandem duplications in Arabidopsis and Oryza. The number of F-box genes in the newly sequenced woody species Vitis (Vitis vinifera; 156) and Carica (Carica papaya; 139) is similar to that in Populus, supporting the hypothesis that the F-box gene family is expanded in herbaceous annual plants relative to woody perennial plants. This study provides insights into the relationship between the structure and composition of the F-box gene family in herbaceous and woody species and their associated developmental and physiological features.     

Differential retention and expansion of the ancestral genes associated with the paleopolyploidies in modern rosid plants,as revealed by analysis of the extensins super-gene family

Background

All modern rosids originated from a common hexapolyploid ancestor, and the genomes of some rosids have undergone one or more cycles of paleopolyploidy. After the duplication of the ancient genome, wholesale gene loss and gene subfunctionalization has occurred. Using the extensin super-gene family as an example, we tracked the differential retention and expansion of ancestral extensin genes in four modern rosids, Arabidopsis, Populus, Vitis and Carica, using several analytical methods.

Results

The majority of extensin genes in each of the modern rosids were found to originate from different ancestral genes. In Arabidopsis and Populus, almost half of the extensins were paralogous duplicates within the genome of each species. By contrast, no paralogous extensins were detected in Vitis and Carica, which have only undergone the common γ-triplication event. It was noteworthy that a group of extensins containing the IPR006706 domain had actively duplicated in Arabidopsis, giving rise to a neo-extensin around every 3 million years. However, such extensins were absent from, or rare in, the other three rosids. A detailed examination revealed that this group of extensins had proliferated significantly in the genomes of a number of species in the Brassicaceae. We propose that this group of extensins might play important roles in the biology and in the evolution of the Brassicaceae. Our analyses also revealed that nearly all of the paralogous and orthologous extensin-pairs have been under strong purifying selection, leading to the strong conservation of the function of extensins duplicated from the same ancestral gene.

Conclusions

Our analyses show that extensins originating from a common ancestor have been differentially retained and expanded among four modern rosids. Our findings suggest that, if Arabidopsis is used as the model plant, we can only learn a limited amount about the functions of a particular gene family. These results also provide an example of how it is essential to learn the origination of a gene when analyzing its function across different plant species.

Electronic supplementary material

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

Identification of shared single copy nuclear genes in Arabidopsis, Populus, Vitis and Oryza and their phylogenetic utility across various taxonomic levels

Background

Although the overwhelming majority of genes found in angiosperms are members of gene families, and both gene- and genome-duplication are pervasive forces in plant genomes, some genes are sufficiently distinct from all other genes in a genome that they can be operationally defined as 'single copy'. Using the gene clustering algorithm MCL-tribe, we have identified a set of 959 single copy genes that are shared single copy genes in the genomes of Arabidopsis thaliana, Populus trichocarpa, Vitis vinifera and Oryza sativa. To characterize these genes, we have performed a number of analyses examining GO annotations, coding sequence length, number of exons, number of domains, presence in distant lineages, such as Selaginella and Physcomitrella, and phylogenetic analysis to estimate copy number in other seed plants and to demonstrate their phylogenetic utility. We then provide examples of how these genes may be used in phylogenetic analyses to reconstruct organismal history, both by using extant coverage in EST databases for seed plants and de novo amplification via RT-PCR in the family Brassicaceae.

Results

There are 959 single copy nuclear genes shared in Arabidopsis, Populus, Vitis and Oryza ["APVO SSC genes"]. The majority of these genes are also present in the Selaginella and Physcomitrella genomes. Public EST sets for 197 species suggest that most of these genes are present across a diverse collection of seed plants, and appear to exist as single or very low copy genes, though exceptions are seen in recently polyploid taxa and in lineages where there is significant evidence for a shared large-scale duplication event. Genes encoding proteins localized in organelles are more commonly single copy than expected by chance, but the evolutionary forces responsible for this bias are unknown. Regardless of the evolutionary mechanisms responsible for the large number of shared single copy genes in diverse flowering plant lineages, these genes are valuable for phylogenetic and comparative analyses. Eighteen of the APVO SSC single copy genes were amplified in the Brassicaceae using RT-PCR and directly sequenced. Alignments of these sequences provide improved resolution of Brassicaceae phylogeny compared to recent studies using plastid and ITS sequences. An analysis of sequences from 13 APVO SSC genes from 69 species of seed plants, derived mainly from public EST databases, yielded a phylogeny that was largely congruent with prior hypotheses based on multiple plastid sequences. Whereas single gene phylogenies that rely on EST sequences have limited bootstrap support as the result of limited sequence information, concatenated alignments result in phylogenetic trees with strong bootstrap support for already established relationships. Overall, these single copy nuclear genes are promising markers for phylogenetics, and contain a greater proportion of phylogenetically-informative sites than commonly used protein-coding sequences from the plastid or mitochondrial genomes.

Conclusions

Putatively orthologous, shared single copy nuclear genes provide a vast source of new evidence for plant phylogenetics, genome mapping, and other applications, as well as a substantial class of genes for which functional characterization is needed. Preliminary evidence indicates that many of the shared single copy nuclear genes identified in this study may be well suited as markers for addressing phylogenetic hypotheses at a variety of taxonomic levels.     

Genome Wide Analysis for Searching Novel Markers to Rapidly Identify Clostridium Strains

Microbial classification is based largely on the 16S rRNA (rrs) gene sequence, which is conserved throughout the prokaryotic domain. The Ribosomal Database Project (RDP) has become a reference point for almost all practical purposes. The use of this gene is limited by the fact that it can be used to identify only to the extent to what has been known and is available in the RDP. In order to identify an organism whose rrs is not present in the RDP database, we need to generate novel markers to place the unknown on the evolutionary map. Here, sequenced genomes of 27 Clostridium strains belonging to 9 species have been used to identify two sets of genes: (1) common to most of the species, and (2) unique to a species. Combinations of genes (recN, dnaJ, secA, mutS, and/or grpE) and their unique restriction endonuclease digestion (AluI, BfaI and/or Tru9I) patterns have been established to rapidly identify Clostridium species. This strategy for identifying novel markers can be extended to all other organisms and diagnostic applications.

Electronic supplementary material

The online version of this article (doi:10.1007/s12088-015-0535-7) contains supplementary material, which is available to authorized users.     

Significant microsynteny with new evolutionary highlights is detected between <Emphasis Type="Italic">Arabidopsis</Emphasis> and legume model plants despite the lack of macrosynteny

The increased amount of data produced by large genome sequencing projects allows scientists to carry out important syntenic studies to a great extent. Detailed genetic maps and entirely or partially sequenced genomes are compared, and macro- and microsyntenic relations can be determined for different species. In our study, the syntenic relationships between key legume plants and two model plants, Arabidopsis thaliana and Populus trichocarpa were investigated. The comparison of the map position of 172 gene-based Medicago sativa markers to the organization of homologous A. thaliana genes could not identify any sign of macrosynteny between the two genomes. A 276 kb long section of chromosome 5 of the model legume Medicago truncatula was used to investigate potential microsynteny with the other legume Lotus japonicus, as well as with Arabidopsis and Populus. Besides the overall correlation found between the legume plants, the comparison revealed several microsyntenic regions in the two more distant plants with significant resemblance. Despite the large phylogenetic distance, clear microsyntenic regions between Medicago and Arabidopsis or Populus were detected unraveling new intragenomic evolutionary relations in Arabidopsis. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Carotenoid biosynthesis genes in rice: structural analysis, genome-wide expression profiling and phylogenetic analysis

Carotenoids, important lipid-soluble antioxidants in photosynthetic tissues, are known to be completely absent in rice endosperm. Many studies, involving transgenic manipulations of carotenoid biosynthesis genes, have been performed to get carotenoid-enriched rice grain. Study of genes involved in their biosynthesis can provide further information regarding the abundance/absence of carotenoids in different tissues. We have identified 16 and 34 carotenoid biosynthesis genes in rice and Populus genomes, respectively. A detailed analysis of the domain structure of carotenoid biosynthesis enzymes in rice, Populus and Arabidopsis has shown that highly conserved catalytic domains, along with other domains, are present in these proteins. Phylogenetic analysis of rice genes with Arabidopsis and other characterized carotenoid biosynthesis genes has revealed that homologous genes exist in these plants, and the duplicated gene copies probably adopt new functions. Expression of rice and Populus genes has been analyzed by full-length cDNA- and EST-based expression profiling. In rice, this analysis was complemented by real-time PCR, microarray and signature-based expression profiling, which reveal that carotenoid biosynthesis genes are highly expressed in light-grown tissues, have differential expression pattern during vegetative/reproductive development and are responsive to stress.     

Genomic structure and evolution of the Pi2/9 locus in wild rice species

Rice blast, caused by the fungal pathogen Magnaporthe oryzae, is a devastating disease of rice worldwide. Among the 85 mapped resistance (R) genes against blast, 13 have been cloned and characterized. However, how these genes originated and how they evolved in the Oryza genus remains unclear. We previously cloned the rice blast R-genes Pi2, Pi9, and Piz-t, and analyzed their genomic structure and evolution in cultivated rice. In this study, we determined the genomic sequences of the Pi2/9 locus in four wild Oryza species representing three genomes (AA, BB and CC). The number of Pi2/9 family members in the four wild species ranges from two copies to 12 copies. Although these genes are conserved in structure and categorized into the same subfamily, sequence duplications and subsequent inversions or uneven crossing overs were observed, suggesting that the locus in different wild species has undergone dynamic changes. Positive selection was found in the leucine-rich repeat region of most members, especially in the largest clade where Pi9 is included. We also provide evidence that the Pi9 gene is more related to its homologues in the recurrent line and other rice cultivars than to those in its alleged donor species O. minuta, indicating a possible origin of the Pi9 gene from O. sativa. Comparative sequence analysis between the four wild Oryza species and the previously established reference sequences in cultivated rice species at the Pi2/9 locus has provided extensive and unique information on the genomic structure and evolution of a complex R-gene cluster in the Oryza genus.     

The complete nucleotide sequence and RNA editing content of the mitochondrial genome of rapeseed (Brassica napus L.): comparative analysis of the mitochondrial genomes of rapeseed and Arabidopsis thaliana

The entire mitochondrial genome of rapeseed (Brassica napus L.) was sequenced and compared with that of Arabidopsis thaliana. The 221 853 bp genome contains 34 protein-coding genes, three rRNA genes and 17 tRNA genes. This gene content is almost identical to that of Arabidopsis. However the rps14 gene, which is a pseudo-gene in Arabidopsis, is intact in rapeseed. On the other hand, five tRNA genes are missing in rapeseed compared to Arabidopsis, although the set of mitochondrially encoded tRNA species is identical in the two Cruciferae. RNA editing events were systematically investigated on the basis of the sequence of the rapeseed mitochondrial genome. A total of 427 C to U conversions were identified in ORFs, which is nearly identical to the number in Arabidopsis (441 sites). The gene sequences and intron structures are mostly conserved (more than 99% similarity for protein-coding regions); however, only 358 editing sites (83% of total editings) are shared by rapeseed and Arabidopsis. Non-coding regions are mostly divergent between the two plants. One-third (about 78.7 kb) and two-thirds (about 223.8 kb) of the rapeseed and Arabidopsis mitochondrial genomes, respectively, cannot be aligned with each other and most of these regions do not show any homology to sequences registered in the DNA databases. The results of the comparative analysis between the rapeseed and Arabidopsis mitochondrial genomes suggest that higher plant mitochondria are extremely conservative with respect to coding sequences and somewhat conservative with respect to RNA editing, but that non-coding parts of plant mitochondrial DNA are extraordinarily dynamic with respect to structural changes, sequence acquisition and/or sequence loss.     

Genome-Wide Analysis in Wild and Cultivated <Emphasis Type="Italic">Oryza</Emphasis> Species Reveals Abundance of NBS Genes in Progenitors of Cultivated Rice

NBS-encoding genes play a critical role in the plant defense system. Wild relatives of crop plants are rich reservoirs of plant defense genes. Here, we performed a stringent genome-wide identification of NBS-encoding genes in three cultivated and eight wild Oryza species, representing three different genomes (AA, BB, and FF) from four continents. A total of 2688 NBS-encoding genes were identified from 11 Oryza genomes. All the three progenitor species of cultivated rice, namely O. barthii, O. rufipogon, and O. nivara, were the richest reservoir of NBS-encoding genes (214, 313, and 307 respectively). Interestingly, the two Asian cultivated species showed a contrasting pattern in the number of NBS-encoding genes. While indica subspecies maintained nearly equal number of NBS genes as its progenitor (309 and 313), the japonica subspecies had retained only two third in the course of evolution (213 and 307). Other major sources for NBS-encoding genes could be (i) O. longistaminata since it had the highest proportion of NBS-encoding genes and (ii) O. glumaepatula as it clustered distinctly away from the rest of the AA genome species. The present study thus revealed that NBS-encoding genes can be exploited from the primary gene pool for disease resistance breeding in rice.     

Structure of two melon regions reveals high microsynteny with sequenced plant species

In this study, two melon bacterial artificial chromosome (BAC) clones have been sequenced and annotated. BAC 1-21-10 spans 92 kb and contains the nsv locus conferring resistance to the Melon Necrotic Spot Virus (MNSV) in melon linkage group 11. BAC 13J4 spans 98 kb and belongs to a BAC contig containing resistance gene homologues, extending a previous sequenced region of 117 kb in linkage group 4. Both regions have microsyntenic relationships to the model plant species Arabidopsis thaliana, and to Medicago truncatula and Populus trichocarpa. The network of synteny found between melon and each of the sequenced genomes reflects the polyploid structure of Arabidopsis, Populus, and Medicago genomes due to whole genome duplications (WGD). A detailed analysis revealed that both melon regions have a lower relative syntenic quality with Arabidopsis (eurosid II) than when compared to Populus and Medicago (eurosid I). Although phylogenetically Cucurbitales seem to be closer to Fabales than to Malphigiales, synteny was higher between both melon regions and Populus. Presented data imply that the recently completed Populus genome sequence could preferentially be used to obtain positional information in melon, based on microsynteny. Wim Deleu and Víctor González contributed equally to this work. The nucleotide sequences of BACs 1-21-10 and 13J4 are available in the DDBJ/EMBL/GenBank databases under the accession numbers EF188258 and EF657230, respectively.     

Phylogeny and Molecular Evolution of miR820 and miR396 microRNA Families in <Emphasis Type="Italic">Oryza</Emphasis> AA Genomes

The phylogeny and evolution of the microRNA families, miR820 and miR396, was analysed across the AA genomes of the Oryza species, the close relatives of domesticated rice. A highly dynamic evolution of the miR820 family was revealed. The number of copies of MIR820 genes, their chromosomal location and the mature microRNA sequence varied greatly with a total of 16 novel miR820 variants being identified. The phylogeny of pre-MIR820 sequences revealed that MIR820 genes of recently evolved Oryza AA genomes may have derived from sequence divergence of one or a few ancestral genes found in wild Australian perennial rice populations, Taxon B (jpn2)-MIR820 genes. Genomic scale duplication played an important role in the evolution of some of the miR396 family genes in AA genome Oryza species. miR396 family contained a MIR396 gene cluster (MIR396a and MIR396c) which was conserved across the cereal genomes. Nucleotide diversity analysis at these two MIR396 loci revealed that domesticated rice has retained less than 10% of the total diversity present in wild species. In contrast, the nucleotide sequence of four MIR396 loci remained almost conserved across domesticated and wild rices, indicating that they were under extreme functional constraint and may be involved in regulating some fundamental processes which are important both for wild and domesticated rices. Expression analysis demonstrated that miR820 variants were expressed in O. glaberrima O. barthi and O. longistaminata genome. These findings pose new challenges to explain the possible role of miR820 variants identified.