Shuli Yang

Whole-genome sequencing of six dog breeds from continuous altitudes reveals adaptation to high-altitude hypoxia

The hypoxic environment imposes severe selective pressure on species living at high altitude. To understand the genetic bases of adaptation to high altitude in dogs, we performed whole-genome sequencing of 60 dogs including five breeds living at continuous altitudes along the Tibetan Plateau from 800 to 5100 m as well as one European breed. More than 150× sequencing coverage for each breed provides us with a comprehensive assessment of the genetic polymorphisms of the dogs, including Tibetan Mastiffs. Comparison of the breeds from different altitudes reveals strong signals of population differentiation at the locus of hypoxia-related genes including endothelial Per-Arnt-Sim (PAS) domain protein 1 (EPAS1) and beta hemoglobin cluster. Notably, four novel nonsynonymous mutations specific to high-altitude dogs are identified at EPAS1, one of which occurred at a quite conserved site in the PAS domain. The association testing between EPAS1 genotypes and blood-related phenotypes on additional high-altitude dogs reveals that the homozygous mutation is associated with decreased blood flow resistance, which may help to improve hemorheologic fitness. Interestingly, EPAS1 was also identified as a selective target in Tibetan highlanders, though no amino acid changes were found. Thus, our results not only indicate parallel evolution of humans and dogs in adaptation to high-altitude hypoxia, but also provide a new opportunity to study the role of EPAS1 in the adaptive processes.

The local origin of the Tibetan pig and additional insights into the origin of Asian pigs.

Background The domestic pig currently indigenous to the Tibetan highlands is supposed to have been introduced during a continuous period of colonization by the ancestors of modern Tibetans. However, there is no direct genetic evidence of either the local origin or exotic migration of the Tibetan pig. Methods and Findings We analyzed mtDNA hypervariable segment I (HVI) variation of 218 individuals from seven Tibetan pig populations and 1,737 reported mtDNA sequences from domestic pigs and wild boars across Asia. The Bayesian consensus tree revealed a main haplogroup M and twelve minor haplogroups, which suggested a large number of small scale in situ domestication episodes. In particular, haplogroups D1 and D6 represented two highly divergent lineages in the Tibetan highlands and Island Southeastern Asia, respectively. Network analysis of haplogroup M further revealed one main subhaplogroup M1 and two minor subhaplogroups M2 and M3. Intriguingly, M2 was mainly distributed in Southeastern Asia, suggesting for a local origin. Similar with haplogroup D6, M3 was mainly restricted in Island Southeastern Asia. This pattern suggested that Island Southeastern Asia, but not Southeastern Asia, might be the center of domestication of the so-called Pacific clade (M3 and D6 here) described in previous studies. Diversity gradient analysis of major subhaplogroup M1 suggested three local origins in Southeastern Asia, the middle and downstream regions of the Yangtze River, and the Tibetan highlands, respectively. Conclusions We identified two new origin centers for domestic pigs in the Tibetan highlands and in the Island Southeastern Asian region.

Genetic diversity and origin of Gayal and cattle in Yunnan revealed by mtDNA control region and SRY gene sequence variation.

There are hump, humpless cattle and gayal distributed in Yunnan province, south-west China, but their genetic background remains unclear. To determine the origin and genetic diversity of Yunnan gayal and cattle (Diqing, Nujiang and Wenshan cattle), we analysed mtDNA control region sequences of 71 samples and SRY gene sequences of 39 samples, together with the available sequences in GenBank. The neighbour-joining phylogeny and the reduced median network analysis showed that Yunnan gayal originated from the hybridization between male Bos frontalis and female Bos taurus or Bos indicus, and that Yunnan cattle mostly originated from B. indicus, also containing some hybrids of male B. indicus and female B. taurus. The phylogenetic pattern of Yunnan cattle was consistent with the recently described cattle matrilineal pool from China and indicated more contribution to the Yunnan cattle from B. indicus than from B. taurus.

Dominant bacterial communities in the rumen of Gayals (Bos frontalis), Yaks (Bos grunniens) and Yunnan Yellow Cattle (Bos taurs) revealed by denaturing gradient gel electrophoresis

The dominant rumen bacteria in Gayals, Yaks and Yunnan Yellow Cattle were investigated using PCR-DGGE approach. The analysis of DGGE profiles, identification of dominant bands and phylogenetic analysis 16S rDNA sequences in DGGE profiles were combined to reveal the dominant bacterial communities and compared the differences between those cattle species. DGGE profiles revealed that Gayals had the most abundant dominant bacteria and the lowest similarity of intraspecies between individuals than other two cattle species. A total of 45 sequences were examined and sequence similarity analysis revealed that Gayals had the most sequences appeared to uncultured bacteria, accounting for 85.0% of the total sequences, Yaks and Yunnan Yellow Cattle had 44.4 and 68.8% uncultured bacterial sequences, respectively. According to phylogenetic analysis, the rumen dominant bacteria of Gayals were mainly phylogenetically placed within phyla firmicutes and bacteroidetes, and the known bacteria were mainly belonged to the genera Lachnospiraceae bacterium, Ruminococcus flavefaciens and Clostridium celerecrescens. Moreover, the dominant bacteria of Yaks were also mainly belonged to phyla firmicutes and bacteroidetes, and the known dominant bacteria were including Ruminococcus flavefaciens, Butyrivibrio fibrisolvens, Pseudobutyrivibrio ruminis, Schwartzia succinivorans and Clostridiales bacterium, most of them are common rumen bacteria. In addition, the dominant bacteria in Yunnan Yellow Cattle were belonged to phyla firmicutes, bacteroidetes and Actinobacteria, and the known dominant bacteria containing Prevotella sp., Staphylococcilentus, Staphylococcus xylosus and Corynebacterium casei. Present study first detected Staphylococcus lentus and Staphylococcus xylosus in the rumen of cattle.

Molecular cloning, sequence identification and tissue expression profile of three novel genes Sfxn1, Snai2 and Cno from Black-boned sheep (Ovis aries)

The complete coding sequences of three of Black-boned sheep (Ovis aries) genes Sfxn1, Snai2 and Cno were amplified using the reverse transcriptase polymerase chain reaction (RT-PCR) according to the conserved sequence information of the cattle or other mammals and known highly homologous sheep ESTs. Black-boned sheep Sfxn1 gene encodes a protein of 322 amino acids which has high homology with the Sfxn1 proteins of five species—cattle 98%, pig 95%, human 95%, rat 93%, and mouse 93%. Black-boned sheep Snai2 gene encodes a protein of 268 amino acids that has high identity with the Snai2 proteins of six species—cattle 99%, pig 94%, human 93%, dog 93%, rat 91%, and mouse 90%. Black-boned sheep Cno gene encodes a protein of 214 amino acids that has high homology with the Cno proteins of four species—cattle 97%, human 75%, mouse 67%, and rat 65%. The phylogenetic tree analysis demonstrated that Black-boned sheep Sfxn1, Snai2 and Cno proteins have close relationship with cattle Sfxn1, Snai2 and Cno proteins. The tissue expression analysis indicated that Black-boned sheep Sfxn1, Snai2 and Cno genes were expressed in a range of tissues including leg muscle, kidney, skin, longissimus dorsi muscle, spleen, heart and liver. Our experiment is the first to provide the primary foundation for further insight into these three sheep genes.

Gene Co-Expression Network Analysis Unraveling Transcriptional Regulation of High-Altitude Adaptation of Tibetan Pig.

Tibetan pigs have survived at high altitude for millennia and they have a suite of adaptive features to tolerate the hypoxic environment. However, the molecular mechanisms underlying the regulation of hypoxia-adaptive phenotypes have not been completely elucidated. In this study, we analyzed differentially expressed genes (DEGs), biological pathways and constructed co-expression regulation networks using whole-transcriptome microarrays from lung tissues of Tibetan and Duroc pigs both at high and low altitude. A total of 3,066 DEGs were identified and this list was over-represented for the ontology terms including metabolic process, catalytic activity, and KEGG pathway including metabolic pathway and PI3K-Akt signaling pathway. The regulatory (RIF) and phenotypic (PIF) impact factor analysis identified several known and several potentially novel regulators of hypoxia adaption, including: IKBKG, KLF6 and RBPJ (RIF1), SF3B1, EFEMP1, HOXB6 and ATF6 (RIF2). These findings provide new details of the regulatory architecture of hypoxia-adaptive genes and also insight into which genes may undergo epigenetic modification for further study in the high-altitude adaptation.

Screening and characterization of a novel ruminal cellulase gene （Umcel-1） from a metagenomic library of gayal （Bos frontalis）

Abstract Gayal is a rare semi-wild bovine species found in the Indo-China. They can graze grasses, including bamboo leaves, as well as reeds and other plant species, and grow to higher mature live weights than Yunnan Yellow cattle maintained in similar harsh environments. The aim of this study was to identify specific cellulase in the gayal rumen. A metagenomic fosmid library was constructed using genomic DNA isolated from the ruminal contents of four adult gayals. This library contained 38 400 clones with an average insert size of 35.5 kb. The Umcel-1 gene was isolated from this library. Investigation of the cellulase activity of 24 random clones led to the identification of the Umcel-1 gene , which exhibited the most potent cellulase activity. Sequencing the Umcel-1 gene revealed that it contained an open reading frame of 942 base pairs that encoded a product of 313 amino acids. The putative gene Umcel-1 product belonged to the glycosyl hydrolase family 5 and showed the highest homology to the cellulase (GenBank accession no. YP_004310852.1) from Clostridium lentocellum DSM 5427, with 44% identity and 62% similarity. The Umcel-1 gene was heterologously expressed in Escherichia coli BL21, and recombinant Umcel-1 was purified. The activity of purified recombinant Umcel-1 was assessed, and the results revealed that it hydrolyzed carboxymethyl cellulose with optimal activity at pH 5.5 and 45°C. To our knowledge, this study provides the first evidence for a cellulase produced by bacteria in gayal rumen.

Natural selection on TMPRSS6 associated with the blunted erythropoiesis and improved blood viscosity in Tibetan Pigs

Tibetan pigs, indigenous to Tibetan plateau, are well adapted to hypoxia. So far, there have been not any definitively described genes and functional sites responsible for hypoxia adaptation for the Tibetan pig. Here we conducted resequencing of the nearly entire genomic region (40.1 kb) of the candidate gene TMPRSS6 (Transmembrane protease, serine 6) associated with hemoglobin concentration (HGB) and red blood cell count (RBC) in 40 domestic pigs and 40 wild boars from five altitudes along the Tea-horse ancient road and identified 708 SNPs, in addition to an indel (CGTG/----) in the intron 10. Both the CGTG deletion frequency and the pairwise r2 linkage disequilibrium showed an increase with elevated altitudes in 838 domestic pigs from five altitudes, suggesting that TMPRSS6 has been under Darwinian positive selection. As the conserved core sequence of hypoxia-response elements (HREs), the deletion of CGTG in Tibetan pigs decreased the expression levels of TMPRSS6 mRNA and protein in the liver revealed by real-time quantitative PCR and western blot, respectively. To explore whether reduced TMPRSS6 expression level could improve blood viscosity, the relationship between CGTG indel and hematologic and hemorheologic parameters in 482 domestic pigs from continuous altitudes was detected and dissected a genetic effect on reducing HGB by 13.25g/L in Gongbo’gyamda Tibetan pigs and decreasing MCV by 4.79 fl in Diqing Tibetan pigs. In conclusion, the CGTG deletion of TMPRSS6 resulted in lower HGB and smaller MCV, thereby blunting erythropoiesis and improving blood viscosity as well as erythrocyte deformability.

Haplotype diversity in mitochondrial DNA reveals the multiple origins of Tibetan horse

The Tibetan horse is a species endemic to the Tibetan plateau, with considerable economic value in the region. However, we currently have little genetic evidence to verify whether the breed originated in Tibet or if it entered the area via an ancient migratory route. In the present study, we analyzed the hypervariable segment I sequences of mitochondrial DNA (mtDNA) in 2,050 horses, including 290 individuals from five Tibetan populations and 1,760 from other areas across Asia. Network analysis revealed multiple maternal lineages in the Tibetan horse. Component analysis of sub-lineage F3 indicated that it decreased in frequency from east to west, a trend reflected both southward and northward from Inner Mongolia. Analysis of population genetics showed that the Deqen horse of eastern Tibet was more closely related to the Ningqiang horse of northern China than to other Tibetan horses or the Yunnan horse. These results indicated that the Tibetan horse migrated first from Central Asia to Mongolia, moved south to eastern Tibet (near Deqen), then finally westward to other regions of Tibet. We also identified a novel lineage K that mainly comprises Tibetan and Yunnan horses, suggesting autochthonous domesticated origin for some Tibetan horse breeds from local wild horses. In conclusion, our study demonstrated that modern Tibetan horse breeds originated from the introgression of local wild horses with exotic domesticated populations outside China.