Host specificity and virulence of Flavobacterium psychrophilum: a comparative study in ayu (Plecoglossus altivelis) and rainbow trout (Oncorhynchus mykiss) hosts.

Flavobacterium psychrophilum, the causative agent of bacterial cold-water disease, is a devastating, worldwide distributed, fish pathogen causing significant economic loss in inland fish farms. Previous epidemiological studies showed that prevalent clonal complexes (CC) differ in fish species affected with disease such as rainbow trout, coho salmon and ayu, indicating significant associations between particular F. psychrophilum genotypes and host species. Yet, whether the population structure is driven by the trade of fish and eggs or by host-specific pathogenicity is uncertain. Notably, all F. psychrophilum isolates retrieved from ayu belong to Type-3 O antigen (O-Ag) whereas only very few strains retrieved from other fish species possess this O-Ag, suggesting a role in outbreaks affecting ayu. Thus, we investigated the links between genotype and pathogenicity by conducting comparative bath infection challenges in two fish hosts, ayu and rainbow trout, for a collection of isolates representing different MLST genotypes and O-Ag. Highly virulent strains in one host species exhibited low to no virulence in the other. F. psychrophilum strains associated with ayu and possessing Type-3 O-Ag demonstrated significant variability in pathogenicity in ayu, ranging from avirulent to highly virulent. Strikingly, F. psychrophilum strains retrieved from rainbow trout and possessing the Type-3 O-Ag were virulent for rainbow trout but not for ayu, indicating that Type-3 O-Ag alone is not sufficient for pathogenicity in ayu, nor does it prevent pathogenicity in rainbow trout. This study revealed that the association between a particular CC and host species partly depends on the pathogen's adaptation to specific host species.

Strengths and limitations of reduced representation bisulfite sequencing (RRBS) in the perspective of DNA methylation analysis in fish: a case-study on rainbow trout spermatozoa.

DNA methylation in CpG dinucleotides is an important epigenetic mark in fish spermatozoa since it has been shown that some sperm methylome features are transmitted to the offspring. Reduced representation bisulfite sequencing (RRBS) is one genome-scale methods developed to assess DNA methylation at CpG sites. It allows the sequencing of a reduced fraction of the genome expected to be enriched in CpGs. The aim of this study is to characterize the extent of the CpG sites that can be identified in the RRBS-reduced sequenced fraction of rainbow trout spermatozoa, in order to evaluate the potential of RRBS for sperm DNA methylation studies. We observed that RRBS did provide a reduced amount of genomic data, the sum of the CpGs analyzed on 12 males spanning 9% of the total genomic CpGs. CpGs were only slightly enriched in the RRBS data (×1.7 times the sequenced nucleotides), the possible causes being linked to trout genome structure and sequenced fragments size. All genomic functional features were represented in our CpG dataset, with a noticeable enrichment in exons but, strikingly, not in promoters. The number of CpGs shared between biological replicates was low, but this proportion reached workable values from six biological replicates (46% of the analyzed cytosines) on. The choices that are to be made regarding fragment size selection and the options during bioinformatic data processing are discussed. In all, RRBS is a relevant first-approach method to scan the CpG DNA methylation status of spermatozoa along rainbow trout genome, although in a very reduced pattern among biological replicates.

Genetic architecture of acute hyperthermia resistance in juvenile rainbow trout (Oncorhynchus mykiss) and genetic correlations with production traits.

BACKGROUND: Selective breeding is a promising solution to reduce the vulnerability of fish farms to heat waves, which are predicted to increase in intensity and frequency. However, limited information about the genetic architecture of acute hyperthermia resistance in fish is available. Two batches of sibs from a rainbow trout commercial line were produced: the first (N = 1382) was phenotyped for acute hyperthermia resistance at nine months of age and the second (N = 1506) was phenotyped for main production traits (growth, body length, muscle fat content and carcass yield) at 20 months of age. Fish were genotyped on a 57 K single nucleotide polymorphism (SNP) array and their genotypes were imputed to high-density based on the parent's genotypes from a 665 K SNP array. RESULTS: The heritability estimate of resistance to acute hyperthermia was 0.29 ± 0.05, confirming the potential of selective breeding for this trait. Since genetic correlations of acute hyperthermia resistance with the main production traits near harvest age were all close to zero, selecting for acute hyperthermia resistance should not impact the main production traits, and vice-versa. A genome-wide association study revealed that resistance to acute hyperthermia is a highly polygenic trait, with six quantitative trait loci (QTL) detected, but explaining less than 5% of the genetic variance. Two of these QTL, including the most significant one, may explain differences in acute hyperthermia resistance across INRAE isogenic lines of rainbow trout. Differences in mean acute hyperthermia resistance phenotypes between homozygotes at the most significant SNP was 69% of the phenotypic standard deviation, showing promising potential for marker-assisted selection. We identified 89 candidate genes within the QTL regions, among which the most convincing functional candidates are dnajc7, hsp70b, nkiras2, cdk12, phb, fkbp10, ddx5, cygb1, enpp7, pdhx and acly. CONCLUSIONS: This study provides valuable insight into the genetic architecture of acute hyperthermia resistance in juvenile rainbow trout. We show that the selection potential for this trait is substantial and selection for this trait should not be too detrimental to improvement of other traits of interest. Identified functional candidate genes provide new knowledge on the physiological mechanisms involved in acute hyperthermia resistance, such as protein chaperoning, oxidative stress response, homeostasis maintenance and cell survival.

Convergence of multiple markers and analysis methods defines the genetic distinctiveness of cryptic pitvipers.

Using multiple markers and multiple analytical approaches is critical for establishing species boundaries reliably, especially so in the case of cryptic species. Despite development of new and powerful analytical methods, most studies continue to adopt a few, with the choice often being subjective. One such example is routine analysis of Amplified Fragment Length Polymorphism (AFLP) data using population genetic models despite disparity between method assumptions and data properties. The application of newly developed methods for analyzing this dominant marker may not be entirely clear in the context of species delimitation. In this study, we use AFLPs and mtDNA to investigate cryptic speciation in the Trimeresurus macrops complex that belongs to a taxonomically difficult lineage of Asian pitvipers. We analyze AFLPs using population genetic, phylogenetic, multivariate statistical, and Bayes Factor Delimitation methods. A gene tree from three mtDNA markers provided additional evidence. Our results show that the inferences about species boundaries that can be derived from population genetic analysis of AFLPs have certain limitations. In contrast, four multivariate statistical analyses produced clear clusters that are consistent with each other, as well as with Bayes Factor Delimitation results, and with mtDNA and total evidence phylogenies. Furthermore, our results concur with allopatric distributions and patterns of variation in individual morphological characters previously identified in the three proposed species: T. macrops sensu stricto, T. cardamomensis, and T. rubeus. Our study provides evidence for reproductive isolation and genetic distinctiveness that define these taxa as full species. In addition, we re-emphasize the importance of examining congruence of results from multiple methods of AFLP analysis for inferring species diversity.

Development of a High-Density 665 K SNP Array for Rainbow Trout Genome-Wide Genotyping.

Single nucleotide polymorphism (SNP) arrays, also named « SNP chips ¼, enable very large numbers of individuals to be genotyped at a targeted set of thousands of genome-wide identified markers. We used preexisting variant datasets from USDA, a French commercial line and 30X-coverage whole genome sequencing of INRAE isogenic lines to develop an Affymetrix 665 K SNP array (HD chip) for rainbow trout. In total, we identified 32,372,492 SNPs that were polymorphic in the USDA or INRAE databases. A subset of identified SNPs were selected for inclusion on the chip, prioritizing SNPs whose flanking sequence uniquely aligned to the Swanson reference genome, with homogenous repartition over the genome and the highest Minimum Allele Frequency in both USDA and French databases. Of the 664,531 SNPs which passed the Affymetrix quality filters and were manufactured on the HD chip, 65.3% and 60.9% passed filtering metrics and were polymorphic in two other distinct French commercial populations in which, respectively, 288 and 175 sampled fish were genotyped. Only 576,118 SNPs mapped uniquely on both Swanson and Arlee reference genomes, and 12,071 SNPs did not map at all on the Arlee reference genome. Among those 576,118 SNPs, 38,948 SNPs were kept from the commercially available medium-density 57 K SNP chip. We demonstrate the utility of the HD chip by describing the high rates of linkage disequilibrium at 2-10 kb in the rainbow trout genome in comparison to the linkage disequilibrium observed at 50-100 kb which are usual distances between markers of the medium-density chip.

Sources of variation of DNA methylation in rainbow trout: combined effects of temperature and genetic background.

Phenotypic plasticity is a key component of the ability of organisms to respond to changing environmental conditions. In this study, we aimed to study the establishment of DNA methylation marks in response to an environmental stress in rainbow trout and to assess whether these marks depend on the genetic background. The environmental stress chosen here was temperature, a known induction factor of epigenetic marks in fish. To disentangle the role of epigenetic mechanisms such as DNA methylation in generating phenotypic variations, nine rainbow trout isogenic lines with no genetic variability within a line were used. For each line, half of the eggs were incubated at standard temperature (11°C) and the other half at high temperature (16°C), from eyed-stage to hatching. In order to gain a first insight into the establishment of DNA methylation marks in response to an early temperature regime (control 11°C vs. heated 16°C), we have studied the expression of 8 dnmt3 (DNA methyltransferase) genes, potentially involved in de novo methylation, and analysed global DNA methylation in the different rainbow trout isogenic lines using LUMA (LUminometric Methylation Assay). Finally, finer investigation of genome-wide methylation patterns was performed using EpiRADseq, a reduced-representation library approach based on the ddRADseq (Double Digest Restriction Associated DNA) protocol, for six rainbow trout isogenic lines. We have demonstrated that thermal history during embryonic development alters patterns of DNA methylation, but to a greater or lesser extent depending on the genetic background.

Combining Multiple Approaches and Models to Dissect the Genetic Architecture of Resistance to Infections in Fish.

Infectious diseases represent a major threat for the sustainable development of fish farming. Efficient vaccines are not available against all diseases, and growing antibiotics resistance limits the use of antimicrobial drugs in aquaculture. It is therefore important to understand the basis of fish natural resistance to infections to help genetic selection and to develop new approaches against infectious diseases. However, the identification of the main mechanisms determining the resistance or susceptibility of a host to a pathogenic microbe is challenging, integrating the complexity of the variation of host genetics, the variability of pathogens, and their capacity of fast evolution and adaptation. Multiple approaches have been used for this purpose: (i) genetic approaches, QTL (quantitative trait loci) mapping or GWAS (genome-wide association study) analysis, to dissect the genetic architecture of disease resistance, and (ii) transcriptomics and functional assays to link the genetic constitution of a fish to the molecular mechanisms involved in its interactions with pathogens. To date, many studies in a wide range of fish species have investigated the genetic determinism of resistance to many diseases using QTL mapping or GWAS analyses. A few of these studies pointed mainly toward adaptive mechanisms of resistance/susceptibility to infections; others pointed toward innate or intrinsic mechanisms. However, in the majority of studies, underlying mechanisms remain unknown. By comparing gene expression profiles between resistant and susceptible genetic backgrounds, transcriptomics studies have contributed to build a framework of gene pathways determining fish responsiveness to a number of pathogens. Adding functional assays to expression and genetic approaches has led to a better understanding of resistance mechanisms in some cases. The development of knock-out approaches will complement these analyses and help to validate putative candidate genes critical for resistance to infections. In this review, we highlight fish isogenic lines as a unique biological material to unravel the complexity of host response to different pathogens. In the future, combining multiple approaches will lead to a better understanding of the dynamics of interaction between the pathogen and the host immune response, and contribute to the identification of potential targets of selection for improved resistance.

Standardized IMGT® Nomenclature of Salmonidae IGH Genes, the Paradigm of Atlantic Salmon and Rainbow Trout: From Genomics to Repertoires.

In teleost fish as in mammals, humoral adaptive immunity is based on B lymphocytes expressing highly diverse immunoglobulins (IG). During B cell differentiation, IG loci are subjected to genomic rearrangements of V, D, and J genes, producing a unique antigen receptor expressed on the surface of each lymphocyte. During the course of an immune response to infections or immunizations, B cell clones specific of epitopes from the immunogen are expanded and activated, leading to production of specific antibodies. Among teleost fish, salmonids comprise key species for aquaculture. Rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar) are especially important from a commercial point of view and have emerged as critical models for fish immunology. The growing interest to capture accurate and comprehensive antibody responses against common pathogens and vaccines has resulted in recent efforts to sequence the IG repertoire in these species. In this context, a unified and standardized nomenclature of salmonid IG heavy chain (IGH) genes is urgently required, to improve accuracy of annotation of adaptive immune receptor repertoire dataset generated by high-throughput sequencing (AIRRseq) and facilitate comparisons between studies and species. Interestingly, the assembly of salmonids IGH genomic sequences is challenging due to the presence of two large size duplicated IGH loci and high numbers of IG genes and pseudogenes. We used data available for Atlantic salmon to establish an IMGT standardized nomenclature of IGH genes in this species and then applied the IMGT rules to the rainbow trout IGH loci to set up a nomenclature, which takes into account the specificities of Salmonid loci. This unique, consistent nomenclature for Salmonid IGH genes was then used to construct IMGT sequence reference directories allowing accurate annotation of AIRRseq data. The complex issues raised by the genetic diversity of salmon and trout strains are discussed in the context of IG repertoire annotation.

Divergent national-scale trends of microbial and animal biodiversity revealed across diverse temperate soil ecosystems.

Soil biota accounts for ~25% of global biodiversity and is vital to nutrient cycling and primary production. There is growing momentum to study total belowground biodiversity across large ecological scales to understand how habitat and soil properties shape belowground communities. Microbial and animal components of belowground communities follow divergent responses to soil properties and land use intensification; however, it is unclear whether this extends across heterogeneous ecosystems. Here, a national-scale metabarcoding analysis of 436 locations across 7 different temperate ecosystems shows that belowground animal and microbial (bacteria, archaea, fungi, and protists) richness follow divergent trends, whereas ß-diversity does not. Animal richness is governed by intensive land use and unaffected by soil properties, while microbial richness was driven by environmental properties across land uses. Our findings demonstrate that established divergent patterns of belowground microbial and animal diversity are consistent across heterogeneous land uses and are detectable using a standardised metabarcoding approach.

Sample size effects on the assessment of eukaryotic diversity and community structure in aquatic sediments using high-throughput sequencing.

Understanding how biodiversity changes in time and space is vital to assess the effects of environmental change on benthic ecosystems. Due to the limitations of morphological methods, there has been a rapid expansion in the application of high-throughput sequencing methods to study benthic eukaryotic communities. However, the effect of sample size and small-scale spatial variation on the assessment of benthic eukaryotic diversity is still not well understood. Here, we investigate the effect of different sample volumes in the genetic assessment of benthic metazoan and non-metazoan eukaryotic community composition. Accordingly, DNA was extracted from five different cumulative sediment volumes comprising 100% of the top 2 cm of five benthic sampling cores, and used as template for Ilumina MiSeq sequencing of 18 S rRNA amplicons. Sample volumes strongly impacted diversity metrics for both metazoans and non-metazoan eukaryotes. Beta-diversity of treatments using smaller sample volumes was significantly different from the beta-diversity of the 100% sampled area. Overall our findings indicate that sample volumes of 0.2 g (1% of the sampled area) are insufficient to account for spatial heterogeneity at small spatial scales, and that relatively large percentages of sediment core samples are needed for obtaining robust diversity measurement of both metazoan and non-metazoan eukaryotes.

Annual time-series analysis of aqueous eDNA reveals ecologically relevant dynamics of lake ecosystem biodiversity.

The use of environmental DNA (eDNA) in biodiversity assessments offers a step-change in sensitivity, throughput and simultaneous measures of ecosystem diversity and function. There remains, however, a need to examine eDNA persistence in the wild through simultaneous temporal measures of eDNA and biota. Here, we use metabarcoding of two markers of different lengths, derived from an annual time series of aqueous lake eDNA to examine temporal shifts in ecosystem biodiversity and in an ecologically important group of macroinvertebrates (Diptera: Chironomidae). The analyses allow different levels of detection and validation of taxon richness and community composition (ß-diversity) through time, with shorter eDNA fragments dominating the eDNA community. Comparisons between eDNA, community DNA, taxonomy and UK species abundance data further show significant relationships between diversity estimates derived across the disparate methodologies. Our results reveal the temporal dynamics of eDNA and validate the utility of eDNA metabarcoding for tracking seasonal diversity at the ecosystem scale.

Genetic variability of environmental sensitivity revealed by phenotypic variation in body weight and (its) correlations to physiological and behavioral traits.

Adaptive phenotypic plasticity is a key component of the ability of organisms to cope with changing environmental conditions. Fish have been shown to exhibit a substantial level of phenotypic plasticity in response to abiotic and biotic factors. In the present study, we investigate the link between environmental sensitivity assessed globally (revealed by phenotypic variation in body weight) and more targeted physiological and behavioral indicators that are generally used to assess the sensitivity of a fish to environmental stressors. We took advantage of original biological material, the rainbow trout isogenic lines, which allowed the disentangling of the genetic and environmental parts of the phenotypic variance. Ten lines were characterized for the changes of body weight variability (weight measurements taken every month during 18 months), the plasma cortisol response to confinement stress (3 challenges) and a set of selected behavioral indicators. This study unambiguously demonstrated the existence of genetic determinism of environmental sensitivity, with some lines being particularly sensitive to environmental fluctuations and others rather insensitive. Correlations between coefficient of variation (CV) for body weight and behavioral and physiological traits were observed. This confirmed that CV for body weight could be used as an indicator of environmental sensitivity. As the relationship between indicators (CV weight, risk-taking, exploration and cortisol) was shown to be likely depending on the nature and intensity of the stressor, the joint use of several indicators should help to investigate the biological complexity of environmental sensitivity.

Metabarcoding Marine Sediments: Preparation of Amplicon Libraries.

The accurate assessment of community composition and ultimately species identification is of utmost importance in any ecological and evolutionary study. Advances in sequencing technologies have allowed the unraveling of levels of biodiversity never imagined before when applied to large-scale environmental DNA studies (also termed metabarcoding/metagenetics/metasystematics/environmental barcoding). Here, we describe a detailed protocol to assess eukaryotic biodiversity in marine sediments, identifying key steps that should not be neglected when preparing Next-Generation Sequencing (NGS) amplicon libraries: DNA extraction, multiple PCR amplification of DNA barcode markers with index/ tag-primers, and final Illumina MiSeq sequencing library preparation.

Environmental metabarcoding reveals heterogeneous drivers of microbial eukaryote diversity in contrasting estuarine ecosystems.

Assessing how natural environmental drivers affect biodiversity underpins our understanding of the relationships between complex biotic and ecological factors in natural ecosystems. Of all ecosystems, anthropogenically important estuaries represent a 'melting pot' of environmental stressors, typified by extreme salinity variations and associated biological complexity. Although existing models attempt to predict macroorganismal diversity over estuarine salinity gradients, attempts to model microbial biodiversity are limited for eukaryotes. Although diatoms commonly feature as bioindicator species, additional microbial eukaryotes represent a huge resource for assessing ecosystem health. Of these, meiofaunal communities may represent the optimal compromise between functional diversity that can be assessed using morphology and phenotype-environment interactions as compared with smaller life fractions. Here, using 454 Roche sequencing of the 18S nSSU barcode we investigate which of the local natural drivers are most strongly associated with microbial metazoan and sampled protist diversity across the full salinity gradient of the estuarine ecosystem. In order to investigate potential variation at the ecosystem scale, we compare two geographically proximate estuaries (Thames and Mersey, UK) with contrasting histories of anthropogenic stress. The data show that although community turnover is likely to be predictable, taxa are likely to respond to different environmental drivers and, in particular, hydrodynamics, salinity range and granulometry, according to varied life-history characteristics. At the ecosystem level, communities exhibited patterns of estuary-specific similarity within different salinity range habitats, highlighting the environmental sequencing biomonitoring potential of meiofauna, dispersal effects or both.

Characterization of 27 microsatellite loci in the European flat oyster Ostrea edulis.

The flat oyster Ostrea edulis is native to Europe and populations have been severely depleted by the parasite Bonamia ostreae since the 1980s. Additional genetic markers are required to improve population genetics study and linkage map development for selection for B. ostrea-resistance in this species. Here, we characterized 27 novel microsatellite loci for O. edulis. Number of alleles per locus ranged from 6 to 25 and observed heterozygosity between 0.375 and 1. Null alleles were suggested at a few loci but most loci were in Hardy-Weinberg agreement enabling their reliable use in further population and mapping genetics approaches.