A note on the rationale for estimating genealogical coancestry from molecular markers

Background

Genetic relatedness or similarity between individuals is a key concept in population, quantitative and conservation genetics. When the pedigree of a population is available and assuming a founder population from which the genealogical records start, genetic relatedness between individuals can be estimated by the coancestry coefficient. If pedigree data is lacking or incomplete, estimation of the genetic similarity between individuals relies on molecular markers, using either molecular coancestry or molecular covariance. Some relationships between genealogical and molecular coancestries and covariances have already been described in the literature.

Methods

We show how the expected values of the empirical measures of similarity based on molecular marker data are functions of the genealogical coancestry. From these formulas, it is easy to derive estimators of genealogical coancestry from molecular data. We include variation of allelic frequencies in the estimators.

Results

The estimators are illustrated with simulated examples and with a real dataset from dairy cattle. In general, estimators are accurate and only slightly biased. From the real data set, estimators based on covariances are more compatible with genealogical coancestries than those based on molecular coancestries. A frequently used estimator based on the average of estimated coancestries produced inflated coancestries and numerical instability. The consequences of unknown gene frequencies in the founder population are briefly discussed, along with alternatives to overcome this limitation.

Conclusions

Estimators of genealogical coancestry based on molecular data are easy to derive. Estimators based on molecular covariance are more accurate than those based on identity by state. A correction considering the random distribution of allelic frequencies improves accuracy of these estimators, especially for populations with very strong drift.     

Maintaining genetic diversity using molecular coancestry: the effect of marker density and effective population size

Background

The most efficient method to maintain genetic diversity in populations under conservation programmes is to optimize, for each potential parent, the number of offspring left to the next generation by minimizing the global coancestry. Coancestry is usually calculated from genealogical data but molecular markers can be used to replace genealogical coancestry with molecular coancestry. Recent studies showed that optimizing contributions based on coancestry calculated from a large number of SNP markers can maintain higher levels of diversity than optimizing contributions based on genealogical data. In this study, we investigated how SNP density and effective population size impact the use of molecular coancestry to maintain diversity.

Results

At low SNP densities, the genetic diversity maintained using genealogical coancestry for optimization was higher than that maintained using molecular coancestry. The performance of molecular coancestry improved with increasing marker density, and, for the scenarios evaluated, it was as efficient as genealogical coancestry if SNP density reached at least 3 times the effective population size.However, increasing SNP density resulted in reduced returns in terms of maintained diversity. While a benefit of 12% was achieved when marker density increased from 10 to 100 SNP/Morgan, the benefit was only 2% when it increased from 100 to 500 SNP/Morgan.

Conclusions

The marker density of most SNP chips already available for farm animals is sufficient for molecular coancestry to outperform genealogical coancestry in conservation programmes aimed at maintaining genetic diversity. For the purpose of effectively maintaining genetic diversity, a marker density of around 500 SNPs/Morgan can be considered as the most cost effective density when developing SNP chips for new species. Since the costs to develop SNP chips are decreasing, chips with 500 SNPs/Morgan should become available in a short-term horizon for non domestic species.     

Genome-wide insights into population structure and genetic history of tunisian local cattle using the illumina bovinesnp50 beadchip

Background

Tunisian local cattle populations are at risk of extinction as they were massively crossed with imported breeds. Preservation of indigenous livestock populations is important because each of them comprises a unique set of genes resulting from a local environment-driven selection that occurred over hundreds of years. The diversity and genetic structure of Tunisian local cattle populations are poorly understood. However, such information is crucial to the conservation and sustainable use of genetic resources.In addition, comparing the genomic structure of population sets from different parts of the world could help yield insight into their origin and history.In the present study, we provide a detailed assessment of the population structure of the three Tunisian local cattle populations using various methods, and we highlight their origin and history by investigating approximately ~38,000 SNPs in a broad panel of 878 individuals from 37 worldwide cattle breeds representative of African, European and indicine populations.

Results

Our study revealed a low level of divergence and high genetic diversity in Tunisian local cattle reflecting low levels of genetic drift. A Comparison with the worldwide cattle panel pinpointed the admixed origin of the genome of the three Tunisian populations with the two main European and African ancestries. Our results were in agreement with previous historical and archaeological reports about the past gene flow that existed between North African and South European breeds, in particular with Iberian cattle. We also detected a low-level indicine introgression in the three Tunisian populations and we inferred that indicine ancestry was inherited via African ancestors.

Conclusions

Our results represent the first study providing genetic evidence about the origin and history of Tunisian local cattle. The information provided by the fine-scale genetic characterization of our study will enhance the establishment of a national conservation strategy for these populations. These results may enable the identification of genetic variants involved in adaptation to harsh environmental conditions.

Electronic supplementary material

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

Removing exogenous information using pedigree data

Management of certain populations requires the preservation of its pure genetic background. When, for different reasons, undesired alleles are introduced, the original genetic conformation must be recovered. The present study tested, through computer simulations, the power of recovery (the ability for removing the foreign information) from genealogical data. Simulated scenarios comprised different numbers of exogenous individuals taking part of the founder population and different numbers of unmanaged generations before the removal program started. Strategies were based on variables arising from classical pedigree analyses such as founders’ contribution and partial coancestry. The efficiency of the different strategies was measured as the proportion of native genetic information remaining in the population. Consequences on the inbreeding and coancestry levels of the population were also evaluated. Minimisation of the exogenous founders’ contributions was the most powerful method, removing the largest amount of genetic information in just one generation. However, as a side effect, it led to the highest values of inbreeding. Scenarios with a large amount of initial exogenous alleles (i.e. high percentage of non native founders), or many generations of mixing became very difficult to recover, pointing out the importance of being careful about introgression events in populations where these are undesired.     

Effect of non-random mating on genomic and BLUP selection schemes

Background

The risk of long-term unequal contribution of mating pairs to the gene pool is that deleterious recessive genes can be expressed. Such consequences could be alleviated by appropriately designing and optimizing breeding schemes i.e. by improving selection and mating procedures.

Methods

We studied the effect of mating designs, random, minimum coancestry and minimum covariance of ancestral contributions on rate of inbreeding and genetic gain for schemes with different information sources, i.e. sib test or own performance records, different genetic evaluation methods, i.e. BLUP or genomic selection, and different family structures, i.e. factorial or pair-wise.

Results

Results showed that substantial differences in rates of inbreeding due to mating design were present under schemes with a pair-wise family structure, for which minimum coancestry turned out to be more effective to generate lower rates of inbreeding. Specifically, substantial reductions in rates of inbreeding were observed in schemes using sib test records and BLUP evaluation. However, with a factorial family structure, differences in rates of inbreeding due mating designs were minor. Moreover, non-random mating had only a small effect in breeding schemes that used genomic evaluation, regardless of the information source.

Conclusions

It was concluded that minimum coancestry remains an efficient mating design when BLUP is used for genetic evaluation or when the size of the population is small, whereas the effect of non-random mating is smaller in schemes using genomic evaluation.     

Combined detection and introgression of QTL in outbred populations

Background

Detecting a QTL is only the first step in genetic improvement programs. When a QTL with desirable characteristics is found, e.g. in a wild or unimproved population, it may be interesting to introgress the detected QTL into the commercial population. One approach to shorten the time needed for introgression is to combine both QTL identification and introgression, into a single step. This combines the strengths of fine mapping and backcrossing and paves the way for introgression of desirable but unknown QTL into recipient animal and plant lines.

Methods

The method consisting in combining QTL mapping and gene introgression has been extended from inbred to outbred populations in which QTL allele frequencies vary both in recipient and donor lines in different scenarios and for which polygenic effects are included in order to model background genes. The effectiveness of the combined QTL detection and introgression procedure was evaluated by simulation through four backcross generations.

Results

The allele substitution effect is underestimated when the favourable QTL allele is not fixed in the donor line. This underestimation is proportional to the frequency differences of the favourable QTL allele between the lines. In most scenarios, the estimates of the QTL location are unbiased and accurate. The retained donor chromosome segment and linkage drag are similar to expected values from other published studies.

Conclusions

In general, our results show that it is possible to combine QTL detection and introgression even in outbred species. Separating QTL mapping and introgression processes is often thought to be longer and more costly. However, using a combined process saves at least one generation. With respect to the linkage drag and obligatory drag, the results of the combined detection and introgression scheme are very similar to those of traditional introgression schemes.     

Methods to estimate effective population size using pedigree data: Examples in dog,sheep, cattle and horse

Background

Effective population sizes of 140 populations (including 60 dog breeds, 40 sheep breeds, 20 cattle breeds and 20 horse breeds) were computed using pedigree information and six different computation methods. Simple demographical information (number of breeding males and females), variance of progeny size, or evolution of identity by descent probabilities based on coancestry or inbreeding were used as well as identity by descent rate between two successive generations or individual identity by descent rate.

Results

Depending on breed and method, effective population sizes ranged from 15 to 133 056, computation method and interaction between computation method and species showing a significant effect on effective population size (P < 0.0001). On average, methods based on number of breeding males and females and variance of progeny size produced larger values (4425 and 356, respectively), than those based on identity by descent probabilities (average values between 93 and 203). Since breeding practices and genetic substructure within dog breeds increased inbreeding, methods taking into account the evolution of inbreeding produced lower effective population sizes than those taking into account evolution of coancestry. The correlation level between the simplest method (number of breeding males and females, requiring no genealogical information) and the most sophisticated one ranged from 0.44 to 0.60 according to species.

Conclusions

When choosing a method to compute effective population size, particular attention should be paid to the species and the specific genetic structure of the population studied.     

Genome-wide prediction methods for detecting genetic effects of donor chromosome segments in introgression populations

Background

Introgression populations are used to make the genetic variation of unadapted germplasm or wild relatives of crops available for plant breeding. They consist of introgression lines that carry small chromosome segments from an exotic donor in the genetic background of an elite line. The goal of our study was to investigate the detection of favorable donor chromosome segments in introgression lines with statistical methods developed for genome-wide prediction.

Results

Computer simulations showed that genome-wide prediction employing heteroscedastic marker variances had a greater power and a lower false positive rate compared with homoscedastic marker variances when the phenotypic difference between the donor and recipient lines was controlled by few genes. The simulations helped to interpret the analyses of glycosinolate and linolenic acid content in a rapeseed introgression population and plant height in a rye introgression population. These analyses support the superiority of genome-wide prediction approaches that use heteroscedastic marker variances.

Conclusions

We conclude that genome-wide prediction methods in combination with permutation tests can be employed for analysis of introgression populations. They are particularly useful when introgression lines carry several donor segments or when the donor segments of different introgression lines are overlapping.     

Genetic variation and relationships of seven sturgeon species and ten interspecific hybrids

Background

Sturgeon cultivation is important for both industry and aquaculture in China. To date, more than 17 species or strains have been farmed for fillets and caviar production. Crossbreeding among different sturgeon species is frequent and the F2 hybrids are fertile. However, large-scale farming can have negative impacts on wild populations i.e. escape of exotic sturgeons and must be taken into consideration. Escape of exotic sturgeons can cause severe ecological problems, including threatening native sturgeon species once the exotic varieties become established or hybridize with native individuals. However, little is known about their genetic resources and variation.

Methods

Genetic diversity and introgression of seven sturgeon species were analyzed using mitochondrial DNA cytochrome oxidase subunit I (COI) and nine microsatellite markers. This study included 189 individuals from seven sturgeon species and 277 individuals from ten lineages of F2 hybrid strains.

Results

MtDNA COI sequences (632 bp long) were generated from 91 individuals across the 17 sturgeon strains and produced 23 different haplotypes. Haplotype diversity was high (h = 0.915 ± 0.015) and nucleotide diversity was low (π = 0.03680 ± 0.00153) in the seven sturgeon species and ten interspecific hybrids. Phylogenetic analyses resulted in almost identical tree topologies, and different haplotype structures were mainly related with sturgeons of different female parents. Analysis of molecular variance revealed that 81.73% of the genetic variance was due to matrilineal differences, while 9.40% resulted from strain variation. Pairwise Fst values obtained with POLYSAT software, were high among strains and ranged from 0.031 to 0.164. Admixture analysis assigned seven distinct groups and ten genotypes of admixed clusters composed of hybrid strains using STRUCTURE when assuming K = 7.

Conclusions

The interspecific mtDNA gene tree corresponded to the expected taxonomic divisions. These relationships were also supported by the results from the microsatellite analysis and contributed to unambiguously identify seven sturgeon species and ten F2 hybrid strains from sturgeon farms in China. Moreover, we found that introgressive hybridization is pervasive, exists in both purebred and hybrid sturgeons, and may reflect widespread mismanagement in sturgeon breeding in China.     

Genetic evidence of tiger population structure and migration within an isolated and fragmented landscape in Northwest India

Background

Majority of the tiger habitat in Indian subcontinent lies within high human density landscapes and is highly sensitive to surrounding pressures. These forests are unable to sustain healthy tiger populations within a tiger-hostile matrix, despite considerable conservation efforts. Ranthambore Tiger Reserve (RTR) in Northwest India is one such isolated forest which is rapidly losing its links with other tiger territories in the Central Indian landscape. Non-invasive genetic sampling for individual identification is a potent technique to understand the relationships between threatened tiger populations in degraded habitats. This study is an attempt to establish tiger movement across a fragmented landscape between RTR and its neighboring forests, Kuno-Palpur Wildlife Sanctuary (KPWLS) and Madhav National Park (MNP) based on non-invasively obtained genetic data.

Methods

Data from twelve microsatellite loci was used to define population structure and also to identify first generation migrants and admixed individuals in the above forests.

Results

Population structure was consistent with the Central Indian landscape and we could determine significant gene flow between RTR and MNP. We could identify individuals of admixed ancestry in both these forests, as well as first generation migrants from RTR to KPWLS and MNP.

Conclusions

Our results indicate reproductive mixing between animals of RTR and MNP in the recent past and migration of animals even today, despite fragmentation and poaching risk, from RTR towards MNP. Substantial conservation efforts should be made to maintain connectivity between these two subpopulations and also higher protection status should be conferred on Madhav National Park.     

Strategies for implementing genomic selection in family-based aquaculture breeding schemes: double haploid sib test populations

Background

Simulation studies have shown that accuracy and genetic gain are increased in genomic selection schemes compared to traditional aquaculture sib-based schemes. In genomic selection, accuracy of selection can be maximized by increasing the precision of the estimation of SNP effects and by maximizing the relationships between test sibs and candidate sibs. Another means of increasing the accuracy of the estimation of SNP effects is to create individuals in the test population with extreme genotypes. The latter approach was studied here with creation of double haploids and use of non-random mating designs.

Methods

Six alternative breeding schemes were simulated in which the design of the test population was varied: test sibs inherited maternal (Mat), paternal (Pat) or a mixture of maternal and paternal (MatPat) double haploid genomes or test sibs were obtained by maximum coancestry mating (MaxC), minimum coancestry mating (MinC), or random (RAND) mating. Three thousand test sibs and 3000 candidate sibs were genotyped. The test sibs were recorded for a trait that could not be measured on the candidates and were used to estimate SNP effects. Selection was done by truncation on genome-wide estimated breeding values and 100 individuals were selected as parents each generation, equally divided between both sexes.

Results

Results showed a 7 to 19% increase in selection accuracy and a 6 to 22% increase in genetic gain in the MatPat scheme compared to the RAND scheme. These increases were greater with lower heritabilities. Among all other scenarios, i.e. Mat, Pat, MaxC, and MinC, no substantial differences in selection accuracy and genetic gain were observed.

Conclusions

In conclusion, a test population designed with a mixture of paternal and maternal double haploids, i.e. the MatPat scheme, increases substantially the accuracy of selection and genetic gain. This will be particularly interesting for traits that cannot be recorded on the selection candidates and require the use of sib tests, such as disease resistance and meat quality.     

Recent genomic heritage in Scotland

Background

The Generation Scotland Scottish Family Health Study (GS:SFHS) includes 23,960 participants from across Scotland with records for many health-related traits and environmental covariates. Genotypes at ~700 K SNPs are currently available for 10,000 participants. The cohort was designed as a resource for genetic and health related research and the study of complex traits. In this study we developed a suite of analyses to disentangle the genomic differentiation within GS:SFHS individuals to describe and optimise the sample and methods for future analyses.

Results

We combined the genotypic information of GS:SFHS with 1092 individuals from the 1000 Genomes project and estimated their genomic relationships. Then, we performed Principal Component Analyses of the resulting relationships to investigate the genomic origin of different groups. We characterised two groups of individuals: those with a few sparse rare markers in the genome, and those with several large rare haplotypes which might represent relatively recent exogenous ancestors. We identified some individuals with likely Italian ancestry and a group with some potential African/Asian ancestry. An analysis of homozygosity in the GS:SFHS sample revealed a very similar pattern to other European populations. We also identified an individual carrying a chromosome 1 uniparental disomy. We found evidence of local geographic stratification within the population having impact on the genomic structure.

Conclusions

These findings illuminate the history of the Scottish population and have implications for further analyses such as the study of the contributions of common and rare variants to trait heritabilities and the evaluation of genomic and phenotypic prediction of disease.

Electronic supplementary material

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

Genetic differentiation of European grayling (Thymallus thymallus) populations in Serbia,based on mitochondrial and nuclear DNA analyses

Background

The structure and diversity of grayling (Thymallus thymallus) populations have been well studied in most of its native habitat; however the southernmost populations of the Balkan Peninsula remain largely unexplored. The purpose of this study was to assess the genetic diversity of Serbian grayling populations, detect the impact of stocking and provide guidelines for conservation and management.

Methods

Eighty grayling individuals were collected from four rivers (Ibar, Lim, Drina and Rzav). The mitochondrial DNA control region (CR; 595 bp of the 3''end and 74 bp of flanking tRNA) and the ATP6 gene (630 bp fragment) were sequenced for 20 individuals (five from each locality). In addition, all individuals were genotyped with 12 microsatellite loci. The diversity and structure of the populations as well as the recent and ancient population declines were studied using specialized software.

Results

We detected three new haplotypes in the mtDNA CR and four haplotypes in the ATP6 gene of which three had not been described before. Previously, one CR haplotype and two ATP6 gene haplotypes had been identified as allochthonous, originating from Slovenia. Reconstruction of phylogenetic relations placed the remaining two CR haplotypes from the River Danube drainage of Serbia into a new clade, which is related to the previously described sister Slovenian clade. These two clades form a new Balkan clade. Microsatellite marker analysis showed that all four populations are genetically distinct from each other without any sign of intra-population structure, although stocking of the most diverse population (Drina River) was confirmed by mtDNA analysis. Recent and historical population declines of Serbian grayling do not differ from those of other European populations.

Conclusions

Our study shows that (1) the Ibar, Lim and Drina Rivers grayling populations are genetically distinct from populations outside of Serbia and thus should be managed as native populations in spite of some introgression in the Drina River population and (2) the Rzav River population is not appropriate for further stocking activities since it originates from stocked Slovenian grayling. However, the Rzav River population does not represent an immediate danger to other populations because it is physically isolated from these.     

No difference between the sexes in fine-scale spatial genetic structure of roe deer

Background

Data on spatial genetic patterns may provide information about the ecological and behavioural mechanisms underlying population structure. Indeed, social organization and dispersal patterns of species may be reflected by the pattern of genetic structure within a population.

Methodology/Principal Findings

We investigated the fine-scale spatial genetic structure of a roe deer (Capreolus capreolus) population in Trois-Fontaines (France) using 12 microsatellite loci. The roe deer is weakly polygynous and highly sedentary, and can form matrilineal clans. We show that relatedness among individuals was negatively correlated with geographic distance, indicating that spatially proximate individuals are also genetically close. More unusually for a large mammalian herbivore, the link between relatedness and distance did not differ between the sexes, which is consistent with the lack of sex-biased dispersal and the weakly polygynous mating system of roe deer.

Conclusions/Significance

Our results contrast with previous reports on highly polygynous species with male-biased dispersal, such as red deer, where local genetic structure was detected in females only. This divergence between species highlights the importance of socio-spatial organization in determining local genetic structure of vertebrate populations.     

European Invasion of North American Pinus strobus at Large and Fine Scales: High Genetic Diversity and Fine-Scale Genetic Clustering over Time in the Adventive Range

Background

North American Pinus strobus is a highly invasive tree species in Central Europe. Using ten polymorphic microsatellite loci we compared various aspects of the large-scale genetic diversity of individuals from 30 sites in the native distribution range with those from 30 sites in the European adventive distribution range. To investigate the ascertained pattern of genetic diversity of this intercontinental comparison further, we surveyed fine-scale genetic diversity patterns and changes over time within four highly invasive populations in the adventive range.

Results

Our data show that at the large scale the genetic diversity found within the relatively small adventive range in Central Europe, surprisingly, equals the diversity found within the sampled area in the native range, which is about thirty times larger. Bayesian assignment grouped individuals into two genetic clusters separating North American native populations from the European, non-native populations, without any strong genetic structure shown over either range. In the case of the fine scale, our comparison of genetic diversity parameters among the localities and age classes yielded no evidence of genetic diversity increase over time. We found that SGS differed across age classes within the populations under study. Old trees in general completely lacked any SGS, which increased over time and reached its maximum in the sapling stage.

Conclusions

Based on (1) the absence of difference in genetic diversity between the native and adventive ranges, together with the lack of structure in the native range, and (2) the lack of any evidence of any temporal increase in genetic diversity at four highly invasive populations in the adventive range, we conclude that population amalgamation probably first happened in the native range, prior to introduction. In such case, there would have been no need for multiple introductions from previously isolated populations, but only several introductions from genetically diverse populations.     

Estimation of relatedness among non-pedigreed Yakutian cryo-bank bulls using molecular data: implications for conservation and breed management

Background

Yakutian cattle, the last remaining native cattle breed in Siberia, are well adapted to the extreme sub-arctic conditions. Nowadays only ca. 1200 purebred animals are left in Yakutia. The semen of six Yakutian bulls was stored in a cryo-bank without any pedigree documentation because of the traditional free herding style of the population.

Methods

To clarify the genetic relatedness between these bulls and to provide recommendations to use their semen in future conservation and breed management programs, we have analysed 30 autosomal microsatellites and mitochondrial DNA sequences in 60 individuals including the six for which semen has been stored. Four relatedness estimators were calculated. In addition, we assessed the value of the cryo-bank bulls for the preservation of genetic variation of the contemporary Yakutian cattle by calculating allelic and gene diversity estimates and mean molecular coancestries.

Results

On the basis of microsatellite variability, including the Yakutian cryo-bank bulls increases the allelic variation in the contemporary population by 3% and in the male subpopulation by 13%. In terms of the mean molecular coancestries, they are less related to the contemporary cow population than the breeding bulls and therefore could be used to reduce inbreeding in the living population. Although 30 loci are insufficient to resolve definitely their relatedness categories, the data suggest four pairs of cryo-bank bulls as possible half-sibs.

Conclusions

Our results show that even relatively limited cryo-bank storage of semen can carry allelic variation through a bottleneck. We propose a breeding scheme based on the rotation of breeding females and the division of cryo-bank bulls into three groups. Thus, if molecular data (e.g. autosomal microsatellite genotypes) for the contemporary population are available and based on relatively small-scale laboratory analyses, it is possible to avoid serious mistakes in their use for breeding applications. The approach suggested here based on the use of Yakutian cryo-bank semen can be easily extended to cryo-bank materials of other animals in future breeding programs.     

Ancient and current gene flow between two distantly related Mediterranean oak species, Quercus suber and Q. ilex

Background and Aims

Quercus suber and Q. ilex are distantly related and their distributions partially overlap. They hybridize occasionally, but the complete replacement of Q. suber chloroplast DNA (cpDNA) by that of Q. ilex was identified in two specific geographical areas. The objective of this study was to determine whether the contrasting situation reflected current or recent geographical interspecies gene flow variation or was the result of ancient introgression.

Methods

cpDNA PCR-RFLPs (restriction fragment length polymorphisms) and variation at ten nuclear microsatellite loci were analysed in populations of each species, in 16 morphologically intermediate individuals and the progeny of several of them. Interspecies nuclear introgression was based on individual admixture rates using a Bayesian approach with no a priori species assignment, and on a maximum-likelihood (ML) method, using allele frequencies in the allopatric populations of each species as controls. Gene flow was compared specifically between populations located within and outside the specific areas.

Key Results

High interspecies nuclear genetic differentiation was observed, with twice the number of alleles in Q. ilex than in Q. suber. According to Bayesian assignment, approx. 1 % of individuals had a high probability of being F₁ hybrids, and bidirectional nuclear introgression affected approx. 4 % of individuals in each species. Hybrid and introgressed individuals were identified predominantly in mixed stands and may have a recent origin. Higher proportions including allospecific genes recovered from past hybridization were obtained using the ML method. Similar rates of hybridization and of nuclear introgression, partially independent of cpDNA interspecies transfer suggestive of gene filtering, were obtained in the populations located within and outside the areas of complete cpDNA replacement.

Conclusions

The results did not provide evidence for geographical variation in interspecies gene flow. In contrast, historical introgression is supported by palynological records and constitutes the more reliable origin of cpDNA replacement in specific regions.Key words: cpDNA PCR-RFLPs, nuclear microsatellite (nSSR) variation, hybridization, interspecies genetic introgression, Quercus suber, Quercus ilex     

Sex, subdivision, and domestic dispersal of Trypanosoma cruzi lineage I in southern Ecuador

Background

Molecular epidemiology at the community level has an important guiding role in zoonotic disease control programmes where genetic markers are suitably variable to unravel the dynamics of local transmission. We evaluated the molecular diversity of Trypanosoma cruzi, the etiological agent of Chagas disease, in southern Ecuador (Loja Province). This kinetoplastid parasite has traditionally been a paradigm for clonal population structure in pathogenic organisms. However, the presence of naturally occurring hybrids, mitochondrial introgression, and evidence of genetic exchange in the laboratory question this dogma.

Methodology/Principal Findings

Eighty-one parasite isolates from domiciliary, peridomiciliary, and sylvatic triatomines and mammals were genotyped across 10 variable microsatellite loci. Two discrete parasite populations were defined: one predominantly composed of isolates from domestic and peridomestic foci, and another predominantly composed of isolates from sylvatic foci. Spatial genetic variation was absent from the former, suggesting rapid parasite dispersal across our study area. Furthermore, linkage equilibrium between loci, Hardy-Weinberg allele frequencies at individual loci, and a lack of repeated genotypes are indicative of frequent genetic exchange among individuals in the domestic/peridomestic population.

Conclusions/Significance

These data represent novel population-level evidence of an extant capacity for sex among natural cycles of T. cruzi transmission. As such they have dramatic implications for our understanding of the fundamental genetics of this parasite. Our data also elucidate local disease transmission, whereby passive anthropogenic domestic mammal and triatomine dispersal across our study area is likely to account for the rapid domestic/peridomestic spread of the parasite. Finally we discuss how this, and the observed subdivision between sympatric sylvatic and domestic/peridomestic foci, can inform efforts at Chagas disease control in Ecuador.     

Genetic variants of diabetes risk and incident cardiovascular events in chronic coronary artery disease

Objective

To determine whether information from genetic risk variants for diabetes is associated with cardiovascular events incidence.

Methods

From the about 30 known genes associated with diabetes, we genotyped single-nucleotide polymorphisms at the 10 loci most associated with type-2 diabetes in 425 subjects from the MASS-II Study, a randomized study in patients with multi-vessel coronary artery disease. The combined genetic information was evaluated by number of risk alleles for diabetes. Performance of genetic models relative to major cardiovascular events incidence was analyzed through Kaplan-Meier curve comparison and Cox Hazard Models and the discriminatory ability of models was assessed for cardiovascular events by calculating the area under the ROC curve.

Results

Genetic information was able to predict 5-year incidence of major cardiovascular events and overall-mortality in non-diabetic individuals, even after adjustment for potential confounders including fasting glycemia. Non-diabetic individuals with high genetic risk had a similar incidence of events then diabetic individuals (cumulative hazard of 33.0 versus 35.1% of diabetic subjects). The addition of combined genetic information to clinical predictors significantly improved the AUC for cardiovascular events incidence (AUC = 0.641 versus 0.610).

Conclusions

Combined information of genetic variants for diabetes risk is associated to major cardiovascular events incidence, including overall mortality, in non-diabetic individuals with coronary artery disease.

Clinical Trial Registration Information

Medicine, Angioplasty, or Surgery Study (MASS II). Unique identifier: ISRCTN66068876 URL.     

On the Origin of Indonesian Cattle

Background

Two bovine species contribute to the Indonesian livestock, zebu (Bos indicus) and banteng (Bos javanicus), respectively. Although male hybrid offspring of these species is not fertile, Indonesian cattle breeds are supposed to be of mixed species origin. However, this has not been documented and is so far only supported by preliminary molecular analysis.

Methods and Findings

Analysis of mitochondrial, Y-chromosomal and microsatellite DNA showed a banteng introgression of 10–16% in Indonesian zebu breeds. East-Javanese Madura and Galekan cattle have higher levels of autosomal banteng introgression (20–30%) and combine a zebu paternal lineage with a predominant (Madura) or even complete (Galekan) maternal banteng origin. Two Madura bulls carried taurine Y-chromosomal haplotypes, presumably of French Limousin origin. In contrast, we did not find evidence for zebu introgression in five populations of the Bali cattle, a domestic form of the banteng.

Conclusions

Because of their unique species composition Indonesian cattle represent a valuable genetic resource, which potentially may also be exploited in other tropical regions.