Mapping of the transcription start site (TSS) and identification of SNPs in the bovine neuropeptide Y (NPY) gene

Background

The domestic goat is one of the important livestock species of India. In the present study we assess genetic diversity of Indian goats using 17 microsatellite markers. Breeds were sampled from their natural habitat, covering different agroclimatic zones.

Results

The mean number of alleles per locus (NA) ranged from 8.1 in Barbari to 9.7 in Jakhrana goats. The mean expected heterozygosity (He) ranged from 0.739 in Barbari to 0.783 in Jakhrana goats. Deviations from Hardy-Weinberg Equilibrium (HWE) were statistically significant (P < 0.05) for 5 loci breed combinations. The D_A measure of genetic distance between pairs of breeds indicated that the lowest distance was between Marwari and Sirohi (0.135). The highest distance was between Pashmina and Black Bengal. An analysis of molecular variance indicated that 6.59% of variance exists among the Indian goat breeds. Both a phylogenetic tree and Principal Component Analysis showed the distribution of breeds in two major clusters with respect to their geographic distribution.

Conclusion

Our study concludes that Indian goat populations can be classified into distinct genetic groups or breeds based on the microsatellites as well as mtDNA information.     

Polymorphism of β-lacto globulin gene in Indian goats and its effect on milk yield

Polymorphism in the beta-LG gene in the Indian goat was investigated by the SDS-PAGE and PCR-RFLP method. SDS-PAGE was carried out in 1098 samples belonging to 8 different breeds of Indian goats. The electrophoretic pattern in the beta-LG locus showed the presence of AA and AB genotypes with frequencies of 0.81 and 0.19, 0.89 and 0.11, 0.50 and 0.50, 0.80 and 0.20, 0.84 and 0.16, 1.00 and 0.00, 0.98 and 0.02 and 0.950 and 0.050 in Jamunapari, Barbari, Marwari, Sirohi, Jakhrana, Beetal, Local UP and Local MP goats, respectively. A total of 358 individuals belonging to 13 different genetic groups were analyzed by the PCR-RFLP method. The amplified product was observed as 426 bp and the restriction digestion with SacII revealed three genotypes, namely S1S1, S1S2 and S2S2 at the beta-LG locus. The frequency of the S2S2 genotype ranged from 0.42 to 1.00 in the population. An analysis was carried out in Jamunapari and Barbari goats to observe the effect of the beta-LG genotype on 90-day milk yield. Least squares analysis of data showed that beta-LG AA animals had higher milk yield than the beta-LG AB genotype in both breeds (P < 0.01).     

Evaluation of adaptability to different seasons in goat breeds of semi-arid region in India through differential expression pattern of heat shock protein genes

The present study was conducted to examine differential expression pattern of HSP genes and adaptability in Indian goat breeds of semi-arid region. The study was conducted in five animals from each breed viz. Barbari, Sirohi, and Jhakrana during winter, thermo-neutral and summer seasons. The respiratory rate (RR) and rectal temperature (RT) of the goats were recorded at 09:00 h during the study period. The blood samples were collected for RNA isolation, cDNA synthesis, and quantitative analysis of HSP genes expression by quantitative RT-PCR. The RR increased significantly (p < 0.01) during summer as compared to winter and thermo-neutral season however, RT did not change (p > 0.05) during different seasons. The expression of HSP genes was significantly (p < 0.01) increased during summer (high THI) as compared to thermo-neutral season in all the goat breeds. Among HSPs, only HSP90 was upregulated (p < 0.01) in Jhakrana goats during winter as compared to thermo-neutral season. The deviation in expression of HSP genes during summer and winter with respect to thermo-neutral season was minimum in Barbari goats. Therefore, it can be concluded that Barbari goats possessed better adaptability during summer and winter as compared to Sirohi and Jhakrana goats in semi-arid climatic conditions of India.     

Analysis of 22 heterologous microsatellite markers for genetic variability in Indian goats

The genetic variability of 22 heterologous microsatellite markers was analyzed in two Indian goat breeds, namely Bengal and Chegu. The heterozygosity, polymorphism information content (PIC), and probability of identity of two individuals were calculated for all microsatellite loci in both the breeds. The observed number of alleles varied between 4 and 13 at the studied microsatellite loci. The evaluated microsatellite loci exhibited high mean heterozygosity of 0.69 +/- 0.11 and 0.66 +/- 0.07 in Bengal and Chegu goats, respectively. The mean PIC values of the studied loci in these breeds were 0.79 +/- 0.08 and 0.78 +/- 0.05, respectively. The probability of identity of two random individuals from different breeds, taking into account, all the 22 microsatellite loci was as low as 5.523 x 10(-40). On the basis of these results, we propose that these microsatellite markers may be used with reliability for studying genetic diversity and for identification of individuals in Indian goat breeds.     

Genetic Structuring of Nine Indian Domestic Goat Breeds Based on SNPs Identified in IGF-1 Gene

The caprine Insulin like Growth Factor1 (IGF1) gene was analyzed for identification of single nucleotide polymorphisms (SNPs) and genetic structuring of Indian goat breeds. A panel of 80 samples belonging to nine Indian goat breeds (Capra hircus) including three large sized breeds (Jamunapari, Beetal and Jakhrana); three medium sized breeds (Sirohi, Barbari, and Osmanabadi) and three small sized breeds (Black Bengal, Changthangi, and Gaddi) were screened for SNP identification and diversity analysis. The comparative gene sequence analysis of all the nine goat breeds studied revealed a total of 18 SNPs in IGF1 gene. All the nucleotide changes were found to be synonymous. The mean observed heterozygosity was found to be maximum (0.074) in Sirohi, Beetal, Osmanabadi, and Gaddi breeds of goat, whereas it is found to be minimum (0.019) in Black Bengal breed of goat. The rest of the breeds were intermediate in terms of heterozygosity. The same has been confirmed by allele frequency distribution across the studied loci. Barbari and Gaddi were found to be more differentiated (0.0123), Changthangi and Jamunapari were least differentiated (0.00110) based on Nei's genetic distance.     

The analysis of genetic diversity and differentiation of six Chinese cattle populations using microsatellite markers

A total of 321 individuals from six cattle populations of four species in a bovine subfamily in China were studied using 12 pairs of microsatellite markers. The genetic diversities within and between populations were calculated. The phylogenetic trees were constructed by （δμ）^2 and DA distances, and the divergence times between populations were estimated by （δμ）^2. Altogether, 144 microsatellite alleles were detected including 24 private alleles and nine shared alleles. Chinese Holstein had the largest number of private alleles （10）, whereas Bohai black and Buffalo had the smallest number of private alleles （2）. Chinese Holstein showed the highest genetic variability. Its observed number of alleles （Na）, mean effective number of alleles （MNA）, and mean heterozygosity （He） were 7.7500, 4.9722, and 0.7719, respectively, whereas, the Buffalo and Yak showed low genetic variability. In the phylogenetic trees, Luxi and Holstein grouped first, followed by Bohai and Minnan. Yak branched next and buffalo emerged as the most divergent population from other cattle populations. Luxi and Bohai were estimated to have diverged 0.039-0.105 million years ago （MYA）, however, buffalo and Holstein diverged 0.501-1.337 MYA. The divergence time of Yak versus Minnan, Holstein and buffalo was 0.136-0.363, 0.273-0.729, and 0.326-0.600 MYA, respectively.     

Rapid genetic divergence in postglacial populations of threespine stickleback (Gasterosteus aculeatus): the role of habitat type, drainage and geographical proximity

The goal of this study was to evaluate the role of common ancestry, and of geographical or reproductive isolation, in genetic divergence in populations of threespine sticklebacks (Gasterosteus aculeatus). Using seven DNA microsatellite loci we compared the effects of habitat type, drainage system and geographical proximity on genetic distance among 16 populations situated in an area in Schleswig-Holstein (Germany) that became deglaciated approximately 12 000 years ago. Stickleback population structure correlated only weakly with drainage system, whereas the primary divergence was among habitat types. Phylogenetic analysis revealed that lake (n = 7) and river (n = 5) populations formed two distinct clades (Cavalli-Sforza's and Edwards' chord distance, 82-100% bootstrap support) at approximately equal genetic distances to a third clade, comprising putative estuarine (n = 4) ancestors. Allele frequencies in lake and river populations represented different subsets of the genetically more diverse estuarine populations. In nested amovas approximately twice the genetic variance was distributed among lake vs. river vs. estuarine populations as compared with the combined effects of drainage system and geographical distance. Limited gene flow between habitat types must have been established after postglacial colonization, suggesting ecological hybrid inferiority or behavioural mating barriers between ecotypes. Within estuarine and lake populations, population differentiation followed an isolation-by-distance model. Given the high observed heterozygosities within the 16 study populations (HO = 0.65-0.87), the mean divergence between lake and river population pairs (FST = 0.18 +/- 0.007) would be reached after 300-6000 generations in a stepwise mutation model, depending on the size of N(e). This demonstrates both the utility of hypervariable microsatellites for detecting recent population divergences and the danger of operating at temporal or spatial scales which are beyond their resolution.     

中国部分家养山羊mtDNA D环区遗传多样性与进化

对我国10个家养山羊品种140只个体的mtDNA D-loop区进行测序分析, 测得结果: 整个D-loop区为 1 211~1 213 bp, 共检测到84种单倍型, 171个多态位点。其中核苷酸多样性(Nucleotide diversity, Pi)为: 0.02063 ± 0.00225, 单倍型多样性(Haplotype gene diversity, Hd)为: 0.988 ± 0.003, 平均核苷酸差异数(Average number of nucleotide differences) k为: 24.896。表明我国家养山羊品种遗传多样性丰富。通过构建NJ网络进化树, 得出中国家养山羊主要分为两大支系, 并且其中一支与角骨羊(Capra aegagrus)聚在一起, 而旋角野山羊(Capra fal-coneri)单独聚为一支, 说明角骨羊对中国家养山羊贡献较大。     

Cross-species amplification of Bovidae microsatellites and low diversity of the endangered Korean goral

The Korean goral (Nemorhaedus caudatus) is an endangered species of wild goat. The conservation and management of this species could benefit from a better understanding of its genetic diversity and structure. Cross-species amplification of 34 Bovidae microsatellite loci was tested on a panel of 6 Korean gorals and 10 domestic goats. After polymerase chain reaction (PCR) optimization, 29 (85.3%) microsatellite loci amplified successfully for the Korean gorals and 27 (79.4%) for the domestic goats. Of the amplified products, 16 (55.2%) were polymorphic in the Korean goral and 22 (81.5%) in domestic goats. Nei's unbiased mean heterozygosity and mean allele number per locus were, respectively, 0.356 and 2.6 in the Korean goral and 0.636 and 4.8 in domestic goats. Low genetic diversity in the Korean gorals observed in this preliminary microsatellite survey suggests an urgent need for further detailed study of genetic diversity in Korean goral populations and a population management strategy based on these studies. Current results of cross-species amplification of domestic Bovidae microsatellites could be employed for the necessary population genetic studies on the Korean goral and other endangered Caprinae species.     

Transcriptomic diversity in longissimus thoracis muscles of Barbari and Changthangi goat breeds of India

The present study is an attempt to examine the differential expression of genes in longissimus thoracis muscles between meat and wool type Indian goat breeds. Barbari goat is considered the best meat breed while Changthangi is famous for its fine fibre quality. RNA sequencing data was generated from four biological replicates of longissimus thoracis muscles of Barbari and Changthangi goats. A clear demarcation could be observed between the breeds in terms of expression of genes associated with lipid metabolism (FASN, SCD, THRSP, DGAT2 and FABP3). Most significant genes with high connectivity identified by gene co-expression network analysis were associated with triacylglycerol biosynthesis pathway in Barbari goat. Highly interactive genes identified in Changthangi goat were mainly associated with muscle fibre type. This study provides an insight into the differential expression of genes in longissimus thoracis muscles between Barbari and Changthangi goats that are adapted to and reared in different agro-climatic regions.     

Genomic divergences between humans and other hominoids and the effective population size of the common ancestor of humans and chimpanzees

To study the genomic divergences among hominoids and to estimate the effective population size of the common ancestor of humans and chimpanzees, we selected 53 autosomal intergenic nonrepetitive DNA segments from the human genome and sequenced them in a human, a chimpanzee, a gorilla, and an orangutan. The average sequence divergence was only 1.24% +/- 0.07% for the human-chimpanzee pair, 1.62% +/- 0.08% for the human-gorilla pair, and 1.63% +/- 0.08% for the chimpanzee-gorilla pair. These estimates, which were confirmed by additional data from GenBank, are substantially lower than previous ones, which included repetitive sequences and might have been based on less-accurate sequence data. The average sequence divergences between orangutans and humans, chimpanzees, and gorillas were 3.08% +/- 0.11%, 3.12% +/- 0.11%, and 3.09% +/- 0.11%, respectively, which also are substantially lower than previous estimates. The sequence divergences in other regions between hominoids were estimated from extensive data in GenBank and the literature, and Alus showed the highest divergence, followed in order by Y-linked noncoding regions, pseudogenes, autosomal intergenic regions, X-linked noncoding regions, synonymous sites, introns, and nonsynonymous sites. The neighbor-joining tree derived from the concatenated sequence of the 53 segments--24,234 bp in length--supports the Homo-Pan clade with a 100% bootstrap value. However, when each segment is analyzed separately, 22 of the 53 segments (approximately 42%) give a tree that is incongruent with the species tree, suggesting a large effective population size (N(e)) of the common ancestor of Homo and Pan. Indeed, a parsimony analysis of the 53 segments and 37 protein-coding genes leads to an estimate of N(e) = 52,000 to 96,000. As this estimate is 5 to 9 times larger than the long-term effective population size of humans (approximately 10,000) estimated from various genetic polymorphism data, the human lineage apparently had experienced a large reduction in effective population size after its separation from the chimpanzee lineage. Our analysis assumes a molecular clock, which is in fact supported by the sequence data used. Taking the orangutan speciation date as 12 to 16 million years ago, we obtain an estimate of 4.6 to 6.2 million years for the Homo-Pan divergence and an estimate of 6.2 to 8.4 million years for the gorilla speciation date, suggesting that the gorilla lineage branched off 1.6 to 2.2 million years earlier than did the human-chimpanzee divergence.     

Post-ice age recolonization and differentiation of Fucus serratus L. (Phaeophyceae; Fucaceae) populations in Northern Europe

The seaweed Fucus serratus is hypothesized to have evolved in the North Atlantic and present populations are thought to reflect recolonization from a southern refugium since the last glacial maximum 18 000-20 000 years bp. We examined genetic structure across several spatial scales by analysing seven microsatellite loci in populations collected from 21 localities throughout the species' range. Spatial auto-correlation analysis of seven microsatellite loci revealed no evidence for spatial clustering of alleles on a scale of 100 m despite limited gamete dispersal in F. serratus of approximately 2 m from parental individuals. Pairwise theta analysis suggested that the minimal panmictic unit for F. serratus was between 0.5 and 2 km. Isolation by distance was significant along some contiguous coastlines. Population differentiation was strong within the Skagerrak-Kattegat-Baltic Seas (SKB) (global theta= 0.17) despite a short history of approximately 7500 years. A neighbour-joining tree based on Reynold's distances computed from the microsatellite data revealed a central assemblage of populations on the Brittany Peninsula surrounded by four well-supported clusters consisting of the SKB, the North Sea (Ireland, Helgoland), and two populations from the northern Spanish coast. Samples from Iceland and Nova Scotia were most closely aligned with northwest Sweden and Brittany, respectively. When sample sizes were standardized (N = 41), allelic diversity was twofold higher for Brittany populations than for populations to the north and threefold higher than southern populations. The Brittany region may be a refugium or a recolonized area, whereas the Spanish populations most likely reflect present-day edge populations that have undergone repeated bottlenecks as a consequence of thermally induced cycles of recolonization and extinction.     

微卫星DNA标记探讨镜鲤的种群结构与遗传变异

采用30个微卫星分子标记, 对5个镜鲤群体的观测杂合度(^H_o)、期望杂合度(^H_e)、多态信息含量(PIC)和有效等位基因数(^A_e)等进行了遗传检测, 根据基因频率计算遗传相似系数和Nei氏标准遗传距离, 以c2检验估计Hardy-Weinberg平衡, 以近交系数(FST)和基因流(Nm)分析群体的遗传分化。同时, 使用PHYLIP3.63软件绘制基于Nei氏标准遗传距离的UPGMA聚类图, 并进行bootstrap自举检验验证进化树的可靠性。在德国镜鲤选育系(Scattered Cyprinus carpio L.)和来自4个不同养殖场(松浦、东岗、奉城和辽中)的德国镜鲤群体中共检测到7 083个扩增片段, 长度在102 ~ 446 bp之间, 在群体内扩增出等位基因1~16个不等, 共计356个等位基因。结果表明: (1)5个群体检测的有效等位基因数在1.07~12.30个不等, 平均多态信息含量为0.74、0.74、0.69、0.75和0.75, 无偏期望杂合度的平均值为0.74、0.78、0.70、0.76和0.78, 说明这几个群体属于高度多态, 遗传多样性水平较高。(2)群体间相似系数在0.52以上, 相似性较高。聚类分析显示, 东岗、奉城和辽中3个养殖场的德国镜鲤群体聚类成一个分支, 而德国镜鲤选育系与松浦群体聚类成另一分支。聚类的先后与它们在地理分布上距离远近有一定的相关性。(3)在与功能基因相关的多个微卫星基因座位上, 扩增产物呈现不同程度的缺失现象, 这些无效等位基因的产生可能与结构基因在育种中受到人工选择的影响较大有关。     

Use of Chronogest implants and Folltropin for estrus synchronization and superovulation in goats.

Six Barbari goats each were assigned randomly to treatments 1,2 or 3, comprising im injections of FSH (folltropin) at 12, 14 or 16 mg dose level respectively. Estrus was synchronized with intravaginal sponge impregnated with flugestone acetate (30 mg; chronogest) inserted for 12 days and cloprostenol (125 micrograms) im at the insertion as well as at removal of sponge. FSH treatment started 48 hr before the sponge removal as 4-day declining dose scheme. Estrus could be effectively synchronized in all goats under the study, with significant difference (P less than 0.05) in the onset of estrus between the treatment groups. All goats were administered with 750 IU hCG i.v. at estrus. Recording of ovarian response and embryo recovery was done 45 hr after the onset of estrus. The prime aim of superovulation was effectively achieved in Barbari goats with the use of chronogest implants and folltropin. There was no difference (P greater than 0.05) between the treatment groups in recovery of transferable embryos, however, 14 mg folltropin appeared to be near optimal dose. There was no adverse effect on the quality of recovered embryos with high doses of folltropin.     

异源四倍体鲫鲤及其原始亲本遗传变异的微卫星标记分析

采用从鲤中分离出来的32对微卫星DNA标记,对异源四倍体鲫鲤、红鲫和野鲤的基因组DNA进行了研究。在筛选出的15对微卫星引物中,随引物不同,各等位基因数为1～8个,大小在100～420bp之间。从3个不同群体内部的遗传相似系数来看,异源四倍体鲫鲤个体之间的遗传相似系数最大,说明异源四倍体鲫鲤群体内部的遗传变异程度最低,已经形成了一个遗传性状稳定的群体。从3个不同群体之间的遗传相似系数来看,异源四倍体鲫鲤和红鲫遗传相似系数为0．5625,和野鲤的遗传相似系数为0．5125,说明异源四倍体鲫鲤接受原始母本的遗传物质比原始父本野鲤要多一些。微卫星标记与以前报道的RAPD标记的检测结果是相似的,然而由微卫星标记获得的种群内和种群间的遗传距离均大于RAPD,说明微卫星标记比RAPD标记显示出更高的个体多态性。     

Phylogeography of the rare Balkan endemic Martino's vole, Dinaromys bogdanovi, reveals strong differentiation within the western Balkan Peninsula

The spatial genetic structure of Martino's vole, a rare palaeoendemic species of the western Balkans, was investigated using DNA isolated from archived museum samples. The study was based on partial sequencing (555 bp) of the mitochondrial cytochrome b gene for 63 specimens from 20 different localities throughout the species' range. Three highly divergent allopatric phylogenetic lineages (Northwestern, Central and Southeastern) were recognized among 47 haplotypes, suggesting three independent glacial differentiation centres within the western Balkans. The Northwestern lineage, which showed the highest divergence from all other samples (mean sequence divergence of 6.64% +/- 1.10), comprised samples collected from northwest of the Neretva River in Croatia, western Bosnia and Herzegovina. The Central and Southeastern lineages were separated by an average sequence divergence of 2.95% +/- 0.66 and were geographically divided by the Drim River (the Kosovo Basin in Serbia). Overall, haplotype diversity decreased from the Northwestern lineage to the Central and subsequently the Southeastern lineage, in a geographical pattern consistent with a stepping stone colonization. The observed distribution indicates a gradual southerly expansion with subsequent allopatry across the Neretva River and Drim River approximately 1 and 0.3 million years ago, respectively. Such a scenario is concordant with palaeontological evidence. Several highly divergent sublineages within the Northwestern and Central lineages showed no significant geographical structuring, suggesting secondary contact of allopatrically evolved lineages. We hypothesize that the topographical complexity of the Balkans promoted allopatry and isolation on a small geographical scale during interglacial periods, with secondary contact during glacial maxima. Furthermore, the three main lineages should be regarded as evolutionary significant units with important implications for conservation. Ecological data show that the Northwestern lineage in particular fulfils all criteria for a highly endangered, evolutionarily significant unit.     

Conservation within artiodactyls of an AATA interrupt in the IGF-I microsatellite for 19-35 million years

Occurrence of an AATA interrupt in the IGF-I microsatellite was investigated in a number of Artiodactyl species, namely pigs, camels, deer, cattle, goats, and sheep. Comparison of DNA sequences in the 5' flank of the microsatellite in these species revealed that the interrupt within the microsatellite is conserved in deer, cattle, sheep, and goats but is absent from camels and pigs. The interrupt was introduced into the Artiodactyl phylogeny after the divergence of the Camelidae but before the divergence of the Cervidae, and thus its time of origin can be estimated to be 19-35 MYA. In contrast to the repeat units which are hypermutable, the interrupt has been conserved for a very long time and may even have suppressed microsatellite variation by inhibiting replication slippage. A 12-bp deletion in the 5' flank of the microsatellite in camels corresponds to a consensus reversed repeat in deer, cattle, sheep, and goats with unknown functional significance. Apart from this deletion, the 5' flank of the microsatellite is highly conserved in Artiodactyl species.     

Genetic divergence of Far Eastern dace species belonging to the genus Tribolodon (Pisces,Cyprinidae) and closely related taxa

Based on a biochemical-genetic approach, heterozygosity and divergence of structural genes of 30 enzyme loci were analyzed in six dace species. In addition, intra- and interspecific divergence of gene expression was analyzed based on a sample of 12 to 15 loci. Mean heterozygosities per individual varied as follows: Tribolodon species, Hobs = 0.007 +/- 0.007 and Hexp = 0.007 +/- 0.007; T. ezoe, Hobs = 0.045 +/- 0.016 and Hexp = 0.067 +/- 0.029. Several variants of genetic distances were estimated. Standard Nei's distances (DN) varied from 0.145 to 0.284 in four dace species studied. As related to Tribolodon dace species, the following genetic distances were obtained for two members of other genera: Pseudaspius leptocephalus, DN = 0.269; Leuciscus waleckii, DN = 0.769. Based on the distance matrices, different clustering algorithms were realized. The main feature shared by different dendrograms was a separate position of the cluster joining Far-Eastern dace species, to which P. leptocephalus and L. waleckii are successively added. Among the species studied, the proportion of loci similar by expression (E) varied from 87 to 100%. The greatest difference was found between anadromous and nonanadromous ecotypes of T. hakonensis, E = 67%. The following conclusions can be made: (1) Four studied species of the genus Tribolodon are rather well genetically differentiated. Diagnostic loci are available. (2) A nominal dace species, T. species, should be considered the fourth isolated species of this genus, which is confirmed by its recent zoological acceptance of this species. (3) The origin and divergence of dace species belonging to the genus Tribolodon are relatively late (1 to 3 Myr ago) historical events. (4) Taxonomically, the genus Tribolodon belong to the tribe Pseudaspinini together with P. leptocephalus, which is confirmed by genetic data. (5) Data on heterozygosity and the divergence of structural and regulatory elements of genome, along with the proposed scheme of speciation types, suggest the following speciation modes for the species studied: for four species, adaptive divergence and for two species, genetic transformation.