Global climate change has a significant effect on extreme environments and a profound influence on species survival. However, little is known of the genome-wide pattern of livestock adaptations to extreme environments over a short time frame following domestication. Sheep (Ovis aries) have become well adapted to a diverse range of agroecological zones, including certain extreme environments (e.g., plateaus and deserts), during their post-domestication (approximately 8–9 kya) migration and differentiation. Here, we generated whole-genome sequences from 77 native sheep, with an average effective sequencing depth of ∼5× for 75 samples and ∼42× for 2 samples. Comparative genomic analyses among sheep in contrasting environments, that is, plateau (>4,000 m above sea level) versus lowland (<100 m), high-altitude region (>1500 m) versus low-altitude region (<1300 m), desert (<10 mm average annual precipitation) versus highly humid region (>600 mm), and arid zone (<400 mm) versus humid zone (>400 mm), detected a novel set of candidate genes as well as pathways and GO categories that are putatively associated with hypoxia responses at high altitudes and water reabsorption in arid environments. In addition, candidate genes and GO terms functionally related to energy metabolism and body size variations were identified. This study offers novel insights into rapid genomic adaptations to extreme environments in sheep and other animals, and provides a valuable resource for future research on livestock breeding in response to climate change.

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.

SUMMARYInterleukin 17A (IL-17A) is implicated in the pathogenesis of several neuroimmunologic diseases. We aimed to evaluate the relationship between IL-17A and seizure severity in patients with epilepsy. Seventy patients with probable symptomatic epilepsy and 68 healthy controls were included. Interictal serum IL-17A and related cytokine (IL-23, IL-6, IL-1b, interferon gamma (IFN-c), and IL-10) levels were measured. The relationship between seizure severity and cytokine concentrations was assessed by Spearman correlation and multivariate linear regression test. IL-17A levels in the cerebrospinal fluid (CSF) were tested in 30 additional patients with epilepsy, either in the postictal or interictal period and 15 patients with idiopathic inflammatory demyelinating diseases (IIDDs). Interictal serum IL-17A levels were significantly elevated in patients with epilepsy compared to controls. IL-6, IFN-c, and IL-1b levels were also markedly elevated. Spearman correlation analysis revealed positive correlation between IL-17A, IL-6 levels and Veterans Administration Seizures Frequency and Severity Rating Scale score (VA score); IFN-c, IL-10 levels, and National Hospital Seizure Severity Scale (NHS3) score. In addition, IL-17A levels correlated significantly with seizure frequency. Multivariate linear regression test showed that only IL-17A levels were independently positively correlated with VA scores (B = 0.288, p = 0.027). Postictal IL-17A levels in the CSF were significantly elevated compared to interictal patients and patients with IIDDs. Our results suggest that interictal IL-17A levels correlated highly with seizure severity.

PGES occurs more commonly after GTCS in sleep than in wakefulness but, in the epilepsy-monitoring unit (EMU), a patient's state at seizure onset does not affect the degree of respiratory dysfunction or duration of postictal immobility. In sleep, outside the hospital setting, GTCS are likely to go unnoticed. Postictal immobility in prone patients prevents head repositioning and unimpeded air exchange. A positive feedback cycle ensues with increasing respiratory distress, potentiating postictal immobility and PGES and eventually leading to asystole. Our findings suggest that the high incidence of nocturnal SUDEP may be related to the unsupervised environment during sleep rather than the severity of sleep-related respiratory dysfunction or PGES duration in the immediate postictal period.

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.

Genome-wide association studies (GWASs) have been widely applied in livestock to identify genes associated with traits of economic interest. Here, we conducted the first GWAS of the supernumerary nipple phenotype in Wadi sheep, a native Chinese sheep breed, based on Ovine Infinium HD SNP BeadChip genotypes in a total of 144 ewes (75 cases with four teats, including two normal and two supernumerary teats, and 69 control cases with two teats). We detected 63 significant SNPs at the chromosome-wise threshold. Additionally, one candidate region (chr1: 170.723-170.734 Mb) was identified by haplotype-based association tests, with one SNP (rs413490006) surrounding functional genes BBX and CD47 on chromosome 1 being commonly identified as significant by the two mentioned analyses. Moreover, Gene Ontology enrichment for the significant SNPs identified by the GWAS analysis was functionally clustered into the categories of receptor activity and synaptic membrane. In addition, pathway mapping revealed four promising pathways (Wnt, oxytocin, MAPK and axon guidance) involved in the development of the supernumerary nipple phenotype. Our results provide novel and important insights into the genetic mechanisms underlying the phenotype of supernumerary nipples in mammals, including humans. These findings may be useful for future breeding and genetics in sheep and other livestock.

Horns are a cranial appendage found exclusively in Bovidae, and play important roles in accessing resources and mates. In sheep (Ovies aries), horns vary from polled to six-horned, and human have been selecting polled animals in farming and breeding. Here, we conducted a genome-wide association study on 24 two-horned versus 22 four-horned phenotypes in a native Chinese breed of Sishui Fur sheep. Together with linkage disequilibrium (LD) analyses and haplotype-based association tests, we identified a genomic region comprising 132.0–133.1 Mb on chromosome 2 that contained the top 10 SNPs (including 4 significant SNPs) and 5 most significant haplotypes associated with the polycerate phenotype. In humans and mice, this genomic region contains the HOXD gene cluster and adjacent functional genes EVX2 and KIAA1715, which have a close association with the formation of limbs and genital buds. Our results provide new insights into the genetic basis underlying variable numbers of horns and represent a new resource for use in sheep genetics and breeding.

No abstract