Optimal nutrient exchange and immune responses operate in partner specificity in the cnidarian-dinoflagellate symbiosis.

The relationship between corals and dinoflagellates of the genus Symbiodinium is fundamental to the functioning of coral ecosystems. It has been suggested that reef corals may adapt to climate change by changing their dominant symbiont type to a more thermally tolerant one, although the capacity for such a shift is potentially hindered by the compatibility of different host-symbiont pairings. Here we combined transcriptomic and metabolomic analyses to characterize the molecular, cellular, and physiological processes that underlie this compatibility, with a particular focus on Symbiodinium trenchii, an opportunistic, thermally tolerant symbiont that flourishes in coral tissues after bleaching events. Symbiont-free individuals of the sea anemone Exaiptasia pallida (commonly referred to as Aiptasia), an established model system for the study of the cnidarian-dinoflagellate symbiosis, were colonized with the "normal" (homologous) symbiont Symbiodinium minutum and the heterologous S. trenchii Analysis of the host gene and metabolite expression profiles revealed that heterologous symbionts induced an expression pattern intermediate between the typical symbiotic state and the aposymbiotic state. Furthermore, integrated pathway analysis revealed that increased catabolism of fixed carbon stores, metabolic signaling, and immune processes occurred in response to the heterologous symbiont type. Our data suggest that both nutritional provisioning and the immune response induced by the foreign "invader" are important factors in determining the capacity of corals to adapt to climate change through the establishment of novel symbioses.

Eye development and photoreceptor differentiation in the cephalopod Doryteuthis pealeii.

Photoreception is a ubiquitous sensory ability found across the Metazoa, and photoreceptive organs are intricate and diverse in their structure. Although the morphology of the compound eye in Drosophila and the single-chambered eye in vertebrates have elaborated independently, the amount of conservation within the 'eye' gene regulatory network remains controversial, with few taxa studied. To better understand the evolution of photoreceptive organs, we established the cephalopod Doryteuthis pealeii as a lophotrochozoan model for eye development. Utilizing histological, transcriptomic and molecular assays, we characterize eye formation in Doryteuthis pealeii Through lineage tracing and gene expression analyses, we demonstrate that cells expressing Pax and Six genes incorporate into the lens, cornea and iris, and the eye placode is the sole source of retinal tissue. Functional assays demonstrate that Notch signaling is required for photoreceptor cell differentiation and retinal organization. This comparative approach places the canon of eye research in traditional models into perspective, highlighting complexity as a result of both conserved and convergent mechanisms.

Heritable variation in bleaching responses and its functional genomic basis in reef-building corals (Orbicella faveolata).

Reef-building corals are highly sensitive to rising ocean temperatures, and substantial adaptation will be required for corals and the ecosystems they support to persist in changing ocean conditions. Genetic variation that might support adaptive responses has been measured in larval stages of some corals, but these estimates remain unavailable for adult corals and the functional basis of this variation remains unclear. In this study, we focused on the potential for adaptation in Orbicella faveolata, a dominant reef-builder in the Caribbean. We conducted thermal stress experiments using corals collected from natural populations in Bocas del Toro, Panama, and used multilocus SNP genotypes to estimate genetic relatedness among samples. This allowed us to estimate narrow-sense heritability of variation in bleaching responses, revealing that variation in these responses was highly heritable (h2  = 0.58). This suggests substantial potential for adaptive responses to warming by natural populations of O. faveolata in this region. We further investigated the functional basis for this variation using genomic and transcriptomic approaches. We used a publicly available genetic linkage map and genome assembly to map markers associated with bleaching responses, identifying twelve markers associated with variation in bleaching responses. We also profiled gene expression in corals with contrasting bleaching phenotypes, uncovering substantial differences in transcriptional stress responses between heat-tolerant and heat-susceptible corals. Together, our findings contribute to the growing body of evidence that natural populations of corals possess genetic variation in thermal stress responses that may potentially support adaptive responses to rising ocean temperatures.

Disparate responses to salinity across species and organizational levels in anchialine shrimps.

Environmentally induced plasticity in gene expression is one of the underlying mechanisms of adaptation to habitats with variable environments. For example, euryhaline crustaceans show predictable changes in the expression of ion-transporter genes during salinity transfers, although studies have typically been limited to specific genes, taxa and ecosystems of interest. Here, we investigated responses to salinity change at multiple organizational levels in five species of shrimp representing at least three independent invasions of the anchialine ecosystem, defined as habitats with marine and freshwater influences with spatial and temporal fluctuations in salinity. Although all five species were generally strong osmoregulators, salinity-induced changes in gill physiology and gene expression were highly species specific. While some species exhibited patterns similar to those of previously studied euryhaline crustaceans, instances of distinct and atypical patterns were recovered from closely related species. Species-specific patterns were found when examining: (1) numbers and identities of differentially expressed genes, (2) salinity-induced expression of genes predicted a priori to play a role in osmoregulation, and (3) salinity-induced expression of orthologs shared among all species. Notably, ion transport genes were unchanged in the atyid Halocaridina rubra while genes normally associated with vision and light perception were among those most highly upregulated. Potential reasons for species-specific patterns are discussed, including variation among anchialine habitats in salinity regimes and divergent evolution in anchialine taxa. Underexplored mechanisms of osmoregulation in crustaceans revealed here by the application of transcriptomic approaches to ecologically and taxonomically understudied systems are also explored.

Genomic and transcriptomic signals of thermal tolerance in heat-tolerant corals (Platygyra daedalea) of the Arabian/Persian Gulf.

Scleractinian corals occur in tropical regions near their upper thermal limits and are severely threatened by rising ocean temperatures. However, several recent studies have shown coral populations can harbour genetic variation in thermal tolerance. Here, we have extended these approaches to study heat tolerance of corals in the Persian/Arabian Gulf, where heat-tolerant local populations experience extreme summer temperatures (up to 36°C). To evaluate whether selection has depleted genetic variation in thermal tolerance, estimate potential future adaptive responses and understand the functional basis for these corals' unusual heat tolerance, we conducted controlled crosses in the Gulf coral Platygyra daedalea. Heat tolerance is highly heritable in this population (h 2  = 0.487-0.748), suggesting substantial potential for adaptive responses to selection for elevated temperatures. To identify genetic markers associated with this variation, we conducted genomewide SNP genotyping in parental corals and tested for relationships between paternal genotype and offspring thermal tolerance. Resulting multilocus SNP genotypes explained a large fraction of variation in thermal tolerance in these crosses (69%). To investigate the functional basis of these differences in thermal tolerance, we profiled transcriptional responses in tolerant and susceptible families, revealing substantial sire effects on transcriptional responses to thermal stress. We also studied sequence variation in these expressed sequences, identifying alleles and functional groups of differentially expressed genes associated with thermal tolerance. Our findings demonstrate that corals in this population harbour extensive genetic variation in thermal tolerance, and heat-tolerant phenotypes differ in both gene sequences and transcriptional stress responses from their susceptible counterparts.

Impacts of temperature and lunar day on gene expression profiles during a monthly reproductive cycle in the brooding coral Pocillopora damicornis.

Reproductive timing in brooding corals has been correlated to temperature and lunar irradiance, but the mechanisms by which corals transduce these environmental variables into molecular signals are unknown. To gain insight into these processes, global gene expression profiles in the coral Pocillopora damicornis were examined (via RNA-Seq) across lunar phases and between temperature treatments, during a monthly planulation cycle. The interaction of temperature and lunar day together had the largest influence on gene expression. Mean timing of planulation, which occurred at lunar days 7.4 and 12.5 for 28- and 23°C-treated corals, respectively, was associated with an upregulation of transcripts in individual temperature treatments. Expression profiles of planulation-associated genes were compared between temperature treatments, revealing that elevated temperatures disrupted expression profiles associated with planulation. Gene functions inferred from homologous matches to online databases suggest complex neuropeptide signalling, with calcium as a central mediator, acting through tyrosine kinase and G protein-coupled receptor pathways. This work contributes to our understanding of coral reproductive physiology and the impacts of environmental variables on coral reproductive pathways.

Host adaptation and unexpected symbiont partners enable reef-building corals to tolerate extreme temperatures.

Understanding the potential for coral adaptation to warming seas is complicated by interactions between symbiotic partners that define stress responses and the difficulties of tracking selection in natural populations. To overcome these challenges, we characterized the contribution of both animal host and symbiotic algae to thermal tolerance in corals that have already experienced considerable warming on par with end-of-century projections for most coral reefs. Thermal responses in Platygyra daedalea corals from the hot Persian Gulf where summer temperatures reach 36°C were compared with conspecifics from the milder Sea of Oman. Persian Gulf corals had higher rates of survival at elevated temperatures (33 and 36°C) in both the nonsymbiotic larval stage (32-49% higher) and the symbiotic adult life stage (51% higher). Additionally, Persian Gulf hosts had fixed greater potential to mitigate oxidative stress (31-49% higher) and their Symbiodinium partners had better retention of photosynthetic performance under elevated temperature (up to 161% higher). Superior thermal tolerance of Persian Gulf vs. Sea of Oman corals was maintained after 6-month acclimatization to a common ambient environment and was underpinned by genetic divergence in both the coral host and symbiotic algae. In P. daedalea host samples, genomewide SNP variation clustered into two discrete groups corresponding with Persian Gulf and Sea of Oman sites. Symbiodinium within host tissues predominantly belonged to ITS2 rDNA type C3 in the Persian Gulf and type D1a in the Sea of Oman contradicting patterns of Symbiodinium thermal tolerance from other regions. Our findings provide evidence that genetic adaptation of both host and Symbiodinium has enabled corals to cope with extreme temperatures in the Persian Gulf. Thus, the persistence of coral populations under continued warming will likely be determined by evolutionary rates in both, rather than single, symbiotic partners.

Gene expression associated with white syndromes in a reef building coral, Acropora hyacinthus.

BACKGROUND: Corals are capable of launching diverse immune defenses at the site of direct contact with pathogens, but the molecular mechanisms of this activity and the colony-wide effects of such stressors remain poorly understood. Here we compared gene expression profiles in eight healthy Acropora hyacinthus colonies against eight colonies exhibiting tissue loss commonly associated with white syndromes, all collected from a natural reef environment near Palau. Two types of tissues were sampled from diseased corals: visibly affected and apparently healthy. RESULTS: Tag-based RNA-Seq followed by weighted gene co-expression network analysis identified groups of co-regulated differentially expressed genes between all health states (disease lesion, apparently healthy tissues of diseased colonies, and fully healthy). Differences between healthy and diseased tissues indicate activation of several innate immunity and tissue repair pathways accompanied by reduced calcification and the switch towards metabolic reliance on stored lipids. Unaffected parts of diseased colonies, although displaying a trend towards these changes, were not significantly different from fully healthy samples. Still, network analysis identified a group of genes, suggestive of altered immunity state, that were specifically up-regulated in unaffected parts of diseased colonies. CONCLUSIONS: Similarity of fully healthy samples to apparently healthy parts of diseased colonies indicates that systemic effects of white syndromes on A. hyacinthus are weak, which implies that the coral colony is largely able to sustain its physiological performance despite disease. The genes specifically up-regulated in unaffected parts of diseased colonies, instead of being the consequence of disease, might be related to the originally higher susceptibility of these colonies to naturally occurring white syndromes.

2b-RAD: a simple and flexible method for genome-wide genotyping.

We describe 2b-RAD, a streamlined restriction site-associated DNA (RAD) genotyping method based on sequencing the uniform fragments produced by type IIB restriction endonucleases. Well-studied accessions of Arabidopsis thaliana were genotyped to validate the method's accuracy and to demonstrate fine-tuning of marker density as needed. The simplicity of the 2b-RAD protocol makes it particularly suitable for high-throughput genotyping as required for linkage mapping and profiling genetic variation in natural populations.

The genetics of divergence and reproductive isolation between ecotypes of Panicum hallii.

The process of plant speciation often involves the evolution of divergent ecotypes in response to differences in soil water availability between habitats. While the same set of traits is frequently associated with xeric/mesic ecotype divergence, it is unknown whether those traits evolve independently or if they evolve in tandem as a result of genetic colocalization either by pleiotropy or genetic linkage. The self-fertilizing C4 grass species Panicum hallii includes two major ecotypes found in xeric (var. hallii) or mesic (var. filipes) habitats. We constructed the first linkage map for P. hallii by genotyping a reduced representation genomic library of an F2 population derived from an intercross of var. hallii and filipes. We then evaluated the genetic architecture of divergence between these ecotypes through quantitative trait locus (QTL) mapping. Overall, we mapped QTLs for nine morphological traits that are involved in the divergence between the ecotypes. QTLs for five key ecotype-differentiating traits all colocalized to the same region of linkage group five. Leaf physiological traits were less divergent between ecotypes, but we still mapped five physiological QTLs. We also discovered a two-locus Dobzhansky-Muller hybrid incompatibility. Our study suggests that ecotype-differentiating traits may evolve in tandem as a result of genetic colocalization.

Integrating transcriptional, metabolomic, and physiological responses to drought stress and recovery in switchgrass (Panicum virgatum L.).

BACKGROUND: In light of the changes in precipitation and soil water availability expected with climate change, understanding the mechanisms underlying plant responses to water deficit is essential. Toward that end we have conducted an integrative analysis of responses to drought stress in the perennial C4 grass and biofuel crop, Panicum virgatum (switchgrass). Responses to soil drying and re-watering were measured at transcriptional, physiological, and metabolomic levels. To assess the interaction of soil moisture with diel light: dark cycles, we profiled gene expression in drought and control treatments under pre-dawn and mid-day conditions. RESULTS: Soil drying resulted in reduced leaf water potential, gas exchange, and chlorophyll fluorescence along with differential expression of a large fraction of the transcriptome (37%). Many transcripts responded differently depending on time of day (e.g. up-regulation pre-dawn and down-regulation mid-day). Genes associated with C4 photosynthesis were down-regulated during drought, while C4 metabolic intermediates accumulated. Rapid changes in gene expression were observed during recovery from drought, along with increased water use efficiency and chlorophyll fluorescence. CONCLUSIONS: Our findings demonstrate that drought responsive gene expression depends strongly on time of day and that gene expression is extensively modified during the first few hours of drought recovery. Analysis of covariation in gene expression, metabolite abundance, and physiology among plants revealed non-linear relationships that suggest critical thresholds in drought stress responses. Future studies may benefit from evaluating these thresholds among diverse accessions of switchgrass and other C4 grasses.

Transcriptome analysis and gene expression atlas for Panicum hallii var. filipes, a diploid model for biofuel research.

Panicum hallii is an emerging model for genetic studies of agronomic traits in Panicum, presenting a tractable diploid alternative study system to the tetra- or octaploid biofuel crop switchgrass (Panicum virgatum). To characterize the gene complement in P. hallii var. filipes and enable gene expression analysis in this system we sequenced, assembled, and annotated the transcriptome. Over 300 Mb of normalized cDNA prepared from multiple tissues and treatments was sequenced using 454-Titanium, producing an annotated assembly including 15 422 unique gene names. Comparison with other grass genomes identified putative P. hallii homologs for >14 000 previously characterized genes. We also developed an atlas of gene expression across tissues and stages using RNA-Seq (the quantitative analysis of short cDNA reads). SOLiD sequencing and quantitative analysis of more than 40 million cDNA tags identified substantial variation in expression profiles among tissues, consistent with known functional differences. Putative homologs were found for all enzymes in the phenylpropanoid pathway leading to lignin biosynthesis, including genes with known effects on biomass conversion efficiency. The resources developed here will enable studies of the genes underlying variation in cell wall composition, drought tolerance, and biomass production in Panicum.

Development of a next-generation NIL library in Arabidopsis thaliana for dissecting complex traits.

BACKGROUND: The identification of the loci and specific alleles underlying variation in quantitative traits is an important goal for evolutionary biologists and breeders. Despite major advancements in genomics technology, moving from QTL to causal alleles remains a major challenge in genetics research. Near-isogenic lines are the ideal raw material for QTL validation, refinement of QTL location and, ultimately, gene discovery. RESULTS: In this study, a population of 75 Arabidopsis thaliana near-isogenic lines was developed from an existing recombinant inbred line (RIL) population derived from a cross between physiologically divergent accessions Kas-1 and Tsu-1. First, a novel algorithm was developed to utilize genome-wide marker data in selecting RILs fully isogenic to Kas-1 for a single chromosome. Seven such RILs were used in 2 generations of crossing to Tsu-1 to create BC1 seed. BC1 plants were genotyped with SSR markers so that lines could be selected that carried Kas-1 introgressions, resulting in a population carrying chromosomal introgressions spanning the genome. BC1 lines were genotyped with 48 genome-wide SSRs to identify lines with a targeted Kas-1 introgression and the fewest genomic introgressions elsewhere. 75 such lines were selected and genotyped at an additional 41 SNP loci and another 930 tags using 2b-RAD genotyping by sequencing. The final population carried an average of 1.35 homozygous and 2.49 heterozygous introgressions per line with average introgression sizes of 5.32 and 5.16 Mb, respectively. In a simple case study, we demonstrate the advantage of maintaining heterozygotes in our library whereby fine-mapping efforts are conducted simply by self-pollination. Crossovers in the heterozygous interval during this single selfing generation break the introgression into smaller, homozygous fragments (sub-NILs). Additionally, we utilize a homozygous NIL for validation of a QTL underlying stomatal conductance, a low heritability trait. CONCLUSIONS: The present results introduce a new and valuable resource to the Brassicaceae research community that enables rapid fine-mapping of candidate loci in parallel with QTL validation. These attributes along with dense marker coverage and genome-wide chromosomal introgressions make this population an ideal starting point for discovery of genes underlying important complex traits of agricultural and ecological significance.

Microsatellite markers for the native Texas perennial grass, Panicum hallii (Poaceae).

PREMISE OF THE STUDY: We developed microsatellites for Panicum hallii for studies of gene flow, population structure, breeding experiments, and genetic mapping. METHODS AND RESULTS: Next-generation (454) genomic sequence data were used to design markers. Eighteen robust markers were discovered, 15 of which were polymorphic across six accessions of P. hallii var. hallii. Fourteen of the markers cross-amplified in a P. capillare accession. For the 15 polymorphic markers, the total number of alleles per locus ranged from two to 26 (mean: 11.0) across six populations (11-19 individuals per population). Observed heterozygosity (mean: 0.031) was 13.7 times lower than the expected heterozygosity (mean: 0.426). CONCLUSIONS: The deficit of heterozygous individuals is consistent with P. hallii having a high rate of self-fertilization. These markers will be useful for studies in P. hallii and related species.

Multi-domain GFP-like proteins from two species of marine hydrozoans.

Proteins homologous to Green Fluorescent Protein (GFP) are widely used as genetically encoded fluorescent labels. Many developments of this technology were spurred by discoveries of novel types of GFP-like proteins (FPs) in nature. Here we report two proteins displaying primary structures never before encountered in natural FPs: they consist of multiple GFP-like domains repeated within the same polypeptide chain. A two-domain green FP (abeGFP) and a four-domain orange-fluorescent FP (Ember) were isolated from the siphonophore Abylopsis eschscholtzii and an unidentified juvenile jellyfish (order Anthoathecata), respectively. Only the most evolutionary ancient domain of Ember is able to synthesize an orange-emitting chromophore (emission at 571 nm), while the other three are purely green (emission at 520 nm) and putatively serve to maintain the stability and solubility of the multidomain protein. When expressed individually, two of the green Ember domains form dimers and the third one exists as a monomer. The low propensity for oligomerization of these domains would simplify their adoption as in vivo labels. Our results reveal a previously unrecognized direction in which natural FPs have diversified, suggesting new avenues to look for FPs with novel and potentially useful features.

Bacterial diversity in Solenopsis invicta and Solenopsis geminata ant colonies characterized by 16S amplicon 454 pyrosequencing.

Social insects harbor diverse assemblages of bacterial microbes, which may play a crucial role in the success or failure of biological invasions. The invasive fire ant Solenopsis invicta (Formicidae, Hymenoptera) is a model system for understanding the dynamics of invasive social insects and their biological control. However, little is known about microbes as biotic factors influencing the success or failure of ant invasions. This pilot study is the first attempt to characterize and compare microbial communities associated with the introduced S. invicta and the native Solenopsis geminata in the USA. Using 16S amplicon 454 pyrosequencing, bacterial communities of workers, brood, and soil from nest walls were compared between neighboring S. invicta and S. geminata colonies at Brackenridge Field Laboratory, Austin, Texas, with the aim of identifying potential pathogenic, commensal, or mutualistic microbial associates. Two samples of S. geminata workers showed high counts of Spiroplasma bacteria, a known pathogen or mutualist of other insects. A subsequent analysis using PCR and sequencing confirmed the presence of Spiroplasma in additional colonies of both Solenopsis species. Wolbachia was found in one alate sample of S. geminata, while one brood sample of S. invicta had a high count of Lactococcus. As expected, ant samples from both species showed much lower microbial diversity than the surrounding soil. Both ant species had similar overall bacterial diversities, although little overlap in specific microbes. To properly characterize a single bacterial community associated with a Solenopsis ant sample, rarefaction analyses indicate that it is necessary to obtain 5,000-10,000 sequences. Overall, 16S amplicon 454 pyrosequencing appears to be a cost-effective approach to screen whole microbial diversity associated with invasive ant species.

Larval development in the Pacific oyster and the impacts of ocean acidification: Differential genetic effects in wild and domesticated stocks.

The adaptive capacity of marine calcifiers to ocean acidification (OA) is a topic of great interest to evolutionary biologists and ecologists. Previous studies have provided evidence to suggest that larval resilience to high pCO2 seawater for these species is a trait with a genetic basis and variability in natural populations. To date, however, it remains unclear how the selective effects of OA occur within the context of complex genetic interactions underpinning larval development in many of the most vulnerable taxa. Here we evaluated phenotypic and genetic changes during larval development of Pacific oysters (Crassostrea gigas) reared in ambient (~400 µatm) and high (~1600 µatm) pCO2 conditions, both in domesticated and naturalized "wild" oysters from the Pacific Northwest, USA. Using pooled DNA samples, we determined changes in allele frequencies across larval development, from early "D-stage" larvae to metamorphosed juveniles (spat), in both groups and environments. Domesticated larvae had ~26% fewer loci with changing allele frequencies across developmental stages and <50% as many loci affected by acidified culture conditions, compared to larvae from wild broodstock. Functional enrichment analyses of genetic markers with significant changes in allele frequency revealed that the structure and function of cellular membranes were disproportionately affected by high pCO2 conditions in both groups. These results indicate the potential for a rapid adaptive response of oyster populations to OA conditions; however, underlying genetic changes associated with larval development differ between these wild and domesticated oyster stocks and influence their adaptive responses to OA conditions.

Enhancing the heat tolerance of reef-building corals to future warming.

Reef-building corals thriving in extreme thermal environments may provide genetic variation that can assist the evolution of populations to rapid climate warming. However, the feasibility and scale of genetic improvements remain untested despite ongoing population declines from recurrent thermal stress events. Here, we show that corals from the hottest reefs in the world transfer sufficient heat tolerance to a naïve population sufficient to withstand end-of-century warming projections. Heat survival increased up to 84% when naïve mothers were selectively bred with fathers from the hottest reefs because of strong heritable genetic effects. We identified genomic loci associated with tolerance variation that were enriched for heat shock proteins, oxidative stress, and immune functions. Unexpectedly, several coral families exhibited survival rates and genomic associations deviating from origin predictions, including a few naïve purebreds with exceptionally high heat tolerance. Our findings highlight previously uncharacterized enhanced and intrinsic potential of coral populations to adapt to climate warming.

Transcriptome-Wide Comparisons and Virulence Gene Polymorphisms of Host-Associated Genotypes of the Cnidarian Parasite Ceratonova shasta in Salmonids.

Ceratonova shasta is an important myxozoan pathogen affecting the health of salmonid fishes in the Pacific Northwest of North America. Ceratonova shasta exists as a complex of host-specific genotypes, some with low to moderate virulence, and one that causes a profound, lethal infection in susceptible hosts. High throughput sequencing methods are powerful tools for discovering the genetic basis of these host/virulence differences, but deep sequencing of myxozoans has been challenging due to extremely fast molecular evolution of this group, yielding strongly divergent sequences that are difficult to identify, and unavoidable host contamination. We designed and optimized different bioinformatic pipelines to address these challenges. We obtained a unique set of comprehensive, host-free myxozoan RNA-seq data from C. shasta genotypes of varying virulence from different salmonid hosts. Analyses of transcriptome-wide genetic distances and maximum likelihood multigene phylogenies elucidated the evolutionary relationship between lineages and demonstrated the limited resolution of the established Internal Transcribed Spacer marker for C. shasta genotype identification, as this marker fails to differentiate between biologically distinct genotype II lineages from coho salmon and rainbow trout. We further analyzed the data sets based on polymorphisms in two gene groups related to virulence: cell migration and proteolytic enzymes including their inhibitors. The developed single-nucleotide polymorphism-calling pipeline identified polymorphisms between genotypes and demonstrated that variations in both motility and protease genes were associated with different levels of virulence of C. shasta in its salmonid hosts. The prospective use of proteolytic enzymes as promising candidates for targeted interventions against myxozoans in aquaculture is discussed. We developed host-free transcriptomes of a myxozoan model organism from strains that exhibited different degrees of virulence, as a unique source of data that will foster functional gene analyses and serve as a base for the development of potential therapeutics for efficient control of these parasites.

Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx.

BACKGROUND: New methods are needed for genomic-scale analysis of emerging model organisms that exemplify important biological questions but lack fully sequenced genomes. For example, there is an urgent need to understand the potential for corals to adapt to climate change, but few molecular resources are available for studying these processes in reef-building corals. To facilitate genomics studies in corals and other non-model systems, we describe methods for transcriptome sequencing using 454, as well as strategies for assembling a useful catalog of genes from the output. We have applied these methods to sequence the transcriptome of planulae larvae from the coral Acropora millepora. RESULTS: More than 600,000 reads produced in a single 454 sequencing run were assembled into approximately 40,000 contigs with five-fold average sequencing coverage. Based on sequence similarity with known proteins, these analyses identified approximately 11,000 different genes expressed in a range of conditions including thermal stress and settlement induction. Assembled sequences were annotated with gene names, conserved domains, and Gene Ontology terms. Targeted searches using these annotations identified the majority of genes associated with essential metabolic pathways and conserved signaling pathways, as well as novel candidate genes for stress-related processes. Comparisons with the genome of the anemone Nematostella vectensis revealed approximately 8,500 pairs of orthologs and approximately 100 candidate coral-specific genes. More than 30,000 SNPs were detected in the coral sequences, and a subset of these validated by re-sequencing. CONCLUSION: The methods described here for deep sequencing of the transcriptome should be widely applicable to generate catalogs of genes and genetic markers in emerging model organisms. Our data provide the most comprehensive sequence resource currently available for reef-building corals, and include an extensive collection of potential genetic markers for association and population connectivity studies. The characterization of the larval transcriptome for this widely-studied coral will enable research into the biological processes underlying stress responses in corals and evolutionary adaptation to global climate change.