Evidence for adaptive morphological plasticity in the Caribbean coral, Acropora cervicornis.

Genotype-by-environment interactions (GxE) indicate that variation in organismal traits cannot be explained by fixed effects of genetics or site-specific plastic responses alone. For tropical coral reefs experiencing dramatic environmental change, identifying the contributions of genotype, environment, and GxE on coral performance will be vital for both predicting persistence and developing restoration strategies. We quantified the impacts of G, E, and GxE on the morphology and survival of the endangered coral, Acropora cervicornis, through an in situ transplant experiment exposing common garden (nursery)-raised clones of ten genotypes to nine reef sites in the Florida Keys. By fate-tracking outplants over one year with colony-level 3D photogrammetry, we uncovered significant GxE on coral size, shape, and survivorship, indicating that no universal winner exists in terms of colony performance. Rather than differences in mean trait values, we found that individual-level morphological plasticity is adaptive in that the most plastic individuals also exhibited the fastest growth and highest survival. This indicates that adaptive morphological plasticity may continue to evolve, influencing the success of A. cervicornis and resulting reef communities in a changing climate. As focal reefs are active restoration sites, the knowledge that variation in phenotype is an important predictor of performance can be directly applied to restoration planning. Taken together, these results establish A. cervicornis as a system for studying the ecoevolutionary dynamics of phenotypic plasticity that also can inform genetic- and environment-based strategies for coral restoration.

Synergistic response to climate stressors in coral is associated with genotypic variation in baseline expression.

As environments are rapidly reshaped due to climate change, phenotypic plasticity plays an important role in the ability of organisms to persist and is considered an especially important acclimatization mechanism for long-lived sessile organisms such as reef-building corals. Often, this ability of a single genotype to display multiple phenotypes depending on the environment is modulated by changes in gene expression, which can vary in response to environmental changes via two mechanisms: baseline expression and expression plasticity. We used transcriptome-wide expression profiling of eleven genotypes of common-gardened Acropora cervicornis to explore genotypic variation in the expression response to thermal and acidification stress, both individually and in combination. We show that the combination of these two stressors elicits a synergistic gene expression response, and that both baseline expression and expression plasticity in response to stress show genotypic variation. Additionally, we demonstrate that frontloading of a large module of coexpressed genes is associated with greater retention of algal symbionts under combined stress. These results illustrate that variation in the gene expression response of individuals to climate change stressors can persist even when individuals have shared environmental histories, affecting their performance under future climate change scenarios.

Rebuild the Academy: Supporting academic mothers during COVID-19 and beyond.

The issues facing academic mothers have been discussed for decades. Coronavirus Disease 2019 (COVID-19) is further exposing these inequalities as womxn scientists who are parenting while also engaging in a combination of academic related duties are falling behind. These inequities can be solved by investing strategically in solutions. Here we describe strategies that would ensure a more equitable academy for working mothers now and in the future. While the data are clear that mothers are being disproportionately impacted by COVID-19, many groups could benefit from these strategies. Rather than rebuilding what we once knew, let us be the architects of a new world.

Performance of Orbicella faveolata larval cohorts does not align with previously observed thermal tolerance of adult source populations.

Orbicella faveolata, commonly known as the mountainous star coral, is a dominant reef-building species in the Caribbean, but populations have suffered sharp declines since the 1980s due to repeated bleaching and disease-driven mortality. Prior research has shown that inshore adult O. faveolata populations in the Florida Keys are able to maintain high coral cover and recover from bleaching faster than their offshore counterparts. However, whether this origin-specific variation in thermal resistance is heritable remains unclear. To address this knowledge gap, we produced purebred and hybrid larval crosses from O. faveolata gametes collected at two distinct reefs in the Upper Florida Keys, a nearshore site (Cheeca Rocks, CR) and an offshore site (Horseshoe Reef, HR), in two different years (2019, 2021). We then subjected these aposymbiotic larvae to severe (36°C) and moderate (32°C) heat challenges to quantify their thermal tolerance. Contrary to our expectation based on patterns of adult thermal tolerance, HR purebred larvae survived better and exhibited gene expression profiles that were less driven by stress response under elevated temperature compared to purebred CR and hybrid larvae. One potential explanation could be the compromised reproductive output of CR adult colonies due to repeated summer bleaching events in 2018 and 2019, as gametes originating from CR in 2019 contained less storage lipids than those from HR. These findings provide an important counter-example to the current selective breeding paradigm, that more tolerant parents will yield more tolerant offspring, and highlight the importance of adopting a holistic approach when evaluating larval quality for conservation and restoration purposes.

Heritable variation and lack of tradeoffs suggest adaptive capacity in Acropora cervicornis despite negative synergism under climate change scenarios.

Knowledge of multi-stressor interactions and the potential for tradeoffs among tolerance traits is essential for developing intervention strategies for the conservation and restoration of reef ecosystems in a changing climate. Thermal extremes and acidification are two major co-occurring stresses predicted to limit the recovery of vital Caribbean reef-building corals. Here, we conducted an aquarium-based experiment to quantify the effects of increased water temperatures and pCO2 individually and in concert on 12 genotypes of the endangered branching coral Acropora cervicornis, currently being reared and outplanted for large-scale coral restoration. Quantification of 12 host, symbiont and holobiont traits throughout the two-month-long experiment showed several synergistic negative effects, where the combined stress treatment often caused a greater reduction in physiological function than the individual stressors alone. However, we found significant genetic variation for most traits and positive trait correlations among treatments indicating an apparent lack of tradeoffs, suggesting that adaptive evolution will not be constrained. Our results suggest that it may be possible to incorporate climate-resistant coral genotypes into restoration and selective breeding programmes, potentially accelerating adaptation.

Global gene expression patterns in Porites white patch syndrome: Disentangling symbiont loss from the thermal stress response in reef-building coral.

The mechanisms resulting in the breakdown of the coral symbiosis once the process of bleaching has been initiated remain unclear. Distinguishing the process of symbiont loss from the thermal stress response may shed light on the cellular and molecular pathways involved in each process. This study examined physiological changes and global gene expression patterns associated with white patch syndrome (WPS) in Porites lobata, which manifests in localized bleaching independent of thermal stress. In addition, a meta-analysis of global gene expression studies in other corals and anemones was used to contrast differential regulation as a result of disease and thermal stress from patterns correlated with symbiotic state. Symbiont density, chlorophyll a content, holobiont productivity, instant calcification rate, and total host protein content were uniformly reduced in WPS relative to healthy tissue. While expression patterns associated with WPS were secondary to fixed effects of source colony, specific functional enrichments combined with a lack of immune regulation suggest that the viral infection putatively giving rise to this condition affects symbiont rather than host cells. Expression in response to WPS also clustered independently of patterns in white syndrome impacted A. hyacinthus, further supporting a distinct aetiology of this syndrome. Expression patterns in WPS-affected tissues were significantly correlated with prior studies that examined short-term thermal stress responses independent of symbiotic state, suggesting that the majority of expression changes reflect a nonspecific stress response. Across studies, the magnitude and direction of expression change among particular functional enrichments suggests unique responses to stressor duration and highlights distinct responses to bleaching in an anemone model.

Gene correlation networks reveal the transcriptomic response to elevated nitrogen in a photosynthetic sponge.

Nutrient levels in coastal environments have been increasing globally due to elevated inputs of sewage and terrigenous sediments carrying fertilizers. Yet, despite their immense filtering capacities, marine sponges appear to be less affected by elevated nutrients than sympatric benthic organisms, such as corals. While the molecular-level stress response of sponges to elevated seawater temperatures and other toxicants has been defined, this study represents the first global gene expression analysis of how sponges respond to elevated nitrogen. Gene correlation network analysis revealed that sponge gene modules, coded by colours, became either highly upregulated (Blue) or downregulated (Turquoise, Black, Brown) as nitrogen treatment levels increased. Gene Ontology enrichment analysis of the different modules revealed genes involved in cell signalling, immune response and flagella motility were affected by increasing nitrogen levels. Notably, a decrease in the regulation of NF-kappaB signalling and an increase in protein degradation was identified, which is comparable to metabolic pathways associated with the sponge thermal stress response. These results highlight that Cymbastela stipitata can rapidly respond to changes in the external environment and identifies pathways that probably contribute to the ability of C. stipitata to tolerate short-term nutrient pulses.

Symbiont community diversity is more variable in corals that respond poorly to stress.

Coral reefs are declining globally as climate change and local water quality press environmental conditions beyond the physiological tolerances of holobionts-the collective of the host and its microbial symbionts. To assess the relationship between symbiont composition and holobiont stress tolerance, community diversity metrics were quantified for dinoflagellate endosymbionts (Family: Symbiodiniaceae) from eight Acropora millepora genets that thrived under or responded poorly to various stressors. These eight selected genets represent the upper and lower tails of the response distribution of 40 coral genets that were exposed to four stress treatments (and control conditions) in a 10-day experiment. Specifically, four 'best performer' coral genets were analyzed at the end of the experiment because they survived high temperature, high pCO2 , bacterial exposure, or combined stressors, whereas four 'worst performer' genets were characterized because they experienced substantial mortality under these stressors. At the end of the experiment, seven of eight coral genets mainly hosted Cladocopium symbionts, whereas the eighth genet was dominated by both Cladocopium and Durusdinium symbionts. Symbiodiniaceae alpha and beta diversity were higher in worst performing genets than in best performing genets. Symbiont communities in worst performers also differed more after stress exposure relative to their controls (based on normalized proportional differences in beta diversity), than did best performers. A generalized joint attribute model estimated the influence of host genet and treatment on Symbiodiniaceae community composition and identified strong associations among particular symbionts and host genet performance, as well as weaker associations with treatment. Although dominant symbiont physiology and function contribute to host performance, these findings emphasize the importance of symbiont community diversity and stochasticity as components of host performance. Our findings also suggest that symbiont community diversity metrics may function as indicators of resilience and have potential applications in diverse disciplines from climate change adaptation to agriculture and medicine.

Molecular convergence and positive selection associated with the evolution of symbiont transmission mode in stony corals.

Heritable symbioses have been critical for the evolution of life. The genetic consequences of evolving a heritable symbiosis from the perspective of the symbiont are well established, but concomitant changes in the host remain unresolved. In stony corals, heritable, vertical transmission has evolved repeatedly, providing a unique opportunity to investigate the genomic basis of this complex trait. We conducted a comparative analysis of 25 coral transcriptomes to identify orthologous genes exhibiting signatures of positive selection and convergent amino acid substitutions in vertically transmitting lineages. The frequency of convergence events tends to be higher among vertically transmitting lineages, consistent with the proposed role of selection in driving the evolution of convergent transmission mode phenotypes. Of 10 774 orthologous genes, 403 exhibited at least one molecular convergence event and evidence of positive selection in at least one vertically transmitting lineage. Functional enrichments among these top candidate genes include processes previously implicated in symbiosis including endocytosis, immune response, cytoskeletal protein binding and cytoplasmic membrane-bounded vesicles. Finally, several novel candidates were identified among 100 genes showing evidence of positive selection at the particular convergence event, highlighting the value of our approach for generating new insight into host mechanisms associated with the evolution of heritable symbioses.

Positive genetic associations among fitness traits support evolvability of a reef-building coral under multiple stressors.

Climate change threatens organisms in a variety of interactive ways that requires simultaneous adaptation of multiple traits. Predicting evolutionary responses requires an understanding of the potential for interactions among stressors and the genetic variance and covariance among fitness-related traits that may reinforce or constrain an adaptive response. Here we investigate the capacity of Acropora millepora, a reef-building coral, to adapt to multiple environmental stressors: rising sea surface temperature, ocean acidification, and increased prevalence of infectious diseases. We measured growth rates (weight gain), coral color (a proxy for Symbiodiniaceae density), and survival, in addition to nine physiological indicators of coral and algal health in 40 coral genets exposed to each of these three stressors singly and combined. Individual stressors resulted in predicted responses (e.g., corals developed lesions after bacterial challenge and bleached under thermal stress). However, corals did not suffer substantially more when all three stressors were combined. Nor were trade-offs observed between tolerances to different stressors; instead, individuals performing well under one stressor also tended to perform well under every other stressor. An analysis of genetic correlations between traits revealed positive covariances, suggesting that selection to multiple stressors will reinforce rather than constrain the simultaneous evolution of traits related to holobiont health (e.g., weight gain and algal density). These findings support the potential for rapid coral adaptation under climate change and emphasize the importance of accounting for corals' adaptive capacity when predicting the future of coral reefs.

Considerations for maximizing the adaptive potential of restored coral populations in the western Atlantic.

Active coral restoration typically involves two interventions: crossing gametes to facilitate sexual larval propagation; and fragmenting, growing, and outplanting adult colonies to enhance asexual propagation. From an evolutionary perspective, the goal of these efforts is to establish self-sustaining, sexually reproducing coral populations that have sufficient genetic and phenotypic variation to adapt to changing environments. Here, we provide concrete guidelines to help restoration practitioners meet this goal for most Caribbean species of interest. To enable the persistence of coral populations exposed to severe selection pressure from many stressors, a mixed provenance strategy is suggested: genetically unique colonies (genets) should be sourced both locally as well as from more distant, environmentally distinct sites. Sourcing three to four genets per reef along environmental gradients should be sufficient to capture a majority of intraspecies genetic diversity. It is best for practitioners to propagate genets with one or more phenotypic traits that are predicted to be valuable in the future, such as low partial mortality, high wound healing rate, high skeletal growth rate, bleaching resilience, infectious disease resilience, and high sexual reproductive output. Some effort should also be reserved for underperforming genets because colonies that grow poorly in nurseries sometimes thrive once returned to the reef and may harbor genetic variants with as yet unrecognized value. Outplants should be clustered in groups of four to six genets to enable successful fertilization upon maturation. Current evidence indicates that translocating genets among distant reefs is unlikely to be problematic from a population genetic perspective but will likely provide substantial adaptive benefits. Similarly, inbreeding depression is not a concern given that current practices only raise first-generation offspring. Thus, proceeding with the proposed management strategies even in the absence of a detailed population genetic analysis of the focal species at sites targeted for restoration is the best course of action. These basic guidelines should help maximize the adaptive potential of reef-building corals facing a rapidly changing environment.

Family matters: Variation in the physiology of brooded Porites astreoides larvae is driven by parent colony effects.

The planktonic larval phase of scleractinian coral life-history represents a crucial stage when dispersal takes place and genetic diversity among populations is maintained. Understanding the dynamics influencing larval survival is especially relevant in the context of climate change, as larvae may be more vulnerable to environmental disturbances than adults. Several physiological parameters of coral larvae have been shown to vary by release time and past environmental history. However, the contribution of parental or genetic effects is largely unknown. To investigate these potential familial effects, we collected adult Porites astreoides colonies in April 2018 from two reef zones in the lower Florida Keys and quantified physiological traits and thermal tolerance of the newly released larvae. Family accounted for more variation than day of release and reef origin, with >60% of the variation in chlorophyll a and protein content explained by family. The survivorship of larvae under 36 °C acute temperature stress was also tightly linked to what parent colony they were released from. During a 32 °C moderate temperature stress experiment, inshore larvae tended to bleach less than offshore larvae, mirroring the enhanced bleaching resistance previously observed in inshore adult coral populations. The significant familial effects identified in the present study suggest that researchers should be cautious when interpreting results of studies which pool larvae among families, and that future studies should take care to account for this variation.

Functional genomic analysis of corals from natural CO2 -seeps reveals core molecular responses involved in acclimatization to ocean acidification.

Little is known about the potential for acclimatization or adaptation of corals to ocean acidification and even less about the molecular mechanisms underpinning these processes. Here, we examine global gene expression patterns in corals and their intracellular algal symbionts from two replicate population pairs in Papua New Guinea that have undergone long-term acclimatization to natural variation in pCO2 . In the coral host, only 61 genes were differentially expressed in response to pCO2 environment, but the pattern of change was highly consistent between replicate populations, likely reflecting the core expression homeostasis response to ocean acidification. Functional annotations highlight lipid metabolism and a change in the stress response capacity of corals as key parts of this process. Specifically, constitutive downregulation of molecular chaperones was observed, which may impact response to combined climate change-related stressors. Elevated CO2 has been hypothesized to benefit photosynthetic organisms but expression changes of in hospite Symbiodinium in response to acidification were greater and less consistent among reef populations. This population-specific response suggests hosts may need to adapt not only to an acidified environment, but also to changes in their Symbiodinium populations that may not be consistent among environments, adding another challenging dimension to the physiological process of coping with climate change.

Fine-scale environmental specialization of reef-building corals might be limiting reef recovery in the Florida Keys.

Despite decades of monitoring global reef decline, we are still largely unable to explain patterns of reef deterioration at local scales, which precludes the development of effective management strategies. Offshore reefs of the Florida Keys, USA, experience milder temperatures and lower nutrient loads in comparison to inshore reefs yet remain considerably more degraded than nearshore patch reefs. A year-long reciprocal transplant experiment of the mustard hill coral (Porites astreoides) involving four source and eight transplant locations reveals that corals adapt and/or acclimatize to their local habitat on a < 10-km scale. Surprisingly, transplantation to putatively similar environmental types (e.g., offshore corals moved to a novel offshore site, or along-shore transplantation) resulted in greater reductions in fitness proxies, such as coral growth, than cross-channel transplantation between inshore and offshore reefs. The only abiotic factor showing significantly greater differences between along-shore sites was daily temperature range extremes (rather than the absolute high or low temperatures reached), providing a possible explanation for this pattern. Offshore-origin corals exhibited significant growth reductions at sites with greater daily temperature ranges, which explained up to 39% of the variation in their mass gain. In contrast, daily temperature range explained at most 9% of growth variation in inshore-origin corals, suggesting that inshore corals are more tolerant of high-frequency temperature fluctuations. Finally, corals incur trade-offs when specializing to their native reef. Across reef locations the coefficient of selection against coral transplants was 0.07 ± 0.02 (mean ± SE). This selection against immigrants could hinder the ability of corals to recolonize devastated reefs, whether through assisted migration efforts or natural recruitment events, providing a unifying explanation for observed patterns of coral decline in this reef system.

Microhabitat acclimatization alters sea anemone-algal symbiosis and thermal tolerance across the intertidal zone.

Contemporary symbioses in extreme environments can give an insight into mechanisms that stabilize species interactions during environmental change. The intertidal sea anemone, Anthopleura elegantissima, engages in a nutritional symbiosis with microalgae similar to tropical coral, but withstands more intense environmental fluctuations during tidal inundations. In this study, we compare baseline symbiotic traits and their sensitivity to thermal stress within and among anemone aggregations across the intertidal using a laboratory-based tank experiment to better understand how fixed genotypic and plastic environmental effects contribute to the successful maintenance of this symbiosis in extreme habitats. High intertidal anemones had lower baseline symbiont-to-host cell ratios under control conditions, but their symbionts had higher baseline photosynthetic efficiency compared to low intertidal anemone symbionts. Symbiont communities were identical across all samples, suggesting that shifts in symbiont density and photosynthetic performance could be an acclimatory mechanism to maintain symbiosis in different environments. Despite lower baseline symbiont-to-host cell ratios, high intertidal anemones maintained greater symbiont-to-host cell ratios under heat stress compared with low intertidal anemones, suggesting greater thermal tolerance of high intertidal holobionts. However, the thermal tolerance of clonal anemones acclimatized to different zones was not explained by tidal height alone, indicating additional environmental variables contribute to physiological differences. Host genotype significantly influenced anemone weight, but only explained a minor proportion of variation among symbiotic traits and their response to thermal stress, further implicating environmental history as the primary driver of holobiont tolerance. These results indicate that this symbiosis is highly plastic and may be able to acclimatize to climate change over ecological timescales, defying the convention that symbiotic organisms are more susceptible to environmental stress.

Abundance of Oligoflexales bacteria is associated with algal symbiont density, independent of thermal stress in Aiptasia anemones.

Many multicellular organisms, such as humans, plants, and invertebrates, depend on symbioses with microbes for metabolic cooperation and exchange. Reef-building corals, an ecologically important order of invertebrates, are particularly vulnerable to environmental stress in part because of their nutritive symbiosis with dinoflagellate algae, and yet also benefit from these and other microbial associations. While coral microbiomes remain difficult to study because of their complexity, the anemone Aiptasia is emerging as a simplified model. Research has demonstrated co-occurrences between microbiome composition and the abundance and type of algal symbionts in cnidarians. However, whether these patterns are the result of general stress-induced shifts or depletions of algal-associated bacteria remains unclear. Our study aimed to distinguish the effect of changes in symbiont density and thermal stress on the microbiome of symbiotic Aiptasia strain CC7 by comparing them with aposymbiotic anemones, depleted of their native symbiont, Symbiodinium linucheae. Our analysis indicated that overall thermal stress had the greatest impact on disrupting the microbiome. We found that three bacterial classes made up most of the relative abundance (60%-85%) in all samples, but the rare microbiome fluctuated between symbiotic states and following thermal stress. We also observed that S. linucheae density correlated with abundance of Oligoflexales, suggesting these bacteria may be primary symbionts of the dinoflagellate algae. The findings of this study help expand knowledge on prospective multipartite symbioses in the cnidarian holobiont and how they respond to environmental disturbance.

Complete Genome Sequence of Roseibium sp. Strain Sym1, a Bacterial Associate of Symbiodinium linucheae, the Microalgal Symbiont of the Anemone Aiptasia.

We sequenced the genome of Roseibium sp. strain Sym1, a strain isolated from a monoculture of a Symbiodiniaceae marine dinoflagellate, Symbiodinium linucheae, a microalgal symbiont of cnidarians. The completed genome consists of one circular chromosome of 6,694,563 bp and four plasmids of 192,102 bp, 160,136 bp, 120,881 bp, and 89,413 bp.

Long-term maintenance of a heterologous symbiont association in Acropora palmata on natural reefs.

The sensitivity of reef-building coral to elevated temperature is a function of their symbiosis with dinoflagellate algae in the family Symbiodiniaceae. Changes in the composition of the endosymbiont community in response to thermal stress can increase coral thermal tolerance. Consequently, this mechanism is being investigated as a human-assisted intervention for rapid acclimation of coral in the face of climate change. Successful establishment of novel symbioses that increase coral thermal tolerance have been demonstrated in laboratory conditions; however, it is unclear how long these heterologous relationships persist in nature. Here, we test the persistence of a novel symbiosis between Acropora palmata and Durusdinium spp. from Mote Marine Laboratory's ex situ nursery by outplanting clonal replicates (ramets) of five A. palmata host genotypes to natural reefs in the lower Florida Keys. Amplicon sequencing analysis of ITS2-type profiles revealed that the majority of surviving ramets remained dominated by Durusdinium spp. two years after transplantation. However, 15% of ramets, including representatives of all genotypes, exhibited some degree of symbiont shuffling or switching at six of eight sites, including complete takeover by site-specific strains of the native symbiont, Symbiodinium fitti. The predominant long-term stability of the novel symbiosis supports the potential effectiveness of symbiont modification as a management tool. Although, the finding that 6-7 year-old coral can alter symbiont community composition in the absence of bleaching indicates that Symbiodiniaceae communities are indeed capable of great flexibility under ambient conditions.

Building consensus around the assessment and interpretation of Symbiodiniaceae diversity.

Within microeukaryotes, genetic variation and functional variation sometimes accumulate more quickly than morphological differences. To understand the evolutionary history and ecology of such lineages, it is key to examine diversity at multiple levels of organization. In the dinoflagellate family Symbiodiniaceae, which can form endosymbioses with cnidarians (e.g., corals, octocorals, sea anemones, jellyfish), other marine invertebrates (e.g., sponges, molluscs, flatworms), and protists (e.g., foraminifera), molecular data have been used extensively over the past three decades to describe phenotypes and to make evolutionary and ecological inferences. Despite advances in Symbiodiniaceae genomics, a lack of consensus among researchers with respect to interpreting genetic data has slowed progress in the field and acted as a barrier to reconciling observations. Here, we identify key challenges regarding the assessment and interpretation of Symbiodiniaceae genetic diversity across three levels: species, populations, and communities. We summarize areas of agreement and highlight techniques and approaches that are broadly accepted. In areas where debate remains, we identify unresolved issues and discuss technologies and approaches that can help to fill knowledge gaps related to genetic and phenotypic diversity. We also discuss ways to stimulate progress, in particular by fostering a more inclusive and collaborative research community. We hope that this perspective will inspire and accelerate coral reef science by serving as a resource to those designing experiments, publishing research, and applying for funding related to Symbiodiniaceae and their symbiotic partnerships.

Reduced representation sequencing accurately quantifies relative abundance and reveals population-level variation in Pseudo-nitzschia spp.

Certain species within the genus Pseudo-nitzschia are able to produce the neurotoxin domoic acid (DA), which can cause illness in humans, mass-mortality of marine animals, and closure of commercial and recreational shellfisheries during toxic events. Understanding and forecasting blooms of these harmful species is a primary management goal. However, accurately predicting the onset and severity of bloom events remains difficult, in part because the underlying drivers of bloom formation have not been fully resolved. Furthermore, Pseudo-nitzschia species often co-occur, and recent work suggests that the genetic composition of a Pseudo-nitzschia bloom may be a better predictor of toxicity than prevailing environmental conditions. We developed a novel next-generation sequencing assay using restriction site-associated DNA (2b-RAD) genotyping and applied it to mock Pseudo-nitzschia communities generated by mixing cultures of different species in known abundances. On average, 94% of the variance in observed species abundance was explained by the expected abundance. In addition, the false positive rate was low (0.45% on average) and unrelated to read depth, and false negatives were never observed. Application of this method to environmental DNA samples collected during natural Pseudo-nitzschia spp. bloom events in Southern California revealed that increases in DA were associated with increases in the relative abundance of P. australis. Although the absolute correlation across time-points was weak, an independent species fingerprinting assay (Automated Ribosomal Intergenic Spacer Analysis) supported this and identified other potentially toxic species. Finally, we assessed population-level genomic variation by mining SNPs from the environmental 2bRAD dataset. Consistent shifts in allele frequencies in P. pungens and P. subpacifica were detected between high and low DA years, suggesting that different intraspecific variants may be associated with prevailing environmental conditions or the presence of DA. Taken together, this method presents a potentially cost-effective and high-throughput approach for studies aiming to evaluate both population and species dynamics in mixed samples.