Genetic differentiation and gene flow between wild and cultivated Prunus avium: An analysis of molecular genetic evidence at a regional scale

The presence of conspecific wild-type and cultivar populations has been a common landscape feature for centuries. As orchards generally continue to expand towards the natural forest, two important issues are raised: the potential reduction of cultivar genetic diversity compared to wild populations and the extent of gene flow between the two population types. These questions were addressed in a study of Prunus avium in northern Greece using nine simple sequence repeat loci to analyse genetic variation in 93 wild-type individuals and 21 cultivars representing the local cultivated germplasm. Results showed a significant reduction of genetic diversity parameters in the cultivated germplasm compared to natural populations. Bayesian, frequency-based and Markov chain – Monte Carlo analyses have revealed that the wild and cultivar groups are genetically divergent and that realized between-group gene flow is almost completely absent. This result was further verified by a principal component analysis showing a clear separation of the two groups in low multivariate space after a principal coordinate analysis. The significant disjunction in flowering time and a considerable geographic distance between the two groups could primarily account for the absence of substantial gene flow. These findings indicate that local wild cherry can provide a source of genetic variation for future breeding in the genetically restricted cultivar group.     

Restricted gene flow between resident Oncorhynchus mykiss and an admixed population of anadromous steelhead

The species Oncorhynchus mykiss is characterized by a complex life history that presents a significant challenge for population monitoring and conservation management. Many factors contribute to genetic variation in O. mykiss populations, including sympatry among migratory phenotypes, habitat heterogeneity, hatchery introgression, and immigration (stray) rates. The relative influences of these and other factors are contingent on characteristics of the local environment. The Rock Creek subbasin in the middle Columbia River has no history of hatchery supplementation and no dams or artificial barriers. Limited intervention and minimal management have led to a dearth of information regarding the genetic distinctiveness of the extant O. mykiss population in Rock Creek and its tributaries. We used 192 SNP markers and collections sampled over a 5‐year period to evaluate the temporal and spatial genetic structures of O. mykiss between upper and lower watersheds of the Rock Creek subbasin. We investigated potential limits to gene flow within the lower watershed where the stream is fragmented by seasonally dry stretches of streambed, and between upper and lower watershed regions. We found minor genetic differentiation within the lower watershed occupied by anadromous steelhead (F_ST = 0.004), and evidence that immigrant influences were prevalent and ubiquitous. Populations in the upper watershed above partial natural barriers were highly distinct (F_ST = 0.093) and minimally impacted by apparent introgression. Genetic structure between watersheds paralleled differences in local demographics (e.g., variation in size), migratory restrictions, and habitat discontinuity. The evidence of restricted gene flow between putative remnant resident populations in the upper watershed and the admixed anadromous population in the lower watershed has implications for local steelhead productivity and regional conservation.     

Gene flow in admixed populations and implications for the conservation of the Western honeybee, Apis mellifera

Anthropogenic activity, especially modern apiculture, has considerable impact on the natural distribution of the Western honeybee, Apis mellifera, leading to the spread, replacement and fragmentation of many subspecies. This creates demand for the conservation of some subspecies, in particular, Apis mellifera mellifera, which once was widely distributed in Western Europe and nowadays is endangered through habitat loss and fragmentation. Moreover, A. m. mellifera may be further endangered by hybridisation in populations that now occur in artificial sympatry with other subspecies. Here, we quantify and compare individual hybridisation between sympatric and allopatric honeybee populations of A. m. mellifera and A. m. carnica using microsatellite markers and a Bayesian model-based approach. We had a special focus on pure breeding populations, which are a major tool in honeybee conservation. Our results demonstrate that subspecies are still highly differentiated, but gene flow is not prevented by the current management strategies, creating urgent demand for an improved conservation management of A. m. mellifera. However, the occurrence of a high number of pure individuals might suggest that some sort of hybrid barrier acts against the complete admixture of the two subspecies.     

Cherokee admixture and its estimation by the gene identity method: A critique

A previous application to a Cherokee hybrid population of two methods for measuring admixture proportions by use of gene identity probabilities was examined. The results suggest that neither procedure has been correctly applied, and that underlying statistical difficulties limit their reliability.     

On the relation between gene flow theory and genetic gain

In conventional gene flow theory the rate of genetic gain is calculated as the summed products of genetic selection differential and asymptotic proportion of genes deriving from sex-age groups. Recent studies have shown that asymptotic proportions of genes predicted from conventional gene flow theory may deviate considerably from true proportions. However, the rate of genetic gain predicted from conventional gene flow theory was accurate. The current note shows that the connection between asymptotic proportions of genes and rate of genetic gain that is embodied in conventional gene flow theory is invalid, even though genetic gain may be predicted correctly from it.     

Population analysis of vitamin D receptor polymorphisms and the role of genetic ancestry in an admixed population

The vitamin D receptor (VDR) is an essential protein related to bone metabolism. Some VDR alleles are differentially distributed among ethnic populations and display variable patterns of linkage disequilibrium (LD). In this study, 200 unrelated Brazilians were genotyped using 21 VDR single nucleotide polymorphisms (SNPs) and 28 ancestry informative markers. The patterns of LD and haplotype distribution were compared among Brazilian and the HapMap populations of African (YRI), European (CEU) and Asian (JPT+CHB) origins. Conditional regression and haplotype-specific analysis were performed using estimates of individual genetic ancestry in Brazilians as a quantitative trait. Similar patterns of LD were observed in the 5' and 3' gene regions. However, the frequency distribution of haplotype blocks varied among populations. Conditional regression analysis identified haplotypes associated with European and Amerindian ancestry, but not with the proportion of African ancestry. Individual ancestry estimates were associated with VDR haplotypes. These findings reinforce the need to correct for population stratification when performing genetic association studies in admixed populations.     

Selecting an appropriate genetic evaluation model for selection in a developing dairy sector

This study aimed to identify genetic evaluation models (GEM) to accurately select cattle for milk production when only limited data are available. It is based on a data set from the Pakistani Sahiwal progeny testing programme which includes records from five government herds, each consisting of 100 to 350 animals, with lactation records dating back to 1968. Different types of GEM were compared, namely: (1) multivariate v. repeatability model when using the first three lactations, (2) an animal v. a sire model, (3) different fixed effects models to account for effects such as herd, year and season; and (4) fitting a model with genetic parameters fixed v. estimating the genetic parameters as part of the model fitting process. Two methods were used for the comparison of models. The first method used simulated data based on the Pakistani progeny testing system and compared estimated breeding values with true breeding values. The second method used cross-validation to determine the best model in subsets of actual Australian herd-recorded data. Subsets were chosen to reflect the Pakistani data in terms of herd size and number of herds. Based on the simulation and the cross-validation method, the multivariate animal model using fixed genetic parameters was generally the superior GEM, but problems arise in determining suitable values for fixing the parameters. Using mean square error of prediction, the best fixed effects structure could not be conclusively determined. The simulation method indicated the simplest fixed effects structure to be superior whereas in contrast, the cross-validation method on actual data concluded that the most complex one was the best. In conclusion it is difficult to propose a universally best GEM that can be used in any data set of this size. However, some general recommendations are that it is more appropriate to estimate the genetic parameters when evaluating for selection purposes, the animal model was superior to the sire model and that in the Pakistani situation the repeatability model is more suitable than a multivariate.     

Gene flow and genetic microdifferentiation of a transplanted Tlaxcaltecan Indian population: Saltillo

This study is one portion of a research program on: (1) the genetic effects of transplantation of Mexican Indian populations into diverse environments; (2) the effect of Spanish colonization on the genetics of Mexico. In 1591, 91 Tlaxcaltecan families were transplanted into northern Mexico to aid in the pacification of the Chichimec Indians. The descendants of these migrants are the residents of the City of Saltillo in the state of Coahuila. Since Saltillo had undergone considerable growth due to industrialization, two adjoining barrios–Chamizal and La Minita–were chosen for this study. Chamizal is an old established barrio where the descendants of the Tlaxcaltecan colonialists resided, while La Minita is a small shantytown made up of squatters from the surrounding areas. Blood specimens were collected by venipuncture from 302 participants and analyzed for 16 blood groups, 13 serum protein and red blood cell enzyme systems. The gene frequency distribution suggests that the Saltillo populations are a triracial hybrid made up of Indian, African and Spanish parental con-tributions. Admixture estimates and genetic distance measures are both in agreement with the ethnohistory of the transplanted Tlaxcaltecan community.     

A deterministic model of admixture and genetic introgression: the case of Neanderthal and Cro-Magnon

There is an ongoing debate in the field of human evolution about the possible contribution of Neanderthals to the modern human gene pool. To study how the Neanderthal private alleles may have spread over the genes of Homo sapiens, we propose a deterministic model based on recursive equations and ordinary differential equations. If the Neanderthal population was large compared to the Homo sapiens population at the beginning of the contact period, we show that genetic introgression should have been fast and complete meaning that most of the Neanderthal private alleles should be found in the modern human gene pool in case of ancient admixture. In order to test/reject ancient admixture from genome-wide data, we incorporate the model of genetic introgression into a statistical hypothesis-testing framework. We show that the power to reject ancient admixture increases as the ratio, at the time of putative admixture, of the population size of Homo sapiens over that of Neanderthal decreases. We find that the power to reject ancient admixture might be particularly low if the population size of Homo sapiens was comparable to the Neanderthal population size.     

A method for the estimation of gene flow parameters from a population structure caused by restricted gene flow and genetic drift

Summary A method has been developed which enables the estimation of the plant gene flow parameters _p (pollen dispersal), _s (seed dispersal) and t (outcrossing rate) from a selection-free continuously structured population in equilibrium. The method uses Wright's F-coefficients and introduces a new F-function which describes the genetic similarity as a function of the spatial distance. The method has been elaborated for wind pollinated plant species but can be modified for insect pollination and for animal species. In practice allozymes will provide for the necessary neutral genetic variation. The more loci used and the more intermediate the gene frequencies, the more reliable the results. For the estimation of _p and t together (when the outcrossing rate is not known) at least two chromosomally unlinked loci are required. The method for estimating _s depends on whether the plant species is annual or perennial. The mechanism of selfing has been analysed by the explanation of the value of t by three components: population density (d), pollen flow (_p) and relative fertilization potential of own pollen (Z). The concepts of neighbourhood size and isolation by distance, developed by Wright, who used a single gene flow parameter , have been extended to the situation which is realistic for seed plants, using all three parameters _p, _s and t. When _p is large with respect to _s, _s largely determines the value of the neighbourhood size, whereas _p is the most dominating factor in isolation by distance. The use of local effective population size and mean gene transport per generation instead of neighbourhood size and neighbourhood area, respectively, is proposed to avoid confusion. Computer simulations have been carried out to check the validity and the reliability of the method. Populations of 200 plants, using two or three loci with intermediate allele frequencies, gave good results in the calculation of _p with known value of t and of _s and N_e. With unknown t, especially with lower values of t, larger populations of at least 1,000 plants are necessary to obtain reasonably accurate results for _p and mean gene transport per generation M.Grassland Species Research Group Publication No. 81     

The influence of cultural factors on the demography and pattern of gene flow from the Makiritare to the Yanomama Indians

A single village of Yanomama Indians was found to have frequencies of Di^a of 0.06 and of Ap^a of 0.08, in contrast to 40 other villages whereDi^a was absent and Ap^a quite rare. The source of these genes was identified as a village of Makiritare Indians, but the two allele frequencies were approximately the same or even higher in the Yanomama than in the Makiritare village. Demographic, social and cultural parameters affecting marriage and reproduction in the two tribes explain this. Genealogical relationships and informants' accounts collected in the field, when viewed against the traditional marriage practices, reproductive advantages of headmen, and differential treatment of captured women, indicate that the mating and reproduction parameters inherent in tribal social organization of this kind constitute an essential part of the explanation of the genetic findings. It is argued that mating systems of this sort are such that the probability of a new gene introduced by a captive surviving in the recipient population is a function of the sex of the initial carrier. The implications for tribalization and potentially radical changes in allele frequencies are briefly explored by considering aspects of settlement pattern and population fissioning known to characterize the tribes in question. Finally, it is shown that genetic sampling from a single location can and does result in unrepresentative allele frequencies when this single sample is taken to characterize the tribe as a whole.     

Using parentage analysis to examine gene flow and spatial genetic structure

Numerous approaches have been developed to examine recent and historical gene flow between populations, but few studies have used empirical data sets to compare different approaches. Some methods are expected to perform better under particular scenarios, such as high or low gene flow, but this, too, has rarely been tested. In this issue of Molecular Ecology , Saenz-Agudelo et   al . (2009 ) apply assignment tests and parentage analysis to microsatellite data from five geographically proximal (2–6 km) and one much more distant (1500 km) panda clownfish populations, showing that parentage analysis performed better in situations of high gene flow, while their assignment tests did better with low gene flow. This unusually complete data set is comprised of multiple exhaustively sampled populations, including nearly all adults and large numbers of juveniles, enabling the authors to ask questions that in many systems would be impossible to answer. Their results emphasize the importance of selecting the right analysis to use, based on the underlying model and how well its assumptions are met by the populations to be analysed.     

Population genetic structure and gene flow patterns between populations of the Antarctic icefish Chionodraco rastrospinosus

Aim A lack of genetic structure is predicted for Antarctic fish due to the duration of pelagic larval stages and the strength of the currents in the Southern Ocean, particularly the Antarctic Circumpolar Current. In this study we explored the population structure of the ocellated icefish, Chionodraco rastrospinosus, by means of analysing a total of 394 individuals collected at four geographical areas off the Antarctic Peninsula in the period 1996–2006. Location Elephant Island, southern South Shetlands, Joinville Island and South Orkneys in the Southern Ocean. Methods The spatio‐temporal genetic structure of Chionodraco rastrospinosus was explored using seven microsatellite loci. Existence and direction of gene flow across sampling locations were investigated using the isolation‐by‐migration procedure. Results Microsatellite data showed a lack of genetic structuring in the area studied, with no differences found at both the geographical or temporal level, and an eastward unidirectional gene flow among sites. This suggested a lack of genetic barriers for this species, attributable to larval dispersal following the Antarctic Circumpolar Current, which fits well with the predicted pattern for Antarctic fish. Re‐examination of genetic data of the closely related icefish Chaenocephalus aceratus, with similar larval duration but displaying genetically structured populations, indicated a weak but significant bidirectional gene flow. Main conclusions Our results point to a relationship that is more complex than expected between potential for dispersal and realized gene flow in the marine environment. In addition to ocean circulation and larval dispersal, other major life‐history traits might be driving connectivity, particularly larval retention.     

Matrilineal genetic structure and female-mediated gene flow in red grouse (Lagopus lagopus scoticus): an analysis using mitochondrial DNA

Abstract.— DNA sequence variation at the hypervariable 5'end of the mitochondrial control region was examined in 247 individuals to detect genetic divergence among 14 populations of red grouse ( Lagopus lagopus scoticus ) in northeastern Scotland. Ten haplotypes were resolved, several of which were shared among populations. Analysis of molecular variance, Nei's yST, and a cladistic estimate of the amount of gene flow indicated a lack of overall population differentiation. Patterns of overall panmixia are in stark contrast to previous reports of localized subdivision among the same set of populations detected using hypervariable microsatellite markers. Because grouse cocks are territorial and show extreme natal philopatry and females are the dispersing sex, such discordance could be explained by sexbiased dispersal, with extensive female-mediated gene flow preventing mitochondrial DNA divergence. However, it is difficult to reconcile how effective dispersal of females would not homogenize both mitochondrial and nuclear structure simultaneously. We use a model that examines the spatial and temporal dynamics of diparentally and uniparentally inherited genes to show that, under realistic ecological scenarios and with specific differences in the dispersal of males and females, the local effective size of the nuclear genome can be less than that of the mitochondrial and the patterns of structuring we observe are meaningful.     

Pollination between maize and teosinte: an important determinant of gene flow in Mexico

Gene flow between maize [Zea mays (L.)] and its wild relatives does occur, but at very low frequencies. Experiments were undertaken in Tapachula, Nayarit, Mexico to investigate gene flow between a hybrid maize, landraces of maize and teosinte (Z. mays ssp. mexicana, races Chalco and Central Plateau). Hybridization, flowering synchrony, pollen size and longevity, silk elongation rates, silk and trichome lengths and tassel diameter and morphology were measured. Hybrid and open-pollinated maize ears produced a mean of 8 and 11 seeds per ear, respectively, when hand-pollinated with teosinte pollen, which is approximately 1–2% of the ovules normally produced on a hybrid maize ear. Teosinte ears produced a mean of 0.2–0.3 seeds per ear when pollinated with maize pollen, which is more than one-fold fewer seeds than produced on a maize ear pollinated with teosinte pollen. The pollination rate on a per plant basis was similar in the context of a maize plant with 400–500 seeds and a teosinte plant with 30–40 inflorescences and 9–12 fruitcases per inflorescence. A number of other factors also influenced gene-flow direction: (1) between 90% and 95% of the fruitcases produced on teosinte that was fertilized by maize pollen were sterile; (2) teosinte collections were made in an area where incompatibility systems that limit fertilization are present; (3) silk longevity was much shorter for teosinte than for maize (approx. 4 days vs. approx. 11 days); (4) teosinte produced more pollen on a per plant basis than the landraces and commercial hybrid maize; (5) teosinte frequently produced lateral branches with silks close to a terminal tassel producing pollen. Collectively these factors tend to favor crossing in the direction of teosinte to maize. Our results support the hypothesis that gene flow and the subsequent introgression of maize genes into teosinte populations most probably results from crosses where teosinte first pollinates maize. The resultant hybrids then backcross with teosinte to introgress the maize genes into the teosinte genome. This approach would slow introgression and may help explain why teosinte continues to co-exist as a separate entity even though it normally grows in the vicinity of much larger populations of maize.     

Polyhydroxyalkaloids in the Aroid Tribes Nephthytideae and Aglaonemateae: Phytochemical Support for an Intertribal relationship

Using living and herbarium material, a survey of polyhydroxyalkaloids in species of 52 genera of Araceae revealed the presence of 2,5-dihydroxymethyl-3,4-dihydroxypyrrolidine (DMDP) and -homonojirimycin (HNJ) in leaves of Nephthytis Schott, Anchomanes Schott, Pseudohydrosme Engl. (tribe Nephthytideae), Aglaonema Schott and Aglaodorum Engl. (tribe Aglaonemateae). Levels were high in living plants, ranging from 0.1 to 1% dry weight DMDP and 0.04 to 0.6% HNJ. Isomers of HNJ, such as -3,4-di-epi-homonojirimycin, were also present in the five genera. Seven of the eight Nephthytis species examined also contained deoxymannojirimycin at levels, in fresh material, of 0.1–0.2% dry weight. Lower levels of DMDP (mean 0.009%) and HNJ (mean 0.002%) were detected in species of Amorphophallus Blume & Decne and Pseudodracontium N.E. Br. (tribe Thomsonieae). The similarity in polyhydroxyalkaloid chemistry between Nephthytideae and Aglaonemateae concurs with recently published chloroplast restriction site data in suggesting a relationship between these tribes.     

Association of serum lipid components and obesity with genetic ancestry in an admixed population of elderly women

The prevalence of metabolic disorders varies among ethnic populations and these disorders represent a critical health care issue for elderly women. This study investigated the correlation between genetic ancestry and body composition, metabolic traits and clinical status in a sample of elderly women. Clinical, nutritional and anthropometric data were collected from 176 volunteers. Genetic ancestry was estimated using 23 ancestry-informative markers. Pearsons correlation test was used to examine the relationship between continuous variables and an independent samples t-test was used to compare the means of continuous traits within categorical variables. Overall ancestry was a combination of European (57.49%), Native American (25.78%) and African (16.73%). Significant correlations were found for European ancestry with body mass index (r = 0.165; p = 0.037) and obesity (mean difference (MD) = 5.3%; p = 0.042). African ancestry showed a significant correlation with LDL (r = 0.159, p = 0.035), VLDL (r = −0.185; p = 0.014), hypertriglyceridemia (MD = 6.4%; p = 0.003) and hyperlipidemia (MD = 4.8%; p = 0.026). Amerindian ancestry showed a significant correlation with triglyceride levels (r = 0.150; p = 0.047) and hypertriglyceridemia (MD = 4.5%; p = 0.039). These findings suggest that genetic admixture may influence the etiology of lipid metabolism-related diseases and obesity in elderly women.     

Trade‐offs between growth and reproduction: an analysis of the quantitative genetic evidence

The assumption of a trade‐off between development time and fecundity, resulting from a positive correlation between body size and fecundity and between body size and development time, is a common feature of life history models. The present paper examines the evidence for such a trade‐off as indicated by genetic correlations between traits. The genetic covariances between traits are derived using a model in which maturation occurs when the organism achieves a genetically variable size threshold, and fecundity is an allometric function of body size with one genetically variable parameter (excluding body size itself). This model predicts that the heritabilities of the life history traits (growth rate, development time, fecundity) will not necessarily be less than the heritability of adult size (i.e. morphological traits). It is shown that if growth rate is genetically correlated with adult size then it is not possible, in general, to predict the sign of the genetic correlation between development time and fecundity. For particular cases the signs of the covariances between traits can be predicted. These predictions are tested using data drawn from the literature.     

AFLP analysis of the phenetic organization and genetic diversity of Vigna unguiculataL. Walp. reveals extensive gene flow between wild and domesticated types

Amplified fragment length polymorphisms (AFLPs) were used to evaluate genetic relationships within cowpea [Vigna unguiculata (L.) Walp.] and to assess the organization of its genetic diversity. Nei’s genetic distances were estimated for a total of 117 accessions including 47 domesticated cowpea (ssp. unguiculata var. unguiculata), 52 wild and weedy annuals (ssp. unguiculata var. spontanea), as well as 18 perennial accessions of the wild subspecies pubescens, tenuis and alba. AFLP variation was also used to study genetic variation among and within domesticated and wild accessions based on their geographical origin (western, eastern and southern Africa). Wild annual cowpea (var. spontanea) (H _T =0.175) was more diverse than domesticated cowpea (H _T =0.108). Wild cowpea was more diverse in eastern (H _S =0.168) than in western Africa (H _S =0.129), suggesting an eastern African origin for the wild taxon. The AFLP data were consistent with earlier findings of a unique domestication event in cowpea in the northern part of the continent and suggested that domestication in eastern or southern Africa was unlikely. It did not allow a more precise localization of domestication due to extensive gene flow between wild and domesticated forms that has led to a large crop-weed complex distributed over the entire African continent. In addition, wild materials from northeastern Africa are still lacking. Overall, the superiority of the AFLP technique over isozymes resided in its ability to uncover variation both within domesticated and wild cowpea, and should be a powerful tool once additional wild material becomes available. Received: 11 September 2000 / Accepted: 14 June 2001