Rainbow trout (Oncorhynchus mykiss) resistance to columnaris disease is heritable and favorably correlated with bacterial cold water disease resistance

Potential of low-density genotype imputation for cost-efficient genomic selection for resistance to Flavobacterium columnare in rainbow trout (Oncorhynchus mykiss)

BACKGROUND

Flavobacterium columnare is the pathogen agent of columnaris disease, a major emerging disease that affects rainbow trout aquaculture. Selective breeding using genomic selection has potential to achieve cumulative improvement of the host resistance. However, genomic selection is expensive partly because of the cost of genotyping large numbers of animals using high-density single nucleotide polymorphism (SNP) arrays. The objective of this study was to assess the efficiency of genomic selection for resistance to F. columnare using in silico low-density (LD) panels combined with imputation. After a natural outbreak of columnaris disease, 2874 challenged fish and 469 fish from the parental generation (n = 81 parents) were genotyped with 27,907 SNPs. The efficiency of genomic prediction using LD panels was assessed for 10 panels of different densities, which were created in silico using two sampling methods, random and equally spaced. All LD panels were also imputed to the full 28K HD panel using the parental generation as the reference population, and genomic predictions were re-evaluated. The potential of prioritizing SNPs that are associated with resistance to F. columnare was also tested for the six lower-density panels.

RESULTS

The accuracies of both imputation and genomic predictions were similar with random and equally-spaced sampling of SNPs. Using LD panels of at least 3000 SNPs or lower-density panels (as low as 300 SNPs) combined with imputation resulted in accuracies that were comparable to those of the 28K HD panel and were 11% higher than the pedigree-based predictions.

CONCLUSIONS

Compared to using the commercial HD panel, LD panels combined with imputation may provide a more affordable approach to genomic prediction of breeding values, which supports a more widespread adoption of genomic selection in aquaculture breeding programmes.

Resistance to Saprolegnia parasitica infection: A heritable trait in Atlantic salmon

A controlled Saprolegnia parasitica infection model was used to challenge 1158 fish representing 105 pedigreed Atlantic salmon families to evaluate the possibility of selecting for Saprolegnia resistance in a commercial breeding programme. Fish were infected in five study tanks and observed for 40 days post-infection for lesion score and survival. Survival analysis of the top 10 resistant and bottom 10 susceptible families indicated that the hazard of dying following Saprolegnia infection was 1509% higher in susceptible families. In all fish, a 10 g increase in weight correlated with a 7.8% increase in the hazard of dying while sex did not affect mortality. Resistance to Saprolegnia was estimated to have a heritability of 0.25, indicating that selection is possible. Genetic and phenotypic correlations indicated that the 11-point scoring system, developed in this study to quantify Saprolegnia infection severity, had a high negative correlation with survival as days to mortality at ≥-0.922(±0.005), suggesting that the scoring method could help assess lesion development in studies where mortality is not the primary biological endpoint.

What Can Genetics Do for the Control of Infectious Diseases in Aquaculture?

Simple Summary

Infectious diseases place an economic burden on aquaculture and a limitation to its growth. This state-of-the-art review describes the application of genetics and genomics as novel tools to control infectious disease in aquaculture.

Abstract

Infectious diseases place an economic burden on aquaculture and a limitation to its growth. An innovative approach to mitigate their impact on production is breeding for disease resistance: selection for domestication, family-based selection, marker-assisted selection, and more recently, genomic selection. Advances in genetics and genomics approaches to the control of infectious diseases are key to increasing aquaculture efficiency, profitability, and sustainability and to reducing its environmental footprint. Interaction and co-evolution between a host and pathogen can, however, turn breeding to boost infectious disease resistance into a potential driver of pathogenic change. Parallel molecular characterization of the pathogen and its virulence and antimicrobial resistance genes is therefore essential to understand pathogen evolution over time in response to host immunity, and to apply appropriate mitigation strategies.

Conserved QTL and chromosomal inversion affect resistance to columnaris disease in 2 rainbow trout (Oncorhyncus mykiss) populations

We present a comparative genetic analysis of the quantitative trait loci underlying resistance to warm water columnaris disease in 2 farmed rainbow trout (Oncorhynchus mykiss) populations. We provide evidence for the conservation of a major quantitative trait loci on Omy03, and the putative role played by a chromosomal rearrangement on Omy05. A total of 3,962 individuals from the 2 populations experienced a natural Flavobacterium columnare outbreak. Data for 25,823 genome-wide SNPs were generated for both cases (fatalities) and controls (survivors). F_ST and pairwise additive genetic relationships suggest that, despite being currently kept as separate broodstocks, the 2 populations are closely related. Association analyses identified a major quantitative trait loci on chromosome Omy03 and a second smaller quantitative trait loci on Omy05. Quantitative trait loci on Omy03 consistently explained 3–11% of genetic variation in both populations, whereas quantitative trait loci on Omy05 showed different degree of association across populations and sexes. The quantitative trait loci on Omy05 was found within a naturally occurring, 54.84 cM long inversion which is easy to tag due to a strong linkage disequilibrium between the 375 tagging SNPs. The ancestral haplotype on Omy05 was associated with decreased mortality. Genetic correlation between mortality in the 2 populations was estimated at 0.64, implying that the genetic basis of resistance is partly similar in the 2 populations. Our quantitative trait loci validation identifies markers that can be potentially used to complement breeding value evaluations to increase resistance against columnaris disease, and help to mitigate effects of climate change on aquaculture.

Growth of Flavobacterium columnare Genomovars in the Presence or Lack of Supplemental Cations

Flavobacterium columnare is a problematic pathogen for the aquaculture industry where isolates are classified by genomovars. Suspended growth in a low nutrient media, like tryptone yeast extract salts, is a common method used for laboratory study. The presence of calcium and magnesium is the factor contributing to growth, virulence, and biofilm formation for F. columnare. Exponential growth occurs within 24 h for F. columnare when grown in complete tryptone yeast extract salts medium at 30 ºC. Withholding CaCl₂ and MgSO₄ components from a complete TYES formulation reduced or completely inhibited growth of genomovar I isolates but not the growth of genomovar II, IIB, or III isolates. Only 3 of 20 genomovar I isolates, MS-FC-4, FC-CSF-53, and 023-08-3, could achieve O.D. 540 readings ≥ 0.3 but only after 48-h incubation in cation-restricted TYES. Independently adding CaCl₂ or MgSO₄ to tryptone and yeast extracts did not result in a genomovar-specific growth phenotype, but generally demonstrated increased clumping with individual isolates presenting abnormal growth. Clumping formed filamentous strings that migrated to the top of the culture tube when isolates were grown in TYE+CaCl₂. Several of the F. columnare isolates from all the genomovars exhibited delayed growth when a single cation source was provided. This study demonstrates phenotypic differences between and within genomovars of a single bacterial species when grown under different TYES media conditions.

Virulence variations of Flavobacterium columnare in rainbow trout (Oncorhynchus mykiss) eyed eggs and alevin

Flavobacterium columnare (Fc) is the causative agent for columnaris disease (CD) in several fish species and an emerging problem for rainbow trout aquaculture. We characterize the virulence phenotype of two Fc isolates, CSF-298-10 and MS-FC-4, against trout from two sources, NCCCWA and a production stock (PS), at the eyed egg and alevin life stages. Immersion challenges demonstrated that NCCCWA eyed eggs were susceptible to the Fc isolate MS-FC-4 (>97% mortality) but no mortality was observed against PS eyed eggs. The CSF-298-10 had little effect on any eyed eggs tested and was not highly virulent to any alevin till day six post-hatch, up to 38% for NCCCWA and ~80% PS alevin. The MS-FC-4 strain produced ≥80% mortality any day an immersion challenge occurred post-hatch. Significant difference in CFU counts was recorded between the Fc strains on 2 days post-hatch immersion challenges. Counts for the NCCCWA alevin were 4.4 × 10³  CFU/ml^-1 and 1.8 × 10⁶  CFU/ml^-1 for the CSF-298-10 strain and MS-FC-4 strain, respectively, and for the PS alevin CSF-298-10 measured 9.9 × 10¹  CFU/ml^-1 and 3.8 × 10⁵  CFU/ml^-1 for MS-FC-4. These two Fc isolates present stark differences in virulence phenotypes to both eyed eggs and alevin and present an interesting model system for virulence kinetics and potentially alternative pathogenic pathways.

Genomic selection for heterobothriosis resistance concurrent with body size in the tiger pufferfish, Takifugu rubripes

Parasite resistance traits in aquaculture species often have moderate heritability, indicating the potential for genetic improvements by selective breeding. However, parasite resistance is often synonymous with an undesirable negative correlation with body size. In this study, we first tested the feasibility of genomic selection (GS) on resistance to heterobothriosis, caused by the monogenean parasite Heterobothrium okamotoi, which leads to huge economic losses in aquaculture of the tiger pufferfish Takifugu rubripes. Then, using a simulation study, we tested the possibility of simultaneous improvement of parasite resistance, assessed by parasite counts on host fish (HC), and standard length (SL). Each trait showed moderate heritability (square-root transformed HC: h2 = 0.308 ± 0.123, S.E.; SL: h2 = 0.405 ± 0.131). The predictive abilities of genomic prediction among 12 models, including genomic Best Linear Unbiased Predictor (GBLUP), Bayesian regressions, and machine learning procedures, were also moderate for both transformed HC (0.248‒0.344) and SL (0.340‒0.481). These results confirmed the feasibility of GS for this trait. Although an undesirable genetic correlation was suggested between transformed HC and SL (rg = 0.228), the simulation study suggested the desired gains index can help achieve simultaneous genetic improvements in both traits.

Salmonid Antibacterial Immunity: An Aquaculture Perspective

The aquaculture industry is continuously threatened by infectious diseases, including those of bacterial origin. Regardless of the disease burden, aquaculture is already the main method for producing fish protein, having displaced capture fisheries. One attractive sector within this industry is the culture of salmonids, which are (a) uniquely under pressure due to overfishing and (b) the most valuable finfish per unit of weight. There are still knowledge gaps in the understanding of fish immunity, leading to vaccines that are not as effective as in terrestrial species, thus a common method to combat bacterial disease outbreaks is the use of antibiotics. Though effective, this method increases both the prevalence and risk of generating antibiotic-resistant bacteria. To facilitate vaccine design and/or alternative treatment efforts, a deeper understanding of the teleost immune system is essential. This review highlights the current state of teleost antibacterial immunity in the context of salmonid aquaculture. Additionally, the success of current techniques/methods used to combat bacterial diseases in salmonid aquaculture will be addressed. Filling the immunology knowledge gaps highlighted here will assist in reducing aquaculture losses in the future.

The Type IX Secretion System Is Required for Virulence of the Fish Pathogen Flavobacterium psychrophilum

Flavobacterium psychrophilum causes bacterial cold-water disease in wild and aquaculture-reared fish and is a major problem for salmonid aquaculture. The mechanisms responsible for cold-water disease are not known. It was recently demonstrated that the related fish pathogen, Flavobacterium columnare, requires a functional type IX protein secretion system (T9SS) to cause disease. T9SSs secrete cell surface adhesins, gliding motility proteins, peptidases, and other enzymes, any of which may be virulence factors. The F. psychrophilum genome has genes predicted to encode components of a T9SS. Here, we used a SacB-mediated gene deletion technique recently adapted for use in the Bacteroidetes to delete a core F. psychrophilum T9SS gene, gldN The ΔgldN mutant cells were deficient for secretion of many proteins in comparison to wild-type cells. Complementation of the mutant with wild-type gldN on a plasmid restored secretion. Compared to wild-type and complemented strains, the ΔgldN mutant was deficient in adhesion, gliding motility, and extracellular proteolytic and hemolytic activities. The ΔgldN mutant exhibited reduced virulence in rainbow trout and complementation restored virulence, suggesting that the T9SS plays an important role in the disease.IMPORTANCE Bacterial cold-water disease, caused by F. psychrophilum, is a major problem for salmonid aquaculture. Little is known regarding the virulence factors involved in this disease, and control measures are inadequate. A targeted gene deletion method was adapted to F. psychrophilum and used to demonstrate the importance of the T9SS in virulence. Proteins secreted by this system are likely virulence factors and targets for the development of control measures.

Recirculation versus flow-through rainbow trout laboratory Flavobacterium columnare challenge

Flavobacterium columnare immersion challenges are affected by water-related environmental parameters and thus are difficult to reproduce. Whereas these challenges are typically conducted using flow-through systems, use of a recirculating challenge system to control environmental parameters may improve reproducibility. We compared mortality, bacterial concentration, and environmental parameters between flow-through and recirculating immersion challenge systems under laboratory conditions using 20 rainbow trout families. Despite identical dose concentration (1:75 dilution), duration of challenge, lot of fish, and temperature, average mortality in the recirculating system (42%) was lower (p < 0.01) compared to the flow-through system (77%), and there was low correlation (r = 0.24) of family mortality. Mean days to death (3.25 vs. 2.99 d) and aquaria-to-aquaria variation (9.6 vs. 10.4%) in the recirculating and flow-through systems, respectively, did not differ (p ≥ 0.30). Despite 10-fold lower water replacement rate in the recirculating (0.4 exchanges h-1) compared to flow-through system (4 exchanges h-1), differences in bacterial concentration between the 2 systems were modest (≤0.6 orders of magnitude) and inconsistent throughout the 21 d challenge. Compared to the flow-through system, dissolved oxygen during the 1 h exposure and pH were greater (p ≤ 0.02), and calcium and hardness were lower (p ≤ 0.03), in the recirculating system. Although this study was not designed to test effects of specific environmental parameters on mortality, it demonstrates that the cumulative effects of these parameters result in poor reproducibility. A recirculating immersion challenge model may be warranted to empirically identify and control environmental parameters affecting mortality and thus may serve as a more repeatable laboratory challenge model.

Resistance to Streptococcus iniae and its genetic associations with traits of economic importance in Asian seabass (Lates calcarifer)

Streptococcus iniae is one of the most serious aquatic pathogens, causing significant economic losses in marine and freshwater species, including Asian seabass (Lates calcarifer). Controlling this gram-positive bacterial pathogen has been an issue in aquaculture systems, due to the combined effects of aquaculture intensification and climatic impacts. To date, there have not been any genetic parameter estimates for S. iniae resistance in Asian seabass. The main aim of this study was to examine genetic variation in S. iniae resistance and its genetic correlations with growth and cannibalism in Asian seabass families produced from a breeding programme for high growth in 2016 and 2017. The study included a total of 5,835 individual fish that were offspring of 41 sires and 60 dams (31 half-sib and 34 full-sib families). The experimental fish were challenged by intraperitoneal injection with a volume containing 10⁵  CFU (colony-forming unit)/fish. Resistance to S. iniae was measured as survival rate at 6 hr, 3, 5, 7, 10 and 15 days post-challenge test. There were significant variations in S. iniae resistance among families at different observation periods (ranging from 24.4% to 80%). Restricted maximum-likelihood method and mixed model analysis were applied to estimate heritability for S. iniae resistance. The heritability for S. iniae resistance ranged from 7% to 18% across different statistical models used. The common full-sib effects accounted for 0.1%-2% of the total variation in resistance to S. iniae. Genetic correlations of the S. iniae resistance at 6 hr and 3 days with later post-challenge test periods were low to moderate. However, these estimates for S. iniae resistance between successive measurement times (5, 7, 10 and 15 days) were high and close to 1. The genetic correlations of resistance with body weights at 180, 270 and 360 days post-hatch were not significant as well with cannibalism. It is concluded that there is substantial additive genetic variation in resistance to S. iniae, suggesting there is potential for genetic improvement of Asian seabass for resistance to S. iniae through selective breeding.

Breeding for improved resistance to Edwardsiella ictaluri in striped catfish (Pangasianodon hypophthalmus): Quantitative genetic parameters

Bacillary necrosis of Pangasius (BNP), caused by Edwardsiella ictaluri, is one of the most devastating diseases in striped catfish farming. To date, quantitative genetic inheritance of BNP resistance is not known in striped catfish Pangasianodon hypophthalmus. The main aim of this study was to estimate genetic parameters for BNP resistance in a breeding population of striped catfish undergoing four generations of selection for high growth. Specifically, the study examined whether BNP resistance is heritable to enable family selection and whether genetic improvement for enhanced BNP resistance may have detrimental effects on growth and survival rate. To test these hypotheses, 720 full- and half-sib families were challenged with E. ictaluri pathogen using injection and cohabitation methods over four years, from 2010 to 2012 and 2015. In total, the data included 398,234 animals in the pedigree, from which 18,849 animals had disease challenge test records and 39,103 siblings had growth performance. Both univariate and bivariate sire-dam linear and threshold mixed models were used to estimate (co)variance components for BNP resistance, survivals and growth traits. The estimates of heritability for the BNP resistance recorded as death or survival were low regardless of models used (0.10-0.16), whereas survival time (days post-challenge test) showed moderate heritability (0.35). The survival rate during hapa rearing had medium heritability (0.33-0.52). The genetic correlations of BNP resistance with body weight and survival were all positive (0.03-0.53), suggesting that selection of increased BNP resistance may have positive impacts on growth and survival traits, and these traits could be easily improved simultaneously in the selective breeding programme for striped catfish.

Whole-genome mapping of quantitative trait loci and accuracy of genomic predictions for resistance to columnaris disease in two rainbow trout breeding populations

Background

Columnaris disease (CD) is an emerging problem for the rainbow trout aquaculture industry in the US. The objectives of this study were to: (1) identify common genomic regions that explain a large proportion of the additive genetic variance for resistance to CD in two rainbow trout (Oncorhynchus mykiss) populations; and (2) estimate the gains in prediction accuracy when genomic information is used to evaluate the genetic potential of survival to columnaris infection in each population.

Methods

Two aquaculture populations were investigated: the National Center for Cool and Cold Water Aquaculture (NCCCWA) odd-year line and the Troutlodge, Inc., May odd-year (TLUM) nucleus breeding population. Fish that survived to 21 days post-immersion challenge were recorded as resistant. Single nucleotide polymorphism (SNP) genotypes were available for 1185 and 1137 fish from NCCCWA and TLUM, respectively. SNP effects and variances were estimated using the weighted single-step genomic best linear unbiased prediction (BLUP) for genome-wide association. Genomic regions that explained more than 1% of the additive genetic variance were considered to be associated with resistance to CD. Predictive ability was calculated in a fivefold cross-validation scheme and using a linear regression method.

Results

Validation on adjusted phenotypes provided a prediction accuracy close to zero, due to the binary nature of the trait. Using breeding values computed from the complete data as benchmark improved prediction accuracy of genomic models by about 40% compared to the pedigree-based BLUP. Fourteen windows located on six chromosomes were associated with resistance to CD in the NCCCWA population, of which two windows on chromosome Omy 17 jointly explained more than 10% of the additive genetic variance. Twenty-six windows located on 13 chromosomes were associated with resistance to CD in the TLUM population. Only four associated genomic regions overlapped with quantitative trait loci (QTL) between both populations.

Conclusions

Our results suggest that genome-wide selection for resistance to CD in rainbow trout has greater potential than selection for a few target genomic regions that were found to be associated to resistance to CD due to the polygenic architecture of this trait, and because the QTL associated with resistance to CD are not sufficiently informative for selection decisions across populations.

Electronic supplementary material

The online version of this article (10.1186/s12711-019-0484-4) contains supplementary material, which is available to authorized users.

Variance and covariance estimates for resistance to bacterial cold water disease and columnaris disease in two rainbow trout breeding populations1

Family-based selective breeding can be an effective strategy for controlling diseases in aquaculture. This study aimed to estimate (co)variance components for resistance to bacterial cold water disease (BCWD) and columnaris disease (CD) in two unrelated rainbow trout nucleus breeding populations: the USDA, ARS, National Center for Cool and Cold Water Aquaculture odd-year line (ARS-Fp-R), which has been subjected to five generations of selection for improved resistance to BCWD, and the Troutlodge, Inc., May-spawning odd-year line (TLUM), which has been selected for improved growth performance but not for disease resistance. A total of 46,805 and 27,821 pedigree records were available from both populations, respectively. Between 44 and 138 families per generation and population were evaluated under controlled BCWD and CD challenges, providing 32,311 and 17,861 phenotypic records for BCWD resistance, and 13,603 and 9,413 for CD resistance, in the ARS-Fp-R and TLUM populations, respectively. A two-trait animal threshold model assuming an underlying normal distribution for the binary survival phenotypes was used to estimate (co)variance components separately for each population. Resistance to BCWD (h2 = 0.27 ± 0.04 and 0.43 ± 0.08) and CD (h2 = 0.23 ± 0.07 and 0.34 ± 0.09) was moderately heritable in the ARS-Fp-R and TLUM populations, respectively. The genetic correlation between the resistance to BCWD and CD was favorably positive in the ARS-Fp-R (0.40 ± 0.17) and TLUM (0.39 ± 0.18) populations. These findings suggest that both disease resistance traits can be improved simultaneously even if genetic selection pressure is applied to only one of the two traits.

An 8-Year Breeding Program for Asian Seabass Lates calcarifer: Genetic Evaluation, Experiences, and Challenges

Selective breeding for marine finfish is challenging due to difficulties in reproduction, larval rearing, and on-growth in captive environments. The farming of Asian seabass (Lates calcarifer) has all these problems and our knowledge of the quantitative genetic information (heritability and correlations) of traits necessary for commercial exploitation is poor. The present study was conducted to address this knowledge gap and to provide information that can be applied to sea bass and other aquaculture species. We carried out a comprehensive genetic evaluation for three traits (body weight, total length, and survival) collected from a breeding population for Asian seabass over an eight-year period from 2010 to 2017. Statistical analysis was carried out on 4,567 adult fish at 105, 180, 270, 360, 450, and 570 days post-hatch (dph). The heritabilities (h²) estimated for body weight and length using linear mixed model were moderate to high (0.12 to 0.78 and 0.41 to 0.85, respectively) and they differed between the measurement periods. Survival during grow-out phase was analyzed using threshold logistic and probit models. The heritability estimates for survival rate on the underlying liability scale ( hL2 ) varied from 0.05 to 0.21. When the observed heritability obtained from the linear mixed model was back-transformed to the liability scale, they were similar but not significant. In addition, we examined effects of genotype by environment (G × E) interaction on body traits. The genetic correlation for body weight between tank and sea cage cultures were high (0.91-0.94) in the first and second rearing periods (180 and 270 dph) but the correlation was decreased to 0.59 ± 0.33 at 360 dph. This suggests that the genotype by environment interaction is important for body traits in this population. Furthermore, the genetic correlations of body weights between different measurement periods were moderate but different from one. This suggests that body weights measured at different time points may be different traits and selection for improved early weight may not capture all genetic expressions in subsequent rearing periods in Asian seabass. Selection of the nucleus in sea cages may produce genotypes that do not perform equally well in tanks, although this deserves further studies to determine a suitable selection environment and optimize the breeding program. This paper discusses challenges encountered during implementation of the selection program for L. calcarifer.

Identification of Four Distinct Phylogenetic Groups in Flavobacterium columnare With Fish Host Associations

Columnaris disease, caused by the Gram-negative bacterium Flavobacterium columnare, is one of the most prevalent fish diseases worldwide. An exceptionally high level of genetic diversity among isolates of F. columnare has long been recognized, whereby six established genomovars have been described to date. However, little has been done to quantify or characterize this diversity further in a systematic fashion. The objective of this research was to perform phylogenetic analyses of 16S rRNA and housekeeping gene sequences to decipher the genetic diversity of F. columnare. Fifty isolates and/or genomes of F. columnare, originating from diverse years, geographic locations, fish hosts, and representative of the six genomovars were analyzed in this study. A multilocus phylogenetic analysis (MLPA) of the 16S rRNA and six housekeeping genes supported four distinct F. columnare genetic groups. There were associations between genomovar and genetic group, but these relationships were imperfect indicating that genomovar assignment does not accurately reflect F. columnare genetic diversity. To expand the dataset, an additional 90 16S rRNA gene sequences were retrieved from GenBank and a phylogenetic analysis of this larger dataset also supported the establishment of four genetic groups. Examination of isolate historical data indicated biological relevance to the identified genetic diversity, with some genetic groups isolated preferentially from specific fish species or families. It is proposed that F. columnare isolates be assigned to the four genetic groups defined in this study rather than genomovar in order to facilitate a standard nomenclature across the scientific community. An increased understanding of which genetic groups are most prevalent in different regions and/or aquaculture industries may allow for the development of improved targeted control and treatment measures for columnaris disease.

Genetic Characterization of the Fish Piaractus brachypomus by Microsatellites Derived from Transcriptome Sequencing

The pirapitinga, Piaractus brachypomus (Characiformes, Serrasalmidae), is a fish from the Amazon basin and is considered to be one of the main native species used in aquaculture production in South America. The objectives of this study were: (1) to perform liver transcriptome sequencing of pirapitinga through NGS and then validate a set of microsatellite markers for this species; and (2) to use polymorphic microsatellites for analysis of genetic variability in farmed stocks. The transcriptome sequencing was carried out through the Roche/454 technology, which resulted in 3,696 non-redundant contigs. Of this total, 2,568 contigs had similarity in the non-redundant (nr) protein database (Genbank) and 2,075 sequences were characterized in the categories of Gene Ontology (GO). After the validation process of 30 microsatellite loci, eight markers showed polymorphism. The analysis of these polymorphic markers in farmed stocks revealed that fish farms from North Brazil had a higher genetic diversity than fish farms from Southeast Brazil. AMOVA demonstrated that the highest proportion of variation was presented within the populations. However, when comparing different groups (1: Wild; 2: North fish farms; 3: Southeast fish farms), a considerable variation between the groups was observed. The F_ST values showed the occurrence of genetic structure among the broodstocks from different regions of Brazil. The transcriptome sequencing in pirapitinga provided important genetic resources for biological studies in this non-model species, and microsatellite data can be used as the framework for the genetic management of breeding stocks in Brazil, which might provide a basis for a genetic pre-breeding programme.

Microsatellites Associated with Growth Performance and Analysis of Resistance to Aeromonas hydrophila in Tambaqui Colossoma macropomum

Tambaqui, Colossoma macropomum, is the main native fish species produced in Brazil, and is an important species for genetic improvement in aquaculture. In addition, breeding studies on this species can be optimized with the use of molecular markers associated with productive phenotypes. The objective of the present study was to test the performance of growth traits and resistance to the bacteria, Aeromonas hydrophila, in association with microsatellite markers in C. macropomum. In this study, three full-sib families were subjected to bacterial challenge and morphometric growth assessments. Tambaqui families subjected to the bacterial challenge differed significantly in death time and mortality rate. There was, however, no association between resistance to bacteria and microsatellite markers. In relation to growth traits, we observed a marker/phenotype association in two microsatellites. The marker in the 6b isoform x5 gene (TNCRC6b) was associated with length, whereas an anonymous marker was associated with height. The present study highlighted the evaluation of molecular markers associated with growth traits, and can serve as the basis for future marker-assisted selection (MAS) of tambaqui.

Draft Genome Sequence of the Fish Pathogen Flavobacterium columnare Strain CSF-298-10

We announce here the draft genome assembly of Flavobacterium columnare CSF-298-10, a strain isolated from an outbreak of columnaris disease at a commercial trout farm in Hagerman Valley, Idaho, USA. The complete genome consists of 13 contigs totaling 3,284,579 bp, with an average G+C content of 31.5% and 2,933 predicted coding genes.

Virulence of Flavobacterium columnare genomovars in rainbow trout Oncorhynchus mykiss

Flavobacterium columnare is the causative agent of columnaris disease and is responsible for significant economic losses in aquaculture. F. columnare is a Gram-negative bacterium, and 5 genetic types or genomovars have been described based on restriction fragment length polymorphism of the 16S rRNA gene. Previous research has suggested that genomovar II isolates are more virulent than genomovar I isolates to multiple species of fish, including rainbow trout Oncorhynchus mykiss. In addition, improved genotyping methods have shown that some isolates previously classified as genomovar I, and used in challenge experiments, were in fact genomovar III. Our objective was to confirm previous results with respect to genomovar II virulence, and to determine the susceptibility of rainbow trout to other genomovars. The virulence of 8 genomovar I, 4 genomovar II, 3 genomovar II-B, and 5 genomovar III isolates originating from various sources was determined through 3 independent challenges in rainbow trout using an immersion challenge model. Mean cumulative percent mortality (CPM) of ~49% for genomovar I isolates, ~1% for genomovar II, ~5% for the II-B isolates, and ~7% for the III isolates was observed. The inability of genomovar II isolates to produce mortalities in rainbow trout was unanticipated based on previous studies, but may be due to a number of factors including rainbow trout source and water chemistry. The source of fish and/or the presence of sub-optimal environment may influence the susceptibility of rainbow trout to different F. columnare genomovars.