Genome-Wide Analysis of the Expansin Gene Superfamily Reveals Grapevine-Specific Structural and Functional Characteristics

Background

Expansins are proteins that loosen plant cell walls in a pH-dependent manner, probably by increasing the relative movement among polymers thus causing irreversible expansion. The expansin superfamily (EXP) comprises four distinct families: expansin A (EXPA), expansin B (EXPB), expansin-like A (EXLA) and expansin-like B (EXLB). There is experimental evidence that EXPA and EXPB proteins are required for cell expansion and developmental processes involving cell wall modification, whereas the exact functions of EXLA and EXLB remain unclear. The complete grapevine (Vitis vinifera) genome sequence has allowed the characterization of many gene families, but an exhaustive genome-wide analysis of expansin gene expression has not been attempted thus far.

Methodology/Principal Findings

We identified 29 EXP superfamily genes in the grapevine genome, representing all four EXP families. Members of the same EXP family shared the same exon–intron structure, and phylogenetic analysis confirmed a closer relationship between EXP genes from woody species, i.e. grapevine and poplar (Populus trichocarpa), compared to those from Arabidopsis thaliana and rice (Oryza sativa). We also identified grapevine-specific duplication events involving the EXLB family. Global gene expression analysis confirmed a strong correlation among EXP genes expressed in mature and green/vegetative samples, respectively, as reported for other gene families in the recently-published grapevine gene expression atlas. We also observed the specific co-expression of EXLB genes in woody organs, and the involvement of certain grapevine EXP genes in berry development and post-harvest withering.

Conclusion

Our comprehensive analysis of the grapevine EXP superfamily confirmed and extended current knowledge about the structural and functional characteristics of this gene family, and also identified properties that are currently unique to grapevine expansin genes. Our data provide a model for the functional characterization of grapevine gene families by combining phylogenetic analysis with global gene expression profiling.     

Expansin gene loss is a common occurrence during adaptation to an aquatic environment

Expansins comprise a superfamily of plant cell wall loosening proteins that can be divided into four individual families (EXPA, EXPB, EXLA and EXLB). Aside from inferred roles in a variety of plant growth and developmental traits, little is known regarding the function of specific expansin clades, for which there are at least 16 in flowering plants (angiosperms); however, there is evidence to suggest that some expansins have cell‐specific functions, in root hair and pollen tube development, for example. Recently, two duckweed genomes have been sequenced (Spirodela polyrhiza strains 7498 and 9509), revealing significantly reduced superfamily sizes. We hypothesized that there would be a correlation between expansin loss and morphological reductions seen among highly adapted aquatic species. In order to provide an answer to this question, we characterized the expansin superfamilies of the greater duckweed Spirodela, the marine eelgrass Zostera marina and the bladderwort Utricularia gibba. We discovered rampant expansin gene and clade loss among the three, including a complete absence of the EXLB family and EXPA‐VII. The most convincing correlation between morphological reduction and expansin loss was seen for Utricularia and Spirodela, which both lack root hairs and the root hair expansin clade EXPA‐X. Contrary to the pattern observed in other species, four Utricularia expansins failed to branch within any clade, suggesting that they may be the result of neofunctionalization. Last, an expansin clade previously discovered only in eudicots was identified in Spirodela, allowing us to conclude that the last common ancestor of monocots and eudicots contained a minimum of 17 expansins.     

Genome histories clarify evolution of the expansin superfamily: new insights from the poplar genome and pine ESTs

Expansins comprise a superfamily of plant cell wall-loosening proteins that has been divided into four distinct families, EXPA, EXPB, EXLA and EXLB. In a recent analysis of Arabidopsis thaliana and Oryza sativa expansins, we proposed a further subdivision of the families into 17 clades, representing independent lineages in the last common ancestor of monocots and eudicots. This division was based on both traditional sequence-based phylogenetic trees and on position-based trees, in which genomic locations and dated segmental duplications were used to reconstruct gene phylogeny. In this article we review recent work concerning the patterns of expansin evolution in angiosperms and include additional insights gained from the genome of a second eudicot species, Populus trichocarpa, which includes at least 36 expansin genes. All of the previously proposed monocot-eudicot orthologous groups, but no additional ones, are represented in this species. The results also confirm that all of these clades are truly independent lineages. Furthermore, we have used position-based phylogeny to clarify the history of clades EXPA-II and EXPA-IV. Most of the growth of the expansin superfamily in the poplar lineage is likely due to a recent polyploidy event. Finally, some monocot-eudicot clades are shown to have diverged before the separation of the angiosperm and gymnosperm lineages. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

毛竹扩展蛋白基因的鉴定及其表达分析

为全面了解毛竹中扩展蛋白的分子特征和表达模式,本研究利用生物信息学方法在毛竹基因组中共鉴定出43个扩展蛋白基因家族成员,属于4个亚家族（EXPA、EXPB、EXLA和EXLB）,分别包含18、17、7和1个成员,分布在37个Scaffold上。除PeEXPA1没有内含子和PeEXLB1含有11个内含子外,其它毛竹扩展蛋白基因的内含子为1~5个。毛竹扩展蛋白基因编码蛋白长度为91~508个氨基酸,所有的氨基酸都具有高频密码子,大部分蛋白为碱性亲水性蛋白。大部分毛竹扩展蛋白二级结构中β转角占比例最少,而β折叠占比例最大,各亚家族多数成员具有类似的三级结构。qRT-PCR结果表明,18个EXPA亚家族成员在不同组织表达存在明显差异,除PeEXPA2、PeEXPA6外其它基因表达的最高值均出现在叶片中,表明它们可能在叶片生长过程中发挥着重要作用。     

Expansins: expanding importance in plant growth and development

Expansins were originally identified as cell wall-loosening proteins. The existence and various roles of expansins have been discovered in many plants. Expansins are encoded by a superfamily of genes comprised of subfamilies that evolved from a common ancestor and encode the α-expansins (EXPAs), the β-expansins (EXPBs), the expansin-like A (EXLA), and expansin-like B (EXLB) proteins. Several expansin-like genes have also been identified in non-plant organisms (e.g. a slime mold, fungi, nematodes, and a mollusk). Localization of EXPA and EXPB in the cell wall was confirmed by immunogold electron microscopy. Studies using transgenic plants provided evidence for a broad range of biological roles of expansins in diverse aspects of plant growth and development, such as cell wall extension, fruit softening, abscission, floral organ development, symbiosis, and the response to environmental stresses.     

Evolutionary divergence of β–expansin structure and function in grasses parallels emergence of distinctive primary cell wall traits

Expansins are wall‐loosening proteins that promote the extension of primary cell walls without the hydrolysis of major structural components. Previously, proteins from the EXPA (α–expansin) family were found to loosen eudicot cell walls but to be less effective on grass cell walls, whereas the reverse pattern was found for EXPB (β–expansin) proteins obtained from grass pollen. To understand the evolutionary and structural bases for the selectivity of EXPB action, we assessed the extension (creep) response of cell walls from diverse monocot families to EXPA and EXPB treatments. Cell walls from Cyperaceae and Juncaceae (families closely related to grasses) displayed a typical grass response (‘β–response’). Walls from more distant monocots, including some species that share with grasses high levels of arabinoxylan, responded preferentially to α–expansins (‘α–response’), behaving in this regard like eudicots. An expansin with selective activity for grass cell walls was detected in Cyperaceae pollen, coinciding with the expression of genes from the divergent EXPB–I branch that includes grass pollen β–expansins. The evolutionary origin of this branch was located within Poales on the basis of phylogenetic analyses and its association with the ‘sigma’ whole‐genome duplication. Accelerated evolution in this branch has remodeled the protein surface in contact with the substrate, potentially for binding highly substituted arabinoxylan. We propose that the evolution of the divergent EXPB–I group made a fundamental change in the target and mechanism of wall loosening in the grass lineage possible, involving a new structural role for xylans and the expansins that target them.     

小立碗藓扩展蛋白基因家族的鉴定与生物信息学分析

扩展蛋白(Expansins,EXP)是一类基因家族,几乎参与了植物发育的全过程,从种子萌发到果实成熟都有扩展蛋白的参与。该研究利用生物信息学的方法对小立碗藓(Physcomitrella patens) Expansin基因家族成员进行鉴定,分析了其基因结构、染色体定位以及系统发生关系。结果表明:小立碗藓基因组中含有Expansin A(EXPA) 32个、Expansin-like A(EXLA) 6个,并未发现Expansin-like B(EXLB)及Expansin B(EXPB)。扩展蛋白氨基酸序列长度在228～290 aa之间,编码蛋白质具有两个保守的结构域Pollen_allerg_1和DPBB_1。蛋白质亚细胞定位预测结果表明:运用CELLO在线工具预测发现小立碗藓中约4/5的EXP家族基因定位于细胞外;而Euk-mPLoc预测结果则显示小立碗藓EXP基因家族成员全定位于细胞外。基因结构分析表明,小立碗藓中约68%Expansin基因有含有1～3个内含子。以上结果可为深入研究小立碗藓扩展蛋白基因的分子进化与生物学功能奠定基础。     

玉米扩展蛋白Expansin基因家族定位及基因表达模式分析

扩展蛋白（Expansin）是一类能够使植物细胞壁松弛的活性蛋白,在植物生长发育过程中起着重要作用。利用PLAZA、NCBI、MaizeGDB、Uniprot、PLEXdb等基因组数据库,获得玉米Expansin家族的基因序列、染色体基因座位、蛋白质序列以及长度,构建玉米Expansin基因家族系统进化树,进行基因组织表达谱的分析。结果表明,玉米基因组中含有93个Expansin基因,分布于玉米的9条染色体上;多数Expansin具有250~300个氨基酸;玉米Expansin基因家族有40个Expansin A（EXPA）、47个Expansin B（EXPB）、6个Expansin-like A（EXLA）,未发现Expansin-like B（EXLB）;44个玉米Expansin基因在不同玉米组织中特异表达。该研究结果不仅为玉米扩展蛋白的深入研究奠定了基础,而且为其他研究人员对基因信息的获取提供了参考。     

The mitochondrial genome of the moss Physcomitrella patens sheds new light on mitochondrial evolution in land plants

The phylogenetic positions of bryophytes and charophytes, together with their genome features, are important for understanding early land plant evolution. Here we report the complete nucleotide sequence (105,340 bp) of the circular-mapping mitochondrial DNA of the moss Physcomitrella patens. Available evidence suggests that the multipartite structure of the mitochondrial genome in flowering plants does not occur in Physcomitrella. It contains genes for 3 rRNAs (rnl, rns, and rrn5), 24 tRNAs, and 42 conserved mitochondrial proteins (14 ribosomal proteins, 4 ccm proteins, 9 nicotinamide adenine dinucleotide dehydrogenase subunits, 5 ATPase subunits, 2 succinate dehydrogenase subunits, apocytochrome b, 3 cytochrome oxidase subunits, and 4 other proteins). We estimate that 5 tRNA genes are missing that might be encoded by the nuclear genome. The overall mitochondrial genome structure is similar in Physcomitrella, Chara vulgaris, Chaetosphaeridium globosum, and Marchantia polymorpha, with easily identifiable inversions and translocations. Significant synteny with angiosperm and chlorophyte mitochondrial genomes was not detected. Phylogenetic analysis of 18 conserved proteins suggests that the moss-liverwort clade is sister to angiosperms, which is consistent with a previous analysis of chloroplast genes but is not consistent with some analyses using mitochondrial sequences. In Physcomitrella, 27 introns are present within 16 genes. Nine of its intron positions are shared with angiosperms and 4 with Marchantia, which in turn shares only one intron position with angiosperms. The phylogenetic analysis as well as the syntenic structure suggest that the mitochondrial genomes of Physcomitrella and Marchantia retain prototype features among land plant mitochondrial genomes.     

Phylogenomic analysis of gene co‐expression networks reveals the evolution of functional modules

Molecular evolutionary studies correlate genomic and phylogenetic information with the emergence of new traits of organisms. These traits are, however, the consequence of dynamic gene networks composed of functional modules, which might not be captured by genomic analyses. Here, we established a method that combines large‐scale genomic and phylogenetic data with gene co‐expression networks to extensively study the evolutionary make‐up of modules in the moss Physcomitrella patens, and in the angiosperms Arabidopsis thaliana and Oryza sativa (rice). We first show that younger genes are less annotated than older genes. By mapping genomic data onto the co‐expression networks, we found that genes from the same evolutionary period tend to be connected, whereas old and young genes tend to be disconnected. Consequently, the analysis revealed modules that emerged at a specific time in plant evolution. To uncover the evolutionary relationships of the modules that are conserved across the plant kingdom, we added phylogenetic information that revealed duplication and speciation events on the module level. This combined analysis revealed an independent duplication of cell wall modules in bryophytes and angiosperms, suggesting a parallel evolution of cell wall pathways in land plants. We provide an online tool allowing plant researchers to perform these analyses at http://www.gene2function.de .     

Conservation and diversification of HAIRY MERISTEM gene family in land plants

The shoot apical meristems (SAMs) of land plants are crucial for plant growth and organ formation. In several angiosperms, the HAIRY MERISTEM (HAM) genes function as key regulators that control meristem development and stem cell homeostasis. To date, the origin and evolutionary history of the HAM family in land plants remains unclear. Potentially shared and divergent functions of HAM family members from angiosperms and non-angiosperms are also not known. In constructing a comprehensive phylogeny of the HAM family, we show that HAM proteins are widely present in land plants and that HAM proteins originated prior to the divergence of bryophytes. The HAM family was duplicated in a common ancestor of angiosperms, leading to two distinct groups: type I and type II. Type-II HAM members are widely present in angiosperms, whereas type-I HAM members were independently lost in different orders of monocots. Furthermore, HAM members from angiosperms and non-angiosperms (including bryophytes, lycophytes, ferns and gymnosperms) are able to replace the role of the type-II HAM genes in Arabidopsis, maintaining established SAMs and promoting the initiation of new stem cell niches. Our results uncover the conserved functions of HAM family members and reveal the conserved regulatory mechanisms underlying HAM expression patterning in meristems, providing insight into the evolution of key stem cell regulators in land plants.     

The mitochondrial genome of the gymnosperm Cycas taitungensis contains a novel family of short interspersed elements, Bpu sequences, and abundant RNA editing sites

The mtDNA of Cycas taitungensis is a circular molecule of 414,903 bp, making it 2- to 6-fold larger than the known mtDNAs of charophytes and bryophytes, but similar to the average of 7 elucidated angiosperm mtDNAs. It is characterized by abundant RNA editing sites (1,084), more than twice the number found in the angiosperm mtDNAs. The A + T content of Cycas mtDNA is 53.1%, the lowest among known land plants. About 5% of the Cycas mtDNA is composed of a novel family of mobile elements, which we designated as "Bpu sequences." They share a consensus sequence of 36 bp with 2 terminal direct repeats (AAGG) and a recognition site for the Bpu 10I restriction endonuclease (CCTGAAGC). Comparison of the Cycas mtDNA with other plant mtDNAs revealed many new insights into the biology and evolution of land plant mtDNAs. For example, the noncoding sequences in mtDNAs have drastically expanded as land plants have evolved, with abrupt increases appearing in the bryophytes, and then in the seed plants. As a result, the genomic organizations of seed plant mtDNAs are much less compact than in other plants. Also, the Cycas mtDNA appears to have been exempted from the frequent gene loss observed in angiosperm mtDNAs. Similar to the angiosperms, the 3 Cycas genes nad1, nad2, and nad5 are disrupted by 5 group II intron squences, which have brought the genes into trans-splicing arrangements. The evolutionary origin and invasion/duplication mechanism of the Bpu sequences in Cycas mtDNA are hypothesized and discussed.     

Phylogenetic analysis of the plant-specific zinc finger-homeobox and mini zinc finger gene families

Zinc finger-homeodomain proteins (ZHD) are present in many plants; however, the evolutionary history of the ZHD gene family remains largely unknown. We show here that ZHD genes are plant-specific, nearly all intronless, and related to MINI ZINC FINGER ( MIF ) genes that possess only the zinc finger. Phylogenetic analyses of ZHD genes from representative land plants suggest that non-seed plant ZHD genes occupy basal positions and angiosperm homologs form seven distinct clades. Several clades contain genes from two or more major angiosperm groups, including eudicots, monocots, magnoliids, and other basal angiosperms, indicating that several duplications occurred before the diversification of flowering plants. In addition, specific lineages have experienced more recent duplications. Unlike the ZHD genes, MIF s are found only from seed plants, possibly derived from ZHD s by loss of the homeodomain before the divergence of seed plants. Moreover, the MIF genes have also undergone relatively recent gene duplications. Finally, genome duplication might have contributed substantially to the expansion of family size in angiosperms and caused a high level of functional redundancy/overlap in these genes.     

Species-specific size expansion and molecular evolution of the oleosins in angiosperms

Oleosins are hydrophobic plant proteins thought to be important for the formation of oil bodies, which supply energy for seed germination and subsequent seedling growth. To better understand the evolutionary history and diversity of the oleosin gene family in plants, especially angiosperms, we systematically investigated the molecular evolution of this family using eight representative angiosperm species. A total of 73 oleosin members were identified, with six members in each of four monocot species and a greater but variable number in the four eudicots. A phylogenetic analysis revealed that the angiosperm oleosin genes belonged to three monophyletic lineages. Species-specific gene duplications, caused mainly by segmental duplication, led to the great expansion of oleosin genes and occurred frequently in eudicots after the monocot–eudicot divergence. Functional divergence analyses indicate that significant amino acid site-specific selective constraints acted on the different clades of oleosins. Adaptive evolution analyses demonstrate that oleosin genes were subject to strong purifying selection after their species-specific duplications and that rapid evolution occurred with a high degree of evolutionary dynamics in the pollen-specific oleosin genes. In conclusion, this study serves as a foundation for genome-wide analyses of the oleosins. These findings provide insight into the function and evolution of this gene family in angiosperms and pave the way for studies in other plants.