Mikel Iriondo

Gene flow within the M evolutionary lineage of Apis mellifera: role of the Pyrenees, isolation by distance and post-glacial re-colonization routes in the western Europe

We present a population genetic study focused on the two subspecies of the M evolutionary lineage, A. m. mellifera and A. m. iberiensis. Nuclear and mtDNA variation was analysed in 27 bee populations from the Iberian Peninsula, France and Belgium. Microsatellite data provides compelling evidence of a barrier to neutral gene flow at the Pyrenees. In addition, they suggest isolation by distance between populations of the M lineage. Mitochondrial data support the hypothesis that the Iberian Peninsula served as glacial refugia for the honeybees of western Europe. They show two paths of post-glacial re-colonization in the extremes of the Pyrenees and suggest that the western path was more significant in the post-glacial re-colonization process. Thus, we report here on three main factors for mellifera and iberiensis subspecies differentiation: the Pyrenean barrier, isolation by distance and the post-glacial re-colonization process.ZusammenfassungDie Unterarten der Honigbiene (Apis mellifera) werden in die 5 evolutionäre Linien A (African), C (northern Mediterranean), M (western Europe), O (Oriental) and Y (Yemenitica) gruppiert. Die in dieser Studie untersuchte evolutionäre Linie M enthält zwei Unterarten: A. m. mellifera wird von Frankreich bis zu den Bergen des Ural gefunden und A. m. iberiensis ist auf der iberischen Halbinsel verbreitet. Obwohl allgemein angenommen wird, dass die Pyrenäen ein bedeutendes Hindernis für den Genfluss zwischen A. m. mellifera und A. m. iberiensis darstellt, wurde ein solcher Barriereeffekt bislang nicht nachgewiesen. Die Existenz genetischer Gefälle vom Süden der Iberischen Halbinsel bis zum nördlichen Europa und eine ungewisse taxonomische Zuordnung einiger Populationen in den Pyrenäen tragen zu der Unsicherheit bezüglich der Rolle der Pyrenäen als genetische Barriere bei. Andererseits wird seit längerem angenommen, dass die Iberische Halbinsel während der Eiszeit als Refugium für die westliche Honigbiene diente (Ruttner, 1952, 1988) und eine nacheiszeitliche Wiederbesiedlung von Nordeuropa wird von mehreren Autoren unterstützt (Garnery et al., 1998a,b; Franck et al., 1998, 2000b). Gegenstand dieser Untersuchung war es, einen potentiell isolierenden Effekt der Pyrenäen nachzuweisen, neue Daten zu dem Differenzierungsprozess der zwei Unterarten beizutragen und den Ablauf der Wiederbesiedlung zu untersuchen. Wir untersuchten 1398 Völker aus 27 Populationen der Iberischen Halbinsel sowie aus Frankreich und Belgien auf Variation an 10 Mikrosatellitenloci und der COI-COII intergenischen Region der mtDNA. Wir verwendeten verschiedene Arten statistischer Analysen wie die DA genetische Distanzmatrix, neighbor-joining trees, Korrelationen erster Ordnung und partielle Korrelationen, AMOVA, Analyse räumlicher Autokorrelationen und COCOPAN für Mikrosatellitendaten. Die Ergebnisse zeigten eine Isolation zwischen den verschiedenen Populationen der M Linie durch die Entfernung auf und lieferten sehr deutliche Hinweise auf eine Barriere für den neutralen Genfluß bei den Pyrenäen. Die Verteilung der mtDNA Haplotypen bestätigte das Iberische Refugium der Westeuropäischen Honigbiene in der Eiszeit. Wir konnten zwei verschiedene Wege der nacheiszeitlichen Wiederbesiedlung von der Iberischen Halbinsel aus an den beiden Enden der Pyrenäen ableiten. Es gab deutliche Unterschiede in der Verteilung der Mitotypen zwischen den westlichen und östlichen Enden der Pyrenäen, diese legten nahe, dass der westliche Weg für den nacheiszeitlichen Widerbesiedlungsprozess wichtiger war. Nach der in dieser Untersuchung beobachteten hohen Variabilität der M Mitotypen südlich der Pyrenäen könnten diese eine nützliche genetische Ressource für die Konservation der Westeuropäischen Honigbienen darstellen. Der nacheiszeitliche Wiederbesiedlungsverlauf, die Isolation durch die Entfernung und die von den Pyrenäen gebildete Verbreitungsbarriere sind Einflüsse, die zu der Ausbildung der Unterarten A. m. mellifera und A. m. iberiensis beigetragen haben.

Both geometric morphometric and microsatellite data consistently support the differentiation of the Apis mellifera M evolutionary branch

Traditional morphometrics, allozymes, and mitochondrial data have supported a close relationship between the M branch subspecies A. m. iberiensis and the North African subspecies (A branch). However, studies using nuclear DNA markers have revealed a clear distinction between the latter and the two European M branch subspecies. In help resolve this paradox, we analyzed 663 colonies from six European and African subspecies. A geometric morphometrics approach was applied to the analysis of wing shape, and the results were compared with data of six microsatellite loci. Both data sets were found to be highly consistent and corroborated a marked divergence of West European subspecies from North African ones. This supports the hypothesis that the presence of the African lineage mitotype in Iberian honey bee populations is likely the consequence of secondary introductions, with a minimal African influence within the current Iberian genetic background. Wing geometric morphometrics appears more appropriate than mitochondrial DNA analysis or traditional morphometrics in the screening and identification of the Africanization process.

Multiple SNP Markers Reveal Fine-Scale Population and Deep Phylogeographic Structure in European Anchovy (Engraulis encrasicolus L.)

Geographic surveys of allozymes, microsatellites, nuclear DNA (nDNA) and mitochondrial DNA (mtDNA) have detected several genetic subdivisions among European anchovy populations. However, these studies have been limited in their power to detect some aspects of population structure by the use of a single or a few molecular markers, or by limited geographic sampling. We use a multi-marker approach, 47 nDNA and 15 mtDNA single nucleotide polymorphisms (SNPs), to analyze 626 European anchovies from the whole range of the species to resolve shallow and deep levels of population structure. Nuclear SNPs define 10 genetic entities within two larger genetically distinctive groups associated with oceanic variables and different life-history traits. MtDNA SNPs define two deep phylogroups that reflect ancient dispersals and colonizations. These markers define two ecological groups. One major group of Iberian-Atlantic populations is associated with upwelling areas on narrow continental shelves and includes populations spawning and overwintering in coastal areas. A second major group includes northern populations in the North East (NE) Atlantic (including the Bay of Biscay) and the Mediterranean and is associated with wide continental shelves with local larval retention currents. This group tends to spawn and overwinter in oceanic areas. These two groups encompass ten populations that differ from previously defined management stocks in the Alboran Sea, Iberian-Atlantic and Bay of Biscay regions. In addition, a new North Sea-English Channel stock is defined. SNPs indicate that some populations in the Bay of Biscay are genetically closer to North Western (NW) Mediterranean populations than to other populations in the NE Atlantic, likely due to colonizations of the Bay of Biscay and NW Mediterranean by migrants from a common ancestral population. Northern NE Atlantic populations were subsequently established by migrants from the Bay of Biscay. Populations along the Iberian-Atlantic coast appear to have been founded by secondary waves of migrants from a southern refuge.

High-density SNP genotyping detects homogeneity of Spanish and French Basques, and confirms their genomic distinctiveness from other European populations

A recent study reported that Basques do not constitute a genetically distinct population, and that Basques from Spanish and French provinces do not show significant genetic similarity. These conclusions disagree with numerous previous studies, and are not consistent with the historical and linguistic evidence that supports the distinctiveness of Basques. In order to further investigate this controversy, we have genotyped 83 Spanish Basque individuals and used these data to infer population structure based on more than 60,000 single nucleotide polymorphisms of several European populations. Here, we present the first high-throughput analysis including Basques from Spanish and French provinces, and show that all Basques constitute a homogeneous group that can be clearly differentiated from other European populations.

Genetic association between bovine NOD2 polymorphisms and infection by Mycobacterium avium subsp. paratuberculosis in Holstein‐Friesian cattle

Nucleotide-Binding Oligomerization Domain 2 (NOD2) has been reported to be a candidate gene for Mycobacterium avium subsp. paratuberculosis (MAP) infection in a Bos taurus × Bos indicus mixed breed based on a genetic association with the c.2197T>C single nucleotide polymorphism (SNP). Nevertheless, this SNP has also been reported to be monomorphic in the B. taurus species. In the present work, 18 SNPs spanning the bovine NOD2 gene have been analysed in a genetic association study of two independent populations of Holstein-Friesian cattle. We found that the C allele of SNP c.*1908C>T, located in the 3′-UTR region of the gene, is significantly more frequent in infected animals than in healthy ones, which supports the idea that the bovine NOD2 gene plays a role in susceptibility to MAP infection. However, in silico analyses of the NOD2 nucleotide sequence did not yield definitive data about a possible direct effect of SNP c.*1908C>T on susceptibility to infection and led us to consider its linkage disequilibrium with the causative variant. A more exhaustive genetic association study including all putative, functional SNPs from this gene and subsequent functional analyses needs to be conducted to achieve a more complete understanding of how different variants of NOD2 may affect susceptibility to MAP infection in cattle.

Genetic polymorphisms of the Basques from Gipuzkoa: genetic heterogeneity of the Basque population

A random sample of 586 Basque individuals from the province of Gipuzkoa was studied for 16 genetic systems: A1A2B0, Rh, MNSs, P, Lewis, Duffy, Kell, GC, TF, AAT, ACP, AK, ADA, ESD, HP and PGM1. The results of this study indicate that the Basque population of Gipuzkoa presents certain differential values with respect to other Basque series, such as maximum values for RH*cde, AK*2 and PGM1*2+ and minimum for PGM1*1-, while the remaining alleles are located within the range of values found in the Basque population to date. It is suggested that there is intraprovincial heterogeneity, as described for Bizkaia by Aguirre et al. in 1991, and the existence of heterogeneity within the Basque population on an inter-provincial level, backing up previous studies in this respect (by Aguirre et al. in 1989 and Manzano et al. 1993).

SP110 as a novel susceptibility gene for Mycobacterium avium subspecies paratuberculosis infection in cattle

The intracellular pathogen resistance 1 (Ipr1) gene has been reported to play a role in mediating innate immunity in a mouse model of Mycobacterium tuberculosis infection, and polymorphisms of its human ortholog, SP110 nuclear body protein, have been suggested to be associated with tuberculosis. Thus, the bovine SP110 gene was considered to be a promising candidate for a genetic association study of bovine paratuberculosis, or Johnes disease, a chronic granulomatous enteritis caused by Mycobacterium avium ssp. paratuberculosis (MAP). Initially, single nucleotide polymorphisms (SNP) within the bovine SP110 gene were identified, and subsequently a population-based genetic association study was carried out. Seventeen new SNP along the SP110 gene were identified in Holstein-Friesian cattle, and 6 more were compiled from public databases. A total of 14 SNP were included in the association study of 2 independent populations. The SNP c.587A>G was found to be significantly associated with MAP infection, with the major allele A appearing to confer greater disease susceptibility in one of the analyzed populations. In addition, 2 haplotypes containing this SNP were also found to be associated with infection in the same population. The SNP c.587A>G is a nonsynonymous mutation that causes an amino acid change in codon 196 from asparagine to serine. In silico analyses point to SNP c.587A>G as a putative causal variant for susceptibility to MAP infection. The elucidation of the precise mechanism by which this SNP can exert its effect in the protein and, as a result, in the risk of infection, requires future functional analyses. Likewise, the absence of genetic association in one of the analyzed populations renders it necessary to carry out this study in other independent populations, with the aim of substantiating the repeatability of the present results. Nevertheless, the present results deepen our understanding of the genetic basis of susceptibility and resistance mechanisms related to MAP infection in cattle and, in turn, constitute a step forward toward the implementation of marker-assisted selection in breeding programs aimed at controlling paratuberculosis.

Structuring the Genetic Heterogeneity of the Basque Population: A View from Classical Polymorphisms

In this study we analyze 18 classical polymorphisms (ABO, Rh, MNSs, Lewis, P, Duffy, Kell, ADA, ESD, PGM1, PGD, AK1, ACP1, GLO1, HP, GC, TF, and PI) in over 2000 autochthonous individuals from 14 natural districts in three provinces of the Basque Country (Alava, Guipuzcoa, and Biscay). Heterogeneity analysis via the 2test and a calculation of indicate that there is significant genetic heterogeneity between the Basque districts. The R> matrix informs us that this heterogeneity is not significantly concentrated in a single district or in the districts of a single province, but is rather distributed among several districts belonging to the three provinces analyzed. We undertake to assess the influence of various historical, geographical, and cultural factors on the genetic structure of the Basque population. Analysis suggests that allele distribution is geographically patterned in the Basque Country. The gradient distributions observed in the case of some alleles (ABO*O, RH*cDE, RH*cde, MNS*MS, and ACP1*C) on the basis of Morans autocorrelation coefficient I, along with the influence of the two main travel routes through the Basque Country (western route through Bilbao and eastern route through Vitoria), suggest that the gene flow tends toward the coast. As regards other factors considered (administrative division, repopulation processes, linguistic heterogeneity, and north vs. south cultural heterogeneity), we detected only a certain influence exerted by an old tribal differentiation (2000 B.P<.), which would diminish with the passage of time.

New Nuclear SNP Markers Unravel the Genetic Structure and Effective Population Size of Albacore Tuna (Thunnus alalunga)

In the present study we have investigated the population genetic structure of albacore (Thunnus alalunga, Bonnaterre 1788) and assessed the loss of genetic diversity, likely due to overfishing, of albacore population in the North Atlantic Ocean. For this purpose, 1,331 individuals from 26 worldwide locations were analyzed by genotyping 75 novel nuclear SNPs. Our results indicated the existence of four genetically homogeneous populations delimited within the Mediterranean Sea, the Atlantic Ocean, the Indian Ocean and the Pacific Ocean. Current definition of stocks allows the sustainable management of albacore since no stock includes more than one genetic entity. In addition, short- and long-term effective population sizes were estimated for the North Atlantic Ocean albacore population, and results showed no historical decline for this population. Therefore, the genetic diversity and, consequently, the adaptive potential of this population have not been significantly affected by overfishing.

Thrombotic Antiphospholipid Syndrome Shows Strong Haplotypic Association with SH2B3-ATXN2 Locus

Background Thrombotic antiphospholipid syndrome is defined as a complex form of thrombophilia that is developed by a fraction of antiphospholipid antibody (aPLA) carriers. Little is known about the genetic risk factors involved in thrombosis development among aPLA carriers. Methods To identify new loci conferring susceptibility to thrombotic antiphospholipid syndrome, a two-stage genotyping strategy was performed. In stage one, 19,000 CNV loci were genotyped in 14 thrombotic aPLA+ patients and 14 healthy controls by array-CGH. In stage two, significant CNV loci were fine-mapped in a larger cohort (85 thrombotic aPLA+, 100 non-thrombotic aPLA+ and 569 healthy controls). Results Array-CGH and fine-mapping analysis led to the identification of 12q24.12 locus as a new susceptibility locus for thrombotic APS. Within this region, a TAC risk haplotype comprising one SNP in SH2B3 gene (rs3184504) and two SNPs in ATXN2 gene (rs10774625 and rs653178) exhibited the strongest association with thrombotic antiphospholipid syndrome (p-value = 5,9 × 10−4 OR 95% CI 1.84 (1.32–2.55)). Conclusion The presence of a TAC risk haplotype in ATXN2-SH2B3 locus may contribute to increased thrombotic risk in aPLA carriers.