Data of whole-genome sequencing of Karakul, Zel, and Kermani sheep breeds.

OBJECTIVE: The data provided herein represent the whole-genome sequencing data associated with three sheep breeds of Iranian native breeds. Sheep are the first domesticated animals that, during the long path of the evolution process, have created gene variants with desirable phenotypic effects, so they can be suitable models for biomedical research. In addition, sheep have a vital role in providing protein to a notable part of the human population around the world. DATA DESCRIPTION: Ten blood samples were taken from three Iranian native sheep breeds, the Zel, Karakul, and Kermani kinds. Blood samples genomes were extracted using the salting-out technique. The Illumina NovaSeq 6000 platform was used to carry out sequencing of the whole genome in a laboratory in China. All sequence information is available through the NCBI database in the sequence read archive (SRA) format under the accession number PRJNA904537. The dataset presented here can provide a useful resource for genome analysis of livestock breeds adapted to hot and dry regions.

Whole-genome resequencing reveals the genomic diversity and signatures of selection in Romanov sheep.

Romanov sheep are adapted to the extremely cold and harsh environment and display a distinctive grey color. Herein, we analyzed the population structure, genetic diversity, and selection signatures of Romanov sheep based on whole-genome sequencing data of 17 Romanov sheep, 114 individuals from other 10 European breeds. The results of PCA, ADMIXTURE, and NJ-tree showed that the Romanov sheep was closely related to other northern European breeds. A relative high level of genetic diversity, low inbreeding coefficient, and large effective population size was observed in Romanov sheep when compared with other European breeds. We then screened the genomic selection signatures of Romanov sheep using FST, XP-XLP, and XP-EHH methods. The most significant region under selection (CHR14:14.2 to 14.3 Mb) harbored a haplotype that contained MC1R gene. Furthermore, this haplotype was also found in other grey-body breeds including Gotland sheep, Grey Tronder Sheep, and German grey heath sheep, suggesting that it was associated with the unique coat color of these breeds. We also found one region (CHR10:40.8Mb- 41.0Mb) harboring PCDH9 gene which was potentially associated with cold environmental adaptation. In summary, this study identified candidate genes that were associated with the unique grey color and environmental adaptation in Romanov sheep, which provided a basis for understanding the genetic background and utilization of this breed.

Romanov sheep is one of the most famous sheep breeds in the word, characterized by adaptability to harsh environment, high fertility, and unique coat color. Understanding its genetic architecture and signatures is of great value for its conservation and utilization. In this study, we analyzed whole-genome sequences of Romanov sheep as compared with 11 other European breeds, to explore for the population structure, genetic diversity, and selection signatures. We discovered a series of candidate genes that likely play a role in the grey coat color and cold adaptation of the Romanov sheep. In particular, we identified MC1R as a strong candidate gene that determines the grey coat color.

A Chinese indicine pangenome reveals a wealth of novel structural variants introgressed from other Bos species.

Chinese indicine cattle harbor a much higher genetic diversity compared with other domestic cattle, but their genome architecture remains uninvestigated. Using PacBio HiFi sequencing data from 10 Chinese indicine cattle across southern China, we assembled 20 high-quality partially phased genomes and integrated them into a multiassembly graph containing 148.5 Mb (5.6%) of novel sequence. We identified 156,009 high-confidence nonredundant structural variants (SVs) and 206 SV hotspots spanning ∼195 Mb of gene-rich sequence. We detected 34,249 archaic introgressed fragments in Chinese indicine cattle covering 1.93 Gb (73.3%) of the genome. We inferred an average of 3.8%, 3.2%, 1.4%, and 0.5% of introgressed sequence originating, respectively, from banteng-like, kouprey-like, gayal-like, and gaur-like Bos species, as well as 0.6% of unknown origin. Introgression from multiple donors might have contributed to the genetic diversity of Chinese indicine cattle. Altogether, this study highlights the contribution of interspecies introgression to the genomic architecture of an important livestock population and shows how exotic genomic elements can contribute to the genetic variation available for selection.

A sheep pangenome reveals the spectrum of structural variations and their effects on tail phenotypes.

Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here we generated high-quality genome assemblies for 15 individuals from genetically diverse sheep breeds using Pacific Biosciences (PacBio) high-fidelity sequencing, discovering 130.3 Mb nonreference sequences, from which 588 genes were annotated. A total of 149,158 biallelic insertions/deletions, 6531 divergent alleles, and 14,707 multiallelic variations with precise breakpoints were discovered. The SV spectrum is characterized by an excess of derived insertions compared to deletions (94,422 vs. 33,571), suggesting recent active LINE expansions in sheep. Nearly half of the SVs display low to moderate linkage disequilibrium with surrounding single-nucleotide polymorphisms (SNPs) and most SVs cannot be tagged by SNP probes from the widely used ovine 50K SNP chip. We identified 865 population-stratified SVs including 122 SVs possibly derived in the domestication process among 690 individuals from sheep breeds worldwide. A novel 168-bp insertion in the 5' untranslated region (5' UTR) of HOXB13 is found at high frequency in long-tailed sheep. Further genome-wide association study and gene expression analyses suggest that this mutation is causative for the long-tail trait. In summary, we have developed a panel of high-quality de novo assemblies and present a catalog of structural variations in sheep. Our data capture abundant candidate functional variations that were previously unexplored and provide a fundamental resource for understanding trait biology in sheep.

RGD v2.0: a major update of the ruminant functional and evolutionary genomics database.

Ruminant Genome Database (RGD; http://animal.nwsuaf.edu.cn/RGD) provides visualization and analysis tools for ruminant comparative genomics and functional annotations. As more high-quality ruminant genome assemblies have become available, we have redesigned the user interface, integrated and expanded multi-omics data, and developed novel features to improve the database. The new version, RGD v2.0, houses 78 ruminant genomes; 110-species synteny alignments for major livestock (including cattle, sheep, goat) and wild ungulates; 21 012 orthologous gene clusters with Gene Ontology and pathway annotation; ∼8 600 000 conserved elements; and ∼1 000 000 cis-regulatory elements by utilizing 1053 epigenomic data sets. The transcriptome data in RGD v2.0 has nearly doubled, currently with 1936 RNA-seq data sets, and 155 174 phenotypic data sets have been newly added. New and updated features include: (i) The UCSC Genome Browser, BLAT, BLAST and Table Browser tools were updated for six available ruminant livestock species. (ii) The LiftOver tool was newly introduced into our browser to allow coordinate conversion between different ruminant assemblies. And (iii) tissue specificity index, tau, was calculated to facilitate batch screening of specifically expressed genes. The enhanced genome annotations and improved functionality in RGD v2.0 will be useful for study of genome evolution, environmental adaption, livestock breeding and biomedicine.

Genome-wide scan for selective footprints and genes related to cold tolerance in Chantecler chickens.

The Chantecler chicken, a unique Canadian indigenous breed, is well adapted to extremely cold environments. However, its genetic characteristics have not been well studied. Here, we analyzed the whole genomes of 10 Chantecler chickens and 121 worldwide chickens, which indicated that Chantecler chickens were derived from commercial chickens and exhibit a high level of inbreeding. Based on a genome-wide scan, we identified two vital candidate regions containing ME3 and ZNF536, which are related to fat metabolism and nervous system in cold adaptation, respectively. We also found that the genetic mechanism of cold adaptation in Chantecler chickens differed from that of chickens from other cold regions, such as northern China. Our study indicated that specialized commercial chickens in the early 20th century contained sufficient genetic diversity to adapt to extreme cold environments over a very short time. These findings enrich our understanding of the adaptive potential of commercial species.

Genomic Analyses of Unveil Helmeted Guinea Fowl (Numida meleagris) Domestication in West Africa.

Domestication of the helmeted guinea fowl (HGF; Numida meleagris) in Africa remains elusive. Here we report a high-quality de novo genome assembly for domestic HGF generated by long- and short-reads sequencing together with optical and chromatin interaction mapping. Using this assembly as the reference, we performed population genomic analyses for newly sequenced whole-genomes for 129 birds from Africa, Asia, and Europe, including domestic animals (n = 89), wild progenitors (n = 34), and their closely related wild species (n = 6). Our results reveal domestication of HGF in West Africa around 1,300-5,500 years ago. Scanning for selective signals characterized the functional genes in behavior and locomotion changes involved in domestication of HGF. The pleiotropy and linkage in genes affecting plumage color and fertility were revealed in the recent breeding of Italian domestic HGF. In addition to presenting a missing piece to the jigsaw puzzle of domestication in poultry, our study provides valuable genetic resources for researchers and breeders to improve production in this species.

A Hu sheep genome with the first ovine Y chromosome reveal introgression history after sheep domestication.

The Y chromosome plays key roles in male fertility and reflects the evolutionary history of paternal lineages. Here, we present a de novo genome assembly of the Hu sheep with the first draft assembly of ovine Y chromosome (oMSY), using nanopore sequencing and Hi-C technologies. The oMSY that we generated spans 10.6 Mb from which 775 Y-SNPs were identified by applying a large panel of whole genome sequences from worldwide sheep and wild Iranian mouflons. Three major paternal lineages (HY1a, HY1b and HY2) were defined across domestic sheep, of which HY2 was newly detected. Surprisingly, HY2 forms a monophyletic clade with the Iranian mouflons and is highly divergent from both HY1a and HY1b. Demographic analysis of Y chromosomes, mitochondrial and nuclear genomes confirmed that HY2 and the maternal counterpart of lineage C represented a distinct wild mouflon population in Iran that diverge from the direct ancestor of domestic sheep, the wild mouflons in Southeastern Anatolia. Our results suggest that wild Iranian mouflons had introgressed into domestic sheep and thereby introduced this Iranian mouflon specific lineage carrying HY2 to both East Asian and Africa sheep populations.

Comparative Transcriptome Profiling Analysis Uncovers Novel Heterosis-Related Candidate Genes Associated with Muscular Endurance in Mules.

Heterosis has been widely exploited in animal and plant breeding programs to enhance the productive traits of hybrid progeny from two breeds or species. However, its underlying genetic mechanisms remain enigmatic. Transcriptome profiling analysis can be used as a method for exploring the mechanism of heterosis. Here, we performed genome-wide gene expression and alternative splicing (AS) analyses in different tissues (muscle, brain, and skin) from crosses between donkeys and horses. Our results indicated that 86.1% of the differentially expressed genes (DEGs) and 87.2% of the differential alternative splicing (DAS) genes showed over-dominance and dominance in muscle. Further analysis showed that the "muscle contraction" pathway was significantly enriched for both the DEGs and DAS genes in mule muscle tissue. Taken together, these DEGs and DAS genes could provide an index for future studies of the genetic and molecular mechanism of heterosis in the hybrids of donkey and horse.

Whole genomes and transcriptomes reveal adaptation and domestication of pistachio.

BACKGROUND: Pistachio (Pistacia vera), one of the most important commercial nut crops worldwide, is highly adaptable to abiotic stresses and is tolerant to drought and salt stresses. RESULTS: Here, we provide a draft de novo genome of pistachio as well as large-scale genome resequencing. Comparative genomic analyses reveal stress adaptation of pistachio is likely attributable to the expanded cytochrome P450 and chitinase gene families. Particularly, a comparative transcriptomic analysis shows that the jasmonic acid (JA) biosynthetic pathway plays an important role in salt tolerance in pistachio. Moreover, we resequence 93 cultivars and 14 wild P. vera genomes and 35 closely related wild Pistacia genomes, to provide insights into population structure, genetic diversity, and domestication. We find that frequent genetic admixture occurred among the different wild Pistacia species. Comparative population genomic analyses reveal that pistachio was domesticated about 8000 years ago and suggest that key genes for domestication related to tree and seed size experienced artificial selection. CONCLUSIONS: Our study provides insight into genetic underpinning of local adaptation and domestication of pistachio. The Pistacia genome sequences should facilitate future studies to understand the genetic basis of agronomically and environmentally related traits of desert crops.

Out of Southern East Asia of the Brown Rat Revealed by Large-Scale Genome Sequencing.

The geographic origin and migration of the brown rat (Rattus norvegicus) remain subjects of considerable debate. In this study, we sequenced whole genomes of 110 wild brown rats with a diverse world-wide representation. We reveal that brown rats migrated out of southern East Asia, rather than northern Asia as formerly suggested, into the Middle East and then to Europe and Africa, thousands of years ago. Comparison of genomes from different geographical populations reveals that many genes involved in the immune system experienced positive selection in the wild brown rat.