西藏牦牛SRAP遗传多样性及分类进化研究

用4对SRAP分子标记引物对西藏11个牦牛类群和四川麦洼牦牛的DNA进行扩增,研究其遗传多样性和分类关系。结果表明,在335头牦牛中,共得到29个基因位点,其中有19个多态位点,多态率占65.52%。12个牦牛群体间的Nei’s遗传多样性和Shannon多样性指数分别为0.048 2和0.073 4,遗传相似系数在0.781 1-0.989 1。巴青牦牛和康布牦牛的遗传多样性指数较其他类群高,分别为0.095 5和0.090 0;桑日牦牛类群的遗传多样性指数最低,仅为0.008 5。这些结果表明,12个牦牛类群的SRAP遗传多样性较低。根据Nei’s遗传距离,利用UPGMAM构建聚类关系图结果显示,嘉黎牦牛、帕里牦牛、类乌齐牦牛、桑桑牦牛、康布牦牛、巴青牦牛、丁青牦牛、斯布牦牛和麦洼牦牛聚为一大类,然后依次才与桑日牦牛、工布江达牦牛和江达牦牛相聚在一起,显示在SRAP分子遗传标记所反映的牦牛基因组的遗传结构中,江达牦牛、工布江达牦牛和桑日牦牛与其他牦牛群体间的亲缘关系较远,牦牛的这种亲缘关系与其地理分布也不一致,说明这12个牦牛的起源、演化关系较复杂,有待于进一步研究分析。     

Genetic diversity and structure in fragmented populations of the endangered species Ranunculus cabrerensis (Ranunculaceae): implications for conservation

Ranunculus cabrerensis is an endemic and endangered species of the Northwestern Iberian Peninsula. The molecular markers AFLP and ISSR were used to investigate the genetic diversity and population structure of four populations across its known distribution. Fifteen selective primer combinations of AFLP and seventeen ISSR primer combinations produced a total of 2830 and 103 unambiguously repeatable fragments respectively, of which 97.57 and 81.38% were polymorphic for both markers. The genetic diversity of R. cabrerensis at species level was high (H _E = 0.294 by ISSR and H _E = 0.191 by AFLP) and differentiation between sampled locations was also relatively high (G _ST = 0.316 and 0.158 by ISSR and AFLP analysis respectively) compared to other studies of endangered and rare species using the same techniques. The analysis of molecular variance (AMOVA) indicated that the main genetic variation was within sampled locations (73% by AFLP; 52% by ISSR), even though the variation among locations was also significant. Principal Coordinates, NeighborNet and Bayesian analyses revealed a weak but significant relationship between the genetic structures of different populations in R. cabrerensis, with gene flow acting as a homogenizing force that prevents stronger differentiation of populations. Finally, suggestions for conservation strategies to preserve the genetic resources of this species are outlined.     

Assessment of cattle genetic introgression into domestic yak populations using mitochondrial and microsatellite DNA markers

Hybridization between yak Poephagus grunniens and taurine Bos taurus or indicine B. indicus cattle has been widely practiced throughout the yak geographical range, and gene flow is expected to have occurred between these species. To assess the impact of cattle admixture on domestic yak, we examined 1076 domestic yak from 29 populations collected in China, Bhutan, Nepal, India, Pakistan, Kyrgyzstan, Mongolia and Russia using mitochondrial DNA and 17 autosomal microsatellite loci. A cattle diagnostic marker‐based analysis reveals cattle‐specific mtDNA and/or autosomal microsatellite allele introgression in 127 yak individuals from 22 populations. The mean level of cattle admixture across the populations, calculated using allelic information at 17 autosomal microsatellite loci, remains relatively low (mY_cattle = 2.66 ± 0.53% and Q_cattle = 0.69 ± 2.58%), although it varies a lot across populations as well as among individuals within population. Although the level of cattle admixture shows a clear geographical structure, with higher levels of admixture in the Qinghai‐Tibetan Plateau and Mongolian and Russian regions, and lower levels in the Himalayan and Pamir Plateau region, our results indicate that the level of cattle admixture is not significantly correlated with the altitude across geographical regions as well as within geographical region. Although yak‐cattle hybridization is primarily driven to produce F₁ hybrids, our results show that the subsequent gene flow between yak and cattle took place and has affected contemporary genetic make‐up of domestic yak. To protect yak genetic integrity, hybridization between yak and cattle should be tightly controlled.     

Characterization of the Russian beef cattle breed gene pools using inter simple sequence repeat DNA analysis (ISSR analysis)

The gene pools of beef cattle breeds bred in Russia were characterized on the basis of inter simple sequence repeat DNA analysis (ISSR analysis). Samples of Aberdeen Angus, Kalmyk, and Kazakh Whitehead breeds from Russia, as well as of Hereford breed, hybrids of Kazakh Whitehead and Hereford breeds, and Kazakh Whitehead breed from the Republic of Kazakhstan, were examined. In the examined breeds, 27 AG-ISSR fragments were identified, 25 of which were polymorphic. The examined breeds were different both in the fragment profiles (the presence/absence of individual ISSR fragments) and in their frequencies. It was demonstrated that the hybrid animals lacked some ISSR fragments that were present with high frequencies in parental forms, suggesting considerable genome rearrangement in the hybrid animals (at the regions of microsatellite localization) in crossings of the individuals from different breeds. The level of genetic diversity in Russian beef breeds was consistent with the values typical of farmed populations (breeds). The genetic diversity parameters assessed by applying Nei’s gene diversity index and the Shannon index varied from 0.0218 to 0.0605 and from 0.0225 to 0.0819, respectively. The highest Shannon index value was detected in the Kalmyk breed (0.0837) and Kazakh Whitehead breed from Russia (0.0819), and the highest level of Nei’s gene diversity index was found in the Kalmyk breed (0.0562) and in both populations of the Kazakh Whitehead breed (0.0509 and 0.0605). The high level of genetic similarity (according to Nei) was revealed between Russian beef cattle breeds and Hereford cattle: 0.839 (for the Kazakh Whitehead breed from Russia) and 0.769 (for the Kalmyk breed).     

Novel Y‐chromosome polymorphisms in Chinese domestic yak

Y‐chromosome‐specific haplotypes (Y‐haplotypes) constructed using single nucleotide polymorphisms (Y‐SNPs) in the MSY (male‐specific region of the Y‐chromosome) are valuable in population genetic studies. But sequence variants in the yak MSY region have been poorly characterized so far. In this study, we screened a total of 16 Y‐chromosome‐specific gene segments from the ZFY, SRY, UTY, USP9Y, AMELY and OFD1Y genes to identify Y‐SNPs in domestic yaks. Six novel Y‐SNPs distributed in the USP9Y (g.223C>T), UTY19 (g.158A>C and g.169C>T), AMELY2 (g.261C>T), OFD1Y9 (g.165A>G) and SRY4 (g.104G>A) loci, which can define three Y‐haplotypes (YH1, YH2 and YH3) in yaks, were discovered. YH1 was the dominant and presumably most ancient haplotype based on the comparison of UTY19 locus with other bovid species. Interestingly, we found informative UTY19 markers (g.158A>C and g.169C>T) that can effectively distinguish the three yak Y‐haplotypes. The nucleotide diversity was 1.7 × 10^?4 ± 0.3 × 10^?4, indicating rich Y‐chromosome diversity in yaks. We identified two highly divergent lineages (YH1 and YH2 vs. YH3) that share similar frequencies (YH1 +  YH2: 0.82–0.89, YH3: 0.11–0.18) among all three populations. In agreement with previous mtDNA studies, we supported the hypothesis that the two highly divergent lineages (YH1 and YH2 vs. YH3) derived from a single gene pool, which can be explained by the reunion of at least two paternal populations with the divergent lineages already accumulated before domestication. We estimated a divergence time of 408 110 years between the two divergent lineages, which is consistent with the data from mitochondrial DNA in yaks.     

Genetic variability of MHC class II DQB exon 2 alleles in yak (Bos grunniens)

The major histocompatibility class (MHC) DQ molecules are dimeric glycoproteins revealing antigen presentation to CD⁴⁺ T cells. In the present study, the exon 2 of the MHC class II DQB gene from 32 yaks (Bos grunniens) was cloned, sequenced and compared with previously reported patterns for other bovidae. It was revealed by sequence analyses that there are 25 DQB exon 2 alleles among 32 yaks, all alleles are found to belong to DQB1 loci. These alleles exhibited a high degree of nucleotide and amino acid polymorphisms with most amino acid variations occurring at positions forming the peptide-binding sites. The DQB loci were analyzed for patterns of synonymous (d _S) and non-synonymous (d _N) substitution. The yak was observed to be under strong positive selection in the DQB exon 2 peptide-binding sites (d _N = 0.15, P < 0.001). It appears that this variability among yaks confers the ability to mount immune responses to a wide variety of peptides or pathogens.     

Comparison of ISSR polymorphism among cattle breeds

Polymorphism analysis of DNA fragments flanked by (AG)₉C and (GA)₉C inverted dinucleotide microsatellite repeats in 766 animals of 19 cattle breeds and one breeding type revealed 66 fragments, of which 64 were polymorphic. The breeds proved to differ in the frequency and presence or absence of amplified DNA fragments at the genomic level, indicating that ISSR fingerprinting is informative for differentiating the PCR product spectra and cattle breeds. Multilocus ISSR polymorphism analysis identified the group of fragments that can be used as Bos taurus and B. indicus species markers to describe the standards of breeds, their genetic profiles, and breed-specific patterns. Based on ISSR polymorphism, a prototypal gene pool of cattle was constructed and the breeds closest to it were identified. Genetic diversity analysis made it possible to assume that an optimal mean heterozygosity is characteristic of cattle breeds and that deviations from this optimum are indicative of various processes occurring in the population (breed).     

Isolation and characterization of polymorphic microsatellites in the genome of Yak (Bos grunniens)

The researches on yak genetics and breeding were extremely restricted due to lacking of reliable DNA molecular markers. The microsatellites with repeat motif (AC)_n/(GT)_n in yak genome were enriched by Dynal magnetic beads and the gene libraries containing (AC)_n/(GT)_n were constructed. Among the 92 identified and sequenced positive clones, 40 contained perfect repeats (43.48 %), 41 contained imperfect repeats (44.57 %) and 11 contained compound repeats (11.96 %). As compared with the percentage of perfect repeats, no significant increases of imperfect repeats were observed in yak genome, which indicated that the level of adaptive evolution of the ability to repair damaged genomic DNA for yaks were high enough to endure the natural pressure of nucleotide substitution resulted from ultraviolet irradiation in high-altitude areas. Totally 19 polymorphic microsatellite loci were screened and genotyped on the basis of electropherograms on an ABI 310 Genetic Analyzer. All the loci exhibited moderate to high-level polymorphisms in a test population of Bos grunniens and the polymorphic information content ranged from 0.299 to 0.861 (mean 0.678). The newly isolated (AC)_n/(GT)_n repeats from yak genome will display their potential values in examining intra-population genetic structure and inter-population relationships, and also in investigating molecular markers for production and adaptive traits of individual/population.     

Population Genetic Study of Fagopyrum tataricum from Western Himalaya Using ISSR Markers

Inter-simple sequence repeat (ISSR) markers were used to analyze genetic diversity and relatedness of 15 germplasms of Fagopyrum tataricum. Samples representing 75 individuals were collected from a range of altitudes in the Western Himalaya. The 13 ISSR primers revealed 98.1% polymorphism among populations, whereas average polymorphism was extremely low (2.18%) within populations. The coefficient of population differentiation was 0.9750, with limited gene flow (N _m) of 0.0128. The average PIC value of the ISSR markers was high (0.812), with a marker ratio of 0.65 and marker index of 6.66. The genetic diversity of F. tataricum significantly correlated with altitude and gene diversity, Shannon’s index, and the percentage of polymorphic bands. The genetic diversity among populations showed broad genetic base and provided a developmental strategy for crop improvement.     

贺兰山牦牛冬春季的生境选择

在2009年12月－2010年1月和2010年4－5月,采用样线法结合直接观察法对贺兰山牦牛的冬春季生境选择进行了研究。结果表明,牦牛冬季对11种生境因子有选择性,偏爱山地针叶林带,海拔小于2 000 m,优势乔木为灰榆,坡度小于10,下坡位,距水源距离大于1 200 m,人为干扰距离2 000－4 000 m,隐蔽级大于70 %;春季牦牛对13种利用生境生态因子有选择性,偏爱于亚高山灌丛和草甸带,海拔大于3 000 m,乔木密度小于1株,乔木高度小于3 m,乔木距离大于3 m,灌木密度大于4 0株,灌木距离小于1 m,植被盖度大于7 0 %,上坡位,距水源的距离小于8 00 m,人为干扰距离大于4 000 m,隐蔽级大于7 0 %。冬春季牦牛在海拔、植被类型、地形特征、优势乔木、灌木种类、坡位、坡向、人为干扰距离、距水源距离上存在显著差异。主成分分析表明冬季第一主成分的贡献率21.100 %,其中绝对值较大的相关系数是乔木距离、优势乔木、乔木高度和乔木密度;春季第一主成分的贡献率是31.247 %,其中绝对值较大的相关系数是植被类型、海拔高度、地形特征和灌木密度。与其他分布地区的牦牛相比,贺兰山地区的牦牛能适应当地的地理特征和气候环境。     

Endangered wild yak (Bos grunniens) in the Tibetan plateau and adjacent regions: Population size,distribution, conservation perspectives and its relation to the domestic subspecies

The yak (Bos grunniens) is a long-haired bovid, endemic to the Tibetan Plateau and the adjacent high-altitude regions. The domesticated subspecies of yak (B. grunniens grunniens) are abundant and closely associated with the livelihoods of herders, while the wild subspecies of yak (B. grunniens mutus) are endangered due primarily to anthropogenic effects. The endangered status of wild yaks calls for consideration, if we are to secure its long term survival, hence this study. Here we hope to provide baseline information necessary for further research and protection of the wild yak resources. We use published data to discuss their evolution, their characteristics as well as their distribution in the Tibetan Plateau and the adjacent high-altitude regions. We were able to come up with a world wild yak distribution map, which may be useful for establishing protected areas, as well as updating the species IUCN Red List Status. From the data available, we were also able to provide an estimate of the wild yak population in China (∼22,000 wild yaks living in China), corresponding to 90% of the total world population. We further discuss the major threats to yaks, and we give some suggestions for future and sustainable conservation.     

Application of bovine microsatellite markers for genetic diversity analysis of Swiss yak (Poephagus grunniens)

In order to assess the applicability of bovine microsatellite markers for population genetic studies in Swiss yak, 131 bovine microsatellite markers were tested on a panel of 10 animals. Efficient amplification was observed for 124 markers (94.6%) with a total of 476 alleles, of which 117 markers (94.3%) were polymorphic. The number of alleles per locus among the polymorphic markers ranged from two to nine. Seven loci (ILSTS005, BMS424B, BMS1825, BMS672, BM1314, ETH123 and BM6017) failed to amplify yak genomic DNA. Two cattle Y-chromosome specific microsatellite markers (INRA126 and BM861) amplified genomic DNA from both male and female yaks. However, two additional markers on cattle Y-chromosome (INRA124 and INRA189) amplified DNA from only males. Of the polymorphic markers, 24 microsatellites proposed by CaDBase for within- and cross-species comparisons and two additional highly polymorphic markers (MHCII and TGLA73) were used to investigate the genetic variability and the population structure of a Swiss yak herd that included 51 additional animals. The polymorphic information content ranged from 0.355 to 0.752, while observed heterozygosity (HO) ranged from 0.348 to 0.823. Furthermore, a set of 13 markers, organized into three multiplex polymerase chain reactions, was evaluated for routine parentage testing. This set provided an exclusion probability in a family of four yaks (both parents and two offspring) of 0.995. These microsatellites serve as useful tools for genetic characterization of the yak, which continues to be an important domestic livestock species.     

The extent of clonality and genetic diversity in the rare Caldesia grandis (Alismataceae): Comparative results for RAPD and ISSR markers

Genetic variation and clonal diversity of three natural populations of the rare, highly clonal marsh herb Caldesia grandis Samuelsson were investigated using random amplified polymorphic DNA (RAPD) and inter-simple sequence repeat (ISSR) markers. Both of the markers worked effectively in clone identification of C. grandis. RAPD markers detected more diversity than ISSR markers in the three populations examined. Of the 60 RAPD primers screened, seven produced highly reproducible bands. Using these primers, a total of 61 DNA fragments were generated with 52 (85.25%) being polymorphic indicating considerable genetic variation at the species level. Analysis of molecular variance (AMOVA) showed that a large proportion of genetic variation (81.5%) resided within populations, while only a small proportion (18.5%) resided among populations. With the use of 52 polymorphic RAPD markers, we were able to identify 127 genets among 342 samples from three populations. The proportion of distinguishable genets (PD: mean 0.37), Simpson's diversity index (D: mean 0.91), and evenness (E: mean 0.78) exhibited high levels of clonal diversity compared to other clonal plants. These results imply that sexual reproduction has played an important role at some time during the history of these populations. Nevertheless, the high level of diversity could have been also partially generated from somatic mutations, although this is unlikely to account for the high diversity generally found among C. grandis genets.     

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.     

Genomic diversity in cattle breeds assessed using polymorphism of intermicrosatellite markers

Using (AG)₉C and (GA)₉C primers and analysis of multilocus profile of intermicrosatellite polymorphism (ISSR-PCR), in 19 cattle breeds, a total of 66 PCR products were revealed. The statistical processing of the data obtained using standard software programs enabled the identification of the ISSR markers, the polymorphism of which made the greatest contribution to the total genomic diversity of both individual breeds and the species Bos taurus in general. These DNA fragments deserve special interest for the further analysis of their localization and functions.     

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.     

Genetic diversity of natural Miscanthus sinensis populations in China revealed by ISSR markers

Miscanthus sinensis is a typical perennial C₄ grass in Asia. In the present study, ISSR markers were used to evaluate the genetic variation within and among grass populations that were sampled from the Zhejiang and Guangdong Provinces of China. Based on nine ISSR primers, 206 clear and reproducible DNA fragments were generated. Relatively low levels of genetic diversity were determined among the populations. The coefficient of genetic differentiation among the populations was 0.52, which indicates that 52.3% of the total molecular variance exists among the populations. Such a high level of divergence present among the populations might be caused by the complex topography and variable climatic conditions present in the two provinces. AMOVA revealed that the majority of the genetic variation was within the populations (50.6%). The application of a novel method, which combines geographical coordinates and genetic differentiation to detect barriers for gene flow, allowed us to identify two zones of lowered gene flow.     

西藏牦牛mtDNA D-loop区的遗传多样性及其遗传分化

通过测定和分析西藏11个牦牛类群114个个体的mtDNA D-loop区全序列,对西藏牦牛的遗传多样性、类群间的亲缘关系及其遗传分化进行了研究。结果表明:①西藏牦牛mtDNA D-loop区全序列长度为890—896 bp,4种核苷酸T、C、A、G的平均比例分别为28.5%、25.3%、32.4%、13.8%,西藏牦牛mtDNA D-loop区富含碱基A+T,表现出一定的碱基偏好性。②共检测到130个变异位点,占分析总位点数的14.33%;其中单一多态位点85个,占多态位点总数的65.38%,简约信息位点45个,占多态位点总数的34.62%。序列变异中碱基缺失、插入和碱基替换等均有,其中碱基替换变异类型中转换114次,颠换12次,在转换变异类型中以A/G、T/C为主,占95.61%,在颠换变异类型中以A/T为主,占75%。③在114个个体中鉴定出90种单倍型,单倍型多样性为0.981±0.008,核苷酸多样性为0.01056±0.00701,均说明西藏牦牛具有丰富的单倍型类型。④90种单倍型分为2个聚类簇(Ⅰ、Ⅱ),聚类簇Ⅰ包含80种单倍型,占全部单倍型的88.89%,涵盖本研究中所有的西藏牦牛类群;聚类簇Ⅱ中有10种单倍型,占单倍型总数的11.11%,涉及的类群有工布江达、帕里、丁青、巴青、江达、类乌齐、桑桑、桑日、斯布,说明西藏牦牛可能有2个母系起源。⑤西藏牦牛类群间核苷酸分歧度(Dxy)在0.503%—1.416%之间,聚类分析和AMOVA分析显示西藏牦牛可分为两大类,康布牦牛、嘉黎牦牛为一类,其余的牦牛类群为另一类。     

SNPs detected in the yak MC4R gene and their association with growth traits

MC4R (melanocortin 4 receptor) is expressed in the appetite-regulating areas of the brain and takes part in leptin signaling pathways. Sequencing of the coding region of the MC4R gene for 354 yaks identified the following five single nucleotide polymorphisms (SNPs): SNP1 (273C>T), SNP2 (321 G>T), SNP3 (864 C>A), SNP4 (1069G>C) and SNP5 (1206 G>C). SNP1, SNP2 and SNP3 were synonymous mutations, whereas SNP4 and SNP5 were missense mutations resulting in amino acid substitutions (V286L and R331S). Pairwise linkage disequilibrium (LD) analysis indicated that two pairs of SNPs, SNP2 and SNP5 (r²=0.81027) and SNP4 and SNP5 (r²=0.53816), exhibited higher degrees of LD. CC genotype of SNP4, CGACG and CTCCC haplotypes for all SNPs were associated with increased BW of animals that were 18 months old and with the average daily gain. The secondary structure and transmembrane region prediction of the yak MC4R protein suggested that SNP4 was correlated with influential changes in the seventh transmembrane domain of the MC4R protein and with the functional deterioration or even incapacitation of MC4R, which may contribute to the increased feed intake, BW and average daily gain of the yaks with CC genotypes. The data from this study suggested that 1069G>C SNP of the MC4R gene could be used in marker-assisted selection of growth traits in the Maiwa yak breed.