Unequal and genotype-dependent expression of mitochondrial genes in larvae of the pacific oyster Crassostrea gigas

Comparative transcriptomics identifies genes underlying growth performance of the Pacific black-lipped pearl oyster Pinctada margaritifera

BACKGROUND

In bivalves, the rate at which organisms grow is a major functional trait underlying many aspects of their commercial production. Growth is a highly polygenic trait, which is typically regulated by many genes with small to moderate effects. Due to its complexity, growth variability in such shellfish remains poorly understood. In this study, we aimed to investigate differential gene expression among spat of the pearl oyster Pinctada margaritifera with distinct growth phenotypes.

RESULTS

We selected two groups of P. margaritifera spat belonging to the same F2 cohort based on their growth performance at 5.5 months old. Transcriptome profile analysis identified a total of 394 differentially expressed genes between these Fast-growing (F) and Slow-growing (S) phenotypes. According to functional enrichment analysis, S oysters overexpressed genes associated with stress-pathways and regulation of innate immune responses. In contrast, F oysters up-regulated genes associated with cytoskeleton activity, cell proliferation, and apoptosis. Analysis of genome polymorphism identified 16 single nucleotide polymorphisms (SNPs) significantly associated with the growth phenotypes. SNP effect categorization revealed one SNP identified for high effect and annotated for a stop codon gained mutation. Interestingly, this SNP is located within a gene annotated for scavenger receptor class F member 1 (SRF1), which is known to modulate apoptosis. Our analyses also revealed that all F oysters showed up-regulation for this gene and were homozygous for the stop-codon mutation. Conversely, S oysters had a heterozygous genotype and a reduced expression of this gene.

CONCLUSIONS

Altogether, our findings suggest that differences in growth among the same oyster cohort may be explained by contrasted metabolic allocation between regulatory pathways for growth and the immune system. This study provides a valuable contribution towards our understanding of the molecular components associated with growth performance in the pearl oyster P. margaritifera and bivalves in general.

Mitochondrial evolution in the Demospongiae (Porifera): Phylogeny, divergence time, and genome biology

The sponge class Demospongiae is the most speciose and morphologically diverse in the phylum Porifera, and the species within it are vital components of a range of ecosystems worldwide. Despite their ubiquity, a number of recalcitrant problems still remain to be solved regarding their phylogenetic inter-relationships, the timing of their appearance, and their mitochondrial biology, the latter of which is only beginning to be investigated. Here we generated 14 new demosponge mitochondrial genomes which, alongside previously published mitochondrial resources, were used to address these issues. In addition to phylogenomic analysis, we have used syntenic data and analysis of coding regions to forge a framework for understanding the inter-relationships between Demospongiae sub-classes and orders. We have also leveraged our new resources to study the mitochondrial biology of these clades in terms of codon usage, optimisation and gene expression, to understand how these vital cellular components may have contributed to the success of the Porifera. Our results strongly support a sister relationship between Keratosa and (Verongimorpha + Heteroscleromorpha), contradicting previous studies using nuclear markers. Our study includes one species of Clionaida, and show for the first time support for a grouping of Suberitida+(Clionaida+(Tethyida + Poecilosclerida). The findings of our phylogenetic analyses are supported by in-depth examination of structural and coding-level evidence from our mitochondrial data. A time-calibrated phylogeny estimated the origin of Demospongiae in the Cambrian (~529 Mya), and suggests that most demosponge order crown-groups emerged in the Mesozoic. This work therefore provides a robust basis for considering demosponge phylogenetic relationships, as well as essential mitochondrial data for understanding the biological basis for their success and diversity.

Biochemical bases of growth variation during development: a study of protein turnover in pedigreed families of bivalve larvae (Crassostrea gigas)

Animal size is a highly variable trait regulated by complex interactions between biological and environmental processes. Despite the importance of understanding the mechanistic bases of growth, predicting size variation in early stages of development remains challenging. Pedigreed lines of the Pacific oyster (Crassostrea gigas) were crossed to produce contrasting growth phenotypes to analyze the metabolic bases of growth variation in larval stages. Under controlled environmental conditions, substantial growth variation of up to 430% in shell length occurred among 12 larval families. Protein was the major biochemical constituent in larvae, with an average protein-to-lipid content ratio of 2.8. On average, 86% of protein synthesized was turned over (i.e. only 14% retained as protein accreted), with a regulatory shift in depositional efficiency resulting in increased protein accretion during later larval growth. Variation in protein depositional efficiency among families did not explain the range in larval growth rates. Instead, changes in protein synthesis rates predicted 72% of growth variation. High rates of protein synthesis to support faster growth, in turn, necessitated greater allocation of the total ATP pool to protein synthesis. An ATP allocation model is presented for larvae of C. gigas that includes the major components (82%) of energy demand: protein synthesis (45%), ion pump activity (20%), shell formation (14%) and protein degradation (3%). The metabolic trade-offs between faster growth and the need for higher ATP allocation to protein synthesis could be a major determinant of fitness for larvae of different genotypes responding to the stress of environmental change.

Oyster reproduction is compromised by acidification experienced seasonally in coastal regions

Atmospheric carbon dioxide concentrations have been rising during the past century, leading to ocean acidification (OA). Coastal and estuarine habitats experience annual pH variability that vastly exceeds the magnitude of long-term projections in open ocean regions. Eastern oyster (Crassostrea virginica) reproduction season coincides with periods of low pH occurrence in estuaries, thus we investigated effects of moderate (pH 7.5, pCO₂ 2260 µatm) and severe OA (pH 7.1, pCO₂ 5584 µatm; and 6.7, pCO₂ 18480 µatm) on oyster gametogenesis, fertilization, and early larval development successes. Exposure at severe OA during gametogenesis caused disruption in oyster reproduction. Oogenesis appeared to be more sensitive compared to spermatogenesis. However, Eastern oyster reproduction was resilient to moderate OA projected for the near-future. In the context of projected climate change exacerbating seasonal acidification, OA of coastal habitats could represent a significant bottleneck for oyster reproduction which may have profound negative implications for coastal ecosystems reliant on this keystone species.

Metabolic Cost of Protein Synthesis in Larvae of the Pacific Oyster (Crassostrea gigas) Is Fixed Across Genotype, Phenotype, and Environmental Temperature

The energy made available through catabolism of specific biochemical reserves is constant using standard thermodynamic conversion equivalents (e.g., 24.0 J mg protein(-1)). In contrast, measurements reported for the energy cost of synthesis of specific biochemical constituents are highly variable. In this study, we measured the metabolic cost of protein synthesis and determined whether this cost was influenced by genotype, phenotype, or environment. We focused on larval stages of the Pacific oyster Crassostrea gigas, a species that offers several experimental advantages: availability of genetically pedigreed lines, manipulation of ploidy, and tractability of larval forms for in vivo studies of physiological processes. The cost of protein synthesis was measured in larvae of C. gigas for 1) multiple genotypes, 2) phenotypes with different growth rates, and 3) different environmental temperatures. For all treatments, the cost of protein synthesis was within a narrow range--near the theoretical minimum--with a fixed cost (mean ± one standard error, n = 21) of 2.1 ± 0.2 J (mg protein synthesized)(-1) We conclude that there is no genetic variation in the metabolic cost of protein synthesis, thereby simplifying bioenergetic models. Protein synthesis is a major component of larval metabolism in C. gigas, accounting for more than half the metabolic rate in diploid (59%) and triploid larvae (54%). These results provide measurements of metabolic cost of protein synthesis in larvae of C. gigas, an indicator species for impacts of ocean change, and provide a quantitative basis for evaluating the cost of resilience.

Physiological response and resilience of early life-stage Eastern oysters (Crassostrea virginica) to past, present and future ocean acidification

The Eastern oyster, Crassostrea virginica (Gmelin, 1791), is the second most valuable bivalve fishery in the USA and is sensitive to high levels of partial pressure of CO2 (pCO2). Here we present experiments that comprehensively examined how the ocean's past, present and projected (21st and 22nd centuries) CO2 concentrations impact the growth and physiology of larval stages of C. virginica. Crassostrea virginica larvae grown in present-day pCO2 concentrations (380 μatm) displayed higher growth and survival than individuals grown at both lower (250 μatm) and higher pCO2 levels (750 and 1500 μatm). Crassostrea virginica larvae manifested calcification rates, sizes, shell thicknesses, metamorphosis, RNA:DNA ratios and lipid contents that paralleled trends in survival, with maximal values for larvae grown at 380 μatm pCO2 and reduced performance in higher and lower pCO2 levels. While some physiological differences among oysters could be attributed to CO2-induced changes in size or calcification rates, the RNA:DNA ratios at ambient pCO2 levels were elevated, independent of these factors. Likewise, the lipid contents of individuals exposed to high pCO2 levels were depressed even when differences in calcification rates were considered. These findings reveal the cascading, interdependent impact that high CO2 can have on oyster physiology. Crassostrea virginica larvae are significantly more resistant to elevated pCO2 than other North Atlantic bivalves, such as Mercenaria mercenaria and Argopecten irradians, a finding that may be related to the biogeography and/or evolutionary history of these species and may have important implications for future bivalve restoration and aquaculture efforts.

Evidence for a fourteenth mtDNA-encoded protein in the female-transmitted mtDNA of marine Mussels (Bivalvia: Mytilidae)

Background

A novel feature for animal mitochondrial genomes has been recently established: i.e., the presence of additional, lineage-specific, mtDNA-encoded proteins with functional significance. This feature has been observed in freshwater mussels with doubly uniparental inheritance of mtDNA (DUI). The latter unique system of mtDNA transmission, which also exists in some marine mussels and marine clams, is characterized by one mt genome inherited from the female parent (F mtDNA) and one mt genome inherited from the male parent (M mtDNA). In freshwater mussels, the novel mtDNA-encoded proteins have been shown to be mt genome-specific (i.e., one novel protein for F genomes and one novel protein for M genomes). It has been hypothesized that these novel, F- and M-specific, mtDNA-encoded proteins (and/or other F- and/or M-specific mtDNA sequences) could be responsible for the different modes of mtDNA transmission in bivalves but this remains to be demonstrated.

Methodology/Principal Findings

We investigated all complete (or nearly complete) female- and male-transmitted marine mussel mtDNAs previously sequenced for the presence of ORFs that could have functional importance in these bivalves. Our results confirm the presence of a novel F genome-specific mt ORF, of significant length (>100aa) and located in the control region, that most likely has functional significance in marine mussels. The identification of this ORF in five Mytilus species suggests that it has been maintained in the mytilid lineage (subfamily Mytilinae) for ∼13 million years. Furthermore, this ORF likely has a homologue in the F mt genome of Musculista senhousia, a DUI-containing mytilid species in the subfamily Crenellinae. We present evidence supporting the functionality of this F-specific ORF at the transcriptional, amino acid and nucleotide levels.

Conclusions/Significance

Our results offer support for the hypothesis that “novel F genome-specific mitochondrial genes” are involved in key biological functions in bivalve species with DUI.