Genetic differentiation and zoogeography of Asian Suncus murinus (Mammalia: Soricidae)

Protein electrophoresis was used to assess the phylogenetic relationships of populations of the phenotypically variable Asian house shrew Suncus murinus. These populations represent a sample of both commensal and wild forms. They were compared to another taxon, S. montanus, which was formerly considered conspecific with S. murinus. Suncus dayi was used as an outgroup in all phylogenetic reconstructions. Within the S. murinus lineage, the allozyme data show very low levels of genetic differentiation among both wild and commensal Southeast Asian and Japanese samples when compared to the Indian populations. This pattern is consistent with the classical hypothesis of a recent introduction by man in Eastern Asia. The higher genetic diversity found within S. murinus from India, as well as previous mitochondrial and karyological results suggest that this area is the probable centre of origin for the species. Although the lack of gene flow between S. murinus and S. montanus is clearly established in an area of sympatry in Southern India, one Asian house shrew sampled in Nepal was more closely related to S. montanus. This could either reflect the retention of an ancestral polymorphism, or result from a hybridization episode between S. murinus and S. montanus. Similar conclusions were also suggested in mitochondrial DNA studies dealing with animals sampled in the Northern parts of the Indian subcontinent. Clearly, further data on Suncus from this area are needed in order to assess these hypotheses. (C) 1995 The Linnean Society of London

GWAS of 126,559 individuals identifies genetic variants associated with educational attainment.

A genome-wide association study (GWAS) of educational attainment was conducted in a discovery sample of 101,069 individuals and a replication sample of 25,490. Three independent single-nucleotide polymorphisms (SNPs) are genome-wide significant (rs9320913, rs11584700, rs4851266), and all three replicate. Estimated effects sizes are small (coefficient of determination R(2) ≈ 0.02%), approximately 1 month of schooling per allele. A linear polygenic score from all measured SNPs accounts for ≈2% of the variance in both educational attainment and cognitive function. Genes in the region of the loci have previously been associated with health, cognitive, and central nervous system phenotypes, and bioinformatics analyses suggest the involvement of the anterior caudate nucleus. These findings provide promising candidate SNPs for follow-up work, and our effect size estimates can anchor power analyses in social-science genetics.

Maintenance of polymorphism : a theoretical approach using multilocus systems

Résumé La thématique de cette thèse peut être résumée par le célèbre paradoxe de biologie évolutive sur le maintien du polymorphisme face à la sélection et par l'équation du changement de fréquence gamétique au cours du temps dû, à la sélection. La fréquence d'un gamète xi à la génération (t + 1) est: !!!Equation tronquée!!! Cette équation est utilisée pour générer des données utlisée tout au long de ce travail pour 2, 3 et 4 locus dialléliques. Le potentiel de l'avantage de l'hétérozygote pour le maintien du polymorphisme est le sujet de la première partie. La définition commune de l'avantage de l'hétérozygote n'etant applicable qu'a un locus ayant 2 allèles, cet avantage est redéfini pour un système multilocus sur les bases de précédentes études. En utilisant 5 définitions différentes de l'avantage de l'hétérozygote, je montre que cet avantage ne peut être un mécanisme général dans le maintien du polymorphisme sous sélection. L'étude de l'influence de locus non-détectés sur les processus évolutifs, seconde partie de cette thèse, est motivée par les travaux moléculaires ayant pour but de découvrir le nombre de locus codant pour un trait. La plupart de ces études sous-estiment le nombre de locus. Je montre que des locus non-détectés augmentent la probabilité d'observer du polymorphisme sous sélection. De plus, les conclusions sur les facteurs de maintien du polymorphisme peuvent être trompeuses si tous les locus ne sont pas détectés. Dans la troisième partie, je m'intéresse à la valeur attendue de variance additive après un goulot d'étranglement pour des traits sélectionés. Une études précédente montre que le niveau de variance additive après goulot d'étranglement augmente avec le nombre de loci. Je montre que le niveau de variance additive après un goulot d'étranglement augmente (comparé à des traits neutres), mais indépendamment du nombre de loci. Par contre, le taux de recombinaison a une forte influence, entre autre en regénérant les gamètes disparus suite au goulot d'étranglement. La dernière partie de ce travail de thèse décrit un programme pour le logiciel de statistique R. Ce programme permet d'itérer l'équation ci-dessus en variant les paramètres de sélection, recombinaison et de taille de populations pour 2, 3 et 4 locus dialléliques. Cette thèse montre qu'utiliser un système multilocus permet d'obtenir des résultats non-conformes à ceux issus de systèmes rnonolocus (la référence en génétique des populations). Ce programme ouvre donc d'intéressantes perspectives en génétique des populations. Abstract The subject of this PhD thesis can be summarized by one famous paradox of evolu-tionary biology: the maintenance of polymorphism in the face of selection, and one classical equation of theoretical population genetics: the changes in gametic frequencies due to selection and recombination. The frequency of gamete xi at generation (t + 1) is given by: !!! Truncated equation!!! This equation is used to generate data on selection at two, three, and four diallelic loci for the different parts of this work. The first part focuses on the potential of heterozygote advantage to maintain genetic polymorphism. Results of previous studies are used to (re)define heterozygote advantage for multilocus systems, since the classical definition is for one diallelic locus. I use 5 different definitions of heterozygote advantage. And for these five definitions, I show that heterozygote advantage is not a general mechanism for the maintenance of polymorphism. The study of the influence of undetected loci on evolutionary processes (second part of this work) is motivated by molecular works which aim at discovering the loci coding for a trait. For most of these works, some coding loci remains undetected. I show that undetected loci increases the probability of maintaining polymorphism under selection. In addition, conclusions about the factor that maintain polymorphism can be misleading if not all loci are considered. This is, therefore, only when all loci are detected that exact conclusions on the level of maintained polymorphism or on the factor(s) that maintain(s) polymorphism could be drawn. In the third part, the focus is on the expected release of additive genetic variance after bottleneck for selected traits. A previous study shows that the expected release of additive variance increases with an increase in the number of loci. I show that the expected release of additive variance after bottleneck increases for selected traits (compared with neutral), but this increase is not a function of the number of loci, but function of the recombination rate. Finally, the last part of this PhD thesis is a description of a package for the statistical software R that implements the Equation given above. It allows to generate data for different scenario regarding selection, recombination, and population size. This package opens perspectives for the theoretical population genetics that mainly focuses on one locus, while this work shows that increasing the number of loci leads not necessarily to straightforward results.

Genetic variation in GIPR influences the glucose and insulin responses to an oral glucose challenge.

Glucose levels 2 h after an oral glucose challenge are a clinical measure of glucose tolerance used in the diagnosis of type 2 diabetes. We report a meta-analysis of nine genome-wide association studies (n = 15,234 nondiabetic individuals) and a follow-up of 29 independent loci (n = 6,958-30,620). We identify variants at the GIPR locus associated with 2-h glucose level (rs10423928, beta (s.e.m.) = 0.09 (0.01) mmol/l per A allele, P = 2.0 x 10(-15)). The GIPR A-allele carriers also showed decreased insulin secretion (n = 22,492; insulinogenic index, P = 1.0 x 10(-17); ratio of insulin to glucose area under the curve, P = 1.3 x 10(-16)) and diminished incretin effect (n = 804; P = 4.3 x 10(-4)). We also identified variants at ADCY5 (rs2877716, P = 4.2 x 10(-16)), VPS13C (rs17271305, P = 4.1 x 10(-8)), GCKR (rs1260326, P = 7.1 x 10(-11)) and TCF7L2 (rs7903146, P = 4.2 x 10(-10)) associated with 2-h glucose. Of the three newly implicated loci (GIPR, ADCY5 and VPS13C), only ADCY5 was found to be associated with type 2 diabetes in collaborating studies (n = 35,869 cases, 89,798 controls, OR = 1.12, 95% CI 1.09-1.15, P = 4.8 x 10(-18)).

Genetic variants influencing circulating lipid levels and risk of coronary artery disease.

OBJECTIVE: Genetic studies might provide new insights into the biological mechanisms underlying lipid metabolism and risk of CAD. We therefore conducted a genome-wide association study to identify novel genetic determinants of low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and triglycerides. METHODS AND RESULTS: We combined genome-wide association data from 8 studies, comprising up to 17 723 participants with information on circulating lipid concentrations. We did independent replication studies in up to 37 774 participants from 8 populations and also in a population of Indian Asian descent. We also assessed the association between single-nucleotide polymorphisms (SNPs) at lipid loci and risk of CAD in up to 9 633 cases and 38 684 controls. We identified 4 novel genetic loci that showed reproducible associations with lipids (probability values, 1.6×10(-8) to 3.1×10(-10)). These include a potentially functional SNP in the SLC39A8 gene for HDL-C, an SNP near the MYLIP/GMPR and PPP1R3B genes for LDL-C, and at the AFF1 gene for triglycerides. SNPs showing strong statistical association with 1 or more lipid traits at the CELSR2, APOB, APOE-C1-C4-C2 cluster, LPL, ZNF259-APOA5-A4-C3-A1 cluster and TRIB1 loci were also associated with CAD risk (probability values, 1.1×10(-3) to 1.2×10(-9)). CONCLUSIONS: We have identified 4 novel loci associated with circulating lipids. We also show that in addition to those that are largely associated with LDL-C, genetic loci mainly associated with circulating triglycerides and HDL-C are also associated with risk of CAD. These findings potentially provide new insights into the biological mechanisms underlying lipid metabolism and CAD risk.

Variable queen number in ant colonies: no impact on queen turnover, inbreeding, and population genetic differentiation in the ant Formica selysi.

Variation in queen number alters the genetic structure of social insect colonies, which in turn affects patterns of kin-selected conflict and cooperation. Theory suggests that shifts from single- to multiple-queen colonies are often associated with other changes in the breeding system, such as higher queen turnover, more local mating, and restricted dispersal. These changes may restrict gene flow between the two types of colonies and it has been suggested that this might ultimately lead to sympatric speciation. We performed a detailed microsatellite analysis of a large population of the ant Formica selysi, which revealed extensive variation in social structure, with 71 colonies headed by a single queen and 41 by multiple queens. This polymorphism in social structure appeared stable over time, since little change in the number of queens per colony was detected over a five-year period. Apart from queen number, single- and multiple-queen colonies had very similar breeding systems. Queen turnover was absent or very low in both types of colonies. Single- and multiple-queen colonies exhibited very small but significant levels of inbreeding, which indicates a slight deviation from random mating at a local scale and suggests that a small proportion of queens mate with related males. For both types of colonies, there was very little genetic structuring above the level of the nest, with no sign of isolation by distance. These similarities in the breeding systems were associated with a complete lack of genetic differentiation between single- and multiple-queen colonies, which provides no support for the hypothesis that change in queen number leads to restricted gene flow between social forms. Overall, this study suggests that the higher rates of queen turnover, local mating, and population structuring that are often associated with multiple-queen colonies do not appear when single- and multiple-queen colonies still coexist within the same population, but build up over time in populations consisting mostly of multiple-queen colonies.

APOBEC3G genetic variants and their influence on the progression to AIDS.

The cytosine deaminase APOBEC3G, in the absence of the human immunodeficiency virus type 1 (HIV-1) accessory gene HIV-1 viral infectivity factor (vif), inhibits viral replication by introducing G-->A hypermutation in the newly synthesized HIV-1 DNA negative strand. We tested the hypothesis that genetic variants of APOBEC3G may modify HIV-1 transmission and disease progression. Single nucleotide polymorphisms were identified in the promoter region (three), introns (two), and exons (two). Genotypes were determined for 3,073 study participants enrolled in six HIV-AIDS prospective cohorts. One codon-changing variant, H186R in exon 4, was polymorphic in African Americans (AA) (f = 37%) and rare in European Americans (f < 3%) or Europeans (f = 5%). For AA, the variant allele 186R was strongly associated with decline in CD4 T cells (CD4 slope on square root scale: -1.86, P = 0.009), The 186R allele was also associated with accelerated progression to AIDS-defining conditions in AA. The in vitro antiviral activity of the 186R enzyme was not inferior to that of the common H186 variant. These studies suggest that there may be a modifying role of variants of APOBEC3G on HIV-1 disease progression that warrants further investigation.

Developmental critical period in genetic redox dysregulation: animal and human studies in schizophrenia

Converging evidence favors an abnormal susceptibility to oxidative stress in schizophrenia. Decreased levels of glutathione (GSH), the major cellular antioxidant and redox regulator, was observed in cerebrospinal-fluid and prefrontal cortex of patients. Importantly, abnormal GSH synthesis of genetic origin was observed: Two case-control studies showed an association with a GAG trinucleotide repeat (TNR) polymorphism in the GSH key synthesizing enzyme glutamate-cysteine-ligase (GCL) catalytic subunit (GCLC) gene. The most common TNR genotype 7/7 was more frequent in controls, whereas the rarest TNR genotype 8/8 was three times more frequent in patients. The disease associated genotypes (35% of patients) correlated with decreased GCLC protein, GCL activity and GSH content. Similar GSH system anomalies were observed in early psychosis patients. Such redox dysregulation combined with environmental stressors at specific developmental stages could underlie structural and functional connectivity anomalies. In pharmacological and knock-out (KO) models, GSH deficit induces anomalies analogous to those reported in patients. (a) morphology: spine density and GABA-parvalbumine immunoreactivity (PV-I) were decreased in anterior cingulate cortex. KO mice showed delayed cortical PV-I at PD10. This effect is exacerbated in mice with increased DA from PD5-10. KO mice exhibit cortical impairment in myelin and perineuronal net known to modulate PV connectivity. (b) physiology: In cultured neurons, NMDA response are depressed by D2 activation. In hippocampus, NMDA-dependent synaptic plasticity is impaired and kainate induced g-oscillations are reduced in parallel to PV-I. (c) cognition: low GSH models show increased sensitivity to stress, hyperactivity, abnormal object recognition, olfactory integration and social behavior. In a clinical study, GSH precursor N-acetyl cysteine (NAC) as add on therapy, improves the negative symptoms and decreases the side effects of antipsychotics. In an auditory oddball paradigm, NAC improves the mismatched negativity, an evoked potential related to pre-attention and to NMDA receptors function. In summary, clinical and experimental evidence converge to demonstrate that a genetically induced dysregulation of GSH synthesis combined with environmental insults in early development represent a major risk factor contributing to the development of schizophrenia

FTO genetic variants, dietary intake and body mass index: insights from 177,330 individuals.

FTO is the strongest known genetic susceptibility locus for obesity. Experimental studies in animals suggest the potential roles of FTO in regulating food intake. The interactive relation among FTO variants, dietary intake and body mass index (BMI) is complex and results from previous often small-scale studies in humans are highly inconsistent. We performed large-scale analyses based on data from 177,330 adults (154 439 Whites, 5776 African Americans and 17 115 Asians) from 40 studies to examine: (i) the association between the FTO-rs9939609 variant (or a proxy single-nucleotide polymorphism) and total energy and macronutrient intake and (ii) the interaction between the FTO variant and dietary intake on BMI. The minor allele (A-allele) of the FTO-rs9939609 variant was associated with higher BMI in Whites (effect per allele = 0.34 [0.31, 0.37] kg/m(2), P = 1.9 × 10(-105)), and all participants (0.30 [0.30, 0.35] kg/m(2), P = 3.6 × 10(-107)). The BMI-increasing allele of the FTO variant showed a significant association with higher dietary protein intake (effect per allele = 0.08 [0.06, 0.10] %, P = 2.4 × 10(-16)), and relative weak associations with lower total energy intake (-6.4 [-10.1, -2.6] kcal/day, P = 0.001) and lower dietary carbohydrate intake (-0.07 [-0.11, -0.02] %, P = 0.004). The associations with protein (P = 7.5 × 10(-9)) and total energy (P = 0.002) were attenuated but remained significant after adjustment for BMI. We did not find significant interactions between the FTO variant and dietary intake of total energy, protein, carbohydrate or fat on BMI. Our findings suggest a positive association between the BMI-increasing allele of FTO variant and higher dietary protein intake and offer insight into potential link between FTO, dietary protein intake and adiposity.

Sex-linked inheritance, genetic correlations and sexual dimorphism in three melanin-based colour traits in the barn owl.

Theory states that genes on the sex chromosomes have stronger effects on sexual dimorphism than genes on the autosomes. Although empirical data are not necessarily consistent with this theory, this situation may prevail because the relative role of sex-linked and autosomally inherited genes on sexual dimorphism has rarely been evaluated. We estimated the quantitative genetics of three sexually dimorphic melanin-based traits in the barn owl (Tyto alba), in which females are on average darker reddish pheomelanic and display more and larger black eumelanic feather spots than males. The plumage traits with higher sex-linked inheritance showed lower heritability and genetic correlations, but contrary to prediction, these traits showed less pronounced sexual dimorphism. Strong offspring sexual dimorphism primarily resulted from daughters not expressing malelike melanin-based traits and from sons expressing femalelike traits to similar degrees as their sisters. We conclude that in the barn owl, polymorphism at autosomal genes rather than at sex-linked genes generate variation in sexual dimorphism in melanin-based traits.

Morph-specific genetic and environmental variation in innate and acquired immune response in a color polymorphic raptor.

Genetic color polymorphism is widespread in nature. There is an increasing interest in understanding the adaptive value of heritable color variation and trade-off resolution by differently colored individuals. Melanin-based pigmentation is often associated with variation in many different life history traits. These associations have recently been suggested to be the outcome of pleiotropic effects of the melanocortin system. Although pharmacological research supports that MC1R, a gene with a major role in vertebrate pigmentation, has important immunomodulatory effects, evidence regarding pleiotropy at MC1R in natural populations is still under debate. We experimentally assessed whether MC1R-based pigmentation covaries with both inflammatory and humoral immune responses in the color polymorphic Eleonora's falcon. By means of a cross-fostering experiment, we disentangled potential genetic effects from environmental effects on the covariation between coloration and immunity. Variation in both immune responses was primarily due to genetic factors via the nestlings' MC1R-related color genotype/phenotype, although environmental effects via the color morph of the foster father also had an influence. Overall, dark nestlings had lower immune responses than pale ones. The effect of the color morph of the foster father was also high, but in the opposite direction, and nestlings raised by dark eumelanic foster fathers had higher immune responses than those raised by pale foster fathers. Although we cannot completely discard alternative explanations, our results suggest that MC1R might influence immunity in this species. Morph-specific variation in immunity as well as pathogen pressure may therefore contribute to the long-term maintenance of genetic color polymorphism in natural populations.

The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study.

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHRadjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHRadjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHRadjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHRadjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHRadjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHRadjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

Multiplicative interaction of functional inflammasome genetic variants in determining the risk of gout.

INTRODUCTION: The acute gout flare results from a localised self-limiting innate immune response to monosodium urate (MSU) crystals deposited in joints in hyperuricaemic individuals. Activation of the caspase recruitment domain-containing protein 8 (CARD8) NOD-like receptor pyrin-containing 3 (NLRP3) inflammasome by MSU crystals and production of mature interleukin-1β (IL-1β) is central to acute gouty arthritis. However very little is known about genetic control of the innate immune response involved in acute gouty arthritis. Therefore our aim was to test functional single nucleotide polymorphism (SNP) variants in the toll-like receptor (TLR)-inflammasome-IL-1β axis for association with gout. METHODS: 1,494 gout cases of European and 863 gout cases of New Zealand (NZ) Polynesian (Māori and Pacific Island) ancestry were included. Gout was diagnosed by the 1977 ARA gout classification criteria. There were 1,030 Polynesian controls and 10,942 European controls including from the publicly-available Atherosclerosis Risk in Communities (ARIC) and Framingham Heart (FHS) studies. The ten SNPs were either genotyped by Sequenom MassArray or by Affymetrix SNP array or imputed in the ARIC and FHS datasets. Allelic association was done by logistic regression adjusting by age and sex with European and Polynesian data combined by meta-analysis. Sample sets were pooled for multiplicative interaction analysis, which was also adjusted by sample set. RESULTS: Eleven SNPs were tested in the TLR2, CD14, IL1B, CARD8, NLRP3, MYD88, P2RX7, DAPK1 and TNXIP genes. Nominally significant (P < 0.05) associations with gout were detected at CARD8 rs2043211 (OR = 1.12, P = 0.007), IL1B rs1143623 (OR = 1.10, P = 0.020) and CD14 rs2569190 (OR = 1.08; P = 0.036). There was significant multiplicative interaction between CARD8 and IL1B (P = 0.005), with the IL1B risk genotype amplifying the risk effect of CARD8. CONCLUSION: There is evidence for association of gout with functional variants in CARD8, IL1B and CD14. The gout-associated allele of IL1B increases expression of IL-1β - the multiplicative interaction with CARD8 would be consistent with a synergy of greater inflammasome activity (resulting from reduced CARD8) combined with higher levels of pre-IL-1β expression leading to increased production of mature IL-1β in gout.

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.

The genetic basis of color-related local adaptation in a ring-like colonization around the Mediterranean.

Uncovering the genetic basis of phenotypic variation and the population history under which it established is key to understand the trajectories along which local adaptation evolves. Here, we investigated the genetic basis and evolutionary history of a clinal plumage color polymorphism in European barn owls (Tyto alba). Our results suggest that barn owls colonized the Western Palearctic in a ring-like manner around the Mediterranean and meet in secondary contact in Greece. Rufous coloration appears to be linked to a recently evolved nonsynonymous-derived variant of the melanocortin 1 receptor (MC1R) gene, which according to quantitative genetic analyses evolved under local adaptation during or following the colonization of Central Europe. Admixture patterns and linkage disequilibrium between the neutral genetic background and color found exclusively within the secondary contact zone suggest limited introgression at secondary contact. These results from a system reminiscent of ring species provide a striking example of how local adaptation can evolve from derived genetic variation.