11个绵羊品种MSTN基因非翻译区的变异

利用PCR-RFLP技术对特克塞尔羊、夏洛莱羊、小尾寒羊、蒙古羊、乌珠穆沁羊、阿勒泰羊、呼伦贝尔羊、塔什库尔干羊、多浪羊、湖羊和岗巴羊11个品种的345个个体的肌肉生长抑制素(Myostatin, MSTN)基因非翻译区(UTR)的变异进行了多态性分析。结果表明大小为271 bp和1 003 bp的扩增片段经限制性内切酶MboⅡ和BsaⅠ酶切后表现多态, 经卡方检验所有品种在该基因座位均处于平衡状态(P>0.05), 3种基因型在11个绵羊品种中的分布差异极显著(P<0.01)。通过限制性内切酶HpyCH4Ⅳ 酶切实验, 证明我国9个地方绵羊品种不存在特克塞尔绵羊中发现的导致肌肉发达的SNP位点, 并在3′UTR区发现了个别碱基突变位点能够形成miRNA作用的靶基序, 测序表明3′UTR区的突变频率较高。     

Identification and Characterization of the miRNA Transcriptome of Ovis aries

The discovery and identification of Ovis aries (sheep) miRNAs will further promote the study of miRNA functions and gene regulatory mechanisms. To explore the microRNAome (miRNAome) of sheep in depth, samples were collected that included eight developmental stages: the longissimus dorsi muscles of Texel fetuses at 70, 85, 100, 120, and 135 days, and the longissimus dorsi muscles of Ujumqin fetuses at 70, 85, 100, 120, and 135 d, and lambs at 0 (birth), 35, and 70 d. These samples covered all of the representative periods of Ovis aries growth and development throughout gestation (about 150 d) and 70 d after birth. Texel and Ujumqin libraries were separately subjected to Solexa deep sequencing; 35,700,772 raw reads were obtained overall. We used ACGT101-miR v4.2 to analyze the sequence data. Following meticulous comparisons with mammalian mature miRNAs, precursor hairpins (pre-miRNAs), and the latest sheep genome, we substantially extended the Ovis aries miRNAome. The list of pre-miRNAs was extended to 2,319, expressing 2,914 mature miRNAs. Among those, 1,879 were genome mapped to unique miRNAs, representing 2,436 genome locations, and 1,754 pre-miRNAs were mapped to chromosomes. Furthermore, the Ovis aries miRNAome was processed using an elaborate bioinformatic analysis that examined multiple end sequence variation in miRNAs, precursors, chromosomal localizations, species-specific expressions, and conservative properties. Taken together, this study provides the most comprehensive and accurate exploration of the sheep miRNAome, and draws conclusions about numerous characteristics of Ovis aries miRNAs, including miRNAs and isomiRs.     

Diversity of copy number variation in a worldwide population of sheep

Copy number variation (CNV) represents a major source of genomic variation. We investigated the diversity of CNV distribution using SNP array data collected from a comprehensive collection of geographically dispersed sheep breeds. We identified 24,558 putative CNVs, which can be merged into 619 CNV regions, spanning 197 Mb of total length and corresponding to ~ 6.9% of the sheep genome. Our results reveal a population differentiation in CNV between different geographical areas, including Africa, America, Asia, Southwestern Asia, Central Europe, Northern Europe and Southwestern Europe. We observed clear distinctions in CNV prevalence between diverse groups, possibly reflecting the population history of different sheep breeds. We sought to determine the gene content of CNV, and found several important CNV-overlapping genes (BTG3, PTGS1 and PSPH) which were involved in fetal muscle development, prostaglandin (PG) synthesis, and bone color. Our study generates a comprehensive CNV map, which may contribute to genome annotation in sheep.     

Investigation of Prolific Sheep from UK and Ireland for Evidence on Origin of the Mutations in BMP15 (FecXG,FecXB) and GDF9 (FecGH) in Belclare and Cambridge Sheep

This paper concerns the likely origin of three mutations with large effects on ovulation rate identified in the Belclare and Cambridge sheep breeds; two in the BMP15 gene (FecX^G and FecX^B) and the third (FecG^H) in GDF9. All three mutations segregate in Belclare sheep while one, FecX^B, has not been found in the Cambridge. Both Belclare and Cambridge breeds are relatively recently developed composites that have common ancestry through the use of genetic material from the Finnish Landrace and Lleyn breeds. The development of both composites also involved major contributions from exceptionally prolific ewes screened from flocks in Ireland (Belclare) and Britain (Cambridge) during the 1960s. The objective of the current study was to establish the likely origin of the mutations (FecX^G, FecX^B and FecG^H) through analysis of DNA from Finnish Landrace and Lleyn sheep, and Galway and Texel breeds which contributed to the development of the Belclare breed. Ewes with exceptionally high prolificacy (hyper-prolific ewes) in current flocks on Irish farms were identified to simulate the screening of ewes from Irish flocks in the 1960s. DNA was obtained from: prolific ewes in extant flocks of Lleyn sheep (n = 44) on the Lleyn peninsula in Wales; hyper-prolific ewes (n = 41); prolific Galway (n = 41) ewes; Finnish Landrace (n = 124) and Texel (n = 19) ewes. The FecX^G mutation was identified in Lleyn but not in Finnish Landrace, Galway or Texel sheep; FecX^B was only found among the hyper-prolific ewes. The FecG^H mutation was identified in the sample of Lleyn sheep. It was concluded from these findings that the Lleyn breed was the most likely source of the FecX^G and FecG^H mutations in Belclare and Cambridge sheep and that the FecX^B mutation came from the High Fertility line that was developed using prolific ewes selected from commercial flocks in Ireland in the 1960′s and subsequently used in the genesis of the Belclare.     

Characterization of OAR1 and OAR18 QTL associated with muscle depth in British commercial terminal sire sheep

This study aimed at verifying previously identified QTL affecting growth and carcass traits on ovine chromosome 18 (OAR18) in Texel sheep (n = 1844), and on OAR1 in Charollais (n = 851) and Suffolk (n = 998) sheep. The QTL were investigated using regression and variance component mapping (VCA) of body weight, muscle and fat depth measurements. In addition, the mode of inheritance of the Texel OAR18 QTL was explored, using data from 4376 Texel sheep, fitting VCA models testing for additive and imprinting effects. We also simulated a 480‐sheep population with different QTL imprinting models and various available levels of marker information to understand the behaviour of the VCA results under different assumed genetic models. In summary, the previously identified QTL were successfully verified using both interval mapping and VCA in the three breeds. We propose a polar overdominance mode of inheritance for the OAR18 QTL in Texel sheep, and we present methods to dissect the QTL mode of inheritance, using the Texel OAR18 QTL as an example.     

Polymorphisms of KiSS-1 and GPR54 genes and their relationships with litter size in sheep

The KiSS-1 and GPR54 genes were studied as candidate genes for the prolificacy in sheep. Four pairs of primers were designed to detect single nucleotide polymorphisms of exon 1 of KiSS-1 gene and exon 1, exon 2 and partial exon 5 of GPR54 gene in high fecundity breeds (Small Tail Han and Hu sheep) and low fecundity breeds (Dorset, Texel and Corriedale sheep) by PCR-SSCP. Polymorphisms in exon 1 of KiSS-1 gene were detected in prolific Small Tail Han sheep (AA, AB and BB genotypes) and Hu sheep (AA and CC genotypes), no polymorphism was found in low fecundity sheep breeds (only AA genotype). Polymorphisms in exon 2 of GPR54 gene were detected in prolific Hu sheep (DD and EE genotypes) and no polymorphism was found in prolific Small Tail Han sheep and low fecundity sheep breeds (only DD genotype). No polymorphism was detected in exon 1 and partial exon 5 of GPR54 gene in five sheep breeds. The polymorphic genotypes were sequenced. While compared the BB genotype with the AA genotype, one nucleotide mutation (G1035A) was detected, which resulted in amino acid change, Val25Met. Five nucleotide mutations were detected from AA to CC genotype (C981T, C996T, T997C, C1034G, C1039T), and among them four caused amino acid changes, that is, Arg7Trp, Phe12Leu, Asn24Lys, Ala26Val. While compared the EE genotype with the DD genotype, two nucleotide mutations (T2360C, A2411C) were detected, which gave rise to amino acid changes, Met90Thr and Asp107Ala, respectively. Genotype frequencies of AA, BB and AB were 0.62, 0.05 and 0.33 in Small Tail Han sheep, respectively. The Small Tail Han sheep ewes with genotype BB or AB had 0.88 (P?<?0.05) or 0.51 (P?<?0.05) lambs more than those with genotype AA; the Small Tail Han sheep ewes with genotype BB had 0.37 (P?>?0.05) lambs more than those with genotype AB. These results preliminarily indicated that the KiSS-1 gene may have some association with prolificacy in sheep.     

Genetic diversity of Chinese indigenous sheep breeds inferred from microsatellite markers

To determine the genetic diversity and phylogenetic relationships among Chinese sheep, 10 indigenous breeds and one introduced breed were genotyped for 19 microsatellite loci. The mean number of alleles per breed ranged from 5.44 (Guide Black Fur sheep) to 9.13 (Ujumqin sheep and Hulunbeier sheep), the expected heterozygosity varied from 0.623 (Guide Black Fur sheep) to 0.737 (Zhaotong sheep), and the allelic richness ranged from 5.169 (Guide Black Fur sheep) to 7.610 (Zhaotong sheep). The deviation from Hardy–Weinberg equilibrium (HWE) was statistically significant (P < 0.05) at three loci (SRCRSP5, OarAE129 and DYMS1) in most of the breeds. Chinese sheep breeds had maintained a high level of within-population genetic differentiation (95.23%), with the remainder explained by differentiation among populations (4.77%). The genetic differentiation pattern and genetic relationships among Chinese sheep breeds displayed a high consistency with the traditional classification. Both the Bayesian cluster and principal component analyses showed a reliable clustering pattern, which revealed three major clusters in Chinese indigenous sheep (Mongolian sheep, Kazakh sheep and Tibetan sheep), except Zhaotong and Guide Black Fur sheep. There were probably caused by different breeding history, geography isolation and different levels of inbreeding. This study will help to interpret the genetic characters of Chinese indigenous sheep and benefit to the future conservation programs.     

Genome-Wide Specific Selection in Three Domestic Sheep Breeds

Background

Commercial sheep raised for mutton grow faster than traditional Chinese sheep breeds. Here, we aimed to evaluate genetic selection among three different types of sheep breed: two well-known commercial mutton breeds and one indigenous Chinese breed.

Results

We first combined locus-specific branch lengths and d_i statistical methods to detect candidate regions targeted by selection in the three different populations. The results showed that the genetic distances reached at least medium divergence for each pairwise combination. We found these two methods were highly correlated, and identified many growth-related candidate genes undergoing artificial selection. For production traits, APOBR and FTO are associated with body mass index. For meat traits, ALDOA, STK32B and FAM190A are related to marbling. For reproduction traits, CCNB2 and SLC8A3 affect oocyte development. We also found two well-known genes, GHR (which affects meat production and quality) and EDAR (associated with hair thickness) were associated with German mutton merino sheep. Furthermore, four genes (POL, RPL7, MSL1 and SHISA9) were associated with pre-weaning gain in our previous genome-wide association study.

Conclusions

Our results indicated that combine locus-specific branch lengths and d_i statistical approaches can reduce the searching ranges for specific selection. And we got many credible candidate genes which not only confirm the results of previous reports, but also provide a suite of novel candidate genes in defined breeds to guide hybridization breeding.     

A Transcriptome-Wide Screen for mRNAs Enriched in Fetal Leydig Cells: CRHR1 Agonism Stimulates Rat and Mouse Fetal Testis Steroidogenesis

Fetal testis steroidogenesis plays an important role in the reproductive development of the male fetus. While regulators of certain aspects of steroidogenesis are known, the initial driver of steroidogenesis in the human and rodent fetal testis is unclear. Through comparative analysis of rodent fetal testis microarray datasets, 54 candidate fetal Leydig cell-specific genes were identified. Fetal mouse testis interstitial expression of a subset of these genes with unknown expression (Crhr1, Gramd1b, Itih5, Vgll3, and Vsnl1) was verified by whole-mount in situ hybridization. Among the candidate fetal Leydig cell-specific factors, three receptors (CRHR1, PRLR, and PROKR2) were tested for a steroidogenic function using ex vivo fetal testes treated with receptor agonists (CRH, PRL, and PROK2). While PRL and PROK2 had no effect, CRH, at low (approximately 1 to 10) nM concentration, increased expression of the steroidogenic genes Cyp11a1, Cyp17a1, Scarb1, and Star in GD15 mouse and GD17 rat testes, and in conjunction, testosterone production was increased. Exposure of GD15 fetal mouse testis to a specific CRHR1 antagonist blunted the CRH-induced steroidogenic gene expression and testosterone responses. Similar to ex vivo rodent fetal testes, ≥10 nM CRH exposure of MA-10 Leydig cells increased steroidogenic pathway mRNA and progesterone levels, showing CRH can enhance steroidogenesis by directly targeting Leydig cells. Crh mRNA expression was observed in rodent fetal hypothalamus, and CRH peptide was detected in rodent amniotic fluid. Together, these data provide a resource for discovering factors controlling fetal Leydig cell biology and suggest that CRHR1 activation by CRH stimulates rat and mouse fetal Leydig cell steroidogenesis in vivo.     

New insights into mitogenomic phylogeny and copy number in eight indigenous sheep populations based on the ATP synthase and cytochrome c oxidase genes

The origins and phylogeny of different sheep breeds has been widely studied using polymorphisms within the mitochondrial hypervariable region. However, little is known about the mitochondrial DNA (mtDNA) content and phylogeny based on mtDNA protein-coding genes. In this study, we assessed the phylogeny and copy number of the mtDNA in eight indigenous (population size, n=184) and three introduced (n=66) sheep breeds in China based on five mitochondrial coding genes (COX1, COX2, ATP8, ATP6 and COX3). The mean haplotype and nucleotide diversities were 0.944 and 0.00322, respectively. We identified a correlation between the lineages distribution and the genetic distance, whereby Valley-type Tibetan sheep had a closer genetic relationship with introduced breeds (Dorper, Poll Dorset and Suffolk) than with other indigenous breeds. Similarly, the Median-joining profile of haplotypes revealed the distribution of clusters according to genetic differences. Moreover, copy number analysis based on the five mitochondrial coding genes was affected by the genetic distance combining with genetic phylogeny; we also identified obvious non-synonymous mutations in ATP6 between the different levels of copy number expressions. These results imply that differences in mitogenomic compositions resulting from geographical separation lead to differences in mitochondrial function.     

IGF2 DNA methylation is a modulator of newborn’s fetal growth and development

The insulin-like growth factor 2 (IGF2) gene, located within a cluster of imprinted genes on chromosome 11p15, encodes a fetal and placental growth factor affecting birth weight. DNA methylation variability at the IGF2 gene locus has been previously reported but its consequences on fetal growth and development are still mostly unknown in normal pediatric population. We collected one hundred placenta biopsies from 50 women with corresponding maternal and cord blood samples and measured anthropometric indices, blood pressure and metabolic phenotypes using standardized procedures. IGF2/H19 DNA methylation and IGF2 circulating levels were assessed using sodium bisulfite pyrosequencing and ELISA, respectively. Placental IGF2 (DMR0 and DMR2) DNA methylation levels were correlated with newborn’s fetal growth indices, such as weight, and with maternal IGF2 circulating concentration at the third trimester of pregnancy, whereas H19 (DMR) DNA methylation levels were correlated with IGF2 levels in cord blood. The maternal genotype of a known IGF2/H19 polymorphism (rs2107425) was associated with birth weight. Taken together, we showed that IGF2/H19 epigenotype and genotypes independently account for 31% of the newborn’s weight variance. No association was observed with maternal diabetic status, glucose concentrations or prenatal maternal body mass index. This is the first study showing that DNA methylation at the IGF2/H19 genes locus may act as a modulator of IGF2 newborn’s fetal growth and development within normal range. IGF2/H19 DNA methylation could represent a cornerstone in linking birth weight and fetal metabolic programming of late onset obesity.     

Genome-wide analysis reveals population structure and selection in Chinese indigenous sheep breeds

Background

Traditionally, Chinese indigenous sheep were classified geographically and morphologically into three groups: Mongolian, Kazakh and Tibetan. Herein, we aimed to evaluate the population structure and genome selection among 140 individuals from ten representative Chinese indigenous sheep breeds: Ujimqin, Hu, Tong, Large-Tailed Han and Lop breed (Mongolian group); Duolang and Kazakh (Kazakh group); and Diqing, Plateau-type Tibetan, and Valley-type Tibetan breed (Tibetan group).

Results

We analyzed the population using principal component analysis (PCA), STRUCTURE and a Neighbor-Joining (NJ)-tree. In PCA plot, the Tibetan and Mongolian groups were clustered as expected; however, Duolang and Kazakh (Kazakh group) were segregated. STRUCTURE analyses suggested two subpopulations: one from North China (Kazakh and Mongolian groups) and the other from the Southwest (Tibetan group). In the NJ-tree, the Tibetan group formed an independent branch and the Kazakh and Mongolian groups were mixed. We then used the d_i statistic approach to reveal selection in Chinese indigenous sheep breeds. Among the 599 genome sequence windows analyzed, sixteen (2.7%) exhibited signatures of selection in four or more breeds. We detected three strong selection windows involving three functional genes: RXFP2, PPP1CC and PDGFD. PDGFD, one of the four subfamilies of PDGF, which promotes proliferation and inhibits differentiation of preadipocytes, was significantly selected in fat type breeds by the Rsb (across pairs of populations) approach. Two consecutive selection regions in Duolang sheep were obviously different to other breeds. One region was in OAR2 including three genes (NPR2, SPAG8 and HINT2) the influence growth traits. The other region was in OAR 6 including four genes (PKD2, SPP1, MEPE, and IBSP) associated with a milk production quantitative trait locus. We also identified known candidate genes such as BMPR1B, MSRB3, and three genes (KIT, MC1R, and FRY) that influence lambing percentage, ear size and coat phenotypes, respectively.

Conclusions

Based on the results presented here, we propose that Chinese native sheep can be divided into two genetic groups: the thin type (Tibetan group), and the fat type (Mongolian and Kazakh group). We also identified important genes that drive valuable phenotypes in Chinese indigenous sheep, especially PDGFD, which may influence fat deposition in fat type sheep.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1384-9) contains supplementary material, which is available to authorized users.     

Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds

A dataset consisting of 787 animals with high‐density SNP chip genotypes (346 774 SNPs) and 939 animals with medium‐density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between F_PED and F_ROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean F_ROH estimates from medium‐density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make‐up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds.