Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Current genetic data do not improve the prediction of type 2 diabetes mellitus: the CoLaus study.

CONTEXT: Several genetic risk scores to identify asymptomatic subjects at high risk of developing type 2 diabetes mellitus (T2DM) have been proposed, but it is unclear whether they add extra information to risk scores based on clinical and biological data. OBJECTIVE: The objective of the study was to assess the extra clinical value of genetic risk scores in predicting the occurrence of T2DM. DESIGN: This was a prospective study, with a mean follow-up time of 5 yr. SETTING AND SUBJECTS: The study included 2824 nondiabetic participants (1548 women, 52 ± 10 yr). MAIN OUTCOME MEASURE: Six genetic risk scores for T2DM were tested. Four were derived from the literature and two were created combining all (n = 24) or shared (n = 9) single-nucleotide polymorphisms of the previous scores. A previously validated clinic + biological risk score for T2DM was used as reference. RESULTS: Two hundred seven participants (7.3%) developed T2DM during follow-up. On bivariate analysis, no differences were found for all but one genetic score between nondiabetic and diabetic participants. After adjusting for the validated clinic + biological risk score, none of the genetic scores improved discrimination, as assessed by changes in the area under the receiver-operating characteristic curve (range -0.4 to -0.1%), sensitivity (-2.9 to -1.0%), specificity (0.0-0.1%), and positive (-6.6 to +0.7%) and negative (-0.2 to 0.0%) predictive values. Similarly, no improvement in T2DM risk prediction was found: net reclassification index ranging from -5.3 to -1.6% and nonsignificant (P ≥ 0.49) integrated discrimination improvement. CONCLUSIONS: In this study, adding genetic information to a previously validated clinic + biological score does not seem to improve the prediction of T2DM.

Quantitative genetic modeling and inference in the presence of nonignorable missing data.

Natural selection is typically exerted at some specific life stages. If natural selection takes place before a trait can be measured, using conventional models can cause wrong inference about population parameters. When the missing data process relates to the trait of interest, a valid inference requires explicit modeling of the missing process. We propose a joint modeling approach, a shared parameter model, to account for nonrandom missing data. It consists of an animal model for the phenotypic data and a logistic model for the missing process, linked by the additive genetic effects. A Bayesian approach is taken and inference is made using integrated nested Laplace approximations. From a simulation study we find that wrongly assuming that missing data are missing at random can result in severely biased estimates of additive genetic variance. Using real data from a wild population of Swiss barn owls Tyto alba, our model indicates that the missing individuals would display large black spots; and we conclude that genes affecting this trait are already under selection before it is expressed. Our model is a tool to correctly estimate the magnitude of both natural selection and additive genetic variance.

Forecasting changes in population genetic structure of alpine plants in response to global warming.

Species range shifts in response to climate and land use change are commonly forecasted with species distribution models based on species occurrence or abundance data. Although appealing, these models ignore the genetic structure of species, and the fact that different populations might respond in different ways because of adaptation to their environment. Here, we introduced ancestry distribution models, that is, statistical models of the spatial distribution of ancestry proportions, for forecasting intra-specific changes based on genetic admixture instead of species occurrence data. Using multi-locus genotypes and extensive geographic coverage of distribution data across the European Alps, we applied this approach to 20 alpine plant species considering a global increase in temperature from 0.25 to 4 °C. We forecasted the magnitudes of displacement of contact zones between plant populations potentially adapted to warmer environments and other populations. While a global trend of movement in a north-east direction was predicted, the magnitude of displacement was species-specific. For a temperature increase of 2 °C, contact zones were predicted to move by 92 km on average (minimum of 5 km, maximum of 212 km) and by 188 km for an increase of 4 °C (minimum of 11 km, maximum of 393 km). Intra-specific turnover-measuring the extent of change in global population genetic structure-was generally found to be moderate for 2 °C of temperature warming. For 4 °C of warming, however, the models indicated substantial intra-specific turnover for ten species. These results illustrate that, in spite of unavoidable simplifications, ancestry distribution models open new perspectives to forecast population genetic changes within species and complement more traditional distribution-based approaches.

Low dispersing species : population genetic processes in space and time : the genus "Trochulus" (Gastropoda: Hygromiidae) as a case study

RÉSUMÉ Une espèce est rarement composée d'une population unique. Parce que les individus ont des capacités de dispersion limitées et que les paysages sont des mosaïques d'habitats, la plupart des espèces sont plutôt composées de sous-populations connectées par la migration. Cette variation spatiale influence directement la distribution de la variabilité génétique dans et entre les populations. Durant ce travail, nous avons abordé certains des processus populationnels qui ont joué un rôle supposé dans l'apparition de nouvelles espèces au sein du genre Trochulus. Plus précisément, nous avons tenté d'évaluer les impacts respectifs de l'isolement passé (facteurs historiques) et présent (facteurs locaux). Nous avons d'abord pu montrer que les faibles capacités de dispersion des escargots terrestres ont directement influencé leur histoire évolutive à toutes les échelles spatiales et temporelles. En réduisant l'effet homogénéisant de la migration, une faible dispersion maintient dans les populations les traces génétiques d'évènements passés. A l'échelle de la distribution globale de Trochulus villosus, ces traces ont permis de reconstruire une histoire faite d'isolements et d'expansions de populations. En combinant des données génétiques avec une modélisation de la niche climatique passée, il a été possible de proposer un scénario significativement meilleur que toutes les hypothèses alternatives que nous avons testées. A l'échelle locale par contre, l'héritage historique est difficile à distinguer de la dynamique actuelle. Ce fut le cas des lignées mitochondriales du complexe sericeus-hispidus : les deux principales lignées étaient phylogénétiquement éloignées, avaient eu des démographies passées différentes et corrélaient avec des différences morphologiques. D'un autre côté, le flux de gène nucléaire était fort, contredisant l'idée de deux espèces cryptiques isolées reproductivement. Pour pouvoir conclure à la présence ou non de deux espèces, il nous a manqué des informations locales sur la dynamique des populations et les conditions écologiques que l'on trouve dans la région d'étude. Enfin, nous avons pu souligner que la connectivité entre populations d'escargots est soumise à la qualité des habitats et à leur organisation spatiale. Les escargots sont dépendants d'un habitat et s'y adaptent, comme l'indiquent la présence de «poils » uniquement sur la coquille d'espèces vivant dans des habitats humides ou la corrélation entre morphologie et habitat au sein du complexe sericeus-hispidus. Logiquement donc, les escargots migrent préférentiellement au travers d'habitats favorables comme l'a montré la réduction de flux de gènes au travers des prairies chez T. villosus (une espèce forestière). De ces données, nous pouvons supposer que les populations d'escargots en particulier, et des espèces à faible dispersion en général, ont de fortes chances d'être affectées par les changements climatiques, avec de probables implications pour leurs histoires évolutives. SUMMARY : Species rarely consists in a single population. Because individuals have limited dispersal abilities, because landscapes are habitat patchworks, most species are made of several subpopulations connected by migration. This spatial variation has consequences on the distribution of genetic diversity within and between populations, creating a structure among the populations. During the present work, we investigated some of the population processes assumed to have played an important role on the speciation within the genus Trochulus. More specifically, we questioned the respective impacts of past (historical factors) or present (local factors) population isolations. We first could show that the poor dispersal abilities of land snails have had profound impacts on their evolutionary histories at all spatial and temporal scales. Low dispersal maintains a strong signature of past events in the populations by minimising the homogenising effects of geneflow. At the scale of Trochulus villosus global distribution, they allowed to retrieve the detailed history of this species population isolations and expansions. Combining a large genetic dataset with paleo-climatic niche modelling ended up with a historical scenario significantly better than all traditional alternatives we tested. At local scale on the contrary, past events become difficult to tease apart from ongoing processes. This was the case for the divergent mitochondria) lineages within the sericeus-hispidus complex: the two principal lineages appeared to be phylogenetically distant, to have experienced different demographic histories and to correlate with morphological differences. On the other hand, nuclear (present day) geneflow was high, contradicting the idea of two reproductively isolated cryptic species. Information on the local population dynamics and environmental conditions are lacking to be able to decide whether past isolation has indeed resulted here in new species. Finally, we emphasised the importance of the habitat types present in a landscape as well as their spatial organisation for the population connectivity of land snails. These species are tightly dependent on a habitat and adapt to it as shown by thé occurrence of hair-like structures only in species living in humid environments or by the correlation between shell morphology and habitat in the sericeus-hispidus complex. As a result, land snails preferentially migrate through favourable habitats: Trochulus villosus, a forest species, had its geneflow significantly reduced across meadows. From these data, we can hypothesise that the populations of land snails in particular and of low dispersing species in general are likely to be strongly affected by the ongoing climate changes, with potential major consequences on their evolutionary histories.

Population pharmacokinetic study of memantine: effects of clinical and genetic factors.

BACKGROUND AND OBJECTIVE: Memantine, a frequently prescribed anti-dementia drug, is mainly eliminated unchanged by the kidneys, partly via tubular secretion. Considerable inter-individual variability in plasma concentrations has been reported. We aimed to investigate clinical and genetic factors influencing memantine disposition. METHODS: A population pharmacokinetic study was performed including data from 108 patients recruited in a naturalistic setting. Patients were genotyped for common polymorphisms in renal cation transporters (SLC22A1/2/5, SLC47A1, ABCB1) and nuclear receptors (NR1I2, NR1I3, RXR, PPAR) involved in transporter expression. RESULTS: The average clearance was 5.2 L/h with a 27 % inter-individual variability (percentage coefficient of variation). Glomerular filtration rate (p = 0.007) and sex (p = 0.001) markedly influenced memantine clearance. NR1I2 rs1523130 was identified as the unique significant genetic covariate for memantine clearance (p = 0.006), with carriers of the NR1I2 rs1523130 CT/TT genotypes presenting a 16 % slower memantine elimination than carriers of the CC genotype. CONCLUSION: The better understanding of inter-individual variability of memantine disposition might be beneficial in the context of individual dose optimization.

Genetic association tests: a method for the joint analysis of family and case-control data.

With the trend in molecular epidemiology towards both genome-wide association studies and complex modelling, the need for large sample sizes to detect small effects and to allow for the estimation of many parameters within a model continues to increase. Unfortunately, most methods of association analysis have been restricted to either a family-based or a case-control design, resulting in the lack of synthesis of data from multiple studies. Transmission disequilibrium-type methods for detecting linkage disequilibrium from family data were developed as an effective way of preventing the detection of association due to population stratification. Because these methods condition on parental genotype, however, they have precluded the joint analysis of family and case-control data, although methods for case-control data may not protect against population stratification and do not allow for familial correlations. We present here an extension of a family-based association analysis method for continuous traits that will simultaneously test for, and if necessary control for, population stratification. We further extend this method to analyse binary traits (and therefore family and case-control data together) and accurately to estimate genetic effects in the population, even when using an ascertained family sample. Finally, we present the power of this binary extension for both family-only and joint family and case-control data, and demonstrate the accuracy of the association parameter and variance components in an ascertained family sample.

Spatio-temporal population genetic structure of the parasitic mite Spinturnix bechsteini is shaped by its own demography and the social system of its bat host.

Information about the population genetic structures of parasites is important for an understanding of parasite transmission pathways and ultimately the co-evolution with their hosts. If parasites cannot disperse independently of their hosts, a parasite's population structure will depend upon the host's spatial distribution. Geographical barriers affecting host dispersal can therefore lead to structured parasite populations. However, how the host's social system affects the genetic structure of parasite populations is largely unknown. We used mitochondrial DNA (mtDNA) to describe the spatio-temporal population structure of a contact-transmitted parasitic wing mite (Spinturnix bechsteini) and compared it to that of its social host, the Bechstein's bat (Myotis bechsteinii). We observed no genetic differentiation between mites living on different bats within a colony. This suggests that mites can move freely among bats of the same colony. As expected in case of restricted inter-colony dispersal, we observed a strong genetic differentiation of mites among demographically isolated bat colonies. In contrast, we found a strong genetic turnover between years when we investigated the temporal variation of mite haplotypes within colonies. This can be explained with mite dispersal occuring between colonies and bottlenecks of mite populations within colonies. The observed absence of isolation by distance could be the result from genetic drift and/or from mites dispersing even between remote bat colonies, whose members may meet at mating sites in autumn or in hibernacula in winter. Our data show that the population structure of this parasitic wing mite is influenced by its own demography and the peculiar social system of its bat host.

Genetic bottlenecks driven by population disconnection.

Connectivity among populations plays a crucial role in maintaining genetic variation at a local scale, especially in small populations affected strongly by genetic drift. The negative consequences of population disconnection on allelic richness and gene diversity (heterozygosity) are well recognized and empirically established. It is not well recognized, however, that a sudden drop in local effective population size induced by such disconnection produces a temporary disequilibrium in allelic frequency distributions that is akin to the genetic signature of a demographic bottleneck. To document this effect, we used individual-based simulations and empirical data on allelic richness and gene diversity in six pairs of isolated versus well-connected (core) populations of European tree frogs. In our simulations, population disconnection depressed allelic richness more than heterozygosity and thus resulted in a temporary excess in gene diversity relative to mutation drift equilibrium (i.e., signature of a genetic bottleneck). We observed a similar excess in gene diversity in isolated populations of tree frogs. Our results show that population disconnection can create a genetic bottleneck in the absence of demographic collapse.

Decades of population genetic research reveal the need for harmonization of molecular markers: the grey wolf Canis lupus as a case study

Following protection measures implemented since the 1970s, large carnivores are currently increasing in number and returning to areas from which they were absent for decades or even centuries. Monitoring programmes for these species rely extensively on non-invasive sampling and genotyping. However, attempts to connect results of such studies at larger spatial or temporal scales often suffer from the incompatibility of genetic markers implemented by researchers in different laboratories. This is particularly critical for long-distance dispersers, revealing the need for harmonized monitoring schemes that would enable the understanding of gene flow and dispersal dynamics. Based on a review of genetic studies on grey wolves Canis lupus from Europe, we provide an overview of the genetic markers currently in use, and identify opportunities and hurdles for studies based on continent-scale datasets. Our results highlight an urgent need for harmonization of methods to enable transnational research based on data that have already been collected, and to allow these data to be linked to material collected in the future. We suggest timely standardization of newly developed genotyping approaches, and propose that action is directed towards the establishment of shared single nucleotide polymorphism panels, next-generation sequencing of microsatellites, a common reference sample collection and an online database for data exchange. Enhanced cooperation among genetic researchers dealing with large carnivores in consortia would facilitate streamlining of methods, their faster and wider adoption, and production of results at the large spatial scales that ultimately matter for the conservation of these charismatic species.

Association between parental history and genetic risk scores for coronary heart disease prediction: The population-based CoLaus study.

BACKGROUND AND AIMS: Parental history (PH) and genetic risk scores (GRSs) are separately associated with coronary heart disease (CHD), but evidence regarding their combined effects is lacking. We aimed to evaluate the joint associations and predictive ability of PH and GRSs for incident CHD. METHODS: Data for 4283 Caucasians were obtained from the population-based CoLaus Study, over median follow-up time of 5.6 years. CHD was defined as incident myocardial infarction, angina, percutaneous coronary revascularization or bypass grafting. Single nucleotide polymorphisms for CHD identified by genome-wide association studies were used to construct unweighted and weighted versions of three GRSs, comprising of 38, 53 and 153 SNPs respectively. RESULTS: PH was associated with higher values of all weighted GRSs. After adjustment for age, sex, smoking, diabetes, systolic blood pressure, low and high density lipoprotein cholesterol, PH was significantly associated with CHD [HR 2.61, 95% CI (1.47-4.66)] and further adjustment for GRSs did not change this estimate. Similarly, one standard deviation change of the weighted 153-SNPs GRS was significantly associated with CHD [HR 1.50, 95% CI (1.26-1.80)] and remained so, after further adjustment for PH. The weighted, 153-SNPs GRS, but not PH, modestly improved discrimination [(C-index improvement, 0.016), p = 0.048] and reclassification [(NRI improvement, 8.6%), p = 0.027] beyond cardiovascular risk factors. After including both the GRS and PH, model performance improved further [(C-index improvement, 0.022), p = 0.006]. CONCLUSION: After adjustment for cardiovascular risk factors, PH and a weighted, polygenic GRS were jointly associated with CHD and provided additive information for coronary events prediction.

Ecological change predicts population dynamics and genetic diversity over 120 000 years.

While ecological effects on short-term population dynamics are well understood, their effects over millennia are difficult to demonstrate and convincing evidence is scant. Using coalescent methods, we analysed past population dynamics of three lizard species (Psammodromus hispanicus, P. edwardsianus, P. occidentalis) and linked the results with climate change data covering the same temporal horizon (120 000 years). An increase in population size over time was observed in two species, and in P. occidentalis, no change was observed. Temporal changes in temperature seasonality and the maximum temperature of the warmest month were congruent with changes in population dynamics observed for the three species and both variables affected population density, either directly or indirectly (via a life-history trait). These results constitute the first solid link between ecological change and long-term population dynamics. The results moreover suggest that ecological change leaves genetic signatures that can be retrospectively traced, providing evidence that ecological change is a crucial driver of genetic diversity and speciation.

Risk prediction of prevalent diabetes in a Swiss population using a weighted genetic score--the CoLaus Study.

AIMS/HYPOTHESIS: Several susceptibility genes for type 2 diabetes have been discovered recently. Individually, these genes increase the disease risk only minimally. The goals of the present study were to determine, at the population level, the risk of diabetes in individuals who carry risk alleles within several susceptibility genes for the disease and the added value of this genetic information over the clinical predictors. METHODS: We constructed an additive genetic score using the most replicated single-nucleotide polymorphisms (SNPs) within 15 type 2 diabetes-susceptibility genes, weighting each SNP with its reported effect. We tested this score in the extensively phenotyped population-based cross-sectional CoLaus Study in Lausanne, Switzerland (n = 5,360), involving 356 diabetic individuals. RESULTS: The clinical predictors of prevalent diabetes were age, BMI, family history of diabetes, WHR, and triacylglycerol/HDL-cholesterol ratio. After adjustment for these variables, the risk of diabetes was 2.7 (95% CI 1.8-4.0, p = 0.000006) for individuals with a genetic score within the top quintile, compared with the bottom quintile. Adding the genetic score to the clinical covariates improved the area under the receiver operating characteristic curve slightly (from 0.86 to 0.87), yet significantly (p = 0.002). BMI was similar in these two extreme quintiles. CONCLUSIONS/INTERPRETATION: In this population, a simple weighted 15 SNP-based genetic score provides additional information over clinical predictors of prevalent diabetes. At this stage, however, the clinical benefit of this genetic information is limited.

Rare chromosome abnormalities, prevalence and prenatal diagnosis rates from population-based congenital anomaly registers in Europe.

The aim of this study is to quantify the prevalence and types of rare chromosome abnormalities (RCAs) in Europe for 2000-2006 inclusive, and to describe prenatal diagnosis rates and pregnancy outcome. Data held by the European Surveillance of Congenital Anomalies database were analysed on all the cases from 16 population-based registries in 11 European countries diagnosed prenatally or before 1 year of age, and delivered between 2000 and 2006. Cases were all unbalanced chromosome abnormalities and included live births, fetal deaths from 20 weeks gestation and terminations of pregnancy for fetal anomaly. There were 10,323 cases with a chromosome abnormality, giving a total birth prevalence rate of 43.8/10,000 births. Of these, 7335 cases had trisomy 21,18 or 13, giving individual prevalence rates of 23.0, 5.9 and 2.3/10,000 births, respectively (53, 13 and 5% of all reported chromosome errors, respectively). In all, 473 cases (5%) had a sex chromosome trisomy, and 778 (8%) had 45,X, giving prevalence rates of 2.0 and 3.3/10,000 births, respectively. There were 1,737 RCA cases (17%), giving a prevalence of 7.4/10,000 births. These included triploidy, other trisomies, marker chromosomes, unbalanced translocations, deletions and duplications. There was a wide variation between the registers in both the overall prenatal diagnosis rate of RCA, an average of 65% (range 5-92%) and the prevalence of RCA (range 2.4-12.9/10,000 births). In all, 49% were liveborn. The data provide the prevalence of families currently requiring specialised genetic counselling services in the perinatal period for these conditions and, for some, long-term care.

Efficacy and safety of certolizumab pegol in an unselected crohn's disease population: 26-week data of the FACTS II survey.

BACKGROUND: Certolizumab pegol (Cimzia, CZP) was approved for the treatment of Crohn's disease (CD) patients in 2007 in Switzerland as the first country worldwide. This prospective phase IV study aimed to evaluate the efficacy and safety of CZP over 26 weeks in a multicenter cohort of practice-based patients. METHODS: Evaluation questionnaires at baseline, week 6, and week 26 were completed by gastroenterologists in hospitals and private practices. Adverse events were evaluated according to World Health Organization (WHO) guidelines. RESULTS: Sixty patients (38F/22M) were included; 53% had complicated disease (stricturing or penetrating), 45% had undergone prior CD-related surgery. All patients had prior exposure to systemic steroids, 96% to immunomodulators, 73% to infliximab, and 43% to adalimumab. A significant decrease of the Harvey-Bradshaw Index (HBI) was observed under CZP therapy (12.2 ± 4.9 at week 0 versus 6.3 ± 4.7 at week 6 and 6.7 ± 5.3 at week 26, both P < 0.001). Response and remission rates were 70% and 40% (week 6) and 67% and 36%, respectively (week 26). The complete perianal fistula closure rate was 36% at week 6 and 55% at week 26. The frequency of adverse drug reactions attributed to CZP was 5%. CZP was continued in 88% of patients beyond week 6 and in 67% beyond week 26. CONCLUSIONS: In a population of CD patients with predominantly complicated disease behavior, CZP proved to be effective in induction and maintenance of response and remission. This series provides the first evidence of CZP's effectiveness in perianal fistulizing CD in clinical practice.