Genomic analyses of wild argali, domestic sheep, and their hybrids provide insights into chromosome evolution, phenotypic variation, and germplasm innovation.

Understanding the genetic mechanisms of phenotypic variation in hybrids between domestic animals and their wild relatives may aid germplasm innovation. Here, we report the high-quality genome assemblies of a male Pamir argali (O ammon polii, 2n = 56), a female Tibetan sheep (O aries, 2n = 54), and a male hybrid of Pamir argali and domestic sheep, and the high-throughput sequencing of 425 ovine animals, including the hybrids of argali and domestic sheep. We detected genomic synteny between Chromosome 2 of sheep and two acrocentric chromosomes of argali. We revealed consistent satellite repeats around the chromosome breakpoints, which could have resulted in chromosome fusion. We observed many more hybrids with karyotype 2n = 54 than with 2n = 55, which could be explained by the selfish centromeres, the possible decreased rate of normal/balanced sperm, and the increased incidence of early pregnancy loss in the aneuploid ewes or rams. We identified genes and variants associated with important morphological and production traits (e.g., body weight, cannon circumference, hip height, and tail length) that show significant variations. We revealed a strong selective signature at the mutation (c.334C > A, p.G112W) in TBXT and confirmed its association with tail length among sheep populations of wide geographic and genetic origins. We produced an intercross population of 110 F2 offspring with varied number of vertebrae and validated the causal mutation by whole-genome association analysis. We verified its function using CRISPR-Cas9 genome editing. Our results provide insights into chromosomal speciation and phenotypic evolution and a foundation of genetic variants for the breeding of sheep and other animals.

Historical Introgression from Wild Relatives Enhanced Climatic Adaptation and Resistance to Pneumonia in Sheep.

How animals, particularly livestock, adapt to various climates and environments over short evolutionary time is of fundamental biological interest. Further, understanding the genetic mechanisms of adaptation in indigenous livestock populations is important for designing appropriate breeding programs to cope with the impacts of changing climate. Here, we conducted a comprehensive genomic analysis of diversity, interspecies introgression, and climate-mediated selective signatures in a global sample of sheep and their wild relatives. By examining 600K and 50K genome-wide single nucleotide polymorphism data from 3,447 samples representing 111 domestic sheep populations and 403 samples from all their seven wild relatives (argali, Asiatic mouflon, European mouflon, urial, snow sheep, bighorn, and thinhorn sheep), coupled with 88 whole-genome sequences, we detected clear signals of common introgression from wild relatives into sympatric domestic populations, thereby increasing their genomic diversities. The introgressions provided beneficial genetic variants in native populations, which were significantly associated with local climatic adaptation. We observed common introgression signals of alleles in olfactory-related genes (e.g., ADCY3 and TRPV1) and the PADI gene family including in particular PADI2, which is associated with antibacterial innate immunity. Further analyses of whole-genome sequences showed that the introgressed alleles in a specific region of PADI2 (chr2: 248,302,667-248,306,614) correlate with resistance to pneumonia. We conclude that wild introgression enhanced climatic adaptation and resistance to pneumonia in sheep. This has enabled them to adapt to varying climatic and environmental conditions after domestication.

Comparative mRNA and miRNA expression in European mouflon (Ovis musimon) and sheep (Ovis aries) provides novel insights into the genetic mechanisms for female reproductive success.

Prolific breeds of domestic sheep (Ovis aries) are important genetic resources due to their reproductive performance, which is characterized by multiple lambs per birth and out-of-season breeding. However, the lack of a comprehensive understanding of the genetic mechanisms underlying the important reproductive traits, particularly from the evolutionary genomics perspective, has impeded the efficient advancement of sheep breeding. Here, for the first time, by performing RNA-sequencing we built a de novo transcriptome assembly of ovarian and endometrial tissues in European mouflon (Ovis musimon) and performed an mRNA-miRNA integrated expression profiling analysis of the wild species and a highly prolific domestic sheep breed, the Finnsheep. We identified several novel genes with differentially expressed mRNAs (e.g., EREG, INHBA, SPP1, AMH, TDRD5, and ZP2) between the wild and domestic sheep, which are functionally involved in oocyte and follicle development and fertilization, and are significantly (adjusted P-value < 0.05) enriched in the Gene Ontology (GO) terms of various reproductive process, including the regulation of fertilization, oogenesis, ovarian follicle development, and sperm-egg recognition. Additionally, we characterized 58 differentially expressed miRNAs and 210 associated target genes that are essential for the regulation of female reproduction cycles through specific regulatory networks [e.g., (miR-136, miR-374a, miR-9-5p)-(EREG, INHBA)]. Furthermore, our integrated mRNA and miRNA expression profiling analysis elucidated novel direct and indirect miRNA/mRNA causal regulatory relationships related to the reproductive traits of the Ovis species. This study provides in-depth insights into the genomic evolution underlying the reproductive traits of the Ovis species and valuable resources for ovine genomics.

[Analysis of lorentzian line shape function broadened by non-sinusoidal wavelength modulation].

In the present work, the Fourier analysis of Lorentzian line shape broadened by non-sinusoidal wavelength modulation was investigated, in which the third order and above harmonic items were ignored. The analytical expression of n-order Fourier coefficient was brought out, where a variable K named harmonic distortion to characterize the ratio of the second harmonic to the first harmonic was introduced. Numerical simulations based on the cases of K > 0.01 and K < 0.01 were carried out, and the result shows: non-sinusoidal modulation has little effect compared with the sinusoidal modulation when K value is less than 0.01, however, if K value is about 0.1 or higher, the center of the Fourier amplitude curve would deviate from the origin of coordinates. With the increase in the harmonic distortion, the deviation of the curve grows, and high order harmonics are more sensitive to the non-sinusoidal modulation compared with the low order harmonics. In addition, when harmonic distortion cannot be ignored, for example K > 0.01, the effect of different depths of modulation on the odd and even order harmonic amplitude curve is significant. And the numerical simulation shows there exists an optimum value of modulation depth which could minimize the impact of the harmonic distortion, and both large K value and small K value would cause a great error. The conclusion of this work could be applied in error analysis of wavelength modulation spectroscopy system And the results are helpful to deepening understanding of WMS and would be the important reference for some kind of frequency stabilization technology in laser instrument.

Characterization of Gut Microbiota Compositions along the Intestinal Tract in CD163/pAPN Double Knockout Piglets and Their Potential Roles in Iron Absorption.

The gut microbiota is known to play a role in regulating host metabolism, yet the mechanisms underlying this regulation are not well elucidated. Our study aimed to characterize the differences in gut microbiota compositions and their roles in iron absorption between wild-type (WT) and CD163/pAPN double-gene-knockout (DKO) weaned piglets. A total of 58 samples along the entire digestive tract were analyzed for microbial community using 16S rRNA gene sequencing. The colonic microbiota and their metabolites were determined by metagenomic sequencing and untargeted liquid chromatography-mass spectrometry (LC-MS), respectively. Our results showed that no alterations in microbial community structure and composition were observed between DKO and WT weaned piglets, with the exception of colonic microbiota. Interestingly, the DKO piglets had selectively increased the relative abundance of the Leeia genus belonging to the Neisseriaceae family and decreased the Ruminococcaceae_UCG_014 genus abundance. Functional capacity analysis showed that organic acid metabolism was enriched in the colon in DKO piglets. In addition, the DKO piglets showed increased iron levels in important tissues compared with WT piglets without any pathological changes. Pearson's correlation coefficient indicated that the specific bacteria such as Leeia and Ruminococcaceae_UCG_014 genus played a key role in host iron absorption. Moreover, the iron levels had significantly (P < 0.05) positive correlation with microbial metabolites, particularly carboxylic acids and their derivatives, which might increase iron absorption by preventing iron precipitation. Overall, this study reveals an interaction between colonic microbiota and host metabolism and has potential significance for alleviating piglet iron deficiency. IMPORTANCE Iron deficiency is a major risk factor for iron deficiency anemia, which is among the most common nutritional disorders in piglets. However, it remains unclear how the gut microbiota interacts with host iron absorption. The current report provides the first insight into iron absorption-microbiome connection in CD163/pAPN double knockout piglets. The present results showed that carboxylic acids and their derivatives contributed to the absorption of nonheme iron by preventing ferric iron precipitation.

Changes in the Mucosa-Associated Microbiome and Transcriptome across Gut Segments Are Associated with Obesity in a Metabolic Syndrome Porcine Model.

Several studies have suggested a role for gut mucosa-associated microbiota in the development of obesity, but the mechanisms involved are poorly defined. Here, the impact of the gut mucosa-associated microbiota on obesity and related metabolic disorders was evaluated in a metabolic syndrome (MetS) porcine model. Body composition was determined among male Wuzhishan minipigs consuming a high-energy diet (HED) and compared to that of those consuming a normal diet (ND), and gut segments (duodenum, jejunum, ileum, cecum, colon, and rectum) were sampled for paired analysis of mucosa-associated microbiota and transcriptome signatures with 16S rRNA gene and RNA sequencing, respectively. Our data indicated that long-term HED feeding significantly increased body weight and visceral fat deposition and aggravated metabolic disorders. Specially, HED feeding induced mucosa-associated microbiota dysbiosis and selectively increased the abundance of the families Enterobacteriaceae, Moraxellaceae, and Lachnospiraceae in the upper intestine. The association analysis indicated that specific bacteria play key roles in adiposity, e.g., Lactobacillus johnsonii in the duodenum, Actinobacillus indolicus in the jejunum, Acinetobacter johnsonii in the ileum, Clostridium butyricum in the cecum, Haemophilus parasuis in the colon, and bacterium NLAEzlP808, Halomonas taeheungii, and Shewanella sp. JNUH029 in the rectum. Transcriptome data further revealed intestinal lipid metabolism and immune dysfunction in the MetS individuals, which may be associated with obesity and related metabolic disorders. Our results indicated that gut mucosa-associated microbiota dysbiosis has the potential to exacerbate obesity, partially through modulating systemic inflammatory responses. IMPORTANCE Obesity is a major risk factor for metabolic syndrome, which is the most common cause of death worldwide, especially in developed countries. The link between obesity and gut mucosa-associated microbiota is unclear due to challenges associated with the collection of intestinal samples from humans. The current report provides the first insight into obesity-microbiome-gut immunity connections in a metabolic syndrome (MetS) porcine model. The present results show that dysbiosis of mucosal microbiota along the entire digestive tract play a critical role in the proinflammatory response in the host-microbial metabolism axis, resulting in obesity and related metabolic disorders in the MetS model.

Multi-Omic Analysis in a Metabolic Syndrome Porcine Model Implicates Arachidonic Acid Metabolism Disorder as a Risk Factor for Atherosclerosis.

Background: The diet-induced gut microbiota dysbiosis has been suggested as a major risk factor for atherothrombosis, however, the detailed mechanism linking these conditions is yet to be fully understood. Methods: We established a long-term excessive-energy diet-induced metabolic syndrome (MetS) inbred Wuzhishan minipig model, which is characterized by its genetic stability, small size, and human-like physiology. The metabolic parameters, atherosclerotic lesions, gut microbiome, and host transcriptome were analyzed. Metabolomics profiling revealed a linkage between gut microbiota and atherothrombosis. Results: We showed that white atheromatous plaque was clearly visible on abdominal aorta in the MetS model. Furthermore, using metagenome and metatranscriptome sequencing, we discovered that the long-term excessive energy intake altered the local intestinal microbiota composition and transcriptional profile, which was most dramatically illustrated by the reduced abundance of SCFAs-producing bacteria including Bacteroides, Lachnospiraceae, and Ruminococcaceae in the MetS model. Liver and abdominal aorta transcriptomes in the MetS model indicate that the diet-induced gut microbiota dysbiosis activated host chronic inflammatory responses and significantly upregulated the expression of genes related to arachidonic acid-dependent signaling pathways. Notably, metabolomics profiling further revealed an intimate linkage between arachidonic acid metabolism and atherothrombosis in the host-gut microbial metabolism axis. Conclusions: These findings provide new insights into the relationship between atherothrombosis and regulation of gut microbiota via host metabolomes and will be of potential value for the treatment of cardiovascular diseases in MetS.

Whole-Genome Selective Scans Detect Genes Associated With Important Phenotypic Traits in Sheep (Ovis aries).

Sheep (Ovis aries) is one of the important livestock with diverse phenotypic traits. However, little is known about the molecular mechanism of diverse phenotypic traits in domestic sheep. Using the genome-wide high-density SNP data (600K) in 253 samples from 13 populations, we conducted the tests of selective sweeps (i.e., pairwise F ST and XP-CLR) associated with several important phenotypic traits (e.g., tail types, horn morphology, prolificacy, coat pigmentation, ear size, milk production, meat production, body size and wool fineness). We identified strong selective signatures in previously reported (e.g., T, RXFP2, BMPR1B, TYRP1, MSRB3, TF, CEBPA, GPR21 and HOXC8) and novel genes associated with the traits, such as CERS6, BTG1, RYR3, SLC6A4, NNAT and OGT for fat deposition in the tails, FOXO4 for fertility, PTCH1 and EMX2 for ear size, and RMI1 and SCD5 for body size. Further gene annotation analysis showed that these genes were identified to be the most probable genes accounting for the diverse phenotypic traits. Our results provide novel insights into the genetic mechanisms underlying the traits and also new genetic markers for genetic improvement in sheep and other livestock.

Genome-Wide Detection of Copy Number Variations and Their Association With Distinct Phenotypes in the World's Sheep.

Copy number variations (CNVs) are a major source of structural variation in mammalian genomes. Here, we characterized the genome-wide CNV in 2059 sheep from 67 populations all over the world using the Ovine Infinium HD (600K) SNP BeadChip. We tested their associations with distinct phenotypic traits by conducting multiple independent genome-wide tests. In total, we detected 7547 unique CNVs and 18,152 CNV events in 1217 non-redundant CNV regions (CNVRs), covering 245 Mb (∼10%) of the whole sheep genome. We identified seven CNVRs with frequencies correlating to geographical origins and 107 CNVRs overlapping 53 known quantitative trait loci (QTLs). Gene ontology and pathway enrichment analyses of CNV-overlapping genes revealed their common involvement in energy metabolism, endocrine regulation, nervous system development, cell proliferation, immune, and reproduction. For the phenotypic traits, we detected significantly associated (adjusted P < 0.05) CNVRs harboring functional candidate genes, such as SBNO2 for polycerate; PPP1R11 and GABBR1 for tail weight; AKT1 for supernumerary nipple; CSRP1, WNT7B, HMX1, and FGFR3 for ear size; and NOS3 and FILIP1 in Wadi sheep; SNRPD3, KHDRBS2, and SDCCAG3 in Hu sheep; NOS3, BMP1, and SLC19A1 in Icelandic; CDK2 in Finnsheep; MICA in Romanov; and REEP4 in Texel sheep for litter size. These CNVs and associated genes are important markers for molecular breeding of sheep and other livestock species.

Whole-genome resequencing of wild and domestic sheep identifies genes associated with morphological and agronomic traits.

Understanding the genetic changes underlying phenotypic variation in sheep (Ovis aries) may facilitate our efforts towards further improvement. Here, we report the deep resequencing of 248 sheep including the wild ancestor (O. orientalis), landraces, and improved breeds. We explored the sheep variome and selection signatures. We detected genomic regions harboring genes associated with distinct morphological and agronomic traits, which may be past and potential future targets of domestication, breeding, and selection. Furthermore, we found non-synonymous mutations in a set of plausible candidate genes and significant differences in their allele frequency distributions across breeds. We identified PDGFD as a likely causal gene for fat deposition in the tails of sheep through transcriptome, RT-PCR, qPCR, and Western blot analyses. Our results provide insights into the demographic history of sheep and a valuable genomic resource for future genetic studies and improved genome-assisted breeding of sheep and other domestic animals.

Genome-Wide Association Analyses Highlight the Potential for Different Genetic Mechanisms for Litter Size Among Sheep Breeds.

Reproduction is an important trait in sheep breeding as well as in other livestock. However, despite its importance the genetic mechanisms of litter size in domestic sheep (Ovis aries) are still poorly understood. To explore genetic mechanisms underlying the variation in litter size, we conducted multiple independent genome-wide association studies in five sheep breeds of high prolificacy (Wadi, Hu, Icelandic, Finnsheep, and Romanov) and one low prolificacy (Texel) using the Ovine Infinium HD BeadChip, respectively. We identified different sets of candidate genes associated with litter size in different breeds: BMPR1B, FBN1, and MMP2 in Wadi; GRIA2, SMAD1, and CTNNB1 in Hu; NCOA1 in Icelandic; INHBB, NF1, FLT1, PTGS2, and PLCB3 in Finnsheep; ESR2 in Romanov and ESR1, GHR, ETS1, MMP15, FLI1, and SPP1 in Texel. Further annotation of genes and bioinformatics analyses revealed that different biological pathways could be involved in the variation in litter size of females: hormone secretion (FSH and LH) in Wadi and Hu, placenta and embryonic lethality in Icelandic, folliculogenesis and LH signaling in Finnsheep, ovulation and preovulatory follicle maturation in Romanov, and estrogen and follicular growth in Texel. Taken together, our results provide new insights into the genetic mechanisms underlying the prolificacy trait in sheep and other mammals, suggesting targets for selection where the aim is to increase prolificacy in breeding projects.

Polymorphisms of FST gene and their association with wool quality traits in Chinese Merino sheep.

Follistatin (FST) is involved in hair follicle morphogenesis. However, its effects on hair traits are not clear. This study was designed to investigate the effects of FST gene single nucleotide polymorphisms (SNP) on wool quality traits in Chinese Merino sheep (Junken Type). We performed gene expression analysis, SNP detection, and association analysis of FST gene with sheep wool quality traits. The real-time RT-PCR analysis showed that FST gene was differentially expressed in adult skin between Chinese Merino sheep (Junken Type) and Suffolk sheep. Immunostaining showed that FST was localized in inner root sheath (IRS) and matrix of hair follicle (HF) in both SF and Suffolk sheep. Sequencing analysis identified a total of seven SNPs (termed SNPs 1-7) in the FST gene in Chinese Merino sheep (Junken Type). Association analysis showed that SNP2 (Chr 16. 25,633,662 G>A) was significantly associated with average wool fiber diameter, wool fineness SD, and wool crimp (P < 0.05). SNP4 (Chr 16. 25,633,569 C>T) was significantly associated with wool fineness SD and CV of fiber diameter (P < 0.05). Similarly, the haplotypes derived from these seven identified SNPs were also significantly associated with average wool fiber diameter, wool fineness SD, CV of fiber diameter, and wool crimp (P < 0.05). Our results suggest that FST influences wool quality traits and its SNPs 2 and 4 might be useful markers for marker-assisted selection and sheep breeding.