Transcriptomics of mRNA and egg quality in farmed fish: Some recent developments and future directions.

Maternal mRNA transcripts deposited in growing oocytes regulate early development and are under intensive investigation as determinants of egg quality. The research has evolved from single gene studies to microarray and now RNA-Seq analyses in which mRNA expression by virtually every gene can be assessed and related to gamete quality. Such studies have mainly focused on genes changing two- to several-fold in expression between biological states, and have identified scores of candidate genes and a few gene networks whose functioning is related to successful development. However, ever-increasing yields of information from high throughput methods for detecting transcript abundance have far outpaced progress in methods for analyzing the massive quantities of gene expression data, and especially for meaningful relation of whole transcriptome profiles to gamete quality. We have developed a new approach to this problem employing artificial neural networks and supervised machine learning with other novel bioinformatics procedures to discover a previously unknown level of ovarian transcriptome function at which minute changes in expression of a few hundred genes is highly predictive of egg quality. In this paper, we briefly review the progress in transcriptomics of fish egg quality and discuss some future directions for this field of study.

Dynamics of the striped bass (Morone saxatilis) ovary proteome reveal a complex network of the translasome.

We evaluated changes in the striped bass (Morone saxatilis) ovary proteome during the annual reproductive cycle using label-free quantitative mass spectrometry and a novel machine learning analysis based on K-means clustering and support vector machines. Modulated modularity clustering was used to group co-variable proteins into expression modules and Gene Ontology (GO) biological process and KEGG pathway enrichment analyses were conducted for proteins within those modules. We discovered that components of the ribosome along with translation initiation and elongation factors generally decrease as the annual ovarian cycle progresses toward ovulation, concomitant with a slight increase in components of the 26S-proteasome. Co-variation within more than one expression module of components from these two multi-protein complexes suggests that they are not only co-regulated, but that co-regulation occurs through more than one sub-network. These components also co-vary with subunits of the TCP-1 chaperonin system and enzymes of intermediary metabolic pathways, suggesting that protein folding and cellular bioenergetic state play important roles in protein synthesis and degradation. We provide further evidence to suggest that protein synthesis and degradation are intimately linked, and our results support function of a proteasome-ribosome supercomplex known as the translasome.

Transcriptomic responses of juvenile Pacific whiteleg shrimp, Litopenaeus vannamei, to hypoxia and hypercapnic hypoxia.

Estuarine crustaceans are often exposed to low dissolved O2 (hypoxia) accompanied by elevated CO2 (hypercapnia), which lowers water pH. Acclimatory responses to hypoxia have been widely characterized; responses to hypercapnia in combination with hypoxia (hypercapnic hypoxia) are less well known. Here we used oligonucleotide microarrays to characterize changes in global gene expression in the hepatopancreas of Pacific whiteleg shrimp, Litopenaeus vannamei, exposed to hypoxia or hypercapnic hypoxia for 4 or 24 h, compared with time-matched animals held in air-saturated water (normoxia). Unigenes whose expressions were significantly impacted by treatment and/or time were used to build artificial neural networks (ANNs) to identify genes with the greatest sensitivity in pairwise discriminations between treatments at each time point and between times for each treatment. ANN gene sets that discriminated hypoxia or hypercapnic hypoxia from normoxia shared functions of translation, mitochondrial energetics, and cellular defense. GO terms protein modification/phosphorylation/cellular protein metabolism and RNA processing/apoptosis/cell cycling occurred at highest frequency in discriminating hypercapnic hypoxia from hypoxia at 4 and 24 h, respectively. For 75.4% of the annotated ANN genes, exposure to hypercapnic hypoxia for 24 h reduced or reversed the transcriptional response to hypoxia alone. These results suggest that high CO2/low pH may interfere with transcriptionally based acclimation to hypoxia or elicit physiological or biochemical responses that relieve internal hypoxia. Whether these data reflect resilience or sensitivity of L. vannamei in the face of expanding hypoxic zones and rising levels of atmospheric CO2 may be important to understanding the survival of this and other estuarine species.

Contaminants and impoSEX: transcriptomics of contaminant-induced sex change.

Reproduction is the goal of living organisms, and environmental conditions that influence sexual development are therefore critical to understanding adaptation in natural populations. It is not surprising that so much attention has been devoted to the impacts of the physical and chemical environment on this process (Vandenberg et al.2012). Chemicals of concern include a variety of endocrine disruptors (EDs) including oestrogen and oestrogen mimics that directly lead to malformation of the gonad. On the molecular side, the impact that EDs have on genes directly involved in the feminization or masculinization of the gonad such as Cyp 19A (or aromatase), foxl2,Sox9, Dmrt1 and NrOb1, has received considerable attention due to their direct involvement in the regulation of oestrogen and testosterone. In this issue of Molecular Ecology, Pascoal et al. (2013) examine the impact of a known endocrine disruptor (tributyltin or TBT) on the transcriptome of the dog whelk, Nucella lapillus (Fig. 1),in relation to the formation of imposex individuals (masculinized females). They conclude that TBT mimics the endogenous ligand of the nuclear retinoid X receptor (RXR) and/or peroxisome profilerator-activated receptor (PPAR) disrupting pathways.

Using machine learning tools to model complex toxic interactions with limited sampling regimes.

A major impediment to understanding the impact of environmental stress, including toxins and other pollutants, on organisms, is that organisms are rarely challenged by one or a few stressors in natural systems. Thus, linking laboratory experiments that are limited by practical considerations to a few stressors and a few levels of these stressors to real world conditions is constrained. In addition, while the existence of complex interactions among stressors can be identified by current statistical methods, these methods do not provide a means to construct mathematical models of these interactions. In this paper, we offer a two-step process by which complex interactions of stressors on biological systems can be modeled in an experimental design that is within the limits of practicality. We begin with the notion that environment conditions circumscribe an n-dimensional hyperspace within which biological processes or end points are embedded. We then randomly sample this hyperspace to establish experimental conditions that span the range of the relevant parameters and conduct the experiment(s) based upon these selected conditions. Models of the complex interactions of the parameters are then extracted using machine learning tools, specifically artificial neural networks. This approach can rapidly generate highly accurate models of biological responses to complex interactions among environmentally relevant toxins, identify critical subspaces where nonlinear responses exist, and provide an expedient means of designing traditional experiments to test the impact of complex mixtures on biological responses. Further, this can be accomplished with an astonishingly small sample size.

The transcriptomic responses of the eastern oyster, Crassostrea virginica, to environmental conditions.

Understanding the mechanisms by which organisms adapt to environmental conditions is a fundamental question for ecology and evolution. In this study, we evaluate changes in gene expression of a marine mollusc, the eastern oyster Crassostrea virginica, associated with the physico-chemical conditions and the levels of metals and other contaminants in their environment. The results indicate that transcript signatures can effectively disentangle the complex interactive gene expression responses to the environment and are also capable of disentangling the complex dynamic effects of environmental factors on gene expression. In this context, the mapping of environment to gene and gene to environment is reciprocal and mutually reinforcing. In general, the response of transcripts to the environment is driven by major factors known to affect oyster physiology such as temperature, pH, salinity, and dissolved oxygen, with pollutant levels playing a relatively small role, at least within the range of concentrations found in the studied oyster habitats. Further, the two environmental factors that dominate these effects (temperature and pH) interact in a dynamic and nonlinear fashion to impact gene expression. Transcriptomic data obtained in our study provide insights into the mechanisms of physiological responses to temperature and pH in oysters that are consistent with the known effects of these factors on physiological functions of ectotherms and indicate important linkages between transcriptomics and physiological outcomes. Should these linkages hold in further studies and in other organisms, they may provide a novel integrated approach for assessing the impacts of climate change, ocean acidification and anthropogenic contaminants on aquatic organisms via relatively inexpensive microarray platforms.

Modelling interactions of acid-base balance and respiratory status in the toxicity of metal mixtures in the American oyster Crassostrea virginica.

Heavy metals, such as copper, zinc and cadmium, represent some of the most common and serious pollutants in coastal estuaries. In the present study, we used a combination of linear and artificial neural network (ANN) modelling to detect and explore interactions among low-dose mixtures of these heavy metals and their impacts on fundamental physiological processes in tissues of the Eastern oyster, Crassostrea virginica. Animals were exposed to Cd (0.001-0.400 microM), Zn (0.001-3.059 microM) or Cu (0.002-0.787 microM), either alone or in combination for 1 to 27 days. We measured indicators of acid-base balance (hemolymph pH and total CO(2)), gas exchange (Po(2)), immunocompetence (total hemocyte counts, numbers of invasive bacteria), antioxidant status (glutathione, GSH), oxidative damage (lipid peroxidation; LPx), and metal accumulation in the gill and the hepatopancreas. Linear analysis showed that oxidative membrane damage from tissue accumulation of environmental metals was correlated with impaired acid-base balance in oysters. ANN analysis revealed interactions of metals with hemolymph acid-base chemistry in predicting oxidative damage that were not evident from linear analyses. These results highlight the usefulness of machine learning approaches, such as ANNs, for improving our ability to recognize and understand the effects of sub-acute exposure to contaminant mixtures.

A transcriptomic analysis of land-use impacts on the oyster, Crassostrea virginica, in the South Atlantic bight.

Increasing utilization and human population density in the coastal zone is widely believed to place increasing stresses on the resident biota, but confirmation of this belief is somewhat lacking. While we have solid evidence that highly disturbed estuarine systems have dramatic changes in the resident biota (black and white if you will), we lack tools that distinguish the shades of grey. In part, this lack of ability to distinguish shades of grey stems from the analytical tools that have been applied to studies of estuarine systems, and perhaps more important, is the insensitivity of the biological end points that we have used to assess these impacts. In this study, we will present data on the phenotypic adjustments as measured by transcriptomic signatures of a resilient organism (oysters) to land-use practices in the surrounding watershed using advanced machine-learning algorithms. We will demonstrate that such an approach can reveal subtle and meaningful shifts in oyster gene expression in response to land use. Further, the data show that gill tissues are far more responsive and provide superior discrimination of land-use classes than hepatopancreas and that transcripts encoding proteins involved in energy production, protein synthesis and basic metabolism are more robust indicators of land use than classic biomarkers such as metallothioneins, GST and cytochrome P-450.

Contributions of functional genomics and proteomics to the study of immune responses in the Pacific white leg shrimp Litopenaeus vannamei.

The need for better control of infectious diseases in shrimp aquaculture and the ecological importance of crustacea in marine ecosystems have prompted interest in the study of crustacean immune systems, particularly those of shrimp. As shrimp and other crustacea are poorly understood from the immunological point of view, functional genomic and proteomic approaches have been applied as a means of quickly obtaining molecular information regarding immune responses in these organisms. In this article, a series of results derived from transcriptomic and proteomic studies in shrimp (Litopenaeus vannamei) are discussed. Expressed Sequence Tag analysis, differential expression cloning through Suppression Subtractive Hybridization, expression profiling using microarrays, and proteomic studies using mass spectrometry, have provided a wealth of useful data and opportunities for new avenues of research. Examples of new research directions arising from these studies in shrimp include the molecular diversity of antimicrobial effectors, the role of double stranded RNA as an inducer of antiviral immunity, and the possible overlap between antibacterial and antiviral responses in the shrimp.

Insights into the immune transcriptome of the shrimp Litopenaeus vannamei: tissue-specific expression profiles and transcriptomic responses to immune challenge.

Infectious disease constitutes a major obstacle to the sustainability of shrimp aquaculture worldwide and a significant threat to natural populations of shrimp and other crustacea. The study of the shrimp immune system, including the response to viral infection, has been hampered by a relative lack of molecular genetic information and of tools suitable for high-throughput assessment of gene expression. In this report, the generation of a cDNA microarray encompassing 2,469 putative unigenes expressed in gills, circulating hemocytes, and hepatopancreas of Litopenaeus vannamei is described. The unigenes printed on the microarray were derived from the analyses of 7,021 expressed sequence tags obtained from standard cDNA libraries as well as from libraries generated by suppression subtractive hybridization, after challenging shrimp with a variety of immune stimuli. The general utility of the cDNA microarray was demonstrated by interrogating the array with labeled RNA from four different shrimp tissues (gills, hemocytes, hepatopancreas, and muscle) and by analyzing the transcriptomic response of shrimp to a lethal challenge with white spot syndrome virus. Our results indicate that white spot syndrome virus infection upregulates (in the hepatopancreas) genes encoding known and potential antimicrobial effectors, while some genes involved in protection from oxidative stress were found to be downregulated by the virus.

Double-stranded RNA and antiviral immunity in marine shrimp: inducible host mechanisms and evidence for the evolution of viral counter-responses.

Double-stranded RNA (dsRNA) is a common virus-associated molecular pattern and a potent inducer of antiviral responses in many organisms. While it is clear that the specific RNA interference (RNAi) response, a phenomenon triggered by dsRNA, serves antiviral functions in invertebrates, innate (non-specific) antiviral immune reactions induced by dsRNA (e.g. the Interferon response) have long been thought to be restricted to vertebrates. Recent work in an underappreciated experimental model, the penaeid shrimp, is challenging these traditional distinctions, by demonstrating the existence of both innate (non sequence-specific) and RNAi-related (sequence-specific) antiviral phenomena in crustacea. Here we discuss the evidence for this bivalent role of dsRNA in the initiation of antiviral responses in shrimp, and present new data that suggest that the antiviral functions of the shrimp RNAi machinery have imposed selective pressures on an evolving viral pathogen. These findings open the door for the discovery of novel mechanisms of innate immunity, and provide a basis for the future development of strategies to control viral diseases in the commercially important penaeid shrimp.

A dolphin peripheral blood leukocyte cDNA microarray for studies of immune function and stress reactions.

A microarray focused on stress response and immune function genes of the bottlenosed dolphin has been developed. Random expressed sequence tags (ESTs) were isolated and sequenced from two dolphin peripheral blood leukocyte (PBL) cDNA libraries biased towards T- and B-cell gene expression by stimulation with IL-2 and LPS, respectively. A total of 2784 clones were sequenced and contig analysis yielded 1343 unigenes (archived and annotated at ). In addition, 52 dolphin genes known to be important in innate and adaptive immune function and stress responses of terrestrial mammals were specifically targeted, cloned and added to the unigene collection. The set of dolphin sequences printed on a cDNA microarray comprised the 1343 unigenes, the 52 targeted genes and 2305 randomly selected (but unsequenced) EST clones. This set was printed in duplicate spots, side by side, and in two replicates per slide, such that the total number of features per microarray slide was 19,200, including controls. The dolphin arrays were validated and transcriptomic profiles were generated using PBL from a wild dolphin, a captive dolphin and dolphin skin cells. The results demonstrate that the array is a reproducible and informative tool for assessing differential gene expression in dolphin PBL and in other tissues.

Characterization of the immunoglobulin A heavy chain gene of the Atlantic bottlenose dolphin (Tursiops truncatus).

Immunoglobulin constant region heavy chain genes of the dolphin (Tursiops truncatus) have been described for IgM and IgG but not for IgA. Here, the heavy chain sequence of dolphin IgA has been cloned and sequenced as cDNA. RT-PCR amplification from blood peripheral lymphocytes was carried out using degenerate primers and a single sequence was detected. The inferred heavy chain structure shows conserved features typical of mammalian IgA heavy chains, including three constant (C) regions, a hinge region between constant region domain 1 (C1) and constant region domain 2 (C2), and conserved residues for interaction with the Fc alpha R1 and N-glycosylation sites. Comparisons of the deduced amino acid sequences of the IgA heavy chain for the dolphin and the evolutionarily related artiodactyl species showed high similarity. In cattle and sheep, as in dolphins, a single IgA subclass has been identified. Southern blot analysis as well as genomic PCR confirmed the presence of multiple IGHA sequences suggesting that IGHA pseudogenes may be present in the dolphin genome.

A cDNA microarray for Crassostrea virginica and C. gigas.

The eastern oyster, Crassostrea virginica, and the Pacific oyster, C. gigas, are species of global economic significance as well as important components of estuarine ecosystems and models for genetic and environmental studies. To enhance the molecular tools available for oyster research, an international group of collaborators has constructed a 27,496-feature cDNA microarray containing 4460 sequences derived from C. virginica, 2320 from C. gigas, and 16 non-oyster DNAs serving as positive and negative controls. The performance of the array was assessed by gene expression profiling using gill and digestive gland RNA derived from both C. gigas and C. virginica, and digestive gland RNA from C. ariakensis. The utility of the microarray for detection of homologous genes by cross-hybridization between species was also assessed and the correlation between hybridization intensity and sequence homology for selected genes determined. The oyster cDNA microarray is publicly available to the research community on a cost-recovery basis.

Regulation of metallothionein genes in the American oyster (Crassostrea virginica): ontogeny and differential expression in response to different stressors.

Metallothioneins (MTs) are typically low molecular weight (6-7 kDa), metal-binding proteins with characteristic repeating cysteine motifs (Cys-X-Cys or Cys-Xn-Cys) and a prolate ellipsoid shape containing single alpha- and beta-domains. While functionally diverse, they play important roles in metals homeostasis, detoxification and the stress response. The present study, combined with previous observations (e.g., Jenny et al., Eur. J. Biochem. 2005; 271:1702-1712) defines an unprecedented diversity of MT primary structure and domain organization in the American oyster, Crassostrea virginica. Two novel molluscan MT families are described. One of these (CvMT-III) is characterized by the presence of two beta-domains and the absence of alpha-domains. This family exhibits constitutive expression during larval development and is the dominant CvMT isoform expressed in larvae. CvMT-III displays low basal levels of expression in adult tissues and only moderate responsiveness to metal challenges in both larvae and adults. A second novel MT isoform (CvMT-IV) was isolated from hemocytes by subtractive hybridization techniques following a 4-hour immune challenge with heat-killed bacteria (Vibrio, Bacillus, Micrococcus spp. mixture). Based on conservation of the cysteine motifs, this isoform appears to be a sub-family related to the molluscan alphabeta-domain MTs. A series of amino acid substitutions has resulted in four additional cysteines which give rise to a Cys-Cys motif and three Cys-Cys-Cys motifs. Northern blot analyses demonstrate that CvMT-IV is down-regulated upon sterile wounding and immune challenge, displays moderate expression in larvae and adults and differential gene induction in response to metals exposure.

New resources for marine genomics: bacterial artificial chromosome libraries for the Eastern and Pacific oysters (Crassostrea virginica and C. gigas).

Large-insert genomic bacterial artificial chromosome (BAC) libraries of two culturally and economically important oyster species, Crassostrea virginica and C. gigas, have been developed as part of an international effort to develop tools and reagents that will advance our ability to conduct genetic and genomic research. A total of 73,728 C. gigas clones with an average insert size of 152 kb were picked and arrayed representing an 11.8-fold genome coverage. A total of 55,296 clones with an average insert size of 150 kb were picked and arrayed for C. virginica, also representing an 11.8-fold genome coverage. The C. gigas and C. virginica libraries were screened with probes derived from selected oyster genes using high-density BAC colony filter arrays. The probes identified 4 to 25 clones per gene for C. virginica and 5 to 50 clones per gene for C. gigas. We conducted a preliminary analysis of genetic polymorphism represented in the C. gigas library. The results suggest that the degree of divergence among similar sequences is highly variable and concentrated in intronic regions. Evidence supporting allelic polymorphism is reported for two genes and allelic and/or locus specific polymorphism for several others. Classical inheritance studies are needed to confirm the nature of these polymorphisms. The oyster BAC libraries are publicly available to the research community on a cost-recovery basis at (www.genome.clemson.edu).

The Immunoglobulin G Heavy Chain (IGHG) genes of the Atlantic bottlenose dolphin, Tursiops truncatus.

Dolphin Immunoglobulin G Heavy Chain (IGHG) sequences were obtained by PCR amplification of cDNA from peripheral blood leukocytes using degenerate primers. Analysis of full-length sequences indicated the presence of two expressed isotypes, IGHG1 and IGHG2 that differ mainly in the hinge region of the molecule. Genomic Southern blot analysis indicated that the IGHG1 and IGHG2 genes are most likely present in single copies. The inferred amino acid sequences show greatest similarity between the dolphin and other closely related artiodactyl species. The genetic structure of the IGHG genes were deduced through genomic PCR and revealed that the hinge regions of both IGHG1 and IGHG2 are encoded by a single exon. The transmembrane region of the dolphin IGHG chain shows similarity to the transmembrane region of other mammalian IGHG chains with a canonical CART motif. This is in contrast to the unusual Ser to Gly substitution previously found in the dolphin IGHM transmembrane region, and the functional significance of this variation for B cell antigen-receptor dimer activation remains unknown.

Predictive Blood Chemistry Parameters for Pansteatitis-Affected Mozambique Tilapia (Oreochromis mossambicus).

One of the largest river systems in South Africa, the Olifants River, has experienced significant changes in water quality due to anthropogenic activities. Since 2005, there have been various "outbreaks" of the inflammatory disease pansteatitis in several vertebrate species. Large-scale pansteatitis-related mortality events have decimated the crocodile population at Lake Loskop and decreased the population at Kruger National Park. Most pansteatitis-related diagnoses within the region are conducted post-mortem by either gross pathology or histology. The application of a non-lethal approach to assess the prevalence and pervasiveness of pansteatitis in the Olifants River region would be of great importance for the development of a management plan for this disease. In this study, several plasma-based biomarkers accurately classified pansteatitis in Mozambique tilapia (Oreochromis mossambicus) collected from Lake Loskop using a commercially available benchtop blood chemistry analyzer combined with data interpretation via artificial neural network analysis. According to the model, four blood chemistry parameters (calcium, sodium, total protein and albumin), in combination with total length, diagnose pansteatitis to a predictive accuracy of 92 percent. In addition, several morphometric traits (total length, age, weight) were also associated with pansteatitis. On-going research will focus on further evaluating the use of blood chemistry to classify pansteatitis across different species, trophic levels, and within different sites along the Olifants River.

Marine genomics: a clearing-house for genomic and transcriptomic data of marine organisms.

BACKGROUND: The Marine Genomics project is a functional genomics initiative developed to provide a pipeline for the curation of Expressed Sequence Tags (ESTs) and gene expression microarray data for marine organisms. It provides a unique clearing-house for marine specific EST and microarray data and is currently available at http://www.marinegenomics.org. DESCRIPTION: The Marine Genomics pipeline automates the processing, maintenance, storage and analysis of EST and microarray data for an increasing number of marine species. It currently contains 19 species databases (over 46,000 EST sequences) that are maintained by registered users from local and remote locations in Europe and South America in addition to the USA. A collection of analysis tools are implemented. These include a pipeline upload tool for EST FASTA file, sequence trace file and microarray data, an annotative text search, automated sequence trimming, sequence quality control (QA/QC) editing, sequence BLAST capabilities and a tool for interactive submission to GenBank. Another feature of this resource is the integration with a scientific computing analysis environment implemented by MATLAB. CONCLUSION: The conglomeration of multiple marine organisms with integrated analysis tools enables users to focus on the comprehensive descriptions of transcriptomic responses to typical marine stresses. This cross species data comparison and integration enables users to contain their research within a marine-oriented data management and analysis environment.