Qianjun Zhao

Mitogenomic meta-analysis identifies two phases of migration in the history of eastern Eurasian sheep

Despite much attention, history of sheep (Ovis aries) evolution, including its dating, demographic trajectory and geographic spread, remains controversial. To address these questions, we generated 45 complete and 875 partial mitogenomic sequences, and performed a meta-analysis of these and published ovine mitochondrial DNA sequences (n = 3,229) across Eurasia. We inferred that O. orientalis and O. musimon share the most recent female ancestor with O. aries at approximately 0.790 Ma (95% CI: 0.637–0.934 Ma) during the Middle Pleistocene, substantially predating the domestication event (∼8–11 ka). By reconstructing historical variations in effective population size, we found evidence of a rapid population increase approximately 20–60 ka, immediately before the Last Glacial Maximum. Analyses of lineage expansions showed two sheep migratory waves at approximately 4.5–6.8 ka (lineages A and B: ∼6.4–6.8 ka; C: ∼4.5 ka) across eastern Eurasia, which could have been influenced by prehistoric West–East commercial trade and deliberate mating of domestic and wild sheep, respectively. A continent-scale examination of lineage diversity and approximate Bayesian computation analyses indicated that the Mongolian Plateau region was a secondary center of dispersal, acting as a “transportation hub” in eastern Eurasia: Sheep from the Middle Eastern domestication center were inferred to have migrated through the Caucasus and Central Asia, and arrived in North and Southwest China (lineages A, B, and C) and the Indian subcontinent (lineages B and C) through this region. Our results provide new insights into sheep domestication, particularly with respect to origins and migrations to and from eastern Eurasia.

Microsatellite analysis revealed genetic diversity and population structure among Chinese cashmere goats

Most cashmere goats are found in northern China and Mongolia. They are regarded as precious resources for their production of high quality natural fibre for the textile industry. It was the first time that the genetic diversity and population structure of nine Chinese cashmere populations has been assessed using 14 ISAG/FAO microsatellite markers. In addition, two Iranian populations and one West African goat population were genotyped for comparison. Results indicated that the genetic diversity of Chinese cashmere goats was rich, but less than those of the Iranian goat populations. All pairwise F(ST) values between the Chinese cashmere goat populations reached a highly significant level (P < 0.001), suggesting that they should all be considered as separate breeds. Finally, clustering analysis divided Chinese cashmere goats into at least two clusters, with the Tibetan Hegu goats alone in one cluster. An extensive admixture was detected among the Chinese goat breeds (except the Hegu), which have important implications for breeding management.

Expression profiling and functional characterization of miR-192 throughout sheep skeletal muscle development

MicroRNAs (miRNAs) are evolutionarily conserved, small, non-coding RNAs that have emerged as key regulators of myogenesis. Here, we examined the miRNA expression profiles of developing sheep skeletal muscle using a deep sequencing approach. We detected 2,396 miRNAs in the sheep skeletal muscle tissues. Of these, miR-192 was found to be up-regulated in prenatal skeletal muscle, but was down-regulated postnatally. MiR-192 expression also decreased during the myogenic differentiation of sheep satellite cells (SCs). MiR-192 overexpression significantly attenuated SCs myogenic differentiation but promoted SCs proliferation, whereas miR-192 inhibition enhanced SCs differentiation but suppressed SCs proliferation. We found that miR-192 targeted retinoblastoma 1 (RB1), a known regulator of myogenesis. Furthermore, knockdown of RB1 in cultured cells significantly inhibited SCs myogenic differentiation but accelerated SCs proliferation, confirming the role of RB1 in myogenesis. Taken together, our findings enrich the ovine miRNA database, and outline the miRNA transcriptome of sheep during skeletal muscle development. Moreover, we show that miR-192 affects SCs proliferation and myogenic differentiation via down-regulation of RB1.

Tracing genetic differentiation of Chinese Mongolian sheep using microsatellites

The genetic consequences of population differentiation and isolation have been the subject of conservation biology. In this study, we analysed the genetic diversity and structure of Mongolian sheep in China. These animals belong to a traditional local breed with high production, extensive adaption, early maturity and roughage resistance. For this purpose, 26 microsatellites were genotyped for five Mongolian sheep populations. The Bayesian clustering indicated five clusters as the most probable genetic structure of the populations investigated. In addition, a clear genetic structure was revealed in three populations distributed at large geographical scales, while the other cluster encompassed UQ and HLBR sheep that displayed no clear differentiation, probably due to their close and small geographical distributions. Overall, our results are helpful in understanding the interplay of population dynamics in these close genetic lineages of Mongolian sheep.

Genomic analysis identified a potential novel molecular mechanism for high-altitude adaptation in sheep at the Himalayas.

Sheep has successfully adapted to the extreme high-altitude Himalayan region. To identify genes underlying such adaptation, we genotyped genome-wide single nucleotide polymorphisms (SNPs) of four major sheep breeds living at different altitudes in Nepal and downloaded SNP array data from additional Asian and Middle East breeds. Using a di value-based genomic comparison between four high-altitude and eight lowland Asian breeds, we discovered the most differentiated variants at the locus of FGF-7 (Keratinocyte growth factor-7), which was previously reported as a good protective candidate for pulmonary injuries. We further found a SNP upstream of FGF-7 that appears to contribute to the divergence signature. First, the SNP occurred at an extremely conserved site. Second, the SNP showed an increasing allele frequency with the elevated altitude in Nepalese sheep. Third, the electrophoretic mobility shift assays (EMSA) analysis using human lung cancer cells revealed the allele-specific DNA-protein interactions. We thus hypothesized that FGF-7 gene potentially enhances lung function by regulating its expression level in high-altitude sheep through altering its binding of specific transcription factors. Especially, FGF-7 gene was not implicated in previous studies of other high-altitude species, suggesting a potential novel adaptive mechanism to high altitude in sheep at the Himalayas.

Population Variation Reveals Independent Selection toward Small Body Size in Chinese Debao Pony

Body size, one of the most important quantitative traits under evolutionary scrutiny, varies considerably among species and among populations within species. Revealing the genetic basis underlying this variation is very important, particularly in humans where there is a close relationship with diseases and in domestic animals as the selective patterns are associated with improvements in production traits. The Debao pony is a horse breed with small body size that is unique to China; however, it is unknown whether the size-related candidate genes identified in Western breeds also account for the small body size of the Debao pony. Here, we compared individual horses from the Debao population with other two Chinese horse populations using single nucleotide polymorphisms (SNPs) identified with the Equine SNP 65 Bead Chip. The previously reported size-related candidate gene HMGA2 showed a significant signature for selection, consistent with its role observed in human populations. More interestingly, we found a candidate gene TBX3, which had not been observed in previous studies on horse body size that displayed the highest differentiation and most significant association, and thus likely is the dominating factor for the small stature of the Debao pony. Further comparison between the Debao pony and other breeds of horses from around the world demonstrated that TBX3 was selected independently in the Debao pony, suggesting that there were multiple origins of small stature in the horse.

Genome-wide detection of copy number variation in Chinese indigenous sheep using an ovine high-density 600 K SNP array

Copy number variants (CNVs) represent a form of genomic structural variation underlying phenotypic diversity. In this study, we used the Illumina Ovine SNP 600 K BeadChip array for genome-wide detection of CNVs in 48 Chinese Tan sheep. A total of 1,296 CNV regions (CNVRs), ranging from 1.2 kb to 2.3 Mb in length, were detected, representing approximately 4.7% of the entire ovine genome (Oar_v3.1). We combined our findings with five existing CNVR reports to generate a composite genome-wide dataset of 4,321 CNVRs, which revealed 556 (43%) novel CNVRs. Subsequently, ten novel CNVRs were randomly chosen for further quantitative real-time PCR (qPCR) confirmation, and eight were successfully validated. Gene functional enrichment revealed that these CNVRs cluster into Gene Ontology (GO) categories of homeobox and embryonic skeletal system morphogenesis. One CNVR overlapping with the homeobox transcription factor DLX3 and previously shown to be associated with curly hair in sheep was identified as the candidate CNV for the special curly fleece phenotype in Tan sheep. We constructed a Chinese indigenous sheep genomic CNV map based on the Illumina Ovine SNP 600 K BeadChip array, providing an important addition to published sheep CNVs, which will be helpful for future investigations of the genomic structural variations underlying traits of interest in sheep.

Detecting signatures of selection within the Tibetan sheep mitochondrial genome

Abstract Tibetan sheep, a Chinese indigenous breed, are mainly distributed in plateau and mountain-valley areas at a terrestrial elevation between 2260 and 4100 m. The herd is genetically distinct from the other domestic sheep and undergoes acclimatization to adapt to the hypoxic environment. To date, whether the mitochondrial DNA modification of Tibetan sheep shares the same feature as the other domestic breed remains unknown. In this study, we compared the whole mitogenome sequences from 32 Tibetan sheep, 22 domestic sheep and 24 commercial sheep to identify the selection signatures of hypoxic-tolerant in Tibetan sheep. Nucleotide diversity analysis using the sliding window method showed that the highest level of nucleotide diversity was observed in the control region with a peak value of π = 0.05215, while the lowest π value was detected in the tRNAs region. qPCR results showed that the relative mtDNA copy number in Tibetan sheep was significantly lower than that in Suffolk sheep. None of the mutations in 12S rRNA were fixed in Tibetan sheep, which indicated that there has been less artificial selection in this herd than the other domestic and commercial breeds. Although one site (1277G) might undergo the purifying selection, it was not identified as the breed-specific allele in Tibetan sheep. We proposed that nature selection was the main drive during the domestication of Tibetan sheep and single mutation (or locus) could not reveal the signature of selection as for the high diversity in the mitogenome of Tibetan sheep.

High occurrence of length heteroplasmy in domestic Bactrian camel (Camelus bactrianus)

Abstract Heteroplasmy is the presence of more than one mitochondrial DNA (mtDNA) variant within a cell, tissue, or individual. In this study, sequence variation was investigated in the control region (CR) of mitochondrial DNA (mtDNA) from 135 individuals belonging to five primary domestic Bactrian camel breeds in China and Mongolia. Due to variation of the repeat unit G(T/C)(AC)n, length heteroplasmy was detected within each camel by direct sequencing and fragment analysis. A high occurrence of mtDNA heteroplasmy, up to 100 percentages was observed in five camel populations. Our study provides the first evidence of extensive length heteroplasmy in Bactrian camels.

Identification of selection signals by large-scale whole-genome resequencing of cashmere goats

Inner Mongolia and Liaoning cashmere goats are two outstanding Chinese multipurpose breeds that adapt well to the semi-arid temperate grassland. These two breeds are characterized by their soft cashmere fibers, thus making them great models to identify genomic regions that are associated with cashmere fiber traits. Whole-genome sequencing of 70 cashmere goats produced more than 5.52 million single-nucleotide polymorphisms and 710,600 short insertions and deletions. Further analysis of these genetic variants showed some population-specific molecular markers for the two cashmere goat breeds that are otherwise phenotypically similar. By analyzing FST and θπ outlier values, we identified 135 genomic regions that were associated with cashmere fiber traits within the cashmere goat populations. These selected genomic regions contained genes, which are potential involved in the production of cashmere fiber, such as FGF5, SGK3, IGFBP7, OXTR, and ROCK1. Gene ontology enrichment analysis of identified short insertions and deletions also showed enrichment in keratinocyte differentiation and epidermal cell differentiation. These findings demonstrate that this genomic resource will facilitate the breeding of cashmere goat and other Capra species in future.