Genetic structure, relationships and admixture with wild relatives in native pig breeds from Iberia and its islands

Background: Native pig breeds in the Iberian Peninsula are broadly classified as belonging to either the Celtic or the Mediterranean breed groups, but there are other local populations that do not fit into any of these groups. Most of the native pig breeds in Iberia are in danger of extinction, and the assessment of their genetic diversity and population structure, relationships and possible admixture between breeds, and the appraisal of conservation alternatives are crucial to adopt appropriate management strategies. Methods: A panel of 24 microsatellite markers was used to genotype 844 animals representing the 17 most important native swine breeds and wild populations existing in Portugal and Spain and various statistical tools were applied to analyze the results. Results: Genetic diversity was high in the breeds studied, with an overall mean of 13.6 alleles per locus and an average expected heterozygosity of 0.80. Signs of genetic bottlenecks were observed in breeds with a small census size, and population substructure was present in some of the breeds with larger census sizes. Variability among breeds accounted for about 20% of the total genetic diversity, and was explained mostly by differences among the Celtic, Mediterranean and Basque breed groups, rather than by differences between domestic and wild pigs. Breeds clustered closely according to group, and proximity was detected between wild pigs and the Mediterranean cluster of breeds. Most breeds had their own structure and identity, with very little evidence of admixture, except for the Retinto and Entrepelado varieties of the Mediterranean group, which are very similar. Genetic influence of the identified breed clusters extends beyond the specific geographical areas across borders throughout the Iberian Peninsula, with a very sharp transition from one breed group to another. Analysis of conservation priorities confirms that the ranking of a breed for conservation depends on the emphasis placed on its contribution to the betweenand within-breed components of genetic diversity. Conclusions: Native pig breeds in Iberia reveal high levels of genetic diversity, a solid breed structure and a clear organization in well-defined clusters.

Population genetics of C. Maenas : oceanography and larval dispersal

Decifrar a complexa interacção entre os ciclos de vida de espécies marinhas e a oceanografia revela-se fundamental para a compreensão do fluxo genético e da conectividade no meio marinho. Nas espécies marinhas com desenvolvimento indirecto o fluxo de genes entre populações depende da distância que separa as populações, bem como da interacção entre a duração do desenvolvimento larvar, do comportamento das larvas e dos padrões de circulação oceânica. A conectividade larvar influencia uma variedade de processos como a dinâmica de stocks e de populações, a distribuição e limites geográficos das espécies, a estrutura genética das populações e a dispersão de espécies invasivas e reveste-se consequentemente de uma importância fundamental na identificação das unidades populacionais evolucionariamente relevantes e para a gestão e conservação marinhas. Os marcadores genéticos e os Modelos Individuais Acoplados a Modelos Físico-Biológicos (“ICPBMs”) são actualmente ferramentas fundamentais para o estudo dos padrões de dispersão larvar e para avaliar o nível de conectividade populacional. A presente tese respeita à avaliação das escalas espaciais de conectividade de populações de uma espécie costeira, o caranguejo Carcinus maenas, e utiliza conjuntamente informação de marcadores genéticos, análise de séries temporais de fornecimento de larvas e um modelo numérico de circulação oceânica. O primeiro capítulo introduz a temática da conectividade em espécies marinhas e inclui algumas referências aos métodos moleculares, analíticos e de modelação seguidos ao longo da tese. Através da utilização de múltiplas ferramentas – avaliação da estrutura genética geográfica de C. maenas na sua distribuição nativa com recurso a marcadores de DNA (microssatélites) (Capítulo 2), avaliação da estrutura genética temporal das larvas que formam os eventos de fornecimento larvar à Ria de Aveiro, NW Portugal (Capítulo 3), descrição da variabilidade inter-anual do fornecimento larvar à Ria de Aveiro, NW Portugal (Capítulo 4) e validação de um modelo ICPBM que descreve os padrões observados de fornecimento (Capítulo 5) – esta tese espera poder contribuir para uma melhor compreensão dos mecanismos que regulam o fluxo de genes e a conectividade entre populações de organismos marinhos. No Capítulo 6 são apresentadas as principais conclusões da investigação. A análise genética com recurso a microssatélites indicou que as populações de C. maenas são geneticamente homogéneas ao longo de várias centenas de km, dentro da distribuição nativa da espécie. Paralelamente, não foram encontrados indícios da existência de reprodução por “sweepstakes” em C. maenas de populações da costa oeste da Península Ibérica, visto que não se obtiveram diferenças genéticas significativas entre os eventos larvares. Também não se encontrou qualquer estrutura familiar entre as larvas que formam cada episódio de fornecimento, e não houve nenhuma redução significativa da variabilidade genética das larvas quando comparada com a de caranguejos adultos. A análise de séries temporais de suprimento de larvas na Ria de Aveiro em cinco anos estudados indica que este é um fenómeno episódico e variável, sendo os maiores episódios de fornecimento coincidentes com as marés vivas e acentuados por fortes ventos de sul. O modelo ICPBM foi validado com sucesso e parece fornecer uma estimativa realística das escalas espaciais e temporais de dispersão larvar, de acordo com as observações da estrutura genética e da ausência de reprodução por “sweepstake” em C. maenas da costa oeste da Península Ibérica

Genetic structure at range edge: Low diversity and high inbreeding in Southeast Asian mangrove (Avicennia marina) populations

Understanding the genetic composition and mating systems of edge populations provides important insights into the environmental and demographic factors shaping species’ distribution ranges. We analysed samples of the mangrove Avicennia marina from Vietnam, northern Philippines and Australia, with microsatellite markers. We compared genetic diversity and structure in edge (Southeast Asia, and Southern Australia) and core (North and Eastern Australia) populations, and also compared our results with previously published data from core and southern edge populations. Comparisons highlighted significantly reduced gene diversity and higher genetic structure in both margins compared to core populations, which can be attributed to very low effective population size, pollinator scarcity and high environmental pressure at distribution margins. The estimated level of inbreeding was significantly higher in northeastern populations compared to core and southern populations. This suggests that despite the high genetic load usually associated with inbreeding, inbreeding or even selfing may be advantageous in margin habitats due to the possible advantages of reproductive assurance, or local adaptation. The very high level of genetic structure and inbreeding show that populations of A. marina are functioning as independent evolutionary units more than as components of a metapopulation system connected by gene flow. The combinations of those characteristics make these peripheral populations likely to develop local adaptations and therefore to be of particular interest for conservation strategies as well as for adaptation to possible future environmental changes.

Population and landscape genetics of the stone marten and red fox in Portugal : implications for conservation management of common carnivores

Tese de doutoramento, Biologia (Biologia da Conservação), Universidade de Lisboa, Faculdade de Ciências, 2015

Proximity-dependent pollen performance in Silene vulgaris.

BACKGROUND AND AIMS: Pollen and seed dispersal in herbaceous insect-pollinated plants are often restricted, inducing strong population structure. To what extent this influences mating within and among patches is poorly understood. This study investigates the influence of population structure on pollen performance using controlled pollinations and genetic markers. METHODS: Population structure was investigated in a patchily distributed population of gynodioecious Silene vulgaris in Switzerland using polymorphic microsatellite markers. Experimental pollinations were performed on 21 hermaphrodite recipients using pollen donors at three spatial scales: (a) self-pollination; (b) within-patch cross-pollinations; and (c) between-patch cross-pollinations. Pollen performance was then compared with respect to crossing distance. KEY RESULTS: The population of S. vulgaris was characterized by a high degree of genetic sub-structure, with neighbouring plants more related to one another than to distant individuals. Inbreeding probably results from both selfing and biparental inbreeding. Pollen performance increased with distance between mates. Between-patch pollen performed significantly better than both self- and within-patch pollen donors. However, no significant difference was detected between self- and within-patch pollen donors. CONCLUSIONS: The results suggest that population structure in animal-pollinated plants is likely to influence mating patterns by favouring cross-pollinations between unrelated plants. However, the extent to which this mechanism could be effective as a pre-zygotic barrier preventing inbred mating depends on the patterns of pollinator foraging and their influence on pollen dispersal.

Clinical and genetic investigation of a large Tunisian family with complete achromatopsia: identification of a new nonsense mutation in GNAT2 gene.

Complete achromatopsia is a rare autosomal recessive disease associated with CNGA3, CNGB3, GNAT2 and PDE6C mutations. This retinal disorder is characterized by complete loss of color discrimination due to the absence or alteration of the cones function. The purpose of the present study was the clinical and the genetic characterization of achromatopsia in a large consanguineous Tunisian family. Ophthalmic evaluation included a full clinical examination, color vision testing and electroretinography. Linkage analysis using microsatellite markers flanking CNGA3, CNGB3, GNAT2 and PDE6C genes was performed. Mutations were screened by direct sequencing. A total of 12 individuals were diagnosed with congenital complete achromatopsia. They are members of six nuclear consanguineous families belonging to the same large consanguineous family. Linkage analysis revealed linkage to GNAT2. Mutational screening of GNAT2 revealed three intronic variations c.119-69G>C, c.161+66A>T and c.875-31G>C that co-segregated with a novel mutation p.R313X. An identical GNAT2 haplotype segregating with this mutation was identified, indicating a founder mutation. All patients were homozygous for the p.R313X mutation. This is the first report of the clinical and genetic investigation of complete achromatopsia in North Africa and the largest family with recessive achromatopsia involving GNAT2; thus, providing a unique opportunity for genotype-phenotype correlation for this extremely rare condition.

Microsatellites can be misleading: an empirical and simulation study.

It has been long recognized that highly polymorphic genetic markers can lead to underestimation of divergence between populations when migration is low. Microsatellite loci, which are characterized by extremely high mutation rates, are particularly likely to be affected. Here, we report genetic differentiation estimates in a contact zone between two chromosome races of the common shrew (Sorex araneus), based on 10 autosomal microsatellites, a newly developed Y-chromosome microsatellite, and mitochondrial DNA. These results are compared to previous data on proteins and karyotypes. Estimates of genetic differentiation based on F- and R-statistics are much lower for autosomal microsatellites than for all other genetic markers. We show by simulations that this discrepancy stems mainly from the high mutation rate of microsatellite markers for F-statistics and from deviations from a single-step mutation model for R-statistics. The sex-linked genetic markers show that all gene exchange between races is mediated by females. The absence of male-mediated gene flow most likely results from male hybrid sterility.

Study of a Swiss dopa-responsive dystonia family with a deletion in GCH1: redefining DYT14 as DYT5.

OBJECTIVE: To report the study of a multigenerational Swiss family with dopa-responsive dystonia (DRD). METHODS: Clinical investigation was made of available family members, including historical and chart reviews. Subject examinations were video recorded. Genetic analysis included a genome-wide linkage study with microsatellite markers (STR), GTP cyclohydrolase I (GCH1) gene sequencing, and dosage analysis. RESULTS: We evaluated 32 individuals, of whom 6 were clinically diagnosed with DRD, with childhood-onset progressive foot dystonia, later generalizing, followed by parkinsonism in the two older patients. The response to levodopa was very good. Two additional patients had late onset dopa-responsive parkinsonism. Three other subjects had DRD symptoms on historical grounds. We found suggestive linkage to the previously reported DYT14 locus, which excluded GCH1. However, further study with more stringent criteria for disease status attribution showed linkage to a larger region, which included GCH1. No mutation was found in GCH1 by gene sequencing but dosage methods identified a novel heterozygous deletion of exons 3 to 6 of GCH1. The mutation was found in seven subjects. One of the patients with dystonia represented a phenocopy. CONCLUSIONS: This study rules out the previously reported DYT14 locus as a cause of disease, as a novel multiexonic deletion was identified in GCH1. This work highlights the necessity of an accurate clinical diagnosis in linkage studies as well as the need for appropriate allele frequencies, penetrance, and phenocopy estimates. Comprehensive sequencing and dosage analysis of known genes is recommended prior to genome-wide linkage analysis.

Male body size and breeding tubercles are both linked to intra-sexual dominance and reproductive success in the minnow

Male dominance hierarchies are usually linked to relative body size and to weapon size, that is, to determinants of fighting ability. Secondary sexual characters that are not directly used as weapons could still be linked to dominance if they reveal determination or overall health and vigour and hence, indirectly, fighting ability. We studied the mating behaviour of the minnow, Phoxinus phoxinus, a cyprinid fish in which males develop breeding tubercles during the spawning season. The function of these breeding tubercles is still not clear. Using microsatellite markers, we determined male reproductive success under controlled conditions. The minnows were territorial and quickly established a dominance hierarchy at the beginning of the spawning season. Dominance was strongly and positively linked to fertilization success. Although body size and number of breeding tubercles were not significantly correlated in our sample, both large males and males with many breeding tubercles were more dominant and achieved higher fertilization success than small males or males with few tubercles. We found multimale fertilization in most clutches, suggesting that sperm competition is important in this species. Females showed behaviour that may be linked to spawning decision, that is, male dominance might not be the only determinant of male reproductive success in minnows

Persistent unequal sex ratio in a population of grayling (salmonidae) and possible role of temperature increase.

In some fishes, water chemistry or temperature affects sex determination or creates sex-specific selection pressures. The resulting population sex ratios are hard to predict from laboratory studies if the environmental triggers interact with other factors, whereas in field studies, singular observations of unusual sex ratios may be particularly prone to selective reporting. Long-term monitoring largely avoids these problems. We studied a population of grayling (Thymallus thymallus) in Lake Thun, Switzerland, that has been monitored since 1948. Samples of spawning fish have been caught about 3 times/week around spawning season, and water temperature at the spawning site has been continuously recorded since 1970. We used scale samples collected in different years to determine the average age of spawners (for life-stage specific analyses) and to identify the cohort born in 2003 (an extraordinarily warm year). Recent tissue samples were genotyped on microsatellite markers to test for genetic bottlenecks in the past and to estimate the genetically effective population size (N(e) ). Operational sex ratios changed from approximately 65% males before 1993 to approximately 85% males from 1993 to 2011. Sex ratios correlated with the water temperatures the fish experienced in their first year of life. Sex ratios were best explained by the average temperature juvenile fish experienced during their first summer. Grayling abundance is declining, but we found no evidence of a strong genetic bottleneck that would explain the apparent lack of evolutionary response to the unequal sex ratio. Results of other studies show no evidence of endocrine disruptors in the study area. Our findings suggest temperature affects population sex ratio and thereby contributes to population decline. Persistencia de Proporción de Sexos Desigual en una Población de Tímalos (Salmonidae) y el Posible Papel del Incremento de la Temperatura.

Haplotype analysis suggest common founders in carriers of the recurrent BRCA2 mutation, 3398delAAAAG, in French Canadian hereditary breast and/ovarian cancer families

Affiliation: Faculté de médicine, Université de Montréal

Genome wide search for genetic determinants of habitual alcohol, tobacco and coffee use, obesity-related traits, response to mental and physical stress and hemodynamic traits

Les habitudes de consommation de substances psychoactives, le stress, l’obésité et les traits cardiovasculaires associés seraient en partie reliés aux mêmes facteurs génétiques. Afin d’explorer cette hypothèse, nous avons effectué, chez 119 familles multi-générationnelles québécoises de la région du Saguenay-Lac-St-Jean, des études d’association et de liaison pangénomiques pour les composantes génétiques : de la consommation usuelle d’alcool, de tabac et de café, de la réponse au stress physique et psychologique, des traits anthropométriques reliés à l’obésité, ainsi que des mesures du rythme cardiaque (RC) et de la pression artérielle (PA). 58000 SNPs et 437 marqueurs microsatellites ont été utilisés et l’annotation fonctionnelle des gènes candidats identifiés a ensuite été réalisée. Nous avons détecté des corrélations phénotypiques significatives entre les substances psychoactives, le stress, l’obésité et les traits hémodynamiques. Par exemple, les consommateurs d’alcool et de tabac ont montré un RC significativement diminué en réponse au stress psychologique. De plus, les consommateurs de tabac avaient des PA plus basses que les non-consommateurs. Aussi, les hypertendus présentaient des RC et PA systoliques accrus en réponse au stress psychologique et un indice de masse corporelle (IMC) élevé, comparativement aux normotendus. D’autre part, l’utilisation de tabac augmenterait les taux corporels d’épinéphrine, et des niveaux élevés d’épinéphrine ont été associés à des IMC diminués. Ainsi, en accord avec les corrélations inter-phénotypiques, nous avons identifié plusieurs gènes associés/liés à la consommation de substances psychoactives, à la réponse au stress physique et psychologique, aux traits reliés à l’obésité et aux traits hémodynamiques incluant CAMK4, CNTN4, DLG2, DAG1, FHIT, GRID2, ITPR2, NOVA1, NRG3 et PRKCE. Ces gènes codent pour des protéines constituant un réseau d’interactions, impliquées dans la plasticité synaptique, et hautement exprimées dans le cerveau et ses tissus associés. De plus, l’analyse des sentiers de signalisation pour les gènes identifiés (P = 0,03) a révélé une induction de mécanismes de Potentialisation à Long Terme. Les variations des traits étudiés seraient en grande partie liées au sexe et au statut d’hypertension. Pour la consommation de tabac, nous avons noté que le degré et le sens des corrélations avec l’obésité, les traits hémodynamiques et le stress sont spécifiques au sexe et à la pression artérielle. Par exemple, si des variations ont été détectées entre les hommes fumeurs et non-fumeurs (anciens et jamais), aucune différence n’a été observée chez les femmes. Nous avons aussi identifié de nombreux traits reliés à l’obésité dont la corrélation avec la consommation de tabac apparaît essentiellement plus liée à des facteurs génétiques qu’au fait de fumer en lui-même. Pour le sexe et l’hypertension, des différences dans l’héritabilité de nombreux traits ont également été observées. En effet, des analyses génétiques sur des sous-groupes spécifiques ont révélé des gènes additionnels partageant des fonctions synaptiques : CAMK4, CNTN5, DNM3, KCNAB1 (spécifique à l’hypertension), CNTN4, DNM3, FHIT, ITPR1 and NRXN3 (spécifique au sexe). Ces gènes codent pour des protéines interagissant avec les protéines de gènes détectés dans l’analyse générale. De plus, pour les gènes des sous-groupes, les résultats des analyses des sentiers de signalisation et des profils d’expression des gènes ont montré des caractéristiques similaires à celles de l’analyse générale. La convergence substantielle entre les déterminants génétiques des substances psychoactives, du stress, de l’obésité et des traits hémodynamiques soutiennent la notion selon laquelle les variations génétiques des voies de plasticité synaptique constitueraient une interface commune avec les différences génétiques liées au sexe et à l’hypertension. Nous pensons, également, que la plasticité synaptique interviendrait dans de nombreux phénotypes complexes influencés par le mode de vie. En définitive, ces résultats indiquent que des approches basées sur des sous-groupes et des réseaux amélioreraient la compréhension de la nature polygénique des phénotypes complexes, et des processus moléculaires communs qui les définissent.

Évidence génétique de populations sympatriques de la perchaude (Perca flavescens) dans le Lac Saint-Pierre

La perchaude (Perca flavescens) constitue une ressource socioéconomique de grande importance dans le Lac Saint-Pierre (Québec, Canada). Bien que ce lac fluvial soit désigné réserve de la biosphère par l’UNESCO, le statut de la perchaude est préoccupant. Afin de permettre à l’espèce de persister en fonction des diverses pressions anthropiques, il est important de comprendre sa dynamique populationnelle et les mécanismes qui en sont responsables. La perchaude est connue pour sa philopatrie ; le fait de toujours se reproduire à son site de naissance peut entraîner la subdivision d’une espèce en de multiples populations, où chacune pourra être pourvue d’adaptations locales distinctes. Il est possible d’étudier ces processus à l’aide des signaux génétiques associés à la reproduction des individus. Toutefois, une faible différentiation génétique entre les populations du Lac Saint-Pierre est envisagée en raison de la colonisation récente du système (moins de 8000 ans). L’objectif de cette étude est de déterminer s’il existe plusieurs populations de perchaude dans le Lac Saint-Pierre. Les simulations réalisées ont révélé que l’utilisation de marqueurs AFLP (amplified fragment length polymorphism), permettant une analyse globale du génome, affiche une meilleure détection de la différentiation des populations que celle des marqueurs microsatellites. Afin d’associer les individus à leur site de naissance, la méthode d’AFLP et des microsatellites ont été utilisées sur des larves capturées suite à l’éclosion des oeufs. Trois analyses distinctes d’AFLP ont indiqué une corrélation entre la composition génétique des individus et des sites géographiques, confirmant ainsi la présence de plusieurs populations sympatriques dans le Lac Saint-Pierre, découlant vraisemblablement de la philopatrie de l’espèce. L’absence de différentiation génétique relatée par les marqueurs microsatellites vient confirmer l’importance du choix des marqueurs génétiques. Bien que la différentiation génétique observée soit relativement faible, la gestion de la perchaude devrait tenir compte de la dynamique des populations distinctes dans ce système.

Caractérisation génétique d’une forme d’ataxie tardive.

Les ataxies forment un groupe de maladies neurodégénératives qui sont caractérisées par un manque de coordination des mouvements volontaires. Mes travaux ont porté sur une forme d'ataxie à début tardif (LOCA), après l’âge de 50 ans. Les principales caractéristiques cliniques sont: atrophie cérébelleuse à l’IRM (88%), dysarthrie (81%), atrophie du lobe frontal (50%) et nystagmus (52%). La ségrégation dans les familles de cette ataxie est en faveur d’une transmission récessive. Afin d'identifier le gène responsable de LOCA, nous avons recruté 38 patients affectés d'une forme tardive d'ataxie, issus du SLSJ, des Cantons de l’Est ou d’autres régions du Québec. Un premier criblage du génome a été effectué avec des marqueurs microsatellites sur une famille clé. Une analyse de liaison paramétrique nous a suggéré une liaison au chromosome 13 (4.4Mb). Une recherche d’un haplotype partagé entre 17 familles LOCA a diminué la taille de l'intervalle candidat à 1.6Mb, mais l’haplotype s’est avéré fréquent dans la population canadienne-française. Un second criblage du génome avec des marqueurs SNP nous a permis d’évaluer par cartographie d’homozygotie la possibilité qu’une mutation fondatrice partagée dans des sous-groupes de malades. Plusieurs stratégies d'analyse ont été effectuées, entre autre par regroupement régional. Aucun loci candidats ne fut identifié avec confiance. Nous avons donc combiné les données de génotypage avec le séquençage exomique afin d'identifier le gène responsable. L'analyse de six individus atteints nous a permis d'obtenir une liste de variants rare contenant quatre gènes potentiels. Cette analyse doit se poursuivre pour identifier le gène responsable de LOCA.

Variations dans la réponse de la diversité génétique de populations de couleuvres insulaires faisant face à la perte d’habitat

Projet de recherche réalisé avec Bernard Angers comme directeur de maîtrise, Denis Réale en tant que co-directeur et grâce à la collaboration active d'Emmanuel Milot.