Gene networks and the evolution of plant morphology

Herbarium specimens reveal links between Capsella bursa-pastoris leaf shape and climate

Studies into the evolution and development of leaf shape have connected variation in plant form, function, and fitness. For species with consistent leaf margin features, patterns in leaf architecture are related to both biotic and abiotic factors. However, for species with inconsistent leaf margin features, quantifying leaf shape variation and the effects of environmental factors on leaf shape has proven challenging. To investigate leaf shape variation in species with inconsistent shapes, we analyzed approximately 500 digitized Capsella bursa-pastoris specimens collected throughout the continental U.S. over a 100-year period with geometric morphometric modeling and deterministic techniques. We generated a morphospace of C. bursa-pastoris leaf shapes and modeled leaf shape as a function of environment and time. Our results suggest C. bursa-pastoris leaf shape variation is strongly associated with temperature over the C. bursa-pastoris growing season, with lobing decreasing as temperature increases. While we expected to see changes in variation over time, our results show that level of leaf shape variation is consistent over the 100-year period. Our findings showed that species with inconsistent leaf shape variation can be quantified using geometric morphometric modeling techniques and that temperature is the main environmental factor influencing leaf shape variation.

How to utilize comparative transcriptomics to dissect morphological diversity in plants

Comparative transcriptomics has emerged as a powerful approach that allows us to unravel the genetic basis of organ morphogenesis and its diversification processes during evolution. However, the application of comparative transcriptomics in studying plant morphological diversity addresses challenges such as identifying homologous gene pairs, selecting appropriate developmental stages for comparison, and extracting biologically meaningful networks. Methods such as phylostratigraphy, clustering, and gene co-expression networks are explored to identify functionally equivalent genes, align developmental stages, and uncover gene regulatory relationships. In the current review, we highlight the importance of these approaches in overcoming the complexity of plant genomes, the impact of heterochrony on stage alignment, and the integration of gene networks with additional data for a comprehensive understanding of morphological evolution.

Genome-wide identification of bZIP gene family and expression analysis of BhbZIP58 under heat stress in wax gourd

BACKGROUND

The basic leucine zipper (bZIP) transcription factor family is one of the most abundant and evolutionarily conserved gene families in plants. It assumes crucial functions in the life cycle of plants, including pathogen defense, secondary metabolism, stress response, seed maturation, and flower development. Although the genome of wax gourd has been published, little is known about the functions, evolutionary background, and gene expression patterns of the bZIP gene family, which limits its utilization.

RESULTS

A total of 61 bZIP genes (BhbZIPs) were identified from wax gourd (Benincasa hispida) genome and divided into 12 subgroups. Whole-genome duplication (WGD) and dispersed duplication (DSD) were the main driving forces of bZIP gene family expansion in wax gourd, and this family may have undergone intense purifying selection pressure during the evolutionary process. We selected BhbZIP58, only one in the member of subgroup B, to study its expression patterns under different stresses, including heat, salt, drought, cold stress, and ABA treatment. Surprisingly, BhbZIP58 had a dramatic response under heat stress. BhbZIP58 showed the highest expression level in the root compared with leaves, stem, stamen, pistil, and ovary. In addition, BhbZIP58 protein was located in the nucleus and had transcriptional activation activity. Overexpression of BhbZIP58 in Arabidopsis enhanced their heat tolerance.

CONCLUSIONS

In this study, bZIP gene family is systematically bioinformatically in wax gourd for the first time. Particularly, BhbZIP58 may have an important role in heat stress. It will facilitate further research on the bZIP gene family regarding their evolutionary history and biological functions.

The role of Class Ⅱ KNOX family in controlling compound leaf patterning in Medicago truncatula

Compound leaf development requires the coordination of genetic factors, hormones, and other signals. In this study, we explored the functions of Class Ⅱ KNOTTED-like homeobox (KNOXII) genes in the model leguminous plant Medicago truncatula. Phenotypic and genetic analyses suggest that MtKNOX4, 5 are able to repress leaflet formation, while MtKNOX3, 9, 10 are not involved in this developmental process. Further investigations have shown that MtKNOX4 represses the CK signal transduction, which is downstream of MtKNOXⅠ-mediated CK biosynthesis. Additionally, two boundary genes, FUSED COMPOUND LEAF1 (orthologue of Arabidopsis Class M KNOX) and NO APICAL MERISTEM (orthologue of Arabidopsis CUP-SHAPED COTYLEDON), are necessary for MtKNOX4-mediated compound leaf formation. These findings suggest, that among the members of MtKNOXⅡ, MtKNOX4 plays a crucial role in integrating the CK pathway and boundary regulators, providing new insights into the roles of MtKNOXⅡ in regulating the elaboration of compound leaves in M. truncatula.

Integration of eQTL and machine learning to dissect causal genes with pleiotropic effects in genetic regulation networks of seed cotton yield

The dissection of a gene regulatory network (GRN) that complements the genome-wide association study (GWAS) locus and the crosstalk underlying multiple agronomical traits remains a major challenge. In this study, we generate 558 transcriptional profiles of lint-bearing ovules at one day post-anthesis from a selective core cotton germplasm, from which 12,207 expression quantitative trait loci (eQTLs) are identified. Sixty-six known phenotypic GWAS loci are colocalized with 1,090 eQTLs, forming 38 functional GRNs associated predominantly with seed yield. Of the eGenes, 34 exhibit pleiotropic effects. Combining the eQTLs within the seed yield GRNs significantly increases the portion of narrow-sense heritability. The extreme gradient boosting (XGBoost) machine learning approach is applied to predict seed cotton yield phenotypes on the basis of gene expression. Top-ranking eGenes (NF-YB3, FLA2, and GRDP1) derived with pleiotropic effects on yield traits are validated, along with their potential roles by correlation analysis, domestication selection analysis, and transgenic plants.

Fundamentals of Plant Morphology and Plant Evo-Devo (Evolutionary Developmental Morphology)

Morphological concepts are used in plant evo-devo (evolutionary developmental biology) and other disciplines of plant biology, and therefore plant morphology is relevant to all of these disciplines. Many plant biologists still rely on classical morphology, according to which there are only three mutually exclusive organ categories in vascular plants such as flowering plants: root, stem (caulome), and leaf (phyllome). Continuum morphology recognizes a continuum between these organ categories. Instead of Aristotelian identity and either/or logic, it is based on fuzzy logic, according to which membership in a category is a matter of degree. Hence, an organ in flowering plants may be a root, stem, or leaf to some degree. Homology then also becomes a matter of degree. Process morphology supersedes structure/process dualism. Hence, structures do not have processes, they are processes, which means they are process combinations. These process combinations may change during ontogeny and phylogeny. Although classical morphology on the one hand and continuum and process morphology on the other use different kinds of logic, they can be considered complementary and thus together they present a more inclusive picture of the diversity of plant form than any one of the three alone. However, continuum and process morphology are more comprehensive than classical morphology. Insights gained from continuum and process morphology can inspire research in plant morphology and plant evo-devo, especially MorphoEvoDevo.

Flower-like meristem conditions and spatial constraints shape architecture of floral pseudanthia in Apioideae

BACKGROUND

Pseudanthia are multiflowered units that resemble single flowers, frequently by association with pseudocorollas formed by enlarged peripheral florets (ray flowers). Such resemblance is not only superficial, because numerous pseudanthia originate from peculiar reproductive meristems with flower-like characteristics, i.e. floral unit meristems (FUMs). Complex FUM-derived pseudanthia with ray flowers are especially common in Apiaceae, but our knowledge about their patterning is limited. In this paper, we aimed to investigate both the genetic and morphological basis of their development.

RESULTS

We analysed umbel morphogenesis with SEM in six species representing four clades of Apiaceae subfamily Apioideae with independently acquired floral pseudanthia. Additionally, using in situ hybridization, we investigated expression patterns of LEAFY (LFY), UNUSUAL FLORAL ORGANS (UFO), and CYCLOIDEA (CYC) during umbel development in carrot (Daucus carota subsp. carota). Here, we show that initial differences in size and shape of umbel meristems influence the position of ray flower formation, whereas an interplay between peripheral promotion and spatial constraints in umbellet meristems take part in the establishment of specific patterns of zygomorphy in ray flowers of Apiaceae. This space-dependent patterning results from flower-like morphogenetic traits of the umbel which are also visible at the molecular level. Transcripts of DcLFY are uniformly distributed in the incipient umbel, umbellet and flower meristems, while DcCYC shows divergent expression in central and peripheral florets.

CONCLUSIONS

Our results indicate that umbels develop from determinate reproductive meristems with flower-like characteristics, which supports their recognition as floral units. The great architectural diversity and complexity of pseudanthia in Apiaceae can be explained by the unique conditions of FUMs-an interplay between expression of regulatory genes, specific spatio-temporal ontogenetic constraints and morphogenetic gradients arising during expansion and repetitive fractionation. Alongside Asteraceae, umbellifers constitute an interesting model for investigation of patterning in complex pseudanthia.

Differential growth dynamics control aerial organ geometry

How gene activities and biomechanics together direct organ shapes is poorly understood. Plant leaf and floral organs develop from highly similar initial structures and share similar gene expression patterns, yet they gain drastically different shapes later-flat and bilateral leaf primordia and radially symmetric floral primordia, respectively. We analyzed cellular growth patterns and gene expression in young leaves and flowers of Arabidopsis thaliana and found significant differences in cell growth rates, which correlate with convergence sites of phytohormone auxin that require polar auxin transport. In leaf primordia, the PRESSED-FLOWER-expressing middle domain grows faster than adjacent adaxial domain and coincides with auxin convergence. In contrast, in floral primordia, the LEAFY-expressing domain shows accelerated growth rates and pronounced auxin convergence. This distinct cell growth dynamics between leaf and flower requires changes in levels of cell-wall pectin de-methyl-esterification and mechanical properties of the cell wall. Data-driven computer model simulations at organ and cellular levels demonstrate that growth differences are central to obtaining distinct organ shape, corroborating in planta observations. Together, our study provides a mechanistic basis for the establishment of early aerial organ symmetries through local modulation of differential growth patterns with auxin and biomechanics.

Co-option of a carotenoid cleavage dioxygenase gene (CCD4a) into the floral symmetry gene regulatory network contributes to the polymorphic floral shape-color combinations in Chrysanthemum sensu lato

Morphological novelties, including formation of trait combinations, may result from de novo gene origination and/or co-option of existing genes into other developmental contexts. A variety of shape-color combinations of capitular florets occur in Chrysanthemum and its allies. We hypothesized that co-option of a carotenoid cleavage dioxygenase gene into the floral symmetry gene network would generate a white zygomorphic ray floret. We tested this hypothesis in an evolutionary context using species in Chrysanthemum sensu lato, a monophyletic group with diverse floral shape-color combinations, based on morphological investigation, interspecific crossing, molecular interaction and transgenic experiments. Our results showed that white color was significantly associated with floret zygomorphy. Specific expression of the carotenoid cleavage dioxygenase gene CCD4a in marginal florets resulted in white color. Crossing experiments between Chrysanthemum lavandulifolium and Ajania pacifica indicated that expression of CCD4a is trans-regulated. The floral symmetry regulator CYC2g can activate expression of CCD4a with a dependence on TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING (TCP) binding element 8 on the CCD4a promoter. Based on all experimental findings, we propose that gene co-option of carotenoid degradation into floral symmetry regulation, and the subsequent dysfunction or loss of either CYC2g or CCD4a, may have led to evolution of capitular shape-color patterning in Chrysanthemum sensu lato.

A seed-like proteome in oil-rich tubers

There are numerous examples of plant organs or developmental stages that are desiccation-tolerant and can withstand extended periods of severe water loss. One prime example are seeds and pollen of many spermatophytes. However, in some plants, also vegetative organs can be desiccation-tolerant. One example are the tubers of yellow nutsedge (Cyperus esculentus), which also store large amounts of lipids similar to seeds. Interestingly, the closest known relative, purple nutsedge (Cyperus rotundus), generates tubers that do not accumulate oil and are not desiccation-tolerant. We generated nanoLC-MS/MS-based proteomes of yellow nutsedge in five replicates of four stages of tuber development and compared them to the proteomes of roots and leaves, yielding 2257 distinct protein groups. Our data reveal a striking upregulation of hallmark proteins of seeds in the tubers. A deeper comparison to the tuber proteome of the close relative purple nutsedge (C. rotundus) and a previously published proteome of Arabidopsis seeds and seedlings indicates that indeed a seed-like proteome was found in yellow but not purple nutsedge. This was further supported by an analysis of the proteome of a lipid droplet-enriched fraction of yellow nutsedge, which also displayed seed-like characteristics. One reason for the differences between the two nutsedge species might be the expression of certain transcription factors homologous to ABSCISIC ACID INSENSITIVE3, WRINKLED1, and LEAFY COTYLEDON1 that drive gene expression in Arabidopsis seed embryos.

Spatially and Temporally Distributed Complexity-A Refreshed Framework for the Study of GRN Evolution

Irrespective of the heuristic value of interpretations of developmental processes in terms of gene regulatory networks (GRNs), larger-angle views often suffer from: (i) an inadequate understanding of the relationship between genotype and phenotype; (ii) a predominantly zoocentric vision; and (iii) overconfidence in a putatively hierarchical organization of animal body plans. Here, we constructively criticize these assumptions. First, developmental biology is pervaded by adultocentrism, but development is not necessarily egg to adult. Second, during development, many unicells undergo transcriptomic profile transitions that are comparable to those recorded in pluricellular organisms; thus, their study should not be neglected from the GRN perspective. Third, the putatively hierarchical nature of the animal body is mirrored in the GRN logic, but in relating genotype to phenotype, independent assessments of the dynamics of the regulatory machinery and the animal’s architecture are required, better served by a combinatorial than by a hierarchical approach. The trade-offs between spatial and temporal aspects of regulation, as well as their evolutionary consequences, are also discussed. Multicellularity may derive from a unicell’s sequential phenotypes turned into different but coexisting, spatially arranged cell types. In turn, polyphenism may have been a crucial mechanism involved in the origin of complex life cycles.

Genomic and Meiotic Changes Accompanying Polyploidization

Hybridization and polyploidy have been considered as significant evolutionary forces in adaptation and speciation, especially among plants. Interspecific gene flow generates novel genetic variants adaptable to different environments, but it is also a gene introgression mechanism in crops to increase their agronomical yield. An estimate of 9% of interspecific hybridization has been reported although the frequency varies among taxa. Homoploid hybrid speciation is rare compared to allopolyploidy. Chromosome doubling after hybridization is the result of cellular defects produced mainly during meiosis. Unreduced gametes, which are formed at an average frequency of 2.52% across species, are the result of altered spindle organization or orientation, disturbed kinetochore functioning, abnormal cytokinesis, or loss of any meiotic division. Meiotic changes and their genetic basis, leading to the cytological diploidization of allopolyploids, are just beginning to be understood especially in wheat. However, the nature and mode of action of homoeologous recombination suppressor genes are poorly understood in other allopolyploids. The merger of two independent genomes causes a deep modification of their architecture, gene expression, and molecular interactions leading to the phenotype. We provide an overview of genomic changes and transcriptomic modifications that particularly occur at the early stages of allopolyploid formation.

Gene expression evolution in pattern-triggered immunity within Arabidopsis thaliana and across Brassicaceae species

Plants recognize surrounding microbes by sensing microbe-associated molecular patterns (MAMPs) to activate pattern-triggered immunity (PTI). Despite their significance for microbial control, the evolution of PTI responses remains largely uncharacterized. Here, by employing comparative transcriptomics of six Arabidopsis thaliana accessions and three additional Brassicaceae species to investigate PTI responses, we identified a set of genes that commonly respond to the MAMP flg22 and genes that exhibit species-specific expression signatures. Variation in flg22-triggered transcriptome responses across Brassicaceae species was incongruent with their phylogeny, while expression changes were strongly conserved within A. thaliana. We found the enrichment of WRKY transcription factor binding sites in the 5'-regulatory regions of conserved and species-specific responsive genes, linking the emergence of WRKY-binding sites with the evolution of gene expression patterns during PTI. Our findings advance our understanding of the evolution of the transcriptome during biotic stress.

An Ontology-Driven Personalized Faceted Search for Exploring Knowledge Bases of Capsicum

Capsicum is a genus of flowering plants in the Solanaceae family in which the members are well known to have a high economic value. The Capsicum fruits, which are popularly known as peppers or chili, have been widely used by people worldwide. It serves as a spice and raw material for many products such as sauce, food coloring, and medicine. For many years, scientists have studied this plant to optimize its production. A tremendous amount of knowledge has been obtained and shared, as reflected in multiple knowledge-based systems, databases, or information systems. An approach to knowledge-sharing is through the adoption of a common ontology to eliminate knowledge understanding discrepancy. Unfortunately, most of the knowledge-sharing solutions are intended for scientists who are familiar with the subject. On the other hand, there are groups of potential users that could benefit from such systems but have minimal knowledge of the subject. For these non-expert users, finding relevant information from a less familiar knowledge base would be daunting. More than that, users have various degrees of understanding of the available content in the knowledge base. This understanding discrepancy raises a personalization problem. In this paper, we introduce a solution to overcome this challenge. First, we developed an ontology to facilitate knowledge-sharing about Capsicum to non-expert users. Second, we developed a personalized faceted search algorithm that provides multiple structured ways to explore the knowledge base. The algorithm addresses the personalization problem by identifying the degree of understanding about the subject from each user. In this way, non-expert users could explore a knowledge base of Capsicum efficiently. Our solution characterized users into four groups. As a result, our faceted search algorithm defines four types of matching mechanisms, including three ranking mechanisms as the core of our solution. In order to evaluate the proposed method, we measured the predictability degree of produced list of facets. Our findings indicated that the proposed matching mechanisms could tolerate various query types, and a high degree of predictability can be achieved by combining multiple ranking mechanisms. Furthermore, it demonstrates that our approach has a high potential contribution to biodiversity science in general, where many knowledge-based systems have been developed with limited access to users outside of the domain.

Deploying Big Data to Crack the Genotype to Phenotype Code

Mechanistically connecting genotypes to phenotypes is a longstanding and central mission of biology. Deciphering these connections will unite questions and datasets across all scales from molecules to ecosystems. Although high-throughput sequencing has provided a rich platform on which to launch this effort, tools for deciphering mechanisms further along the genome to phenome pipeline remain limited. Machine learning approaches and other emerging computational tools hold the promise of augmenting human efforts to overcome these obstacles. This vision paper is the result of a Reintegrating Biology Workshop, bringing together the perspectives of integrative and comparative biologists to survey challenges and opportunities in cracking the genotype to phenotype code and thereby generating predictive frameworks across biological scales. Key recommendations include promoting the development of minimum "best practices" for the experimental design and collection of data; fostering sustained and long-term data repositories; promoting programs that recruit, train, and retain a diversity of talent; and providing funding to effectively support these highly cross-disciplinary efforts. We follow this discussion by highlighting a few specific transformative research opportunities that will be advanced by these efforts.

A crosstalk between auxin and brassinosteroid regulates leaf shape by modulating growth anisotropy

Leaf shape is highly variable within and among plant species, ranging from slender to oval shaped. This is largely determined by the proximodistal axis of growth. However, little is known about how proximal-distal growth is controlled to determine leaf shape. Here, we show that Arabidopsis leaf and sepal proximodistal growth is tuned by two phytohormones. Two class A AUXIN RESPONSE FACTORs (ARFs), ARF6 and ARF8, activate the transcription of DWARF4, which encodes a key brassinosteroid (BR) biosynthetic enzyme. At the cellular level, the phytohormones promote more directional cell expansion along the proximodistal axis, as well as final cell sizes. BRs promote the demethyl-esterification of cell wall pectins, leading to isotropic in-plane cell wall loosening. Notably, numerical simulation showed that isotropic cell wall loosening could lead to directional cell and organ growth along the proximodistal axis. Taken together, we show that auxin acts through biosynthesis of BRs to determine cell wall mechanics and directional cell growth to generate leaves of variable roundness.

Polygenic Adaptation: Integrating Population Genetics and Gene Regulatory Networks

The adaptation of populations to local environments often relies on the selection of optimal values for polygenic traits. Here, we first summarize the results obtained from different quantitative genetics and population genetics models, about the genetic architecture of polygenic traits and their response to directional selection. We then highlight the contribution of systems biology to the understanding of the molecular bases of polygenic traits and the evolution of gene regulatory networks involved in these traits. Finally, we discuss the need for a unifying framework merging the fields of population genetics, quantitative genetics and systems biology to better understand the molecular bases of polygenic traits adaptation.

Zinc Finger-Homeodomain Transcriptional Factors (ZF-HDs) in Wheat (Triticum aestivum L.): Identification, Evolution, Expression Analysis and Response to Abiotic Stresses

Zinc finger-homeodomain transcriptional factors (ZF-HDs), a kind of plant-specific transcription factor, play important roles in plant growth, development and various stress responses. In this study, the genome-wide analysis of the ZF-HD gene family was performed in wheat. A total of 37 TaZF-HD genes were identified in T. aestivum and classified into six groups. The results of a synteny analysis showed that gene replication events contributed to the expansion of the TaZF-HD gene family. The TaZF-HD paralogous gene pairs with similar chromosomal locations in different subgenomes had similar expression patterns. TaZF-HDs were highly induced under PEG (polyethylene glycol), NaCl and cold stress but not induced under heat stress. Gene ontology (GO) annotation and protein-protein interactions suggested that TaZF-HD proteins may participate in various biological processes of plants. These results increase our understanding of ZF-HD genes and provide robust candidate genes for future functional investigations aimed at crop improvement.

In Vitro Plant Regeneration in Conifers: The Role of WOX and KNOX Gene Families

Conifers are a group of woody plants with an enormous economic and ecological importance. Breeding programs are necessary to select superior varieties for planting, but they have many limitations due to the biological characteristics of conifers. Somatic embryogenesis (SE) and de novo organogenesis (DNO) from in vitro cultured tissues are two ways of plant mass propagation that help to overcome this problem. Although both processes are difficult to achieve in conifers, they offer advantages like a great efficiency, the possibilities to cryopreserve the embryogenic lines, and the ability of multiplying adult trees (the main bottleneck in conifer cloning) through DNO. Moreover, SE and DNO represent appropriate experimental systems to study the molecular bases of developmental processes in conifers such as embryogenesis and shoot apical meristem (SAM) establishment. Some of the key genes regulating these processes belong to the WOX and KNOX homeobox gene families, whose function has been widely described in Arabidopsis thaliana. The sequences and roles of these genes in conifers are similar to those found in angiosperms, but some particularities exist, like the presence of WOXX, a gene that putatively participates in the establishment of SAM in somatic embryos and plantlets of Pinus pinaster.

Characterization of the KNOTTED1-LIKE HOMEOBOX (KNOX) gene family in Pinus pinaster Ait

KNOTTED1-LIKE HOMEOBOX (KNOX) genes are a family of plant-specific homeobox transcription factors with important roles in plant development that have been classified into two subfamilies with differential expression domains and functions. Studies in angiosperms have shown that class I members are related to the maintenance of meristem homeostasis and leaf development, whereas class II members promote differentiation of tissues and organs. However, little is known about its diversification and function in gymnosperms. By combining PCR-based detection and transcriptome data analysis, we identified four class I and two class II KNOX genes in Pinus pinaster. Expression analyses showed that class I members were mainly expressed in meristematic regions and differentiating tissues, with practically no expression in lateral organs, whereas expression of class II members was restricted to lateral organs. Furthermore, overexpression of P. pinaster KNOX genes in Arabidopsis thaliana caused similar phenotypic effects to those described for their angiosperms counterparts. This is the first time to our knowledge that functional analyses of class II members are reported in a conifer species. These results suggest a high conservation of the KNOX gene family throughout seed plants, as the functional differentiation of both subfamilies observed in angiosperms might be partially conserved in gymnosperms.

Genomics of Evolutionary Novelty in Hybrids and Polyploids

It has long been recognized that hybridization and polyploidy are prominent processes in plant evolution. Although classically recognized as significant in speciation and adaptation, recognition of the importance of interspecific gene flow has dramatically increased during the genomics era, concomitant with an unending flood of empirical examples, with or without genome doubling. Interspecific gene flow is thus increasingly thought to lead to evolutionary innovation and diversification, via adaptive introgression, homoploid hybrid speciation and allopolyploid speciation. Less well understood, however, are the suite of genetic and genomic mechanisms set in motion by the merger of differentiated genomes, and the temporal scale over which recombinational complexity mediated by gene flow might be expressed and exposed to natural selection. We focus on these issues here, considering the types of molecular genetic and genomic processes that might be set in motion by the saltational event of genome merger between two diverged species, either with or without genome doubling, and how these various processes can contribute to novel phenotypes. Genetic mechanisms include the infusion of new alleles and the genesis of novel structural variation including translocations and inversions, homoeologous exchanges, transposable element mobilization and novel insertional effects, presence-absence variation and copy number variation. Polyploidy generates massive transcriptomic and regulatory alteration, presumably set in motion by disrupted stoichiometries of regulatory factors, small RNAs and other genome interactions that cascade from single-gene expression change up through entire networks of transformed regulatory modules. We highlight both these novel combinatorial possibilities and the range of temporal scales over which such complexity might be generated, and thus exposed to natural selection and drift.

Temporal flexibility of gene regulatory network underlies a novel wing pattern in flies

Developmental genes can be coopted to generate evolutionary novelties by changing their spatial regulation. However, developmental genes seldom act independently, but rather work in a gene regulatory network (GRN). How is it possible to recruit a single gene from a whole GRN? What are the properties that allow parallel cooptions of the same genes during evolution? Here, we show that a novel engrailed gene expression underlies a novel wing color pattern in flies. We show that cooption is facilitated 1) because of GRN flexibility over development and 2) because every single gene of the GRN has its own functional time window. We suggest these two temporal properties could explain why the same gene can be independently recruited several times during evolution.

Organisms have evolved endless morphological, physiological, and behavioral novel traits during the course of evolution. Novel traits were proposed to evolve mainly by orchestration of preexisting genes. Over the past two decades, biologists have shown that cooption of gene regulatory networks (GRNs) indeed underlies numerous evolutionary novelties. However, very little is known about the actual GRN properties that allow such redeployment. Here we have investigated the generation and evolution of the complex wing pattern of the fly Samoaia leonensis. We show that the transcription factor Engrailed is recruited independently from the other players of the anterior–posterior specification network to generate a new wing pattern. We argue that partial cooption is made possible because 1) the anterior–posterior specification GRN is flexible over time in the developing wing and 2) this flexibility results from the fact that every single gene of the GRN possesses its own functional time window. We propose that the temporal flexibility of a GRN is a general prerequisite for its possible cooption during the course of evolution.

Pleiotropic changes revealed by in situ recovery of the semi-dwarf gene sd1 in rice

The "Green Revolution" that dramatically reduced cultivar heights and sharply boosted rice production mid-century was achieved in large part through introgression of defective alleles of Semi-Dwarf 1 (SD1), which encodes a GA20ox oxidase involved in the final steps of the synthesis of bioactive gibberellin in rice. Here, we ask whether converting the defective sd1 version in a modern semi-dwarf cultivar back to wild-type SD1 in situ recovers ancestral plant traits, and more broadly, what it reveals about pleiotropic effects of this gene. We assess these effects of SD1 restoration in three independent recombinant lines recovered from F₂ progeny of a cross between 93-11 and PA64s. We then used RNA-seq to dissect gene network changes that accompanied SD1 restoration. We report that this in situ restoration of wild-type SD1 nearly doubles plant height, increases total grain yield per panicle, and elongates the second-leaf length. Comparison of expression profiles reveals changes in key nodes of the gibberellin pathway, such as OsKO1 and OsGA2ox3, and more broadly in genes related to metabolic networks, defense response, and catabolic processes. Two JA-induced genes, RIR1b and OsPR1b, are extremely down-regulated after SD1 restoration, suggesting that SD1 restoration alters the balance between GA and JA to plant growth, at the cost of degrading the defense response. This in situ approach at the SD1 locus also provides a model example that is applicable to other systems and will further understanding of gene networks underlying high-yield traits in crops.

Contrasted evolutionary trajectories of plant transcription factors

Because of their prominent roles in plant development, transcription factors (TF) play central roles as drivers of innovation in the evolution of the green lineage (viridiplantae). The advent of massive sequencing combined with comparative genetics/genomics allows a rigorous investigation of how TF families have contributed to plant diversification from charophyte algae to bryophytes to angiosperms. Here, we review recent progress on TF family reconstruction and the identification of distantly related TFs present throughout the evolutionary timeline from algae to angiosperms. These data provide examples of contrasting evolutionary trajectories of TF families and illustrate how conserved TFs adopt diverse roles over the course of evolution.

Taxonomic Notes on the 'Mahat' (Artocarpus lacucha and A. thailandicus, Moraceae) Species Complex in Thailand

‘Mahat’ is a well-known medicinal plant utilized in Thailand. The Thai name ‘Mahat’ has been used in many scientific articles for years. However, it is, unpredictably, a homonym of two scientific names in Flora of Thailand, i.e., A. lacucha and A. thailandicus. Additionally, both species are complex due to their high morphological variation. This causes difficulties in species identification especially when this Thai name is referred to as the scientific name for research publication, quality control of pharmaceutical raw materials, and registration of pharmaceutical products. In this study, we scrutinized the taxonomy of ‘Mahat’ by detailed examination of its morphology and distribution, including molecular and qualitative phytochemical studies. Leaf surfaces were inspected using scanning electron microscopy. The phylogeny of both species was studied using DNA sequences of nuclear and plastid regions. Chromatographic fingerprints, focusing on the major active compound oxyresveratrol, were identified using high-performance liquid chromatography. According to our current study, phylogenetic evidence showed that some samples of both species were clustered together in the same clade and phytochemical fingerprints were almost identical. These results are valuable data for taxonomic revision in the near future and reveal the possible utilization of A. thailandicus as a new material source of oxyresveratrol in the pharmaceutical industry.

Transcriptional Landscapes of Divergent Sporophyte Development in Two Mosses, Physcomitrium (Physcomitrella) patens and Funaria hygrometrica

Understanding the molecular basis of morphological shifts is a fundamental question of evolutionary biology. New morphologies may arise through the birth/death of genes (gene gain/loss) or by reutilizing existing gene sets. Yet, the relative contribution of these two processes to radical morphological shifts is still poorly understood. Here, we use the model system of two mosses, Funaria hygrometrica and Physcomitrium (Physcomitrella) patens, to investigate the molecular mechanisms underlying contrasting sporophyte architectures. We used comparative analysis of time-series expression data for four stages of sporophyte development in both species to address this question in detail. We found that large-scale differences in sporophytic architecture are mainly governed by orthologous (i.e., shared) genes frequently experiencing temporal gene expression shifts between the two species. While the absolute number of species-specific genes expressed during sporophyte development is somewhat smaller, we observed a significant increase of their proportion in preferentially sporophyte expressed genes, suggesting a fundamental role in the sporophyte phase. However, further functional studies are necessary to determine their contribution to diverging sporophyte morphologies. Our results add to the growing set of studies suggesting that radical changes in morphology may rely on the heterochronic expression of conserved regulators.

Evolution of Lineage-Specific Gene Networks Underlying the Considerable Fruit Shape Diversity in Persimmon

The shapes of plant organs reflect the evolution of each lineage and have been diversified according to lineage-specific adaptations to environment. Research on the molecular pathways responsible for organ shapes has traditionally been focused mainly on leaves or flowers. Thus, little is known about the pathways controlling fruit shapes, despite their diversity in some plant species. In this study, we analyzed oriental persimmon (Diospyros kaki), which exhibits considerable diversity in fruit shapes among cultivars, to elucidate the underlying molecular mechanism using transcriptomic data and quantitative evaluation. First, to filter the candidate genes associated with persimmon fruit shapes, the whole gene expression patterns obtained using mRNA-Seq analysis from 100 individuals, including a segregated population and various cultivars, were assessed to detect correlations with principal component scores for fruit shapes characterized with elliptic Fourier descriptors. Next, a gene co-expression network analysis with weighted gene co-expression network analysis (WGCNA) package revealed that class 1 KNOX family genes and SEEDSTICK function as integrators along with some phytohormone-related genes, to regulate the fruit shape diversity. On the other hand, the OVATE family genes also contribute to fruit shape diversity, of which pathway would be potentially shared with other plant species. Evolutionary aspects suggest that acquisition of a high lineage-specific and variable expression of class 1 KNOX gene, knotted-like homeobox of Arabidopsis thaliana 1 (KNAT1), in young fruit is important for establishing the persimmon-specific mechanism that determines fruit shape diversity.

Genome-wide analysis of MpBHLH12, a IIIf basic helix-loop-helix transcription factor of Marchantia polymorpha

The evolution of plants on land required adaptation to UV radiation and dry environments, and involved the appearance and/or rewiring of genetic connections, known as gene regulatory networks (GRNs), which consist of one or more transcription factors (TFs). The liverwort, Marchantia polymorpha, is a basal land plant, with a recently sequenced genome. The number of genes encoding basic helix-loop-helix (bHLH) family members is considerably higher in M. polymorpha than in charophyte green algae, suggesting the contribution of bHLH proteins to the evolution of GRNs associated with the adaptation of plants to land. Although an understanding of the evolutionary aspects of GRNs is fundamental for elucidating the mechanisms of environmental adaptation, the evolution of GRNs that led to land adaptation in plants remains poorly understood. In this study, we isolated a single gene encoding a IIIf bHLH TF from M. polymorpha, MpBHLH12. Transgenic M. polymorpha constitutively overexpressing MpBHLH12 showed smaller and fewer gemma cups than wild type, suggesting that MpBHLH12 is involved in the regulation of morphological development. Transcriptomic analysis of MpBHLH12 overexpressor (MpBHLH12ox) lines revealed an overlap with the GRN of MpMYB14, which regulates the biosynthesis of anthocyanins and phenolic compounds. However, MpBHLH12ox did not show anthocyanin accumulation. Results of the transient reporter assay suggest that MpBHLH12 could function in repression rather than activation. Our findings suggest that although the IIIf bHLH MpBHLH12 shows highest amino acid similarity with IIIf bHLH clade and is involved in developmental process and partly biosynthesis of phenolic compounds in M. polymorpha like Arabidopsis IIIf bHLH, the GRN involving MpBHLH12 would be distinct one from those of the IIIf bHLH TFs of seed plants.

Perspectives and applications of machine learning for evolutionary developmental biology

Evolutionary Developmental Biology (Evo-Devo) is an ever-expanding field that aims to understand how development was modulated by the evolutionary process. In this sense, "omic" studies emerged as a powerful ally to unravel the molecular mechanisms underlying development. In this scenario, bioinformatics tools become necessary to analyze the growing amount of information. Among computational approaches, machine learning stands out as a promising field to generate knowledge and trace new research perspectives for bioinformatics. In this review, we aim to expose the current advances of machine learning applied to evolution and development. We draw clear perspectives and argue how evolution impacted machine learning techniques.