Lingyang Xu

Analysis of copy number variations in the sheep genome using 50K SNP BeadChip array

BackgroundIn recent years, genome-wide association studies have successfully uncovered single-nucleotide polymorphisms (SNPs) associated with complex traits such as diseases and quantitative phenotypes. These variations account for a small proportion of heritability. With the development of high throughput techniques, abundant submicroscopic structural variations have been found in organisms, of which the main variations are copy number variations (CNVs). Therefore, CNVs are increasingly recognized as an important and abundant source of genetic variation and phenotypic diversity.ResultsAnalyses of CNVs in the genomes of three sheep breeds were performed using the Ovine SNP50 BeadChip array. A total of 238 CNV regions (CNVRs) were identified, including 219 losses, 13 gains, and six with both events (losses and gains), which cover 60.35 Mb of the sheep genomic sequence and correspond to 2.27% of the autosomal genome sequence. The length of the CNVRs on autosomes range from 13.66 kb to 1.30 Mb with a mean size of 253.57 kb, and 75 CNVRs events had a frequency > 3%. Among these CNVRs, 47 CNVRs identified by the PennCNV overlapped with the CNVpartition. Functional analysis indicated that most genes in the CNVRs were significantly enriched for involvement in the environmental response. Furthermore, 10 CNVRs were selected for validation and 6 CNVRs were further experimentally confirmed by qPCR. In addition, there were 57 CNVRs overlapped in our new dataset and other published ruminant CNV studies.ConclusionsIn this study, we firstly constructed a sheep CNV map based on the Ovine SNP50 array. Our results demonstrated the differences of two detection tools and integration of multiple algorithms can enhance the detection of sheep genomic structure variations. Furthermore, our findings would be of help for understanding the sheep genome and provide preliminary foundation for carrying out the CNVs association studies with economically important phenotypes of sheep in the future.

Genome-Wide Association Studies for Growth and Meat Production Traits in Sheep

Background Growth and meat production traits are significant economic traits in sheep. The aim of the study is to identify candidate genes affecting growth and meat production traits at genome level with high throughput single nucleotide polymorphisms (SNP) genotyping technologies. Methodology and Results Using Illumina OvineSNP50 BeadChip, we performed a GWA study in 329 purebred sheep for 11 growth and meat production traits (birth weight, weaning weight, 6-month weight, eye muscle area, fat thickness, pre-weaning gain, post-weaning gain, daily weight gain, height at withers, chest girth, and shin circumference). After quality control, 319 sheep and 48,198 SNPs were analyzed by TASSEL program in a mixed linear model (MLM). 36 significant SNPs were identified for 7 traits, and 10 of them reached genome-wise significance level for post-weaning gain. Gene annotation was implemented with the latest sheep genome Ovis_aries_v3.1 (released October 2012). More than one-third SNPs (14 out of 36) were located within ovine genes, others were located close to ovine genes (878bp-398,165bp apart). The strongest new finding is 5 genes were thought to be the most crucial candidate genes associated with post-weaning gain: s58995.1 was located within the ovine genes MEF2B and RFXANK, OAR3_84073899.1, OAR3_115712045.1 and OAR9_91721507.1 were located within CAMKMT, TRHDE, and RIPK2 respectively. GRM1, POL, MBD5, UBR2, RPL7 and SMC2 were thought to be the important candidate genes affecting post-weaning gain too. Additionally, 25 genes at chromosome-wise significance level were also forecasted to be the promising genes that influencing sheep growth and meat production traits. Conclusions The results will contribute to the similar studies and facilitate the potential utilization of genes involved in growth and meat production traits in sheep in future.

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.

Genetic structure of Chinese indigenous goats and the special geographical structure in the Southwest China as a geographic barrier driving the fragmentation of a large population.

Background China has numerous native domestic goat breeds, however, extensive studies are focused on the genetic diversity within the fewer breeds and limited regions, the population demograogic history and origin of Chinese goats are still unclear. The roles of geographical structure have not been analyzed in Chinese goat domestic process. In this study, the genetic relationships of Chinese indigenous goat populations were evaluated using 30 microsatellite markers. Methodology/Principal Findings Forty Chinese indigenous populations containing 2078 goats were sampled from different geographic regions of China. Moderate genetic diversity at the population level (HS of 0.644) and high population diversity at the species level (HT value of 0.737) were estimated. Significant moderate population differentiation was detected (FST value of 0.129). Significant excess homozygosity (FIS of 0.105) and recent population bottlenecks were detected in thirty-six populations. Neighbour-joining tree, principal components analysis and Bayesian clusters all revealed that Chinese goat populations could be subdivided into at least four genetic clusters: Southwest China, South China, Northwest China and East China. It was observed that the genetic diversity of Northern China goats was highest among these clusters. The results here suggested that the goat populations in Southwest China might be the earliest domestic goats in China. Conclusions/Significance Our results suggested that the current genetic structure of Chinese goats were resulted from the special geographical structure, especially in the Western China, and the Western goat populations had been separated by the geographic structure (Hengduan Mountains and Qinling Mountains-Huaihe River Line) into two clusters: the Southwest and Northwest. It also indicated that the current genetic structure was caused by the geographical origin mainly, in close accordance with the human’s migration history throughout China. This study provides a fundamental genetic profile for the conservation of these populations and better to understand the domestication process and origin of Chinese goats.

Identification of the crucial molecular events during the large-scale myoblast fusion in sheep

It is well known that in sheep most myofibers are formed before birth; however, the crucial myogenic stage and the cellular and molecular mechanisms underpinning phenotypic variation of fetal muscle development remain to be ascertained. We used histological, microarray, and quantitative real-time PCR (qPCR) methods to examine the developmental characteristics of fetal muscle at 70, 85, 100, 120, and 135 days of gestation in sheep. We show that day 100 is an important checkpoint for change in muscle transcriptome and histomorphology in fetal sheep and that the period of 85-100 days is the vital developmental stage for large-scale myoblast fusion. Furthermore, we identified the cis-regulatory motifs for E2F1 or MEF2A in a list of decreasingly or increasingly expressed genes between 85 and 100 days, respectively. Further analysis demonstrated that the mRNA and phosphorylated protein levels of E2F1 and MEF2A significantly declined with myogenic progression in vivo and in vitro. qRT-PCR analysis indicated that PI3K and FST, as targets of E2F1, may be involved in myoblast differentiation and fusion and that downregulation of MEF2A contributes to transition of myofiber types by differential regulation of the target genes involved at the stage of 85-100 days. We clarify for the first time the timing of myofiber proliferation and development during gestation in sheep, which would be beneficial to meat sheep production. Our findings present a repertoire of gene expression in muscle during large-scale myoblast fusion at transcriptome-wide level, which contributes to elucidate the regulatory network of myogenic differentiation.