De novo Sequencing and Comparative Transcriptomics of Floral Development of the Distylous Species Lithospermum multiflorum

Probing of plant transcriptomes reveals the hidden genetic diversity of the family Secoviridae

Secoviruses are single-stranded RNA viruses that infect plants. In the present study, we identified 61 putative novel secoviral genomes in various plant species by mining publicly available plant transcriptome data. These viral sequences represent the genomes of 13 monopartite and 48 bipartite secovirids. The genome sequences of 52 secovirids were coding-complete, and nine were partial. Except for small open reading frames (ORFs) determined in waikaviral genomes and RNA2 of torradoviruses, all of the recovered genomes/genome segments contained a large ORF encoding a polyprotein. Based on genome organization and phylogeny, all but three of the novel secoviruses were assigned to different genera. The genome organization of two identified waika-like viruses resembled that of the recently identified waika-like virus Triticum aestivum secovirus. Phylogenetic analysis revealed a pattern of host-virus co-evolution in a few waika- and waika-like viruses and increased phylogenetic diversity of nepoviruses. The study provides a basis for further investigation of the biological properties of these novel secoviruses.

Expanding the Repertoire of the Plant-Infecting Ophioviruses through Metatranscriptomics Data

Ophioviruses (genus Ophiovirus, family Aspiviridae) are plant-infecting viruses with non-enveloped, filamentous, naked nucleocapsid virions. Members of the genus Ophiovirus have a segmented single-stranded negative-sense RNA genome (ca. 11.3–12.5 kb), encompassing three or four linear segments. In total, these segments encode four to seven proteins in the sense and antisense orientation, both in the viral and complementary strands. The genus Ophiovirus includes seven species with viruses infecting both monocots and dicots, mostly trees, shrubs and some ornamentals. From a genomic perspective, as of today, there are complete genomes available for only four species. Here, by exploring large publicly available metatranscriptomics datasets, we report the identification and molecular characterization of 33 novel viruses with genetic and evolutionary cues of ophioviruses. Genetic distance and evolutionary insights suggest that all the detected viruses could correspond to members of novel species, which expand the current diversity of ophioviruses ca. 4.5-fold. The detected viruses increase the tentative host range of ophioviruses for the first time to mosses, liverwort and ferns. In addition, the viruses were linked to several Asteraceae, Orchidaceae and Poaceae crops/ornamental plants. Phylogenetic analyses showed a novel clade of mosses, liverworts and fern ophioviruses, characterized by long branches, suggesting that there is still plenty of unsampled hidden diversity within the genus. This study represents a significant expansion of the genomics of ophioviruses, opening the door to future works on the molecular and evolutionary peculiarity of this virus genus.

Morphological and physiological mechanism of cytoplasmic inheritance stigma exsertion trait expression in tobacco (Nicotiana tabacum)

Stigma exsertion is an essential outcrossing trait that can improve hybrid seed production efficiencies. In this study, the morphological and physiological mechanisms of cytoplasmic inheritance stigma exsertion trait expression in a tobacco line (MSK326SE) which generated from a stigma exsertion tobacco mutant through continuous backcross were investigated. Compared with its homonuclear-heteroplasmic lines (MSK326 and K326 with inserted stigmas), the exserted stigma phenotype of MSK326SE was mainly caused by corolla shortening, while was stable under different environmental temperature. The different responses of mainly endogenous hormones and expression of cell division- and expansion-related genes caused the differences in cell division and expansion in different flower organs, which further determined the lengths of the corolla. Furthermore, the significant decrease of MSK326SE corolla epidermal cell size caused corolla shortening and finally resulting in stigma exsertion. Exogenous JA could shorten the corolla and more effective increased stigma exsertion degree of MSK326SE, suggesting a potential relationship between stigma exsertion and high JA levels during early bud development. The hybrid seed production efficiency could be improved in tobacco. Our results provide a basis for elucidating the cytoplasmic inheritance stigma exsertion trait expression in tobacco while helping to improve hybrid seed production efficiency.

Comparative transcriptomics reveals commonalities and differences in the genetic underpinnings of a floral dimorphism

Distyly, a floral dimorphism associated with heteromorphic self-incompatibility and controlled by the S-locus supergene, evolved independently multiple times. Comparative analyses of the first transcriptome atlas for the main distyly model, Primula veris, with other distylous species produced the following findings. A set of 53 constitutively expressed genes in P. veris did not include any of the housekeeping genes commonly used to normalize gene expression in qPCR experiments. The S-locus gene CYP^T acquired its role in controlling style elongation via a change in expression profile. Comparison of genes differentially expressed between floral morphs revealed that brassinosteroids and auxin are the main hormones controlling style elongation in P. veris and Fagopyrum esculentum, respectively. Furthermore, shared biochemical pathways might underlie the expression of distyly in the distantly related P. veris, F. esculentum and Turnera subulata, suggesting a degree of correspondence between evolutionary convergence at phenotypic and molecular levels. Finally, we provide the first evidence supporting the previously proposed hypothesis that distyly supergenes of distantly related species evolved via the recruitment of genes related to the phytochrome-interacting factor (PIF) signaling network. To conclude, this is the first study that discovered homologous genes involved in the control of distyly in distantly related taxa.

The S-Gene YUC6 Pleiotropically Determines Male Mating Type and Pollen Size in Heterostylous Turnera (Passifloraceae): A Novel Neofunctionalization of the YUCCA Gene Family

In heterostylous, self-incompatible Turnera species, a member of the YUCCA gene family, YUC6, resides at the S-locus and has been hypothesized to determine the male mating type. YUCCA gene family members synthesize the auxin, indole-3-acetic acid, via a two-step process involving the TAA gene family. Consequently, it has been speculated that differences in auxin concentration in developing anthers are the biochemical basis underlying the male mating type. Here, we provide empirical evidence that supports this hypothesis. Using a transgenic knockdown approach, we show that YUC6 acts pleiotropically to control both the male physiological mating type and pollen size, but not the filament length dimorphism associated with heterostyly in Turnera. Using qPCR to assess YUC6 expression in different transgenic lines, we demonstrate that the level of YUC6 knockdown correlates with the degree of change observed in the male mating type. Further assessment of YUC6 expression through anther development, in the knockdown lines, suggests that the male mating type is irreversibly determined during a specific developmental window prior to microsporogenesis, which is consistent with the genetically sporophytic nature of this self-incompatibility system. These results represent the first gene controlling male mating type to be characterized in any species with heterostyly.

Phylogenomics, floral evolution, and biogeography of Lithospermum L. (Boraginaceae)

Lithospermum (Boraginaceae), a geographically cosmopolitan medium-sized genus, includes diverse floral morphology, with variation in corolla size and shape and in breeding system. Over the past decade, multiple studies have examined the evolutionary history of Lithospermum, with most utilizing DNA regions from the plastid genome and/or the nuclear ribosomal internal transcribed spacer. These studies have, in general, not resulted in well-resolved and well-supported phylogenies. In the present study, 298 nuclear DNA regions, amplified via target sequence capture, were utilized for phylogenetic reconstruction for Lithospermum and relatives in Boraginaceae, and patterns of floral evolution, species diversification, and biogeography were examined. Based on multiple phylogenetic methods, Lithospermum is resolved as monophyletic, and the New World species of the genus are also monophyletic. While minimal phylogenetic incongruence is resolved within the nuclear genome, incongruence between the nuclear and plastid genomes is recovered. This is likely due to incomplete lineage sorting during early diversification of the genus in the Americas approximately 7.8 million years ago. At least four shifts to longer corollas are identified throughout Lithospermum, and this may be due to selection for hummingbird-pollinated flowers, particularly for species in Mexico and the southwestern United States. In the New World, one clade of species of the genus diversified primarily across the United States and Canada, and another radiated throughout the mountains of Mexico.

Differential Gene Expression with an Emphasis on Floral Organ Size Differences in Natural and Synthetic Polyploids of Nicotiana tabacum (Solanaceae)

Floral organ size, especially the size of the corolla, plays an important role in plant reproduction by facilitating pollination efficiency. Previous studies have outlined a hypothesized organ size pathway. However, the expression and function of many of the genes in the pathway have only been investigated in model diploid species; therefore, it is unknown how these genes interact in polyploid species. Although correlations between ploidy and cell size have been shown in many systems, it is unclear whether there is a difference in cell size between naturally occurring and synthetic polyploids. To address these questions comparing floral organ size and cell size across ploidy, we use natural and synthetic polyploids of Nicotiana tabacum (Solanaceae) as well as their known diploid progenitors. We employ a comparative transcriptomics approach to perform analyses of differential gene expression, focusing on candidate genes that may be involved in floral organ size, both across developmental stages and across accessions. We see differential expression of several known floral organ candidate genes including ARF2, BIG BROTHER, and GASA/GAST1. Results from linear models show that ploidy, cell width, and cell number positively influence corolla tube circumference; however, the effect of cell width varies by ploidy, and diploids have a significantly steeper slope than both natural and synthetic polyploids. These results demonstrate that polyploids have wider cells and that polyploidy significantly increases corolla tube circumference.

Transcriptome and Network Analyses of Heterostyly in Turnera subulata Provide Mechanistic Insights: Are S-Loci a Red-Light for Pistil Elongation?

Heterostyly employs distinct hermaphroditic floral morphs to enforce outbreeding. Morphs differ structurally in stigma/anther positioning, promoting cross-pollination, and physiologically blocking self-fertilization. Heterostyly is controlled by a self-incompatibility (S)-locus of a small number of linked S-genes specific to short-styled morph genomes. Turnera possesses three S-genes, namely TsBAHD (controlling pistil characters), TsYUC6, and TsSPH1 (controlling stamen characters). Here, we compare pistil and stamen transcriptomes of floral morphs of T. subulata to investigate hypothesized S-gene function(s) and whether hormonal differences might contribute to physiological incompatibility. We then use network analyses to identify genetic networks underpinning heterostyly. We found a depletion of brassinosteroid-regulated genes in short styled (S)-morph pistils, consistent with hypothesized brassinosteroid-inactivating activity of TsBAHD. In S-morph anthers, auxin-regulated genes were enriched, consistent with hypothesized auxin biosynthesis activity of TsYUC6. Evidence was found for auxin elevation and brassinosteroid reduction in both pistils and stamens of S- relative to long styled (L)-morph flowers, consistent with reciprocal hormonal differences contributing to physiological incompatibility. Additional hormone pathways were also affected, however, suggesting S-gene activities intersect with a signaling hub. Interestingly, distinct S-genes controlling pistil length, from three species with independently evolved heterostyly, potentially intersect with phytochrome interacting factor (PIF) network hubs which mediate red/far-red light signaling. We propose that modification of the activities of PIF hubs by the S-locus could be a common theme in the evolution of heterostyly.

How to build distylous flowers: comparative floral development and evolution of distylous species across the angiosperms

PREMISE

Distyly, a plant breeding system characterized by two floral morphs that have reciprocal positioning of anthers and stigmas, is known from at least 27 angiosperm families, making it an excellent example of convergent evolution. The various manners in which patterns of floral development produce the distinct anther and stigma heights in each morph remain largely unexplored from developmental and evolutionary perspectives.

METHODS

In 15 species representing at least 12 origins of distyly, heights and lengths of floral organs in each morph throughout development were examined using light microscopy. Patterns of floral organ development were determined and compared among species. Family-level phylogenies of distylous species and relatives were reconstructed, and patterns of ancestral herkogamy were resolved.

RESULTS

Differences in floral development between morphs resulted in 12 patterns leading to the anther and stigma positions characterizing distyly. Distylous species evolved from ancestors with different types of herkogamy, with approach herkogamy and lack of herkogamy resolved most frequently.

CONCLUSIONS

Seven of the 12 patterns of floral development are known from only one species, with three other patterns described among pairs of close relatives. The most common pattern of floral development, described from at least seven genera, involves for anther heights, distinct intermorph growth rates and for stigma heights, growth rates that differ between morphs only during later development. This pattern is common among subclass Lamiidae, suggesting canalized development within the taxon. Among distylous species, the same type of ancestral herkogamy can give rise to different patterns of floral development.

Global transcriptome and gene co-expression network analyses on the development of distyly in Primula oreodoxa

Distyly is a genetically controlled flower polymorphism that has intrigued both botanists and evolutionary biologists ever since Darwin's time. Despite extensive reports on the pollination and evolution of distylous systems, the genetic basis and mechanism of molecular regulation remain unclear. In the present study, comparative transcriptome profiling was conducted in primrose (Primula oreodoxa), the prime research model for heterostyly. Thirty-six transcriptomes were sequenced for styles at different stages and corolla tube in the three morphs of P. oreodoxa. Large numbers of differentially expressed genes (DEGs) were detected in the transcriptomes of styles across different morphs. Several transcription factors (TFs) and phytohormone metabolism-related genes were highlighted in S-morphs. A growing number of genes showed differential expression patterns along with the development of styles, suggesting that the genetic control of distyly may be more complicated than ever expected. Analysis of co-expression networks and module-trait relationships identified modules significantly associated with style development. CYP734A50, a key S-locus gene whose products degrade brassinosteroids, was co-expressed with many genes in the module and showed significant negative association with style length. In addition, crucial TFs involved in phytohormone signaling pathways were found to be connected with CYP734A50 in the co-expression module. Our global transcriptomic analysis has identified DEGs that are potentially involved in regulation of style length in P. oreodoxa, and may shed light on the evolution and broad biological processes of heterostyly.

Unravelling a stearidonic acid-rich triacylglycerol biosynthetic pathway in the developing seeds of Buglossoides arvensis: A transcriptomic landscape

Buglossoides arvensis is an emerging oilseed crop that is rich in stearidonic acid (SDA) and has several potential applications in human health and nutrition. The molecular basis of SDA biosynthesis in this plant remains unknown due to lack of genomic information. To unravel key genes involved in SDA-rich triacylglycerol (TAG) biosynthesis, we performed transcriptome sequencing of pooled mRNA from five different developmental stages of B. arvensis seeds using Illumina NextSeq platform. De novo transcriptome assembly generated 102,888 clustered transcripts from 39.83 million high-quality reads. Of these, 62.1% and 55.54% of transcripts were functionally annotated using Uniprot-Viridiplantae and KOG databases, respectively. A total of 10,021 SSR-containing sequences were identified using the MISA tool. Deep mining of transcriptome assembly using in silico tools led to the identification of genes involved in fatty acid and TAG biosynthesis. Expression profiling of 17 key transcripts involved in fatty acid desaturation and TAG biosynthesis showed expression patterns specific to the development stage that positively correlated with polyunsaturated fatty acid accumulation in the developing seeds. This first comprehensive transcriptome analysis provides the basis for future research on understanding molecular mechanisms of SDA-rich TAG accumulation in B. arvensis and aids in biotechnological production of SDA in other oilseed crops.