RESUMO
Land plants evolved from charophytic algae, among which Charophyceae possess the most complex body plans. We present the genome of Chara braunii; comparison of the genome to those of land plants identified evolutionary novelties for plant terrestrialization and land plant heritage genes. C. braunii employs unique xylan synthases for cell wall biosynthesis, a phragmoplast (cell separation) mechanism similar to that of land plants, and many phytohormones. C. braunii plastids are controlled via land-plant-like retrograde signaling, and transcriptional regulation is more elaborate than in other algae. The morphological complexity of this organism may result from expanded gene families, with three cases of particular note: genes effecting tolerance to reactive oxygen species (ROS), LysM receptor-like kinases, and transcription factors (TFs). Transcriptomic analysis of sexual reproductive structures reveals intricate control by TFs, activity of the ROS gene network, and the ancestral use of plant-like storage and stress protection proteins in the zygote.
Assuntos
Chara/genética , Genoma de Planta , Evolução Biológica , Parede Celular/metabolismo , Chara/crescimento & desenvolvimento , Embriófitas/genética , Redes Reguladoras de Genes , Pentosiltransferases/genética , Reguladores de Crescimento de Plantas/metabolismo , Proteínas de Plantas/genética , Proteínas de Plantas/metabolismo , Proteínas Serina-Treonina Quinases/genética , Proteínas Serina-Treonina Quinases/metabolismo , Espécies Reativas de Oxigênio/metabolismo , Fatores de Transcrição/genética , Fatores de Transcrição/metabolismo , TranscriptomaRESUMO
The genome of the liverwort Marchantia polymorpha is an important step toward development of a new plant model system (Bowman et al., 2017). Liverworts may be the sister taxon to all other land plants, and the genome shows features that illuminate the ancestor of all land plants and give insights into how plant systems function and evolved.
Assuntos
Embriófitas , Marchantia/genética , PlantasRESUMO
Green plants, broadly defined as green algae and the land plants (together, Viridiplantae), constitute the primary eukaryotic lineage that successfully colonized Earth's emergent landscape. Members of various clades of green plants have independently made the transition from fully aquatic to subaerial habitats many times throughout Earth's history. The transition, from unicells or simple filaments to complex multicellular plant bodies with functionally differentiated tissues and organs, was accompanied by innovations built upon a genetic and phenotypic toolkit that have served aquatic green phototrophs successfully for at least a billion years. These innovations opened an enormous array of new, drier places to live on the planet and resulted in a huge diversity of land plants that have dominated terrestrial ecosystems over the past 500 million years. This review examines the greening of the land from several perspectives, from paleontology to phylogenomics, to water stress responses and the genetic toolkit shared by green algae and plants, to the genomic evolution of the sporophyte generation. We summarize advances on disparate fronts in elucidating this important event in the evolution of the biosphere and the lacunae in our understanding of it. We present the process not as a step-by-step advancement from primitive green cells to an inevitable success of embryophytes, but rather as a process of adaptations and exaptations that allowed multiple clades of green plants, with various combinations of morphological and physiological terrestrialized traits, to become diverse and successful inhabitants of the land habitats of Earth.
Assuntos
Clorófitas , Embriófitas , Evolução Biológica , Ecossistema , Embriófitas/genética , Filogenia , Plantas/genética , Clorófitas/genética , Evolução MolecularRESUMO
The Neoproterozoic Era records the transition from a largely bacterial to a predominantly eukaryotic phototrophic world, creating the foundation for the complex benthic ecosystems that have sustained Metazoa from the Ediacaran Period onward. This study focuses on the evolutionary origins of green seaweeds, which play an important ecological role in the benthos of modern sunlit oceans and likely played a crucial part in the evolution of early animals by structuring benthic habitats and providing novel niches. By applying a phylogenomic approach, we resolve deep relationships of the core Chlorophyta (Ulvophyceae or green seaweeds, and freshwater or terrestrial Chlorophyceae and Trebouxiophyceae) and unveil a rapid radiation of Chlorophyceae and the principal lineages of the Ulvophyceae late in the Neoproterozoic Era. Our time-calibrated tree points to an origin and early diversification of green seaweeds in the late Tonian and Cryogenian periods, an interval marked by two global glaciations with strong consequent changes in the amount of available marine benthic habitat. We hypothesize that unicellular and simple multicellular ancestors of green seaweeds survived these extreme climate events in isolated refugia, and diversified in benthic environments that became increasingly available as ice retreated. An increased supply of nutrients and biotic interactions, such as grazing pressure, likely triggered the independent evolution of macroscopic growth via different strategies, including true multicellularity, and multiple types of giant-celled forms.
Assuntos
Clorófitas/crescimento & desenvolvimento , Evolução Molecular , Alga Marinha/crescimento & desenvolvimento , Clorófitas/classificação , Ecossistema , Filogenia , Alga Marinha/classificaçãoRESUMO
SARS-CoV-2 is a member of the subgenus Sarbecovirus and thus contains the genetic element s2m. We have extensively mined nucleotide data in GenBank in order to obtain a comprehensive list of s2m sequences both in the four virus families where s2m has previously been described and in other groups of organisms. Surprisingly, there seems to be a xenologue of s2m in a large number of insect species. The function of s2m is unknown, but our data show a very high degree of sequence conservation both in insects and in viruses and that the version of s2m found in SARS-CoV-2 has unique features, not seen in any other virus or insect strains.
Assuntos
Transferência Genética Horizontal , Genoma Viral , Insetos/genética , SARS-CoV-2/genética , Animais , Mineração de Dados , Filogenia , RNA Viral/genéticaRESUMO
The primarily freshwater genus Chara is comprised of many species that exhibit a wide range of salinity tolerance. The range of salt tolerance provides a good platform for investigating the role of transport mechanisms in response to salt stress, and the close evolutionary relationship between Charophytes and land plants can provide broader insights. We investigated the response to salt stress of previously identified transport mechanisms in two species of Chara, Chara longifolia (salt-tolerant), and Chara australis (salt-sensitive): a cation transporter (HKT), a Na+ /H+ antiport (NHX), H+ -ATPase (AHA), and a Na+ -ATPase (ENA). The presence of these candidate genes has been confirmed in both species of Chara, with the exception of the Na+ -ATPase, which is present only in salt-tolerant Chara longifolia. Time-course Illumina transcriptomes were created using RNA from multiple time points (0, 6, 12, 24 and 48 h) after freshwater cultures for each species were exposed to salt stress. These transcriptomes verified our hypotheses of these mechanisms conferring salt tolerance in the two species examined and also aided in identification of specific transcripts representing our genes of interest in both species. The expression of these transcripts was validated through use of qPCR, in a similar experimental set-up used for the RNAseq data described above. The RNAseq and qPCR data showed significant changes of expression mechanisms in C. longifolia (respectively), a down-regulation of HKT and a substantial up-regulation of ENA. Significant responses to salt stress in salt-sensitive C. australis show up-regulation of NHX and AHA.
Assuntos
Chara , Salinidade , Adenosina Trifosfatases , Expressão Gênica , Tolerância ao Sal/genéticaRESUMO
Species within the genus Chara have variable salinity tolerance. Their close evolutionary relationship with embryophytes makes their study crucial to understanding the evolution of salt tolerance and key evolutionary processes shared among the phyla. We examined salt-tolerant Chara longifolia and salt-sensitive Chara australis for mechanisms of salt tolerance and their potential role in adaptation to salt. We hypothesize that there are shared mechanisms similar to those in embryophytes, which assist in conferring salt tolerance in Chara, including a cation transporter (HKT), a Na+ /H+ antiport (NHX), a H+ -ATPase (AHA), and a Na+ -ATPase (ENA). Illumina transcriptomes were created using cultures grown in freshwater and exposed to salt stress. The presence of these candidate genes, identified by comparing with genes known from embryophytes, has been confirmed in both species of Chara, with the exception of ENA, present only in salt-tolerant C. longifolia. These transcriptomes provide evidence for the contribution of these mechanisms to differences in salt tolerance in the two species and for the independent evolution of the Na+ -ATPase. We also examined genes that may have played a role in important evolutionary processes, suggested by previous work on the Chara braunii genome. Among the genes examined, cellulose synthase protein (GT43) and response regulator (RRB) were confirmed in both species. Genes absent from all three Chara species were members of the GRAS family, microtubule-binding protein (TANGLED1), and auxin synthesizers (YUCCA, TAA). Results from this study shed light on the evolutionary relationship between Chara and embryophytes through confirmation of shared salt tolerance mechanisms, as well as unique mechanisms that do not occur in angiosperms.
Assuntos
Chara , Carofíceas , Proteínas de Transporte , Transporte de Íons , Tolerância ao SalRESUMO
Despite their diversity and ecological importance, many areas of the SAR-Stramenopila, Alveolata, and Rhizaria-clade are poorly understood as the majority (90%) of SAR species lack molecular data and only 5% of species are from well-sampled families. Here, we review and summarize the state of knowledge about the three major clades of SAR, describing the diversity within each clade and identifying synapomorphies when possible. We also assess the "dark area" of SAR: the morphologically described species that are missing molecular data. The majority of molecular data for SAR lineages are characterized from marine samples and vertebrate hosts, highlighting the need for additional research effort in areas such as freshwater and terrestrial habitats and "non-vertebrate" hosts. We also describe the paucity of data on the biogeography of SAR species, and point to opportunities to illuminate diversity in this major eukaryotic clade. See also the video abstract here: https://youtu.be/_VUXqaX19Rw.
Assuntos
Alveolados/fisiologia , Escuridão , Rhizaria/fisiologia , Estramenópilas/fisiologia , Ecologia , Ecossistema , Eucariotos/fisiologia , Filogenia , Análise de Sequência de DNARESUMO
Dinoflagellates are key species in marine environments, but they remain poorly understood in part because of their large, complex genomes, unique molecular biology, and unresolved in-group relationships. We created a taxonomically representative dataset of dinoflagellate transcriptomes and used this to infer a strongly supported phylogeny to map major morphological and molecular transitions in dinoflagellate evolution. Our results show an early-branching position of Noctiluca, monophyly of thecate (plate-bearing) dinoflagellates, and paraphyly of athecate ones. This represents unambiguous phylogenetic evidence for a single origin of the group's cellulosic theca, which we show coincided with a radiation of cellulases implicated in cell division. By integrating dinoflagellate molecular, fossil, and biogeochemical evidence, we propose a revised model for the evolution of thecal tabulations and suggest that the late acquisition of dinosterol in the group is inconsistent with dinoflagellates being the source of this biomarker in pre-Mesozoic strata. Three distantly related, fundamentally nonphotosynthetic dinoflagellates, Noctiluca, Oxyrrhis, and Dinophysis, contain cryptic plastidial metabolisms and lack alternative cytosolic pathways, suggesting that all free-living dinoflagellates are metabolically dependent on plastids. This finding led us to propose general mechanisms of dependency on plastid organelles in eukaryotes that have lost photosynthesis; it also suggests that the evolutionary origin of bioluminescence in nonphotosynthetic dinoflagellates may be linked to plastidic tetrapyrrole biosynthesis. Finally, we use our phylogenetic framework to show that dinoflagellate nuclei have recruited DNA-binding proteins in three distinct evolutionary waves, which included two independent acquisitions of bacterial histone-like proteins.
Assuntos
Dinoflagellida/genética , Evolução Molecular , Filogenia , Plastídeos , RNA de Protozoário/genética , Análise de Sequência de RNA , TranscriptomaRESUMO
Coccolithophores have influenced the global climate for over 200 million years. These marine phytoplankton can account for 20 per cent of total carbon fixation in some systems. They form blooms that can occupy hundreds of thousands of square kilometres and are distinguished by their elegantly sculpted calcium carbonate exoskeletons (coccoliths), rendering them visible from space. Although coccolithophores export carbon in the form of organic matter and calcite to the sea floor, they also release CO2 in the calcification process. Hence, they have a complex influence on the carbon cycle, driving either CO2 production or uptake, sequestration and export to the deep ocean. Here we report the first haptophyte reference genome, from the coccolithophore Emiliania huxleyi strain CCMP1516, and sequences from 13 additional isolates. Our analyses reveal a pan genome (core genes plus genes distributed variably between strains) probably supported by an atypical complement of repetitive sequence in the genome. Comparisons across strains demonstrate that E. huxleyi, which has long been considered a single species, harbours extensive genome variability reflected in different metabolic repertoires. Genome variability within this species complex seems to underpin its capacity both to thrive in habitats ranging from the equator to the subarctic and to form large-scale episodic blooms under a wide variety of environmental conditions.
Assuntos
Genoma/genética , Haptófitas/genética , Haptófitas/isolamento & purificação , Fitoplâncton/genética , Calcificação Fisiológica , Cálcio/metabolismo , Anidrases Carbônicas/genética , Anidrases Carbônicas/metabolismo , Ecossistema , Haptófitas/classificação , Haptófitas/metabolismo , Oceanos e Mares , Filogenia , Proteoma/genética , Água do MarRESUMO
It is well known that ethylene regulates a diverse set of developmental and stress-related processes in angiosperms, yet its roles in early-diverging embryophytes and algae are poorly understood. Recently, it was shown that ethylene functions as a hormone in the charophyte green alga Spirogyra pratensis Since land plants evolved from charophytes, this implies conservation of ethylene as a hormone in green plants for at least 450 million years. However, the physiological role of ethylene in charophyte algae has remained unknown. To gain insight into ethylene responses in Spirogyra, we used mRNA sequencing to measure changes in gene expression over time in Spirogyra filaments in response to an ethylene treatment. Our analyses show that at the transcriptional level, ethylene predominantly regulates three processes in Spirogyra: (1) modification of the cell wall matrix by expansins and xyloglucan endotransglucosylases/hydrolases, (2) down-regulation of chlorophyll biosynthesis and photosynthesis, and (3) activation of abiotic stress responses. We confirmed that the photosynthetic capacity and chlorophyll content were reduced by an ethylene treatment and that several abiotic stress conditions could stimulate cell elongation in an ethylene-dependent manner. We also found that the Spirogyra transcriptome harbors only 10 ethylene-responsive transcription factor (ERF) homologs, several of which are regulated by ethylene. These results provide an initial understanding of the hormonal responses induced by ethylene in Spirogyra and help to reconstruct the role of ethylene in ancestral charophytes prior to the origin of land plants.
Assuntos
Parede Celular/efeitos dos fármacos , Etilenos/farmacologia , Perfilação da Expressão Gênica/métodos , Fotossíntese/efeitos dos fármacos , Spirogyra/efeitos dos fármacos , Proteínas de Algas/genética , Proteínas de Algas/metabolismo , Sequência de Aminoácidos , Parede Celular/genética , Parede Celular/metabolismo , Análise por Conglomerados , Expressão Gênica/efeitos dos fármacos , Ontologia Genética , Luz , Fotossíntese/genética , Reação em Cadeia da Polimerase Via Transcriptase Reversa , Homologia de Sequência de Aminoácidos , Cloreto de Sódio/farmacologia , Spirogyra/genética , Spirogyra/metabolismo , Estresse Fisiológico/efeitos dos fármacos , Estresse Fisiológico/genética , TemperaturaRESUMO
Isoprenoid metabolism occupies a central position in the anabolic metabolism of all living cells. In plastid-bearing organisms, two pathways may be present for de novo isoprenoid synthesis, the cytosolic mevalonate pathway (MVA) and nuclear-encoded, plastid-targeted nonmevalonate pathway (DOXP). Using transcriptomic data we find that dinoflagellates apparently make exclusive use of the DOXP pathway. Using phylogenetic analyses of all DOXP genes we inferred the evolutionary origins of DOXP genes in dinoflagellates. Plastid replacements led to a DOXP pathway of multiple evolutionary origins. Dinoflagellates commonly referred to as dinotoms due to their relatively recent acquisition of a diatom plastid, express two completely redundant DOXP pathways. Dinoflagellates with a tertiary plastid of haptophyte origin, by contrast, express a hybrid pathway of dual evolutionary origin. Here, changes in the targeting motif of signal/transit peptide likely allow for targeting the new plastid by the proteins of core isoprenoid metabolism proteins. Parasitic dinoflagellates of the Amoebophyra species complex appear to have lost the DOXP pathway, suggesting that they may rely on their host for sterol synthesis.
Assuntos
Dinoflagellida/genética , Dinoflagellida/metabolismo , Redes e Vias Metabólicas/genética , Terpenos/metabolismo , Perfilação da Expressão Gênica , Ácido Mevalônico/metabolismo , Filogenia , Plastídeos , Esteróis/biossínteseRESUMO
BACKGROUND: Clustering protein sequences according to inferred homology is a fundamental step in the analysis of many large data sets. Since the publication of the Markov Clustering (MCL) algorithm in 2002, it has been the centerpiece of several popular applications. Each of these approaches generates an undirected graph that represents sequences as nodes connected to each other by edges weighted with a BLAST-based metric. MCL is then used to infer clusters of homologous proteins by analyzing these graphs. The various approaches differ only by how they weight the edges, yet there has been very little direct examination of the relative performance of alternative edge-weighting metrics. This study compares the performance of four BLAST-based edge-weighting metrics: the bit score, bit score ratio (BSR), bit score over anchored length (BAL), and negative common log of the expectation value (NLE). Performance is tested using the Extended CEGMA KOGs (ECK) database, which we introduce here. RESULTS: All metrics performed similarly when analyzing full-length sequences, but dramatic differences emerged as progressively larger fractions of the test sequences were split into fragments. The BSR and BAL successfully rescued subsets of clusters by strengthening certain types of alignments between fragmented sequences, but also shifted the largest correct scores down near the range of scores generated from spurious alignments. This penalty outweighed the benefits in most test cases, and was greatly exacerbated by increasing the MCL inflation parameter, making these metrics less robust than the bit score or the more popular NLE. Notably, the bit score performed as well or better than the other three metrics in all scenarios. CONCLUSIONS: The results provide a strong case for use of the bit score, which appears to offer equivalent or superior performance to the more popular NLE. The insight that MCL-based clustering methods can be improved using a more tractable edge-weighting metric will greatly simplify future implementations. We demonstrate this with our own minimalist Python implementation: Porthos, which uses only standard libraries and can process a graph with 25 m + edges connecting the 60 k + KOG sequences in half a minute using less than half a gigabyte of memory.
Assuntos
Algoritmos , Análise por Conglomerados , Biologia Computacional/métodos , Cadeias de Markov , Proteínas/química , Alinhamento de Sequência/métodos , Análise de Sequência de Proteína/métodos , Sequência de Aminoácidos , Bases de Dados Factuais , Humanos , Dados de Sequência Molecular , Proteínas/metabolismo , Homologia de Sequência de Aminoácidos , SoftwareRESUMO
The competition-relatedness hypothesis (CRH) predicts that the strength of competition is the strongest among closely related species and decreases as species become less related. This hypothesis is based on the assumption that common ancestry causes close relatives to share biological traits that lead to greater ecological similarity. Although intuitively appealing, the extent to which phylogeny can predict competition and co-occurrence among species has only recently been rigorously tested, with mixed results. When studies have failed to support the CRH, critics have pointed out at least three limitations: (i) the use of data poor phylogenies that provide inaccurate estimates of species relatedness, (ii) the use of inappropriate statistical models that fail to detect relationships between relatedness and species interactions amidst nonlinearities and heteroskedastic variances, and (iii) overly simplified laboratory conditions that fail to allow eco-evolutionary relationships to emerge. Here, we address these limitations and find they do not explain why evolutionary relatedness fails to predict the strength of species interactions or probabilities of coexistence among freshwater green algae. First, we construct a new data-rich, transcriptome-based phylogeny of common freshwater green algae that are commonly cultured and used for laboratory experiments. Using this new phylogeny, we re-analyse ecological data from three previously published laboratory experiments. After accounting for the possibility of nonlinearities and heterogeneity of variances across levels of relatedness, we find no relationship between phylogenetic distance and ecological traits. In addition, we show that communities of North American green algae are randomly composed with respect to their evolutionary relationships in 99% of 1077 lakes spanning the continental United States. Together, these analyses result in one of the most comprehensive case studies of how evolutionary history influences species interactions and community assembly in both natural and experimental systems. Our results challenge the generality of the CRH and suggest it may be time to re-evaluate the validity and assumptions of this hypothesis.
Assuntos
Evolução Biológica , Clorófitas/fisiologia , Filogenia , Clorófitas/genética , Ecossistema , Dinâmica Populacional , Especificidade da Espécie , TranscriptomaRESUMO
DEFECTIVE KERNEL1 (DEK1) of higher plants plays an essential role in position-dependent signaling and consists of a large transmembrane domain (MEM) linked to a protease catalytic domain and a regulatory domain. Here, we show that the postulated sensory Loop of the MEM domain plays an important role in the developmental regulation of DEK1 activity in the moss Physcomitrella patens. Compared with P. patens lacking DEK1 (∆dek1), the dek1∆loop mutant correctly positions the division plane in the bud apical cell. In contrast with an early developmental arrest of ∆dek1 buds, dek1∆loop develops aberrant gametophores lacking expanded phyllids resulting from misregulation of mitotic activity. In contrast with the highly conserved sequence of the protease catalytic domain, the Loop is highly variable in land plants. Functionally, the sequence from Marchantia polymorpha fully complements the dek1∆loop phenotype, whereas sequences from maize (Zea mays) and Arabidopsis (Arabidopsis thaliana) give phenotypes with retarded growth and affected phyllid development. Bioinformatic analysis identifies MEM as a member of the Major Facilitator Superfamily, membrane transporters reacting to stimuli from the external environment. Transcriptome analysis comparing wild-type and ∆dek1 tissues identifies an effect on two groups of transcripts connected to dek1 mutant phenotypes: transcripts related to cell wall remodeling and regulation of the AINTEGUMENTA, PLETHORA, and BABY BOOM2 (APB2) and APB3 transcription factors known to regulate bud initiation. Finally, sequence data support the hypothesis that the advanced charophyte algae that evolved into ancestral land plants lost cytosolic calpains, retaining DEK1 as the sole calpain in the evolving land plant lineage.
Assuntos
Padronização Corporal , Bryopsida/genética , Genes de Plantas , Proteínas de Plantas/genética , Sequência de Aminoácidos , Dados de Sequência Molecular , Proteínas de Plantas/química , Proteínas de Plantas/fisiologia , Homologia de Sequência de AminoácidosRESUMO
Dinoflagellates are one of the last major lineages of eukaryotes for which little is known about genome structure and organization. We report here the sequence and gene structure of a clone isolated from a cosmid library which, to our knowledge, represents the largest contiguously sequenced, dinoflagellate genomic, tandem gene array. These data, combined with information from a large transcriptomic library, allowed a high level of confidence of every base pair call. This degree of confidence is not possible with PCR-based contigs. The sequence contains an intron-rich set of five highly expressed gene repeats arranged in tandem. One of the tandem repeat gene members contains an intron 26,372 bp long. This study characterizes a splice site consensus sequence for dinoflagellate introns. Two to nine base pairs around the 3' splice site are repeated by an identical two to nine base pairs around the 5' splice site. The 5' and 3' splice sites are in the same locations within each repeat so that the repeat is found only once in the mature mRNA. This identically repeated intron boundary sequence might be useful in gene modeling and annotation of genomes.
Assuntos
Dinoflagellida/genética , Genoma de Protozoário , Genômica/métodos , Íntrons , Splicing de RNA , Sequência de Aminoácidos , Sequência de Bases , Cosmídeos , Biblioteca Gênica , Genoma de Protozoário/genética , Dados de Sequência Molecular , Sequências de Repetição em TandemRESUMO
The alveolates are composed of three major lineages, the ciliates, dinoflagellates, and apicomplexans. Together these 'protist' taxa play key roles in primary production and ecology, as well as in illness of humans and other animals. The interface between the dinoflagellate and apicomplexan clades has been an area of recent discovery, blurring the distinction between these two clades. Moreover, phylogenetic analysis has yet to determine the position of basal dinoflagellate clades hence the deepest branches of the dinoflagellate tree currently remain unresolved. Large-scale mRNA sequencing was applied to 11 species of dinoflagellates, including strains of the syndinean genera Hematodinium and Amoebophrya, parasites of crustaceans and dinoflagellates, respectively, to optimize and update the dinoflagellate tree. From the transcriptome-scale data a total of 73 ribosomal protein-coding genes were selected for phylogeny. After individual gene orthology assessment, the genes were concatenated into a >15,000 amino acid alignment with 76 taxa from dinoflagellates, apicomplexans, ciliates, and the outgroup heterokonts. Overall the tree was well resolved and supported, when the data was subsampled with gblocks or constraint trees were tested with the approximately unbiased test. The deepest branches of the dinoflagellate tree can now be resolved with strong support, and provides a clearer view of the evolution of the distinctive traits of dinoflagellates.
Assuntos
Dinoflagellida/genética , Filogenia , Proteínas Ribossômicas/genética , Animais , Análise de Sequência de DNA , TranscriptomaRESUMO
Metagenomic methods provide a powerful means to investigate complex ecological phenomena. Developed originally for study of Bacteria and Archaea, the application of these methods to eukaryotic microorganisms is yet to be fully realized. Most prior environmental molecular studies of eukaryotes have relied heavily on PCR amplification with eukaryote-specific primers. Here we apply high throughput short-read sequencing of poly-A selected RNA to capture the metatranscriptome of an estuarine dinoflagellate bloom. To validate the metatranscriptome assembly process we simulated metatranscriptomic datasets using short-read sequencing data from clonal cultures of four algae of varying phylogenetic distance. We find that the proportion of chimeric transcripts reconstructed from community transcriptome sequencing is low, suggesting that metatranscriptomic sequencing can be used to accurately reconstruct the transcripts expressed by bloom-forming communities of eukaryotes. To further validate the bloom metatransciptome assembly we compared it to a transcriptomic assembly from a cultured, clonal isolate of the dominant bloom-causing alga and found that the two assemblies are highly similar. Eukaryote-wide phylogenetic analyses reveal the taxonomic composition of the bloom community, which is comprised of several dinoflagellates, ciliates, animals, and fungi. The assembled metatranscriptome reveals the functional genomic composition of a metabolically active community. Highlighting the potential power of these methods, we found that relative transcript abundance patterns suggest that the dominant dinoflagellate might be expressing toxin biosynthesis related genes at a higher level in the presence of competitors, predators and prey compared to it growing in monoculture.