Analysis of the non‐recombining Y chromosome defines polymorphisms in domestic pig breeds: ancestral bases identified by comparative sequencing

Sequences from 20 amplicons representing nine different loci and 11369bp from the short arm of the pig Y chromosome were compared using pools of DNA from different European and Chinese breeds. A total of 33 polymorphic sites were identified, including five indels and 28 single nucleotide polymorphisms (SNPs). Three high frequency SNPs within the coding regions of SRY were further analysed across 889 males representing 25 European and 25 Asian breeds or Lines, plus a European Line of Meishan. Two haplotypes seen to be associated with ‘European’ or ‘Chinese’ origin in the initial SNP discovery phase were found to be the most common in their respective groups of breeds in a more detailed genotyping study. Two further SRY haplotypes are relatively rare. One was found exclusively within Tamworth, at low frequency in Retinto, and in three Chinese breeds (Huai, Sahwutou and Xiaomeishan). The other uncommon haplotype is found exclusively in Bamajiang, two further Chinese breeds (Hangjiang Black and Longling) and two European rare breeds (Mangalica and Linderödssvin), but appears based on comparison with other suids to represent an ancestral sequence.     

TYR基因外显子1的序列变异

为了分析家猪与野猪的遗传多样性及起源,测定了来自１２个中国地方家猪品种、３个欧洲引进猪品种以及8个中国野猪和2个越南野猪共36个个体的酪氨酸酶基因（TYR）外显子1的序列,共检出６个单核苷酸多态性位点（SNPs）,且这6个位点的变异均为同义突变,根据这些变异可将酪氨酸酶基因DNA序列归结为４种单倍型。结合已发表的数据,构建了简约中介网络图。 在网络图中,单倍型TYR*2主要为欧洲家猪与欧洲野猪和三条亚洲家猪染色体。大部分亚洲家猪和野猪共享单倍型TYR*1,表明这是一个亚洲类型的单倍型;同时也有部分欧洲家猪与野猪携带这一单倍型。 而单倍型TYR*3和TYR*4为本研究检测到的稀有单倍型,这两种单倍型主要由中国家猪与亚洲野猪组成。这种网络图结构支持家猪的欧洲和亚洲独立起源学说,同时也表明相当部分的欧洲家猪品种受到亚洲猪的基因渗透,而少量中国家猪和日本野猪也受到了欧洲猪的基因渗透。     

Mitochondrial diversity in European and Chinese pigs is consistent with population expansions that occurred prior to domestication

Mitochondrial DNA (mtDNA) diversity in European and Asian pigs was assessed using 1536 samples representing 45 European and 21 Chinese breeds. Diagnostic nucleotide differences in the cytochrome b (Cytb) gene between the European and Asian mtDNA variants were determined by pyrosequencing as a rapid screening method. Subsequently, 637bp of the hypervariable control region was sequenced to further characterize mtDNA diversity. All sequences belonged to the D1 and D2 clusters of pig mtDNA originating from ancestral wild boar populations in Europe and Asia, respectively. The average frequency of Asian mtDNA haplotypes was 29% across European breeds, but varied from 0 to 100% within individual breeds. A neighbour-joining (NJ) tree of control region sequences showed that European and Asian haplotypes form distinct clusters consistent with the independent domestication of pigs in Asia and Europe. The Asian haplotypes found in the European pigs were identical or closely related to those found in domestic pigs from Southeast China. The star-like pattern detected by network analysis for both the European and Asian haplotypes was consistent with a previous demographic expansion. Mismatch analysis supported this notion and suggested that the expansion was initiated before domestication.     

Molecular anatomy of the cytoplasmic domain of bovine growth hormone receptor, a quantitative trait locus

Quantitative trait loci (QTL) studies have indicated growth hormone receptor (GHR) as a candidate gene affecting cattle milk yield and composition. In order to characterize genetic variation at GHR in cattle, we studied European and East African breeds with different histories of selection, and Bos grunniens, Ovis aries, Sus scrofa, Bison bison and Rangifer tarandus as references. We sequenced most of the cytoplasmic domain (900 bp of exon 10), 89 bp of exon 8, including the putative causative mutation for the QTL effect, and 390 bp of intron 8 for comparison. In the cytoplasmic domain, seven synonymous and seven non-synonymous single nucleotide polymorphisms (SNP) were identified in cattle. Three non-synonymous SNPs were found in sheep and one synonymous SNP in yak, while other studied species were monomorphic. Three major haplotypes were observed, one unique to African breeds, one unique to European breeds and one shared. Bison and yak haplotypes are derivatives of the European haplotype lineage. Most of the exon 10 non-synonymous cattle SNPs appear at phylogenetically highly conserved sites. The polymorphisms in exon 10 cluster around a ruminant-specific tyrosine residue, suggesting that this site may act as an additional signalling domain of GHR in ruminants. Alternative explanations for the persistent polymorphism include balancing selection, hitch-hiking, pleiotropic or sexually antagonistic fitness effects or relaxed functional constraints.     

Mitochondrial haplotypes reveal a strong genetic structure for three Indian sheep breeds

This survey represents the first characterization of mitochondrial DNA diversity within three breeds of Indian sheep (two strains of the Deccani breed, as well as the Bannur and Garole breeds) from different geographic regions and with divergent phenotypic characteristics. A 1061-bp fragment of the mitochondrial genome spanning the control region, a portion of the 12S rRNA gene and the complete phenyl tRNA gene, was sequenced from 73 animals and compared with the corresponding published sequence from European and Asian breeds and the European Mouflon (Ovis musimon). Analysis of all 156 sequences revealed 73 haplotypes, 52 of which belonged to the Indian breeds. The three Indian breeds had no haplotypes in common, but one Indian haplotype was shared with European and other Asian breeds. The highest nucleotide and haplotype diversity was observed in the Bannur breed (0.00355 and 0.981 respectively), while the minimum was in the Sangamneri strain of the Deccani breed (0.00167 and 0.882 respectively). All 52 Indian haplotypes belonged to mitochondrial lineage A. Therefore, these Indian sheep are distinct from other Asian and European breeds studied so far. The relationships among the haplotypes showed strong breed structure and almost no introgression among these Indian breeds, consistent with Indian sheep husbandry, which discourages genetic exchange between breeds. These results have implications for the conservation of India's ovine biodiversity and suggest a common origin for the breeds investigated.     

PKLR- GBA region shows almost complete linkage disequilibrium over 70 kb in a set of worldwide populations

Haplotype diversity in a genomic region of approximately 70 kb in 1q21 between genes PKLR and GBA was characterized by typing one single nucleotide polymorphism (SNP) in PKLR, two SNPs in GBA and one short tandem repeat polymorphism (STRP) in PKLR in 1792 chromosomes from 17 worldwide populations. Two other SNPs in GBA were typed in three African populations. Most chromosomes carried one of either two phylogenetically distinct haplotypes with different alleles at each site. Allele diversity at the STRP was tightly linked to haplotype background. Linkage disequilibrium (LD) was highly significant for all SNP pairs in all populations, although it was, on average, slightly higher in non-African populations than in sub-Saharan Africans. Variation at PKLR-GBA was also tightly linked to that at the GBA pseudogene, 16 kb downstream from GBA. Thus, a 90 kb-long LD block was observed, which points to a low recombination rate in this region. Detailed haplotype phylogeny suggests that the chimpanzee GBA haplotype is not one of the two most frequent haplotypes. Based on variability at the PKLR STRP and on the geographical distribution of LD, the expansion of the two main haplotypes may have predated the "Out of Africa" expansion of anatomically modern humans. LD and STRP variability in non-Africans are approximately 87% of those in Africans, in contrast with other loci; this implies that the "out of Africa" bottleneck may have had a broad distribution of effects across loci.     

Extensive sharing of chloroplast haplotypes among East Asian Cerris oaks: The imprints of shared ancestral polymorphism and introgression

Shared ancestral polymorphism and introgression are two main causes of chloroplast DNA (cpDNA) haplotype sharing among closely related angiosperms. In this study, we explored the roles of these two processes in shaping the phylogeographic patterns of East Asian Cerris oaks by examining the geographic distributions of randomly and locally distributed shared haplotypes, which coincide with the expectations of shared ancestry and introgression, respectively. We sequenced 1340 bp of non‐coding cpDNA from Quercus acutissima (n = 418) and Q. chenii (n = 183) and compiled previously published sequence data of Q. variabilis (n = 439). The phylogenetic relationships among haplotypes were examined using a median‐joining network. The geographic patterns of interspecifically shared haplotypes were assessed to test whether nearby populations have a higher degree of interspecific cpDNA sharing than distant ones. We identified a total of 27 haplotypes that were grouped into three non‐species‐specific lineages with overlapping distributions. Ancestral haplotypes were extensively shared and randomly distributed across populations of the three species. Some young haplotypes were locally shared in mountainous areas that may have been shared refugia. The local exchange of cpDNA resulted in an excess of similar haplotypes between nearby populations. Our study demonstrated that the haplotype sharing pattern among East Asian Cerris oaks reflected the imprints of both shared ancestral polymorphism and introgression. This pattern was also associated with the relatively stable climates and complex landscapes in East Asia, which not only allowed the long‐term persistence of ancestral lineages but also connected the survived populations across refugia.     

Haplotype diversity in the endosperm specific β-amylase gene Bmy1 of cultivated barley (Hordeum vulgare L.)

Five single nucleotide polymorphism (SNP) sites corresponding to substitutions in the protein sequence of the β-amylase gene Bmy1 at amino acid (AA) positions 115, 165, 233, 347 and 430 were genotyped in 493 cultivated barley accessions by Pyrosequencing and a CAPS assay. A total of 6 different haplotypes for the Bmy1 gene were discovered of which 4 haplotypes were identified as previously described alleles Bmy1-Sd1, Bmy1-Sd2L, Bmy1-Sd2H and Bmy1-Sd3, while 2 haplotypes were new. A broad spectrum of haplotypes was found in spring barleys, while the winter barleys were dominated by the newly described haplotype Bmy1-Sd4. Individual haplotype frequencies varied between the geographic regions.Three pairs of SNP loci within the gene showed highly significant (P<0.0001) elevated values of linkage disequilibrium (LD) with r ² > 0.6. In the European and Asian subpopulations different loci were in linkage disequilibrium due to the differences in haplotype frequency distributions. By applying LD data to select haplotype tagging SNPs, three SNP sites corresponding to AA positions 115, 233 and 347 were identified that allowed to discriminate 4 haplotypes and to capture 91.6% of the available diversity by distinguishing 452 out of 493 accessions. In a subset of 2-rowed German spring barley varieties 4 SNPs and 2 haplotypes had a significant association with the malting quality parameter final attenuation limit which is related to the total amylolytic enzymatic activity.     

Further resolution of porcine phylogeny in Southeast Asia by Thai mtDNA haplotypes

The mtDNA control region of 72 Thai native pigs and 11 Thai wild boars indigenous to Northern Thailand was comparatively sequenced. In total, 36 nucleotide variations that accounted for 24 haplotypes have been described (TNH01 to TNH20 and TWH01 to TWH04). These haplotypes and further publicly available mtDNA haplotypes were used to assess phylogenetic relationships. Twenty-three of the 24 haplotypes became integrated into the Asian clade of the phylogenetic tree and eight of them recapitulated another major cluster of haplotypes within this clade (Thai haplogroup, THG). Only haplotype TNH01 fit in with the European clade of the phylogenetic tree. An additional analysis using 510 bp of the mtDNA incorporated the THG haplotypes in to clade MTSEA (mountainous and Southeast Asian distribution) to form haplogroup MTSEA-THG. Recently, MTSEA was renamed in MC3. MC3 contains only signatures of pigs scattered across the Indo-Burma Biodiversity Hotspot (IBBH), a region including Thailand to the Kra Isthmus. Here we propose a putative independent porcine domestication event in South-east Asia (SEA). All haplotypes of haplogroup MTSEA-THG have revealed unique and previously unknown nucleotide signatures at positions 24 (nucleotide A) and at positions 183 (nucleotide C) that differentiate them from all other porcine mtDNA haplotypes.     

Worldwide genetic analysis of the CFTR region

Mutations at the cystic fibrosis transmembrane conductance regulator gene (CFTR) cause cystic fibrosis, the most prevalent severe genetic disorder in individuals of European descent. We have analyzed normal allele and haplotype variation at four short tandem repeat polymorphisms (STRPs) and two single-nucleotide polymorphisms (SNPs) in CFTR in 18 worldwide population samples, comprising a total of 1,944 chromosomes. The rooted phylogeny of the SNP haplotypes was established by typing ape samples. STRP variation within SNP haplotype backgrounds was highest in most ancestral haplotypes-although, when STRP allele sizes were taken into account, differences among haplotypes became smaller. Haplotype background determines STRP diversity to a greater extent than populations do, which indicates that haplotype backgrounds are older than populations. Heterogeneity among STRPs can be understood as the outcome of differences in mutation rate and pattern. STRP sites had higher heterozygosities in Africans, although, when whole haplotypes were considered, no significant differences remained. Linkage disequilibrium (LD) shows a complex pattern not easily related to physical distance. The analysis of the fraction of possible different haplotypes not found may circumvent some of the methodological difficulties of LD measure. LD analysis showed a positive correlation with locus polymorphism, which could partly explain the unusual pattern of similar LD between Africans and non-Africans. The low values found in non-Africans may imply that the size of the modern human population that emerged "Out of Africa" may be larger than what previous LD studies suggested.     

Phylogeographic study of brown trout from Serbia,based on mitochondrial DNA control region analysis

In order to illuminate the phylogeography of brown trout (Salmo trutta) populations in the Balkan state of Serbia, the 561 bp 5''-end of mtDNA control region of 101 individuals originating from upland tributaries of the Danubian, Aegean and Adriatic drainages were sequenced and compared to corresponding brown trout sequences obtained in previous studies. Among 15 haplotypes found, 14 were considered native, representing the Danubian and Adriatic lineages of the brown trout, while one haplotype (ATcs1), found only in two individuals originating from two stocked rivers, corresponded to the Atlantic lineage and was considered introduced. Native haplotypes exhibited a strong geographic pattern of distribution: the Danubian haplotypes were strictly confined to the Danubian drainage, while the Adriatic haplotypes dominated in the Aegean and Adriatic drainages; most of the total molecular variance (69%) was attributed to differences among the drainages. Phylogenetic reconstruction, supplemented with seven haplotypes newly described in this study, suggested a sister position of the Atlantic-Danubian and Adriatic-Mediterranean-marmoratus ("southern") phylogenetic group, and pointed to the existence of a distinct clade, detected within the "southern" group. The data obtained confirmed our expectation of the existence of high genetic diversity in Balkan trout populations, and we recommend more widespread surveys covering trout stocks from the region.     

Nucleotide Diversity and Molecular Evolution of the <Emphasis Type="Italic">ALK</Emphasis> Gene in Cultivated Rice and its Wild Relatives

Gelatinization temperature (GT), an important parameter for rice cooking quality, is mainly regulated by the ALK gene encoding starch synthase IIa. Here, we reported the nucleotide diversity of the ALK gene in 122 cultivated accessions and 199 wild rice accessions that were collected around the Pearl River Basin in China. A total of 93 single nucleotide ploymorphisms (SNPs) were identified, with an average of one SNP per 40 bp. Tajima D statistics revealed that the DNA sequences covering the last exon have probably evolved under balancing selection. Based on two functional SNPs (an A to G substitution at 4198 bp and a GC to TT dinucleotide substitution at 4330/4331 bp), three haplotypes, G/GC, G/TT, and A/GC, were identified in both wild and cultivated accessions, with the G/GC haplotype being predominant. Interestingly, the A/GC haplotype was exclusively found in the wild accessions from Guangdong province, while the G/TT haplotype was only present in the wild accessions from Jiangxi province and Hainan Island. This suggests that the G/TT and A/GC variants may have arisen independently and undergone balancing selection on separate haplotypes in multiple populations. Our result supports earlier hypothesis that cultivated rice was independently domesticated from multiple domestication events in China. Our study aids in the understanding of the domestication process that led to the improvement of rice grain quality.     

The phylogeography of the asian sockeye salmon <Emphasis Type="Italic">Oncorhynchus nerka</Emphasis>, inferred from the data on the variability of mitochondrial SNP loci: Analysis of scenarios for post-glacial expansion of the species over the asian coast of the pacific ocean

The variability of three single nucleotide polymorphism (SNP) loci in the mitochondrial DNA (mtSNP) is analyzed for sockeye salmon populations over a major part of the species range, from Chukotka to the Kuril Islands. Two basic mtSNP haplotypes, GCC and GTT, have been revealed in 20 sockeye samples from 15 lake–river systems on the Asian coast of the Pacific Ocean. In most of the samples, the ratio of the haplotypes is approximately equal. The GTT haplotype dominates the populations from the Kuril Islands (except Shumshu Island); only the GCC haplotype has been found in the sample from the Commander Islands. This geographic pattern of haplotype distribution was presumably caused by the historico-demographic events related to the formation of the Asian sockeye range in the Middle–Late Pleistocene, viz., fragmentation of the species range and subsequent secondary contact between previously diverged populations. These data provide a basis to consider different scenarios for the formation of the modern diversity of sockeye mtSNP haplotypes, all of which suggest multiple expansions of the species to Asian waters during the periods of oceanic transgression after the Pleistocene glaciations.     

Strong signature of natural selection within an FHIT intron implicated in prostate cancer risk

Previously, a candidate gene linkage approach on brother pairs affected with prostate cancer identified a locus of prostate cancer susceptibility at D3S1234 within the fragile histidine triad gene (FHIT), a tumor suppressor that induces apoptosis. Subsequent association tests on 16 SNPs spanning approximately 381 kb surrounding D3S1234 in Americans of European descent revealed significant evidence of association for a single SNP within intron 5 of FHIT. In the current study, re-sequencing and genotyping within a 28.5 kb region surrounding this SNP further delineated the association with prostate cancer risk to a 15 kb region. Multiple SNPs in sequences under evolutionary constraint within intron 5 of FHIT defined several related haplotypes with an increased risk of prostate cancer in European-Americans. Strong associations were detected for a risk haplotype defined by SNPs 138543, 142413, and 152494 in all cases (Pearson's chi(2) = 12.34, df 1, P = 0.00045) and for the homozygous risk haplotype defined by SNPs 144716, 142413, and 148444 in cases that shared 2 alleles identical by descent with their affected brothers (Pearson's chi(2) = 11.50, df 1, P = 0.00070). In addition to highly conserved sequences encompassing SNPs 148444 and 152413, population studies revealed strong signatures of natural selection for a 1 kb window covering the SNP 144716 in two human populations, the European American (pi = 0.0072, Tajima's D = 3.31, 14 SNPs) and the Japanese (pi = 0.0049, Fay & Wu's H = 8.05, 14 SNPs), as well as in chimpanzees (Fay & Wu's H = 8.62, 12 SNPs). These results strongly support the involvement of the FHIT intronic region in an increased risk of prostate cancer.     

Sequence polymorphism at the human apolipoprotein AII gene (<Emphasis Type="Italic">APOA2</Emphasis>): unexpected deficit of variation in an African-American sample

A 3.3-kb region, encompassing the APOA2 gene and 2 kb of 5' and 3' flanking DNA, was re-sequenced in a "core" sample of 24 individuals, sampled without regard to the health from each of three populations: African-Americans from Jackson (Miss., USA), Europeans from North Karelia (Finland), and non-Hispanic European-Americans from Rochester, (Minn., USA). Fifteen variable sites were identified (14 SNPs and one multi-allelic microsatellite, all silent), and these sites segregated as 18 sequence haplotypes (or nine, if SNPs only are considered). The haplotype distribution in the core African-American sample was unusual, with a deficit of particular haplotypes compared with those found in the other two samples, and a significantly (P<0.05) low level of nucleotide diversity relative to patterns of polymorphism and divergence at other human loci. Six of the 14 SNPs, whose variation captured the haplotype structure of the core data, were then genotyped by oligonucleotide ligation assay in an additional 2183 individuals from the same three populations (n=843, n=452, and n=888, respectively). All six sites varied in each of the larger "epidemiological" samples, and together, they defined 19 SNP haplotypes, seven with relative frequencies greater than 1% in the total sample; all of these common haplotypes had been identified earlier in the core re-sequencing survey. Here also, the African-American sample showed significantly lower SNP heterozygosity and haplotype diversity than the other two samples. The deficit of polymorphism is consistent with a population-specific non-neutral increase in the relative frequency of several haplotypes in Jackson.     

Association analysis of the RGS4 gene in Han Chinese and Scottish populations with schizophrenia

We investigated the RGS4 as a susceptibility gene for schizophrenia in Chinese Han (184 trios and 138 sibling pairs, a total of 322 families) and Scottish (580 cases and 620 controls) populations using both a family trio and case-control design. Both the samples had statistical power greater than 70% to detect a heterozygote genotype relative risk of >1.2 for frequent RGS4-risk alleles. We genotyped four single nucleotide polymorphisms (SNPs) which have previously been associated with schizophrenia as either individually or part of haplotypes. Allele frequencies and linkage disequilibrium between the SNPs was similar in the two populations. In the Chinese sample, no individual SNPs or any of their haplotypes were associated with schizophrenia. In the Scottish population, one SNP (SNP7) was significantly over-represented in the cases compared with the controls (0.44 vs. 0.38; A allele; chi(2) 7.08, P = 0.011 after correction for correlation between markers by permutation testing). One two-marker haplotype, composed of alleles T and A of SNP4 and SNP7, respectively, showed individual significance after correction by permutation testing (chi(2) 6.8; P = 0.04). None of the full four-marker haplotypes showed association, including the G-G-G-G haplotype previously associated with schizophrenia in more than one sample and the A-T-A-A haplotype. Thus, our data do not directly replicate previous associations of RGS4, but association with SNP 7 in the Scottish population provides some support for a role in schizophrenia susceptibility. We cannot conclusively exclude RGS4, as associated haplotypes are likely to be surrogates for unknown causative alleles, whose relationship with overlying haplotypes may differ between the population groups. Differences in the association seen across the two populations could result from methodological factors such as diagnostic differences but most likely result from ethnic differences in haplotype structures within RGS4.     

Low genetic diversity and significant structuring of the common hamster populations <Emphasis Type="Italic">Cricetus cricetus</Emphasis> in Poland revealed by the mtDNA control region sequence variation

The genetic diversity of 12 populations in the present range of the common hamster Cricetus cricetus (Linnaeus, 1758) in Poland was established. The 366 bp of the mtDNA control region was sequenced for 195 individuals. As few as seven haplotypes were found and their distribution was geographically structured. The large geographic areas were fixed or almost fixed for a single haplotype and three groups of populations, that do not share any haplotypes, have been defined. Proportions of genetic diversity attributable to variation between groups of populations, between populations within groups and within populations were 93.64, 1.92 and 4.45% (SAMOVA: p < 0.001 for all estimates), respectively. Such pattern of variation is most probably the result of historical, postglacial bottlenecks and present genetic drift after the population decline in the last few decades.     

Genomewide distribution of high-frequency, completely mismatching SNP haplotype pairs observed to be common across human populations

Knowledge of human haplotype structure has important implications for strategies of disease-gene mapping and for understanding human evolutionary history. Many attributes of SNPs and haplotypes appear to exhibit highly nonrandom behavior, suggesting past operation of selection or other nonneutral forces. We report the exceptional abundance of a particular haplotype pattern in which two high-frequency haplotypes have different alleles at every SNP site (hence the name "yin yang haplotypes"). Analysis of common haplotypes in 62 random genomic loci and 85 gene coding regions in humans shows that the proportion of the genome spanned by yin yang haplotypes is 75%-85%. Population data of 28 genomic loci in Drosophila melanogaster reveal a similar pattern. The high recurrence (>/=85%) of these haplotype patterns in four distinct human populations suggests that the yin yang haplotypes are likely to predate the African diaspora. The pattern initially appeared to suggest deep population splitting or maintenance of ancient lineages by selection; however, coalescent simulation reveals that the yin yang phenomenon can be explained by strictly neutral evolution in a well-mixed population.