How accurate is the current picture of human genetic variation?

Our understanding of the distribution of worldwide human genomic diversity has greatly increased over recent years thanks to the availability of large data sets derived from short tandem repeats (STRs), insertion deletion polymorphisms (indels) and single nucleotide polymorphisms (SNPs). A concern, however, is that the current picture of worldwide human genomic diversity may be inaccurate because of biases in the selection process of genetic markers (so-called 'ascertainment bias'). To evaluate this problem, we first compared the distribution of genomic diversity between these three types of genetic markers in the populations from the HGDP-CEPH panel for evidence of bias or incongruities. In a second step, using a very relaxed set of criteria to prevent the intrusion of bias, we developed a new set of unbiased STR markers and compared the results against those from available panels. Contrarily to recent claims, our results show that the STR markers suffer from no discernible bias, and can thus be used as a baseline reference for human genetic diversity and population differentiation. The bias on SNPs is moderate compared to that on the set of indels analysed, which we recommend should be avoided for work describing the distribution of human genetic diversity or making inference on human settlement history.

Phylogeographical structure, postglacial recolonization and barriers to gene flow in the distinctive Valais chromosome race of the common shrew (Sorex araneus).

Using one male-inherited and eight biparentally inherited microsatellite markers, we investigate the population genetic structure of the Valais chromosome race of the common shrew (Sorex araneus) in the Central Alps of Europe. Unexpectedly, the Y-chromosome microsatellite suggests nearly complete absence of male gene flow among populations from the St-Bernard and Simplon regions (Switzerland). Autosomal markers also show significant genetic structuring among these two geographical areas. Isolation by distance is significant and possible barriers to gene flow exist in the study area. Two different approaches are used to better understand the geographical patterns and the causes of this structuring. Using a principal component analysis for which testing procedure exists, and partial Mantel tests, we show that the St-Bernard pass does not represent a significant barrier to gene flow although it culminates at 2469 m, close to the highest altitudinal record for this species. Similar results are found for the Simplon pass, indicating that both passes represented potential postglacial recolonization routes into Switzerland from Italian refugia after the last Pleistocene glaciations. In contrast with the weak effect of these mountain passes, the Rhône valley lowlands significantly reduce gene flow in this species. Natural obstacles (the large Rhône river) and unsuitable habitats (dry slopes) are both present in the valley. Moreover, anthropogenic changes to landscape structures are likely to have strongly reduced available habitats for this shrew in the lowlands, thereby promoting genetic differentiation of populations found on opposite sides of the Rhône valley.

Infantile spasms is associated with deletion of the MAGI2 gene on chromosome 7q11.23-q21.11.

Infantile spasms (IS) is the most severe and common form of epilepsy occurring in the first year of life. At least half of IS cases are idiopathic in origin, with others presumed to arise because of brain insult or malformation. Here, we identify a locus for IS by high-resolution mapping of 7q11.23-q21.1 interstitial deletions in patients. The breakpoints delineate a 500 kb interval within the MAGI2 gene (1.4 Mb in size) that is hemizygously disrupted in 15 of 16 participants with IS or childhood epilepsy, but remains intact in 11 of 12 participants with no seizure history. MAGI2 encodes the synaptic scaffolding protein membrane-associated guanylate kinase inverted-2 that interacts with Stargazin, a protein also associated with epilepsy in the stargazer mouse.

Use of forensic investigations in anti-doping.

The fight against doping is mainly focused on direct detection, using analytical methods for the detection of doping agents in biological samples. However, the World Anti-Doping Code also defines doping as possession, administration or attempted administration of prohibited substances or methods, trafficking or attempted trafficking in any prohibited substance or methods. As these issues correspond to criminal investigation, a forensic approach can help assessing potential violation of these rules. In the context of a rowing competition, genetic analyses were conducted on biological samples collected in infusion apparatus, bags and tubing in order to obtain DNA profiles. As no database of athletes' DNA profiles was available, the use of information from the location detection as well as contextual information were key to determine a population of suspected athletes and to obtain reference DNA profiles for comparison. Analysis of samples from infusion systems provided 8 different DNA profiles. The comparison between these profiles and 8 reference profiles from suspected athletes could not be distinguished. This case-study is one of the first where a forensic approach was applied for anti-doping purposes. Based on this investigation, the International Rowing Federation authorities decided to ban not only the incriminated athletes, but also the coaches and officials for 2 years.

Flexible social organization and high incidence of drifting in the sweat bee, Halictus scabiosae.

The very diverse social systems of sweat bees make them interesting models to study social evolution. Here we focus on the dispersal behaviour and social organization of Halictus scabiosae, a common yet poorly known species of Europe. By combining field observations and genetic data, we show that females have multiple reproductive strategies, which generates a large diversity in the social structure of nests. A detailed microsatellite analysis of 60 nests revealed that 55% of the nests contained the offspring of a single female, whereas the rest had more complex social structures, with three clear cases of multiple females reproducing in the same nest and frequent occurrence of unrelated individuals. Drifting among nests was surprisingly common, as 16% of the 122 nests in the overall sample and 44% of the nests with complex social structure contained females that had genotypes consistent with being full-sisters of females sampled in other nests of the population. Drifters originated from nests with an above-average productivity and were unrelated to their nestmates, suggesting that drifting might be a strategy to avoid competition among related females. The sex-specific comparison of genetic differentiation indicated that dispersal was male-biased, which would reinforce local resource competition among females. The pattern of genetic differentiation among populations was consistent with a dynamic process of patch colonization and extinction, as expected from the unstable, anthropogenic habitat of this species. Overall, our data show that H. scabiosae varies greatly in dispersal behaviour and social organization. The surprisingly high frequency of drifters echoes recent findings in wasps and bees, calling for further investigation of the adaptive basis of drifting in the social insects.

A sex-specific marker reveals male heterogamety in European tree frogs.

Most amphibians examined so far show undifferentiated sex chromosomes. The heterogametic sex's identity, usually revealed through indirect means, often varies among closely related species or even populations (as do sex-linkage groups), suggesting great evolutionary instability of the sex-determining genes. Here we take advantage of a sex-specific marker that amplifies in several related species of European tree frogs (Hyla arborea group) to disclose a homogeneous pattern of male heterogamety. Besides relevance for evolutionary studies of sex determination in amphibians, our results have potential for addressing practical issues in conservation biology because sex reversal by anthropogenic endocrine disruptors is considered one possible cause of amphibian decline.

Genetic differentiation of continental and island populations of Bombus terrestris (Hymenoptera: Apidae) in Europe.

Ten microsatellite loci and a partial sequence of the COII mitochondrial gene were used to investigate genetic differentiation in B. terrestris, a bumble bee of interest for its high-value crop pollination. The analysis included eight populations from the European continent, five from Mediterranean islands (six subspecies altogether) and one from Tenerife (initially described as a colour form of B. terrestris but recently considered as a separate species, B. canariensis). Eight of the 10 microsatellite loci displayed high levels of polymorphism in most populations. In B. terrestris populations, the total number of alleles detected per polymorphic locus ranged from 3 to 16, with observed allelic diversity from 3.8 +/- 0.5 to 6.5 +/- 1.4 and average calculated heterozygosities from 0.41 +/- 0.09 to 0.65 +/- 0.07. B. canariensis showed a significantly lower average calculated heterozygosity (0.12 +/- 0.08) and observed allelic diversity (1.5 +/- 0.04) as compared to both continental and island populations of B. terrestris. No significant differentiation was found among populations of B. terrestris from the European continent. In contrast, island populations were all significantly and most of them strongly differentiated from continental populations. B. terrestris mitochondrial DNA is characterized by a low nucleotide diversity: 0.18% +/- 0.07%, 0.20% +/- 0.04% and 0.27% +/- 0.04% for the continental populations, the island populations and all populations together, respectively. The only haplotype found in the Tenerife population differs by a single nucleotide substitution from the most common continental haplotype of B. terrestris. This situation, identical to that of Tyrrhenian islands populations and quite different from that of B. lucorum (15 substitutions between terrestris and lucorum mtDNA) casts doubts on the species status of B. canariensis. The large genetic distance between the Tenerife and B. terrestris populations estimated from microsatellite data result, most probably, from a severe bottleneck in the Canary island population. Microsatellite and mitochondrial DNA data call for the protection of the island populations of B. terrestris against importation of bumble bees of foreign origin which are used as crop pollinators.

Postglacial recolonization at a snail's pace (Trochulus villosus): confronting competing refugia hypotheses using model selection.

The localization of Last Glacial Maximum (LGM) refugia is crucial information to understand a species' history and predict its reaction to future climate changes. However, many phylogeographical studies often lack sampling designs intensive enough to precisely localize these refugia. The hairy land snail Trochulus villosus has a small range centred on Switzerland, which could be intensively covered by sampling 455 individuals from 52 populations. Based on mitochondrial DNA sequences (COI and 16S), we identified two divergent lineages with distinct geographical distributions. Bayesian skyline plots suggested that both lineages expanded at the end of the LGM. To find where the origin populations were located, we applied the principles of ancestral character reconstruction and identified a candidate refugium for each mtDNA lineage: the French Jura and Central Switzerland, both ice-free during the LGM. Additional refugia, however, could not be excluded, as suggested by the microsatellite analysis of a population subset. Modelling the LGM niche of T. villosus, we showed that suitable climatic conditions were expected in the inferred refugia, but potentially also in the nunataks of the alpine ice shield. In a model selection approach, we compared several alternative recolonization scenarios by estimating the Akaike information criterion for their respective maximum-likelihood migration rates. The 'two refugia' scenario received by far the best support given the distribution of genetic diversity in T. villosus populations. Provided that fine-scale sampling designs and various analytical approaches are combined, it is possible to refine our necessary understanding of species responses to environmental changes.

Phenotypic and molecular assessment of seven patients with 6p25 deletion syndrome: relevance to ocular dysgenesis and hearing impairment.

BACKGROUND: Thirty-nine patients have been described with deletions involving chromosome 6p25. However, relatively few of these deletions have had molecular characterization. Common phenotypes of 6p25 deletion syndrome patients include hydrocephalus, hearing loss, and ocular, craniofacial, skeletal, cardiac, and renal malformations. Molecular characterization of deletions can identify genes that are responsible for these phenotypes. METHODS: We report the clinical phenotype of seven patients with terminal deletions of chromosome 6p25 and compare them to previously reported patients. Molecular characterization of the deletions was performed using polymorphic marker analysis to determine the extents of the deletions in these seven 6p25 deletion syndrome patients. RESULTS: Our results, and previous data, show that ocular dysgenesis and hearing impairment are the two most highly penetrant phenotypes of the 6p25 deletion syndrome. While deletion of the forkhead box C1 gene (FOXC1) probably underlies the ocular dysgenesis, no gene in this region is known to be involved in hearing impairment. CONCLUSIONS: Ocular dysgenesis and hearing impairment are the two most common phenotypes of 6p25 deletion syndrome. We conclude that a locus for dominant hearing loss is present at 6p25 and that this locus is restricted to a region distal to D6S1617. Molecular characterization of more 6p25 deletion patients will aid in refinement of this locus and the identification of a gene involved in dominant hearing loss.

Estimating sex-specific dispersal rates with autosomal markers in hierarchically structured populations.

A recent study suggests that sex-specific dispersal rates can be quantitatively estimated on the basis of sex- and state-specific (pre- vs. postdispersal) F-statistics. In the present paper, we extend this approach to account for the hierarchical structure of natural populations, and we validate it through individual-based simulations. The model is applied to an empirical data set consisting of 536 individuals (males, females, and predispersal juveniles) of greater white-toothed shrews (Crocidura russula), sampled according to a hierarchical design and typed for seven autosomal microsatellite loci. From this dataset, dispersal is significantly female biased at the local scale (breeding-group level), but not at the larger scale (among local populations). We argue that selective pressures on dispersal are likely to depend on the spatial scale considered, and that short-distance dispersal should mainly respond to kin interactions (inbreeding or kin competition avoidance), which exert differential pressure on males and females.

Characterization and distribution of Pogonomyrmex harvester ant lineages with genetic caste determination.

Genetic caste determination has been described in two populations of Pogonomyrmex harvester ants, each comprising a pair of interbreeding lineages. Queens mate with males of their own and of the alternate lineage and produce two types of diploid offspring, those fertilized by males of the queens' lineage which develop into queens and those fertilized by males of the other lineage which develop into workers. Each of the lineages has been shown to be itself of hybrid origin between the species Pogonomyrmex barbatus and Pogonomyrmex rugosus, which both have typical, environmentally determined caste differentiation. In a large scale genetic survey across 35 sites in Arizona, New Mexico and Texas, we found that genetic caste determination associated with pairs of interbreeding lineages occurred frequently (in 26 out of the 35 sites). Overall, we identified eight lineages with genetic caste determination that always co-occurred in the same complementary lineage pairs. Three of the four lineage pairs appear to have a common origin while their relationship with the fourth remains unclear. The level of genetic differentiation among these eight lineages was significantly higher than the differentiation between P. rugosus and P. barbatus, which questions the appropriate taxonomic status of these genetic lineages. In addition to being genetically isolated from one another, all lineages with genetic caste determination were genetically distinct from P. rugosus and P. barbatus, even when colonies of interbreeding lineages co-occurred with colonies of either putative parent at the same site. Such nearly complete reproductive isolation between the lineages and the species with environmental caste determination might prevent the genetic caste determination system to be swept away by gene flow.

From the Apennines to the Alps: colonization genetics of the naturally expanding Italian wolf (Canis lupus) population

Wolves in Italy strongly declined in the past and were confined south of the Alps since the turn of the last century, reduced in the 1970s to approximately 100 individuals surviving in two fragmented subpopulations in the central-southern Apennines. The Italian wolves are presently expanding in the Apennines, and started to recolonize the western Alps in Italy, France and Switzerland about 16 years ago. In this study, we used a population genetic approach to elucidate some aspects of the wolf recolonization process. DNA extracted from 3068 tissue and scat samples collected in the Apennines (the source populations) and in the Alps (the colony), were genotyped at 12 microsatellite loci aiming to assess (i) the strength of the bottleneck and founder effects during the onset of colonization; (ii) the rates of gene flow between source and colony; and (iii) the minimum number of colonizers that are needed to explain the genetic variability observed in the colony. We identified a total of 435 distinct wolf genotypes, which showed that wolves in the Alps: (i) have significantly lower genetic diversity (heterozygosity, allelic richness, number of private alleles) than wolves in the Apennines; (ii) are genetically distinct using pairwise F(ST) values, population assignment test and Bayesian clustering; (iii) are not in genetic equilibrium (significant bottleneck test). Spatial autocorrelations are significant among samples separated up to c. 230 km, roughly correspondent to the apparent gap in permanent wolf presence between the Alps and north Apennines. The estimated number of first-generation migrants indicates that migration has been unidirectional and male-biased, from the Apennines to the Alps, and that wolves in southern Italy did not contribute to the Alpine population. These results suggest that: (i) the Alps were colonized by a few long-range migrating wolves originating in the north Apennine subpopulation; (ii) during the colonization process there has been a moderate bottleneck; and (iii) gene flow between sources and colonies was moderate (corresponding to 1.25-2.50 wolves per generation), despite high potential for dispersal. Bottleneck simulations showed that a total of c. 8-16 effective founders are needed to explain the genetic diversity observed in the Alps. Levels of genetic diversity in the expanding Alpine wolf population, and the permanence of genetic structuring, will depend on the future rates of gene flow among distinct wolf subpopulation fragments.

Genetic and karyotypic structure in the shrews of the Sorex araneus group: are they independent?

The species of the common shrew (Sorex araneus) group are morphologically very similar but exhibit high levels of karyotypic variation. Here we used genetic variation at 10 microsatellite markers in a data set of 212 individuals mostly sampled in the western Alps and composed of five karyotypic taxa (Sorex coronatus, Sorex antinorii and the S. araneus chromosome races Cordon, Bretolet and Vaud) to investigate the concordance between genetic and karyotypic structure. Bayesian analysis confirmed the taxonomic status of the three sampled species since individuals consistently grouped according to their taxonomical status. However, introgression can still be detected between S. antinorii and the race Cordon of S. araneus. This observation is consistent with the expected low karyotypic complexity of hybrids between these two taxa. Geographically based cryptic substructure was discovered within S. antinorii, a pattern consistent with the different postglaciation recolonization routes of this species. Additionally, we detected two genetic groups within S. araneus notwithstanding the presence of three chromosome races. This pattern can be explained by the probable hybrid status of the Bretolet race but also suggests a relatively low impact of chromosomal differences on genetic structure compared to historical factors. Finally, we propose that the current data set (available at http://www.unil.ch/dee/page7010_en.html#1) could be used as a reference by those wanting to identify Sorex individuals sampled in the western Alps.

Is the Gibraltar strait a barrier to gene flow for the bat Myotis myotis (Chiroptera: Vespertilionidae)?

Because of their role in limiting gene flow, geographical barriers like mountains or seas often coincide with intraspecific genetic discontinuities. Although the Strait of Gibraltar represents such a potential barrier for both plants and animals, few studies have been conducted on its impact on gene flow. Here we test this effect on a bat species (Myotis myotis) which is apparently distributed on both sides of the strait. Six colonies of 20 Myotis myotis each were sampled in southern Spain and northern Morocco along a linear transect of 1350 km. Results based on six nuclear microsatellite loci reveal no significant population structure within regions, but a complete isolation between bats sampled on each side of the strait. Variability at 600 bp of a mitochondrial gene (cytochrome b) confirms the existence of two genetically distinct and perfectly segregating clades, which diverged several million years ago. Despite the narrowness of the Gibraltar Strait (14 km), these molecular data suggest that neither males, nor females from either region have ever reproduced on the opposite side of the strait. Comparisons of molecular divergence with bats from a closely related species (M. blythii) suggest that the North African clade is possibly a distinct taxon warranting full species rank. We provisionally refer to it as Myotis cf punicus Felten 1977, but a definitive systematic understanding of the whole Mouse-eared bat species complex awaits further genetic sampling, especially in the Eastern Mediterranean areas.

Suggestive linkage to chromosome 1q for bone mineral apparent density in Brazilian sister adolescents.

OBJECTIVE: To investigate linkage to chromosome 1q and 11q region for lumbar spine, femoral neck and total body BMD and volumetric BMD in Brazilian sister adolescents aged 10-20-year-old and 57 mothers. METHODS: We evaluated 161 sister pairs (n=329) aged 10-20 years old and 57 of their mothers in this study. Physical traits and lifestyle factors were collected as covariates for lumbar spine (LS), femoral neck (FN) and total body (TB) BMD and bone mineral apparent density (BMAD). We selected nine microsatellite markers in chromosome 1q region (spanning nearly 33cM) and eight in chromosome 11q region (spanning nearly 34cM) to perform linkage analysis. RESULTS: The highest LOD score values obtained from our data were in sister pairs LS BMAD analysis. Their values were: 1.32 (P<0.006), 2.61 (P<0.0002) and 2.44 (P<0.0004) in D1S218, D1S2640 and D1S2623 markers, respectively. No significant LOD score was found with LS and FN BMD/BMAD in chromosome 11q region. Only TB BMD showed significant linkage higher than 1.0 for chromosome 11q region in the markers D11S4191 and D11S937. DISCUSSION/CONCLUSIONS: Our results provided suggestive linkage for LS BMAD at D1S2640 marker in adolescent sister pairs and suggest a possible candidate gene (LHX4) related to adolescent LS BMAD in this region. These results reinforce chromosome 1q21-23 as a candidate region to harbor one or more bone formation/maintenance gene. In the other hand, it did not repeat for chromosome 11q12-13 in our population.