Genetic Structure and Differentiation of Three Chinese Indigenous Cattle Populations

Levels of genetic differentiation, gene flow, and genetic structure of three indigenous cattle populations (Luxi, Bohai, and Minnan) and two reference cattle populations (Chinese Holstein and Qinhai yak) in China were estimated using the information from 12 microsatellites, and 141 microsatellite alleles were identified. The mean number of alleles per locus ranged from 2.9005 in yak to 4.9722 in Holstein. The observed heterozygosity ranged from 0.5325 (yak) to 0.7719 (Holstein); 29 private alleles were detected. The global heterozygote deficit across all populations amounted to 58.5% (p < 0.001). The overall significant (p < 0.001) deficit of heterozygotes because of inbreeding within breeds amounted to 43.2%. The five cattle populations were highly differentiated (F _st = 26.9%, p < 0.001) at all loci. The heterozygote deficit within the population was highest in Luxi cattle and lowest in yak. The average number of effective migrants exchanged per generation was highest (1.149) between Luxi and Holstein, and lowest (0.509) between Luxi and yak. With the application of prior population information, cluster analysis achieved posterior probabilities from 91% to 98% of correctly assigning individuals to populations. Combining the information of cluster analysis, gene flow, and Structure analysis, the five cattle populations belong to three genetic clusters, a taurine (Luxi and Chinese Holstein), a zebu (Bohai and Minnan), and a yak cluster. This indicates that Bohai black is closer to Bos indicus than Luxi cattle. The evolution and development of three indigenous cattle populations are discussed.     

Comparative analysis of ISSR marker polymorphism in populations of yak (Bos mutus) and in F1 hybrids between yak and cattle in the Sayan-Altai region

The genetic variability in seven yak populations from the Sayan-Altai region and in F₁ hybrids between yak and cattle (khainags) was investigated with the help of a technique that involves the use of inter simple sequence repeat (ISSR) markers generated with PCR primers (AG)⁹C and (GA)⁹C. Samples for the analysis were collected in Mongolia, Tuva, and Altai from 2008 through 2012. The examined yak populations differed in in the presence/absence of ISSR fragments, as well as in their frequency. In total, 46 ISSR fragments were identified using two marker systems; the proportion of polymorphic loci constituted 76% and 90% for the AG-ISSR and GA-ISSR markers, respectively. For the total sample of yaks, total genetic diversity (H _t), within-population diversity (H _s), and interpopulation diversity (G _st) constituted 0.081, 0.044, and 0.459 for the AG-ISSR and 0.137, 0.057, and 0.582 for the GA-ISSR markers, respectively. Based on ISSR finger-printing, species- and breed-specific DNA patterns were described for the three groups of animals (yaks, cattle, khainags). For the domestic yak, the species-specific profile was represented by eight ISSR fragments. Genetic relationships between the yak populations, cattle breeds, and khainags were examined with the help of four different approaches used in the analysis of population structure: estimation of phylogenetic similarity, multidimensional scaling, principal component analysis, and cluster analysis. Clear evidence on the differentiation of the populations examined at the interspecific, as well as at intraspecific, level were obtained. Similar (relative), as well as remote (isolated), yak populations were identified. Khainags occupy an intermediate position between yak and cattle. However, the data on the ISSR-PCR marker polymorphism (genome polymorphism, population structure) indicate that part of the analyzed khainag genome was more similar to the yak genome than to the cattle genome.     

Genetic admixture despite ecological segregation in a North African sparrow hybrid zone (Aves,Passeriformes, Passer domesticus × Passer hispaniolensis)

Under different environmental conditions, hybridization between the same species might result in different patterns of genetic admixture. Particularly, species pairs with large distribution ranges and long evolutionary history may have experienced several independent hybridization events over time in different zones of overlap. In birds, the diverse hybrid populations of the house sparrow (Passer domesticus) and the Spanish sparrow (Passer hispaniolensis) provide a striking example. Throughout their range of sympatry, these two species do not regularly interbreed; however, a stabilized hybrid form (Passer italiae) exists on the Italian Peninsula and on several Mediterranean islands. The spatial distribution pattern on the Eurasian continent strongly contrasts the situation in North Africa, where house sparrows and Spanish sparrows occur in close vicinity of phenotypically intermediate populations across a broad mosaic hybrid zone. In this study, we investigate patterns of divergence and admixture among the two parental species, stabilized and nonstabilized hybrid populations in Italy and Algeria based on a mitochondrial marker, a sex chromosomal marker, and 12 microsatellite loci. In Algeria, despite strong spatial and temporal separation of urban early‐breeding house sparrows and hybrids and rural late‐breeding Spanish sparrows, we found strong genetic admixture of mitochondrial and nuclear markers across all study populations and phenotypes. That pattern of admixture in the North African hybrid zone is strikingly different from i) the Iberian area of sympatry where we observed only weak asymmetrical introgression of Spanish sparrow nuclear alleles into local house sparrow populations and ii) the very homogenous Italian sparrow population where the mitogenome of one parent (P. domesticus) and the Z‐chromosomal marker of the other parent (P. hispaniolensis) are fixed. The North African sparrow hybrids provide a further example of enhanced hybridization along with recent urbanization and anthropogenic land‐use changes in a mosaic landscape.     

Influence of niche similarity on hybridization between Myriophyllum sibiricum and M. spicatum

The impact of ecological factors on natural hybridization is of widespread interest. Here, we asked whether climate niche influences hybridization between the two closely related plant species Myriophyllum sibiricum and M. spicatum. Eight microsatellite loci and two chloroplast fragments were used to investigate the occurrence of hybridization between these two species in two co‐occurring regions: north‐east China (NEC) and the Qinghai‐Tibetan Plateau (QTP). The climate niches of the species were quantified by principal component analysis with bioclimatic data, and niche comparisons were performed between the two species in each region. Reciprocal hybridization was observed, and M. sibiricum was favoured as the maternal species. Furthermore, hybrids were rare in NEC but common in the QTP. Accordingly, in NEC, the two species were climatically distinct, and hybrids only occurred in the narrow geographical or ecological transition zone, whereas in the QTP, obvious niche overlaps were found for the two species, and hybrids occurred in multiple contact zones. This association between hybridization pattern and climate niche similarity suggests that the level of hybridization was promoted by niche overlap. Compared with the parental species, similar climate niches were found for the hybrid populations in the QTP, indicating that other environmental factors rather than climate were important for hybrid persistence. Our findings highlight the significance of climate niche with respect to hybridization patterns in plants.     

Population genetic structure of Bombus terrestris in Europe: Isolation and genetic differentiation of Irish and British populations

The genetic structure of the earth bumblebee (Bombus terrestris L.) was examined across 22 wild populations and two commercially reared populations using eight microsatellite loci and two mitochondrial genes. Our study included wild bumblebee samples from six populations in Ireland, one from the Isle of Man, four from Britain and 11 from mainland Europe. A further sample was acquired from New Zealand. Observed levels of genetic variability and heterozygosity were low in Ireland and the Isle of Man, but relatively high in continental Europe and among commercial populations. Estimates of F_st revealed significant genetic differentiation among populations. Bayesian cluster analysis indicated that Irish populations were highly differentiated from British and continental populations, the latter two showing higher levels of admixture. The data suggest that the Irish Sea and prevailing south westerly winds act as a considerable geographical barrier to gene flow between populations in Ireland and Britain; however, some immigration from the Isle of Man to Ireland was detected. The results are discussed in the context of the recent commercialization of bumblebees for the European horticultural industry.     

Detecting adaptive trait loci in nonmodel systems: divergence or admixture mapping?

Mapping adaptive trait loci (ATL) underlying ecological divergence is an essential step towards understanding the processes that generate phenotypic diversity. Technological advances have made it possible to sequence exomes in nonmodel systems, providing an efficient means of analysing functional genetic variants. Divergence scans of genetic markers for outlier loci, or ‘divergence mapping’, have been used to map locally adapted genes, but this approach is likely to be underpowered when background divergence is elevated. Genotype–phenotype association tests in admixed populations, or ‘admixture mapping’, may provide a useful approach for mapping locally adapted loci when neutral divergence is high. To determine the power and limits of divergence mapping, we simulated exomes containing a single ATL across two parental populations of varying neutral divergence, estimated divergence and quantified the power to identify the ATL. We found that divergence mapping had very high power when background F_ST is <0.2, but decreased dramatically above this level. To evaluate the utility of admixture mapping, we simulated exomes from admixed populations, then simulated phenotypes, conducted genotype–phenotype association tests and found that even two generations of random mating after admixture could provide high mapping power in scenarios where pure divergence mapping was ineffective (F_ST = 0.35). Moreover, admixture mapping had high power across all levels of divergence after 20 generations since admixture. Together with high‐throughput exome sequencing, admixture mapping could be used to map ATL in systems such as Heliconius butterflies or Gryllus crickets when experimental design and analytical approach are chosen accordingly.     

Phylogeny of Y chromosomes from bovine species

The Bovini tribe contains domestic and wild cattle‐like species, several of which are cross‐fertile. We present a completely resolved Y‐chromosomal phylogeny, which is better in agreement with autosomal phylogeny, morphological data, cross‐fertility and estimates of divergence times than mitochondrial data. The tree links Bos grunniens (yak) to Bison, so the commonly accepted Bos genus is not monophyletic. Therefore, we advocate the term Poephagus as designation of yak. This work illustrates the resolving power of Y‐chromosomal variation for cladistic studies of closely related species. © The Willi Hennig Society 2008.     

Microgeographic structure of Anopheles gambiae in western Kenya based on mtDNA and microsatellite loci

The population genetic structure of the Anopheles gambiae in western Kenya was studied using length variation at five microsatellite loci and sequence variation in a 648-nt mtDNA fragment. Mosquitoes were collected from houses in villages spanning up to 50 km distance, The following questions were answered, (i) Are mosquitoes in a house more related genetically to each other than mosquitoes between houses? (ii) What degree of genetic differentiation occurs on these geographical scales? (iii) How consistent are the results obtained with both types of genetic markers? At the house level, no differentiation was detected by F_ST and R_ST, and the band sharing index test revealed no significant associations of alleles across loci. Likewise, indices of kinship based on mtDNA haplotypes in houses were even lower than in the pooled sample. Therefore, the hypothesis that mosquitoes in a house are more related genetically was rejected. At increasing geographical scales, microsatellite allele distributions were similar among all population samples and no subdivision of the gene pool was detected by F_ST or R_ST. Likewise, estimates of haplotype divergence of mtDNA between populations were not higher than the within population estimates, and mtDNA-based F_ST values were not significantly different from zero. That sequence variation in mtDNA provided matching results with microsatellite loci (while high genetic variation was observed in all loci), suggested that this pattern represents the whole genome. The minimum area associated with a deme of A. gambiae in western Kenya is therefore larger than 50 km in diameter.     

Novel Y‐chromosome short tandem repeats in Sus scrofa and their variation in European wild boar and domestic pig populations

Y‐chromosome markers are important tools for studying male‐specific gene flow within and between populations, hybridization patterns and kinship. However, their use in non‐human mammals is often hampered by the lack of Y‐specific polymorphic markers. We identified new male‐specific short tandem repeats (STRs) in Sus scrofa using the available genome sequence. We selected four polymorphic loci (5–10 alleles per locus), falling in one duplicated and two single‐copy regions. A total of 32 haplotypes were found by screening 211 individuals from eight wild boar populations across Europe and five domestic pig populations. European wild boar were characterized by significantly higher levels of haplotype diversity compared to European domestic pigs (H_D = 0.904 ± 0.011 and H_D = 0.491 ± 0.077 respectively). Relationships among STR haplotypes were investigated by combining them with single nucleotide polymorphisms at two linked genes (AMELY and UTY) in a network analysis. A differentiation between wild and domestic populations was observed (F_ST = 0.229), with commercial breeds sharing no Y haplotype with the sampled wild boar. Similarly, a certain degree of geographic differentiation was observed across Europe, with a number of local private haplotypes and high diversity in northern populations. The described Y‐chromosome markers can be useful to track male inheritance and gene flow in wild and domestic populations, promising to provide insights into evolutionary and population genetics in Sus scrofa.     

Isolation and characterization of microsatellite loci in Bemisia tabaci

Bemisia tabaci (Hemiptera: Aleyrodidae) is a haplo‐diploid species with a global distribution demonstrating strong geographical structure with eight recognizable genetic groups. Fifteen microsatellite loci (335 alleles, 6–44 alleles per locus) were derived from four of the eight groups and were then screened across 33 populations. These loci clearly differentiate the populations. The microsatellites amplified best in individuals from genetic groups rep‐resenting the Mediterranean, Middle East, Asia (three groups) and Australasia/Oceania and amplified less well with populations from sub‐Saharan Africa and the New World. This differential amplification pattern is a direct result of the relatedness to the microsatellite source material.     

Seed-Mediated Gene Flow Promotes Genetic Diversity of Weedy Rice within Populations: Implications for Weed Management

Increased infestation of weedy rice—a noxious agricultural pest has caused significant reduction of grain yield of cultivated rice (Oryza sativa) worldwide. Knowledge on genetic diversity and structure of weedy rice populations will facilitate the design of effective methods to control this weed by tracing its origins and dispersal patterns in a given region. To generate such knowledge, we studied genetic diversity and structure of 21 weedy rice populations from Sri Lanka based on 23 selected microsatellite (SSR) loci. Results indicated an exceptionally high level of within-population genetic diversity (H_e = 0.62) and limited among-population differentiation (F_st = 0.17) for this predominantly self-pollinating weed. UPGMA analysis showed a loose genetic affinity of the weedy rice populations in relation to their geographical locations, and no obvious genetic structure among populations across the country. This phenomenon was associated with the considerable amount of gene flow between populations. Limited admixture from STRUCTURE analyses suggested a very low level of hybridization (pollen-mediated gene flow) between populations. The abundant within-population genetic diversity coupled with limited population genetic structure and differentiation is likely caused by the considerable seed-mediated gene flow of weedy rice along with the long-distance exchange of farmer-saved rice seeds between weedy-rice contaminated regions in Sri Lanka. In addition to other effective weed management strategies, promoting the application of certified rice seeds with no weedy rice contamination should be the immediate action to significantly reduce the proliferation and infestation of this weed in rice ecosystems in countries with similar rice farming styles as in Sri Lanka.     

The role of recombination on genome‐wide patterns of local ancestry exemplified by supplemented brook charr populations

Assessing the immediate and long‐term evolutionary consequences of human‐mediated hybridization is of major concern for conservation biology. Several studies have documented how selection in interaction with recombination modulates introgression at a genome‐wide scale, but few have considered the dynamics of this process within and among chromosomes. Here, we used an exploited freshwater fish, the brook charr (Salvelinus fontinalis), for which decades of stocking practices have resulted in admixture between wild populations and an introduced domestic strain, to assess both the temporal dynamics and local chromosomal variation in domestic ancestry. We provide a detailed picture of the domestic ancestry patterns across the genome using about 33,000 mapped single nucleotide polymorphisms genotyped in 611 individuals from 24 supplemented populations. For each lake, we distinguished early‐ and late‐generation hybrids using information regarding admixture tracts. To assess the selective outcomes following admixture we then evaluated the relationship between recombination and admixture proportions at three different scales: the whole genome, chromosomes and within 2‐Mb windows. This allowed us to detect a wide range of evolutionary mechanisms varying along the genome, as reflected by the finding of favoured or disfavoured introgression of domestic haplotypes. Among these, the main factor modulating local ancestry was probably the presence of deleterious recessive mutations in the wild populations, which can be efficiently hidden to selection in the presence of long admixture tracts. Overall, our results emphasize the relevance of taking into consideration local ancestry information to assess both the temporal and the chromosomal variation in local admixture ancestry toward better understanding post‐hybridization evolutionary outcomes.     

Associations of single nucleotide polymorphisms in candidate genes with the polled trait in Datong domestic yaks

The domestic yak (Bos grunniens) is an iconic symbol of animal husbandry at high altitudes. Yaks exhibit unique external characteristics including long hair and large horns. However, hornless yaks can be found in different breeds and different populations. The hornless trait is also known as polled, and the POLL locus has been fine‐mapped to chromosome 1 in cattle (Bos taurus), although the underlying genetic basis of the polled trait is still unclear in the yak. Thus, we performed an association study to identify the genetic polymorphisms responsible for the polled trait in the yak. Fifty polled Datong domestic yaks and 51 horned individuals were selected randomly from a huge herd and were used as the case and control groups respectively for the association analysis. Twelve genes located in the candidate region of the POLL locus in cattle were used as references to detect DNA polymorphisms related to yak polledness, which were analyzed by sequencing and a high‐resolution melting test. We applied Fisher's exact test and haplotype analysis to show that a 147‐kb segment that included three protein‐coding genes C1H21orf62, GCFC1 and SYNJ1 was the most likely location of the POLL mutation in domestic yaks.     

Genotyping of microsatellite markers to study genetic structure of the wild striped snakehead Channa striata in Malaysia

Genetic variability and differences in wild striped snakehead Channa striata from Malaysia were analysed by genotyping nine novel nuclear microsatellite loci. Analysis revealed moderate‐to‐high genetic diversity in most of the populations, indicative of large effective population sizes. The highly diversified populations are admixed populations and, therefore, can be recommended as potential candidates for selective breeding and conservation since they each contain most of the alleles found in their particular region. Three homogenous groups of the wild populations were identified, apparently separated by effective barriers, in accordance with contemporary drainage patterns. The highest population pairwise F_ST found between members of the same group reflects the ancient population connectivity; yet prolonged geographical isolation resulted in adaptation of alleles to local contemporary environmental change. A significant relationship between genetic distance and geographical isolation was observed (r = 0·644, P < 0·01). Anthropogenic perturbations indicated apparent genetic proximity between distant populations.     

Hybridization and introgression across different ploidy levels in the Neotropical orchids Epidendrum fulgens and E. puniceoluteum (Orchidaceae)

The hypothesis of gene flow between species with large differences in chromosome numbers has rarely been tested in the wild, mainly because species of different ploidy are commonly assumed to be reproductively isolated from each other because of instantaneous and strong postzygotic barriers. In this study, a broad‐scale survey of molecular variation was carried out between two orchid species with different ploidy levels: Epidendrum fulgens (2n = 2x = 24 chromosomes) and Epidendrum puniceoluteum (2n = 4x = 52 chromosomes). To test the strength of their reproductive barriers, we investigated the distribution of genetic variation in sympatric and allopatric populations of these two species and conducted crossing experiments. Nuclear and plastid microsatellite loci were used to genotype 463 individuals from eight populations across the geographical range of both species along the Brazilian coastal plain. All six sympatric populations analysed presented hybrid zones, indicating that hybridization between E. fulgens and E. puniceoluteum is a common phenomenon. Bayesian assignment analysis detected the presence of F₁ and F₂ individuals and also signs of introgression, demonstrating a high potential for interspecific gene flow. Introgression occurs preferentially from E. fulgens to E. puniceoluteum. Pure parental individuals of both species display strong genotype–habitat associations, indicating that environment‐dependent selection could be acting in all hybrid zones. This study suggests that hybridization and introgression are evolutionary processes playing a role in the diversification of Epidendrum and indicates the importance of investigations of hybrid zones in understanding reproductive barriers and speciation processes in Neotropical orchid species.     

Environmental and anthropogenic correlates of hybridization between westslope cutthroat trout (Oncorhynchus clarkii lewisi) and introduced rainbow trout (O. mykiss)

Hybridization with introduced taxa is one of the major threats to the persistence of native biodiversity. The westslope cutthroat trout (Oncorhynchus clarkii lewisi) is found in southeastern British Columbia and southwestern Alberta, Canada, and adjacent areas of Montana, Idaho, and Washington State, USA. Through much of this area, native populations are threatened by hybridization with introduced rainbow trout (O. mykiss). We surveyed 159 samples comprising over 5,000 fish at 10 microsatellite DNA loci to assess the level of admixture between native westslope cutthroat trout (wsct) and introduced rainbow trout in southwestern Alberta. Admixture levels (q_wsct of 0 = pure rainbow trout, q_wsct of 1.0 = pure westslope cutthroat trout) ranged from <0.01 to 0.99 and averaged from 0.72 to 0.99 across seven drainage areas. Regression tree analyses indicated that water temperature, elevation, distance to the nearest stocking site, and distance to the nearest railway line were significant components of a model that explained 34 % of the variation across sites in q_wsct across 58 localities for which habitat variables were available. Partial dependence plots indicated that admixture with rainbow trout increased with increasing water temperature and distance to the nearest railway line, but decreased with increasing elevation and distance from stocking site to sample site. Our results support the hypothesis that westslope cutthroat trout may be less susceptible to hybridization with rainbow trout in colder, higher elevation streams, and illustrate the interaction between abiotic and anthropogenic factors in influencing hybridization between native and introduced taxa.     

Effect of the Aegean Sea barrier between Europe and Asia on differentiation in Juniperus drupacea (Cupressaceae)

Juniperus drupacea is an eastern Mediterranean mountain tree with a disjunct geographical range. We hypothesized that this disjunct occurrence (the Peloponnese in Europe and the Taurus and Lebanon Mountains in Asia) should be reflected in the patterns of genetic and morphological diversity and differentiation. Nuclear microsatellite markers (nSSR) and biometric variables of the cones and seeds were examined on material sampled from four populations in Europe and eight in Asia. The Asian populations were characterized by a higher level of genetic diversity than the European populations. The genetic differentiation among populations was moderate but significant (F_ST = 0.101, P < 0.001). According to the clustering performed with BAPS, six genetically and geographically groups of populations were found: I and II from the Peloponnese; III from the Taurus Mountains; IV and V from the Anti‐Taurus Mountains; and VI from the Lebanon Mountains. The level of genetic differentiation among these six groups (4.30%, P = 0.012) probably reflects long‐lasting genetic isolation during the Pleistocene, as limited genetic admixture was found. In accordance with genetic analysis, the biometric investigations indicated a high level of morphological divergence between the European and Asian populations of the species, with further differentiation between the populations from the Taurus and Lebanon Mountains.     

Genetic Diversity and Divergence in Chinese Yak (Bos grunniens) Populations Inferred from Blood Protein Electrophoresis

In 6 Chinese yak (Bos. Grunniens) populations including 177 yaks, 34 blood protein loci were studied by horizontal starch gel electrophoresis, four of these loci (AKP, ALB, LDH-1, TF) were found to be polymorphic. The percentage of polymorphic loci(P) is 0.118, the mean individual heterozygosity(H) is 0.015, which means a low level of genetic diversity in the whole Chinese yak population. The coefficient of gene differentiation (G _ST ) is 0.0625, which indicated an almost-indistinguishable divergence among different populations at the level of blood protein electrophoresis.     

The genetics of brown coat color and white spotting in domestic yaks (Bos grunniens)

Domestic yaks (Bos grunniens) exhibit two major coat color variations: a brown vs. wild‐type black pigmentation and a white spotting vs. wild‐type solid color pattern. The genetic basis for these variations in color and distribution remains largely unknown and may be complicated by a breeding history involving hybridization between yaks and cattle. Here, we investigated 92 domestic yaks from China using a candidate gene approach. Sequence variations in MC1R, PMEL and TYRP1 were surveyed in brown yaks; TYRP1 was unassociated with the coloration and excluded. Recessive mutations from MC1R, or p.Gln34*, p.Met73Leu and possibly p.Arg142Pro, are reported in bovids for the first time and accounted for approximately 40% of the brown yaks in this study. The remaining 60% of brown individuals correlated with a cattle‐derived deletion mutation from PMEL (p.Leu18del) in a dominant manner. Degrees of white spotting found in yaks vary from color sidedness and white face, to completely white. After examining the candidate gene KIT, we suggest that color‐sided and all‐white yaks are caused by the serial translations of KIT (Cs₆ or Cs₂₉) as reported for cattle. The white‐faced phenotype in yaks is associated with the KIT haplotype S^wf. All KIT mutations underlying the serial phenotypes of white spotting in yaks are identical to those in cattle, indicating that cattle are the likely source of white spotting in yaks. Our results reveal the complex genetic origins of domestic yak coat color as either native in yaks through evolution and domestication or as introduced from cattle through interspecific hybridization.