Gene Duplication, Shifting Selection, and Dosage Balance of Silicon Transporter Proteins in Marine and Freshwater Diatoms.

Numerous factors shape the evolution of protein-coding genes, including shifts in the strength or type of selection following gene duplications or changes in the environment. Diatoms and other silicifying organisms use a family of silicon transporters (SITs) to import dissolved silicon from the environment. Freshwaters contain higher silicon levels than oceans, and marine diatoms have more efficient uptake kinetics and less silicon in their cell walls, making them better competitors for a scarce resource. We compiled SITs from 37 diatom genomes to characterize shifts in selection following gene duplications and marine-freshwater transitions. A deep gene duplication, which coincided with a whole-genome duplication, gave rise to two gene lineages. One of them (SIT1-2) is present in multiple copies in most species and is known to actively import silicon. These SITs have evolved under strong purifying selection that was relaxed in freshwater taxa. Episodic diversifying selection was detected but not associated with gene duplications or habitat shifts. In contrast, genes in the second SIT lineage (SIT3) were present in just half the species, the result of multiple losses. Despite conservation of SIT3 in some lineages for the past 90-100 million years, repeated losses, relaxed selection, and low expression highlighted the dispensability of SIT3, consistent with a model of deterioration and eventual loss due to relaxed selection on SIT3 expression. The extensive but relatively balanced history of duplications and losses, together with paralog-specific expression patterns, suggest diatoms continuously balance gene dosage and expression dynamics to optimize silicon transport across major environmental gradients.

Resolving Marine-Freshwater Transitions by Diatoms Through a Fog of Gene Tree Discordance.

Despite the obstacles facing marine colonists, most lineages of aquatic organisms have colonized and diversified in freshwaters repeatedly. These transitions can trigger rapid morphological or physiological change and, on longer timescales, lead to increased rates of speciation and extinction. Diatoms are a lineage of ancestrally marine microalgae that have diversified throughout freshwater habitats worldwide. We generated a phylogenomic data set of genomes and transcriptomes for 59 diatom taxa to resolve freshwater transitions in one lineage, the Thalassiosirales. Although most parts of the species tree were consistently resolved with strong support, we had difficulties resolving a Paleocene radiation, which affected the placement of one freshwater lineage. This and other parts of the tree were characterized by high levels of gene tree discordance caused by incomplete lineage sorting and low phylogenetic signal. Despite differences in species trees inferred from concatenation versus summary methods and codons versus amino acids, traditional methods of ancestral state reconstruction supported six transitions into freshwaters, two of which led to subsequent species diversification. Evidence from gene trees, protein alignments, and diatom life history together suggest that habitat transitions were largely the product of homoplasy rather than hemiplasy, a condition where transitions occur on branches in gene trees not shared with the species tree. Nevertheless, we identified a set of putatively hemiplasious genes, many of which have been associated with shifts to low salinity, indicating that hemiplasy played a small but potentially important role in freshwater adaptation. Accounting for differences in evolutionary outcomes, in which some taxa became locked into freshwaters while others were able to return to the ocean or become salinity generalists, might help further distinguish different sources of adaptive mutation in freshwater diatoms.

Phylotranscriptomics reveals the reticulate evolutionary history of a widespread diatom species complex.

In contrast to surveys based on a few genes that often provide limited taxonomic resolution, transcriptomes provide a wealth of genomic loci that can resolve relationships among taxonomically challenging lineages. Diatoms are a diverse group of aquatic microalgae that includes important bioindicator species and many such lineages. One example is Nitzschia palea, a widespread species complex with several morphologically defined taxonomic varieties, some of which are critical pollution indicators. Morphological differences among the varieties are subtle and phylogenetic studies based on a few genes fail to resolve their evolutionary relationships. We conducted morphometric and transcriptome analyses of 10 Nitzschia palea strains to resolve the relationships among strains and taxonomic varieties. Nitzschia palea was resolved into three clades, one of which corresponds to a group of strains with narrow linear-lanceolate valves. The other morphological group recovered in the shape outline analysis was not monophyletic and consisted of two clades. Gene-tree concordance analyses and phylogenetic network estimations revealed patterns of incomplete lineage sorting and gene flow between intraspecific lineages. We detected reticulated evolutionary patterns among lineages with different morphologies, resulting in a putative recent hybrid. Our study shows that phylogenomic analyses of unlinked nuclear loci, complemented with morphometrics, can resolve complex evolutionary histories of recently diverged species complexes.

The genome of a nonphotosynthetic diatom provides insights into the metabolic shift to heterotrophy and constraints on the loss of photosynthesis.

Although most of the tens of thousands of diatom species are photoautotrophs, a small number of heterotrophic species no longer photosynthesize. We sequenced the genome of a nonphotosynthetic diatom, Nitzschia Nitz4, to determine how carbon metabolism was altered in the wake of this trophic shift. Nitzschia Nitz4 has retained its plastid and plastid genome, but changes associated with the transition to heterotrophy were cellular-wide and included losses of photosynthesis-related genes from the nuclear and plastid genomes, elimination of isoprenoid biosynthesis in the plastid, and remodeling of mitochondrial glycolysis to maximize adenosine triphosphte (ATP) yield. The genome contains a ß-ketoadipate pathway that may allow Nitzschia Nitz4 to metabolize lignin-derived compounds. Diatom plastids lack an oxidative pentose phosphate pathway (oPPP), leaving photosynthesis as the primary source of NADPH to support essential biosynthetic pathways in the plastid and, by extension, limiting available sources of NADPH in nonphotosynthetic plastids. The genome revealed similarities between nonphotosynthetic diatoms and apicomplexan parasites for provisioning NADPH in their plastids and highlighted the ancestral absence of a plastid oPPP as a potentially important constraint on loss of photosynthesis, a hypothesis supported by the higher frequency of transitions to parasitism or heterotrophy in lineages that have a plastid oPPP.

Genome Sequences of Bacteria Associated with the Diatom Cyclotella cryptica Strain CCMP332.

We report draft genomes of two bacterial strains in the genera Hyphobacterium and Reichenbachiella, which are associated with the diatom Cyclotella cryptica strain CCMP332. Genomes from these strains were 2,691,501 bp and 3,325,829 bp, respectively, and will be useful for understanding interactions between diatoms and bacteria.

Improved Reference Genome for Cyclotella cryptica CCMP332, a Model for Cell Wall Morphogenesis, Salinity Adaptation, and Lipid Production in Diatoms (Bacillariophyta).

The diatom, Cyclotella cryptica, is a well-established model species for physiological studies and biotechnology applications of diatoms. To further facilitate its use as a model diatom, we report an improved reference genome assembly and annotation for C. cryptica strain CCMP332. We used a combination of long- and short-read sequencing to assemble a high-quality and contaminant-free genome. The genome is 171 Mb in size and consists of 662 scaffolds with a scaffold N50 of 494 kb. This represents a 176-fold decrease in scaffold number and 41-fold increase in scaffold N50 compared to the previous assembly. The genome contains 21,250 predicted genes, 75% of which were assigned putative functions. Repetitive DNA comprises 59% of the genome, and an improved classification of repetitive elements indicated that a historically steady accumulation of transposable elements has contributed to the relatively large size of the C. cryptica genome. The high-quality C. cryptica genome will serve as a valuable reference for ecological, genetic, and biotechnology studies of diatoms.

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.

Phylogenomic analyses reveal extensive gene flow within the magic flowers (Achimenes).

PREMISE OF THE STUDY: The Neotropical Gesneriaceae is a lineage known for its colorful and diverse flowers, as well as an extensive history of intra- and intergeneric hybridization, particularly among Achimenes (the magic flowers) and other members of subtribe Gloxiniinae. Despite numerous studies seeking to elucidate the evolutionary relationships of these lineages, relatively few have sought to infer specific patterns of gene flow despite evidence of widespread hybridization. METHODS: To explore the utility of phylogenomic data for reassessing phylogenetic relationships and inferring patterns of gene flow among species of Achimenes, we sequenced 12 transcriptomes. We used a variety of methods to infer the species tree, examine gene tree discordance, and infer patterns of gene flow. KEY RESULTS: Phylogenomic analyses resolve clade relationships at the crown of the lineage with strong support. In contrast to previous analyses, we recovered strong support for several new relationships despite a significant amount of gene tree discordance. We present evidence for at least two introgression events between two species pairs that share pollinators, and suggest that the species status of Achimenes admirabilis be reexamined. CONCLUSIONS: Our study demonstrates the utility of transcriptome data for phylogenomic analyses, and inferring patterns of gene flow despite gene tree discordance. Moreover, these data provide another example of prevalent interspecific gene flow among Neotropical plants that share pollinators.

Comparative transcriptome analyses of flower development in four species of Achimenes (Gesneriaceae).

BACKGROUND: Flowers have an amazingly diverse display of colors and shapes, and these characteristics often vary significantly among closely related species. The evolution of diverse floral form can be thought of as an adaptive response to pollination and reproduction, but it can also be seen through the lens of morphological and developmental constraints. To explore these interactions, we use RNA-seq across species and development to investigate gene expression and sequence evolution as they relate to the evolution of the diverse flowers in a group of Neotropical plants native to Mexico-magic flowers (Achimenes, Gesneriaceae). RESULTS: The assembled transcriptomes contain between 29,000 and 42,000 genes expressed during development. We combine sequence orthology and coexpression clustering with analyses of protein evolution to identify candidate genes for roles in floral form evolution. Over 25% of transcripts captured were distinctive to Achimenes and overrepresented by genes involved in transcription factor activity. Using a model-based clustering approach we find dynamic, temporal patterns of gene expression among species. Selection tests provide evidence of positive selection in several genes with roles in pigment production, flowering time, and morphology. Combining these approaches to explore genes related to flower color and flower shape, we find distinct patterns that correspond to transitions of floral form among Achimenes species. CONCLUSIONS: The floral transcriptomes developed from four species of Achimenes provide insight into the mechanisms involved in the evolution of diverse floral form among closely related species with different pollinators. We identified several candidate genes that will serve as an important and useful resource for future research. High conservation of sequence structure, patterns of gene coexpression, and detection of positive selection acting on few genes suggests that large phenotypic differences in floral form may be caused by genetic differences in a small set of genes. Our characterized floral transcriptomes provided here should facilitate further analyses into the genomics of flower development and the mechanisms underlying the evolution of diverse flowers in Achimenes and other Neotropical Gesneriaceae.

Distinct Processes Drive Diversification in Different Clades of Gesneriaceae.

Using a time-calibrated phylogenetic hypothesis including 768 Gesneriaceae species (out of [Formula: see text]3300 species) and more than 29,000 aligned bases from 26 gene regions, we test Gesneriaceae for diversification rate shifts and the possible proximal drivers of these shifts: geographic distributions, growth forms, and pollination syndromes. Bayesian Analysis of Macroevolutionary Mixtures analyses found five significant rate shifts in Beslerieae, core Nematanthus, core Columneinae, core Streptocarpus, and Pacific Cyrtandra These rate shifts correspond with shifts in diversification rates, as inferred by Binary State Speciation and Extinction Model and Geographic State Speciation and Extinction model, associated with hummingbird pollination, epiphytism, unifoliate growth, and geographic area. Our results suggest that diversification processes are extremely variable across Gesneriaceae clades with different combinations of characters influencing diversification rates in different clades. Diversification patterns between New and Old World lineages show dramatic differences, suggesting that the processes of diversification in Gesneriaceae are very different in these two geographic regions.