Evolutionary developmental transcriptomics reveals a gene network module regulating interspecific diversity in plant leaf shape

Shaping leaves through TALE homeodomain transcription factors

The first TALE homeodomain transcription factor gene to be described in plants was maize knotted1 (kn1). Dominant mutations in kn1 disrupt leaf development, with abnormal knots of tissue forming in the leaf blade. kn1 was found to be expressed in the shoot meristem but not in a peripheral region that gives rise to leaves. Furthermore, KN1 and closely related proteins were excluded from initiating and developing leaves. These findings were a prelude to a large body of work wherein TALE homeodomain proteins have been identified as vital regulators of meristem homeostasis and organ development in plants. KN1 homologues are widely represented across land plant taxa. Thus, studying the regulation and mechanistic action of this gene class has allowed investigations into the evolution of diverse plant morphologies. This review will focus on the function of TALE homeodomain transcription factors in leaf development in eudicots. Here, we discuss how TALE homeodomain proteins contribute to a spectrum of leaf forms, from the simple leaves of Arabidopsis thaliana to the compound leaves of Cardamine hirsuta and species beyond the Brassicaceae.

Leaf growth - complex regulation of a seemingly simple process

Understanding the underlying mechanisms of plant development is crucial to successfully steer or manipulate plant growth in a targeted manner. Leaves, the primary sites of photosynthesis, are vital organs for many plant species, and leaf growth is controlled by a tight temporal and spatial regulatory network. In this review, we focus on the genetic networks governing leaf cell proliferation, one major contributor to final leaf size. First, we provide an overview of six regulator families of leaf growth in Arabidopsis: DA1, PEAPODs, KLU, GRFs, the SWI/SNF complexes, and DELLAs, together with their surrounding genetic networks. Next, we discuss their evolutionary conservation to highlight similarities and differences among species, because knowledge transfer between species remains a big challenge. Finally, we focus on the increase in knowledge of the interconnectedness between these genetic pathways, the function of the cell cycle machinery as their central convergence point, and other internal and environmental cues.

Ancient allopatry and ecological divergence act together to promote plant diversity in mountainous regions: evidence from comparative phylogeography of two genera in the Sino-Himalayan region

BACKGROUND

How geographical isolation and ecological divergence act together to promote plant diversity in mountainous regions remains largely unknown. In this study, we chose two genera comprising a small number of species distributed in the Sino-Himalayan region, Megacodon (Gentianaceae) and Beesia (Ranunculaceae), which both exhibit a fragmented distribution pattern and are found across a wide range of elevations. By summarizing their common patterns of speciation and/or divergence processes, we aim to understand how environmental changes accelerated lineage diversification in the Sino-Himalayan region through ancient allopatry and ecological divergence.

RESULTS

Using ddRAD-seq, chloroplast genome sequences, and specific molecular markers, we studied the phylogenetic relationships, population structure, and historical biogeography of Beesia and Megacodon. Both genera began to diverge from the late Miocene onwards, with ancient allopatry at lower elevations formed narrow-range species or relict populations. Mantel tests between genetic distance and climatic, elevational, or geographic distance revealed an isolation-by-distance pattern in Beesia and Megacodon stylophorus. Megacodon showed two clades occupying entirely different altitudinal ranges, whereas Beesia calthifolia exhibited a genetic divergence pattern along an elevation gradient. Furthermore, we conducted morphological measurements on Beesia calthifolia and found that different elevational groups had distinct leaf shapes.

CONCLUSIONS

The regional disjunctions of plant groups in the Sino-Himalayan region are drastic and closely related to several biogeographic boundaries. As a consequence of major geological and climate change, ecological divergence when different elevations are colonized often happens simultaneously within plant groups. Although habitat fragmentation and parapatric ecological divergence each spur speciation to different extents, a combined effect of these two factors is a common phenomenon in the Sino-Himalayan region.

Immunity priming uncouples the growth-defense trade-off in tomato

Plants have developed an array of mechanisms to protect themselves against pathogen invasion. The deployment of defense mechanisms is imperative for plant survival, but can come at the expense of plant growth, leading to the 'growth-defense trade-off' phenomenon. Following pathogen exposure, plants can develop resistance to further attack. This is known as induced resistance, or priming. Here, we investigated the growth-defense trade-off, examining how defense priming via systemic acquired resistance (SAR), or induced systemic resistance (ISR), affects tomato development and growth. We found that defense priming can promote, rather than inhibit, plant development, and that defense priming and growth trade-offs can be uncoupled. Cytokinin response was activated during induced resistance, and found to be required for the observed growth and disease resistance resulting from ISR activation. ISR was found to have a stronger effect than SAR on plant development. Our results suggest that growth promotion and induced resistance can be co-dependent, and that, in certain cases, defense priming can drive developmental processes and promote plant yield.

Diversity and genetic architecture of agro-morphological traits in a core collection of European traditional tomato

European traditional tomato varieties have been selected by farmers given their consistent performance and adaptation to local growing conditions. Here we developed a multipurpose core collection, comprising 226 accessions representative of the genotypic, phenotypic, and geographical diversity present in European traditional tomatoes, to investigate the basis of their phenotypic variation, gene×environment interactions, and stability for 33 agro-morphological traits. Comparison of the traditional varieties with a modern reference panel revealed that some traditional varieties displayed excellent agronomic performance and high trait stability, as good as or better than that of their modern counterparts. We conducted genome-wide association and genome-wide environment interaction studies and detected 141 quantitative trait loci (QTLs). Out of those, 47 QTLs were associated with the phenotype mean (meanQTLs), 41 with stability (stbQTLs), and 53 QTL-by-environment interactions (QTIs). Most QTLs displayed additive gene actions, with the exception of stbQTLs, which were mostly recessive and overdominant QTLs. Both common and specific loci controlled the phenotype mean and stability variation in traditional tomato; however, a larger proportion of specific QTLs was observed, indicating that the stability gene regulatory model is the predominant one. Developmental genes tended to map close to meanQTLs, while genes involved in stress response, hormone metabolism, and signalling were found within regions affecting stability. A total of 137 marker-trait associations for phenotypic means and stability were novel, and therefore our study enhances the understanding of the genetic basis of valuable agronomic traits and opens up a new avenue for an exploitation of the allelic diversity available within European traditional tomato germplasm.

Deciphering the Molecular Mechanisms of Autonomic Nervous System Neuron Induction through Integrative Bioinformatics Analysis

In vitro derivation of human neurons in the autonomic nervous system (ANS) is an important technology, given its regulatory roles in maintaining homeostasis in the human body. Although several induction protocols for autonomic lineages have been reported, the regulatory machinery remains largely undefined, primarily due to the absence of a comprehensive understanding of the molecular mechanism regulating human autonomic induction in vitro. In this study, our objective was to pinpoint key regulatory components using integrated bioinformatics analysis. A protein-protein interaction network construction for the proteins encoded by the differentially expressed genes from our RNA sequencing data, and conducting subsequent module analysis, we identified distinct gene clusters and hub genes involved in the induction of autonomic lineages. Moreover, we analyzed the impact of transcription factor (TF) activity on target gene expression, revealing enhanced autonomic TF activity that could lead to the induction of autonomic lineages. The accuracy of this bioinformatics analysis was corroborated by employing calcium imaging to observe specific responses to certain ANS agonists. This investigation offers novel insights into the regulatory machinery in the generation of neurons in the ANS, which would be valuable for further understanding and precise regulation of autonomic induction and differentiation.

Diversity of tomato leaf form provides novel insights into breeding

Tomato (Solanum lycopersicum L.) is cultivated widely globally. The crop exhibits tremendous morphological variations because of its long breeding history. Apart from the commercial tomato varieties, wild species and heirlooms are grown in certain regions of the world. Since the fruit constitutes the edible part, much of the agronomical research is focused on it. However, recent studies have indicated that leaf morphology influences fruit quality. As leaves are specialized photosynthetic organs and the vascular systems transport the photosynthetic products to sink organs, the architectural characteristics of the leaves have a strong influence on the final fruit quality. Therefore, comprehensive research focusing on both the fruit and leaf morphology is required for further tomato breeding. This review summarizes an overview of knowledge of the basic tomato leaf development, morphological diversification, and molecular mechanisms behind them and emphasizes its importance in breeding. Finally, we discuss how these findings and knowledge can be applied to future tomato breeding.

Model-based reconstruction of whole organ growth dynamics reveals invariant patterns in leaf morphogenesis

Plant organ morphogenesis spans several orders of magnitude in time and space. Because of limitations in live-imaging, analysing whole organ growth from initiation to mature stages typically rely on static data sampled from different timepoints and individuals. We introduce a new model-based strategy for dating organs and for reconstructing morphogenetic trajectories over unlimited time windows based on static data. Using this approach, we show that Arabidopsis thaliana leaves are initiated at regular 1-day intervals. Despite contrasted adult morphologies, leaves of different ranks exhibited shared growth dynamics, with linear gradations of growth parameters according to leaf rank. At the sub-organ scale, successive serrations from same or different leaves also followed shared growth dynamics, suggesting that global and local leaf growth patterns are decoupled. Analysing mutants leaves with altered morphology highlighted the decorrelation between adult shapes and morphogenetic trajectories, thus stressing the benefits of our approach in identifying determinants and critical timepoints during organ morphogenesis.

Leaf-size control beyond transcription factors: Compensatory mechanisms

Plant leaves display abundant morphological richness yet grow to characteristic sizes and shapes. Beginning with a small number of undifferentiated founder cells, leaves evolve via a complex interplay of regulatory factors that ultimately influence cell proliferation and subsequent post-mitotic cell enlargement. During their development, a sequence of key events that shape leaves is both robustly executed spatiotemporally following a genomic molecular network and flexibly tuned by a variety of environmental stimuli. Decades of work on Arabidopsis thaliana have revisited the compensatory phenomena that might reflect a general and primary size-regulatory mechanism in leaves. This review focuses on key molecular and cellular events behind the organ-wide scale regulation of compensatory mechanisms. Lastly, emerging novel mechanisms of metabolic and hormonal regulation are discussed, based on recent advances in the field that have provided insights into, among other phenomena, leaf-size regulation.

The transcriptional co-regulators NBCL1 and NBCL2 redundantly coordinate aerial organ development and root nodule identity in legumes

Medicago truncatula NODULE ROOT1 (MtNOOT1) and Pisum sativum COCHLEATA1 (PsCOCH1) are orthologous genes belonging to the NOOT-BOP-COCH-LIKE (NBCL) gene family which encodes key transcriptional co-regulators of plant development. In Mtnoot1 and Pscoch1 mutants, the development of stipules, flowers, and symbiotic nodules is altered. MtNOOT2 and PsCOCH2 represent the single paralogues of MtNOOT1 and PsCOCH1, respectively. In M. truncatula, MtNOOT1 and MtNOOT2 are both required for the establishment and maintenance of symbiotic nodule identity. In legumes, the role of NBCL2 in above-ground development is not known. To better understand the roles of NBCL genes in legumes, we used M. truncatula and P. sativum nbcl mutants, isolated a knockout mutant for the PsCOCH2 locus and generated Pscoch1coch2 double mutants in P. sativum. Our work shows that single Mtnoot2 and Pscoch2 mutants develop wild-type stipules, flowers, and symbiotic nodules. However, the number of flowers was increased and the pods and seeds were smaller compared to the wild type. Furthermore, in comparison to the corresponding nbcl1 single mutants, both the M. truncatula and P. sativum nbcl double mutants show a drastic alteration in stipule, inflorescence, flower, and nodule development. Remarkably, in both M. truncatula and P. sativum nbcl double mutants, stipules are transformed into a range of aberrant leaf-like structures.

Gene expression and DNA methylation altering lead to the high oil content in wild allotetraploid peanut (A. monticola)

INTRODUCTION

The wild allotetraploid peanut Arachis monticola contains a higher oil content than the cultivated allotetraploid Arachis hypogaea. Besides the fact that increasing oil content is the most important peanut breeding objective, a proper understanding of its molecular mechanism controlling oil accumulation is still lacking.

METHODS

We investigated this aspect by performing comparative transcriptomics from developing seeds between three wild and five cultivated peanut varieties.

RESULTS

The analyses not only showed species-specific grouping transcriptional profiles but also detected two gene clusters with divergent expression patterns between two species enriched in lipid metabolism. Further analysis revealed that expression alteration of lipid metabolic genes with co-expressed transcription factors in wild peanut led to enhanced activity of oil biogenesis and retarded the rate of lipid degradation. In addition, bisulfite sequencing was conducted to characterize the variation of DNA methylation between wild allotetraploid (245, WH 10025) and cultivated allotetraploid (Z16, Zhh 7720) genotypes. CG and CHG context methylation was found to antagonistically correlate with gene expression during seed development. Differentially methylated region analysis and transgenic assay further illustrated that variations of DNA methylation between wild and cultivated peanuts could affect the oil content via altering the expression of peroxisomal acyl transporter protein (Araip.H6S1B).

DISCUSSION

From the results, we deduced that DNA methylation may negatively regulate lipid metabolic genes and transcription factors to subtly affect oil accumulation divergence between wild and cultivated peanuts. Our work provided the first glimpse on the regulatory mechanism of gene expression altering for oil accumulation in wild peanut and gene resources for future breeding applications.

Genome-wide characterization of the TALE homeodomain family and the KNOX-BLH interaction network in tomato

Comprehensive yeast and protoplast two-hybrid analyses illustrated the protein-protein interaction network of the TALE homeodomain protein family, KNOX and BLH proteins, in tomato leaf and fruit development. KNOTTED-like (KNOX, KN) proteins and BELL1-like (BLH) proteins, which belong to the same TALE homeodomain family, act together by forming KNOX-BLH heterodimer modules. These modules play crucial roles in regulating multiple developmental processes in plants, like organ differentiation. However, despite the increasing knowledge about individual KNOX and BLH functions, a comprehensive view of their functional protein-protein interaction (PPI) network remains elusive in most plants, including tomato (Solanum lycopersicum), an important model plant to study fruit and leaf development. Here, we characterized eight tomato KNOX genes (SlKN1 to SlKN8) and fourteen tomato BLH genes (SlBLH1 to SlBLH14) by expression profiling, co-expression analysis, and PPI network analysis using two-hybrid techniques in yeasts (Y2H) and protoplasts (P2H). We identified 75 pairwise KNOX-BLH interactions, including ten novel interactors of SlKN2/TKN2, a primary class I KNOX protein, and nine novel interactors of SlKN5, a primary class II KNOX protein. Based on these data, we classified KNOX-BLH modules into several categories, which made us infer the order and combination of the KNOX-BLH modules involved in differentiation processes in leaf and fruit. Notably, the co-expression and interaction of SlKN5 and fruit preferentially expressing BLH1-clade paralogs (SlBLH5/SlBEL11 and SlBLH7) suggest their important roles in regulating fruit differentiation. Furthermore, in silico modeling of the KNOX-BLH modules, sequence analysis, and P2H assay identified several residues and a linker region potentially influencing the affinity of BLHs to KNOXs within their conserved dimerization domains. Together, these findings provide insights into the regulatory mechanism of KNOX-BLH modules underlying tomato organ differentiation.

Understanding of the various aspects of gene regulatory networks related to crop improvement

The hierarchical relationship between transcription factors, associated proteins, and their target genes is defined by a gene regulatory network (GRN). GRNs allow us to understand how the genotype and environment of a plant are incorporated to control the downstream physiological responses. During plant growth or environmental acclimatization, GRNs are diverse and can be differently regulated across tissue types and organs. An overview of recent advances in the development of GRN that speed up basic and applied plant research is given here. Furthermore, the overview of genome and transcriptome involving GRN research along with the exciting advancement and application are discussed. In addition, different approaches to GRN predictions were elucidated. In this review, we also describe the role of GRN in crop improvement, crop plant manipulation, stress responses, speed breeding and identifying genetic variations/locus. Finally, the challenges and prospects of GRN in plant biology are discussed.

Molecular mechanisms underlying leaf development, morphological diversification, and beyond

The basic mechanisms of leaf development have been revealed through a combination of genetics and intense analyses in select model species. The genetic basis for diversity in leaf morphology seen in nature is also being unraveled through recent advances in techniques and technologies related to genomics and transcriptomics, which have had a major impact on these comparative studies. However, this has led to the emergence of new unresolved questions about the mechanisms that generate the diversity of leaf form. Here, we provide a review of the current knowledge of the fundamental molecular genetic mechanisms underlying leaf development with an emphasis on natural variation and conserved gene regulatory networks involved in leaf development. Beyond that, we discuss open questions/enigmas in the area of leaf development, how recent technologies can best be deployed to generate a unified understanding of leaf diversity and its evolution, and what untapped fields lie ahead.

Evolution of miRNA binding sites and regulatory networks in cichlids

The divergence of regulatory regions and gene regulatory network (GRN) rewiring is a key driver of cichlid phenotypic diversity. However, the contribution of miRNA binding site turnover has yet to be linked to GRN evolution across cichlids. Here, we extend our previous studies by analysing the selective constraints driving evolution of miRNA and transcription factor (TF) binding sites of target genes, to infer instances of cichlid GRN rewiring associated with regulatory binding site turnover. Comparative analyses identified increased species-specific networks that are functionally associated to traits of cichlid phenotypic diversity. The evolutionary rewiring is associated with differential models of miRNA and TF binding site turnover, driven by a high proportion of fast-evolving polymorphic sites in adaptive trait genes compared to subsets of random genes. Positive selection acting upon discrete mutations in these regulatory regions is likely to be an important mechanism in rewiring GRNs in rapidly radiating cichlids. Regulatory variants of functionally associated miRNA and TF binding sites of visual opsin genes differentially segregate according to phylogeny and ecology of Lake Malawi species, identifying both rewired e.g. clade-specific and conserved network motifs of adaptive trait associated GRNs. Our approach revealed several novel candidate regulators, regulatory regions and three-node motifs across cichlid genomes with previously reported associations to known adaptive evolutionary traits.

Heinz-resistant tomato cultivars exhibit a lignin-based resistance to field dodder (Cuscuta campestris) parasitism

Cuscuta species (dodders) are agriculturally destructive, parasitic angiosperms. These parasitic plants use haustoria as physiological bridges to extract nutrients and water from hosts. Cuscuta campestris has a broad host range and wide geographical distribution. While some wild tomato relatives are resistant, cultivated tomatoes are generally susceptible to C. campestris infestations. However, some specific Heinz tomato (Solanum lycopersicum) hybrid cultivars exhibit resistance to dodders in the field, but their defense mechanism was previously unknown. Here, we discovered that the stem cortex in these resistant lines responds with local lignification upon C. campestris attachment, preventing parasite entry into the host. Lignin Induction Factor 1 (LIF1, an AP2-like transcription factor), SlMYB55, and Cuscuta R-gene for Lignin-based Resistance 1, a CC-NBS-LRR (CuRLR1) are identified as factors that confer host resistance by regulating lignification. SlWRKY16 is upregulated upon C. campestris infestation and potentially negatively regulates LIF1 function. Intriguingly, CuRLR1 may play a role in signaling or function as an intracellular receptor for receiving Cuscuta signals or effectors, thereby regulating lignification-based resistance. In summary, these four regulators control the lignin-based resistance response in specific Heinz tomato cultivars, preventing C. campestris from parasitizing resistant tomatoes. This discovery provides a foundation for investigating multilayer resistance against Cuscuta species and has potential for application in other essential crops attacked by parasitic plants.

Analyzing the Impact of Greenhouse Planting Strategy and Plant Architecture on Tomato Plant Physiology and Estimated Dry Matter

Determine the level of significance of planting strategy and plant architecture and how they affect plant physiology and dry matter accumulation within greenhouses is essential to actual greenhouse plant management and breeding. We thus analyzed four planting strategies (plant spacing, furrow distance, row orientation, planting pattern) and eight different plant architectural traits (internode length, leaf azimuth angle, leaf elevation angle, leaf length, leaflet curve, leaflet elevation, leaflet number/area ratio, leaflet length/width ratio) with the same plant leaf area using a formerly developed functional-structural model for a Chinese Liaoshen-solar greenhouse and tomato plant, which used to simulate the plant physiology of light interception, temperature, stomatal conductance, photosynthesis, and dry matter. Our study led to the conclusion that the planting strategies have a more significant impact overall on plant radiation, temperature, photosynthesis, and dry matter compared to plant architecture changes. According to our findings, increasing the plant spacing will have the most significant impact to increase light interception. E-W orientation has better total light interception but yet weaker light uniformity. Changes in planting patterns have limited influence on the overall canopy physiology. Increasing the plant leaflet area by leaflet N/A ratio from what we could observe for a rose the total dry matter by 6.6%, which is significantly better than all the other plant architecture traits. An ideal tomato plant architecture which combined all the above optimal architectural traits was also designed to provide guidance on phenotypic traits selection of breeding process. The combined analysis approach described herein established the causal relationship between investigated traits, which could directly apply to provide management and breeding insights on other plant species with different solar greenhouse structures.

Developmental Analysis of Compound Leaf Development in Arachis hypogaea

Leaves are the primary photosynthetic structures, while photosynthesis is the direct motivation of crop yield formation. As a legume plant, peanut (Arachis hypogaea) is one of the most economically essential crops as well as an important source of edible oil and protein. The leaves of A. hypogaea are in the tetrafoliate form, which is different from the trifoliate leaf pattern of Medicago truncatula, a model legume species. In A. hypogaea, an even-pinnate leaf with a pair of proximal and distal leaflets was developed; however, only a single terminal leaflet and a pair of lateral leaflets were formed in the odd-pinnate leaf in M. truncatula. In this study, the development of compound leaf in A. hypogaea was investigated. Transcriptomic profiles revealed that the common and unique differentially expressed genes were identified in a proximal leaflet and a distal leaflet, which provided a research route to understand the leaf development in A. hypogaea. Then, a naturally occurring mutant line with leaf developmental defects in A. hypogaea was obtained, which displayed a pentafoliate form with an extra terminal leaflet. The characterization of the mutant indicated that cytokinin and class I KNOTTED-LIKE HOMEOBOX were involved in the control of compound leaf pattern in A. hypogaea. These results expand our knowledge and provide insights into the molecular mechanism underlying the formation of different compound leaf patterns among species.

Investigating Host and Parasitic Plant Interaction by Tissue-Specific Gene Analyses on Tomato and Cuscuta campestris Interface at Three Haustorial Developmental Stages

Parasitic weeds cause billions of dollars in agricultural losses each year worldwide. Cuscuta campestris (C. campestris), one of the most widespread and destructive parasitic plants in the United States, severely reduces yield in tomato plants. Reducing the spread of parasitic weeds requires understanding the interaction between parasites and hosts. Several studies have identified factors needed for parasitic plant germination and haustorium induction, and genes involved in host defense responses. However, knowledge of the mechanisms underlying the interactions between host and parasitic plants, specifically at the interface between the two organisms, is relatively limited. A detailed investigation of the crosstalk between the host and parasite at the tissue-specific level would enable development of effective parasite control strategies. To focus on the haustorial interface, we used laser-capture microdissection (LCM) with RNA-seq on early, intermediate and mature haustorial stages. In addition, the tomato host tissue that immediately surround the haustoria was collected to obtain tissue- resolution RNA-Seq profiles for C. campestris and tomato at the parasitism interface. After conducting RNA-Seq analysis and constructing gene coexpression networks (GCNs), we identified CcHB7, CcPMEI, and CcERF1 as putative key regulators involved in C. campestris haustorium organogenesis, and three potential regulators, SlPR1, SlCuRe1-like, and SlNLR, in tomatoes that are involved in perceiving signals from the parasite. We used host-induced gene silencing (HIGS) transgenic tomatoes to knock-down the candidate genes in C. campestris and produced CRISPR transgenic tomatoes to knock out candidate genes in tomatoes. The interactions of C. campestris with these transgenic lines were tested and compared with that in wild-type tomatoes. The results of this study reveal the tissue-resolution gene regulatory mechanisms at the parasitic plant-host interface and provide the potential of developing a parasite-resistant system in tomatoes.

Common Themes and Future Challenges in Understanding Gene Regulatory Network Evolution

A major driving force behind the evolution of species-specific traits and novel structures is alterations in gene regulatory networks (GRNs). Comprehending evolution therefore requires an understanding of the nature of changes in GRN structure and the responsible mechanisms. Here, we review two insect pigmentation GRNs in order to examine common themes in GRN evolution and to reveal some of the challenges associated with investigating changes in GRNs across different evolutionary distances at the molecular level. The pigmentation GRN in Drosophila melanogaster and other drosophilids is a well-defined network for which studies from closely related species illuminate the different ways co-option of regulators can occur. The pigmentation GRN for butterflies of the Heliconius species group is less fully detailed but it is emerging as a useful model for exploring important questions about redundancy and modularity in cis-regulatory systems. Both GRNs serve to highlight the ways in which redeployment of trans-acting factors can lead to GRN rewiring and network co-option. To gain insight into GRN evolution, we discuss the importance of defining GRN architecture at multiple levels both within and between species and of utilizing a range of complementary approaches.

Brachypodium distachyon UNICULME4 and LAXATUM-A are redundantly required for development

In cultivated grasses, tillering, leaf, and inflorescence architecture, as well as abscission ability, are major agronomical traits. In barley (Hordeum vulgare), maize (Zea mays), rice (Oryza sativa), and brachypodium (Brachypodium distachyon), NOOT-BOP-COCH-LIKE (NBCL) genes are essential regulators of vegetative and reproductive development. Grass species usually possess 2-4 NBCL copies and until now a single study in O. sativa showed that the disruption of all NBCL genes strongly altered O. sativa leaf development. To improve our understanding of the role of NBCL genes in grasses, we extended the study of the two NBCL paralogs BdUNICULME4 (CUL4) and BdLAXATUM-A (LAXA) in the nondomesticated grass B. distachyon. For this, we applied reversed genetics and generated original B. distachyon single and double nbcl mutants by clustered regularly interspaced short palindromic repeats - CRISPR associated protein 9 (CRISPR-Cas9) approaches and genetic crossing between nbcl targeting induced local lesions in genomes (TILLING) mutants. Through the study of original single laxa CRISPR-Cas9 null alleles, we validated functions previously proposed for LAXA in tillering, leaf patterning, inflorescence, and flower development and also unveiled roles for these genes in seed yield. Furthermore, the characterization of cul4laxa double mutants revealed essential functions for nbcl genes in B. distachyon development, especially in the regulation of tillering, stem cell elongation and secondary cell wall composition as well as for the transition toward the reproductive phase. Our results also highlight recurrent antagonist interactions between NBCLs occurring in multiple aspects of B. distachyon development.

Universal gene co-expression network reveals receptor-like protein genes involved in broad-spectrum resistance in pepper (Capsicum annuum L.)

Receptor-like proteins (RLPs) on plant cells have been implicated in immune responses and developmental processes. Although hundreds of RLP genes have been identified in plants, only a few RLPs have been functionally characterized in a limited number of plant species. Here, we identified RLPs in the pepper (Capsicum annuum) genome and performed comparative transcriptomics coupled with the analysis of conserved gene co-expression networks (GCNs) to reveal the role of core RLP regulators in pepper-pathogen interactions. A total of 102 RNA-seq datasets of pepper plants infected with four pathogens were used to construct CaRLP-targeted GCNs (CaRLP-GCNs). Resistance-responsive CaRLP-GCNs were merged to construct a universal GCN. Fourteen hub CaRLPs, tightly connected with defense-related gene clusters, were identified in eight modules. Based on the CaRLP-GCNs, we evaluated whether hub CaRLPs in the universal GCN are involved in the biotic stress response. Of the nine hub CaRLPs tested by virus-induced gene silencing, three genes (CaRLP264, CaRLP277, and CaRLP351) showed defense suppression with less hypersensitive response-like cell death in race-specific and non-host resistance response to viruses and bacteria, respectively, and consistently enhanced susceptibility to Ralstonia solanacearum and/or Phytophthora capsici. These data suggest that key CaRLPs are involved in the defense response to multiple biotic stresses and can be used to engineer a plant with broad-spectrum resistance. Together, our data show that generating a universal GCN using comprehensive transcriptome datasets can provide important clues to uncover genes involved in various biological processes.

Evolution of polarity protein BASL and the capacity for stomatal lineage asymmetric divisions

Asymmetric and oriented stem cell divisions enable the continued production of patterned tissues. The molecules that guide these divisions include several "polarity proteins" that are localized to discrete plasma membrane domains, are differentially inherited during asymmetric divisions, and whose scaffolding activities can guide division plane orientation and subsequent cell fates. In the stomatal lineages on the surfaces of plant leaves, asymmetric and oriented divisions create distinct cell types in physiologically optimized patterns. The polarity protein BREAKING OF ASYMMETRY IN THE STOMATAL LINEAGE (BASL) is a major regulator of stomatal lineage division and cell fate asymmetries in Arabidopsis, but its role in the stomatal lineages of other plants is unclear. Here, using phylogenetic and functional assays, we demonstrate that BASL is a eudicot-specific polarity protein. Dicot BASL orthologs can polarize in heterologous systems and rescue the Arabidopsis BASL mutant. The more widely distributed BASL-like proteins, although they share BASL's conserved C-terminal domain, are neither polarized nor do they function in asymmetric divisions of the stomatal lineage. Comparison of BASL protein localization and loss of function BASL phenotypes in Arabidopsis and tomato revealed previously unappreciated differences in how asymmetric cell divisions are employed for pattern formation in different species. This multi-species analysis therefore provides insight into the evolution of a unique polarity regulator and into the developmental choices available to cells as they build and pattern tissues.

Leaf form diversification in an ornamental heirloom tomato results from alterations in two different HOMEOBOX genes

Domesticated plants display diverse phenotypic traits. However, the influence of breeding effort on this phenotypic diversity remains unknown. Here, we demonstrate that a single nucleotide deletion in the homeobox motif of BIPINNATA, a BEL-LIKE HOMEODOMAIN gene, led to a highly complex leaf phenotype in an heirloom tomato (Solanum lycopersicum), Silvery Fir Tree (SiFT), which is used as a landscaping and ornamental plant. A comparative gene network analysis revealed that repression of SOLANIFOLIA, the ortholog of WUSCHEL RELATED HOMEOBOX 1, caused the narrow leaflet phenotype seen in SiFT. Comparative genomics indicated that the bip mutation in SiFT likely arose de novo and is unique to SiFT and not introgressed from other tomato genomes. These results provide new insights into the natural variation in phenotypic traits introduced into crops during improvement processes after domestication and establish homeobox genes as evolutionary hotspots.

Leaf form diversification in an ornamental heirloom tomato results from alterations in two different HOMEOBOX genes

Domesticated plants display diverse phenotypic traits. However, the influence of breeding effort on this phenotypic diversity remains unknown. Here, we demonstrate that a single nucleotide deletion in the homeobox motif of BIPINNATA, a BEL-LIKE HOMEODOMAIN gene, led to a highly complex leaf phenotype in an heirloom tomato (Solanum lycopersicum), Silvery Fir Tree (SiFT), which is used as a landscaping and ornamental plant. A comparative gene network analysis revealed that repression of SOLANIFOLIA, the ortholog of WUSCHEL RELATED HOMEOBOX 1, caused the narrow leaflet phenotype seen in SiFT. Comparative genomics indicated that the bip mutation in SiFT likely arose de novo and is unique to SiFT and not introgressed from other tomato genomes. These results provide new insights into the natural variation in phenotypic traits introduced into crops during improvement processes after domestication and establish homeobox genes as evolutionary hotspots.

Meristem transitions and plant architecture-learning from domestication for crop breeding

Genetic networks that regulate meristem transitions were recurrent targets of selection during crop domestication and allow fine-tuning of plant architecture for improved crop productivity.

Digital morphometrics: Application of MorphoLeaf in shape visualization and species delimitation, using Cucurbitaceae leaves as a model

PREMISE

Plant leaves are one of the most important organs for plant identification due to their variability across different taxonomic groups. While traditional morphometrics has contributed tremendously to reducing the problems accompanying plant identification and morphology-based species delimitation, image-analysis digital solutions have made it easy to detect more characters to complement existing leaf data sets.

METHODS

Here, we apply MorphoLeaf to generate a morphometric data set from 140 leaf specimens of seven Cucurbitaceae species via landmark extraction, the reparameterization of leaf contours, and data quantification and normalization. A statistical analysis was performed on the resulting data set.

RESULTS

A principal component analysis revealed that leaf blade area, blade perimeter, tooth area, tooth perimeter, the measure of the distance from tooth position to the tip, and the measure of the distance from tooth position to the base are important and informative landmarks that contribute to the variation within the species studied.

DISCUSSION

MorphoLeaf can be applied to quantitatively track leaf diversity, thereby functionally integrating morphometrics and shape visualization into the digital identification of plants. The success of digital morphometrics in leaf outline analyses presents researchers with opportunities to carry out more accurate image-based research in areas such as plant development, evolution, and phenotyping.

LATERAL ORGAN BOUNDARIES DOMAIN 25 functions as a key regulator of haustorium development in dodders

Parasitic plants reduce crop yield worldwide. Dodder (Cuscuta campestris) is a stem parasite that attaches to its host, using haustoria to extract nutrients and water. We analyzed the transcriptome of six C. campestris tissues and identified a key gene, LATERAL ORGAN BOUNDARIES DOMAIN 25 (CcLBD25), as highly expressed in prehaustoria and haustoria. Gene coexpression networks from different tissue types and laser-capture microdissection RNA-sequencing data indicated that CcLBD25 could be essential for regulating cell wall loosening and organogenesis. We employed host-induced gene silencing by generating transgenic tomato (Solanum lycopersicum) hosts that express hairpin RNAs to target and down-regulate CcLBD25 in the parasite. Our results showed that C. campestris growing on CcLBD25 RNAi transgenic tomatoes transited to the flowering stage earlier and had reduced biomass compared with C. campestris growing on wild-type (WT) hosts, suggesting that parasites growing on transgenic plants were stressed due to insufficient nutrient acquisition. We developed an in vitro haustorium system to assay the number of prehaustoria produced on strands from C. campestris. Cuscuta campestris grown on CcLBD25 RNAi tomatoes produced fewer prehaustoria than those grown on WT tomatoes, indicating that down-regulating CcLBD25 may affect haustorium initiation. Cuscuta campestris haustoria growing on CcLBD25 RNAi tomatoes exhibited reduced pectin digestion and lacked searching hyphae, which interfered with haustorium penetration and formation of vascular connections. The results of this study elucidate the role of CcLBD25 in haustorium development and might contribute to developing parasite-resistant crops.

Evolutionarily conserved hierarchical gene regulatory networks for plant salt stress response

Plant cells constantly alter their gene expression profiles to respond to environmental fluctuations. These continuous adjustments are regulated by multi-hierarchical networks of transcription factors. To understand how such gene regulatory networks (GRNs) have stabilized evolutionarily while allowing for species-specific responses, we compare the GRNs underlying salt response in the early-diverging and late-diverging plants Marchantia polymorpha and Arabidopsis thaliana. Salt-responsive GRNs, constructed on the basis of the temporal transcriptional patterns in the two species, share common trans-regulators but exhibit an evolutionary divergence in cis-regulatory sequences and in the overall network sizes. In both species, WRKY-family transcription factors and their feedback loops serve as central nodes in salt-responsive GRNs. The divergent cis-regulatory sequences of WRKY-target genes are probably associated with the expansion in network size, linking salt stress to tissue-specific developmental and physiological responses. The WRKY modules and highly linked WRKY feedback loops have been preserved widely in other plants, including rice, while keeping their binding-motif sequences mutable. Together, the conserved trans-regulators and the quickly evolving cis-regulatory sequences allow salt-responsive GRNs to adapt over a long evolutionary timescale while maintaining some consistent regulatory structure. This strategy may benefit plants as they adapt to changing environments.

Coordinating the morphogenesis-differentiation balance by tweaking the cytokinin-gibberellin equilibrium

Morphogenesis and differentiation are important stages in organ development and shape determination. However, how they are balanced and tuned during development is not fully understood. In the compound leaved tomato, an extended morphogenesis phase allows for the initiation of leaflets, resulting in the compound form. Maintaining a prolonged morphogenetic phase in early stages of compound-leaf development in tomato is dependent on delayed activity of several factors that promote differentiation, including the CIN-TCP transcription factor (TF) LA, the MYB TF CLAU and the plant hormone Gibberellin (GA), as well as on the morphogenesis-promoting activity of the plant hormone cytokinin (CK). Here, we investigated the genetic regulation of the morphogenesis-differentiation balance by studying the relationship between LA, CLAU, TKN2, CK and GA. Our genetic and molecular examination suggest that LA is expressed earlier and more broadly than CLAU and determines the developmental context of CLAU activity. Genetic interaction analysis indicates that LA and CLAU likely promote differentiation in parallel genetic pathways. These pathways converge downstream on tuning the balance between CK and GA. Comprehensive transcriptomic analyses support the genetic data and provide insights into the broader molecular basis of differentiation and morphogenesis processes in plants.

Morphogenesis and differentiation are crucial steps in the formation and shaping of organs in both plants and animals. A wide array of transcription factors and hormones were shown to act together to support morphogenesis or promote differentiation. However, a comprehensive molecular and genetic understating of how morphogenesis and differentiation are coordinated during development is still missing. We addressed these questions in the context of the development of the tomato compound leaf, for which many regulators have been described. Investigating the coordination among these different actors, we show that several discrete genetic pathways promote differentiation. Downstream of these separate pathways, two important plant hormones, cytokinin and gibberellin, act antagonistically to tweak the morphogenesis-differentiation balance.

Spatial transcriptional signatures define margin morphogenesis along the proximal-distal and medio-lateral axes in tomato (Solanum lycopersicum) leaves

Leaf morphogenesis involves cell division, expansion, and differentiation in the developing leaf, which take place at different rates and at different positions along the medio-lateral and proximal-distal leaf axes. The gene expression changes that control cell fate along these axes remain elusive due to difficulties in precisely isolating tissues. Here, we combined rigorous early leaf characterization, laser capture microdissection, and transcriptomic sequencing to ask how gene expression patterns regulate early leaf morphogenesis in wild-type tomato (Solanum lycopersicum) and the leaf morphogenesis mutant trifoliate. We observed transcriptional regulation of cell differentiation along the proximal-distal axis and identified molecular signatures delineating the classically defined marginal meristem/blastozone region during early leaf development. We describe the role of endoreduplication during leaf development, when and where leaf cells first achieve photosynthetic competency, and the regulation of auxin transport and signaling along the leaf axes. Knockout mutants of BLADE-ON-PETIOLE2 exhibited ectopic shoot apical meristem formation on leaves, highlighting the role of this gene in regulating margin tissue identity. We mapped gene expression signatures in specific leaf domains and evaluated the role of each domain in conferring indeterminacy and permitting blade outgrowth. Finally, we generated a global gene expression atlas of the early developing compound leaf.

Insights into the evo-devo of plant reproduction using next-generation sequencing approaches

The development of plant model organisms has traditionally been analyzed using resource-heavy, tailored applications that are not easily transferable to distantly related non-model taxa. Thus, our understanding of plant development has been limited to a subset of traits, and evolutionary studies conducted most effectively either across very wide [e.g. Arabidopsis thaliana and Oryza sativa (rice)] or narrow (i.e. population level) phylogenetic distances. As plant biologists seek to capitalize on natural diversity for crop improvement, enhance ecosystem functioning, and better understand plant responses to climate change, high-throughput and broadly applicable forms of existing molecular biology assays are becoming an invaluable resource. Next-generation sequencing (NGS) is increasingly becoming a powerful tool in evolutionary developmental biology (evo-devo) studies, particularly through its application to understanding trait evolution at different levels of gene regulation. Here, I review some of the most common and emerging NGS-based methods, using exemplar studies in reproductive plant evo-devo to illustrate their potential.

Genomic Analyses of Phenotypic Differences Between Native and Invasive Populations of Diffuse Knapweed (Centaurea diffusa)

Invasive species represent excellent opportunities to study the evolutionary potential of traits important to success in novel environments. Although some ecologically important traits have been identified in invasive species, little is typically known about the genetic mechanisms that underlie invasion success in non-model species. Here, we use a genome-wide association (GWAS) approach to identify the genetic basis of trait variation in the non-model, invasive, diffuse knapweed [Centaurea diffusa Lam. (Asteraceae)]. To assist with this analysis, we have assembled the first draft genome reference and fully annotated plastome assembly for this species, and one of the first from this large, weedy, genus, which is of major ecological and economic importance. We collected phenotype data from 372 individuals from four native and four invasive populations of C. diffusa grown in a common environment. Using these individuals, we produced reduced-representation genotype-by-sequencing (GBS) libraries and identified 7,058 SNPs. We identify two SNPs associated with leaf width in these populations, a trait which significantly varies between native and invasive populations. In this rosette forming species, increased leaf width is a major component of increased biomass, a common trait in invasive plants correlated with increased fitness. Finally, we use annotations from Arabidopsis thaliana to identify 98 candidate genes that are near the associated SNPs and highlight several good candidates for leaf width variation.

Evolution of regulatory networks associated with traits under selection in cichlids

BACKGROUND

Seminal studies of vertebrate protein evolution speculated that gene regulatory changes can drive anatomical innovations. However, very little is known about gene regulatory network (GRN) evolution associated with phenotypic effect across ecologically diverse species. Here we use a novel approach for comparative GRN analysis in vertebrate species to study GRN evolution in representative species of the most striking examples of adaptive radiations, the East African cichlids. We previously demonstrated how the explosive phenotypic diversification of East African cichlids can be attributed to diverse molecular mechanisms, including accelerated regulatory sequence evolution and gene expression divergence.

RESULTS

To investigate these mechanisms across species at a genome-wide scale, we develop a novel computational pipeline that predicts regulators for co-extant and ancestral co-expression modules along a phylogeny, and candidate regulatory regions associated with traits under selection in cichlids. As a case study, we apply our approach to a well-studied adaptive trait-the visual system-for which we report striking cases of network rewiring for visual opsin genes, identify discrete regulatory variants, and investigate their association with cichlid visual system evolution. In regulatory regions of visual opsin genes, in vitro assays confirm that transcription factor binding site mutations disrupt regulatory edges across species and segregate according to lake species phylogeny and ecology, suggesting GRN rewiring in radiating cichlids.

CONCLUSIONS

Our approach reveals numerous novel potential candidate regulators and regulatory regions across cichlid genomes, including some novel and some previously reported associations to known adaptive evolutionary traits.

Reproductive Proteins Evolve Faster Than Non-reproductive Proteins Among Solanum Species

Elevated rates of evolution in reproductive proteins are commonly observed in animal species, and are thought to be driven by the action of sexual selection and sexual conflict acting specifically on reproductive traits. Whether similar patterns are broadly observed in other biological groups is equivocal. Here, we examine patterns of protein divergence among wild tomato species (Solanum section Lycopersicon), to understand forces shaping the evolution of reproductive genes in this diverse, rapidly evolving plant clade. By comparing rates of molecular evolution among loci expressed in reproductive and non-reproductive tissues, our aims were to test if: (a) reproductive-specific loci evolve more rapidly, on average, than non-reproductive loci; (b) 'male'-specific loci evolve at different rates than 'female'-specific loci; (c) genes expressed exclusively in gametophytic (haploid) tissue evolve differently from genes expressed in sporophytic (diploid) tissue or in both tissue types; and (d) mating system variation (a potential proxy for the expected strength of sexual selection and/or sexual conflict) affects patterns of protein evolution. We observed elevated evolutionary rates in reproductive proteins. However, this pattern was most evident for female- rather than male-specific loci, both broadly and for individual loci inferred to be positively selected. These elevated rates might be facilitated by greater tissue-specificity of reproductive proteins, as faster rates were also associated with more narrow expression domains. In contrast, we found little evidence that evolutionary rates are consistently different in loci experiencing haploid selection (gametophytic-exclusive loci), or in lineages with quantitatively different mating systems. Overall while reproductive protein evolution is generally elevated in this diverse plant group, some specific patterns of evolution are more complex than those reported in other (largely animal) systems, and include a more prominent role for female-specific loci among adaptively evolving genes.

Co-expression Networks From Gene Expression Variability Between Genetically Identical Seedlings Can Reveal Novel Regulatory Relationships

Co-expression networks are a powerful tool to understand gene regulation. They have been used to identify new regulation and function of genes involved in plant development and their response to the environment. Up to now, co-expression networks have been inferred using transcriptomes generated on plants experiencing genetic or environmental perturbation, or from expression time series. We propose a new approach by showing that co-expression networks can be constructed in the absence of genetic and environmental perturbation, for plants at the same developmental stage. For this, we used transcriptomes that were generated from genetically identical individual plants that were grown under the same conditions and for the same amount of time. Twelve time points were used to cover the 24-h light/dark cycle. We used variability in gene expression between individual plants of the same time point to infer a co-expression network. We show that this network is biologically relevant and use it to suggest new gene functions and to identify new targets for the transcriptional regulators GI, PIF4, and PRR5. Moreover, we find different co-regulation in this network based on changes in expression between individual plants, compared to the usual approach requiring environmental perturbation. Our work shows that gene co-expression networks can be identified using variability in gene expression between individual plants, without the need for genetic or environmental perturbations. It will allow further exploration of gene regulation in contexts with subtle differences between plants, which could be closer to what individual plants in a population might face in the wild.

Leaf shape in Populus tremula is a complex, omnigenic trait

Leaf shape is a defining feature of how we recognize and classify plant species. Although there is extensive variation in leaf shape within many species, few studies have disentangled the underlying genetic architecture. We characterized the genetic architecture of leaf shape variation in Eurasian aspen (Populus tremula L.) by performing genome-wide association study (GWAS) for physiognomy traits. To ascertain the roles of identified GWAS candidate genes within the leaf development transcriptional program, we generated RNA-Seq data that we used to perform gene co-expression network analyses from a developmental series, which is publicly available within the PlantGenIE resource. We additionally used existing gene expression measurements across the population to analyze GWAS candidate genes in the context of a population-wide co-expression network and to identify genes that were differentially expressed between groups of individuals with contrasting leaf shapes. These data were integrated with expression GWAS (eQTL) results to define a set of candidate genes associated with leaf shape variation. Our results identified no clear adaptive link to leaf shape variation and indicate that leaf shape traits are genetically complex, likely determined by numerous small-effect variations in gene expression. Genes associated with shape variation were peripheral within the population-wide co-expression network, were not highly connected within the leaf development co-expression network, and exhibited signatures of relaxed selection. As such, our results are consistent with the omnigenic model.

Contributions of Adaptive Plant Architecture to Transgressive Salinity Tolerance in Recombinant Inbred Lines of Rice: Molecular Mechanisms Based on Transcriptional Networks

Genetic novelties are important nucleators of adaptive speciation. Transgressive segregation is a major mechanism that creates genetic novelties with morphological and developmental attributes that confer adaptive advantages in certain environments. This study examined the morpho-developmental and physiological profiles of recombinant inbred lines (RILs) from the salt-sensitive IR29 and salt-tolerant Pokkali rice, representing the total range of salt tolerance including the outliers at both ends of the spectrum. Morpho-developmental and physiological profiles were integrated with a hypothesis-driven interrogation of mRNA and miRNA transcriptomes to uncover the critical genetic networks that have been rewired for novel adaptive architecture. The transgressive super-tolerant FL510 had a characteristic small tiller angle and wider, more erect, sturdier, and darker green leaves. This unique morphology resulted in lower transpiration rate, which also conferred a special ability to retain water more efficiently for osmotic avoidance. The unique ability for water retention conferred by such adaptive morphology appeared to enhance the efficacy of defenses mediated by Na+ exclusion mechanism (SalTol-effects) inherited from Pokkali. The super-tolerant FL510 and super-sensitive FL499 had the smallest proportions of differentially expressed genes with little overlaps. Genes that were steadily upregulated in FL510 comprised a putative cytokinin-regulated genetic network that appeared to maintain robust growth under salt stress through well-orchestrated cell wall biogenesis and cell expansion, likely through major regulatory (OsRR23, OsHK5) and biosynthetic (OsIPT9) genes in the cytokinin signaling pathway. Meanwhile, a constitutively expressed cluster in FL510 prominently featured two transcription factors (OsIBH1, TAC3) that control tiller angle and growth habit through the brassinosteroid signaling pathway. Both the putative cytokinin-mediated and brassinosteroid-mediated clusters appeared to function as highly coordinated network synergies in FL510. In contrast, both networks appeared to be sub-optimal and inferior in the other RILs and parents as they were disjointed and highly fragmented. Transgressively expressed miRNAs (miR169, miR397, miR827) were also identified as prominent signatures of FL510, with functional implications to mechanisms that support robust growth, homeostasis, and osmotic stress avoidance. Results of this study demonstrate how genetic recombination creates novel morphology that complements inducible defenses hence transgressive adaptive phenotypes.

Conservation and Divergence in Duplicated Fiber Coexpression Networks Accompanying Domestication of the Polyploid Gossypium hirsutum L

Gossypium hirsutum L. (Upland cotton) has an evolutionary history involving inter-genomic hybridization, polyploidization, and subsequent domestication. We analyzed the developmental dynamics of the cotton fiber transcriptome accompanying domestication using gene coexpression networks for both joint and homoeologous networks. Remarkably, most genes exhibited expression for at least one homoeolog, confirming previous reports of widespread gene usage in cotton fibers. Most coexpression modules comprising the joint network are preserved in each subgenomic network and are enriched for similar biological processes, showing a general preservation of network modular structure for the two co-resident genomes in the polyploid. Interestingly, only one fifth of homoeologs co-occur in the same module when separated, despite similar modular structures between the joint and homoeologous networks. These results suggest that the genome-wide divergence between homoeologous genes is sufficient to separate their co-expression profiles at the intermodular level, despite conservation of intramodular relationships within each subgenome. Most modules exhibit D-homoeolog expression bias, although specific modules do exhibit A-homoeolog bias. Comparisons between wild and domesticated coexpression networks revealed a much tighter and denser network structure in domesticated fiber, as evidenced by its fewer modules, 13-fold increase in the number of development-related module member genes, and the poor preservation of the wild network topology. These results demonstrate the amazing complexity that underlies the domestication of cotton fiber.

Transcriptomic variation of the flower-fruit transition in Physalis and Solanum

Gene expression variations in response to fertilization between Physalis and Solanum might play essential roles in species divergence and fruit evolution. Fertilization triggers variation in fruit development and morphology. The Chinese lantern, a morphological novelty derived from the calyx, is formed upon fertilization in Physalis but is not observed in Solanum. The underlying genetic variations are largely unknown. Here, we documented the developmental and morphological differences in the flower and fruit between Physalis floridana and Solanum pimpinellifolium and then evaluated both the transcript sequence variation and gene expression at the transcriptomic level at fertilization between the two species. In Physalis transcriptomic analysis, 468 unigenes were identified as differentially expressed genes (DEGs) that were strongly regulated by fertilization across 3 years. In comparison with tomato, 14,536 strict single-copy orthologous gene pairs were identified between P. floridana and S. pimpinellifolium in the flower-fruit transcriptome. Nine types of gene variations with specific GO-enriched patterns were identified, covering 58.82% orthologous gene pairs that were DEGs in either trend or dosage at the flower-fruit transition between the two species, which could adequately distinguish Solanum and Physalis, implying that differential gene expression at fertilization might play essential roles during the divergence and fruit evolution of Solanum-Physalis. Virus-induced gene silencing analyses revealed the developmental roles of some transcription factor genes in fertility, Chinese lantern development, and fruit weight control in Physalis. This study presents the first floral transcriptomic resource of Physalis, and reveals some candidate genetic variations accounting for the early fruit developmental evolution in Physalis in comparison to Solanum.

Roles of BdUNICULME4 and BdLAXATUM-A in the non-domesticated grass Brachypodium distachyon

In cultivated grasses, tillering, spike architecture and seed shattering represent major agronomical traits. In barley, maize and rice, the NOOT-BOP-COCH-LIKE (NBCL) genes play important roles in development, especially in ligule development, tillering and flower identity. However, compared with dicots, the role of grass NBCL genes is underinvestigated. To better understand the role of grass NBCLs and to overcome any effects of domestication that might conceal their original functions, we studied TILLING nbcl mutants in the non-domesticated grass Brachypodium distachyon. In B. distachyon, the NBCL genes BdUNICULME4 (CUL4) and BdLAXATUM-A (LAXA) are orthologous, respectively, to the barley HvUniculme4 and HvLaxatum-a, to the maize Zmtassels replace upper ears1 and Zmtassels replace upper ears2 and to the rice OsBLADE-ON-PETIOLE1 and OsBLADE-ON-PETIOLE2/3. In B. distachyon, our reverse genetics study shows that CUL4 is not essential for the establishment of the blade-sheath boundary but is necessary for the development of the ligule and auricles. We report that CUL4 also exerts a positive role in tillering and a negative role in spikelet meristem activity. On the other hand, we demonstrate that LAXA plays a negative role in tillering, positively participates in spikelet development and contributes to the control of floral organ number and identity. In this work, we functionally characterized two new NBCL genes in a context of non-domesticated grass and highlighted original roles for grass NBCL genes that are related to important agronomical traits.

Co-expression clustering across flower development identifies modules for diverse floral forms in Achimenes (Gesneriaceae)

BACKGROUND

Genetic pathways involved with flower color and shape are thought to play an important role in the development of flowers associated with different pollination syndromes, such as those associated with bee, butterfly, or hummingbird pollination. Because pollination syndromes are complex traits that are orchestrated by multiple genes and pathways, the gene regulatory networks have not been explored. Gene co-expression networks provide a systems level approach to identify important contributors to floral diversification.

METHODS

RNA-sequencing was used to assay gene expression across two stages of flower development (an early bud and an intermediate stage) in 10 species of Achimenes (Gesneriaceae). Two stage-specific co-expression networks were created from 9,503 orthologs and analyzed to identify module hubs and the network periphery. Module association with bee, butterfly, and hummingbird pollination syndromes was tested using phylogenetic mixed models. The relationship between network connectivity and evolutionary rates (d N/d S) was tested using linear models.

RESULTS

Networks contained 65 and 62 modules that were largely preserved between developmental stages and contained few stage-specific modules. Over a third of the modules in both networks were associated with flower color, shape, and pollination syndrome. Within these modules, several hub nodes were identified that related to the production of anthocyanin and carotenoid pigments and the development of flower shape. Evolutionary rates were decreased in highly connected genes and elevated in peripheral genes.

DISCUSSION

This study aids in the understanding of the genetic architecture and network properties underlying the development of floral form and provides valuable candidate modules and genes for future studies.

The remarkable morphological diversity of leaf shape in sweet potato (Ipomoea batatas): the influence of genetics, environment, and G×E

Leaf shape, a spectacularly diverse plant trait, varies across taxonomic levels, geography and in response to environmental differences. However, comprehensive intraspecific analyses of leaf shape variation across variable environments is surprisingly absent. Here, we performed a multilevel analysis of leaf shape using diverse accessions of sweet potato (Ipomoea batatas), and uncovered the role of genetics, environment, and G×E on this important trait. We examined leaf shape using a variety of morphometric analyses, and complement this with a transcriptomic survey to identify gene expression changes associated with shape variation. Additionally, we examined the role of genetics and environment on leaf shape by performing field studies in two geographically separate common gardens. We showed that extensive leaf shape variation exists within I. batatas, and identified promising candidate genes associated with this variation. Interestingly, when considering traditional measures, we found that genetic factors are largely responsible for most of leaf shape variation, but that the environment is highly influential when using more quantitative measures via leaf outlines. This extensive and multilevel examination of leaf shape shows an important role of genetics underlying a potentially important agronomic trait, and highlights that the environment can be a strong influence when using more quantitative measures of leaf shape.

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.

Overexpression of a Novel LcKNOX Transcription Factor from Liriodendron chinense Induces Lobed Leaves in Arabidopsis thaliana

Liriodendron chinense is a common ornamental tree that has attractive leaves, which is a valuable trait for use in landscape architecture. In this work, we aimed to identify the potential genes that control and regulate the development of L. chinense leaf lobes. Based on the transcriptome data for the leaf developmental stages we previously generated, two candidate genes were identified in this study. KNOTTED-LIKE HOMEOBOX(KNOX), encoding homeobox family proteins, play a large role in leaf lobe and leaf complexity regulation. Here, two full length KNOX genes from L. chinense were amplified and named LcKNOX1 and LcKNOX6 according to their sequence similarities with the respective Arabidopsis thaliana KNOX family genes. Overexpression vectors were constructed and subsequently transformed into wild type (WT) A. thaliana. Additionally, LcKNOX6 was expressed in tobacco leaves to examine its subcellular localization, and the 35S::LcKNOX6 transgenic A. thaliana leaf cells were imaged with the use of SEM. The expression of several genes that participate in KNOX gene regulation were validated by quantitative real-time PCR. The results show that LcKNOX1 produces almost the same phenotype as that found in WT A. thaliana. Notably, the LcKNOX6-1 lines presented deep leaf lobes that were similar to L. chinense leaf lobes. Two 35S::LcKNOX6 lines induced an abnormal growth phenotype whose seeds were abortive. In short, these results indicate that the LcKNOX6 gene might affect leaf development in A. thaliana and provide insights into the regulation of L. chinense leaf shaping.

Into the bloom: Molecular response of pelagic tunicates to fluctuating food availability

The planktonic tunicates appendicularians and thaliaceans are highly efficient filter feeders on a wide range of prey size including bacteria and have shorter generation times than any other marine grazers. These traits allow some tunicate species to reach high population densities and ensure their success in a favourable environment. However, there are still few studies focusing on which genes and gene pathways are associated with responses of pelagic tunicates to environmental variability. Herein, we present the effect of food availability increase on tunicate community and gene expression at the Marquesas Islands (South-East Pacific Ocean). By using data from the Tara Oceans expedition, we show that changes in phytoplankton density and composition trigger the success of a dominant larvacean species (an undescribed appendicularian). Transcriptional signature to the autotroph bloom suggests key functions in specific physiological processes, i.e., energy metabolism, muscle contraction, membrane trafficking, and proteostasis. The relative abundance of reverse transcription-related Pfams was lower at bloom conditions, suggesting a link with adaptive genetic diversity in tunicates in natural ecosystems. Downstream of the bloom, pelagic tunicates were outcompeted by copepods. Our work represents the first metaomics study of the biological effects of phytoplankton bloom on a key zooplankton taxon.

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.

Coordination of Leaf Development Across Developmental Axes

Leaves are initiated as lateral outgrowths from shoot apical meristems throughout the vegetative life of the plant. To achieve proper developmental patterning, cell-type specification and growth must occur in an organized fashion along the proximodistal (base-to-tip), mediolateral (central-to-edge), and adaxial–abaxial (top-bottom) axes of the developing leaf. Early studies of mutants with defects in patterning along multiple leaf axes suggested that patterning must be coordinated across developmental axes. Decades later, we now recognize that a highly complex and interconnected transcriptional network of patterning genes and hormones underlies leaf development. Here, we review the molecular genetic mechanisms by which leaf development is coordinated across leaf axes. Such coordination likely plays an important role in ensuring the reproducible phenotypic outcomes of leaf morphogenesis.

IBR5 Regulates Leaf Serrations Development via Modulation of the Expression of PIN1

Biodiversity in plant shape is mainly attributable to the diversity of leaf shape, which is largely determined by the transient morphogenetic activity of the leaf margin that creates leaf serrations. However, the precise mechanism underlying the establishment of this morphogenetic capacity remains poorly understood. We report here that INDOLE-3-BUTYRIC ACID RESPONSE 5 (IBR5), a dual-specificity phosphatase, is a key component of leaf-serration regulatory machinery. Loss-of-function mutants of IBR5 exhibited pronounced serrations due to increased cell area. IBR5 was localized in the nucleus of leaf epidermis and petiole cells. Introducing a C129S mutation within the highly conserved VxVHCx₂GxSRSx₅AYLM motif of IBR5 rendered it unable to rescue the leaf-serration defects of the ibr5-3 mutant. In addition, auxin reporters revealed that the distribution of auxin maxima was expanded ectopically in ibr5-3. Furthermore, we found that the distribution of PIN1 on the plasma membrane of the epidermal and cells around the leaf vein was compromised in ibr5-3. We concluded that IBR5 is essential for the establishment of PIN-FORMED 1 (PIN1)-directed auxin maxima at the tips of leaf serration, which is vital for the elaborated regulation during its formation.

BLADE-ON-PETIOLE genes are not involved in the transition from protonema to gametophore in the moss Physcomitrella patens

The timing of the transition between developmental phases is a critical determinant of plant form. In the moss Physcomitrella patens, the transition from protonema to gametophore is a particularly important step as it results in a change from two-dimensional to three-dimensional growth of the plant body. It is well known that this transition is promoted by cytokinin (CK), however, the underlying mechanisms are poorly understood. Previously, it was reported that P. patens orthologs of BLADE-ON-PETIOLE (BOP) genes (PpBOPs) work downstream of CK to promote the transition to gametophore. To further understand the role of PpBOPs in the control of this transition, we performed functional analyses of PpBOP genes. We simultaneously disrupted the function of all three PpBOP genes in P. patens using CRISPR technology, however, no abnormal phenotypes were observed in the triple mutant during either the gametophytic or the sporophytic growth stages. CK treatment did not alter the phase change in the triple mutant. We conclude that PpBOP genes are unnecessary in the control of P. patens development under normal conditions. We propose that BOP genes are not involved in the control of developmental processes in bryophytes and other basal land plants, but may function in physiological processes such as in the defense response.