A global analysis of CNVs in Chinese indigenous fine-wool sheep populations using whole-genome resequencing

Prenatal diagnosis and genetic counseling of a paternally inherited chromosome 8q24.22q24.23 microdeletion in a Chinese family

BACKGROUND

Copy number variants (CNVs) are an important source of normal and pathogenic genome variations. Especially CNVs identified in prenatal cases need careful considerations and correct interpretation if those are harmless or harmful variants from the norm. The literature on 8q24.22q24.23 microdeletion is rare, which is a challenge for genetic counselling.

CASE PRESENTATION

We have performed prenatal diagnosis and genetic counseling of a paternally inherited 8q24.22q24.23 microdeletion. In this family, father with normal phenotype and fetus with abnormal phenotype have the same microdeletion.

CONCLUSION

Chromosomal microdeletions and microduplications are difficult to detect by conventional cytogenetics, combination of prenatal ultrasound, karyotype analysis, copy number variation sequencing (CNV-seq), whole-exome sequencing (WES) and genetic counseling is helpful for the prenatal diagnosis of chromosomal microdeletions/microduplications.

Leveraging Whole-Genome Resequencing to Uncover Genetic Diversity and Promote Conservation Strategies for Ruminants in Asia

Whole-genome resequencing (WGRS) is a critical branch of whole-genome sequencing (WGS), primarily targeting species with existing reference genomes. By aligning sequencing data to the reference genome, WGRS enables precise detection of genetic variations in individuals or populations. As a core technology in genomic research, WGS holds profound significance in ruminant studies. It not only reveals the intricate structure of ruminant genomes but also provides essential data for deciphering gene function, variation patterns, and evolutionary processes, thereby advancing the exploration of ruminant genetic mechanisms. However, WGS still faces several challenges, such as incomplete and inaccurate genome assembly, as well as the incomplete annotation of numerous unknown genes or gene functions. Although WGS can identify a vast number of genomic variations, the specific relationships between these variations and phenotypes often remain unclear, which limits its potential in functional studies and breeding applications. By performing WGRS on multiple samples, these assembly challenges can be effectively addressed, particularly in regions with high repeat content or complex structural variations. WGRS can accurately identify subtle variations among different individuals or populations and further elucidate their associations with specific traits, thereby overcoming the limitations of WGS and providing more precise genetic information for functional research and breeding applications. This review systematically summarizes the latest applications of WGRS in the analysis of ruminant genetic structures, genetic diversity, economic traits, and adaptive traits, while also discussing the challenges faced by this technology. It aims to provide a scientific foundation for the improvement and conservation of ruminant genetic resources.

Genome wide detection of CNV and their association with body size in Danzhou chickens

Copy number variation (CNV) is a crucial component of genetic diversity in the genome, serving as the foundation for the genetic architecture and phenotypic variability of complex traits. In this study, we examined CNVs in the Danzhou (DZ) chicken, an indigenous breed exclusive to Hainan Province, China. By employing whole-genome resequencing data from 200 DZ chickens, we conducted a comprehensive genome-wide analysis of CNVs using CNVpytor and performed CNV-based genome-wide association studies (GWAS) on 6 body size traits, including body slope length (BSL), keel length (KeL), tibial length (TiL), tibial circumference (TiC), chest width (ChW), and chest depth (ChD) utilizing linear mixed model methods considering a genomic relationship matrix. We identified a total of 144,265 autosomal CNVs among the 200 individuals, comprising 67,818 deletions and 76,447 duplications. After merging these variants together, we obtained 4,824 distinct copy number variant regions, which accounted for approximately 20% of the chicken autosomal genome. Furthermore, we discovered several significantly associated CNV segments with body size traits located proximal to genes such as IHH, WNT6, WNT10A, LPR4, FZD2, WNT7B, and GNAS that have been extensively implicated in skeletal development and growth processes. These findings enhance our understanding of CNVs in chickens and their potential impact on body size traits by revealing candidate genes involved in the regulation of these traits. This establishes a solid framework for future studies and may prove particularly beneficial for exploring genetic structural variation in chickens.

α-mannosidosis diagnosis in Brazilian patients with MPS-like symptoms

BACKGROUND

α-mannosidosis is an inborn error of metabolism caused by the deficiency of the lysosomal enzyme α-mannosidase, which is encoded by the MAN2B1 gene and inherited in an autosomal recessive manner. The impairment of affected individuals is multisystemic and very similar to the observed in some mucopolysaccharidosis (MPS) patients. The aim of this study was to search for α-mannosidosis cases in individuals with clinical suspicion of MPS without a confirmed diagnosis. Biochemical and molecular analysis were standardized by our group for this study. Two hundred and fifty samples from patients with clinical suspicion of MPS, but with inconclusive MPS biochemical and/or molecular analysis, were screened for α-mannosidase activity. Subsequently the MAN2B1 gene was sequenced in samples from 53 patients by the Sanger method.

RESULTS

The measurement of enzymatic activity detected fifty-three samples with abnormal results, suggesting α-mannosidosis. Molecular analysis confirmed three affected families, which presented the nonsense variant p.Ser899Ter. This variant generates a premature stop codon in exon 22, resulting in a truncated protein with no residual enzymatic activity.

CONCLUSION

In conclusion, this work brings data for the beginning of a genetic characterization of α-mannosidosis in the Brazilian population. It also shows that α-mannosidosis cases may be underdiagnosed due to the clinical similarity to MPS and the lack of information about this ultra-rare disease. Based on our data, we strongly recommend to all screening centers to consider α-mannosidosis testing together with screening for MPS as a tool for diagnosis to MPS-like phenotype individuals, since the phenotype similarity between these diseases poses a significant challenge for clinicians worldwide and often leads to the failure of the correct clinical diagnosis and treatment.

Genome-wide selection signatures address trait specific candidate genes in cattle indigenous to arid regions of India

The peculiarity of Indian cattle lies in milk quality, resistance to diseases and stressors as well as adaptability. The investigation addressed selection signatures in Gir and Tharparkar cattle, belonging to arid ecotypes of India. Double digest restriction-site associated DNA sequencing (ddRAD-seq) yielded nearly 26 million high-quality reads from unrelated seven Gir and seven Tharparkar cows. In all, 19,127 high-quality SNPs were processed for selection signature analysis. An approach involving within-population composite likelihood ratio (CLR) statistics and between-population FST statistics was used to capture selection signatures within and between the breeds, respectively. A total of 191 selection signatures were addressed using CLR and FST approaches. Selection signatures overlapping 86 and 73 genes were detected as Gir- and Tharparkar-specific, respectively. Notably, genes related to production (CACNA1D, GHRHR), reproduction (ESR1, RBMS3), immunity (NOSTRIN, IL12B) and adaptation (ADAM22, ASL) were annotated to selection signatures. Gene pathway analysis revealed genes in insulin/IGF pathway for milk production, gonadotropin releasing hormone pathway for reproduction, Wnt signalling pathway and chemokine and cytokine signalling pathway for adaptation. This is the first study where selection signatures are identified using ddRAD-seq in indicine cattle breeds. The study shall help in conservation and leveraging genetic improvements in Gir and Tharparkar cattle.

Genome-wide association study reveals 2 copy number variations associated with the variation of plumage color in the white duck hybrid population

Plumage color is an intuitive external poultry characteristic with rich manifestations and complex genetic mechanisms. In our previous study, we observed that there were more dark variations in plumage color in the F2 population derived from the hybridization of 2 white duck varieties. Therefore, based on the statistics of plumage color of 308 F2 populations, we further used the resequencing data of these individuals to detect copy number variations (CNVs) in the whole genome and conducted genome-wide association studies (GWAS) to determine the genetic basis related to plumage color traits. The CNV detection revealed 9,337 CNVs, with an average length of 15,950 bp and a total length of 142.02 MB, accounting for approximately 12.91% of the reference genome. The CNV distribution on the chromosomes was relatively uniform, and the number of CNVs on each chromosome positively correlated with the length of the chromosome. In the pure black plumage group, 2,101 CNVs were only identified, and 1,714 were specifically identified in the pure white plumage group. Ten CNVs were randomly selected for validation using quantitative real-time PCR, and 9 CNVs had the same CNV types as predicted, with an accuracy of 90%. Based on GWAS, we identified 2 CNVs potentially associated with plumage color variations, with the associated CNV regions covering 9 genes. Enrichment analysis of these 9 candidate genes showed significant enrichment of 3 pathways (ribosome biogenesis in eukaryotes, RNA transport, and protein export) and 17 gene ontology terms. Among these, VWA5A can downregulate MITF by binding to the regulatory factors SOX10. The occurrence of CNV may indirectly contribute to duck plumage color variation by affecting the regulatory factors of the switch gene MITF in the melanogenesis pathway. These findings have improved the understanding of the genetic basis of duck plumage color variation and have been beneficial for developing and using plumage color traits in subsequent poultry breeding.

Analysis of genomic copy number variations through whole-genome scan in Yunling cattle

Yunling cattle is a new breed of beef cattle bred in Yunnan Province, China, which has the advantages of fast growth, excellent meat quality, improved tolerance ability, and important landscape value. Copy number variation (CNV) is a significant source of gene structural variation and plays a crucial role in evolution and phenotypic diversity. Based on the latest reference genome ARS-UCD2.0, this study analyzed the genome-wide distribution of CNVs in Yunling cattle using short-read whole-genome sequencing data (n = 129) and single-molecule long-read sequencing data (n = 1), and a total of 16,507 CNVs were detected. After merging CNVs with overlapping genomic positions, 3,728 CNV regions (CNVRs) were obtained, accounting for 0.61% of the reference genome. The functional analysis indicated significant enrichment of CNVRs in 96 GO terms and 57 KEGG pathways, primarily related to cell adhesion, signal transduction, neuromodulation, and nutritional metabolism. Additionally, 111 CNVRs overlapped with 76 quantitative trait loci (QTLs), including Subcutaneous fat thickness QTL, Longissimus muscle area QTL, and Marbling score QTL. Several CNVR-overlapping genes, including BZW1, AOX1, and LOC100138449, overlap with regions associated with meat color and quality QTLs. Furthermore, Vst analysis showed that PSMB4, ERICH1, SMC2, and PPP4R3A were highly divergent between Yunling and Brahman cattle. In summary, we have constructed the genomic CNV map of Yunling cattle for the first time using whole-genome resequencing. This provides valuable genetic variation resources for the study of the Yunling cattle genome and contributes to the study of economic traits in Yunling cattle.

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.

Advancements in copy number variation screening in herbivorous livestock genomes and their association with phenotypic traits

Copy number variations (CNVs) have garnered increasing attention within the realm of genetics due to their prevalence in human, animal, and plant genomes. These structural genetic variations have demonstrated associations with a broad spectrum of phenotypic diversity, economic traits, environmental adaptations, epidemics, and other essential aspects of both plants and animals. Furthermore, CNVs exhibit extensive sequence variability and encompass a wide array of genomes. The advancement and maturity of microarray and sequencing technologies have catalyzed a surge in research endeavors pertaining to CNVs. This is particularly prominent in the context of livestock breeding, where molecular markers have gained prominence as a valuable tool in comparison to traditional breeding methods. In light of these developments, a contemporary and comprehensive review of existing studies on CNVs becomes imperative. This review serves the purpose of providing a brief elucidation of the fundamental concepts underlying CNVs, their mutational mechanisms, and the diverse array of detection methods employed to identify these structural variations within genomes. Furthermore, it seeks to systematically analyze the recent advancements and findings within the field of CNV research, specifically within the genomes of herbivorous livestock species, including cattle, sheep, horses, and donkeys. The review also highlighted the role of CNVs in shaping various phenotypic traits including growth traits, reproductive traits, pigmentation and disease resistance etc., in herbivorous livestock. The main goal of this review is to furnish readers with an up-to-date compilation of knowledge regarding CNVs in herbivorous livestock genomes. By integrating the latest research findings and insights, it is anticipated that this review will not only offer pertinent information but also stimulate future investigations into the realm of CNVs in livestock. In doing so, it endeavors to contribute to the enhancement of breeding strategies, genomic selection, and the overall improvement of herbivorous livestock production and resistance to diseases.

Comparative analysis of genome-wide copy number variations between Tibetan sheep and White Suffolk sheep

The DNA copy number variations (CNVs) are widely involved in affecting various kinds of biological functions, such as environmental adaptation. Tibetan sheep and White Suffolk sheep are two representative indigenous and exotic breeds raised in Sichuan, China, and both of them have many contrasting biological characteristics. In this study, we employed high-throughput sequencing approach to investigate genome-wide CNVs between the two sheep breeds. A total of 11,135 CNV regions (CNVRs) consisting of 6,488 deletions and 4,647 duplications were detected, whose length ranged from 1,599 bp to 0.56 Mb with the mean of 4,658 bp. There were 281 CNVRs segregated between Tibetan sheep and White Suffolk sheep, and 18 of them have been fixed within both breeds. Functional analyses of candidate genes within the segregating CNVRs revealed the thyroid hormone signaling pathway and CTNNB1 gene that would be responsible for differential biological characteristics of breeds, such as energy metabolism, seasonal reproduction, and litter size. Furthermore, the segregating CNVRs identified in this study were overlapped with many known quantitative trait loci that are associated with growth, testis weight, and reproductive seasonality. In conclusion, these results help us better understanding differential biological characteristics between Tibetan sheep and White Suffolk sheep.

Genome-wide identification of candidate copy number polymorphism genes associated with complex traits of Tibetan-sheep

Copy number variation (CNV) is a genetic structural polymorphism important for phenotypic diversity and important economic traits of livestock breeds, and it plays an important role in the desired genetic variation. This study used whole genome sequencing to detect the CNV variation in the genome of 6 local Tibetan sheep groups. We detected 69,166 CNV events and 7230 copy number variable regions (CNVRs) after merging the overlapping CNVs, accounting for 2.72% of the reference genome. The CNVR length detected ranged from 1.1 to 1693.5 Kb, with a total length of 118.69 Mb and an average length of 16.42 Kb per CNVR. Functional GO cluster analysis showed that the CNVR genes were mainly involved in sensory perception systems, response to stimulus, and signal transduction. Through CNVR-based Vst analysis, we found that the CACNA2D3 and CTBP1 genes related to hypoxia adaptation, the HTR1A gene related to coat color, and the TRNAS-GGA and PIK3C3 genes related to body weight were all strongly selected. The findings of our study will contribute novel insights into the genetic structural variation underlying hypoxia adaptation and economically important traits in Tibetan sheep.

Identification of copy number variation in Tibetan sheep using whole genome resequencing reveals evidence of genomic selection

BACKGROUND

Copy number variation (CNV) is an important source of structural variation in the mammalian genome. CNV assays present a new method to explore the genomic diversity of environmental adaptations in animals and plants and genes associated with complex traits. In this study, the genome-wide CNV distribution characteristics of 20 Tibetan sheep from two breeds (10 Oula sheep and 10 Panou sheep) were analysed using whole-genome resequencing to investigate the variation in the genomic structure of Tibetan sheep during breeding.

RESULTS

CNVs were detected using CNVnator, and the overlapping regions of CNVs between individual sheep were combined. Among them, a total of 60,429 CNV events were detected between the indigenous sheep breed (Oula) and the synthetic sheep breed (Panou). After merging the overlapping CNVs, 4927 CNV regions (CNVRs) were finally obtained. Of these, 4559 CNVRs were shared by two breeds, and there were 368 differential CNVRs. Deletion events have a higher percentage of occurrences than duplication events. Functional enrichment analysis showed that the shared CNVRs were significantly enriched in 163 GO terms and 62 KEGG pathways, which were mainly associated with organ development, neural regulation, immune regulation, digestion and metabolism. In addition, 140 QTLs overlapped with some of the CNVRs at more than 1 kb, such as average daily gain QTL, body weight QTL, and total lambs born QTL. Many of the CNV-overlapping genes such as PPP3CA, SSTR1 and FASN, overlap with the average daily weight gain and carcass weight QTL regions. Moreover, VST analysis showed that XIRP2, ABCB1, CA1, ASPA and EEF2 differed significantly between the synthetic breed and local sheep breed. The duplication of the ABCB1 gene may be closely related to adaptation to the plateau environment in Panou sheep, which deserves further study. Additionally, cluster analysis, based on all individuals, showed that the CNV clustering could be divided into two origins, indicating that some Tibetan sheep CNVs are likely to arise independently in different populations and contribute to population differences.

CONCLUSIONS

Collectively, we demonstrated the genome-wide distribution characteristics of CNVs in Panou sheep by whole genome resequencing. The results provides a valuable genetic variation resource and help to understand the genetic characteristics of Tibetan sheep. This study also provides useful information for the improvement and breeding of Tibetan sheep in the future.

Research progress on the regulation of production traits by gastrointestinal microbiota in dairy cows

The composition and abundance of microorganisms in the gastrointestinal tract of cows are complex and extensive, and they play a crucial role in regulating nutrient digestion, absorption, maintaining digestive tract stability, and promoting the production and health of the host. The fermentation carried out by these microorganisms in the gastrointestinal tract is fundamental to the health and productivity of cows. Rumen microorganisms produce the majority of enzymes required to break down feed substrates, such as cellulose, protein, lipids, and other plant materials, through fermentation. This process provides energy metabolism substrates that satisfy approximately 70% of the host's energy requirements for physiological activities. Gut microorganisms primarily decompose cellulose that is difficult to digest in the rumen, thereby providing heat and energy to the hosts. Additionally, they have an impact on host health and productivity through their role in immune function. Understanding the composition and function of the cow gut microbiota can help regulate dairy cattle breeding traits and improve their health status. As a result, it has become a popular research topic in dairy cattle breeding. This article provides a review of the composition, structure, physiological characteristics, and physiological effects of the cow gut microbiota, serving as a theoretical foundation for future studies that aim to utilize the gut microbiota for dairy cattle breeding or improving production traits. It may also serve as a reference for research on gut microbiota of other ruminants.

Genomic Landscape of Copy Number Variations and Their Associations with Climatic Variables in the World's Sheep

Sheep show characteristics of phenotypic diversity and adaptation to diverse climatic regions. Previous studies indicated associations between copy number variations (CNVs) and climate-driven adaptive evolution in humans and other domestic animals. Here, we constructed a genomic landscape of CNVs (n = 39,145) in 47 old autochthonous populations genotyped at a set of high-density (600 K) SNPs to detect environment-driven signatures of CNVs using a multivariate regression model. We found 136 deletions and 52 duplications that were significantly (Padj. < 0.05) associated with climatic variables. These climate-mediated selective CNVs are involved in functional candidate genes for heat stress and cold climate adaptation (e.g., B3GNTL1, UBE2L3, and TRAF2), coat and wool-related traits (e.g., TMEM9, STRA6, RASGRP2, and PLA2G3), repairing damaged DNA (e.g., HTT), GTPase activity (e.g., COPG), fast metabolism (e.g., LMF2 and LPIN3), fertility and reproduction (e.g., SLC19A1 and CCDC155), growth-related traits (e.g., ADRM1 and IGFALS), and immune response (e.g., BEGAIN and RNF121) in sheep. In particular, we identified significant (Padj. < 0.05) associations between probes in deleted/duplicated CNVs and solar radiation. Enrichment analysis of the gene sets among all the CNVs revealed significant (Padj. < 0.05) enriched gene ontology terms and pathways related to functions such as nucleotide, protein complex, and GTPase activity. Additionally, we observed overlapping between the CNVs and 140 known sheep QTLs. Our findings imply that CNVs can serve as genomic markers for the selection of sheep adapted to specific climatic conditions.

A comprehensive analysis of the genetic diversity and environmental adaptability in worldwide Merino and Merino-derived sheep breeds

BACKGROUND

To enhance and extend the knowledge about the global historical and phylogenetic relationships between Merino and Merino-derived breeds, 19 populations were genotyped with the OvineSNP50 BeadChip specifically for this study, while an additional 23 populations from the publicly available genotypes were retrieved. Three complementary statistical tests, Rsb (extended haplotype homozygosity between-populations), XP-EHH (cross-population extended haplotype homozygosity), and runs of homozygosity (ROH) islands were applied to identify genomic variants with potential impact on the adaptability of Merino genetic type in two contrasting climate zones.

RESULTS

The results indicate that a large part of the Merino's genetic relatedness and admixture patterns are explained by their genetic background and/or geographic origin, followed by local admixture. Multi-dimensional scaling, Neighbor-Net, Admixture, and TREEMIX analyses consistently provided evidence of the role of Australian, Rambouillet and German strains in the extensive gene introgression into the other Merino and Merino-derived breeds. The close relationship between Iberian Merinos and other South-western European breeds is consistent with the Iberian origin of the Merino genetic type, with traces from previous contributions of other Mediterranean stocks. Using Rsb and XP-EHH approaches, signatures of selection were detected spanning four genomic regions located on Ovis aries chromosomes (OAR) 1, 6 and 16, whereas two genomic regions on OAR6, that partially overlapped with the previous ones, were highlighted by ROH islands. Overall, the three approaches identified 106 candidate genes putatively under selection. Among them, genes related to immune response were identified via the gene interaction network. In addition, several candidate genes were found, such as LEKR1, LCORL, GHR, RBPJ, BMPR1B, PPARGC1A, and PRKAA1, related to morphological, growth and reproductive traits, adaptive thermogenesis, and hypoxia responses.

CONCLUSIONS

To the best of our knowledge, this is the first comprehensive dataset that includes most of the Merino and Merino-derived sheep breeds raised in different regions of the world. The results provide an in-depth picture of the genetic makeup of the current Merino and Merino-derived breeds, highlighting the possible selection pressures associated with the combined effect of anthropic and environmental factors. The study underlines the importance of Merino genetic types as invaluable resources of possible adaptive diversity in the context of the occurring climate changes.

Resequencing of a Pekin duck breeding population provides insights into the genomic response to short-term artificial selection

BACKGROUND

Short-term, intense artificial selection drives fast phenotypic changes in domestic animals and leaves imprints on their genomes. However, the genetic basis of this selection response is poorly understood. To better address this, we employed the Pekin duck Z2 pure line, in which the breast muscle weight was increased nearly 3-fold after 10 generations of breeding. We denovo assembled a high-quality reference genome of a female Pekin duck of this line (GCA_003850225.1) and identified 8.60 million genetic variants in 119 individuals among 10 generations of the breeding population.

RESULTS

We identified 53 selected regions between the first and tenth generations, and 93.8% of the identified variations were enriched in regulatory and noncoding regions. Integrating the selection signatures and genome-wide association approach, we found that 2 regions covering 0.36 Mb containing UTP25 and FBRSL1 were most likely to contribute to breast muscle weight improvement. The major allele frequencies of these 2 loci increased gradually with each generation following the same trend. Additionally, we found that a copy number variation region containing the entire EXOC4 gene could explain 1.9% of the variance in breast muscle weight, indicating that the nervous system may play a role in economic trait improvement.

CONCLUSIONS

Our study not only provides insights into genomic dynamics under intense artificial selection but also provides resources for genomics-enabled improvements in duck breeding.

Copy Number Variation Analysis Revealed the Evolutionary Difference between Chinese Indigenous Pigs and Asian Wild Boars

Copy number variation (CNV) has been widely used to study the evolution of different species. We first discovered different CNVs in 24 Anqingliubai pigs and 6 Asian wild boars using next-generation sequencing at the whole-genome level with 10× depth to understand the relationship between genetic evolution and production traits in wild boars and domestic pigs. A total of 97,489 CNVs were identified and divided into 10,429 copy number variation regions (CNVRs), occupying 32.06% of the porcine genome. Chromosome 1 had the most CNVRs, and chromosome 18 had the least. Ninety-six CNVRs were selected using VST 1% based on the signatures of all CNVRs, and sixty-five genes were identified in the selected regions. These genes were strongly correlated with traits distinguishing groups by enrichment in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways, such as growth (CD36), reproduction (CIT, RLN), detoxification (CYP3A29), and fatty acid metabolism (ELOVL6). The QTL overlapping regions were associated with meat traits, growth, and immunity, which was consistent with CNV analysis. Our findings increase the understanding of evolved genome structural variations between wild boars and domestic pigs, and provide new molecular biomarkers to guide breeding and the efficient use of available genetic resources.

In it for the long run: perspectives on exploiting long-read sequencing in livestock for population scale studies of structural variants

Studies have demonstrated that structural variants (SV) play a substantial role in the evolution of species and have an impact on Mendelian traits in the genome. However, unlike small variants (< 50 bp), it has been challenging to accurately identify and genotype SV at the population scale using short-read sequencing. Long-read sequencing technologies are becoming competitively priced and can address several of the disadvantages of short-read sequencing for the discovery and genotyping of SV. In livestock species, analysis of SV at the population scale still faces challenges due to the lack of resources, high costs, technological barriers, and computational limitations. In this review, we summarize recent progress in the characterization of SV in the major livestock species, the obstacles that still need to be overcome, as well as the future directions in this growing field. It seems timely that research communities pool resources to build global population-scale long-read sequencing consortiums for the major livestock species for which the application of genomic tools has become cost-effective.

Genome-wide analysis of CNVs in three populations of Tibetan sheep using whole-genome resequencing

Copy number variation (CNV), an important source of genomic structural variation, can disturb genetic structure, dosage, regulation and expression, and is associated with phenotypic diversity and adaptation to local environments in mammals. In the present study, 24 resequencing datasets were used to characterize CNVs in three ecotypic populations of Tibetan sheep and assess CNVs related to domestication and adaptation in Qinghai-Tibetan Plateau. A total of 87,832 CNV events accounting for 0.3% of the sheep genome were detected. After merging the overlapping CNVs, 2777 CNV regions (CNVRs) were obtained, among which 1098 CNVRs were shared by the three populations. The average length of these CNVRs was more than 3 kb, and duplication events were more frequent than deletions. Functional analysis showed that the shared CNVRs were significantly enriched in 56 GO terms and 18 KEGG pathways that were mainly concerned with ABC transporters, olfactory transduction and oxygen transport. Moreover, 188 CNVRs overlapped with 97 quantitative trait loci (QTLs), such as growth and carcass QTLs, immunoglobulin QTLs, milk yield QTLs and fecal egg counts QTLs. PCDH15, APP and GRID2 overlapped with body weight QTLs. Furthermore, Vst analysis showed that RUNX1, LOC101104348, LOC105604082 and PAG11 were highly divergent between Highland-type Tibetan Sheep (HTS) and Valley-type Tibetan sheep (VTS), and RUNX1 and LOC101111988 were significantly differentiated between VTS and Oura-type Tibetan sheep (OTS). The duplication of RUNX1 may facilitate the hypoxia adaptation of OTS and HTS in Qinghai-Tibetan Plateau, which deserves further research in detail. In conclusion, for the first time, we represented the genome-wide distribution characteristics of CNVs in Tibetan sheep by resequencing, and provided a valuable genetic variation resource, which will facilitate the elucidation of the genetic basis underlying the distinct phenotypic traits and local adaptation of Tibetan sheep.

Genome-wide evaluation of copy gain and loss variations in three Afghan sheep breeds

Copy number variation (CNV) is one of the main sources of variation between different individuals that has recently attracted much researcher interest as a major source for heritable variation in complex traits. The aim of this study was to identify CNVs in Afghan indigenous sheep consisting of three Arab, Baluchi, and Gadik breeds using genomic arrays containing 53,862 single nucleotide polymorphism (SNP) markers. Data were analyzed using the Hidden Markov Model (HMM) of PennCNV software. In this study, out of 45 sheep studied, 97.8% (44 animals) have shown CNVs. In total, 411 CNVs were observed for autosomal chromosomes and the entire sequence length of around 144 Mb was identified across the genome. The average number of CNVs per each sheep was 9.13. The identified CNVs for Arab, Baluchi, and Gadik breeds were 306, 62, and 43, respectively. After merging overlapped regions, a total of 376 copy number variation regions (CNVR) were identified, which are 286, 50, and 40 for Arab, Baluchi, and Gadik breeds, respectively. Bioinformatics analysis was performed to identify the genes and QTLs reported in these regions and the biochemical pathways involved by these genes. The results showed that many of these CNVRs overlapped with the genes or QTLs that are associated with various pathways such as immune system development, growth, reproduction, and environmental adaptions. Furthermore, to determine a genome-wide pattern of selection signatures in Afghan sheep breeds, the unbiased estimates of FST was calculated and the results indicated that 37 of the 376 CNVRs (~ 10%) have been also under selection signature, most of those overlapped with the genes influencing production, reproduction and immune system. Finally, the statistical methods used in this study was applied in an external dataset including 96 individuals of the Iranian sheep breed. The results indicated that 20 of the 114 CNVRs (18%) identified in Iranian sheep breed were also identified in our study, most of those overlapped with the genes influencing production, reproduction and immune system. Overall, this is the first attempts to develop the genomic map of loss and gain variation in the genome of Afghan indigenous sheep breeds, and may be important to shed some light on the genomic regions associated with some economically important traits in these breeds.

Prenatal diagnosis of a novel 7q31.31q31.33 microduplication with a favorable outcome

BACKGROUND

Copy number variants (CNVs) are an important source of normal and pathogenic genome variations. Especially CNVs identified in prenatal cases need careful considerations and correct interpretation if those are harmless or harmful variants from the norm.

CASE PRESENTATION

Herein, we reported a paternally inherited duplication of 7.6 Mb in 7q31.3 with, surprisingly, a favorable outcome. GTG-banding and CMA on the DNA derived from uncultured amniocytes revealed a karyotype: 46,XX.arr[GRCh37] 7q31.31q31.33(118,601,001_126,177,044) × 3. Ultrasound examination showed no dysmorphisms or intrauterine growth restriction in the fetus and the father was clinically normal as well.

CONCLUSION

Prenatal detection of a 7.6 Mb in 7q31.31 to 7q31.33 duplication in a female fetus turned out to be a yet unreported unbalanced chromosome abnormality. This is another example that parental testing and GTG-banding are necessary additional tests to be done in prenatal cases, before a reliable conclusion on the meaning of an aberration can be drawn.

Whole-genome sequence analysis unveils different origins of European and Asiatic mouflon and domestication-related genes in sheep

The domestication and subsequent development of sheep are crucial events in the history of human civilization and the agricultural revolution. However, the impact of interspecific introgression on the genomic regions under domestication and subsequent selection remains unclear. Here, we analyze the whole genomes of domestic sheep and their wild relative species. We found introgression from wild sheep such as the snow sheep and its American relatives (bighorn and thinhorn sheep) into urial, Asiatic and European mouflons. We observed independent events of adaptive introgression from wild sheep into the Asiatic and European mouflons, as well as shared introgressed regions from both snow sheep and argali into Asiatic mouflon before or during the domestication process. We revealed European mouflons might arise through hybridization events between a now extinct sheep in Europe and feral domesticated sheep around 6000-5000 years BP. We also unveiled later introgressions from wild sheep to their sympatric domestic sheep after domestication. Several of the introgression events contain loci with candidate domestication genes (e.g., PAPPA2, NR6A1, SH3GL3, RFX3 and CAMK4), associated with morphological, immune, reproduction or production traits (wool/meat/milk). We also detected introgression events that introduced genes related to nervous response (NEURL1), neurogenesis (PRUNE2), hearing ability (USH2A), and placental viability (PAG11 and PAG3) into domestic sheep and their ancestral wild species from other wild species.