Genome-wide detection of copy number variation in American mink using whole-genome sequencing

Whole-genome sequencing of copy number variation analysis in Ethiopian cattle reveals adaptations to diverse environments

BACKGROUND

Genomic structural variations (GSVs), notably copy number variations (CNVs), significantly shape genetic diversity and facilitate adaptation in cattle populations. Despite their importance, the genome-wide characterization of CNVs in indigenous Ethiopian cattle breeds-Abigar, Fellata, and Gojjam-Highland remains largely unexplored. In this study, we applied a read-depth approach to whole genome sequencing (WGS) data to conduct the first comprehensive analysis of CNVs in these populations.

RESULTS

We identified 3,893 CNV regions (CNVRs) covering 19.15 Mb (0.71% of the cattle genome). These CNVRs ranged from 1.60 kb to 488.0 kb, with an average size of 4.92 kb. These CNVRs included deletions (1713), duplications (1929), and mixed events (251) showing notable differences in distribution among the breeds. Four out of five randomly selected CNVRs were successfully validated using real time polymerase chain reaction (qPCR). Further analyses identified candidate genes associated with high-altitude adaptation (GBE1 and SOD1), heat stress adaptation (HSPA13, DNAJC18, and DNAJC8) and resistance to tick infestations (BoLA and KRT33A). In addition, variance stabilizing transformation (VST) statistics highlighted population-specific CNVRs, emphasizing the unique genetic signatures of high-altitude adaptation in the Gojjam-Highland cattle breed. Among the detected CNVRs, 4.93% (192 out of 3,893) overlapped with 520 quantitative traits loci (QTLs) associated with six economically important trait categories suggesting that these CNVRs may significantly contribute to the genetic variation underlying these traits.

CONCLUSIONS

Our comprehensive analysis reveals significant CNVRs associated with key adaptive traits in Ethiopian cattle breeds highlighting their genetic diversity and resilience. These findings offer valuable insights into the genetic basis of adaptability and can inform sustainable breeding practices and conservation efforts. Future research should prioritize the functional validation of these CNVRs and their integration into breeding programs to enhance traits such as disease resistance and environmental adaptability.

A comprehensive map of copy number variations in dromedary camels based on whole genome sequence data

Copy number variants (CNVs) are structural variants within the eukaryotic genome that vary among individuals of a species. These variants have been associated with different phenotypic traits, making them a valuable consideration as markers for designing breeding programmes. In this study, whole genome sequence data of 60 dromedary camel samples originating from the Arabian Peninsula were analyzed to construct a comprehensive dromedary CNV map. Utilizing four CNV callers employing read-depth, split-read and paired-end mapping approaches, a total of 37,519 CNV events (17,847 deletions and 19,672 duplications) were called on the dromedary autosomes. These CNV events were merged into 2,557 regions, categorized as 1,322 losses, 122 gains, and 1,113 "mixed regions" comprising both types. The cumulative size of the CNV regions amounted to 22.5 Mb, covering roughly 1.16% of the dromedary autosomes. Approximately 32% of the defined CNV regions (comprising 60% losses, 18% gains, and 0.27% mixed regions) were found in ≥ 90% of the dromedary samples, classifying them as prevalent regions. Genes with biological functions related to the different adaptive physiologies of dromedary camels, such as fertility, heat stress, musculoskeletal development, and fat metabolism, were overlapping with or in close proximity to ~ 68% of the defined CNV regions, demonstrating their potential role in dromedaries' physiology. This study presents the first comprehensive CNV map of dromedary camels and builds on the present knowledge in understanding the genetic structure of this species.

Genome-wide copy number variant analysis reveals candidate genes associated with milk production traits in water buffalo (Bubalus bubalis)

Buffaloes are vital contributors to the global dairy industry. Understanding the genetic basis of milk production traits in buffalo populations is essential for breeding programs and improving productivity. In this study, we conducted whole-genome resequencing on 387 buffalo genomes from 29 diverse Asian breeds, including 132 river buffaloes, 129 swamp buffaloes, and 126 crossbred buffaloes. We identified 36,548 copy number variants (CNV) spanning 133.29 Mb of the buffalo genome, resulting in 2,100 CNV regions (CNVR), with 1,993 shared CNVR being found within the studied buffalo types. Analyzing CNVR highlighted distinct genetic differentiation between river and swamp buffalo subspecies, verified by evolutionary tree and principal component analyses. Admixture analysis grouped buffaloes into river and swamp categories, with crossbred buffaloes displaying mixed ancestry. To identify candidate genes associated with milk production traits, we employed 3 approaches. First, we used Vst-based population differentiation, revealing 11 genes within CNVR that exhibited significant divergence between different buffalo breeds, including genes linked to milk production traits. Second, expression quantitative loci analysis revealed differentially expressed CNVR-derived genes (DECG) associated with milk production traits. Notably, known milk production-related genes were among these DECG, validating their relevance. Last, a GWAS identified 3 CNVR significantly linked to peak milk yield. Our study provides comprehensive genomic insights into buffalo populations and identifies candidate genes associated with milk production traits. These findings facilitate genetic breeding programs aimed at increasing milk yield and improving quality in this economically important livestock species.

Genomic Regions Associated with Growth and Reproduction Traits in Pink-Eyed White Mink

In mink breeding, balanced selection for growth and reproductive features is essential because these traits are contradictory. The variables of total number born (TNB), number born alive (NBA), and body weight (BW) are highly valuable in terms of their importance in mink production. A comprehensive understanding of the molecular mechanisms that drive these features could offer vital insights into their genetic compositions. In the present study, the single-nucleotide polymorphism (SNP) genotypes of 219 minks were obtained via double digest restriction-site associated DNA sequencing (ddRAD-seq). Following several rounds of screening, about 2,415,121 high-quality SNPs were selected for a genome-wide association study (GWAS). The GWAS was used to determine BW and reproductive traits in pink-eyed white mink. It was suggested that SLC26A36, STXBP5L, and RPS 29 serve as potential genes for the total number of kits born (TNB), while FSCB, PDPN, NKX 2-1, NFKB 1, NFKBIA, and GABBR1 are key genes for the number born alive (NBA). Moreover, RTTN, PRPF31, MACROD1, and KYAT1 are possible BW genes based on association results and available functional data from gene and mammalian phenotype databases. These results offer essential information about the variety of mink and theoretical principles for applying mink breeds.

Genome-wide association studies for economically important traits in mink using copy number variation

Copy number variations (CNVs) are structural variants consisting of duplications and deletions of DNA segments, which are known to play important roles in the genetics of complex traits in livestock species. However, CNV-based genome-wide association studies (GWAS) have remained unexplored in American mink. Therefore, the purpose of the current study was to investigate the association between CNVs and complex traits in American mink. A CNV-based GWAS was performed with the ParseCNV2 software program using deregressed estimated breeding values of 27 traits as pseudophenotypes, categorized into traits of growth and feed efficiency, reproduction, pelt quality, and Aleutian disease tests. The study identified a total of 10,137 CNVs (6968 duplications and 3169 deletions) using the Affymetrix Mink 70K single nucleotide polymorphism (SNP) array in 2986 American mink. The association analyses identified 250 CNV regions (CNVRs) associated with at least one of the studied traits. These CNVRs overlapped with a total of 320 potential candidate genes, and among them, several genes have been known to be related to the traits such as ARID1B, APPL1, TOX, and GPC5 (growth and feed efficiency traits); GRM1, RNASE10, WNT3, WNT3A, and WNT9B (reproduction traits); MYO10, and LIMS1 (pelt quality traits); and IFNGR2, APEX1, UBE3A, and STX11 (Aleutian disease tests). Overall, the results of the study provide potential candidate genes that may regulate economically important traits and therefore may be used as genetic markers in mink genomic breeding programs.

Read-depth based approach on whole genome resequencing data reveals important insights into the copy number variation (CNV) map of major global buffalo breeds

BACKGROUND

Elucidating genome-wide structural variants including copy number variations (CNVs) have gained increased significance in recent times owing to their contribution to genetic diversity and association with important pathophysiological states. The present study aimed to elucidate the high-resolution CNV map of six different global buffalo breeds using whole genome resequencing data at two coverages (10X and 30X). Post-quality control, the sequence reads were aligned to the latest draft release of the Bubaline genome. The genome-wide CNVs were elucidated using a read-depth approach in CNVnator with different bin sizes. Adjacent CNVs were concatenated into copy number variation regions (CNVRs) in different breeds and their genomic coverage was elucidated.

RESULTS

Overall, the average size of CNVR was lower at 30X coverage, providing finer details. Most of the CNVRs were either deletion or duplication type while the occurrence of mixed events was lesser in number on a comparative basis in all breeds. The average CNVR size was lower at 30X coverage (0.201 Mb) as compared to 10X (0.013 Mb) with the finest variants in Banni buffaloes. The maximum number of CNVs was observed in Murrah (2627) and Pandharpuri (25,688) at 10X and 30X coverages, respectively. Whereas the minimum number of CNVs were scored in Surti at both coverages (2092 and 17,373). On the other hand, the highest and lowest number of CNVRs were scored in Jaffarabadi (833 and 10,179 events) and Surti (783 and 7553 events) at both coverages. Deletion events overnumbered duplications in all breeds at both coverages. Gene profiling of common overlapped genes and longest CNVRs provided important insights into the evolutionary history of these breeds and indicate the genomic regions under selection in respective breeds.

CONCLUSION

The present study is the first of its kind to elucidate the high-resolution CNV map in major buffalo populations using a read-depth approach on whole genome resequencing data. The results revealed important insights into the divergence of major global buffalo breeds along the evolutionary timescale.

Rosemary extract improves egg quality by altering gut barrier function, intestinal microbiota and oviductal gene expressions in late-phase laying hens

BACKGROUND

Rosemary extract (RE) has been reported to exert antioxidant property. However, the application of RE in late-phase laying hens on egg quality, intestinal barrier and microbiota, and oviductal function has not been systematically studied. This study was investigated to detect the potential effects of RE on performance, egg quality, serum parameters, intestinal heath, cecal microbiota and metabolism, and oviductal gene expressions in late-phase laying hens. A total of 210 65-week-old "Jing Tint 6" laying hens were randomly allocated into five treatments with six replicates and seven birds per replicate and fed basal diet (CON) or basal diet supplemented with chlortetracycline at 50 mg/kg (CTC) or RE at 50 mg/kg (RE50), 100 mg/kg (RE100), and 200 mg/kg (RE200).

RESULTS

Our results showed that RE200 improved (P < 0.05) Haugh unit and n-6/n-3 of egg yolk, serum superoxide dismutase (SOD) compared with CON. No significant differences were observed for Haugh unit and n-6/n-3 of egg yolk among CTC, RE50, RE100 and RE200 groups. Compared with CTC and RE50 groups, RE200 increased serum SOD activity on d 28 and 56. Compared with CON, RE supplementation decreased (P < 0.05) total cholesterol (TC) level. CTC, RE100 and RE200 decreased (P < 0.05) serum interleukin-6 (IL-6) content compared with CON. CTC and RE200 increased jejunal mRNA expression of ZO-1 and Occludin compared with CON. The biomarkers of cecal microbiota and metabolite induced by RE 200, including Firmicutes, Eisenbergiella, Paraprevotella, Papillibacter, and butyrate, were closely associated with Haugh unit, n-6/n-3, SOD, IL-6, and TC. PICRUSt2 analysis indicated that RE altered carbohydrate and amino acid metabolism of cecal microbiota and increased butyrate synthesizing enzymes, including 3-oxoacid CoA-transferase and butyrate-acetoacetate CoA-transferase. Moreover, transcriptomic analysis revealed that RE200 improved gene expressions and functional pathways related to immunity and albumen formation in the oviductal magnum.

CONCLUSIONS

Dietary supplementation with 200 mg/kg RE could increase egg quality of late-phase laying hens via modulating intestinal barrier, cecal microbiota and metabolism, and oviductal function. Overall, RE could be used as a promising feed additive to improve egg quality of laying hens at late stage of production.

Population Structure and Selection Signatures Underlying Domestication Inferred from Genome-Wide Copy Number Variations in Chinese Indigenous Pigs

Single nucleotide polymorphism was widely used to perform genetic and evolution research in pigs. However, little is known about the effect of copy number variation (CNV) on characteristics in pigs. This study performed a genome-wide comparison of CNVs between Wannan black pigs (WBP) and Asian wild boars (AWB), using whole genome resequencing data. By using Manta, we detected in total 28,720 CNVs that covered approximately 1.98% of the pig genome length. We identified 288 selected CNVs (top 1%) by performing Fst statistics. Functional enrichment analyses for genes located in selected CNVs were found to be muscle related (NDN, TMOD4, SFRP1, and SMYD3), reproduction related (GJA1, CYP26B1, WNT5A, SRD5A2, PTPN11, SPEF2, and CCNB1), residual feed intake (RFI) related (MAP3K5), and ear size related (WIF1). This study provides essential information on selected CNVs in Wannan black pigs for further research on the genetic basis of the complex phenotypic and provides essential information for direction in the protection and utilization of Wannan black pig.