Construction of Ultradense Linkage Maps with Lep-MAP2 : Stickleback F-2 Recombinant Crosses as an Example

High-density linkage maps are important tools for genome biology and evolutionary genetics by quantifying the extent of recombination, linkage disequilibrium, and chromosomal rearrangements across chromosomes, sexes, and populations. They provide one of the best ways to validate and refine de novo genome assemblies, with the power to identity errors in assemblies increasing with marker density. However, assembly of high-density linkage maps is still challenging due to software limitations. We describe Lep-MAP2, a software for ultradense genome-wide linkage map construction. Lep-MAP2 can handle various family structures and can account for achiasmatic meiosis to gain linkage map accuracy. Simulations show that Lep-MAP2 outperforms other available mapping software both in computational efficiency and accuracy. When applied to two large F-2-generation recombinant crosses between two nine-spined stickleback (Pungitius pungitius) populations, it produced two high-density (similar to 6 markers/cM) linkage maps containing 18,691 and 20,054 single nucleotide polymorphisms. The two maps showed a high degree of synteny, but female maps were 1.5-2 times longer than male maps in all linkage groups, suggesting genome-wide recombination suppression in males. Comparison with the genome sequence of the three-spined stickleback (Gasterosteus aculeatus) revealed a high degree of interspecific synteny with a low frequency (Peer reviewe

Contrasting population genetic responses to migration barriers in two native and an invasive freshwater fish

Habitat fragmentation impacts the distribution of genetic diversity and population genetic structure. Therefore, protecting the evolutionary potential of species, especially in the context of the current rate of human-induced environmental change, is an important goal. In riverine ecosystems, migration barriers affect the genetic structure of native species, while also influencing the spread of invasive species. In this study, we compare genetic patterns of two native and one highly invasive riverine fish species in a Belgian river basin, namely the native three-spined stickleback (Gasterosteus aculeatus) and stone loach (Barbatula barbatula), and the non-native and invasive topmouth gudgeon (Pseudorasbora parva). We aimed to characterize both natural and anthropogenic determinants of genetic diversity and population genetic connectivity. Genetic diversity was highest in topmouth gudgeon, followed by stone loach and three-spined stickleback. The correlation between downstream distance and genetic diversity, a pattern often observed in riverine systems, was only marginally significant in stone loach and three-spined stickleback, while genetic diversity strongly declined with increasing number of barriers in topmouth gudgeon. An Isolation-By-Distance pattern characterizes the population genetic structure of each species. Population differentiation was only associated with migration barriers in the invasive topmouth gudgeon, while genetic composition of all species seemed at least partially determined by the presence of migration barriers. Among the six barrier types considered (watermills, sluices, tunnels, weirs, riverbed obstructions, and others), the presence of watermills was the strongest driver of genetic structure and composition. Our results indicate that conservation and restoration actions, focusing on conserving genetic patterns, cannot be generalized across species. Moreover, measures might target either on restoring connectivity, while risking a rapid spread of the invasive topmouth gudgeon, or not restoring connectivity, while risking native species extinction in upstream populations

Reconciling seascape genetics and fisheries science in three codistributed flatfishes

Uncertainty hampers innovative mixed‐fisheries management by the scales at which connectivity dynamics are relevant to management objectives. The spatial scale of sustainable stock management is species‐specific and depends on ecology, life history and population connectivity. One valuable approach to understand these spatial scales is to determine to what extent population genetic structure correlates with the oceanographic environment. Here, we compare the level of genetic connectivity in three codistributed and commercially exploited demersal flatfish species living in the North East Atlantic Ocean. Population genetic structure was analysed based on 14, 14 and 10 neutral DNA microsatellite markers for turbot, brill and sole, respectively. We then used redundancy analysis (RDA) to attribute the genetic variation to spatial (geographical location), temporal (sampling year) and oceanographic (water column characteristics) components. The genetic structure of turbot was composed of three clusters and correlated with variation in the depth of the pycnocline, in addition to spatial factors. The genetic structure of brill was homogenous, but correlated with average annual stratification and spatial factors. In sole, the genetic structure was composed of three clusters, but was only linked to a temporal factor. We explored whether the management of data poor commercial fisheries, such as in brill and turbot, might benefit from population‐specific information. We conclude that the management of fish stocks has to consider species‐specific genetic structures and may benefit from the documentation of the genetic seascape and life‐history traits.publishedVersionUnit Licence Agreemen

Reconciling seascape genetics and fisheries science in three codistributed flatfishes

Uncertainty hampers innovative mixed‐fisheries management by the scales at which connectivity dynamics are relevant to management objectives. The spatial scale of sustainable stock management is species‐specific and depends on ecology, life history and population connectivity. One valuable approach to understand these spatial scales is to determine to what extent population genetic structure correlates with the oceanographic environment. Here, we compare the level of genetic connectivity in three codistributed and commercially exploited demersal flatfish species living in the North East Atlantic Ocean. Population genetic structure was analysed based on 14, 14 and 10 neutral DNA microsatellite markers for turbot, brill and sole, respectively. We then used redundancy analysis (RDA) to attribute the genetic variation to spatial (geographical location), temporal (sampling year) and oceanographic (water column characteristics) components. The genetic structure of turbot was composed of three clusters and correlated with variation in the depth of the pycnocline, in addition to spatial factors. The genetic structure of brill was homogenous, but correlated with average annual stratification and spatial factors. In sole, the genetic structure was composed of three clusters, but was only linked to a temporal factor. We explored whether the management of data poor commercial fisheries, such as in brill and turbot, might benefit from population‐specific information. We conclude that the management of fish stocks has to consider species‐specific genetic structures and may benefit from the documentation of the genetic seascape and life‐history traits.publishedVersionUnit Licence Agreemen

A High-Quality Assembly of the Nine-Spined Stickleback (Pungitius pungitius) Genome

The Gasterosteidae fish family hosts several species that are important models for eco-evolutionary, genetic, and genomic research. In particular, a wealth of genetic and genomic data has been generated for the three-spined stickleback (Gasterosteus aculeatus), the "ecology's supermodel," whereas the genomic resources for the nine-spined stickleback (Pungitius pungitius) have remained relatively scarce. Here, we report a high-quality chromosome-level genome assembly of P. pungitius consisting of 5,303 contigs (N50 = 1.2Mbp) with a total size of 521 Mbp. These contigs were mapped to 21 linkage groups using a high-density linkage map, yielding a final assembly with 98.5% BUSCO completeness. A total of 25,062 protein-coding genes were annotated, and about 23% of the assembly was found to consist of repetitive elements. A comprehensive analysis of repetitive elements uncovered centromere-specific tandem repeats and provided insights into the evolution of retrotransposons. A multigene phylogenetic analysis inferred a divergence time of about 26 million years ago (Ma) between nine- and three-spined sticklebacks, which is far older than the commonly assumed estimate of 13 Ma. Compared with the three-spined stickleback, we identified an additional duplication of several genes in the hemoglobin cluster. Sequencing data from populations adapted to different environments indicated potential copy number variations in hemoglobin genes. Furthermore, genome-wide synteny comparisons between three- and nine-spined sticklebacks identified chromosomal rearrangements underlying the karyotypic differences between the two species. The high-quality chromosome-scale assembly of the nine-spined stickleback genome obtained with long-read sequencing technology provides a crucial resource for comparative and population genomic investigations of stickleback fishes and teleosts.Peer reviewe

Predicting fish community responses to environmental policy targets

The European Union adopted the Water Framework Directive (WFD) in the year 2000 to tackle the rapid degradation of freshwater systems. However, biological, hydromorphological, and physico-chemical water quality targets are currently not met, and identifying successful policy implementation and management actions is of key importance. We built a joint species distribution model for riverine fish in Flanders (Belgium) to better understand the response of fish communities to current environmental policy goals. Environmental covariates included physico-chemical variables and hydromorphological quality indices, while waterway distances accounted for spatial effects. We detected strong effects of physico-chemistry on fish species' distributions. Evaluation of fish community responses to simulated policy scenarios revealed that targeting a 'good' status, following the WFD, increases average species richness with a fraction of species (0.13-0.69 change in accumulated occurrence probabilities). Targeting a 'very good' status, however, predicted an increase of 0.17-1.38 in average species richness. These simulations indicated that riverbed quality, nitrogen, and conductivity levels should be the focal point of policy. However, the weak response of species to a 'good' quality together with the complexity of nutrient-associated problems, suggest a challenging future for river restoration in Flanders.Peer reviewe

MultiPhen: Joint Model of Multiple Phenotypes Can Increase Discovery in GWAS

The genome-wide association study (GWAS) approach has discovered hundreds of genetic variants associated with diseases and quantitative traits. However, despite clinical overlap and statistical correlation between many phenotypes, GWAS are generally performed one-phenotype-at-a-time. Here we compare the performance of modelling multiple phenotypes jointly with that of the standard univariate approach. We introduce a new method and software, MultiPhen, that models multiple phenotypes simultaneously in a fast and interpretable way. By performing ordinal regression, MultiPhen tests the linear combination of phenotypes most associated with the genotypes at each SNP, and thus potentially captures effects hidden to single phenotype GWAS. We demonstrate via simulation that this approach provides a dramatic increase in power in many scenarios. There is a boost in power for variants that affect multiple phenotypes and for those that affect only one phenotype. While other multivariate methods have similar power gains, we describe several benefits of MultiPhen over these. In particular, we demonstrate that other multivariate methods that assume the genotypes are normally distributed, such as canonical correlation analysis (CCA) and MANOVA, can have highly inflated type-1 error rates when testing case-control or non-normal continuous phenotypes, while MultiPhen produces no such inflation. To test the performance of MultiPhen on real data we applied it to lipid traits in the Northern Finland Birth Cohort 1966 (NFBC1966). In these data MultiPhen discovers 21% more independent SNPs with known associations than the standard univariate GWAS approach, while applying MultiPhen in addition to the standard approach provides 37% increased discovery. The most associated linear combinations of the lipids estimated by MultiPhen at the leading SNPs accurately reflect the Friedewald Formula, suggesting that MultiPhen could be used to refine the definition of existing phenotypes or uncover novel heritable phenotypes

Mannose-Binding Lectin 2 Gene and Risk of Adult Glioma

Collection of specimens in this study, including funding: JMG JM MS SST DJH CAC. Selection of genes and SNPs to include in the analysis: FCFC MES. Manuscript editing: DMB DGC MS SST CAC ATP. Conceived and designed the experiments: DSM ATP. Performed the experiments: AS. Analyzed the data: DSM AS DGC. Wrote the paper: DSM.Background and AimsThe immune system is likely to play a key role in the etiology of gliomas. Genetic polymorphisms in the mannose-binding lectin gene, a key activator in the lectin complement pathway, have been associated with risk of several cancers.MethodsTo examine the role of the lectin complement pathway, we combined data from prospectively collected cohorts with available DNA specimens. Using a nested case-control design, we genotyped 85 single nucleotide polymorphisms (SNPs) in 9 genes in the lectin complement pathway and 3 additional SNPs in MBL2 were tested post hoc). Initial SNPs were selected using tagging SNPs for haplotypes; the second group of SNPs for MBL2 was selected based on functional SNPs related to phenotype. Associations were examined using logistic regression analysis. All statistical tests were two-sided. Nominal p-values are presented and are not corrected for multiple comparisons.ResultsA total of 143 glioma cases and 419 controls were available for this analysis. Statistically significant associations were observed for two SNPs in the mannose-binding lectin 2 (ML2) gene and risk of glioma (rs1982266 and rs1800450, test for trend p = 0.003 and p = 0.04, respectively, using the additive model). One of these SNPs, rs1800450, was associated with a 58% increase in glioma risk among those carrying one or two mutated alleles (odds ratio = 1.58, 95% confidence interval = 0.99–2.54), compared to those homozygous for the wild type allele.ConclusionsOverall, our findings suggest that MBL may play a role in the etiology of glioma. Future studies are needed to confirm these findings which may be due to chance, and if reproduced, to determine mechanisms that link glioma pathogenesis with the MBL complement pathway.Yeshttp://www.plosone.org/static/editorial#pee

A Genome-Wide Association Study of Neuroticism in a Population-Based Sample

Neuroticism is a moderately heritable personality trait considered to be a risk factor for developing major depression, anxiety disorders and dementia. We performed a genome-wide association study in 2,235 participants drawn from a population-based study of neuroticism, making this the largest association study for neuroticism to date. Neuroticism was measured by the Eysenck Personality Questionnaire. After Quality Control, we analysed 430,000 autosomal SNPs together with an additional 1.2 million SNPs imputed with high quality from the Hap Map CEU samples. We found a very small effect of population stratification, corrected using one principal component, and some cryptic kinship that required no correction. NKAIN2 showed suggestive evidence of association with neuroticism as a main effect (p<10−6) and GPC6 showed suggestive evidence for interaction with age (p≈10−7). We found support for one previously-reported association (PDE4D), but failed to replicate other recent reports. These results suggest common SNP variation does not strongly influence neuroticism. Our study was powered to detect almost all SNPs explaining at least 2% of heritability, and so our results effectively exclude the existence of loci having a major effect on neuroticism