Search

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. [ABSTRACT FROM AUTHOR]

Background: Coffin-Siris syndrome is a clinically elusive and rare genetic disease characterized by a wide range of clinical manifestations. This study deeply analyzed and identified the clinical phenotype and genetic variant location in a pediatric patient with Coffin-Siris syndrome, aiming to enhance the understanding of this syndrome and assist in its screening and diagnosis. Methods: A combination of advanced diagnostic tools, including high-throughput whole-exome sequencing (WES) and first-generation sequencing technologies, was employed to ascertain the etiology of the disease in the child. Results: The clinical phenotype was characterized by stunted growth, reduced stature, spina bifida, enuresis, and a ventricular septal defect. WES revealed a de novo variant in the SOX11 gene locus (c.700G > T), identified as pathogenic. It is noteworthy that this variant has not been previously reported. Conclusions: The combination of clinical presentation and genetic testing results supports that the patient suffers from Coffin-Siris syndrome due to a genetic variant in the SOX11 gene. This de novo variant expands our understanding of human gene variation, which is conducive to genetic counseling and screening for early diagnosis of Coffin-Siris syndrome. [ABSTRACT FROM AUTHOR]

The adaptation of marine fish to freshwater environments includes prodigious examples of rapid evolution.Given the scarcity of ω‐3 long‐chain polyunsaturated fatty acids (LC‐PUFA) in freshwater, we expect selection to be strong on fish transitioning to freshwater habitats and yet the underlying ecological causes of genomic and phenotypic differentiation are poorly understood for traits associated with lipid content and composition. Threespine stickleback Gasterosteus aculeatus have repeatedly colonized, and adapted to, freshwater habitats across the Northern Hemisphere. These freshwater populations often show elevated copy number of the fatty acid desaturase 2 gene (FADS2), which increases the biosynthetic capacity of LC‐PUFA. The starkly lower content of LC‐PUFA in freshwater compared to marine prey, especially docosahexaenoic acid (DHA), likely imposes strong positive selection on freshwater fish for either increased biosynthesis or greater dietary acquisition of LC‐PUFA. The recently colonized and relatively undisturbed threespine stickleback populations in postglacial coastal lakes of southern Greenland offer an exceptional opportunity to study how variation in the copy number of FADS2 is related to abiotic and biotic conditions of lakes and their morphometry. As expected, given its position on the stickleback X chromosome, we found strong sexual dimorphism in FADS2 copy number in all populations (19 freshwater, 1 marine and 1 brackish), and an increased dimorphism in some freshwater populations. We also found that FADS2 copy number was negatively correlated, for both males and females, with the abundance of copepods, which are a DHA‐rich food source in the zooplankton community. Overall, our results suggest that the prey community context of lakes might influence the process of metabolic adaptation of marine fish colonizing freshwater ecosystems. [ABSTRACT FROM AUTHOR]

The abundance of Lp(a) protein holds significant implications for the risk of cardiovascular disease (CVD), which is directly impacted by the copy number (CN) of KIV-2, a 5.5 kbp sub-region. KIV-2 is highly polymorphic in the population and accurate analysis is challenging. In this study, we present the DRAGEN KIV-2 CN caller, which utilizes short reads. Data across 166 WGS show that the caller has high accuracy, compared to optical mapping and can further phase approximately 50% of the samples. We compared KIV-2 CN numbers to 24 previously postulated KIV-2 relevant SNVs, revealing that many are ineffective predictors of KIV-2 copy number. Population studies, including USA-based cohorts, showed distinct KIV-2 CN, distributions for European-, African-, and Hispanic-American populations and further underscored the limitations of SNV predictors. We demonstrate that the CN estimates correlate significantly with the available Lp(a) protein levels and that phasing is highly important. [ABSTRACT FROM AUTHOR]

Human salivary amylase secretion increases in response to stress; the activity has been reported to rise significantly with high‐intensity exercise. The human salivary amylase gene (AMY1) has copy number variation, with the copy number correlating with salivary amylase activity. However, the relationship between individual AMY1 copy number and salivary amylase activity in response to exercise remains unclear. In this study, we investigated AMY1 copy number and fluctuations in amylase activity in 42 healthy university students (25 males and 17 females). Participants engaged in intermittent round‐trip interval training on a basketball court. Saliva samples were collected pre‐ and post‐exercise to measure amylase activity. DNA was extracted from the oral mucosa, and AMY1 copy number was quantified using RT‐PCR. Results showed a significant increase in amylase activity postexercise. Additionally, amylase activity pre‐ and post‐exercise was positively correlated with AMY1 copy number. The generalize linear model showed that the exercise‐induced increase in amylase activity per AMY1 gene was negatively related to the AMY1 copy number and aerobic fitness. Gender has no effect on amylase activity. These results suggest a different mechanism for the constitutive and exercise‐induced amylase secretion, while aerobic fitness may be independently involved in the secretion. [ABSTRACT FROM AUTHOR]

Background: Structural variations (SVs) are widespread across genome and have a great impact on evolution, disease, and phenotypic diversity. Despite the development of numerous bioinformatic tools, commonly referred to as SV callers, tailored for detecting SVs using whole genome sequence (WGS) data and employing diverse algorithms, their performance necessitates rigorous evaluation with real data and validated SVs. Moreover, a considerable proportion of these tools have been primarily designed and optimized using human genome data. Consequently, their applicability and performance in Avian species, characterized by smaller genomes and distinct genomic architectures, remain inadequately assessed. Results: We performed a comprehensive assessment of the performance of ten widely used SV callers using population-level real genomic data with the validated five common types of SVs. The performance of SV callers varies with the types and sizes of SVs. As compared with other tools, GRIDSS, Lumpy, Wham, and Manta present better detection accuracy. Pindel can detect more small SVs than others. CNVnator and CNVkit can detect more medium and large copy number variations. Given the poor consistency among different SV callers, the combination calling strategy is not recommended. All tools show poor ability in the detection of insertions (especially with size > 150 bp). At least 50× read depth is required to detect more than 80% of the SVs for most tools. Conclusions: This study highlights the importance and necessity of using real sequencing data, rather than simulated data only, with validated SVs for SV caller evaluation. Some practical guidance and suggestions are provided for SV detection in future researches. [ABSTRACT FROM AUTHOR]

Background: Repetitive genome regions, such as variable number of tandem repeats (VNTR) or short tandem repeats (STR), are major constituents of the uncharted dark genome and evade conventional sequencing approaches. The protein-coding LPA kringle IV type-2 (KIV-2) VNTR (5.6 kb per unit, 1–40 units per allele) is a medically highly relevant example with a particularly intricate structure, multiple haplotypes, intragenic homologies, and an intra-VNTR STR. It is the primary regulator of plasma lipoprotein(a) [Lp(a)] concentrations, an important cardiovascular risk factor. Lp(a) concentrations vary widely between individuals and ancestries. Multiple variants and functional haplotypes in the LPA gene and especially in the KIV-2 VNTR strongly contribute to this variance. Methods: We evaluated the performance of amplicon-based nanopore sequencing with unique molecular identifiers (UMI-ONT-Seq) for SNP detection, haplotype mapping, VNTR unit consensus sequence generation, and copy number estimation via coverage-corrected haplotypes quantification in the KIV-2 VNTR. We used 15 human samples and low-level mixtures (0.5 to 5%) of KIV-2 plasmids as a validation set. We then applied UMI-ONT-Seq to extract KIV-2 VNTR haplotypes in 48 multi-ancestry 1000 Genome samples and analyzed at scale a poorly characterized STR within the KIV-2 VNTR. Results: UMI-ONT-Seq detected KIV-2 SNPs down to 1% variant level with high sensitivity, specificity, and precision (0.977 ± 0.018; 1.000 ± 0.0005; 0.993 ± 0.02) and accurately retrieved the full-length haplotype of each VNTR unit. Human variant levels were highly correlated with next-generation sequencing (R2 = 0.983) without bias across the whole variant level range. Six reads per UMI produced sequences of each KIV-2 unit with Q40 quality. The KIV-2 repeat number determined by coverage-corrected unique haplotype counting was in close agreement with droplet digital PCR (ddPCR), with 70% of the samples falling even within the narrow confidence interval of ddPCR. We then analyzed 62,679 intra-KIV-2 STR sequences and explored KIV-2 SNP haplotype patterns across five ancestries. Conclusions: UMI-ONT-Seq accurately retrieves the SNP haplotype and precisely quantifies the VNTR copy number of each repeat unit of the complex KIV-2 VNTR region across multiple ancestries. This study utilizes the KIV-2 VNTR, presenting a novel and potent tool for comprehensive characterization of medically relevant complex genome regions at scale. [ABSTRACT FROM AUTHOR]

Background: Retinal dystrophies (RDs) are the most common cause of inherited blindness worldwide and are caused by genetic defects in about 300 different genes. While targeted next‐generation sequencing (NGS) has been demonstrated to be a reliable and efficient method to identify RD disease‐causing variants, it doesn't routinely identify pathogenic structural variant as copy number variations (CNVs). Targeted NGS‐based CNV detection has become a crucial step for RDs molecular diagnosis, particularly in cases without identified causative single nucleotide or Indels variants. Herein, we report the exome sequencing (ES) data‐based read‐depth bioinformatic analysis in a group of 30 unrelated Mexican RD patients with a negative or inconclusive genetic result after ES. Methods: CNV detection was performed using ExomeDepth software, an R package designed to detect CNVs using exome data. Bioinformatic validation of identified CNVs was conducted through a commercially available CNV caller. All identified candidate pathogenic CNVs were orthogonally verified through quantitative PCR assays. Results: Pathogenic or likely pathogenic CNVs were identified in 6 out of 30 cases (20%), and of them, a definitive molecular diagnosis was reached in 5 cases, for a final diagnostic rate of ~17%. CNV‐carrying genes included CLN3 (2 cases), ABCA4 (novel deletion), EYS, and RPGRIP1. Conclusions: Our results indicate that bioinformatic analysis of ES data is a reliable method for pathogenic CNV detection and that it should be incorporated in cases with a negative or inconclusive molecular result after ES. [ABSTRACT FROM AUTHOR]

In livestock breeding, single nucleotide polymorphism arrays have become a cornerstone of modern livestock breeding. SNP arrays facilitate the identification of genetic markers linked to economically important traits and provide a powerful tool for predicting breeding values. However, conventional breeding programs often overlook additional genomic features contained in the SNP array data that can provide valuable insights into the genetic diversity, copy number variation, inbreeding levels and potential challenges in breeding lines. Here we present genomeprofile, a tool using SNP array‐based genomic data, offering a comprehensive profile of breeding animals including the identification of copy number variants and runs of homozygosity, and screening for aneuploidy. By integrating these features into the breeding landscape, genomeprofile enables a more comprehensive picture of genomic variation, ultimately enhancing precision breeding strategies. To illustrate the practicality and efficacy of genomeprofile, we applied the tool to a dataset of four pig breeding lines. The genomeprofile tool is a user‐friendly tool that processes genotype data in finalreport or plink ped format efficiently into useful output. The output contains copy number variations, runs of homozygosity, selection signatures, aneuploidy and inbreeding per individual and across populations. This allows breeding companies and researchers to identify unique individuals or regions in the genome of interest based on routinely collected data. [ABSTRACT FROM AUTHOR]

Objective: This study aimed to explore the heterogeneity of tumor endothelial cells (TECs) in hepatocellular carcinoma (HCC) and their role in tumor progression, with the goal of identifying new therapeutic targets and strategies to improve patient prognosis. Methods: Single-cell RNA sequencing data from nine primary liver cancer samples were analyzed, obtained from the Gene Expression Omnibus (GEO) database. Data preprocessing, normalization, dimensionality reduction, and batch effect correction were performed based on the Seurat package. HCC cell types were identified using uniform manifold approximation and projection (UMAP) and cluster analysis, and the different cell types were annotated using the CellMarker database. Pseudotime trajectory analysis was conducted with Monocle to explore the differentiation trajectory of TECs. MAPK signaling pathway activity and copy number variations (CNV) in TECs were analyzed in conjunction with data from The Cancer Genome Atlas (TCGA), the trans-well and wound healing assay was used for cell invasion and migration activity assessment. Results: Two subgroups of TECs (TECs 1 and TECs 2) were identified, exhibiting distinct functional activities and signaling pathways. Specifically, TECs 1 may be involved in tumor cell proliferation and inflammatory responses, whereas TECs 2 is not only involved in cell proliferation pathways, but also enriched in pathways such as metabolic synthesis. Pseudotime analysis revealed dynamic changes in TECs subgroups during HCC progression, correlating specific gene expressions (such as PDGFRB, PGF, JUN, and NR4A1). Subsequently, the JUN gene was predicted by performing binding sites and was shown to act as a transcription factor that may regulate the expression of the PGF gene. CNV analysis highlighted key genes and pathways in TECs that might influence HCC progression, and the PGF as key regulatory factor mediated cell proliferation and migration. Conclusion: The study revealed the heterogeneity of TECs in HCC and their potential roles in tumor progression, offering new perspectives and potential therapeutic targets for HCC molecular mechanisms. The findings emphasize the importance of further exploring TECs heterogeneity for understanding HCC pathogenesis and developing personalized treatment strategies. [ABSTRACT FROM AUTHOR]

Spontaneous forward–reverse mutations were reported by us earlier in clinical samples from various types of cancers and in HeLa cells under normal culture conditions. To investigate the effects of chemical stimulations on such mutation cycles, the present study examined single nucleotide variations (SNVs) and copy number variations (CNVs) in HeLa and A549 cells exposed to wogonin-containing or acidic medium. In wogonin, both cell lines showed a mutation cycle during days 16–18. In acidic medium, both cell lines displayed multiple mutation cycles of different magnitudes. Genomic feature colocalization analysis suggests that CNVs tend to occur in expanded and unstable regions, and near promoters, histones, and non-coding transcription sites. Moreover, phenotypic variations in cell morphology occurred during the forward–reverse mutation cycles under both types of chemical treatments. In conclusion, chemical stresses imposed by wogonin or acidity promoted cyclic forward–reverse mutations in both HeLa and A549 cells to different extents. [ABSTRACT FROM AUTHOR]

Objective: Copy number changes at Chromosomal 16p13.11 have been implicated in a variety of human diseases including congenital cardiac abnormalities. The clinical correlation of copy number variants (CNVs) in this region with developmental abnormalities remains controversial as most of the patients inherit the duplication from an unaffected parent. Methods: We performed CNV analysis on 164 patients with defective left-right (LR) patterning based on whole genome-exome sequencing (WG-ES) followed by multiplex ligation-dependent probe amplification (MLPA) validation. Most cases were accompanied with complex congenital heart disease (CHD). Results: CNVs at 16p13.11 were identified in a total of 21 cases, accounting for 12.80% (21/164) evaluated cases. We observed a marked overrepresentation of chromosome 16p13.11 duplications in cases when compared with healthy controls according to literature reports (15/164, 9.14% versus 0.09% in controls). Notably, in two independent family trios, de novo 16p13.11 microduplications were identified in two patients with laterality defects and CHD. Moreover, 16p13.11 micro-duplication was segregated with the disease in a family trio containing 2 affected individuals. Notably, five coding genes, NOMO1, PKD1P3, NPIPA1, PDXDC1, and NTAN1, were potentially affected by micro-CNV at 16p13.11 in these patients. Conclusion: Our study provides new family-trio based evidences to support 16p13.11 micro-duplications predispose individuals to defective cardiac left-right patterning and laterality disorder. [ABSTRACT FROM AUTHOR]

Background: Tumour dormancy, a resistance mechanism employed by cancer cells, is a significant challenge in cancer treatment, contributing to minimal residual disease (MRD) and potential relapse. Despite its clinical importance, the mechanisms underlying tumour dormancy and MRD remain unclear. In this study, we employed two syngeneic murine models of myeloid leukemia and melanoma to investigate the genetic, epigenetic, transcriptomic and protein signatures associated with tumour dormancy. We used a multiomics approach to elucidate the molecular mechanisms driving MRD and identify potential therapeutic targets. Results: We conducted an in-depth omics analysis encompassing whole-exome sequencing (WES), copy number variation (CNV) analysis, chromatin immunoprecipitation followed by sequencing (ChIP-seq), transcriptome and proteome investigations. WES analysis revealed a modest overlap of gene mutations between melanoma and leukemia dormancy models, with a significant number of mutated genes found exclusively in dormant cells. These exclusive genetic signatures suggest selective pressure during MRD, potentially conferring resistance to the microenvironment or therapies. CNV, histone marks and transcriptomic gene expression signatures combined with Gene Ontology (GO) enrichment analysis highlighted the potential functional roles of the mutated genes, providing insights into the pathways associated with MRD. In addition, we compared "murine MRD genes" profiles to the corresponding human disease through public datasets and highlighted common features according to disease progression. Proteomic analysis combined with multi-omics genetic investigations, revealed a dysregulated proteins signature in dormant cells with minimal genetic mechanism involvement. Pathway enrichment analysis revealed the metabolic, differentiation and cytoskeletal remodeling processes involved in MRD. Finally, we identified 11 common proteins differentially expressed in dormant cells from both pathologies. Conclusions: Our study underscores the complexity of tumour dormancy, implicating both genetic and nongenetic factors. By comparing genomic, transcriptomic, proteomic, and epigenomic datasets, our study provides a comprehensive understanding of the molecular landscape of minimal residual disease. These results provide a robust foundation for forthcoming investigations and offer potential avenues for the advancement of targeted MRD therapies in leukemia and melanoma patients, emphasizing the importance of considering both genetic and nongenetic factors in treatment strategies. [ABSTRACT FROM AUTHOR]

Genetic variation is one of the important factors leading to the incidence and development of cancer. Copy number variation is an important source of genetic diversity, which can be expressed as gene amplification or deletion in structure, and is related to the occurrence and development of different tumors. High-throughput sequencing and gene chip technology can be adopted to detect the variation of copy number, and provide relevant information about tumor molecular characteristics, prognosis and treatment, which is conducive to more accurate diagnosis and treatment decisions for patients in clinical practice. Ovarian cancer is one of the female reproductive system diseases with the highest mortality rate. Understanding its pathogenesis is of significance for improving the survival rate of ovarian cancer. At present, the specific role and mechanism of copy number variation in ovarian cancer are still unclear. In this article, relevant copy number variation in ovarian cancer was reviewed based on the existing research results, aiming to provide novel ideas and methods for the prevention, diagnosis and treatment of ovarian cancer. [ABSTRACT FROM AUTHOR]

Accurate assessment of fragment abundance within a genome is crucial in clinical genomics applications such as the analysis of copy number variation (CNV). However, this task is often hindered by biased coverage in regions with varying guanine–cytosine (GC) content. These biases are particularly exacerbated in hybridization capture sequencing due to GC effects on probe hybridization and polymerase chain reaction (PCR) amplification efficiency. Such GC content–associated variations can exert a negative impact on the fidelity of CNV calling within hybridization capture panels. In this report, we present panelGC, a novel metric, to quantify and monitor GC biases in hybridization capture sequencing data. We establish the efficacy of panelGC, demonstrating its proficiency in identifying and flagging potential procedural anomalies, even in situations where instrument and experimental monitoring data may not be readily accessible. Validation using real-world datasets demonstrates that panelGC enhances the quality control and reliability of hybridization capture panel sequencing. [ABSTRACT FROM AUTHOR]

Purpose: To analyze the copy number variation (CNV) in the X-linked genes BCORL1, POF1B, and USP9X in idiopathic diminished ovarian reserve (DOR). Methods: This case-control study included 47 women, 26 with DOR and 21 in the control group. Age, weight, height, BMI, and FSH level were evaluated, as well as antral follicle count (AFC), oocyte retrieval after controlled ovarian stimulation, and metaphase II (MII) oocytes. The CNVs of BCORL1, USP9X, and POF1B genes were measured by quantitative real time PCR (qPCR) using two reference genes, the HPRT1 (X-linked) and MFN2 (autosomal). Protein–protein interaction network and functional enrichment analysis were performed using the STRING database. Results: The mean age was 36.52 ± 4.75 in DOR women and 35.38 ± 4.14 in control. Anthropometric measures did not differ between the DOR and control groups. DOR women presented higher FSH (p = 0.0025) and lower AFC (p <.0001), oocyte retrieval after COS (p = 0.0004), and MII oocytes (p <.0001) when compared to the control group. BCORL1 and POF1B did not differ in copy number between DOR and control. However, DOR women had more copies of USP9X than the control group (p = 0.028). Conclusion: The increase in the number of copies of the USP9X gene may lead to overexpression in idiopathic DOR and contribute to altered folliculogenesis and oocyte retrieval. [ABSTRACT FROM AUTHOR]

Maturity‐Onset Diabetes of the Young (MODY) is a diabetes mellitus subtype caused by a single gene. The detection rate of the responsible gene is 27% in the United Kingdom, indicating that the causative gene remains unknown in the majority of clinically diagnosed MODY cases. To improve the detection rate, we applied comprehensive genetic testing using whole exome sequencing (WES) followed by Multiplex Ligation‐dependent Probe Amplification (MLPA) and functional analyses. Twenty‐one unrelated Japanese participants with MODY were enrolled in the study. To detect copy number variations (CNVs), WES was performed first, followed by MLPA analysis for participants who were negative on the basis of WES. Undetermined variants were analyzed according to their functional properties. WES identified 7 pathogenic and 3 novel likely pathogenic variants in the 21 participants. Functional analyses revealed that 1 in 3 variants was pathogenic. MLPA analysis applied to the remaining 13 undetermined samples identified 4 cases with pathogenic CNVs: 3 in HNF4A and 1 in HNF1B. Pathogenic variants were identified in 12 participants (12/21, 57.1%) – relatively high rate reported to date. Notably, one‐third of the participants had CNVs in HNF4A or HNF1B, indicating a limitation of WES‐only screening. [ABSTRACT FROM AUTHOR]

The list of standard abbreviations for JDS is available at adsa.org/jds-abbreviations-24. Nonstandard abbreviations are available in the Notes. 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. [ABSTRACT FROM AUTHOR]

Abstract 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 (V ST ) 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.

Abstract Background Coffin-Siris syndrome is a clinically elusive and rare genetic disease characterized by a wide range of clinical manifestations. This study deeply analyzed and identified the clinical phenotype and genetic variant location in a pediatric patient with Coffin-Siris syndrome, aiming to enhance the understanding of this syndrome and assist in its screening and diagnosis. Methods A combination of advanced diagnostic tools, including high-throughput whole-exome sequencing (WES) and first-generation sequencing technologies, was employed to ascertain the etiology of the disease in the child. Results The clinical phenotype was characterized by stunted growth, reduced stature, spina bifida, enuresis, and a ventricular septal defect. WES revealed a de novo variant in the SOX11 gene locus (c.700G > T), identified as pathogenic. It is noteworthy that this variant has not been previously reported. Conclusions The combination of clinical presentation and genetic testing results supports that the patient suffers from Coffin-Siris syndrome due to a genetic variant in the SOX11 gene. This de novo variant expands our understanding of human gene variation, which is conducive to genetic counseling and screening for early diagnosis of Coffin-Siris syndrome.

Abstract The abundance of Lp(a) protein holds significant implications for the risk of cardiovascular disease (CVD), which is directly impacted by the copy number (CN) of KIV-2, a 5.5 kbp sub-region. KIV-2 is highly polymorphic in the population and accurate analysis is challenging. In this study, we present the DRAGEN KIV-2 CN caller, which utilizes short reads. Data across 166 WGS show that the caller has high accuracy, compared to optical mapping and can further phase approximately 50% of the samples. We compared KIV-2 CN numbers to 24 previously postulated KIV-2 relevant SNVs, revealing that many are ineffective predictors of KIV-2 copy number. Population studies, including USA-based cohorts, showed distinct KIV-2 CN, distributions for European-, African-, and Hispanic-American populations and further underscored the limitations of SNV predictors. We demonstrate that the CN estimates correlate significantly with the available Lp(a) protein levels and that phasing is highly important.

Abstract Background Structural variations (SVs) are widespread across genome and have a great impact on evolution, disease, and phenotypic diversity. Despite the development of numerous bioinformatic tools, commonly referred to as SV callers, tailored for detecting SVs using whole genome sequence (WGS) data and employing diverse algorithms, their performance necessitates rigorous evaluation with real data and validated SVs. Moreover, a considerable proportion of these tools have been primarily designed and optimized using human genome data. Consequently, their applicability and performance in Avian species, characterized by smaller genomes and distinct genomic architectures, remain inadequately assessed. Results We performed a comprehensive assessment of the performance of ten widely used SV callers using population-level real genomic data with the validated five common types of SVs. The performance of SV callers varies with the types and sizes of SVs. As compared with other tools, GRIDSS, Lumpy, Wham, and Manta present better detection accuracy. Pindel can detect more small SVs than others. CNVnator and CNVkit can detect more medium and large copy number variations. Given the poor consistency among different SV callers, the combination calling strategy is not recommended. All tools show poor ability in the detection of insertions (especially with size > 150 bp). At least 50× read depth is required to detect more than 80% of the SVs for most tools. Conclusions This study highlights the importance and necessity of using real sequencing data, rather than simulated data only, with validated SVs for SV caller evaluation. Some practical guidance and suggestions are provided for SV detection in future researches.

Liubing Lan,1,2 Dandan Luo,1,2 Jianwen Lian,1 Lingna She,1,3 Bosen Zhang,1,3 Hua Zhong,1 Huaxian Wang,1 Heming Wu1 1Department of Prenatal Diagnostic Center, Meizhou People’s Hospital, Meizhou, People’s Republic of China; 2Department of Obstetrics, Meizhou People’s Hospital, Meizhou, People’s Republic of China; 3Department of Ultrasound, Meizhou People’s Hospital, Meizhou, People’s Republic of ChinaCorrespondence: Heming Wu, Department of Prenatal Diagnostic Center, Meizhou People’s Hospital, Meizhou, People’s Republic of China, Email wuheming1986@126.comObjective: Chromosomal microarray analysis (CMA) is a first-line test to assess the genetic etiology of fetal ultrasound abnormalities. The aim of this study was to evaluate the effectiveness of CMA in detecting chromosomal abnormalities in fetuses with ultrasound abnormalities, including structural abnormalities and non-structural abnormalities.Methods: A retrospective study was conducted on 368 fetuses with abnormal ultrasound who received interventional prenatal diagnosis at Meizhou People’s Hospital from October 2022 to December 2023. Samples of villi, amniotic fluid, and umbilical cord blood were collected according to different gestational weeks, and karyotype and CMA analyses were performed. The detection rate of chromosomal abnormalities in different ultrasonic abnormalities was analyzed.Results: There were 368 fetuses with abnormal ultrasound, including 114 (31.0%) with structural abnormalities, 225 (61.1%) with non-structural abnormalities, and 29 (7.9%) with structural combined with non-structural abnormalities. The detection rate of aneuploidy and pathogenic (P)/likely pathogenic (LP) copy number variations (CNVs) of CMA in fetuses with structural abnormalities was 5.26% (6/114), the detection rate of karyotype was 2.63% (3/114), and the additional diagnosis rate of CMA was 2.63%. In the fetuses with ultrasonic non-structural abnormalities, the detection rate of karyotype was 6.22% (14/225), the detection rate of aneuploidy and P/LP CNVs in fetuses with ultrasonic structural abnormalities was 9.33% (21/225), and the additional diagnosis rate of CMA was 3.11%. There was no significant difference in chromosome abnormality detection rate of CMA among structural abnormality, non-structural abnormality, and structural abnormality combined with non-structural abnormality groups (5.3%, 9.3%, and 13.8%, p = 0.241), also among multiple ultrasonic abnormality and single ultrasonic abnormality groups (14.8%, and 7.3%, p = 0.105).Conclusion: CMA can significantly improve the detection rate of genetic abnormalities in prenatal diagnosis of ultrasonic abnormal fetuses compared with karyotype analysis. CMA is a more effective tool than karyotyping alone in detecting chromosomal abnormalities in fetuses with ultrasound abnormalities.Keywords: chromosomal microarray analysis, copy number variation, abnormal ultrasound fetus, prenatal diagnosis

Abstract Background Repetitive genome regions, such as variable number of tandem repeats (VNTR) or short tandem repeats (STR), are major constituents of the uncharted dark genome and evade conventional sequencing approaches. The protein-coding LPA kringle IV type-2 (KIV-2) VNTR (5.6 kb per unit, 1–40 units per allele) is a medically highly relevant example with a particularly intricate structure, multiple haplotypes, intragenic homologies, and an intra-VNTR STR. It is the primary regulator of plasma lipoprotein(a) [Lp(a)] concentrations, an important cardiovascular risk factor. Lp(a) concentrations vary widely between individuals and ancestries. Multiple variants and functional haplotypes in the LPA gene and especially in the KIV-2 VNTR strongly contribute to this variance. Methods We evaluated the performance of amplicon-based nanopore sequencing with unique molecular identifiers (UMI-ONT-Seq) for SNP detection, haplotype mapping, VNTR unit consensus sequence generation, and copy number estimation via coverage-corrected haplotypes quantification in the KIV-2 VNTR. We used 15 human samples and low-level mixtures (0.5 to 5%) of KIV-2 plasmids as a validation set. We then applied UMI-ONT-Seq to extract KIV-2 VNTR haplotypes in 48 multi-ancestry 1000 Genome samples and analyzed at scale a poorly characterized STR within the KIV-2 VNTR. Results UMI-ONT-Seq detected KIV-2 SNPs down to 1% variant level with high sensitivity, specificity, and precision (0.977 ± 0.018; 1.000 ± 0.0005; 0.993 ± 0.02) and accurately retrieved the full-length haplotype of each VNTR unit. Human variant levels were highly correlated with next-generation sequencing (R 2 = 0.983) without bias across the whole variant level range. Six reads per UMI produced sequences of each KIV-2 unit with Q40 quality. The KIV-2 repeat number determined by coverage-corrected unique haplotype counting was in close agreement with droplet digital PCR (ddPCR), with 70% of the samples falling even within the narrow confidence interval of ddPCR. We then analyzed 62,679 intra-KIV-2 STR sequences and explored KIV-2 SNP haplotype patterns across five ancestries. Conclusions UMI-ONT-Seq accurately retrieves the SNP haplotype and precisely quantifies the VNTR copy number of each repeat unit of the complex KIV-2 VNTR region across multiple ancestries. This study utilizes the KIV-2 VNTR, presenting a novel and potent tool for comprehensive characterization of medically relevant complex genome regions at scale.

Genetic variation is one of the important factors leading to the incidence and development of cancer. Copy number variation is an important source of genetic diversity， which can be expressed as gene amplification or deletion in structure， and is related to the occurrence and development of different tumors. High-throughput sequencing and gene chip technology can be adopted to detect the variation of copy number， and provide relevant information about tumor molecular characteristics， prognosis and treatment， which is conducive to more accurate diagnosis and treatment decisions for patients in clinical practice. Ovarian cancer is one of the female reproductive system diseases with the highest mortality rate. Understanding its pathogenesis is of significance for improving the survival rate of ovarian cancer. At present， the specific role and mechanism of copy number variation in ovarian cancer are still unclear. In this article， relevant copy number variation in ovarian cancer was reviewed based on the existing research results， aiming to provide novel ideas and methods for the prevention， diagnosis and treatment of ovarian cancer.

Purpose of the study. To evaluate the cellular, genomic (gene copy number) and transcriptomic (gene expression) effects of P.hybridus (L.) secondary metabolites when they affect the HeLa cell line.Materials and methods. The isolation of secondary metabolites from plant material and its identification were carried out by preparative chromatography. The composition was determined using mass spectrometric analysis, and the final verification of structural formulas was carried out by nuclear magnetic resonance at the Department of Natural Compounds, the Faculty of Chemistry of the Southern Federal University. The subsequent phase of the study was conducted using both cultural and molecular methods. HeLa cells were cultivated under standard conditions in a MEM medium. Once the confluence level was reached 75–80 %, the nutrient medium was replaced with the introduction of the studied compounds (at a concentration of 4 micrograms/ml) and cultivated for 72 hours. Cell mortality was determined using a NanoEnTek JuliFl counter (Korea) in the presence of 0.4 % trypan blue. The assessment of apoptosis following secondary metabolite exposure was conducted on a BD FACSCanto II flow cytometer using the FITC Annexin V Apoptosis Detection Kit I. The level of replication and expression of the genes responsible for apoptosis was assessed by digital droplet PCR (ddPCR).Results. The following compounds were isolated and verified, and were assigned the following sequence numbers to facilitate their use in the experiment: No. 2 – 2,4-dihydroxy-2,5-dimethylfuran-3(2H)-one, No. 3 – 5-(hydroxymethyl) furan-2-carbaldehyde, No. 5.3 – 2,2,8-trimethyldecahydroazulene-5,6-dicarbaldehyde, P. hybridus (L.) At the stage of cell death assessment, it was found that the greatest effect was achieved in the compound under ordinal No. 2. However, the evaluation of the copy number and expression of the CASP8, CASP9, CASP3, BAX, BCL2, TP53, MDM2, CDKN1B, CDK1, CCND1, CCND3, and RB1 genes by DD-PCR revealed the presence of apoptosis initiation in tumor cells at the molecular level under the action of compounds No. 2 and No. 5.3 obtained from P. hybridus (L.).Conclusion. The outcomes were multifeatured. Only compound 2,4-dihydroxy-2,5-dimethylfuran-3(2H)-one exhibited a pronounced cytostatic effect out of all compounds utilized in the experiment. Concurrently, the compound 2,2,8-trimethyldecahydroazulene-5,6-dicarbaldehyde was found to induce an increase in the expression of the CASP3, CASP8, TP53, and BAX genes.

Abstract Spontaneous forward–reverse mutations were reported by us earlier in clinical samples from various types of cancers and in HeLa cells under normal culture conditions. To investigate the effects of chemical stimulations on such mutation cycles, the present study examined single nucleotide variations (SNVs) and copy number variations (CNVs) in HeLa and A549 cells exposed to wogonin-containing or acidic medium. In wogonin, both cell lines showed a mutation cycle during days 16–18. In acidic medium, both cell lines displayed multiple mutation cycles of different magnitudes. Genomic feature colocalization analysis suggests that CNVs tend to occur in expanded and unstable regions, and near promoters, histones, and non-coding transcription sites. Moreover, phenotypic variations in cell morphology occurred during the forward–reverse mutation cycles under both types of chemical treatments. In conclusion, chemical stresses imposed by wogonin or acidity promoted cyclic forward–reverse mutations in both HeLa and A549 cells to different extents.

Objective: To present the ultrasound imaging and genetic diagnosis of a fetus with prenatal lethal form of Gaucher disease. Case report: A 37-year-old primiparous woman was pregnant at her 23 weeks of gestation and the prenatal fetal ultrasound revealed hydrops fetalis, cerebellum hypoplasia, and fetal immobility. The pregnancy was terminated due to major fetal anomaly, and whole exome sequencing (WES) analysis of fetal tissue and parental blood unveiled a pathogenic variant in exon 10 of the GBA gene (NM_001005741.3: c.1265T > G: p.L422R) originating from the mother. Additionally, a novel CNV (chr1: 155204785–155205635 deletion, 0.85 kb) spanning exon 10–12 in the GBA gene was identified from the father. This compound heterozygosity confirmed the diagnosis of prenatal lethal form of Gaucher disease and was informative for genetic counseling. Conclusion: WES is a powerful tool to detect pathogenic variants among fetuses with nonimmune hydrops fetalis and complex abnormality from prenatal ultrasound. Compound heterozygosity consisted of single nucleotide variants (SNV) and copy number variations (CNVs) may lead rare inherited metabolic disorders including prenatal lethal form of Gaucher disease.

ABSTRACT: 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.

PURPOSE OF REVIEW: To summarize current literature available on sleep in 22q11.2 Deletion Syndrome (22q11.2DS; Velocardiofacial or DiGeorge Syndrome), a neurogenetic disorder caused by a hemizygous deletion in a genomic region critical for neurodevelopment. Due to the greatly increased risk of developmental psychiatric disorders (e.g., autism and schizophrenia) in 22q11.2DS, this review focuses on clinical correlates of sleep disturbances and potential neurobiological underpinnings of these relationships. RECENT FINDINGS: Sleep disturbances are widely prevalent in 22q11.2DS and are associated with worse behavioral, psychiatric, and physical health outcomes. There are reports of sleep architecture and sleep neurophysiology differences, but the literature is limited by logistical challenges posed by objective sleep measures, resulting in small study samples to date. Sleep disturbances in 22q11.2DS are prevalent and have a substantial impact on well-being. Further investigation of sleep in 22q11.2DS utilizing multimodal sleep assessments has the potential to provide new insight into neurobiological mechanisms and a potential trans-diagnostic treatment target in 22q11.2DS.

Copy number variation (CNV) is an essential genetic driving factor of cancer formation and progression, making intelligent classification based on CNV feasible. However, there are a few challenges in the current machine learning and deep learning methods, such as the design of base classifier combination schemes in ensemble methods and the selection of layers of neural networks, which often result in low accuracy. Therefore, an adaptive bilinear dynamic cascade model (Adap-BDCM) is developed to further enhance the accuracy and applicability of these methods for intelligent classification on CNV datasets. In this model, a feature selection module is introduced to mitigate the interference of redundant information, and a bilinear model based on the gated attention mechanism is proposed to extract more beneficial deep fusion features. Furthermore, an adaptive base classifier selection scheme is designed to overcome the difficulty of manually designing base classifier combinations and enhance the applicability of the model. Lastly, a novel feature fusion scheme with an attribute recall submodule is constructed, effectively avoiding getting stuck in local solutions and missing some valuable information. Numerous experiments have demonstrated that our Adap-BDCM model exhibits optimal performance in cancer classification, stage prediction, and recurrence on CNV datasets. This study can assist physicians in making diagnoses faster and better. [ABSTRACT FROM AUTHOR]

Abstract Copy number variation (CNV) tends to occur in genetically enriched regions and is likely associated with a number of complex diseases such as skin aging. In this study, we investigated the genome-wide CNVs in 20 wrinkled skin cases (WSC) of Xiang pigs and 63 controls, and identified 7893 copy number variable regions (CNVRs). We estimated the F-statistic (Fst) at each locus and identified that 93 case-controls stratified CNVRs (Fst > = 0.15) overlapped with 87 known genes. Functional enrichment analysis showed that most of these genes were predominantly enriched in pathways and terms related to the extracellular matrix. Finally, we found that some CNVs were predicted to have high effects on genes such as VCAN, TIMP1 and FOXO1 through transcriptional amplification, transcript ablation and so on. Most of the genes overlapped with those CNVRs have been reported to be related to aging in human or animals. The copy numbers presented the positive correlations with the transcript level of the genes in skins between the cases and controls. Our results suggested that those 22 CNVRs, including 19 CNV losses and 3 CNV gains, were putatively associated with the skin wrinkle of Xiang pigs.

Xiaoqin Chen,1,2 Liubing Lan,1,2 Heming Wu,1 Mei Zeng,1,2 Zhiyuan Zheng,1 Qiuping Zhong,2 Fengdan Lai,2 Yonghe Hu2 1Department of Prenatal Diagnostic Center, Meizhou People’s Hospital, Meizhou Academy of Medical Sciences, Meizhou, People’s Republic of China; 2Department of Obstetrics, Meizhou People’s Hospital, Meizhou Academy of Medical Sciences, Meizhou, People’s Republic of ChinaCorrespondence: Liubing Lan, Department of Prenatal Diagnostic Center, Meizhou People’s Hospital, Meizhou Academy of Medical Sciences, No. 63 Huangtang Road, Meijiang District, Meizhou, People’s Republic of China, Email 13411221666@163.comObjective: To explore and evaluate the value of chromosomal microarray analysis (CMA) in prenatal diagnosis of fetuses with ultrasound abnormalities.Methods: A retrospective analysis was performed on 370 fetuses with ultrasound abnormalities received invasive prenatal diagnosis at Meizhou People’s Hospital from October 2022 to December 2023. Fetal specimens were analyzed by CMA, and the detection rates of aneuploidy and pathogenic (P)/likely pathogenic (LP) copy number variations (CNVs) in ultrasound structural abnormalities (malformations of fetal anatomy) and non-structural abnormalities (abnormalities of fetal nonanatomical structure) were analyzed.Results: There were 114 (30.8%) cases with isolated ultrasound structural abnormalities, 226 (61.1%) cases with isolated non-structural abnormalities (182 isolated ultrasound soft markers abnormalities, 30 isolated fetal growth restriction (FGR), and 8 isolated abnormalities of amniotic fluid volume), and 30 (8.1%) cases with both structural and non-structural abnormalities. The overall detection rate of aneuploidy and P/LP CNVs in isolated ultrasonic structural abnormalities was 5.3%, among which cardiovascular system abnormalities were the highest. In addition, the largest number of fetuses with non-structural abnormalities was nuchal translucency (NT) thickening (n = 81), followed by ventriculomegaly (n = 29), and nasal bone dysplasia (n = 24). The detection rate of chromosomal abnormalities of fetuses with abnormal ultrasound soft markers was 9.9%, and the detection rate in single abnormal ultrasound soft marker, and multiple ultrasound soft markers abnormalities was 9.7% (16/165) and 11.8% (2/17), respectively. Moreover, the detection rate of chromosomal abnormalities of fetuses with FGR and structural abnormalities combined with non-structural abnormalities was 6.7% (2/30), and 13.3% (4/30), respectively.Conclusion: The incidence of chromosomal abnormalities (aneuploidy and P/LP CNVs) varies among different fetal ultrasound abnormalities.Keywords: chromosomal microarray analysis, copy number variation, abnormal ultrasound fetus, prenatal diagnosis

Abstract Background Giardiasis, caused by the protozoan parasite Giardia intestinalis, often presents a treatment challenge, particularly in terms of resistance to metronidazole. Despite extensive research, markers for metronidazole resistance have not yet been identified. Methods This study analysed 28 clinical samples of G. intestinalis from sub-assemblage AII, characterised by varying responses to metronidazole treatment. We focussed on copy number variation (CNV) of the multi-copy flavohemoprotein gene, analysed using digital polymerase chain reaction (dPCR) and next generation sequencing (NGS). Additionally, chromosomal ploidy was tested in 18 of these samples. Flavohemoprotein CNV was also assessed in 17 samples from other sub-assemblages. Results Analyses revealed variable CNVs of the flavohemoprotein gene among the isolates, with no correlation to clinical metronidazole resistance. Discrepancies in CNVs detected from NGS data were attributed to biases linked to the whole genome amplification. However, dPCR helped to clarify these discrepancies by providing more consistent CNV data. Significant differences in flavohemoprotein CNVs were observed across different G. intestinalis sub-assemblages. Notably, Giardia exhibits a propensity for aneuploidy, contributing to genomic variability within and between sub-assemblages. Conclusions The complexity of the clinical metronidazole resistance in Giardia is influenced by multiple genetic factors, including CNVs and aneuploidy. No significant differences in the CNV of the flavohemoprotein gene between isolates from metronidazole-resistant and metronidazole-sensitive cases of giardiasis were found, underscoring the need for further research to identify reliable genetic markers for resistance. We demonstrate that dPCR and NGS are robust methods for analysing CNVs and provide cross-validating results, highlighting their utility in the genetic analyses of this parasite. Graphical Abstract

Abstract Background Teat number is a vital reproductive trait in sows, crucial for providing immunity and nutrition to piglets during lactation. However, “missing heritability” in Single Nucleotide Polymorphism (SNP)-based Genome-Wide Association Studies (GWAS) has led to an increasing focus on structural variations in the genetic analysis of complex biological traits. Results In this study, we generated a comprehensive CNV map in a population of French Yorkshire pigs (n = 644) and identified 429 CNVRs. Notably, 44% (189 CNVRs) of these were detected for the first time. Subsequently, we conducted GWAS for teat number in the French Yorkshire pig population using both 80K chip and its imputed data, as well as a GWAS analysis based on CNV regions (CNVR). Interestingly, 80K chip GWAS identified two SNPs located on Sus scrofa chromosome 5 (SSC5) that were simultaneously associated with Total Teat Number (TTN), Left Teat Number (LTN), and Right Teat Number (RTN). The leading SNP (WU_10.2_5_76130558) explained 3.33%, 2.69%, and 2.67% of the phenotypic variance for TTN, LTN, and RTN, respectively. Moreover, through imputed GWAS, we successfully identified 30 genetic variants associated with TTN located within the 73.22 -73.30 Mb region on SSC5. The two SNPs identified in the 80K chip GWAS were also located in this region. In addition, CNVR-based GWAS revealed three significant CNVRs associated with TTN. Finally, through gene annotation, we pinpointed two candidate genes, TRIM66 and PRICKLE1, which are related to diverse processes such as breast cancer and abnormal vertebral development. Conclusions Our research provides an in-depth analysis of the complex genetic structure underlying teat number, contributing to the genetic enhancement of sows with improved reproductive performance and, ultimately, bolstering the economic benefits of swine production enterprises.

Abstract Background Structural variations play an important role in bacterial genomes. They can mediate genome adaptation quickly in response to the external environment and thus can also play a role in antibiotic resistance. The detection of structural variations in bacteria is challenging, and the recognition of even small rearrangements can be important. Even though most detection tools are aimed at and benchmarked on eukaryotic genomes, they can also be used on prokaryotic genomes. The key features of detection are the ability to detect small rearrangements and support haploid genomes. Because of the limiting performance of a single detection tool, combining the detection abilities of multiple tools can lead to more robust results. There are already available workflows for structural variation detection for long-reads technologies and for the detection of single-nucleotide variation and indels, both aimed at bacteria. Yet we are unaware of structural variations detection workflows for the short-reads sequencing platform. Motivated by this gap we created our workflow. Further, we were interested in increasing the detection performance and providing more robust results. Results We developed an open-source bioinformatics pipeline, ProcaryaSV, for the detection of structural variations in bacterial isolates from paired-end short sequencing reads. Multiple tools, starting with quality control and trimming of sequencing data, alignment to the reference genome, and multiple structural variation detection tools, are integrated. All the partial results are then processed and merged with an in-house merging algorithm. Compared with a single detection approach, ProcaryaSV has improved detection performance and is a reproducible easy-to-use tool. Conclusions The ProcaryaSV pipeline provides an integrative approach to structural variation detection from paired-end next-generation sequencing of bacterial samples. It can be easily installed and used on Linux machines. It is publicly available on GitHub at https://github.com/robinjugas/ProcaryaSV .

Abstract Background The massive structural variations and frequent introgression highly contribute to the genetic diversity of wheat, while the huge and complex genome of polyploid wheat hinders efficient genotyping of abundant varieties towards accurate identification, management, and exploitation of germplasm resources. Results We develop a novel workflow that identifies 1240 high-quality large copy number variation blocks (CNVb) in wheat at the pan-genome level, demonstrating that CNVb can serve as an ideal DNA fingerprinting marker for discriminating massive varieties, with the accuracy validated by PCR assay. We then construct a digitalized genotyping CNVb map across 1599 global wheat accessions. Key CNVb markers are linked with trait-associated introgressions, such as the 1RS·1BL translocation and 2NvS translocation, and the beneficial alleles, such as the end-use quality allele Glu-D1d (Dx5 + Dy10) and the semi-dwarf r-e-z allele. Furthermore, we demonstrate that these tagged CNVb markers promote a stable and cost-effective strategy for evaluating wheat germplasm resources with ultra-low-coverage sequencing data, competing with SNP array for applications such as evaluating new varieties, efficient management of collections in gene banks, and describing wheat germplasm resources in a digitalized manner. We also develop a user-friendly interactive platform, WheatCNVb ( http://wheat.cau.edu.cn/WheatCNVb/ ), for exploring the CNVb profiles over ever-increasing wheat accessions, and also propose a QR-code-like representation of individual digital CNVb fingerprint. This platform also allows uploading new CNVb profiles for comparison with stored varieties. Conclusions The CNVb-based approach provides a low-cost and high-throughput genotyping strategy for enabling digitalized wheat germplasm management and modern breeding with precise and practical decision-making.

The genetic mechanism of congenital heart disease (CHD) is complex and currently lacks a clear understanding. Literature studies on CHD often report the presence of concurrent extracardiac anomalies, but since the majority of CHD cases are isolated, presenting only a single cardiac malformation, the etiological mechanisms remain uncertain, especially regarding the genetic aspects. Furthermore, there is a scarcity of case studies focusing on isolated CHD, resulting in a lack of comprehensive research data. Therefore, elucidating the genetic causes of isolated CHD and providing guidance for its clinical treatment holds significant research value in the field of CHD. This article reviews the known genetic causes and potential genetic mechanisms of isolated CHD, as well as provides recommendations for genetic testing in patients with isolated CHD.

Abstract Background Copy number variation (CNV) is a key genetic characteristic for cancer diagnostics and can be used as a biomarker for the selection of therapeutic treatments. Using data sets established in our previous study, we benchmark the performance of cancer CNV calling by six most recent and commonly used software tools on their detection accuracy, sensitivity, and reproducibility. In comparison to other orthogonal methods, such as microarray and Bionano, we also explore the consistency of CNV calling across different technologies on a challenging genome. Results While consistent results are observed for copy gain, loss, and loss of heterozygosity (LOH) calls across sequencing centers, CNV callers, and different technologies, variation of CNV calls are mostly affected by the determination of genome ploidy. Using consensus results from six CNV callers and confirmation from three orthogonal methods, we establish a high confident CNV call set for the reference cancer cell line (HCC1395). Conclusions NGS technologies and current bioinformatics tools can offer reliable results for detection of copy gain, loss, and LOH. However, when working with a hyper-diploid genome, some software tools can call excessive copy gain or loss due to inaccurate assessment of genome ploidy. With performance matrices on various experimental conditions, this study raises awareness within the cancer research community for the selection of sequencing platforms, sample preparation, sequencing coverage, and the choice of CNV detection tools.

Background: CYP2D6 testing is increasingly used to guide drug therapy and thus, reliable methods are needed to test this complex and polymorphic gene locus. A particular challenge arises from the detection and interpretation of structural variants (SVs) including gene deletions, duplications, and hybrids with the CYP2D7 pseudogene. This study validated the Absolute Q™ platform for digital PCR-based CYP2D6 copy number variation (CNV) determination by comparing results to those obtained with a previously established method using the QX200 platform. In addition, protocols for streamlining CYP2D6 CNV testing were established and validated including the "One-pot" single-step restriction enzyme digestion and a multiplex assay simultaneously targeting the CYP2D6 5'UTR, intron 6, and exon 9 regions. Methods: Genomic DNA (gDNA) samples from Coriell (n = 13) and from blood, saliva, and liver tissue (n = 17) representing 0-6 copies were tested on the Absolute Q and QX200 platforms. Custom TaqMan™ copy number (CN) assays targeting CYP2D6 the 5'UTR, intron 6, and exon 9 regions and a reference gene assay (TERT or RNaseP) were combined for multiplexing by optical channel. In addition, two digestion methods (One-pot digestion and traditional) were assessed. Inconclusive CN values on the Absolute Q were resolved using an alternate reference gene and/or diluting gDNA. Results: Overall, results between the two platforms and digestions methods were consistent. The "One-pot" digestion method and optically multiplexing up to three CYP2D6 regions yielded consistent result across DNA sample types and diverse SVs, reliably detecting up to 6 gene copies. Rare variation in reference genes were found to interfere with results and interpretation, which were resolved by using a different reference. Conclusion: The Absolute Q produced accurate and reliable CYP2D6 copy number results allowing for a streamlined and economical protocol using One-pot digestion and multiplexing three target regions. Protocols are currently being expanded to other pharmacogenes presenting with SVs/CNVs. [ABSTRACT FROM AUTHOR]

Copy number variation (CNV) tends to occur in genetically enriched regions and is likely associated with a number of complex diseases such as skin aging. In this study, we investigated the genome-wide CNVs in 20 wrinkled skin cases (WSC) of Xiang pigs and 63 controls, and identified 7893 copy number variable regions (CNVRs). We estimated the F-statistic (Fst) at each locus and identified that 93 case-controls stratified CNVRs (Fst > = 0.15) overlapped with 87 known genes. Functional enrichment analysis showed that most of these genes were predominantly enriched in pathways and terms related to the extracellular matrix. Finally, we found that some CNVs were predicted to have high effects on genes such as VCAN, TIMP1 and FOXO1 through transcriptional amplification, transcript ablation and so on. Most of the genes overlapped with those CNVRs have been reported to be related to aging in human or animals. The copy numbers presented the positive correlations with the transcript level of the genes in skins between the cases and controls. Our results suggested that those 22 CNVRs, including 19 CNV losses and 3 CNV gains, were putatively associated with the skin wrinkle of Xiang pigs. [ABSTRACT FROM AUTHOR]

Genetic variation in the FCGR3B gene is responsible for different variants of human neutrophil antigen 1 (HNA-1). Laboratory techniques currently utilized for routine HNA-1 genotyping, predominantly PCR-sequence-specific primer (PCR-SSP) and PCR-sequence-based typing (PCR-SBT), lack specificity for FCGR3B. This study compares the capabilities and limitations of existing technologies including an in-house TaqMan PCR, a commercial PCR-SSP test, PCR-SBT and multiplex ligation-dependent probe amplification (MLPA) with those of a long-read nanopore sequencing assay. Testing was performed with both related and unrelated Danish samples with different copy numbers and/or rare alleles. Long-read nanopore sequencing was validated by blind testing of ten English samples. The results showed that FCGR3B copy numbers correlate with a dose-dependent distribution of alleles that complicates genotyping by TaqMan PCR, PCR-SSP and PCR-SBT, due to co-amplification of the homologous FCGR3A gene. MLPA can correctly quantify the dose-dependent distribution but not detect novel variants. Long-read nanopore sequencing showed high specificity for FCGR3B and was able to detect dosage-dependent distribution, and rare and novel variants that were previously not described. Current HNA-1 genotyping methods cannot produce unambiguous allele-level results, whereas long-read nanopore sequencing has shown the potential to resolve observed ambiguities, identify new HNA-1 variants and allow definitive allele assignment. [ABSTRACT FROM AUTHOR]

Human pluripotent stem cells (hPSCs) are pivotal in regenerative medicine, yet their in vitro expansion often leads to genetic abnormalities, raising concerns about their safety in clinical applications. This study analyzed ten human embryonic stem cell lines across multiple passages to elucidate the dynamics of chromosomal abnormalities and single-nucleotide variants (SNVs) in 380 cancer-related genes. Prolonged in vitro culture resulted in 80% of the lines acquiring gains of chromosome 20q or 1q, both known for conferring an in vitro growth advantage. 70% of lines also acquired other copy number variants (CNVs) outside the recurrent set. Additionally, we detected 122 SNVs in 88 genes, with all lines acquiring at least one de novo SNV during culture. Our findings showed higher loads of both CNVs and SNVs at later passages, which were due to the cumulative acquisition of mutations over a longer time in culture, and not to an increased rate of mutagenesis over time. Importantly, we observed that SNVs and rare CNVs followed the acquisition of chromosomal gains in 1q and 20q, while most of the low-passage and genetically balanced samples were devoid of cancer-associated mutations. This suggests that recurrent chromosomal abnormalities are potential drivers for the acquisition of other mutations. [ABSTRACT FROM AUTHOR]

Background: Giardiasis, caused by the protozoan parasite Giardia intestinalis, often presents a treatment challenge, particularly in terms of resistance to metronidazole. Despite extensive research, markers for metronidazole resistance have not yet been identified. Methods: This study analysed 28 clinical samples of G. intestinalis from sub-assemblage AII, characterised by varying responses to metronidazole treatment. We focussed on copy number variation (CNV) of the multi-copy flavohemoprotein gene, analysed using digital polymerase chain reaction (dPCR) and next generation sequencing (NGS). Additionally, chromosomal ploidy was tested in 18 of these samples. Flavohemoprotein CNV was also assessed in 17 samples from other sub-assemblages. Results: Analyses revealed variable CNVs of the flavohemoprotein gene among the isolates, with no correlation to clinical metronidazole resistance. Discrepancies in CNVs detected from NGS data were attributed to biases linked to the whole genome amplification. However, dPCR helped to clarify these discrepancies by providing more consistent CNV data. Significant differences in flavohemoprotein CNVs were observed across different G. intestinalis sub-assemblages. Notably, Giardia exhibits a propensity for aneuploidy, contributing to genomic variability within and between sub-assemblages. Conclusions: The complexity of the clinical metronidazole resistance in Giardia is influenced by multiple genetic factors, including CNVs and aneuploidy. No significant differences in the CNV of the flavohemoprotein gene between isolates from metronidazole-resistant and metronidazole-sensitive cases of giardiasis were found, underscoring the need for further research to identify reliable genetic markers for resistance. We demonstrate that dPCR and NGS are robust methods for analysing CNVs and provide cross-validating results, highlighting their utility in the genetic analyses of this parasite. [ABSTRACT FROM AUTHOR]

We summarize the copy number variations (CNVs) and phenotype spectrum of infantile epileptic spasms syndrome (IESS) in a Chinese cohort. The CNVs were identified by genomic copy number variation sequencing. The CNVs and clinical data were analyzed. 74 IESS children with CNVs were enrolled. 35 kinds of CNVs were identified. There were 11 deletions and 5 duplications not reported previously in IESS, including 2 CNVs not reported in epilepsy. 87.8% were de novo, 9.5% were inherited from mother and 2.7% from father. Mosaicism occurred in one patient with Xq21.31q25 duplication. 16.2% (12/74) were 1p36 deletion, and 20.3% (15/74) were 15q11‐q13 duplication. The age of seizure onset ranged from 17 days to 24 months. Seizure types included epileptic spasms, focal seizures, tonic seizures, and myoclonic seizures. All patients displayed developmental delay. Additional features included craniofacial anomaly, microcephaly, congenital heart defects, and hemangioma. 29.7% of patients were seizure‐free for more than 12 months, and 70.3% still had seizures after trying 2 or more anti‐seizure medications. In conclusion, CNVs is a prominent etiology of IESS. 1p36 deletion and 15q duplication occurred most frequently. CNV detection should be performed in patients with IESS of unknown causes, especially in children with craniofacial anomalies and microcephaly. [ABSTRACT FROM AUTHOR]

Autistic spectrum disorder (ASD) is a neurodevelopmental disability characterised by the impairment of social interaction and communication ability. The alarming increase in its prevalence in children urged researchers to obtain a better understanding of the causes of this disease. Genetic factors are considered to be crucial, as ASD has a tendency to run in families. In recent years, with technological advances, the importance of structural variations (SVs) in ASD began to emerge. Most of these studies, however, focus on the Caucasian population. As a populated ethnicity, ASD shall be a significant health issue in China. This systematic review aims to summarise current case-control studies of SVs associated with ASD in the Chinese population. A list of genes identified in the nine included studies is provided. It also reveals that similar research focusing on other genetic backgrounds is demanded to manifest the disease etiology in different ethnic groups, and assist the development of accurate ethnic-oriented genetic diagnosis. [ABSTRACT FROM AUTHOR]

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 thicknessQTL, Longissimusmuscle areaQTL, andMarbling scoreQTL. 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. [ABSTRACT FROM AUTHOR]

Introduction: Copy number variation (CNV) consist of deletion, insertion, and duplications. It is an important source of genetic variation in organisms and thus influences on the gene expression and phenotypic variation. Copy number variation (CNV) is one of the structural variant with an intermediate size class larger than 50bp which involves unbalanced rearrangements that increase or decrease the amount of DNA (Pirooznia et al 2015, Alkan et al 2011). The size of CNVs is larger than 50bp, while smaller segments are known as insertions or deletions (indels). Thereupon these structural variations comprise more polymorphic than SNPs because of enormity, detection of them and their effect on phenotype has caught the attention of many researchers recently. It has been reported that CNVs changes in gene dosage and regulation as well as in transcript structure, and thus contribute to phenotypic variability (Pirooznia et al 2015, Alkan et al 2011). The pea-comb phenotype is caused by a CNV mapping to intron 1 of the SRY (sex determining region Y)-box 5 (SOX5) gene (Wright et al 2009). Late feathering in chickens is due to incomplete duplication in PRLR and SPEF2 genes (Elfrink et al 2008). In swine, dominant white colour has been related with a duplication of a 450-kb fragment of the KIT gene (Giuffra et al 1999) and a splice mutation causing the skipping of exon 17 (Giuffra et al 1999). In sheep, doubling in the ASIP gene results in the regulation of pigment in body coat (Norris et al 2008). Doubling the 4.6 k base pair into the six introns of the STX17 gene results in a gray body color in the horse with age. Deletion of the intergenic region with a length of 11.7 kbp in the goat genome leads to the removal of horns (Clop et al 2012). Chicken is the most intensively farmed animal on earth and is a major food source with billions of birds used in meat and egg production each year. A big share of chicken CNVs involves protein coding or regulatory sequences. A comprehensive study of chicken CNV can provide valuable information on genetic diversity and assist future analyses of associations between CNV and economically important traits in chickens. Unique chicken genome with macro and micro chromosomes and its biology make it an ideal organism for studies in development and evolution, as well as applications in agriculture and medicine (Burt 2005). In the last several years, there has been an increasing interest in the study of CNVs in the chicken. This study focuses on comparison of CNV between the broilers and layers chicken to find evidence of domestication on the genome using whole genome sequencing Material and methods: we used n=90 female birds of two commercial broiler (n=40) and layer (n=50) chicken. The broilers (BRs) were represented by 20 DNA samples of each of two lines (BRA and BRB) established independently and previously collected as part of the AVIANDIV project. In the layer group (LRs), data from 25 birds each from purebred white (WL) and brown (BL) egg laying populations, sequenced in the frame of the SYNBREED project (http://www.synbreed.tum.de/index.php?id=2(, were included. The paired-end reads with a read length of 101bp were mapped against the current reference genome assembly Galgal6 using the Burrows-Wheeler aligner (bwa, 0.6.2-r126 Version, with default parameters. Duplicate reads were masked during post-processing using the Picard tool set (version 2.9.2, http://picard.sourceforge.net). Finally, Genome Analysis Toolkit-3.3.0 was used to realign reads for correcting errors caused by indels. Using GATK software package and Depth of Coverage function (McKenna et al 2010), the depth of readings was calculated for each sample. Then filter out reads with mapping quality below 20. Because comparing the genomes of individuals in different groups was time consuming and computationally difficult for all parts of the genome, the genomes of each individual were divided into 1000 bp non-overlapping windows and the average reading depth per window was calculated. Then the results were normalized against the BL sample that showed highest average depth. In short, we created a correction factor per population and applied it on the depth of coverage value for each window. For all the contrasts, we performed an analysis of variance (ANOVA) as described (Carneiro et al 2014). For the Broilers-Layers contrast we scanned 935247 windows. 69522windows showed significant by FDR with P < 0.001, with ANOVA using the Benjamini-Hochberg FDR method for multiple corrections (Benjamini and Hochberg 1995). Results and discussion: Mapping sequencing data to galGal6 assembly showed an average 98.61% mapping rate and 11.51 depth. Manhattan plot was plotted for regions of the genome that differed significantly between the two groups (FDR = 0.001). The points above the hypothetical line were identified and examined in a 25 Kbp confidence interval to identify possible genes. 39 regions were identified that half of them does not contain any genes. Although Long noncoding RNAs are under lower selective pressure than protein-coding genes (Batista and Chang 2013), the other 11 regions contained 16 genes related to long non-coding RNAs. Long noncoding RNAs (lncRNAs) play a critical role in organizing the 3-dimensional genome architecture and regulating gene activity in cis or in trans through multiple mechanisms (Zhang et al. 2019, Batista and Chang 2013). 6 Othere regions also contained 12 coding genes. Most of the identified genes were somehow linked to the immune system disease or cancer. Genes such as DEDs and TNFAIP8 are involved in programmed cell death (apoptosis) and two genes NPAL3 and RCAN, which are involved in the immune system, had a copy number variation in the studied samples. In addition, RCAN is involved in Down syndrome. The PFDN gene, located on chromosome 25, is also involved in Alzheimer's and Parkinson's disease. [ABSTRACT FROM AUTHOR]

We developed a generally applicable method, CRISPR/Cas9-targeted long-read sequencing (CTLR-Seq), to resolve, haplotype-specifically, the large and complex regions in the human genome that had been previously impenetrable to sequencing analysis, such as large segmental duplications (SegDups) and their associated genome rearrangements. CTLR-Seq combines in vitro Cas9-mediated cutting of the genome and pulse-field gel electrophoresis to isolate intact large (i.e., up to 2,000 kb) genomic regions that encompass previously unresolvable genomic sequences. These targets are then sequenced (amplification-free) at high on-target coverage using long-read sequencing, allowing for their complete sequence assembly. We applied CTLR-Seq to the SegDup-mediated rearrangements that constitute the boundaries of, and give rise to, the 22q11.2 Deletion Syndrome (22q11DS), the most common human microdeletion disorder. We then performed de novo assembly to resolve, at base-pair resolution, the full sequence rearrangements and exact chromosomal breakpoints of 22q11.2DS (including all common subtypes). Across multiple patients, we found a high degree of variability for both the rearranged SegDup sequences and the exact chromosomal breakpoint locations, which coincide with various transposons within the 22q11.2 SegDups, suggesting that 22q11DS can be driven by transposon-mediated genome recombination. Guided by CTLR-Seq results from two 22q11DS patients, we performed three-dimensional chromosomal folding analysis for the 22q11.2 SegDups from patient-derived neurons and astrocytes and found chromosome interactions anchored within the SegDups to be both cell type-specific and patient-specific. Lastly, we demonstrated that CTLR-Seq enables cell-type specific analysis of DNA methylation patterns within the deletion haplotype of 22q11DS. [ABSTRACT FROM AUTHOR]

Background: Teat number is a vital reproductive trait in sows, crucial for providing immunity and nutrition to piglets during lactation. However, "missing heritability" in Single Nucleotide Polymorphism (SNP)-based Genome-Wide Association Studies (GWAS) has led to an increasing focus on structural variations in the genetic analysis of complex biological traits. Results: In this study, we generated a comprehensive CNV map in a population of French Yorkshire pigs (n = 644) and identified 429 CNVRs. Notably, 44% (189 CNVRs) of these were detected for the first time. Subsequently, we conducted GWAS for teat number in the French Yorkshire pig population using both 80K chip and its imputed data, as well as a GWAS analysis based on CNV regions (CNVR). Interestingly, 80K chip GWAS identified two SNPs located on Sus scrofa chromosome 5 (SSC5) that were simultaneously associated with Total Teat Number (TTN), Left Teat Number (LTN), and Right Teat Number (RTN). The leading SNP (WU_10.2_5_76130558) explained 3.33%, 2.69%, and 2.67% of the phenotypic variance for TTN, LTN, and RTN, respectively. Moreover, through imputed GWAS, we successfully identified 30 genetic variants associated with TTN located within the 73.22 -73.30 Mb region on SSC5. The two SNPs identified in the 80K chip GWAS were also located in this region. In addition, CNVR-based GWAS revealed three significant CNVRs associated with TTN. Finally, through gene annotation, we pinpointed two candidate genes, TRIM66 and PRICKLE1, which are related to diverse processes such as breast cancer and abnormal vertebral development. Conclusions: Our research provides an in-depth analysis of the complex genetic structure underlying teat number, contributing to the genetic enhancement of sows with improved reproductive performance and, ultimately, bolstering the economic benefits of swine production enterprises. [ABSTRACT FROM AUTHOR]