Nonparallelism in MHCIIβ diversity accompanies nonparallelism in pathogen infection of lake whitefish (Coregonus clupeaformis) species pairs as revealed by next-generation sequencing

Template-specific optimization of NGS genotyping pipelines reveals allele-specific variation in MHC gene expression

Using high-throughput sequencing for precise genotyping of multi-locus gene families, such as the major histocompatibility complex (MHC), remains challenging, due to the complexity of the data and difficulties in distinguishing genuine from erroneous variants. Several dedicated genotyping pipelines for data from high-throughput sequencing, such as next-generation sequencing (NGS), have been developed to tackle the ensuing risk of artificially inflated diversity. Here, we thoroughly assess three such multi-locus genotyping pipelines for NGS data, the DOC method, AmpliSAS and ACACIA, using MHC class IIβ data sets of three-spined stickleback gDNA, cDNA and "artificial" plasmid samples with known allelic diversity. We show that genotyping of gDNA and plasmid samples at optimal pipeline parameters was highly accurate and reproducible across methods. However, for cDNA data, the gDNA-optimal parameter configuration yielded decreased overall genotyping precision and consistency between pipelines. Further adjustments of key clustering parameters were required tο account for higher error rates and larger variation in sequencing depth per allele, highlighting the importance of template-specific pipeline optimization for reliable genotyping of multi-locus gene families. Through accurate paired gDNA-cDNA typing and MHC-II haplotype inference, we show that MHC-II allele-specific expression levels correlate negatively with allele number across haplotypes. Lastly, sibship-assisted cDNA-typing of MHC-I revealed novel variants linked in haplotype blocks, and a higher-than-previously-reported individual MHC-I allelic diversity. In conclusion, we provide novel genotyping protocols for the three-spined stickleback MHC-I and -II genes, and evaluate the performance of popular NGS-genotyping pipelines. We also show that fine-tuned genotyping of paired gDNA-cDNA samples facilitates amplification bias-corrected MHC allele expression analysis.

Environmental correlates of adaptive diversification in postglacial freshwater fishes

Determining how environmental conditions contribute to divergence among populations and drive speciation is fundamental to resolving mechanisms and understanding outcomes in evolutionary biology. Postglacial freshwater fish species in the Northern Hemisphere are ideal biological systems to explore the effects of environment on diversification in morphology, ecology, and genetics (ecomorph divergences) within lakes. To date, various environmental factors have been implicated in the presence of multiple ecomorphs within particular lakes or regions. However, concerted evidence for generalizable patterns in environmental variables associated with speciation across geographical regions and across species and genera has been lacking. Here, we aimed to identify key biotic and abiotic factors associated with ecological divergence of postglacial freshwater fish species into multiple sympatric ecomorphs, focusing on species in the well-studied, widespread, and co-distributed genera Gasterosteus, Salvelinus, and Coregonus (stickleback, charr, and whitefish, respectively). We found that the presence of multiple sympatric ecomorphs tended to be associated with increasing lake surface area, maximum depth, and nutrient availability. In addition, predation, competition, and prey availability were suggested to play a role in divergence into multiple ecomorphs, but the effects of biotic factors require further study. Although we identified several environmental factors correlated with the presence of multiple ecomorphs, there were substantial data gaps across species and regions. An improved understanding of these systems may provide insight into both generalizable environmental factors involved in speciation in other systems, and potential ecological and evolutionary responses of species complexes when these variables are altered by environmental change.

Tenets in Microbial Endocrinology: A New Vista in Teleost Reproduction

Climate vulnerability and induced changes in physico-chemical properties of aquatic environment can bring impairment in metabolism, physiology and reproduction in teleost. Variation in environmental stimuli mainly acts on reproduction by interfering with steroidogenesis, gametogenesis and embryogenesis. The control on reproductive function in captivity is essential for the sustainability of aquaculture production. There are more than 3,000 teleost species across the globe having commercial importance; however, adequate quality and quantity of seed production have been the biggest bottleneck. Probiotics are widely used in aquaculture as a growth promoter, stress tolerance, pathogen inhibition, nutrient digestibility and metabolism, reproductive performance and gamete quality. As the gut microbiota exerts various effects on the intestinal milieu which influences distant organs and pathways, therefore it is considered to be a full-fledged endocrine organ. Researches on Gut-Brain-Gonad axis (GBG axis) and its importance on physiology and reproduction have already been highlighted for higher mammals; however, the study on fish physiology and reproduction is limited. While looking into the paucity of information, we have attempted to review the present status of microbiome and its interaction between the brain and gut. This review will address a process of the microbiome physiological mechanism involved in fish reproduction. The gut microbiota influences the BPG axis through a wide variety of compounds, including neuropeptides, neurotransmitter homologs and transmitters. Currently, research is being conducted to determine the precise process by which gut microbial composition influences brain function in fish. The gut-brain bidirectional interaction can influence brain biochemistry such as GABA, serotonin and tryptophan metabolites which play significant roles in CNS regulation. This review summarizes the fact, how microbes from gut, skin and other parts of the body influence fish reproduction through the Gut-Brain-Gonad axis.

Divergent and non-parallel evolution of MHC IIB in the Neotropical Midas cichlid species complex

Background

Ecological diversification is the result of divergent natural selection by contrasting habitat characteristics that favours the evolution of distinct phenotypes. This process can happen in sympatry and in allopatry. Habitat-specific parasite communities have the potential to drive diversification among host populations by imposing selective pressures on their host's immune system. In particular, the hyperdiverse genes of the major histocompatibility complex (MHC) are implicated in parasite-mediated host divergence. Here, we studied the extent of divergence at MHC, and discuss how it may have contributed to the Nicaraguan Midas cichlid species complex diversification, one of the most convincing examples of rapid sympatric parallel speciation.

Results

We genotyped the MHC IIB for individuals from six sympatric Midas cichlid assemblages, each containing species that have adapted to exploit similar habitats. We recovered large allelic and functional diversity within the species complex. While most alleles were rare, functional groups of alleles (supertypes) were common, suggesting that they are key to survival and that they were maintained during colonization and subsequent radiations. We identified lake-specific and habitat-specific signatures for both allelic and functional diversity, but no clear pattern of parallel divergence among ecomorphologically similar phenotypes.

Conclusions

Colonization and demographic effects of the fish could have contributed to MHC evolution in the Midas cichlid in conjunction with habitat-specific selective pressures, such as parasites associated to alternative preys or environmental features. Additional ecological data will help evaluating the role of host–parasite interactions in the Midas cichlid radiations and aid in elucidating the potential role of non-parallel features differentiating crater lake species assemblages.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01997-9.

Rapid adaptation through genomic and epigenomic responses following translocations in an endangered salmonid

Identifying the molecular mechanisms facilitating adaptation to new environments is a key question in evolutionary biology, especially in the face of current rapid and human-induced changes. Translocations have become an important tool for species conservation, but the attendant small population sizes and new ecological pressures might affect phenotypic and genotypic variation and trajectories dramatically and in unknown ways. In Scotland, the European whitefish (Coregonus lavaretus) is native to only two lakes and vulnerable to extirpation. Six new refuge populations were established over the last 30 years as a conservation measure. In this study, we examined whether there is a predictable ecological and evolutionary response of these fishes to translocation. We found eco-morphological differences, as functional traits relating to body shape differed between source and refuge populations. Dual isotopic analyses suggested some ecological release, with the diets in refuge populations being more diverse than in source populations. Analyses of up to 9117 genome-mapped SNPs showed that refuge populations had reduced genetic diversity and elevated inbreeding and relatedness relative to source populations, though genomic differentiation was low (F ST = 0.002-0.030). We identified 14 genomic SNPs that showed shared signals of a selective response to translocations, including some located near or within genes involved in the immune system, nervous system and hepatic functions. Analysis of up to 120,897 epigenomic loci identified a component of consistent differential methylation between source and refuge populations. We found that epigenomic variation and genomic variation were associated with morphological variation, but we were not able to infer an effect of population age because the patterns were also linked with the methodology of the translocations. These results show that conservation-driven translocations affect evolutionary potential by impacting eco-morphological, genomic and epigenomic components of diversity, shedding light on acclimation and adaptation process in these contexts.

A novel workflow to improve genotyping of multigene families in wildlife species: An experimental set-up with a known model system

Genotyping complex multigene families in novel systems is particularly challenging. Target primers frequently amplify simultaneously multiple loci leading to high PCR and sequencing artefacts such as chimeras and allele amplification bias. Most genotyping pipelines have been validated in nonmodel systems whereby the real genotype is unknown and the generation of artefacts may be highly repeatable. Further hindering accurate genotyping, the relationship between artefacts and genotype complexity (i.e. number of alleles per genotype) within a PCR remains poorly described. Here, we investigated the latter by experimentally combining multiple known major histocompatibility complex (MHC) haplotypes of a model organism (chicken, Gallus gallus, 43 artificial genotypes with 2-13 alleles per amplicon). In addition to well-defined 'optimal' primers, we simulated a nonmodel species situation by designing 'cross-species' primers based on sequence data from closely related Galliform species. We applied a novel open-source genotyping pipeline (ACACIA; https://gitlab.com/psc_santos/ACACIA), and compared its performance with another, previously published pipeline (AmpliSAS). Allele calling accuracy was higher when using ACACIA (98.5% versus 97% and 77.8% versus 75% for the 'optimal' and 'cross-species' data sets, respectively). Systematic allele dropout of three alleles owing to primer mismatch in the 'cross-species' data set explained high allele calling repeatability (100% when using ACACIA) despite low accuracy, demonstrating that repeatability can be misleading when evaluating genotyping workflows. Genotype complexity was positively associated with nonchimeric artefacts, chimeric artefacts (nonlinearly by levelling when amplifying more than 4-6 alleles) and allele amplification bias. Our study exemplifies and demonstrates pitfalls researchers should avoid to reliably genotype complex multigene families.

A continuously changing selective context on microbial communities associated with fish, from egg to fork

Fast increase of fish aquaculture production to meet consumer demands is accompanied by important ecological concerns such as disease outbreaks. Meanwhile, food waste is an important concern with fish products since they are highly perishable. Recent aquaculture and fish product microbiology, and more recently, microbiota research, paved the way to a highly integrated approach to understand complex relationships between host fish, product and their associated microbial communities at health/disease and preservation/spoilage frontiers. Microbial manipulation strategies are increasingly validated as promising tools either to replace or to complement traditional veterinary and preservation methods. In this review, we consider evolutionary forces driving fish microbiota assembly, in particular the changes in the selective context along the production chain. We summarize the current knowledge concerning factors governing assembly and dynamics of fish hosts and food microbial communities. Then, we discuss the current microbial community manipulation strategies from an evolutionary standpoint to provide a perspective on the potential for risks, conflict and opportunities. Finally, we conclude that to harness evolutionary forces in the development of sustainable microbiota manipulation applications in the fish industry, an integrated knowledge of the controlling abiotic and especially biotic factors is required.

A way forward with eco evo devo: an extended theory of resource polymorphism with postglacial fishes as model systems

A major goal of evolutionary science is to understand how biological diversity is generated and altered. Despite considerable advances, we still have limited insight into how phenotypic variation arises and is sorted by natural selection. Here we argue that an integrated view, which merges ecology, evolution and developmental biology (eco evo devo) on an equal footing, is needed to understand the multifaceted role of the environment in simultaneously determining the development of the phenotype and the nature of the selective environment, and how organisms in turn affect the environment through eco evo and eco devo feedbacks. To illustrate the usefulness of an integrated eco evo devo perspective, we connect it with the theory of resource polymorphism (i.e. the phenotypic and genetic diversification that occurs in response to variation in available resources). In so doing, we highlight fishes from recently glaciated freshwater systems as exceptionally well-suited model systems for testing predictions of an eco evo devo framework in studies of diversification. Studies on these fishes show that intraspecific diversity can evolve rapidly, and that this process is jointly facilitated by (i) the availability of diverse environments promoting divergent natural selection; (ii) dynamic developmental processes sensitive to environmental and genetic signals; and (iii) eco evo and eco devo feedbacks influencing the selective and developmental environments of the phenotype. We highlight empirical examples and present a conceptual model for the generation of resource polymorphism - emphasizing eco evo devo, and identify current gaps in knowledge.

Polymorphism of MHC class IIB in an acheilognathid species, Rhodeus sinensis shaped by historical selection and recombination

BACKGROUND

Rhodeus sinensis is a bitterling species occurring throughout the numerous freshwater systems on the East Asia. Here, we analyzed the diversity of the MHC class IIB (DAB) genes from this species, which may offer meaningful insights into evolutionary processes in this species as well as other bitterlings.

RESULTS

Using cDNA and gDNA samples from 50 individuals, we discovered classical 140 allelic sequences that could be allocated into either DAB1 (Rhsi-DAB1) or DAB3 (Rhsi-DAB3). DAB sequences completely lacking the intron, but identical or similar to Rhsi-DAB1, were also discovered from our gDNA samples, and this intron loss likely originated from the retrotransposition events of processed mDNA. The β1 domain was the most polymorphic in both Rhsi-DAB1 and -DAB3. Putative peptide biding residues (PBRs) in Rhsi-DAB1, but not in Rhsi-DAB3, exhibited a significant dN/dS, presumably indicating that different selection pressures have acted on those two DABs. Recombination between different alleles seemed to have contributed to the increase of diversity in Rhsi-DABs. Upon phylogenetic analysis, Rhsi-DAB1 and -DAB3 formed independent clusters. Several alleles from other species of Cypriniformes were embedded in the clade of Rhsi-DAB1, whereas Rhsi-DAB3 clustered with alleles from the wider range of taxa (Cyprinodontiformes), indicating that these two Rhsi-DABs have taken different historical paths.

CONCLUSIONS

A great deal of MHC class IIB allelic diversity was found in R. sinensis, and gene duplication, selection and recombination may have contributed to this diversity. Based on our data, it is presumed that such historical processes have commonly or differently acted on the polymorphism of Rhsi-DAB1 and -DAB3.

(Non)Parallel Evolution

Parallel evolution across replicate populations has provided evolutionary biologists with iconic examples of adaptation. When multiple populations colonize seemingly similar habitats, they may evolve similar genes, traits, or functions. Yet, replicated evolution in nature or in the laboratory often yields inconsistent outcomes: Some replicate populations evolve along highly similar trajectories, whereas other replicate populations evolve to different extents or in distinct directions. To understand these heterogeneous outcomes, biologists are increasingly treating parallel evolution not as a binary phenomenon but rather as a quantitative continuum ranging from parallel to nonparallel. By measuring replicate populations’ positions along this (non)parallel continuum, we can test hypotheses about evolutionary and ecological factors that influence the extent of repeatable evolution. We review evidence regarding the manifestation of (non)parallel evolution in the laboratory, in natural populations, and in applied contexts such as cancer. We enumerate the many genetic, ecological, and evolutionary processes that contribute to variation in the extent of parallel evolution.

Holobionts and ecological speciation: the intestinal microbiota of lake whitefish species pairs

BACKGROUND

It is well established that symbionts have considerable impact on their host, yet the investigation of the possible role of the holobiont in the host's speciation process is still in its infancy. In this study, we compared the intestinal microbiota among five sympatric pairs of dwarf (limnetic) and normal (benthic) lake whitefish Coregonus clupeaformis representing a continuum in the early stage of ecological speciation. We sequenced the 16s rRNA gene V3-V4 regions of the intestinal microbiota present in a total of 108 wild sympatric dwarf and normal whitefish as well as the water bacterial community from five lakes to (i) test for differences between the whitefish intestinal microbiota and the water bacterial community and (ii) test for parallelism in the intestinal microbiota of dwarf and normal whitefish.

RESULTS

The water bacterial community was distinct from the intestinal microbiota, indicating that intestinal microbiota did not reflect the environment, but rather the intrinsic properties of the host microbiota. Our results revealed a strong influence of the host (dwarf or normal) on the intestinal microbiota with pronounced conservation of the core intestinal microbiota (mean ~ 44% of shared genera). However, no clear evidence for parallelism was observed, whereby non-parallel differences between dwarf and normal whitefish were observed in three of the lakes while similar taxonomic composition was observed for the two other species pairs.

CONCLUSIONS

This absence of parallelism across dwarf vs. normal whitefish microbiota highlighted the complexity of the holobiont and suggests that the direction of selection could be different between the host and its microbiota.

Applied Hologenomics: Feasibility and Potential in Aquaculture

Aquaculture will play an essential role in feeding a growing human population, but several biological challenges impede sustainable growth of production. Emerging evidence across all areas of life has revealed the importance of the intimate biological interactions between animals and their associated gut microbiota. Based on challenges in aquaculture, we leverage current knowledge in molecular biology and host microbiota interactions to propose an applied holo-omic framework that integrates molecular data including genomes, transcriptomes, epigenomes, proteomes, and metabolomes for analyzing fish and their gut microbiota as interconnected and coregulated systems. With an eye towards aquaculture, we discuss the feasibility and potential of our holo-omic framework to improve growth, health, and sustainability in any area of food production, including livestock and agriculture.

Natural selection on MHC IIβ in parapatric lake and stream stickleback: Balancing, divergent, both or neither?

Major histocompatibility complex (MHC) genes encode proteins that play a central role in vertebrates' adaptive immunity to parasites. MHC loci are among the most polymorphic in vertebrates' genomes, inspiring many studies to identify evolutionary processes driving MHC polymorphism within populations and divergence between populations. Leading hypotheses include balancing selection favouring rare alleles within populations, and spatially divergent selection. These hypotheses do not always produce diagnosably distinct predictions, causing many studies of MHC to yield inconsistent or ambiguous results. We suggest a novel strategy to distinguish balancing vs. divergent selection on MHC, taking advantage of natural admixture between parapatric populations. With divergent selection, individuals with immigrant alleles will be more infected and less fit because they are susceptible to novel parasites in their new habitat. With balancing selection, individuals with locally rare immigrant alleles will be more fit (less infected). We tested these contrasting predictions using three-spine stickleback from three replicate pairs of parapatric lake and stream habitats. We found numerous positive and negative associations between particular MHC IIβ alleles and particular parasite taxa. A few allele-parasite comparisons supported balancing selection, and others supported divergent selection between habitats. But, there was no overall tendency for fish with immigrant MHC alleles to be more or less heavily infected. Instead, locally rare MHC alleles (not necessarily immigrants) were associated with heavier infections. Our results illustrate the complex relationship between MHC IIβ allelic variation and spatially varying multispecies parasite communities: different hypotheses may be concurrently true for different allele-parasite combinations.

The evolution of the major histocompatibility complex in upstream versus downstream river populations of the longnose dace

Populations in upstream versus downstream river locations can be exposed to vastly different environmental and ecological conditions and can thus harbor different genetic resources due to selection and neutral processes. An interesting question is how upstream–downstream directionality in rivers affects the evolution of immune response genes. We used next‐generation amplicon sequencing to identify eight alleles of the major histocompatibility complex (MHC) class II β exon 2 in the cyprinid longnose dace (Rhinichthys cataractae) from three rivers in Alberta, upstream and downstream of municipal and agricultural areas along contaminant gradients. We used these data to test for directional and balancing selection on the MHC. We also genotyped microsatellite loci to examine neutral population processes in this system. We found evidence for balancing selection on the MHC in the form of increased nonsynonymous variation relative to neutral expectations, and selection occurred at more amino acid residues upstream than downstream in two rivers. We found this pattern despite no population structure or isolation by distance, based on microsatellite data, at these sites. Overall, our results suggest that MHC evolution is driven by upstream–downstream directionality in fish inhabiting this system.

Microsatellite and major histocompatibility complex variation in an endangered rattlesnake, the Eastern Massasauga (Sistrurus catenatus)

Genetic diversity is fundamental to maintaining the long-term viability of populations, yet reduced genetic variation is often associated with small, isolated populations. To examine the relationship between demography and genetic variation, variation at hypervariable loci (e.g., microsatellite DNA loci) is often measured. However, these loci are selectively neutral (or near neutral) and may not accurately reflect genomewide variation. Variation at functional trait loci, such as the major histocompatibility complex (MHC), can provide a better assessment of adaptive genetic variation in fragmented populations. We compared patterns of microsatellite and MHC variation across three Eastern Massasauga (Sistrurus catenatus) populations representing a gradient of demographic histories to assess the relative roles of natural selection and genetic drift. Using 454 deep amplicon sequencing, we identified 24 putatively functional MHC IIB exon 2 alleles belonging to a minimum of six loci. Analysis of synonymous and nonsynonymous substitution rates provided evidence of historical positive selection at the nucleotide level, and Tajima's D provided support for balancing selection in each population. As predicted, estimates of microsatellite allelic richness, observed, heterozygosity, and expected heterozygosity varied among populations in a pattern qualitatively consistent with demographic history and abundance. While MHC allelic richness at the population and individual levels revealed similar trends, MHC nucleotide diversity was unexpectedly high in the smallest population. Overall, these results suggest that genetic variation in the Eastern Massasauga populations in Illinois has been shaped by multiple evolutionary mechanisms. Thus, conservation efforts should consider both neutral and functional genetic variation when managing captive and wild Eastern Massasauga populations.

Next-generation genotyping of hypervariable loci in many individuals of a non-model species: technical and theoretical implications

BACKGROUND

Across species, diversity at the Major Histocompatibility Complex (MHC) is critical to disease resistance and population health; however, use of MHC diversity to quantify the genetic health of populations has been hampered by the extreme variation found in MHC genes. Next generation sequencing (NGS) technology generates sufficient data to genotype even the most diverse species, but workflows for distinguishing artifacts from alleles are still under development. We used NGS to evaluate the MHC diversity of over 300 captive and wild ring-tailed lemurs (Lemur catta: Primates: Mammalia). We modified a published workflow to address errors that arise from deep sequencing individuals and tested for evidence of selection at the most diverse MHC genes.

RESULTS

In addition to evaluating the accuracy of 454 Titanium and Ion Torrent PGM for genotyping large populations at hypervariable genes, we suggested modifications to improve current methods of allele calling. Using these modifications, we genotyped 302 out of 319 individuals, obtaining an average sequencing depth of over 1000 reads per amplicon. We identified 55 MHC-DRB alleles, 51 of which were previously undescribed, and provide the first sequences of five additional MHC genes: DOA, DOB, DPA, DQA, and DRA. The additional five MHC genes had one or two alleles each with little sequence variation; however, the 55 MHC-DRB alleles showed a high dN/dS ratio and trans-species polymorphism, indicating a history of positive selection. Because each individual possessed 1-7 MHC-DRB alleles, we suggest that ring-tailed lemurs have four, putatively functional, MHC-DRB copies.

CONCLUSIONS

In the future, accurate genotyping methods for NGS data will be critical to assessing genetic variation in non-model species. We recommend that future NGS studies increase the proportion of replicated samples, both within and across platforms, particularly for hypervariable genes like the MHC. Quantifying MHC diversity within non-model species is the first step to assessing the relationship of genetic diversity at functional loci to individual fitness and population viability. Owing to MHC-DRB diversity and copy number, ring-tailed lemurs may serve as an ideal model for estimating the interaction between genetic diversity, fitness, and environment, especially regarding endangered species.

Characterization of MHC class IIB for four endangered Australian freshwater fishes obtained from ecologically divergent populations

Genetic diversity is an essential aspect of species viability, and assessments of neutral genetic diversity are regularly implemented in captive breeding and conservation programs. Despite their importance, information from adaptive markers is rarely included in such programs. A promising marker of significance in fitness and adaptive potential is the major histocompatibility complex (MHC), a key component of the adaptive immune system. Populations of Australian freshwater fishes are generally declining in numbers due to human impacts and the introduction of exotic species, a scenario of particular concern for members of the family Percichthyidae, several of which are listed as nationally vulnerable or endangered, and hence subject to management plans, captive breeding, and restoration plans. We used a next-generation sequencing approach to characterize the MHC IIB locus and provide a conservative description of its levels of diversity in four endangered percichthyids: Gadopsis marmoratus, Macquaria australasica, Nannoperca australis, and Nannoperca obscura. Evidence is presented for a duplicated MHC IIB locus, positively selected sites and recombination of MHC alleles. Relatively moderate levels of diversity were detected in the four species, as well as in different ecotypes within each species. Phylogenetic analyses revealed genus specific clustering of alleles and no allele sharing among species. There were also no shared alleles observed between two ecotypes within G. marmoratus and within M. australasica, which might be indicative of ecologically-driven divergence and/or long divergence times. This represents the first characterization and assessment of MHC diversity for Percichthyidae, and also for Australian freshwater fishes in general, providing key genetic resources for a vertebrate group of increasing conservation concern.

Among-lake reciprocal transplants induce convergent expression of immune genes in threespine stickleback

Geographic variation in parasite communities can drive evolutionary divergence in host immune genes. However, biotic and abiotic environmental variation can also induce plastic differences in immune function among populations. At present, there is little information concerning the relative magnitudes of heritable vs. induced immune divergence in natural populations. We examined immune gene expression profiles of threespine stickleback (Gasterosteus aculeatus) from six lakes on Vancouver Island, British Columbia. Parasite community composition differs between lake types (large or small, containing limnetic- or benthic-like stickleback) and between watersheds. We observed corresponding differences in immune gene expression profiles among wild-caught stickleback, using a set of seven immune genes representing distinct branches of the immune system. To evaluate the role of environmental effects on this differentiation, we experimentally transplanted wild-caught fish into cages in their native lake, or into a nearby foreign lake. Transplanted individuals' immune gene expression converged on patterns typical of their destination lake, deviating from their native expression profile. Transplant individuals' source population had a much smaller effect, suggesting relatively weak genetic underpinning of population differences in immunity, as viewed through gene expression. This strong environmental regulation of immune gene expression provides a counterpoint to the large emerging literature documenting microevolution and genetic diversification of immune function. Our findings illustrate the value of studying immunity in natural environmental settings where the immune system has evolved and actively functions.

Identifying designatable units for intraspecific conservation prioritization: a hierarchical approach applied to the lake whitefish species complex (Coregonus spp.)

The concept of the designatable unit (DU) affords a practical approach to identifying diversity below the species level for conservation prioritization. However, its suitability for defining conservation units in ecologically diverse, geographically widespread and taxonomically challenging species complexes has not been broadly evaluated. The lake whitefish species complex (Coregonus spp.) is geographically widespread in the Northern Hemisphere, and it contains a great deal of variability in ecology and evolutionary legacy within and among populations, as well as a great deal of taxonomic ambiguity. Here, we employ a set of hierarchical criteria to identify DUs within the Canadian distribution of the lake whitefish species complex. We identified 36 DUs based on (i) reproductive isolation, (ii) phylogeographic groupings, (iii) local adaptation and (iv) biogeographic regions. The identification of DUs is required for clear discussion regarding the conservation prioritization of lake whitefish populations. We suggest conservation priorities among lake whitefish DUs based on biological consequences of extinction, risk of extinction and distinctiveness. Our results exemplify the need for extensive genetic and biogeographic analyses for any species with broad geographic distributions and the need for detailed evaluation of evolutionary history and adaptive ecological divergence when defining intraspecific conservation units.

Spatial and temporal variation at major histocompatibility complex class IIB genes in the endangered Blakiston's fish owl

INTRODUCTION

Quantifying intraspecific genetic variation in functionally important genes, such as those of the major histocompatibility complex (MHC), is important in the establishment of conservation plans for endangered species. The MHC genes play a crucial role in the vertebrate immune system and generally show high levels of diversity, which is likely due to pathogen-driven balancing selection. The endangered Blakiston's fish owl (Bubo blakistoni) has suffered marked population declines on Hokkaido Island, Japan, during the past several decades due to human-induced habitat loss and fragmentation. We investigated the spatial and temporal patterns of genetic diversity in MHC class IIβ genes in Blakiston's fish owl, using massively parallel pyrosequencing.

RESULTS

We found that the Blakiston's fish owl genome contains at least eight MHC class IIβ loci, indicating recent gene duplications. An analysis of sequence polymorphism provided evidence that balancing selection acted in the past. The level of MHC variation, however, was low in the current fish owl populations in Hokkaido: only 19 alleles were identified from 174 individuals. We detected considerable spatial differences in MHC diversity among the geographically isolated populations. We also detected a decline of MHC diversity in some local populations during the past decades.

CONCLUSIONS

Our study demonstrated that the current spatial patterns of MHC variation in Blakiston's fish owl populations have been shaped by loss of variation due to the decline and fragmentation of populations, and that the short-term effects of genetic drift have counteracted the long-term effects of balancing selection.

High intralocus variability and interlocus recombination promote immunological diversity in a minimal major histocompatibility system

Background

The genes of the major histocompatibility complex (MHC/MH) have attracted considerable scientific interest due to their exceptional levels of variability and important function as part of the adaptive immune system. Despite a large number of studies on MH class II diversity of both model and non-model organisms, most research has focused on patterns of genetic variability at individual loci, failing to capture the functional diversity of the biologically active dimeric molecule. Here, we take a systematic approach to the study of MH variation, analyzing patterns of genetic variation at MH class IIα and IIβ loci of the seahorse, which together form the immunologically active peptide binding cleft of the MH class II molecule.

Results

The seahorse carries a minimal class II system, consisting of single copies of both MH class IIα and IIβ, which are physically linked and inherited in a Mendelian fashion. Both genes are ubiquitously expressed and detectible in the brood pouch of male seahorses throughout pregnancy. Genetic variability of the two genes is high, dominated by non-synonymous variation concentrated in their peptide-binding regions. Coding variation outside these regions is negligible, a pattern thought to be driven by intra- and interlocus recombination. Despite the tight physical linkage of MH IIα and IIβ loci, recombination has produced novel composite alleles, increasing functional diversity at sites responsible for antigen recognition.

Conclusions

Antigen recognition by the adaptive immune system of the seahorse is enhanced by high variability at both MH class IIα and IIβ loci. Strong positive selection on sites involved in pathogen recognition, coupled with high levels of intra- and interlocus recombination, produce a patchwork pattern of genetic variation driven by genetic hitchhiking. Studies focusing on variation at individual MH loci may unintentionally overlook an important component of ecologically relevant variation.

Stepwise threshold clustering: a new method for genotyping MHC loci using next-generation sequencing technology

Genes of the vertebrate major histocompatibility complex (MHC) are of great interest to biologists because of their important role in immunity and disease, and their extremely high levels of genetic diversity. Next generation sequencing (NGS) technologies are quickly becoming the method of choice for high-throughput genotyping of multi-locus templates like MHC in non-model organisms. Previous approaches to genotyping MHC genes using NGS technologies suffer from two problems:1) a "gray zone" where low frequency alleles and high frequency artifacts can be difficult to disentangle and 2) a similar sequence problem, where very similar alleles can be difficult to distinguish as two distinct alleles. Here were present a new method for genotyping MHC loci--Stepwise Threshold Clustering (STC)--that addresses these problems by taking full advantage of the increase in sequence data provided by NGS technologies. Unlike previous approaches for genotyping MHC with NGS data that attempt to classify individual sequences as alleles or artifacts, STC uses a quasi-Dirichlet clustering algorithm to cluster similar sequences at increasing levels of sequence similarity. By applying frequency and similarity based criteria to clusters rather than individual sequences, STC is able to successfully identify clusters of sequences that correspond to individual or similar alleles present in the genomes of individual samples. Furthermore, STC does not require duplicate runs of all samples, increasing the number of samples that can be genotyped in a given project. We show how the STC method works using a single sample library. We then apply STC to 295 threespine stickleback (Gasterosteus aculeatus) samples from four populations and show that neighboring populations differ significantly in MHC allele pools. We show that STC is a reliable, accurate, efficient, and flexible method for genotyping MHC that will be of use to biologists interested in a variety of downstream applications.

Neutral and selective processes shape MHC gene diversity and expression in stocked brook charr populations (Salvelinus fontinalis)

The capacity of an individual to battle infection is an important fitness determinant in wild vertebrate populations. The major histocompatibility complex (MHC) genes are crucial for a host's adaptive immune system to detect pathogens. However, anthropogenic activities may disrupt natural cycles of co-evolution between hosts and pathogens. In this study, we investigated the dynamic sequence and expression variation of host parasite interactions in brook charr (Salvelinus fontinalis) in a context of past human disturbance via population supplementation from domestic individuals. To do so, we developed a new method to examine selection shaping MHC diversity within and between populations and found a complex interplay between neutral and selective processes that varied between lakes that were investigated. We provided evidence for a lower introgression rate of domestic alleles and found that parasite infection increased with domestic genomic background of individuals. We also documented an association between individual MHC alleles and parasite taxa. Finally, longer cis-regulatory minisatellites were positively correlated with MHC II down-regulation and domestic admixture, suggesting that inadvertent selection during domestication resulted in a lower immune response capacity, through a trade-off between growth and immunity, which explained the negative selection of domestic alleles at least under certain circumstances.

Teleost microbiomes: the state of the art in their characterization, manipulation and importance in aquaculture and fisheries

Indigenous microbiota play a critical role in the lives of their vertebrate hosts. In human and mouse models it is increasingly clear that innate and adaptive immunity develop in close concert with the commensal microbiome. Furthermore, several aspects of digestion and nutrient metabolism are governed by intestinal microbiota. Research on teleosts has responded relatively slowly to the introduction of massively parallel sequencing procedures in microbiomics. Nonetheless, progress has been made in biotic and gnotobiotic zebrafish models, defining a core microbiome and describing its role in development. However, microbiome research in other teleost species, especially those important from an aquaculture perspective, has been relatively slow. In this review, we examine progress in teleost microbiome research to date. We discuss teleost microbiomes in health and disease, microbiome ontogeny, prospects for successful microbiome manipulation (especially in an aquaculture setting) and attempt to identify important future research themes. We predict an explosion in research in this sector in line with the increasing global demand for fish protein, and the need to find sustainable approaches to improve aquaculture yield. The reduced cost and increasing ease of next generation sequencing technologies provides the technological backing, and the next 10 years will be an exciting time for teleost microbiome research.

Ultra-deep Illumina sequencing accurately identifies MHC class IIb alleles and provides evidence for copy number variation in the guppy (Poecilia reticulata)

We address the bioinformatic issue of accurately separating amplified genes of the major histocompatibility complex (MHC) from artefacts generated during high-throughput sequencing workflows. We fit observed ultra-deep sequencing depths (hundreds to thousands of sequences per amplicon) of allelic variants to expectations from genetic models of copy number variation (CNV). We provide a simple, accurate and repeatable method for genotyping multigene families, evaluating our method via analyses of 209 b of MHC class IIb exon 2 in guppies (Poecilia reticulata). Genotype repeatability for resequenced individuals (N = 49) was high (100%) within the same sequencing run. However, repeatability dropped to 83.7% between independent runs, either because of lower mean amplicon sequencing depth in the initial run or random PCR effects. This highlights the importance of fully independent replicates. Significant improvements in genotyping accuracy were made by greatly reducing type I genotyping error (i.e. accepting an artefact as a true allele), which may occur when using low-depth allele validation thresholds used by previous methods. Only a small amount (4.9%) of type II error (i.e. rejecting a genuine allele as an artefact) was detected through fully independent sequencing runs. We observed 1-6 alleles per individual, and evidence of sharing of alleles across loci. Variation in the total number of MHC class II loci among individuals, both among and within populations was also observed, and some genotypes appeared to be partially hemizygous; total allelic dosage added up to an odd number of allelic copies. Collectively, observations provide evidence of MHC CNV and its complex basis in natural populations.