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.

Y chromosome haplotype diversity of domestic sheep (Ovis aries) in northern Eurasia.

Variation in two SNPs and one microsatellite on the Y chromosome was analyzed in a total of 663 rams representing 59 breeds from a large geographic range in northern Eurasia. SNPA-oY1 showed the highest allele frequency (91.55%) across the breeds, whereas SNPG-oY1 was present in only 56 samples. Combined genotypes established seven haplotypes (H4, H5, H6, H7, H8, H12 and H19). H6 dominated in northern Eurasia, and H8 showed the second-highest frequency. H4, which had been earlier reported to be absent in European breeds, was detected in one European breed (Swiniarka), whereas H7, which had been previously identified to be unique to European breeds, was present in two Chinese breeds (Ninglang Black and Large-tailed Han), one Buryatian (Transbaikal Finewool) and two Russian breeds (North Caucasus Mutton-Wool and Kuibyshev). H12, which had been detected only in Turkish breeds, was also found in Chinese breeds in this work. An overall low level of haplotype diversity (median h = 0.1288) was observed across the breeds with relatively higher median values in breeds from the regions neighboring the Near Eastern domestication center of sheep. H6 is the dominant haplotype in northwestern and eastern China, in which the haplotype distribution could be explained by the historical translocations of the H4 and H8 Y chromosomes to China via the Mongol invasions followed by expansions to northwestern and eastern China. Our findings extend previous results of sheep Y chromosomal genetic variability and indicate probably recent paternal gene flows between sheep breeds from distinct major geographic regions.

The expression and mutation of BMPR1B and its association with litter size in small-tail Han sheep (Ovis aries).

Previous studies have shown that BMPR1B promotes follicular development and ovarian granulosa cell proliferation, thereby affecting ovulation in mammals. In this study, the expression and polymorphism of the BMPR1B gene associated with litter size in small-tail Han (STH) sheep were determined. The expression of BMPR1B was detected in 14 tissues of STH sheep during the follicular phase as well as in the hypothalamic-pituitary-gonadal (HPG) axis of monotocous and polytocous STH sheep during the follicular and luteal phases using quantitative polymerase chain reaction (qPCR). Sequenom MassARRAY® single nucleotide polymorphism (SNP) technology was also used to detect the polymorphism of SNPs in seven sheep breeds. Here, BMPR1B was highly expressed in hypothalamus, ovary, uterus, and oviduct tissue during the follicular phase, and BMPR1B was expressed significantly more in the hypothalamus of polytocous ewes than in monotocous ewes during both the follicular and luteal phases ( P < 0.05 ). For genotyping, we found that genotype and allele frequencies of three loci of the BMPR1B gene were extremely significantly different ( P < 0.01 ) between the monotocous and polytocous groups. Association analysis results showed that the g.29380965A > G locus had significant negative effects on the litter size of STH sheep, and the combination of g.29380965A > G and FecB (Fec - fecundity and B - Booroola; A746G) at the BMPR1B gene showed that the litter size of AG-GG, AA-GG, and GG-GG genotypes was significantly higher compared with other genotypes ( P < 0.05 ). This is the first study to find a new molecular marker affecting litter size and to systematically analyze the expression of BMPR1B in different fecundity and physiological periods of STH sheep.

Novel Insights reveal Anti-microbial Gene Regulation of Piglet Intestine Immune in response to Clostridium perfringens Infection.

LncRNA play important roles in regulation of host immune and inflammation responses in defending bacterial infection. Clostridium perfringens (C. perfringens) type C is one of primary bacteria leading to piglet diarrhea and other intestinal inflammatory diseases. For the differences of host immune capacity, individuals usually show resistance and susceptibility to bacterial infection. However, whether and how lncRNAs involved in modulating host immune resistance have not been reported. We have investigated the expression patterns of ileum lncRNAs of 7-day-old piglets infected by C. perfringens type C through RNA sequencing. A total of 16 lncRNAs and 126 mRNAs were significantly differentially expressed in resistance (IR) and susceptibility (IS) groups. Many lncRNAs and mRNAs were identified to regulate resistance and susceptibility of piglets through immune related pathways. Five lncRNAs may have potential function on regulating the expressions of cytokines, these lncRNAs and cytokines work together to co-regulated piglet immune response to C. perfringens, affecting host resistance and susceptibility. These results provide valuable information for understanding the functions of lncRNA and mRNA in affecting piglet diarrhea resistance of defensing to C. perfringens type C, these lncRNAs and mRNAs may be used as the important biomarkers for decreasing C. perfringens spread and diseases in human and piglets.

The complete mitochondrial genome of the Qinghai Tibetan pig.

The complete mitochondrial genome sequence of the Qinghai Tibetan pig was first determined in this study. The total length of the mitogenome is 16,720 bp. Indicating the an A + T(60.5%)-rich feature, including 2 ribosomal RNA genes, 13 protein-coding genes. 22 transfer RNA genes and 1 non-coding control region. The NJ phylogenetic tree analysis showed that the phylogenetic relationship between Qinghai Tibetan pig and Yimenghei pig was the closest, and the relationship with Chinese northeas wildboar was farthest.

GCH1 plays a role in the high-altitude adaptation of Tibetans.

Tibetans are well adapted to high-altitude hypoxia. Previous genome-wide scans have reported many candidate genes for this adaptation, but only a few have been studied. Here we report on a hypoxia gene ( GCH1, GTP-cyclohydrolase I), involved in maintaining nitric oxide synthetase (NOS) function and normal blood pressure, that harbors many potentially adaptive variants in Tibetans. We resequenced an 80.8 kb fragment covering the entire gene region of GCH1 in 50 unrelated Tibetans. Combined with previously published data, we demonstrated many GCH1 variants showing deep divergence between highlander Tibetans and lowlander Han Chinese. Neutrality tests confirmed a signal of positive Darwinian selection on GCH1 in Tibetans. Moreover, association analysis indicated that the Tibetan version of GCH1 was significantly associated with multiple physiological traits in Tibetans, including blood nitric oxide concentration, blood oxygen saturation, and hemoglobin concentration. Taken together, we propose that GCH1 plays a role in the genetic adaptation of Tibetans to high altitude hypoxia.

EP300 contributes to high-altitude adaptation in Tibetans by regulating nitric oxide production.

The genetic adaptation of Tibetans to high altitude hypoxia likely involves a group of genes in the hypoxic pathway, as suggested by earlier studies. To test the adaptive role of the previously reported candidate gene EP300 (histone acetyltransferase p300), we conducted resequencing of a 108.9 kb gene region of EP300 in 80 unrelated Tibetans. The allele-frequency and haplotype-based neutrality tests detected signals of positive Darwinian selection on EP300 in Tibetans, with a group of variants showing allelic divergence between Tibetans and lowland reference populations, including Han Chinese, Europeans, and Africans. Functional prediction suggested the involvement of multiple EP300 variants in gene expression regulation. More importantly, genetic association tests in 226 Tibetans indicated significant correlation of the adaptive EP300 variants with blood nitric oxide (NO) concentration. Collectively, we propose that EP300 harbors adaptive variants in Tibetans, which might contribute to high-altitude adaptation through regulating NO production.

NADPH oxidase-mediated generation of reactive oxygen species is critically required for survival of undifferentiated human promyelocytic leukemia cell line HL-60.

Nicotinamide adenine dinucleotide phosphate oxidase (NADPH oxidase) mediated generation of reactive oxygen species (ROS) was originally identified as the powerful host defense machinery against microorganism in phagocytes. But recent reports indicated that some non-phagocytic cells also have the NADPH oxidase activity, and the ROS produced by it may act as cell signal molecule. But as far as today, whether the NADPH oxidase also plays similar role in phagocyte has not been paid much attention. Utilizing the undifferentiated HL-60 promyelocytic leukemia cells as a model, the aim of the present study was to determine whether NADPH oxidase plays a role on ROS generation in undifferentiated HL-60, and the ROS mediated by it was essential for cell's survival. For the first time, we verified that the release of ROS in undifferentiated HL-60 was significantly increased by the stimulation with Calcium ionophore or opsonized zymosan, which are known to trigger respiration burst in phagocytes by NADPH oxidase pathway. Diphenylene iodonium (DPI) or apocynin (APO), two inhibitors of NADPH oxidase, significantly suppressed the increasing of ROS caused by opsonized zymosan. Cell survival assay and fluorescence double dyeing with acridine orange and ethidium bromide showed that DPI and APO, as well as superoxide dismutase (SOD) and catalase (CAT) concentration-dependently decreased the viability of undifferentiated HL-60 cells, whereas exogenous H2O2 can rescue the cells from death obviously. Our results suggested that the ROS, generated by NADPH oxidase play an essential role in the survival of undifferentiated HL-60 cells.

Evaluating the immune responses of mice to subcutaneous immunization with Helicobacter pylori urease B subunit.

BACKGROUND: Helicobacter pylori, a gram-negative bacterial pathogen that expresses a strong urease activity, is associated with the development of gastroduodenal disease. Urease B subunit, one of the two structural subunits of urease, was expressed in E. coli BL21 (DE3) strain. The objective of this study was to evaluate the effects of Helicobacter pylori urease B subunit on the immune responses in mice by subcutaneous immunization. METHODS: The mice were immunized and boosted with Helicobacter pylori urease B subunit antigen subcutaneously three times with 2-wk intervals between the immunizations and boosters. The mice in the control group were immunized with PBS. The adjuvant group received PBS containing complete/incomplete freund's adjuvant identical to antigen group without Helicobacter pylori urease B subunit antigen. Four weeks after the final booster, all the mice were sacrificed. Blood was collected on d 0, 14, 28 and 56 before immunization, booster and sacrifice, respectively. Immediately after sacrifice, gastric liquid and spleen were collected for antibody and cytokine analyses. RESULTS: Urease B subunit increased the concentrations of serum and gastric anti-urease B antigen specific IgG, and the levels of interleukin-4 and interferon-γ in splenocytes of the mice (P < 0.05). CONCLUSIONS: This study demonstrated that recombinant urease B subunit can induce systemic and local immune responses in mice by subcutaneous immunization, which might be used as the effective component of vaccine against Helicobacter pylori.