Genetic Basis of Dorper Sheep (Ovis aries) Revealed by Long-Read De Novo Genome Assembly

Genetic advancements and future directions in ruminant livestock breeding: from reference genomes to multiomics innovations

Ruminant livestock provide a rich source of products, such as meat, milk, and wool, and play a critical role in global food security and nutrition. Over the past few decades, genomic studies of ruminant livestock have provided valuable insights into their domestication and the genetic basis of economically important traits, facilitating the breeding of elite varieties. In this review, we summarize the main advancements for domestic ruminants in reference genome assemblies, population genomics, and the identification of functional genes or variants for phenotypic traits. These traits include meat and carcass quality, reproduction, milk production, feed efficiency, wool and cashmere yield, horn development, tail type, coat color, environmental adaptation, and disease resistance. Functional genomic research is entering a new era with the advancements of graphical pangenomics and telomere-to-telomere (T2T) gap-free genome assembly. These advancements promise to improve our understanding of domestication and the molecular mechanisms underlying economically important traits in ruminant livestock. Finally, we provide new perspectives and future directions for genomic research on ruminant genomes. We suggest how ever-increasing multiomics datasets will facilitate future studies and molecular breeding in livestock, including the potential to uncover novel genetic mechanisms underlying phenotypic traits, to enable more accurate genomic prediction models, and to accelerate genetic improvement programs.

Tissue-specific expression, functional analysis, and polymorphism of the KRT2 gene in sheep horn

Horn is a defensive weapon of sheep, consisting of a horny sheath and a bony core. The KRT2 gene is related to keratinization of the epidermis, so it is likely to be one of the contributor genes affecting horn type in sheep. In this study, we first analyzed the species-specific and tissue-specific expression of the KRT2 gene using transcriptome sequencing data. Then, by comparing the protein sequences of 20 species, we identified 28 specific amino acid sites in Artiodactyla animals, constructed a phylogenetic tree of the KRT2 gene, and predicted its three-dimensional protein structure. Finally, whole genome sequencing data was used and mined 4 functional SNP sites of KRT2 gene, and use KASP assay to verify the loci. In addition, we explored the relationship between the KRT2 gene and the evolution of Artiodactyla animals, and predicted the possible mechanism by which the KRT2 gene affects the horn type of sheep.

History and genetic diversity of African sheep: Contrasting phenotypic and genomic diversity

Domesticated sheep have adapted to contrasting and extreme environments and continue to play important roles in local community-based economies throughout Africa. Here we review the Neolithic migrations of thin-tailed sheep and the later introductions of fat-tailed sheep into eastern Africa. According to contemporary pictorial evidence, the latter occurred in Egypt not before the Ptolemaic period (305-25 BCE). We further describe the more recent history of sheep in Egypt, the Maghreb, west and central Africa, central-east Africa, and southern Africa. We also present a comprehensive molecular survey based on the analysis of 50 K SNP genotypes for 59 African breeds contributed by several laboratories. We propose that gene flow and import of fat-tailed sheep have partially overwritten the diversity profile created by the initial migration. We found a genetic contrast between sheep north and south of the Sahara and a west-east contrast of thin- and fat-tailed sheep. There is no close relationship between African and central and east Asian fat-tailed breeds, whereas we observe within Africa only a modest effect of tail types on breed relationships.

Significance of KLK7 expression, polymorphisms, and function in sheep horn growth

BACKGROUND

Sheep horns play a critical role in the survival and reproduction of sheep. Research on sheep horns not only aids in comprehending their biological roles but is also vital for developing hornless breeds. Although previous studies have suggested that KLK7 may be associated with keratin growth, there are few studies that have focused on the role of KLK7 in sheep horns. This study aims to elucidate the relationship between KLK7 and sheep horns by analyzing the expression, genetic polymorphisms, and potential functions of KLK7 in sheep horns.

RESULTS

This study utilized RNA sequencing (RNA-Seq) data to analyze the expression levels of the KLK7 gene across different species, as well as among different breeds, tissue types, and genders in sheep. Potential functional sites of KLK7 were explored using whole-genome sequencing (WGS) data. Allele specific expression sites in the KLK7 gene sequence were identified. Finally, the WGS data were linked with sheep horn length for association analysis, and significantly different single-nucleotide polymorphisms (SNPs) were detected and validated by Kompetitive Allele Specific PCR (KASP) genotyping.

CONCLUSION

Our results demonstrate that KLK7 was highly expressed in soft horn and skin tissues, and its expression was significantly higher in small-horned sheep than it was in large-horned sheep, suggesting that KLK7 may have an inhibitory effect on horn growth. By comparing the amino acid sequence of KLK7 with KLK7 sequences in other species, we discovered eight amino acids at specific positions in the KLK7 protein sequence that may have regulatory functions in determining the size of horns in ruminant animals. Thirteen SNPs with F-statistic value (Fst) > 0.15 were identified. By integrating RNA-Seq and WGS data, we discovered two SNPs (g.56695395 T > C and g.56695484A > C) with allele specific expression between the large- and small-horned sheep. The two SNPs were validated and were found to have significantly different (P < 0.05) effects on horn length. Our findings suggest a strong association between KLK7 and sheep horn length, indicating the potential role of KLK7 in inhibiting horn growth and providing novel insights into the functionality of KLK7.

Identification of ALOX12B Gene Expression, Evolution, and Potential Functional Sites in Horn Development of Sheep

The growth and development of horns are primarily controlled by the skin. The ALOX12B gene is crucial for epidermal barrier function and may have a significant impact on horn growth. The purpose of this study was to investigate the expression of ALOX12B across different sheep breeds and tissues by utilizing RNA sequencing. Additionally, potential functional sites were identified in conjunction with whole genome sequencing. Our findings revealed that ALOX12B was highly expressed in the scurred horn group as opposed to the normal horn group (SHE). ALOX12B expression was also notably high in the skin across several species. Eight loci that may influence horn size were indicated in this study. Through the alignment of the ALOX12B protein sequence from 16 species, 15 amino acid sites were identified specifically expressed in horned animals. In conclusion, this study established a connection between ALOX12B and horn size and identified a series of functional sites that may serve as molecular markers for reducing the presence of horns in Chinese sheep breeds.

Tissue-specific and functional loci analysis of CASP14 gene in the sheep horn

Under the current context of intensive farming, small-horned animals are more suitable for large-scale breeding. The CASP14 gene, closely associated with skin and keratin formation, may influence horn size due to its link with skin development. This study comprehensively analyzed the tissue-specific expression of CASP14 using RNA-Seq data, identified functional sites through whole-genome sequencing (WGS), and investigated allele-specific expression (ASE) validated by KASP assays. Results showed significantly higher CASP14 expression in the scurred group com-pared to the SHE group, with pronounced expression in the skin. Interbreed comparisons also revealed elevated CASP14 levels in the scurred group. Analysis of potential functional sites indicated structural similarities in the CASP14 protein among horned mammals. Five ASE events were discovered, and intersecting these with SNPs and high fixation index loci (Fst > 0.05) identified three potential functional sites: 7941628, 7941817, and 7941830. The SNP site 7944295 was selected for T-test analysis and further validated through KASP assays, corroborating the role of CASP14 in sheep horn phenotypes. Our findings suggest that CASP14 plays a significant role in horn development, offering insights into breeding strategies for small-horned animals.

Structural variant landscapes reveal convergent signatures of evolution in sheep and goats

BACKGROUND

Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution.

RESULTS

We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats.

CONCLUSIONS

Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.

Effects of nonsynonymous single nucleotide polymorphisms of the KIAA1217, SNTA1 and LTBP1 genes on the growth traits of Ujumqin sheep

Sheep body size can directly reflect the growth rates and fattening rates of sheep and is also an important index for measuring the growth performance of meat sheep. In this study, high-resolution resequencing data from four sheep breeds (Dorper sheep, Suffolk sheep, Ouessant sheep, and Shetland sheep) were analyzed. The nonsynonymous single nucleotide polymorphisms of three candidate genes (KIAA1217, SNTA1, and LTBP1) were also genotyped in 642 healthy Ujumqin sheep using MALDI-TOFMS and the genotyping results were associated with growth traits. The results showed that different genotypes of the KIAA1217 g.24429511T>C locus had significant effects on the chest circumferences of Ujumqin sheep. The SNTA1 g.62222626C>A locus had different effects on the chest depths, shoulder widths and rump widths of Ujumqin sheep. This study showed that these two sites can be used for marker-assisted selection, which will be beneficial for future precision molecular breeding.

Expression, Polymorphism, and Potential Functional Sites of the BMPR1A Gene in the Sheep Horn

Sheep horns are composed of bone and sheaths, and the BMPR1A gene is required for cartilage and osteogenic differentiation. Therefore, the BMPR1A gene may have a function related to the sheep horn, but its relationship with the sheep horn remains unclear. In this study, we first utilized RNA sequencing (RNA-seq) data to investigate the expression of the BMPR1A gene in different tissues and breeds of sheep. Second, whole-genome sequencing (WGS) data were used to explore the functional sites of the BMPR1A gene. Lastly, the allele-specific expression of the BMPR1A gene was explored. Our results indicate that BMPR1A gene expression is significantly higher in the normal horn groups than in the scurred groups. Importantly, this trend is consistent across several sheep breeds. Therefore, this finding suggests that the BMPR1A gene may be related to horn type. A total of 43 Single-Nucleotide Polymorphisms (SNPs) (F-statistics > 0.15) and 10 allele-specific expressions (ASEs) exhibited difference between the large and small horn populations. It is probable that these sites significantly impact the size of sheep horns. Compared to other polled species, we discovered ten amino acid sites that could influence horn presence. By combining RNA-seq and WGS functional loci results, we identified a functional site at position 40574836 on chromosome 25 that is both an SNP and exhibits allele-specific expression. In conclusion, we demonstrated that the BMPR1A gene is associated with horn type and identified some important functional sites which can be used as molecular markers in the breeding of sheep horns.

When Livestock Genomes Meet Third-Generation Sequencing Technology: From Opportunities to Applications

Third-generation sequencing technology has found widespread application in the genomic, transcriptomic, and epigenetic research of both human and livestock genetics. This technology offers significant advantages in the sequencing of complex genomic regions, the identification of intricate structural variations, and the production of high-quality genomes. Its attributes, including long sequencing reads, obviation of PCR amplification, and direct determination of DNA/RNA, contribute to its efficacy. This review presents a comprehensive overview of third-generation sequencing technologies, exemplified by single-molecule real-time sequencing (SMRT) and Oxford Nanopore Technology (ONT). Emphasizing the research advancements in livestock genomics, the review delves into genome assembly, structural variation detection, transcriptome sequencing, and epigenetic investigations enabled by third-generation sequencing. A comprehensive analysis is conducted on the application and potential challenges of third-generation sequencing technology for genome detection in livestock. Beyond providing valuable insights into genome structure analysis and the identification of rare genes in livestock, the review ventures into an exploration of the genetic mechanisms underpinning exemplary traits. This review not only contributes to our understanding of the genomic landscape in livestock but also provides fresh perspectives for the advancement of research in this domain.

Pan-Omics in Sheep: Unveiling Genetic Landscapes

Multi-omics-integrated analysis, known as panomics, represents an advanced methodology that harnesses various high-throughput technologies encompassing genomics, epigenomics, transcriptomics, proteomics, and metabolomics. Sheep, playing a pivotal role in agricultural sectors due to their substantial economic importance, have witnessed remarkable advancements in genetic breeding through the amalgamation of multiomics analyses, particularly with the evolution of high-throughput technologies. This integrative approach has established a robust theoretical foundation, enabling a deeper understanding of sheep genetics and fostering improvements in breeding strategies. The comprehensive insights obtained through this approach shed light on diverse facets of sheep development, including growth, reproduction, disease resistance, and the quality of livestock products. This review primarily focuses on the application of principal omics analysis technologies in sheep, emphasizing correlation studies between multiomics data and specific traits such as meat quality, wool characteristics, and reproductive features. Additionally, this paper anticipates forthcoming trends and potential developments in this field.

Dynamic transcriptome analysis of Maiwa yak corpus luteum during the estrous cycle

Maiwa yak is a special breed of animal living on the Qinghai-Tibet Plateau, which has great economic value, but its fertility rate is low. The corpus luteum (CL) is a temporary tissue that plays a crucial role in maintaining the physiological cycle. However, little is known about the transcriptome profile in Maiwa yak CL. In the present study, the transcriptome of Maiwa yak CL at early (EYCL), middle (MYCL) and late-stages (LYCL) was studied employing high-throughput sequencing. A total of 25,922 transcripts were identified, including 22,277 known as well as 3,645 novel ones. Furthermore, 690 and 212 differentially expressed (DE) mRNAs were detected in the EYCL vs. MYCL and MYCL vs. LYCL groups, respectively. KEGG pathway enrichment analysis of DEGs illustrated that the most enriched pathway was PI3K-Akt pathway. Furthermore, twenty-six DEGs were totally found to be associated with different biological processes of CL development. One of these genes, PGRMC1, displayed a dynamical expression trend during the lifespan of yak CL. The knockdown of PGRMC1 in luteinized yak granulosa cells resulted in defective steroidogenesis. In conclusion, this study analyzed the transcriptome profiles in yak CL of different stages, and provided a novel database for analyzing the gene network in yak CL.HIGHLIGHTSThe manuscript analyzed the transcriptome profiles in yak CL during the estrous cycle.Twenty-six DEGs were found to be associated with the development or function of CL.One of the DEGs, PGRMC1, was found to be responsible for steroidogenesis in luteinized yak granulosa cells.

Recent advances in the genomic resources for sheep

Sheep (Ovis aries) provide a vital source of protein and fibre to human populations. In coming decades, as the pressures associated with rapidly changing climates increase, breeding sheep sustainably as well as producing enough protein to feed a growing human population will pose a considerable challenge for sheep production across the globe. High quality reference genomes and other genomic resources can help to meet these challenges by: (1) informing breeding programmes by adding a priori information about the genome, (2) providing tools such as pangenomes for characterising and conserving global genetic diversity, and (3) improving our understanding of fundamental biology using the power of genomic information to link cell, tissue and whole animal scale knowledge. In this review we describe recent advances in the genomic resources available for sheep, discuss how these might help to meet future challenges for sheep production, and provide some insight into what the future might hold.

Chromosome-level assembly of the Rangifer tarandus genome and validation of cervid and bovid evolution insights

BACKGROUND

Genome assembly into chromosomes facilitates several analyses including cytogenetics, genomics and phylogenetics. Despite rapid development in bioinformatics, however, assembly beyond scaffolds remains challenging, especially in species without closely related well-assembled and available reference genomes. So far, four draft genomes of Rangifer tarandus (caribou or reindeer, a circumpolar distributed cervid species) have been published, but none with chromosome-level assembly. This emblematic northern species is of high interest in ecological studies and conservation since most populations are declining.

RESULTS

We have designed specific probes based on Oligopaint FISH technology to upgrade the latest published reindeer and caribou chromosome-level genomes. Using this oligonucleotide-based method, we found six mis-assembled scaffolds and physically mapped 68 of the largest scaffolds representing 78% of the most recent R. tarandus genome assembly. Combining physical mapping and comparative genomics, it was possible to document chromosomal evolution among Cervidae and closely related bovids.

CONCLUSIONS

Our results provide validation for the current chromosome-level genome assembly as well as resources to use chromosome banding in studies of Rangifer tarandus.