Population genetics of Atlantic bluefin tuna,Thunnus thynnus (Linnaeus, 1758), in the Mediterranean: implications for its conservation management

The movement of Atlantic bluefin tuna (ABFT) across international boundaries necessitates traceability strategies that would provide more accurate information needed for stock assessment. The Mediterranean Sea is one of the main contributors to ABFT reproduction and global population genetic diversity. In the present study this genetic variability was investigated using 193 samples of adult bluefin tuna from Spain, Turkey and Malta – a longitudinal distance of 3400 km. Analysed were 13 microsatellite loci (eight of which were newly‐tested) as genetic markers for the population study. Allele richness measured per locus and sampling location varied from 1.89 to 8.88, taking into account rarefaction. ABFT private alleles were detected in each of the three sampling sites. No significant spatial genetic divergence was found between pairs at the studied locations (F_ST values <0.0001; P‐values >0.05). Bayesian clustering analysis corroborated a single and panmictic ABFT population in the Mediterranean Sea. Statistical power analyses indicated a high probability of detecting genetic differentiation and population structure with the sample size and microsatellites used, even at an F_ST value of 0.005. From the results it may be postulated that migrating ABFT during the spawning season are allowing gene flow within the Mediterranean Sea. The complex interplay of movements, including plasticity in the selection of spawning sites with increasing age and environmental conditions, require multiple and new fisheries monitoring and management techniques in order to target the ABFT long‐term conservation effectively.     

Population genetics of Mediterranean and Saharan olives: geographic patterns of differentiation and evidence for early generations of admixture

Background and Aims

The olive (Olea europaea subsp. europaea) was domesticated in the Mediterranean area but its wild relatives are distributed over three continents, from the Mediterranean basin to South Africa and south-western Asia. Recent studies suggested that this crop originated in the Levant while a secondary diversification occurred in most westward areas. A possible contribution of the Saharan subspecies (subsp. laperrinei) has been highlighted, but the data available were too limited to draw definite conclusions. Here, patterns of genetic differentiation in the Mediterranean and Saharan olives are analysed to test for recent admixture between these taxa.

Methods

Nuclear microsatellite and plastid DNA (ptDNA) data were compiled from previous studies and completed for a sample of 470 cultivars, 390 wild Mediterranean trees and 270 Saharan olives. A network was reconstructed for the ptDNA haplotypes, while a Bayesian clustering method was applied to identify the main gene pools in the data set and then simulate and test for early generations of admixture between Mediterranean and Saharan olives.

Key Results

Four lineages of ptDNA haplotypes are recognized: three from the Mediterranean basin and one from the Sahara. Only one haplotype, primarily distributed in the Sahara, is shared between laperrinei and europaea. This haplotype is detected once in ‘Dhokar’, a cultivar from the Maghreb. Nuclear microsatellites show geographic patterns of genetic differentiation in the Mediterranean olive that reflect the primary origins of cultivars in the Levant, and indicate a high genetic differentiation between europaea and laperrinei. No first-generation hybrid between europaea and laperrinei is detected, but recent, reciprocal admixture between Mediterranean and Saharan subspecies is found in a few accessions, including ‘Dhokar’.

Conclusions

This study reports for the first time admixture between Mediterranean and Saharan olives. Although its contribution remains limited, Laperrine''s olive has been involved in the diversification of cultivated olives.     

Ecological Niche Models and Coalescent Analysis of Gene Flow Support Recent Allopatric Isolation of Parasitoid Wasp Populations in the Mediterranean

Background

The integration of multiple complementary approaches is a powerful way to understand the processes of diversification and speciation. The parasitoid wasp Aphidius transcaspicus Telenga (Hymenoptera: Braconidae) is a parasitoid of Hyalopterus aphids across a wide geographic range. This species shows a remarkable degree of genetic structure among western, central, and eastern Mediterranean population clusters. In this paper we attempt to better characterize this genetic structure.

Methodology/Principal Findings

We use a Bayesian coalescent analysis of gene flow under the Isolation with Migration model using mitochondrial and microsatellite markers together with climate-based ecological niche models to better understand the genetic structure of A. transcaspicus in the Mediterranean. The coalescent analysis revealed low levels of migration among western and eastern Mediterranean populations (Nm<1) that were not statistically distinguishable from zero. Niche models showed that localities within population clusters each occupy areas of continuously high environmental suitability, but are separated from each other by large regions of completely unsuitable habitat that could limit dispersal. Overall, environmental characteristics were similar among the population clusters, though significant differences did emerge.

Conclusions/Significance

These results support contemporary allopatric isolation of Mediterranean populations of A. transcaspicus, which together with previous analyses indicating partial behaviorally mediated reproductive isolation, suggest that the early stages of cryptic speciation may be in progress.     

Genetic and Metabolite Diversity of Sardinian Populations of Helichrysum italicum

Background

Helichrysum italicum (Asteraceae) is a small shrub endemic to the Mediterranean Basin, growing in fragmented and diverse habitats. The species has attracted attention due to its secondary metabolite content, but little effort has as yet been dedicated to assessing the genetic and metabolite diversity present in these populations. Here, we describe the diversity of 50 H. italicum populations collected from a range of habitats in Sardinia.

Methods

H. italicum plants were AFLP fingerprinted and the composition of their leaf essential oil characterized by GC-MS. The relationships between the genetic structure of the populations, soil, habitat and climatic variables and the essential oil chemotypes present were evaluated using Bayesian clustering, contingency analyses and AMOVA.

Key results

The Sardinian germplasm could be partitioned into two AFLP-based clades. Populations collected from the southwestern region constituted a homogeneous group which remained virtually intact even at high levels of K. The second, much larger clade was more diverse. A positive correlation between genetic diversity and elevation suggested the action of natural purifying selection. Four main classes of compounds were identified among the essential oils, namely monoterpenes, oxygenated monoterpenes, sesquiterpenes and oxygenated sesquiterpenes. Oxygenated monoterpene levels were significantly correlated with the AFLP-based clade structure, suggesting a correspondence between gene pool and chemical diversity.

Conclusions

The results suggest an association between chemotype, genetic diversity and collection location which is relevant for the planning of future collections aimed at identifying valuable sources of essential oil.     

Spatial genetic structure of Aquilegia taxa endemic to the island of Sardinia

Background and Aims

The Mediterranean Basin is one of the most important regions for the Earth''s plant biodiversity; however, the scarcity of studies on fine scale patterns of genetic variation in this region is striking. Here, an assessment is made of the spatial genetic structure of all known locations of the three Sardinian endemic species of Aquilegia in order to determine the relative roles of gene flow and genetic drift as underlying evolutionary forces canalizing the divergence of Sardinian Aquilegia taxa, and to see if the spatial genetic structure found fits the current taxonomic differentiation of these taxa.

Methods

DNA from 89 individuals from all known locations of Aquilegia across Sardinia was analysed by means of amplified fragment length polymorphism (AFLP) markers. Both principal co-ordinates analysis (PCoA) and Bayesian clustering analyses were used to determine the spatial genetic structure irrespective of any taxonomic affiliation. Historical effects of gene flow and genetic drift were assessed by checking for the existence of isolation-by-distance patterns.

Key Results

STRUCTURE and PCoA analyses revealed a pattern of genetic variation geographically structured into four spatial genetic groups. No migration–drift equilibrium was detected for Aquilegia in Sardinia, when analysed either as a whole or in individual groups. The scenario approached a Case III pattern sensu Hutchinson and Templeton, which is associated with extreme isolation conditions where genetic drift has historically played a dominant role over gene flow.

Conclusions

The pattern of genetic variation of Sardinian taxa of Aquilegia indicates that genetic drift has been historically more influential than gene flow on population structure of Sardinian species of Aquilegia. Limited seed dispersal and divergent selection imposed by habitat conditions have been probably the main causes reinforcing post-Pleistocene geographical isolation of Aquilegia populations. The spatial genetic structure found here is not fully compatible with current taxonomic affiliations of Sardinian Aquilegia taxa. This is probably a consequence of the uncoupling between morphological and genetic patterns of differentiation frequently found in recently radiated taxa.     

Genetic consequences of anagenetic speciation in Acer okamotoanum (Sapindaceae) on Ullung Island, Korea

Background and Aim

Anagenesis (also known as phyletic speciation) is an important process of speciation in endemic species of oceanic islands. We investigated genetic variation in Acer okamotoanum, an anagenetically derived species endemic to Ullung Island, South Korea, to infer genetic consequences of anagenesis in comparison with other groups that have undergone cladogenesis (and adaptive radiation).

Methods

We examined genetic variation based on eight polymorphic microsatellite markers from 145 individuals of A. okamotoanum and 134 individuals of its putative progenitor A. mono. We employed standard population genetic analyses, clustering analyses, Bayesian clustering analyses in STRUCTURE and bottleneck analyses.

Key Results

Based on both the Neighbor–Joining tree and Bayesian clustering analyses, clear genetic distinctions were found between the two species. Genetic diversity in terms of allelic richness and heterozygosity shows slightly lower levels in A. okamotoanum in comparison with A. mono. Bayesian clustering analyses showed a relatively high F-value in the cluster of A. okamotoanum, suggesting a strong episode of genetic drift during colonization and speciation. There was no clear evidence of a bottleneck based on allelic frequency distribution and excess of observed heterozygotes, but the M-ratio indicated a historical bottleneck in several populations of A. okamotoanum. No geographical genetic structure within the island was found, and the genetic variation among populations of A. okamotoanum was quite low.

Conclusions

We hypothesized that genetic consequences of oceanic-endemic plants derived via anagenesis would be quite different from those derived via cladogenesis. Populations of A. okamotoanum form a cluster and are clearly differentiated from A. mono, which suggests a single origin for the anagenetically derived island endemic. No pattern of geographical differentiation of populations occurs in A. okamotoanum, which supports the concept of initial founder populations diverging through time by accumulation of mutations in a relatively uniform environment without further specific differentiation.     

A Single Locus Is Responsible for Salinity Tolerance in a Chinese Landrace Barley (Hordeum vulgare L.)

Introduction

Salinity and waterlogging are two major abiotic stresses severely limiting barley production. The lack of a reliable screening method makes it very hard to improve the tolerance through breeding programs.

Methods

This work used 188 DH lines from a cross between a Chinese landrace variety, TX9425 (waterlogging and salinity tolerant), and a Japanese malting barley, Naso Nijo (waterlogging and salinity sensitive), to identify QTLs associated with the tolerance.

Results

Four QTLs were found for waterlogging tolerance. The salinity tolerance was evaluated with both a hydroponic system and in potting mixture. In the trial with potting mixture, only one major QTL was identified to associate with salinity tolerance. This QTL explained nearly 50% of the phenotypic variation, which makes it possible for further fine mapping and cloning of the gene. This QTL was also identified in the hydroponic experiment for different salt-related traits. The position of this QTL was located at a similar position to one of the major QTLs for waterlogging tolerance, indicating the possibility of similar mechanisms controlling both waterlogging and salinity tolerance.

Conclusion

The markers associated with the QTL provided a unique opportunity in breeding programs for selection of salinity and waterlogging tolerance.     

The genetic population structure of Thunnus thynnus (Linnaeus, 1758) in the Mediterranean Sea,a controversial issue

The Atlantic bluefin tuna (ABFT), Thunnus thynnus (Linnaeus, 1758), is an important commercial species managed as two different stocks, western and eastern Atlantic, with their spawning grounds in the Gulf of Mexico and in the Mediterranean Sea, respectively. The eastern Atlantic stock has been overexploited in the last decades, leading to the application of specific management measures introduced by the International Commission for the Atlantic Tuna (ICCAT). A clear understanding of the genetic structure of ABFT Mediterranean population should be pursued in order to support management decisions. To date the genetic studies on the Mediterranean ABFT, carried out with different molecular markers and sampling procedures, have produced unclear results. Here, we analysed ABFT samples from central and western Mediterranean Sea with mitochondrial sequences and 11 microsatellite loci to investigate, among the others, the area of the Strait of Messina, where environmental conditions seem to support a resident population of ABFT. Furthermore, genetic analyses of mitochondrial sequences were carried out including nucleotide sequences of Adriatic ABFT wild larvae retrieved from GenBank. Among the investigated areas a genetic differentiation was detected between the Strait of Messina and the Tyrrhenian Sea with microsatellite loci according to the exact G test, but not to the Bayesian analyses carried out with STRUCTURE. The analyses with mitochondrial sequences do not reveal any differentiation among sampled areas, however, a highly significant genetic divergence was observed between the Adriatic mitochondrial sequences retrieved from GenBank and the central‐western Mediterranean sequences obtained in the present work. Our results provide some evidence of population structure of Mediterranean ABFT adding pieces to a still unclear picture.     

A Bayesian Outlier Criterion to Detect SNPs under Selection in Large Data Sets

Background

The recent advent of high-throughput SNP genotyping technologies has opened new avenues of research for population genetics. In particular, a growing interest in the identification of footprints of selection, based on genome scans for adaptive differentiation, has emerged.

Methodology/Principal Findings

The purpose of this study is to develop an efficient model-based approach to perform Bayesian exploratory analyses for adaptive differentiation in very large SNP data sets. The basic idea is to start with a very simple model for neutral loci that is easy to implement under a Bayesian framework and to identify selected loci as outliers via Posterior Predictive P-values (PPP-values). Applications of this strategy are considered using two different statistical models. The first one was initially interpreted in the context of populations evolving respectively under pure genetic drift from a common ancestral population while the second one relies on populations under migration-drift equilibrium. Robustness and power of the two resulting Bayesian model-based approaches to detect SNP under selection are further evaluated through extensive simulations. An application to a cattle data set is also provided.

Conclusions/Significance

The procedure described turns out to be much faster than former Bayesian approaches and also reasonably efficient especially to detect loci under positive selection.     

Using QTL mapping to investigate the relationships between abiotic stress tolerance (drought and salinity) and agronomic and physiological traits

Background

Drought and salinity are two major abiotic stresses that severely limit barley production worldwide. Physiological and genetic complexity of these tolerance traits has significantly slowed the progress of developing stress-tolerant cultivars. Marker-assisted selection (MAS) may potentially overcome this problem. In the current research, seventy two double haploid (DH) lines from a cross between TX9425 (a Chinese landrace variety with superior drought and salinity tolerance) and a sensitive variety, Franklin were used to identify quantitative trait loci (QTL) for drought and salinity tolerance, based on a range of developmental and physiological traits.

Results

Two QTL for drought tolerance (leaf wilting under drought stress) and one QTL for salinity tolerance (plant survival under salt stress) were identified from this population. The QTL on 2H for drought tolerance determined 42% of phenotypic variation, based on three independent experiments. This QTL was closely linked with a gene controlling ear emergency. The QTL on 5H for drought tolerance was less affected by agronomic traits and can be effectively used in breeding programs. A candidate gene for this QTL on 5H was identified based on the draft barley genome sequence. The QTL for proline accumulation, under both drought and salinity stresses, were located on different positions to those for drought and salinity tolerance, indicating no relationship with plant tolerance to either of these stresses.

Conclusions

Using QTL mapping, the relationships between QTL for agronomic and physiological traits and plant drought and salinity tolerance were studied. A new QTL for drought tolerance which was not linked to any of the studied traits was identified. This QTL can be effectively used in breeding programs. It was also shown that proline accumulation under stresses was not necessarily linked with drought or salinity tolerance based on methods of phenotyping used in this experiment. The use of proline content in breeding programs can also be limited by the accuracy of phenotyping.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1243-8) contains supplementary material, which is available to authorized users.     

Predictions of native American population structure using linguistic covariates in a hidden regression framework

Background

The mainland of the Americas is home to a remarkable diversity of languages, and the relationships between genes and languages have attracted considerable attention in the past. Here we investigate to which extent geography and languages can predict the genetic structure of Native American populations.

Methodology/Principal Findings

Our approach is based on a Bayesian latent cluster regression model in which cluster membership is explained by geographic and linguistic covariates. After correcting for geographic effects, we find that the inclusion of linguistic information improves the prediction of individual membership to genetic clusters. We further compare the predictive power of Greenberg''s and The Ethnologue classifications of Amerindian languages. We report that The Ethnologue classification provides a better genetic proxy than Greenberg''s classification at the stock and at the group levels. Although high predictive values can be achieved from The Ethnologue classification, we nevertheless emphasize that Choco, Chibchan and Tupi linguistic families do not exhibit a univocal correspondence with genetic clusters.

Conclusions/Significance

The Bayesian latent class regression model described here is efficient at predicting population genetic structure using geographic and linguistic information in Native American populations.     

Analysis of Whole-Brain Resting-State fMRI Data Using Hierarchical Clustering Approach

Background

Previous studies using hierarchical clustering approach to analyze resting-state fMRI data were limited to a few slices or regions-of-interest (ROIs) after substantial data reduction.

Purpose

To develop a framework that can perform voxel-wise hierarchical clustering of whole-brain resting-state fMRI data from a group of subjects.

Materials and Methods

Resting-state fMRI measurements were conducted for 86 adult subjects using a single-shot echo-planar imaging (EPI) technique. After pre-processing and co-registration to a standard template, pair-wise cross-correlation coefficients (CC) were calculated for all voxels inside the brain and translated into absolute Pearson''s distances after imposing a threshold CC≥0.3. The group averages of the Pearson''s distances were then used to perform hierarchical clustering with the developed framework, which entails gray matter masking and an iterative scheme to analyze the dendrogram.

Results

With the hierarchical clustering framework, we identified most of the functional connectivity networks reported previously in the literature, such as the motor, sensory, visual, memory, and the default-mode functional networks (DMN). Furthermore, the DMN and visual system were split into their corresponding hierarchical sub-networks.

Conclusion

It is feasible to use the proposed hierarchical clustering scheme for voxel-wise analysis of whole-brain resting-state fMRI data. The hierarchical clustering result not only confirmed generally the finding in functional connectivity networks identified previously using other data processing techniques, such as ICA, but also revealed directly the hierarchical structure within the functional connectivity networks.     

Unravelling genetics at the top: mountain islands or isolated belts?

Background and Aims

In mountain plant populations, local adaptation has been described as one of the main responses to climate warming, allowing plants to persist under stressful conditions. This is especially the case for marginal populations at their lowest elevation, as they are highly vulnerable. Adequate levels of genetic diversity are required for selection to take place, while high levels of altitudinal gene flow are seen as a major limiting factor potentially precluding local adaptation processes. Thus, a compromise between genetic diversity and gene flow seems necessary to guarantee persistence under oncoming conditions. It is therefore critical to determine if gene flow occurs preferentially between mountains at similar altitudinal belts, promoting local adaptation at the lowest populations, or conversely along altitude within each mountain.

Methods

Microsatellite markers were used to unravel genetic diversity and population structure, inbreeding and gene flow of populations at two nearby altitudinal gradients of Silene ciliata, a Mediterranean high-mountain cushion plant.

Key Results

Genetic diversity and inbreeding coefficients were similar in all populations. Substantial gene flow was found both along altitudinal gradients and horizontally within each elevation belt, although greater values were obtained along altitudinal gradients. Gene flow may be responsible for the homogeneous levels of genetic diversity found among populations. Bayesian cluster analyses also suggested that shifts along altitudinal gradients are the most plausible scenario.

Conclusions

Past population shifts associated with glaciations and interglacial periods in temperate mountains may partially explain current distributions of genetic diversity and population structure. In spite of the predominance of gene flow along the altitudinal gradients, local genetic differentiation of one of the lower populations together with the detection of one outlier locus might support the existence of different selection forces at low altitudes.     

Multigene molecular systematics confirm species status of morphologically convergent Pagurus hermit crabs

Introduction

In spite of contemporary morphological taxonomy appraisals, apparent high morphological similarity raises uncertainty about the species status of certain Pagurus hermit crabs. This is exemplified between two European species, Pagurus excavatus (Herbst, 1791) and Pagurus alatus (Fabricius 1775), whose species status is still difficult to resolve using morphological criteria alone.

Methodology/Principal Findings

To address such ambiguities, we used combinations of Maximum Likelihood (ML) and Bayesian Inference (BI) methods to delineate species boundaries of P. alatus and P. excavatus and formulate an intermediate Pagurus phylogenetic hypothesis, based upon single and concatenated mitochondrial (cytochrome oxidase I [COI]) and nuclear (16S and 28s ribosomal RNA) gene partitions. The molecular data supported the species status of P. excavatus and P. alatus and also clearly resolved two divergent clades within hermit crabs from the Northeast Atlantic Ocean and the Mediterranean Sea.

Conclusions/Significance

Despite the abundance and prominent ecological role of hermit crabs, Pagurus, in North East Atlantic Ocean and Mediterranean Sea ecosystems, many important aspects of their taxonomy, biology, systematics and evolution remain poorly explored. The topologies presented here should be regarded as hypotheses that can be incorporated into the robust and integrated understanding of the systematic relationships within and between species of the genus Pagurus inhabiting the Northeast Atlantic Ocean and the Mediterranean Sea.     

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 between- and 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.     

Origin of the invasive Arundo donax (Poaceae): a trans-Asian expedition in herbaria

Background and Aims

The hypothesis of an ancient introduction, i.e. archaeophyte origin, is one of the most challenging questions in phylogeography. Arundo donax (Poaceae) is currently considered to be one of the worst invasive species globally, but it has also been widely utilzed by man across Eurasia for millennia. Despite a lack of phylogenetic data, recent literature has often speculated on its introduction to the Mediterranean region.

Methods

This study tests the hypothesis of its ancient introduction from Asia to the Mediterranean by using plastid DNA sequencing and morphometric analysis on 127 herbarium specimens collected across sub-tropical Eurasia. In addition, a bioclimatic species distribution model calibrated on 1221 Mediterranean localities was used to identify similar ecological niches in Asia.

Key Results

Despite analysis of several plastid DNA hypervariable sites and the identification of 13 haplotypes, A. donax was represented by a single haplotype from the Mediterranean to the Middle East. This haplotype is shared with invasive samples worldwide, and its nearest phylogenetic relatives are located in the Middle East. Morphometric data characterized this invasive clone by a robust morphotype distinguishable from all other Asian samples. The ecological niche modelling designated the southern Caspian Sea, southern Iran and the Indus Valley as the most suitable regions of origin in Asia for the invasive clone of A. donax.

Conclusions

Using an integrative approach, an ancient dispersion of this robust, polyploid and non-fruiting clone is hypothesized from the Middle East to the west, leading to its invasion throughout the Mediterranean Basin.     

Geographical structuring of genetic diversity across the whole distribution range of Narcissus longispathus, a habitat-specialist, Mediterranean narrow endemic

Background and Aims

High mountain ranges of the Mediterranean Basin harbour a large number of narrowly endemic plants. In this study an investigation is made of the levels and partitioning of genetic diversity in Narcissus longispathus, a narrow endemic of south-eastern Spanish mountains characterized by a naturally fragmented distribution due to extreme specialization on a rare habitat type. By using dense sampling of populations across the species'' whole geographical range, genetic structuring at different geographical scales is also examined.

Methods

Using horizontal starch-gel electrophoresis, allozyme variability was screened at 19 loci for a total of 858 individuals from 27 populations. The data were analysed by means of standard statistical approaches in order to estimate gene diversity and the genetic structure of the populations.

Key Results

Narcissus longispathus displayed high levels of genetic diversity and extensive diversification among populations. At the species level, the percentage of polymorphic loci was 68 %, with average values of 2·1, 0·11 and 0·14 for the number of alleles per locus, observed heterozygosity and expected heterozygosity, respectively. Southern and more isolated populations tended to have less genetic variability than northern and less-isolated populations. A strong spatial patterning of genetic diversity was found at the various spatial scales. Gene flow/drift equilibrium occurred over distances <4 km. Beyond that distance divergence was relatively more influenced by drift. The populations studied seem to derive from three panmictic units or ‘gene pools’, with levels of admixture being greatest in the central and south-eastern portions of the species'' range.

Conclusions

In addition to documenting a case of high genetic diversity in a narrow endemic plant with naturally fragmented populations, the results emphasize the need for dense population sampling and examination of different geographical scales for understanding population genetic structure in habitat specialists restricted to ecological islands.Key words: Allozymes, genetic diversity, geographical scale, habitat isolation, Narcissus longispathus, Mediterranean endemism, mountain range, natural fragmented distribution     

Genetic diversity and networks of exchange: a combined approach to assess intra-breed diversity

Background

Cryopreservation of three endangered Belgian sheep breeds required to characterize their intra-breed genetic diversity. It is assumed that the genetic structure of a livestock breed depends mostly on gene flow due to exchanges between herds. To quantify this relation, molecular data and analyses of the exchanges were combined for three endangered Belgian breeds.

Methods

For each breed, between 91 and 225 sheep were genotyped with 19 microsatellites. Genetic differentiations between breeds and among herds within a breed were evaluated and the genetic structure of the breeds was described using Bayesian clustering (Structure). Exchanges of animals between 20, 46 and 95 herds according to breed were identified via semi-directed interviews and were analyzed using the concepts of the network theory to calculate average degrees and shortest path lengths between herds. Correlation between the Reynolds’ genetic distances and the shortest path lengths between each pair of herds was assessed by a Mantel test approach.

Results

Genetic differentiation between breeds was high (0.16). Overall Fst values among herds were high in each breed (0.17, 0.11 and 0.10). Use of the Bayesian approach made it possible to identify genetic groups of herds within a breed. Significant correlations between the shortest path lengths and the Reynolds’ genetic distances were found in each breed (0.87, 0.33 and 0.41), which demonstrate the influence of exchanges between herds on the genetic diversity. Correlation differences between breeds could be explained by differences in the average degree of the animal exchange networks, which is a measure of the number of exchanges per herd. The two breeds with the highest average degree showed the lowest correlation. Information from the exchange networks was used to assign individuals to the genetic groups when molecular information was incomplete or missing to identify donors for a cryobank.

Conclusions

A fine-scale picture of the population genetic structure at the herd level was obtained for the three breeds. Network analysis made it possible to highlight the influence of exchanges on genetic structure and to complete or replace molecular information in establishing a conservation program.     

Bayesian neural networks for detecting epistasis in genetic association studies

Background

Discovering causal genetic variants from large genetic association studies poses many difficult challenges. Assessing which genetic markers are involved in determining trait status is a computationally demanding task, especially in the presence of gene-gene interactions.

Results

A non-parametric Bayesian approach in the form of a Bayesian neural network is proposed for use in analyzing genetic association studies. Demonstrations on synthetic and real data reveal they are able to efficiently and accurately determine which variants are involved in determining case-control status. By using graphics processing units (GPUs) the time needed to build these models is decreased by several orders of magnitude. In comparison with commonly used approaches for detecting interactions, Bayesian neural networks perform very well across a broad spectrum of possible genetic relationships.

Conclusions

The proposed framework is shown to be a powerful method for detecting causal SNPs while being computationally efficient enough to handle large datasets.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0368-0) contains supplementary material, which is available to authorized users.     

New insights into the Lake Chad Basin population structure revealed by high-throughput genotyping of mitochondrial DNA coding SNPs

Background

Located in the Sudan belt, the Chad Basin forms a remarkable ecosystem, where several unique agricultural and pastoral techniques have been developed. Both from an archaeological and a genetic point of view, this region has been interpreted to be the center of a bidirectional corridor connecting West and East Africa, as well as a meeting point for populations coming from North Africa through the Saharan desert.

Methodology/Principal Findings

Samples from twelve ethnic groups from the Chad Basin (n = 542) have been high-throughput genotyped for 230 coding region mitochondrial DNA (mtDNA) Single Nucleotide Polymorphisms (mtSNPs) using Matrix-Assisted Laser Desorption/Ionization Time-Of-Flight (MALDI-TOF) mass spectrometry. This set of mtSNPs allowed for much better phylogenetic resolution than previous studies of this geographic region, enabling new insights into its population history. Notable haplogroup (hg) heterogeneity has been observed in the Chad Basin mirroring the different demographic histories of these ethnic groups. As estimated using a Bayesian framework, nomadic populations showed negative growth which was not always correlated to their estimated effective population sizes. Nomads also showed lower diversity values than sedentary groups.

Conclusions/Significance

Compared to sedentary population, nomads showed signals of stronger genetic drift occurring in their ancestral populations. These populations, however, retained more haplotype diversity in their hypervariable segments I (HVS-I), but not their mtSNPs, suggesting a more ancestral ethnogenesis. Whereas the nomadic population showed a higher Mediterranean influence signaled mainly by sub-lineages of M1, R0, U6, and U5, the other populations showed a more consistent sub-Saharan pattern. Although lifestyle may have an influence on diversity patterns and hg composition, analysis of molecular variance has not identified these differences. The present study indicates that analysis of mtSNPs at high resolution could be a fast and extensive approach for screening variation in population studies where labor-intensive techniques such as entire genome sequencing remain unfeasible.