Genetic variability, local selection and demographic history: genomic evidence of evolving towards allopatric speciation in Asian seabass

Genomic structural variation in Barramundi Perch Lates calcarifer and potential roles in speciation and adaptation

Advancements in genome sequencing and assembly techniques have increased the documentation of structural variants in wild organisms. Of these variants, chromosomal inversions are especially prominent due to their large size and active recombination suppression between alternative homokaryotypes. This suppression enables the 2 forms of the inversion to be maintained and allows the preservation of locally adapted alleles. The Barramundi Perch (BP; Lates calcarifer) is a widespread species complex with 3 main genetic lineages located in the biogeographic regions of Australia and New Guinea (AUS + NG), Southeast Asia (SEA), and the Indian Subcontinent (IND). BP are typically considered to be a protandrous sequential hermaphrodite species that exhibits catadromy. Freshwater occupancy and intraspecific variation in life history (e.g. partially migratory populations) exist and provide opportunities for strongly divergent selection associated with, for example, salinity tolerance, swimming ability, and marine dispersal. Herein, we utilize genomic data generated from all 3 genetic lineages to identify and describe 3 polymorphic candidate chromosomal inversions. These candidate chromosomal inversions appear to be fixed for ancestral variants in the IND lineage and for inverted versions in the AUS + NG lineage and exhibit variation in all 3 inversions in the SEA lineage. BP have a diverse portfolio of life history options that includes migratory strategy as well as sexual system (i.e. hermaphroditism and gonochorism). We propose that the some of the life history variabilities observed in BP may be linked to inversions and, in doing so, we present genetic data that might be useful in enhancing aquaculture production and population management.

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

A Genomic Quantitative Study on the Contribution of the Ancestral-State Bases Relative to Derived Bases in the Divergence and Local Adaptation of Populus davidiana

Identifying alleles associated with adaptation to new environments will advance our understanding of evolution from the molecular level. Previous studies have found that the Populus davidiana southwest population in East Asia has differentiated from other populations in the range. We aimed to evaluate the contributions of the ancestral-state bases (ASBs) relative to derived bases (DBs) in the local adaptation of P. davidiana in the Yunnan-Guizhou Plateau from a quantitative perspective using whole-genome re-sequencing data from 90 P. davidiana samples from three regions across the species range. Our results showed that the uplift of the Qinghai-Tibet Plateau during the Neogene and associated climate fluctuations during the Middle Pleistocene were likely an important factor in the early divergence of P. davidiana. Highly differentiated genomic regions between populations were inferred to have undergone strong linked natural selection, and ASBs are the chief means by which populations of P. davidiana adapt to novel environmental conditions; however, when adapting to regions with high environmental differences relative to the ancestral range, the proportion of DBs was significantly higher than that of background regions, as ASBs are insufficient to cope with these environments. Finally, a number of genes were identified in the outlier region.

Linked selection, differential introgression and recombination rate variation promote heterogeneous divergence in a pair of yellow croakers

Understanding the mechanisms underlying heterogeneous genomic divergence is of particular interest in evolutionary biology. Highly differentiated genomic regions, known as genomic islands, often evolve between diverging lineages. These genomic islands may be related to selection promoting adaptation or reproductive isolation. Based on whole genome assembly and genome-wide RAD sequencing in a pair of yellow croakers (genus: Larimichthys), we investigated the evolutionary processes shaping genomic landscapes of divergence. Demographic modelling indicated that the two species diverged following a secondary contact scenario, where differential introgression and linked selection were suggested to be involved in heterogeneous genomic divergence. We identified reduced recombination rate in genomic islands and a relatively good conservation of both genetic diversity and recombination landscapes between species, which highlight the roles of linked selection and recombination rate variation in promoting heterogeneous divergence in the common ancestral lineage of the two species. In addition, we found a positive correlation between differentiation (F_ST ) and absolute sequence divergence (D_xy ), and elevated D_xy in genomic islands, indicating that the genomic landscape of divergence was not shaped by linked selection alone. Restricted gene flow in highly differentiated regions has probably remodelled the landscape of heterogeneous genomic divergence. This study highlights that highly differentiated genomic regions can also arise from a combination of linked selection and differential gene flow in interaction with varying recombination rates.

Worldwide Population Structure of the Coffee Rust Fungus Hemileia vastatrix Is Strongly Shaped by Local Adaptation and Breeding History

The devastating disease coffee leaf rust, caused by Hemileia vastatrix, has been a major constraint to worldwide coffee production. Recently, H. vastatrix populations were shown to be structured into three divergent genetic lineages with marked host specialization (C1, C2, and C3). However, there is yet no overall understanding of the population dynamics and adaptation of the most widespread and epidemiological relevant H. vastatrix group (C3). We used restriction site-associated DNA sequencing to generate 13,804 single nucleotide polymorphisms (SNPs) across a worldwide collection of 99 H. vastatrix isolates. Phylogenetic analyses uncovered a well-supported structuring within C3, with three main subgroups (SGs; SGI, SGII, and SGIII), which seem to reflect the historical distribution, breeding, and exchange of coffee cultivars. SGI shows a ladder-like diversification pattern and occurs across all four continents sampled, SGII is mainly restricted to Africa, and SGIII is observed only in Timor, revealing a higher genetic differentiation. Outlier and association tests globally identified 112 SNPs under putative positive selection, which impacted population structure. In particular, 29 overlapping SNPs per se seemed to have an extremely strong effect on H. vastatrix population divergence. We also found exclusive and fixed alleles associated with the SGs supporting local adaptation. Functional annotation revealed that transposable elements may play a role in host adaptation. Our study provides a higher-resolution perspective on the evolutionary history of H. vastatrix on cultivated coffee, showing its strong ability to adapt and the strength of the selective force imposed by coffee hosts, which should be taken into account when designing strategies for pathogen dissemination control and selective breeding.

Restoring Genetic Resource through In Vitro Culturing Testicular Cells from the Cryo-Preserved Tissue of the American Shad (Alosa sapidissima)

Simple Summary

Cryopreservation and in vitro culture of germ cells are key techniques for the genetic resource preservation of the declining population of American shad. Two types of cryopreserved samples, namely testis pieces and testicular cells of American shad, were comparatively analyzed for cell viability. The results showed that the cell viability of the cryopreserved testis pieces was much higher than that of the cryopreserved testicular cells. The viability of the cells from the cryopreserved testis pieces ranged from 65.2 ± 2.2 (%) to 93.8 ± 0.6 (%), whereas the viability of the dissociated cells after cryopreservation was 38.5 ± 0.8 (%) to 87.1 ± 2.6 (%). Moreover, the testicular cells isolated from the post-thaw testicular tissue could be cultured and propagated in vitro. Our findings would benefit further investigations on genetic resource preservation and other manipulations of germ cells in a commercially and ecologically important fish species.

Abstract

Germ cells, as opposed to somatic cells, can transmit heredity information between generations. Cryopreservation and in vitro culture of germ cells are key techniques for genetic resource preservation and cellular engineering breeding. In this study, two types of cryopreserved samples, namely testis pieces and testicular cells of American shad, were comparatively analyzed for cell viability. The results showed that the cell viability of the cryopreserved testis pieces was much higher than that of the cryopreserved testicular cells. The viability of cells from the cryopreserved testis pieces ranged from 65.2 ± 2.2 (%) to 93.8 ± 0.6 (%), whereas the viability of the dissociated cells after cryopreservation was 38.5 ± 0.8 (%) to 87.1 ± 2.6 (%). Intriguingly, the testicular cells from the post-thaw testicular tissue could be cultured in vitro. Likewise, most of the cultured cells exhibited germ cell properties and highly expressed Vasa and PCNA protein. This study is the first attempt to effectively preserve and culture the male germ cells through freezing tissues in the American shad. The findings of this study would benefit further investigations on genetic resource preservation and other manipulations of germ cells in a commercially and ecologically important fish species.

Genome-wide analysis of butterfly bush (Buddleja alternifolia) in three uplands provides insights into biogeography, demography and speciation

Understanding processes that generate and maintain large disjunctions within plant species can provide valuable insights into plant diversity and speciation. The butterfly bush Buddleja alternifolia has an unusual disjunct distribution, occurring in the Himalaya, Hengduan Mountains (HDM) and the Loess Plateau (LP) in China. We generated a high-quality, chromosome-level genome assembly of B. alternifolia, the first within the family Scrophulariaceae. Whole-genome re-sequencing data from 48 populations plus morphological and petal colour reflectance data covering its full distribution range were collected. Three distinct genetic lineages of B. alternifolia were uncovered, corresponding to Himalayan, HDM and LP populations, with the last also differentiated morphologically and phenologically, indicating occurrence of allopatric speciation likely to be facilitated by geographic isolation and divergent adaptation to distinct ecological niches. Moreover, speciation with gene flow between populations from either side of a mountain barrier could be under way within LP. The current disjunctions within B. alternifolia might result from vicariance of a once widespread distribution, followed by several past contraction and expansion events, possibly linked to climate fluctuations promoted by the Kunlun-Yellow river tectonic movement. Several adaptive genes are likely to be either uniformly or diversely selected among regions, providing a footprint of local adaptations. These findings provide new insights into plant biogeography, adaptation and different processes of allopatric speciation.

Complex and divergent histories gave rise to genome-wide divergence patterns amongst European whitefish (Coregonus lavaretus)

Pleistocene glaciations dramatically affected species distribution in regions that were impacted by ice cover and subsequent postglacial range expansion impacted contemporary biodiversity in complex ways. The European whitefish, Coregonus lavaretus, is a widely distributed salmonid fish species on mainland Europe, but in Britain it has only seven native populations, all of which are found on the western extremes of the island. The origins and colonization routes of the species into Britain are unknown but likely contributed to contemporary genetic patterns and regional uniqueness. Here, we used up to 25,751 genome‐wide polymorphic loci to reconstruct the history and to discern the demographic and evolutionary forces underpinning divergence between British populations. Overall, we found lower genetic diversity in Scottish populations but high differentiation (F_ST = 0.433–0.712) from the English/Welsh and other European populations. Differentiation was elevated genome‐wide rather than in particular genomic regions. Demographic modelling supported a postglacial colonization into western Scotland from northern refugia and a separate colonization route for the English/Welsh populations from southern refugia, with these two groups having been separated for more than ca. 50 Ky. We found cyto‐nuclear discordance at a European scale, with the Scottish populations clustering closely with Baltic population in the mtDNA analysis but not in the nuclear data, and with the Norwegian and Alpine populations displaying the same mtDNA haplotype but being distantly related in the nuclear tree. These findings suggest that neutral processes, primarily drift and regionally distinct pre‐glacial evolutionary histories, are important drivers of genomic divergence in British populations of European whitefish. This sheds new light on the establishment of the native British freshwater fauna after the last glacial maximum.

Robust and efficient software for reference-free genomic diversity analysis of genotyping-by-sequencing data on diploid and polyploid species

Genotyping-by-sequencing (GBS) is a widely used and cost-effective technique for obtaining large numbers of genetic markers from populations by sequencing regions adjacent to restriction cut sites. Although a standard reference-based pipeline can be followed to analyse GBS reads, a reference genome is still not available for a large number of species. Hence, reference-free approaches are required to generate the genetic variability information that can be obtained from a GBS experiment. Unfortunately, available tools to perform de novo analysis of GBS reads face issues of usability, accuracy and performance. Furthermore, few available tools are suitable for analysing data sets from polyploid species. In this manuscript, we describe a novel algorithm to perform reference-free variant detection and genotyping from GBS reads. Nonexact searches on a dynamic hash table of consensus sequences allow for efficient read clustering and sorting. This algorithm was integrated in the Next Generation Sequencing Experience Platform (NGSEP) to integrate the state-of-the-art variant detector already implemented in this tool. We performed benchmark experiments with three different empirical data sets of plants and animals with different population structures and ploidies, and sequenced with different GBS protocols at different read depths. These experiments show that NGSEP has comparable and in some cases better accuracy and always better computational efficiency compared to existing solutions. We expect that this new development will be useful for many research groups conducting population genetic studies in a wide variety of species.

Effects of Ocean Acidification on Transcriptomes in Asian Seabass Juveniles

Ocean acidification is changing the fate of marine organisms. It is essential to predict the biological responses and evolutionary processes driven by ocean acidification, to maintain the equilibrium of the marine ecosystem and to facilitate aquaculture. However, how marine organisms, particularly the marine fish species, respond to ocean acidification, is still poorly understood. Consequences of ocean acidification on finfish aquaculture are largely not well known. We studied the effects of ocean acidification for 7 days on growth, behaviour and gene expression profiles in the brain, gill and kidney of Asian seabass juveniles. Results showed that growth and behaviour were not affected by short-term ocean acidification. We found tissue-specific differentially expressed genes (DEGs) involving many molecular processes, such as organ development, growth, muscle development, ion homeostasis and neurogenesis and development, as well as behaviours. Most of the DEGs, which were functionally enriched in ion homeostasis, were related to calcium transport, followed by sodium/potassium channels. We found that genes associated with neurogenesis and development were significantly enriched, implying that ocean acidification has also adversely affected the neural regulatory mechanism. Our results indicate that although the short-term ocean acidification does not cause obvious phenotypic and behavioural changes, it causes substantial changes of gene expressions in all three analysed tissues. All these changes of gene expressions may eventually affect physiological fitness. The DEGs identified here should be further investigated to discover DNA markers associated with adaptability to ocean acidification to improve fish's capability to adapt to ocean acidification.

Whole-genome re-sequencing data to infer historical demography and speciation processes in land snails: the study of two Candidula sister species

Despite the global biodiversity of terrestrial gastropods and their ecological and economic importance, the genomic basis of ecological adaptation and speciation in land snail taxa is still largely unknown. Here, we combined whole-genome re-sequencing with population genomics to evaluate the historical demography and the speciation process of two closely related species of land snails from western Europe, Candidula unifasciata and C. rugosiuscula. Historical demographic analysis indicated fluctuations in the size of ancestral populations, probably driven by Pleistocene climatic fluctuations. Although the current population distributions of both species do not overlap, our approximate Bayesian computation model selection approach on several speciation scenarios suggested that gene flow has occurred throughout the divergence process until recently. Positively selected genes diverging early in the process were associated with intragenomic and cyto-nuclear incompatibilities, respectively, potentially fostering reproductive isolation as well as ecological divergence. Our results suggested that the speciation between species entails complex processes involving both gene flow and ecological speciation, and that further research based on whole-genome data can provide valuable understanding on species divergence. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Mapping and Validating QTL for Fatty Acid Compositions and Growth Traits in Asian Seabass

Asian seabass is an important food fish species. While improving growth, increasing the nutritional value is important, omega-3 fatty acids are indispensable to human health. Identifying and validating DNA markers associated with traits is the first step towards marker-assisted selection (MAS). We quantified 13 different fatty acids and three growth traits in 213 F₂ Asian seabass from a family at the age 270 days post hatch, and screened QTL for these traits. The content of total fatty acids in 100 g flesh was 2.57 ± 0.80 g, while the proportions of docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA) were 16.96 ± 2.20% and 5.42 ± 0.90%, respectively. A linkage map with 2424 SNPs was constructed and used for QTL mapping. For fatty acid compositions, 14 significant QTL were identified on three linkage groups (LG5, LG11 and LG14), with phenotypic variance explained (PVE) from 12.8 to 24.6%. Thirty-nine suggestive QTL were detected on 16 LGs. Two significant QTL for EPA were identified on LG5 and LG14, with PVE of 15.2% and 15.1%, respectively. No significant QTL was identified for DHA. For growth traits, six significant and 13 suggestive QTL were identified on two and seven LGs, respectively. Only a few significant QTL for fatty acids overlapped with previously mapped QTL for these traits, suggesting that most QTL detected in a family are family-specific and could only be used in MAS in the family per se. To facilitate population-wide molecular breeding, more powerful methods (e.g. GWAS) should be used to identify SNPs for genomic selection.

Population Genomic Signatures of Genetic Structure and Environmental Selection in the Catadromous Roughskin Sculpin Trachidermus fasciatus

Understanding the patterns of genetic diversity and adaptation across species' range is crucial to assess its long-term persistence and determine appropriate conservation measures. The impacts of human activities on the genetic diversity and genetic adaptation to heterogeneous environments remain poorly understood in the marine realm. The roughskin sculpin (Trachidermus fasciatus) is a small catadromous fish, and has been listed as a second-class state protected aquatic animal since 1988 in China. To elucidate the underlying mechanism of population genetic structuring and genetic adaptations to local environments, RAD tags were sequenced for 202 individuals in nine populations across the range of T. fasciatus in China. The pairwise FST values over 9,271 filtered SNPs were significant except that between Dongying and Weifang. All the genetic clustering analysis revealed significant population structure with high support for eight distinct genetic clusters. Both the minor allele frequency spectra and Ne estimations suggested extremely small Ne in some populations (e.g., Qinhuangdao, Rongcheng, Wendeng, and Qingdao), which might result from recent population bottleneck. The strong genetic structure can be partly attributed to genetic drift and habitat fragmentation, likely due to the anthropogenic activities. Annotations of candidate adaptive loci suggested that genes involved in metabolism, development, and osmoregulation were critical for adaptation to spatially heterogenous environment of local populations. In the context of anthropogenic activities and environmental change, results of the present population genomic work provided important contributions to the understanding of genetic differentiation and adaptation to changing environments.

Population structure, demographic history and local adaptation of the grass carp

BACKGROUND

Genetic diversity within a species reflects population evolution, ecology, and ability to adapt. Genome-wide population surveys of both natural and introduced populations provide insights into genetic diversity, the evolutionary processes and the genetic basis underlying local adaptation. Grass carp is the most important freshwater foodfish species for food and water weed control. However, there is as yet no overall picture on genetic variations and population structure of this species, which is important for its aquaculture.

RESULTS

We used 43,310 SNPs to infer the population structure, evidence of local adaptation and sources of introduction. The overall genetic differentiation of this species was low. The native populations were differentiated into three genetic clusters, corresponding to the Yangtze, Pearl and Heilongjiang River Systems, respectively. The populations in Malaysia, India and Nepal were introduced from both the Yangtze and Pearl River Systems. Loci and genes involved in putative local selection for native locations were identified. Evidence of both positive and balancing selection was found in the introduced locations. Genes associated with loci under putative selection were involved in many biological functions. Outlier loci were grouped into clusters as genomic islands within some specific genomic regions, which likely agrees with the divergence hitchhiking scenario of divergence-with-gene-flow.

CONCLUSIONS

This study, for the first time, sheds novel insights on the population differentiation of the grass carp, genetics of its strong ability in adaption to diverse environments and sources of some introduced grass carp populations. Our data also suggests that the natural populations of the grass carp have been affected by the aquaculture besides neutral and adaptive forces.

Constructing High-Density Genetic Maps and Developing Sexing Markers in Northern Snakehead (Channa argus)

High-density genetic maps are essential for mapping QTL, improving genome assembly, comparative genomics, and studying sex chromosome evolution. The northern snakehead (Channa argus) is an economically important foodfish species with significant sexual dimorphism, where the males grow much faster and bigger than the females. However, to date, the sex determination pattern is still not clear, limiting identification of sex chromosomes, even sex determination genes and development of monosex populations that are valuable for both sex evolution of vertebrates and aquaculture practices. Here, a sex-averaged map and two sex-specific genetic maps were constructed with 2974, 2323, and 2338 SNPs, respectively. Little difference was observed in the pattern of sex-specific recombination between female- and male-specific genetic maps. Genome scan identified a major locus for sex determination at LG16. Females and males are, respectively, homogametic and heterogametic, suggesting an XY sex determination system for this species. By resequencing genomes, InDels in the sex-associated QTL region were discovered and used for developing sex-specific PCR assays for fast sexing of snakehead. These high-density genetic maps provide useful resources for future genomic studies in snakehead and its related species. The PCR assays for sexing are of importance in developing all male populations for aquaculture.

B Chromosomes of the Asian Seabass (Lates calcarifer) Contribute to Genome Variations at the Level of Individuals and Populations

The Asian seabass (Lates calcarifer) is a bony fish from the Latidae family, which is widely distributed in the tropical Indo-West Pacific region. The karyotype of the Asian seabass contains 24 pairs of A chromosomes and a variable number of AT- and GC-rich B chromosomes (Bchrs or Bs). Dot-like shaped and nucleolus-associated AT-rich Bs were microdissected and sequenced earlier. Here we analyzed DNA fragments from Bs to determine their repeat and gene contents using the Asian seabass genome as a reference. Fragments of 75 genes, including an 18S rRNA gene, were found in the Bs; repeats represented 2% of the Bchr assembly. The 18S rDNA of the standard genome and Bs were similar and enriched with fragments of transposable elements. A higher nuclei DNA content in the male gonad and somatic tissue, compared to the female gonad, was demonstrated by flow cytometry. This variation in DNA content could be associated with the intra-individual variation in the number of Bs. A comparison between the copy number variation among the B-related fragments from whole genome resequencing data of Asian seabass individuals identified similar profiles between those from the South-East Asian/Philippines and Indian region but not the Australian ones. Our results suggest that Bs might cause variations in the genome among the individuals and populations of Asian seabass. A personalized copy number approach for segmental duplication detection offers a suitable tool for population-level analysis across specimens with low coverage genome sequencing.

Population Genomics Reveals Genetic Divergence and Adaptive Differentiation of Chinese Sea Bass (Lateolabrax maculatus)

The marine species usually show high dispersal capabilities accompanied by high levels of gene flow. On the other hand, many physical barriers distribute along the continental marginal seas and may prevent dispersals and increase population divergence. These complexities along the continental margin generate serious challenges to population genetic studies of marine species. Chinese sea bass Lateolabrax maculatus distributes broad latitudinal gradient spanning from the tropical to the mid-temperate zones in the continental margin seas of the Northwest Pacific Ocean. Using the double digest restriction-site-associated DNA tag sequencing (ddRAD) approach, we genotyped 10,297 SNPs for 219 Chinese seabass individuals of 12 populations along the Chinese coast in the Northwest Pacific region. Genetic divergence among these populations was evaluated, and population structure was established. The results suggested that geographically distant populations in the Bohai Gulf and the Beibu Gulf retain significant genetic divergence, which are connected by a series of intermediate populations in between. The results also suggested that Leizhou Peninsula, Hainan Island, and Shandong Peninsula are major physical barriers and substantially block gene flow and genetic admixture of L. maculatus. We also investigated the potential genetic basis of local adaptation correlating with population differentiation of L. maculatus. The sea surface temperature is a significantly differentiated environmental factor for the distribution of L. maculatus. The correlation of water temperature and genetic variations in extensively distributed populations was investigated with Bayesian-based approaches. The candidate genes underlying the local selection in geographically divergent populations were identified and annotated, providing clues to understand the potential mechanisms of adaptive evolution. Overall, our genome scale population genetic analysis provided insight into population divergence and local adaptation of Chinese sea bass in the continental marginal seas along Chinese coast.

Survival and divergence in a small group: The extraordinary genomic history of the endangered Apennine brown bear stragglers

About 100 km east of Rome, in the central Apennine Mountains, a critically endangered population of ∼50 brown bears live in complete isolation. Mating outside this population is prevented by several 100 km of bear-free territories. We exploited this natural experiment to better understand the gene and genomic consequences of surviving at extremely small population size. We found that brown bear populations in Europe lost connectivity since Neolithic times, when farming communities expanded and forest burning was used for land clearance. In central Italy, this resulted in a 40-fold population decline. The overall genomic impact of this decline included the complete loss of variation in the mitochondrial genome and along long stretches of the nuclear genome. Several private and deleterious amino acid changes were fixed by random drift; predicted effects include energy deficit, muscle weakness, anomalies in cranial and skeletal development, and reduced aggressiveness. Despite this extreme loss of diversity, Apennine bear genomes show nonrandom peaks of high variation, possibly maintained by balancing selection, at genomic regions significantly enriched for genes associated with immune and olfactory systems. Challenging the paradigm of increased extinction risk in small populations, we suggest that random fixation of deleterious alleles (i) can be an important driver of divergence in isolation, (ii) can be tolerated when balancing selection prevents random loss of variation at important genes, and (iii) is followed by or results directly in favorable behavioral changes.

QTL Mapping for Resistance to Iridovirus in Asian Seabass Using Genotyping-by-Sequencing

Identifying quantitative trait loci (QTL) for viral disease resistance is of particular importance in selective breeding programs of fish species. Genetic markers linked to QTL can be useful in marker-assisted selection (MAS) for elites resistant to specific pathogens. Here, we conducted a genome scan for QTL associated with Singapore grouper iridovirus (SGIV) resistance in an Asian seabass (Lates calcarifer) family, using a high-density linkage map generated with genotyping-by-sequencing. One genome-wide significant and three suggestive QTL were detected at LG21, LG6, LG13, and LG15, respectively. The phenotypic variation explained (PVE) by the four QTL ranged from 7.5 to 15.6%. The position of the most significant QTL at LG21 was located between 31.88 and 36.81 cM. The SNP marker (SNP130416) nearest to the peak of this QTL was significantly associated with SGIV resistance in an unrelated multifamily population. One candidate gene, MECOM, close to the peak of this QTL region, was predicted. Evidence of alternative splicing was observed for MECOM and one specific category of splicing variants was differentially expressed at 5 days post-SGIV infection. The QTL detected in this study are valuable resources and can be used in the selective breeding programs of Asian seabass with regard to resistance to SGIV.

Genome-Wide Association Study Identifies Loci Associated with Resistance to Viral Nervous Necrosis Disease in Asian Seabass

Viral nervous necrosis disease (VNN), caused by nervous necrosis virus (NNV), is one major threat to mariculture. Identifying loci and understanding the mechanisms associated with resistance to VNN are important in selective breeding programs. We performed a genome-wide association study (GWAS) using genotyping-by-sequencing (GBS) to study the genomic architecture of resistance to NNV infection in Asian seabass. We genotyped 986 individuals from 43 families produced by 15 founders with 44498 bi-allelic genetic variants using GBS. The GWAS identified three genome-wide significant loci on chromosomes 16, 19, and 20, respectively, and six suggestive loci on chromosomes 1, 8, 14, 15, 21, and 24, respectively, associated with resistance to NNV infection measured as binary and quantitative traits. Using the 500 most significant markers in combination with a training population of 800 samples could reach a genomic prediction accuracy of 0.7. Candidate genes significantly associated with resistance to NNV, including lysine-specific demethylase 2A, beta-defensin 1, and cystatin-B, which play important roles in immune responses against virus infection, were identified. Almost all the candidate genes were differentially expressed in different tissues against NNV infection. The significant genetic variants can be used in genomic selection and help understand the mechanism of resistance to VNN. Future studies should use populations of large effective size and whole genome resequencing to identify more useful genetic variants.

Construction of high-resolution recombination maps in Asian seabass

Background

A high-density genetic map is essential for de novo genome assembly, fine mapping QTL for important complex traits, comparative genomic studies and understanding the mechanisms of genome evolution. Although a number of genomic resources are available in Asian seabass (Lates calcarifer), a high-density linkage map is still lacking. To facilitate QTL mapping for marker-assisted selection and genome assembly, and to understand the genome-wide recombination rates, we constructed high density linkage maps using three families and genotyping by sequencing.

Results

A high-density consensus linkage map consisting of 8, 274 markers was constructed based on sex-averaged genetic maps. The genetic maps were then aligned and integrated with the current genome assembly of Asian seabass. More than 90% of the genome contig sequences were anchored onto the consensus genetic map. Evidence of assembly errors in the current genome assembly was identified. A fragment of up to 2.5 Mb belonging to LG14 was assembled into Chr15. The length of family-specific sex-averaged maps ranged from 1348.96 to 1624.65 cM. Female maps were slightly longer than male maps using common markers. Female-to-male ratios were highly variable both across chromosomes within each family and throughout three families for each chromosome. However, the distribution patterns of recombination along chromosomes were similar between sexes across the whole genome. The overall recombination rates were significantly correlated with genome-wide GC content and the correlations were revealed to be stronger in females than in males.

Conclusions

These high-density genetic maps provide not only essential tools for facilitating de novo genome assembly and comparative genomic studies in teleosts, but also critical resources for fine mapping QTL and genome-wide association mapping for economically important traits in Asian seabass.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3462-z) contains supplementary material, which is available to authorized users.