Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Genome-wide association study for Streptococcus iniae in Nile tilapia (Oreochromis niloticus) identifies a significant QTL for disease resistance

Streptococcus iniae is a problematic gram-positive bacterium negatively affecting Nile tilapia (Oreochromis niloticus), one of the main aquacultural species produced worldwide. The aim of this study was to identify the genetic architecture of survival to S. iniae and identify single nucleotide polymorphism (SNPs) linked to quantitative trait loci (QTL) related to survival to S. iniae challenge. With this purpose, Nile tilapia from the Spring Genetics breeding program were sent to a controlled S. iniae challenge test where phenotypes were scored as dead for fish that died during challenge test and survivors for the fish alive at the termination of the test. Additionally, fin-clip samples from all fish in the test were collected for DNA extraction. Out of 1904 fish in the challenge test, tissue samples of 321 fish were sent for genotyping using double digest restriction site associated DNA sequencing (ddRADseq). After quality control and filtering, 9,085 SNPs were used to perform a genome-wide association study (GWAS). A significant signal in LG8 was observed indicating association with survival to S. iniae challenge, with SNPs explaining from 12% to 26% of the genetic variance. To demonstrate the usefulness of marker assisted selection (MAS) to selectively breed fish for survival to S. iniae, offspring of breeding candidates classified as "resistant" and "susceptible" based on haplotypes of the four most significant markers were sent to a controlled S. iniae challenge test. At the end of the test, the differences in mortality between the two groups were strikingly different with a final cumulative percent mortality of less than 1% and 73% for offspring from "resistant" and "susceptible" parents, respectively. These results demonstrate that MAS for improved resistance to S. iniae is feasible.

LSD1 represses a neonatal/reparative gene program in adult intestinal epithelium

Intestinal epithelial homeostasis is maintained by adult intestinal stem cells, which, alongside Paneth cells, appear after birth in the neonatal period. We aimed to identify regulators of neonatal intestinal epithelial development by testing a small library of epigenetic modifier inhibitors in Paneth cell-skewed organoid cultures. We found that lysine-specific demethylase 1A (Kdm1a/Lsd1) is absolutely required for Paneth cell differentiation. Lsd1-deficient crypts, devoid of Paneth cells, are still able to form organoids without a requirement of exogenous or endogenous Wnt. Mechanistically, we find that LSD1 enzymatically represses genes that are normally expressed only in fetal and neonatal epithelium. This gene profile is similar to what is seen in repairing epithelium, and we find that Lsd1-deficient epithelium has superior regenerative capacities after irradiation injury. In summary, we found an important regulator of neonatal intestinal development and identified a druggable target to reprogram intestinal epithelium toward a reparative state.

Differential Effects of Dietary Supplementation of Krill Meal, Soybean Meal, Butyrate, and Bactocell® on the Gene Expression of Atlantic Salmon Head Kidney

The head kidney is a key organ that plays a fundamental role in the regulation of the fish immune response and in the maintenance of endocrine homeostasis. Previous studies indicate that the supplementation of exogenous dietary components, such as krill meal (KM), soybean meal (SM), Bactocell® (BA), and butyrate (BU), can have a significant effect on the immune function of the head kidney. The aim of this study was to investigate the differential effect of these four dietary ingredients on the transcriptional profiles of the head kidney of the Atlantic salmon. This study revealed that just a small number of genes were responsive to the feeding regime after a long-term (12 weeks) treatment, and evidenced that the most significant alterations, both in terms of the number of affected genes and magnitude of changes in gene expression, were detectable in the BU- and KM-fed groups compared with controls, while the SM diet had a nearly negligible effect, and BA had no significant effects at all. Most of the differentially expressed genes were involved in the immune response and, in line with data previously obtained from pyloric caeca, major components of the complement system were significantly affected. These alterations were accompanied by an increase in the density of melanomacrophage centers in the KM- and SM-fed group and their reduction in the BU-fed group. While three types of dietary supplements (BU, KM, and SM) were able to produce a significant modulation of some molecular players of the immune system, the butyrate-rich diet was revealed as the one with the most relevant immune-stimulating properties in the head kidney. These preliminary results suggest that further investigations should be aimed towards the elucidation of the potential beneficial effects of butyrate and krill meal supplementation on farmed salmon health and growth performance.

Selection response for Streptococcus agalactiae resistance in Nile tilapia Oreochromis niloticus

The potential of selection to improve resistance to streptococcosis was evaluated in a commercial population of Nile tilapia in Thailand. The base generation (G0) consisted of offspring from 98 sires and 149 dams using a partly nested design. At 60 days post-hatch, 30 fish from each family were injected intraperitoneally with a Streptococcosis agalactiae solution (1 × 10⁹  CFU/ml) and evaluated for 14 days. Disease resistance was recorded as the number of days from challenge until death (DD) and as a binary (BIN) trait (dead/alive) on day 14. Three models were used for genetic analyses: Cox frailty model for DD; animal model for DD; and animal model for BIN. Age at challenge was fitted as a covariate and contemporary group as fixed or random effect, depending on the model. Fish from the 18 most resistant families were selected to produce the first generation (G1). Heritability estimates for G0 were 0.22, 0.14 ± 0.02 and 0.11 ± 0.02 for the Cox, linear DD and linear BIN models, respectively. Selection response indicated that the risk of death decreased to 54%, survival time increased to 3.4 days and survival rate increased to 21%. These results suggest that genetic improvement is possible for this population.

Novel BRAF and KRAS Mutations in Papillary Thyroid Carcinoma Arising in Struma Ovarii

Papillary carcinomas of thyroid type rarely arise within struma ovarii. There are limited data on the immunohistochemical and molecular features of these tumors. Three cases of papillary carcinoma arising in struma ovarii (PCSO) were identified. The clinicopathological features were reviewed and immunohistochemical staining for HBME-1, cytokeratin (CK) 19, and CD56 was performed. Tumor DNA was sequenced for somatic mutations using a panel of 26 oncogenes, with a particular focus on BRAF and KRAS mutations. The patients were aged 22, 48, and 55 years. All cases were FIGO stage IA. Two tumors were of classical histological type, and one was a follicular variant papillary carcinoma. All tumors expressed HBME-1 and two were positive for CK19. CD56 was negative in all three cases. One tumor demonstrated a BRAF G469A mutation in exon 11, and in a second case, a KRAS Q61K double base mutation in exon 3 was detected. These mutations have not been described previously in PCSO. No mutations were detected in the benign follicular components of the tumors adjacent to the malignant papillary tissue. None of the patients had tumor recurrence on clinical follow-up (range 11 months to 8½ years). HBME-1, CK19, and CD56 are useful immunohistochemical markers of PCSO. Novel BRAF and KRAS mutations were identified in two of three tumors suggesting that mutations in PCSO may differ from those commonly identified in papillary carcinoma of the eutopic thyroid. The clinical significance of these mutations is uncertain but follow-up data in this small series support the generally good prognosis of PCSO.

ClusTrack: feature extraction and similarity measures for clustering of genome-wide data sets

Clustering is a popular technique for explorative analysis of data, as it can reveal subgroupings and similarities between data in an unsupervised manner. While clustering is routinely applied to gene expression data, there is a lack of appropriate general methodology for clustering of sequence-level genomic and epigenomic data, e.g. ChIP-based data. We here introduce a general methodology for clustering data sets of coordinates relative to a genome assembly, i.e. genomic tracks. By defining appropriate feature extraction approaches and similarity measures, we allow biologically meaningful clustering to be performed for genomic tracks using standard clustering algorithms. An implementation of the methodology is provided through a tool, ClusTrack, which allows fine-tuned clustering analyses to be specified through a web-based interface. We apply our methods to the clustering of occupancy of the H3K4me1 histone modification in samples from a range of different cell types. The majority of samples form meaningful subclusters, confirming that the definitions of features and similarity capture biological, rather than technical, variation between the genomic tracks. Input data and results are available, and can be reproduced, through a Galaxy Pages document at http://hyperbrowser.uio.no/hb/u/hb-superuser/p/clustrack. The clustering functionality is available as a Galaxy tool, under the menu option "Specialized analyzis of tracks", and the submenu option "Cluster tracks based on genome level similarity", at the Genomic HyperBrowser server: http://hyperbrowser.uio.no/hb/.

Chromatin states reveal functional associations for globally defined transcription start sites in four human cell lines

BACKGROUND

Deciphering the most common modes by which chromatin regulates transcription, and how this is related to cellular status and processes is an important task for improving our understanding of human cellular biology. The FANTOM5 and ENCODE projects represent two independent large scale efforts to map regulatory and transcriptional features to the human genome. Here we investigate chromatin features around a comprehensive set of transcription start sites in four cell lines by integrating data from these two projects.

RESULTS

Transcription start sites can be distinguished by chromatin states defined by specific combinations of both chromatin mark enrichment and the profile shapes of these chromatin marks. The observed patterns can be associated with cellular functions and processes, and they also show association with expression level, location relative to nearby genes, and CpG content. In particular we find a substantial number of repressed inter- and intra-genic transcription start sites enriched for active chromatin marks and Pol II, and these sites are strongly associated with immediate-early response processes and cell signaling. Associations between start sites with similar chromatin patterns are validated by significant correlations in their global expression profiles.

CONCLUSIONS

The results confirm the link between chromatin state and cellular function for expressed transcripts, and also indicate that active chromatin states at repressed transcripts may poise transcripts for rapid activation during immune response.

Identification of serum microRNA profiles in colon cancer

Background:

microRNAs (miRNAs) exist in blood in an apparently stable form. We have explored whether serum miRNAs can be used as non-invasive early biomarkers of colon cancer.

Methods:

Serum samples from 30 patients with colon cancer stage IV and 10 healthy controls were examined for the expression of 375 cancer-relevant miRNAs. Based on the miRNA profile in this study, 34 selected miRNAs were measured in serum from 40 patients with stage I–II colon cancer and from 10 additional controls.

Results:

Twenty miRNAs were differentially expressed in serum from stage IV patients compared with controls (P<0.01). Unsupervised clustering revealed four subgroups; one corresponding mostly to the control group and the three others to the patient groups. Of the 34 miRNAs measured in the follow-up study of stage I–II patients, 21 showed concordant expression between stage IV and stage I–II patient. Based on the profiles of these 21 miRNAs, a supervised linear regression analysis (Partial Least Squares Regression) was performed. Using this model we correctly assigned stage I–II colon cancer patients based on miRNA profiles of stage IV patients.

Conclusion:

Serum miRNA expression profiling may be utilised in early detection of colon cancer.

Survival, growth and sexual maturation in Atlantic salmon exposed to infectious pancreatic necrosis: a multi-variate mixture model approach

Background

Outbreaks of infectious pancreatic necrosis (IPN) in Atlantic salmon can result in reduced growth rates in a fraction of the surviving fish (runts). Genetic and environmental variation also affects growth rates within different categories of healthy animals and runts, which complicates identification of runts. Mixture models are commonly used to identify the underlying structures in such data, and the aim of this study was to develop Bayesian mixture models for the genetic analysis of health status (runt/healthy) of surviving fish from an IPN outbreak.

Methods

Five statistical models were tested on data consisting of 10 972 fish that died and 3959 survivors with recorded growth data. The most complex models (4 and 5) were multivariate normal-binary mixture models including growth, sexual maturity and field survival traits. Growth rate and liability of sexual maturation were treated as two-component normal mixtures, assuming phenotypes originated from two potentially overlapping distributions, (runt/normal). Runt status was an unobserved binary trait. These models were compared to mixture models with fewer traits (Models 2 and 3) and a classical linear animal model for growth (Model 1).

Results

Assuming growth as a mixture trait improved the predictive ability of the statistical model considerably (Model 2 vs. 1). The final models (4 and 5) yielded the following results: estimated (underlying) heritabilities were moderate for growth in healthy fish (0.32 ± 0.04 and 0.35 ± 0.05), runt status (0.39 ± 0.07 and 0.36 ± 0.08) and sexual maturation (0.33 ± 0.05), and high for field survival (0.47 ± 0.03 and 0.48 ± 0.03). Growth in healthy animals, runt status and survival showed consistent favourable genetic associations. Sexual maturation showed an unfavourable non-significant genetic correlation with runt status, but favourable genetic correlations with other traits. The estimated fraction of healthy fish was 81-85%. The estimated breeding values for runt status and (normal) growth were consistent for the most complex models (4 and 5), but showed imperfect correlations with estimated breeding values from the simpler models.

Conclusions

Modelling growth in IPN survivors as a mixture trait improved the predictive ability of the model compared with a classical linear model. The results indicated considerable genetic variation in health status among survivors. Mixture modelling may be useful for the genetic analysis of diseases detected mainly through indicator traits.

Cell-type specificity of ChIP-predicted transcription factor binding sites

BACKGROUND

Context-dependent transcription factor (TF) binding is one reason for differences in gene expression patterns between different cellular states. Chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) identifies genome-wide TF binding sites for one particular context-the cells used in the experiment. But can such ChIP-seq data predict TF binding in other cellular contexts and is it possible to distinguish context-dependent from ubiquitous TF binding?

RESULTS

We compared ChIP-seq data on TF binding for multiple TFs in two different cell types and found that on average only a third of ChIP-seq peak regions are common to both cell types. Expectedly, common peaks occur more frequently in certain genomic contexts, such as CpG-rich promoters, whereas chromatin differences characterize cell-type specific TF binding. We also find, however, that genotype differences between the cell types can explain differences in binding. Moreover, ChIP-seq signal intensity and peak clustering are the strongest predictors of common peaks. Compared with strong peaks located in regions containing peaks for multiple transcription factors, weak and isolated peaks are less common between the cell types and are less associated with data that indicate regulatory activity.

CONCLUSIONS

Together, the results suggest that experimental noise is prevalent among weak peaks, whereas strong and clustered peaks represent high-confidence binding events that often occur in other cellular contexts. Nevertheless, 30-40% of the strongest and most clustered peaks show context-dependent regulation. We show that by combining signal intensity with additional data-ranging from context independent information such as binding site conservation and position weight matrix scores to context dependent chromatin structure-we can predict whether a ChIP-seq peak is likely to be present in other cellular contexts.

Clustered ChIP-Seq-defined transcription factor binding sites and histone modifications map distinct classes of regulatory elements

Background

Transcription factor binding to DNA requires both an appropriate binding element and suitably open chromatin, which together help to define regulatory elements within the genome. Current methods of identifying regulatory elements, such as promoters or enhancers, typically rely on sequence conservation, existing gene annotations or specific marks, such as histone modifications and p300 binding methods, each of which has its own biases.

Results

Herein we show that an approach based on clustering of transcription factor peaks from high-throughput sequencing coupled with chromatin immunoprecipitation (Chip-Seq) can be used to evaluate markers for regulatory elements. We used 67 data sets for 54 unique transcription factors distributed over two cell lines to create regulatory element clusters. By integrating the clusters from our approach with histone modifications and data for open chromatin, we identified general methylation of lysine 4 on histone H3 (H3K4me) as the most specific marker for transcription factor clusters. Clusters mapping to annotated genes showed distinct patterns in cluster composition related to gene expression and histone modifications. Clusters mapping to intergenic regions fall into two groups either directly involved in transcription, including miRNAs and long noncoding RNAs, or facilitating transcription by long-range interactions. The latter clusters were specifically enriched with H3K4me1, but less with acetylation of lysine 27 on histone 3 or p300 binding.

Conclusion

By integrating genomewide data of transcription factor binding and chromatin structure and using our data-driven approach, we pinpointed the chromatin marks that best explain transcription factor association with different regulatory elements. Our results also indicate that a modest selection of transcription factors may be sufficient to map most regulatory elements in the human genome.

Selection for improved resistance to Aeromonas hydrophila in Indian major carp Labeo rohita: survival and innate immune responses in first generation of resistant and susceptible lines

Selection for disease resistance in fish may be performed directly on basis of survival data obtained in controlled challenge trials, or indirectly using information from immunological or molecular markers linked to differential survival. In the present study, several key innate immune parameters were measured in aeromoniasis resistant and susceptible lines of rohu Labeo rohita to assess their suitability as immune markers for use in indirect selection for increased resistance. Experimental infection with Aeromonas hydrophila (9.55 × 10(6) cfu g(-1) fish) through the intraperitoneal route produced higher survival in the resistant line (73.33%) as compared to the susceptible line (16.67%). Blood and liver tissue samples from both lines were collected to study some of the innate immune parameters and immune-related gene expression. The respiratory burst activity of blood phagocytes, serum myeloperoxidase activity and ceruloplasmin level were significantly (p < 0.05) higher in the resistant line compared to the susceptible line. Lower level of blood glucose and serum natural haemolysin titre were marked in the resistant line as compared to the susceptible line. No significant difference was measured in total serum protein concentration, antiprotease activity and bacterial agglutinin level between two lines, while the expression of transferrin, complement factor C3 and TLR 22-like transcripts were significantly (P < 0.05) higher in liver samples of the susceptible line. However, no such difference was found in β(2)-microglobulin and lysozyme gene expression between lines. The study demonstrated the possibility of using some of the investigated innate immune parameters as indirect marker traits for selection for improved resistance to aeromoniasis in rohu.

Quantitative genetics of Taura syndrome resistance in Pacific white shrimp (Penaeus vannamei): a cure model approach

Background

In aquaculture breeding, resistance against infectious diseases is commonly assessed as time until death under exposure to a pathogen. For some diseases, a fraction of the individuals may appear as "cured" (non-susceptible), and the resulting survival time may thus be a result of two confounded underlying traits, i.e., endurance (individual hazard) and susceptibility (whether at risk or not), which may be accounted for by fitting a cure survival model. We applied a cure model to survival data of Pacific white shrimp (Penaeus vannamei) challenged with the Taura syndrome virus, which is one of the major pathogens of Panaeid shrimp species.

Methods

In total, 15,261 individuals of 513 full-sib families from three generations were challenge-tested in 21 separate tests (tanks). All challenge-tests were run until mortality naturally ceased. Time-until-event data were analyzed with a mixed cure survival model using Gibbs sampling, treating susceptibility and endurance as separate genetic traits.

Results

Overall mortality at the end of test was 28%, while 38% of the population was considered susceptible to the disease. The estimated underlying heritability was high for susceptibility (0.41 ± 0.07), but low for endurance (0.07 ± 0.03). Furthermore, endurance and susceptibility were distinct genetic traits (r_g = 0.22 ± 0.25). Estimated breeding values for endurance and susceptibility were only moderately correlated (0.50), while estimated breeding values from classical models for analysis of challenge-test survival (ignoring the cured fraction) were closely correlated with estimated breeding values for susceptibility, but less correlated with estimated breeding values for endurance.

Conclusions

For Taura syndrome resistance, endurance and susceptibility are apparently distinct genetic traits. However, genetic evaluation of susceptibility based on the cure model showed clear associations with standard genetic evaluations that ignore the cure fraction for these data. Using the current testing design, genetic variation in observed survival time and absolute survival at the end of test were most likely dominated by genetic variation in susceptibility. If the aim is to reduce susceptibility, earlier termination of the challenge-test or back-truncation of the follow-up period should be avoided, as this may shift focus of selection towards endurance rather than susceptibility.

Family association between immune parameters and resistance to Aeromonas hydrophila infection in the Indian major carp, Labeo rohita

Seven innate immune parameters were investigated in 64 full-sib families (the offspring of 64 sires and 45 dams) from two year-classes of farmed rohu carp (Labeo rohita). Survival rates were also available from Aeromonas hydrophila infection (aeromoniasis) recorded in controlled challenge tests on a different sample of individuals from the same families. Due to strong confounding between the animal additive genetic effect and the family effects (common environmental+non-additive genetic), reliable additive (co)variance components and hence heritabilities and genetic correlations could not be obtained for the investigated parameters. Therefore, estimates of the association of challenge test survival with the studied immune parameters were obtained as product moment correlations between family least square means. These correlations revealed statistically significant (p<0.05) negative correlations of survival with bacterial agglutination titre (-0.48), serum haemolysin titre (-0.29) and haemagglutination titre (-0.34); and significant positive correlation with ceruloplasmin level (0.51). The correlations of survival to aeromoniasis with myeloperoxidase activity, superoxide production and lysozyme activity were found to be not significantly different from zero (p>0.05). Assuming that the negatively correlated candidate traits are not favourable as indirect selection criteria, the results suggest that ceruloplasmin level could potentially be a marker for resistance to aeromoniasis in rohu. The use of this immune parameter as an indirect selection criterion for increased resistance to aeromoniasis in rohu will, however, require that the parameter shows significant additive genetic variation and a significant genetic correlation with survival. Further studies are therefore needed to obtain a reliable heritability estimate for ceruloplasmin and its genetic correlation with survival from aeromoniasis.

Pixel-based analysis of multiple images for the identification of changes: a novel approach applied to unravel proteome patterns [corrected] of 2-D electrophoresis gel images

A novel approach for revealing patterns of proteome variation among series of 2-DE gel images is presented. The approach utilises image alignment to ensure that each pixel represents the same information across all gels. Gel images are normalised, and background corrected, followed by unfolding of the images to 1-D pixel vectors and analysing pixel vectors by multivariate data modelling. Information resulting from the data analysis is refolded back to the image domain for visualisation and interpretation. The method is rapid and suitable for automatic routines applied after the gel alignment. The approach is compared with spot volume analysis to illustrate how this approach can solve persistent problems like mismatch of protein spots, erroneous missing values and failure to detect variation in overlapping proteins. The method may also detect variation in the border area of saturated proteins. The approach is given the name pixel-based analysis of multiple images for the identification of changes (PMC). The method can be used for multiple images in general. Effects of pretreatment of the images are discussed.