Algebraic connectivity may explain the evolution of gene regulatory networks

Gene expression is a result of the interplay between the structure, type, kinetics, and specificity of gene regulatory interactions, whose diversity gives rise to the variety of life forms. As the dynamic behavior of gene regulatory networks depends on their structure, here we attempt to determine structural reasons which, despite the similarities in global network properties, may explain the large differences in organismal complexity. We demonstrate that the algebraic connectivity, the smallest non-trivial eigenvalue of the Laplacian, of the directed gene regulatory networks decreases with the increase of organismal complexity, and may therefore explain the difference between the variety of analyzed regulatory networks. In addition, our results point out that, for the species considered in this study, evolution favours decreasing concentration of strategically positioned feed forward loops, so that the network as a whole can increase the specificity towards changing environments. Moreover, contrary to the existing results, we show that the average degree, the length of the longest cascade, and the average cascade length of gene regulatory networks cannot recover the evolutionary relationships between organisms. Whereas the dynamical properties of special subnetworks are relatively well understood, there is still limited knowledge about the evolutionary reasons for the already identified design principles pertaining to these special subnetworks, underlying the global quantitative features of gene regulatory networks of different organisms. The behavior of the algebraic connectivity, which we show valid on gene regulatory networks extracted from curated databases, can serve as an additional evolutionary principle of organism-specific regulatory networks.     

Rethinking convergence in plant parasitism through the lens of molecular and population genetic processes

The autotrophic lifestyle of photosynthetic plants has profoundly shaped their body plan, physiology, and gene repertoire. Shifts to parasitism and heterotrophy have evolved at least 12 times in more than 4000 species, and this transition has consequently left major evolutionary footprints among these parasitic lineages. Features that are otherwise rare at the molecular level and beyond have evolved repetitively, including reduced vegetative bodies, carrion-mimicking during reproduction, and the incorporation of alien genetic material. Here, I propose an integrated conceptual model, referred to as the funnel model, to define the general evolutionary trajectory of parasitic plants and provide a mechanistic explanation for their convergent evolution. This model connects our empirical understanding of gene regulatory networks in flowering plants with classical theories of molecular and population genetics. It emphasizes that the cascading effects brought about by the loss of photosynthesis may be a major force constraining the physiological capacity of parasitic plants and shaping their genomic landscapes. Here I review recent studies on the anatomy, physiology, and genetics of parasitic plants that lend support to this photosynthesis-centered funnel model. Focusing on nonphotosynthetic holoparasites, I elucidate how they may inevitably reach an evolutionary terminal status (i.e., extinction) and highlight the utility of a general, explicitly described and falsifiable model for future studies of parasitic plants.     

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 .     

De Novo Evolution of Complex, Global and Hierarchical Gene Regulatory Mechanisms

Gene regulatory networks exhibit complex, hierarchical features such as global regulation and network motifs. There is much debate about whether the evolutionary origins of such features are the results of adaptation, or the by-products of non-adaptive processes of DNA replication. The lack of availability of gene regulatory networks of ancestor species on evolutionary timescales makes this a particularly difficult problem to resolve. Digital organisms, however, can be used to provide a complete evolutionary record of lineages. We use a biologically realistic evolutionary model that includes gene expression, regulation, metabolism and biosynthesis, to investigate the evolution of complex function in gene regulatory networks. We discover that: (i) network architecture and complexity evolve in response to environmental complexity, (ii) global gene regulation is selected for in complex environments, (iii) complex, inter-connected, hierarchical structures evolve in stages, with energy regulation preceding stress responses, and stress responses preceding growth rate adaptations and (iv) robustness of evolved models to mutations depends on hierarchical level: energy regulation and stress responses tend not to be robust to mutations, whereas growth rate adaptations are more robust and non-lethal when mutated. These results highlight the adaptive and incremental evolution of complex biological networks, and the value and potential of studying realistic in silico evolutionary systems as a way of understanding living systems.     

Genes under weaker stabilizing selection increase network evolvability and rapid regulatory adaptation to an environmental shift

Regulatory networks play a central role in the modulation of gene expression, the control of cellular differentiation, and the emergence of complex phenotypes. Regulatory networks could constrain or facilitate evolutionary adaptation in gene expression levels. Here, we model the adaptation of regulatory networks and gene expression levels to a shift in the environment that alters the optimal expression level of a single gene. Our analyses show signatures of natural selection on regulatory networks that both constrain and facilitate rapid evolution of gene expression level towards new optima. The analyses are interpreted from the standpoint of neutral expectations and illustrate the challenge to making inferences about network adaptation. Furthermore, we examine the consequence of variable stabilizing selection across genes on the strength and direction of interactions in regulatory networks and in their subsequent adaptation. We observe that directional selection on a highly constrained gene previously under strong stabilizing selection was more efficient when the gene was embedded within a network of partners under relaxed stabilizing selection pressure. The observation leads to the expectation that evolutionarily resilient regulatory networks will contain optimal ratios of genes whose expression is under weak and strong stabilizing selection. Altogether, our results suggest that the variable strengths of stabilizing selection across genes within regulatory networks might itself contribute to the long‐term adaptation of complex phenotypes.     

Gene networks controlling the initiation of flower development

The onset of flower formation is a key regulatory event during the life cycle of angiosperm plants, which marks the beginning of the reproductive phase of development. It has been shown that floral initiation is under tight genetic control, and deciphering the underlying molecular mechanisms has been a main area of interest in plant biology for the past two decades. Here, we provide an overview of the developmental and genetic processes that occur during floral initiation. We further review recent studies that have led to the genome-wide identification of target genes of key floral regulators and discuss how they have contributed to an in-depth understanding of the gene regulatory networks controlling early flower development. We focus especially on a master regulator of floral initiation in Arabidopsis thaliana APETALA1 (AP1), but also outline what is known about the AP1 network in other plant species and the evolutionary implications.     

Detecting phylogenetic signal in mutualistic interaction networks using a Markov process model

Ecological interaction networks, such as those describing the mutualistic interactions between plants and their pollinators or between plants and their frugivores, exhibit non‐random structural properties that cannot be explained by simple models of network formation. One factor affecting the formation and eventual structure of such a network is its evolutionary history. We argue that this, in many cases, is closely linked to the evolutionary histories of the species involved in the interactions. Indeed, empirical studies of interaction networks along with the phylogenies of the interacting species have demonstrated significant associations between phylogeny and network structure. To date, however, no generative model explaining the way in which the evolution of individual species affects the evolution of interaction networks has been proposed. We present a model describing the evolution of pairwise interactions as a branching Markov process, drawing on phylogenetic models of molecular evolution. Using knowledge of the phylogenies of the interacting species, our model yielded a significantly better fit to 21% of a set of plant–pollinator and plant–frugivore mutualistic networks. This highlights the importance, in a substantial minority of cases, of inheritance of interaction patterns without excluding the potential role of ecological novelties in forming the current network architecture. We suggest that our model can be used as a null model for controlling evolutionary signals when evaluating the role of other factors in shaping the emergence of ecological networks.     

Convergent occurrence of the developmental hourglass in plant and animal embryogenesis?

Background The remarkable similarity of animal embryos at particular stages of development led to the proposal of a developmental hourglass. In this model, early events in development are less conserved across species but lead to a highly conserved ‘phylotypic period’. Beyond this stage, the model suggests that development once again becomes less conserved, leading to the diversity of forms. Recent comparative studies of gene expression in animal groups have provided strong support for the hourglass model. How and why might such an hourglass pattern be generated? More importantly, how might early acting events in development evolve while still maintaining a later conserved stage?Scope The discovery that an hourglass pattern may also exist in the embryogenesis of plants provides comparative data that may help us explain this phenomenon. Whether the developmental hourglass occurs in plants, and what this means for our understanding of embryogenesis in plants and animals is discussed. Models by which conserved early-acting genes might change their functional role in the evolution of gene networks, how networks buffer these changes, and how that might constrain, or confer diversity, of the body plan are also discused.Conclusions Evidence of a morphological and molecular hourglass in plant and animal embryogenesis suggests convergent evolution. This convergence is likely due to developmental constraints imposed upon embryogenesis by the need to produce a viable embryo with an established body plan, controlled by the architecture of the underlying gene regulatory networks. As the body plan is largely laid down during the middle phases of embryo development in plants and animals, then it is perhaps not surprising this stage represents the narrow waist of the hourglass where the gene regulatory networks are the oldest and most robust and integrated, limiting species diversity and constraining morphological space.     

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.     

Exploring the plant transcriptome through phylogenetic profiling

Publicly available protein sequences represent only a small fraction of the full catalog of genes encoded by the genomes of different plants, such as green algae, mosses, gymnosperms, and angiosperms. By contrast, an enormous amount of expressed sequence tags (ESTs) exists for a wide variety of plant species, representing a substantial part of all transcribed plant genes. Integrating protein and EST sequences in comparative and evolutionary analyses is not straightforward because of the heterogeneous nature of both types of sequence data. By combining information from publicly available EST and protein sequences for 32 different plant species, we identified more than 250,000 plant proteins organized in more than 12,000 gene families. Approximately 60% of the proteins are absent from current sequence databases but provide important new information about plant gene families. Analysis of the distribution of gene families over different plant species through phylogenetic profiling reveals interesting insights into plant gene evolution, and identifies species- and lineage-specific gene families, orphan genes, and conserved core genes across the green plant lineage. We counted a similar number of approximately 9,500 gene families in monocotyledonous and eudicotyledonous plants and found strong evidence for the existence of at least 33,700 genes in rice (Oryza sativa). Interestingly, the larger number of genes in rice compared to Arabidopsis (Arabidopsis thaliana) can partially be explained by a larger amount of species-specific single-copy genes and species-specific gene families. In addition, a majority of large gene families, typically containing more than 50 genes, are bigger in rice than Arabidopsis, whereas the opposite seems true for small gene families.     

种子植物抗细胞凋亡DAD基因的演化

抗细胞凋亡基因(DAD)是一个高度保守的细胞凋亡抑制基因,在植物生长发育中承担重要功能。为全面了解DAD基因在种子植物中的分布和演化规律,该文利用31种植物的全基因组数据,通过生物信息学手段,深入探讨和分析了不同植物类群中DAD基因的拷贝数目、基因结构和染色体定位,并综合另外7种裸子植物的转录组数据探讨了其在种子植物中的演化趋势。结果表明,DAD基因属于低拷贝基因,在不同种子植物中只具有1–3个拷贝;不同DAD基因编码的氨基酸长度在108–170 aa之间变动。同线性和系统发育分析进一步表明,种子植物DAD基因的演化具有明显的谱系特异性。随机复制和染色体大片段复制及其随后的基因丢失可能是其维持低拷贝的重要方式。     

Evolution of plant Ash1 SET genes: structural divergence and functional differentiation

Plant Ash1 SET proteins are involved in H3K36 methylation, and play a key role in plant reproductive development. Genes encoding Ash1 SET proteins constitute a multigene family in which the copy number varies among plant species and functional divergence appears to have occurred repeatedly. To investigate the evolutionary history and functional differentiation of the Ash1 SET gene family, we made a comprehensive evolutionary analysis of this gene family from eleven major representatives of green plants. A novel deep sister relationship grouping previously resolved II-1 and II-2 orthologous groups was identified. The absence of AWS domain in the group II-2 suggests that the independent losses of AWS domain have occurred during evolution. A diversity of gene structures in plant Ash1 SET gene family have been presented since the divergence of Physcomitrella patens (moss) from the other land plants. A small proportion of codons in SET domain regions were detected to be under positive selection along the branches ancestral to land plant and angiosperms, which may have allowed changes of substrate specificity among different evolutionary groups while maintaining the primary function of SET domains. Our predictive subcellular localization and comparative anatomical meta-expression analyses can assort with the structural divergences of Ash1 SET proteins.     

ATTED-II updates: condition-specific gene coexpression to extend coexpression analyses and applications to a broad range of flowering plants

ATTED-II (http://atted.jp) is a gene coexpression database for a wide variety of experimental designs, such as prioritizations of genes for functional identification and analyses of the regulatory relationships among genes. Here, we report updates of ATTED-II focusing on two new features: condition-specific coexpression and homologous coexpression with rice. To analyze a broad range of biological phenomena, it is important to collect data under many diverse experimental conditions, but the meaning of coexpression can become ambiguous under these conditions. One approach to overcome this difficulty is to calculate the coexpression for each set of conditions with a clear biological meaning. With this viewpoint, we prepared five sets of experimental conditions (tissue, abiotic stress, biotic stress, hormones and light conditions), and users can evaluate the coexpression by employing comparative gene lists and switchable gene networks. We also developed an interactive visualization system, using the Cytoscape web system, to improve the network representation. As the second update, rice coexpression is now available. The previous version of ATTED-II was specifically developed for Arabidopsis, and thus coexpression analyses for other useful plants have been difficult. To solve this problem, we extended ATTED-II by including comparison tables between Arabidopsis and rice. This representation will make it possible to analyze the conservation of coexpression among flowering plants. With the ability to investigate condition-specific coexpression and species conservation, ATTED-II can help researchers to clarify the functional and regulatory networks of genes in a broad array of plant species.     

植物查尔酮合成酶超基因家族的分子进化

查尔酮合成酶(CHS)超基因家族又称为植物类型III聚酮合酶超基因家族, 其编码酶通过催化和合成一系列结构多样及生理活性各异的次生代谢物, 在植物生长发育和适应环境的过程中扮演着重要角色。为全面了解CHS超基因家族在植物中的进化规律, 重建其进化历史, 该研究利用14种具有全基因组数据的代表植物, 通过生物信息学手段, 深入挖掘和分析了不同植物类群基因组中查尔酮合成酶超基因家族的成员构成, 推测了其可能的扩增机制和功能分歧, 并探讨了该超基因家族在植物中的总体进化趋势。结果共识别144条具有表达信息的同源序列, 它们全部来自9种陆生植物的基因组, 藻类植物基因组中没有发现相关序列。系统发育和进化分析表明, CHS超基因家族的起源古老, 它们可能为适应复杂的生态环境而出现在早期的陆生植物中, 之后在长期的进化过程中不断发生谱系的特异扩张和拷贝丢失, 最后通过功能分歧的形式在不同植物类群中被分别固定。此外, 进化检验也显示, 尽管CHS超基因家族内部发生了多样的遗传改变, 但整个超基因家族仍处于强烈的纯化选择之下, 并且个体基因中也无任何单氨基酸位点受到正向选择的影响。