ProcaryaSV: structural variation detection pipeline for bacterial genomes using short-read sequencing.

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 .

Tagging large CNV blocks in wheat boosts digitalization of germplasm resources by ultra-low-coverage sequencing.

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.

CDKN2A copy number alteration in bladder cancer: Integrative analysis in patient-derived xenografts and cancer patients.

Bladder cancer (BlCa) is an extensively heterogeneous disease that leads to great variability in tumor evolution scenarios and lifelong patient surveillance, emphasizing the need for modern, minimally invasive precision medicine. Here, we explored the clinical significance of copy number alterations (CNAs) in BlCa. CNA profiling was performed in 15 patient-derived xenografts (PDXs) and validated in The Cancer Genome Atlas BlCa (TCGA-BLCA; n = 408) and Lindgren et al. (n = 143) cohorts. CDKN2A copy number loss was identified as the most frequent CNA in bladder tumors, associated with reduced CDKN2A expression, tumors of a papillary phenotype, and prolonged PDX survival. The study's screening cohort consisted of 243 BlCa patients, and CDKN2A copy number was assessed in genomic DNA and cell-free DNA (cfDNA) from 217 tumors and 189 pre-treatment serum samples, respectively. CDKN2A copy number loss was correlated with superior disease-free and progression-free survival of non-muscle-invasive BlCa (NMIBC) patients. Moreover, a higher CDKN2A index (CDKN2A/LEP ratio) in pre-treatment cfDNA was associated with advanced tumor stage and grade and short-term NMIBC progression to invasive disease, while multivariate models fitted for CDKN2A index in pre-treatment cfDNA offered superior risk stratification of T1/high-grade and EORTC high-risk patients, enhancing prediction of treatment outcome. CDKN2A copy number status could serve as a minimally invasive tool to improve risk stratification and support personalized prognosis in BlCa.

Unveiling novel genetic variants in 370 challenging medically relevant genes using the long read sequencing data of 41 samples from 19 global populations.

BACKGROUND: A large number of challenging medically relevant genes (CMRGs) are situated in complex or highly repetitive regions of the human genome, hindering comprehensive characterization of genetic variants using next-generation sequencing technologies. In this study, we employed long-read sequencing technology, extensively utilized in studying complex genomic regions, to characterize genetic alterations, including short variants (single nucleotide variants and short insertions and deletions) and copy number variations, in 370 CMRGs across 41 individuals from 19 global populations. RESULTS: Our analysis revealed high levels of genetic variants in CMRGs, with 68.73% exhibiting copy number variations and 65.20% containing short variants that may disrupt protein function across individuals. Such variants can influence pharmacogenomics, genetic disease susceptibility, and other clinical outcomes. We observed significant differences in CMRG variation across populations, with individuals of African ancestry harboring the highest number of copy number variants and short variants compared to samples from other continents. Notably, 15.79% to 33.96% of short variants were exclusively detectable through long-read sequencing. While the T2T-CHM13 reference genome significantly improved the assembly of CMRG regions, thereby facilitating variant detection in these regions, some regions still lacked resolution. CONCLUSION: Our results provide an important reference for future clinical and pharmacogenetic studies, highlighting the need for a comprehensive representation of global genetic diversity in the reference genome and improved variant calling techniques to fully resolve medically relevant genes.

A novel copy number variant in the murine Cdh23 gene gives rise to profound deafness and vestibular dysfunction.

Hearing loss is the most common congenital sensory deficit worldwide and exhibits high genetic heterogeneity, making molecular diagnoses elusive for most individuals. Detecting novel mutations that contribute to hearing loss is crucial to providing accurate personalized diagnoses, tailored interventions, and improving prognosis. Copy number variants (CNVs) are structural mutations that are understudied, potential contributors to hearing loss. Here, we present the Abnormal Wobbly Gait (AWG) mouse, the first documented mutant exhibiting waltzer-like locomotor dysfunction, hyperactivity, circling behaviour, and profound deafness caused by a spontaneous CNV deletion in cadherin 23 (Cdh23). We were unable to identify the causative mutation through a conventional whole-genome sequencing (WGS) and variant detection pipeline, but instead found a linked variant in hexokinase 1 (Hk1) that was insufficient to recapitulate the AWG phenotype when introduced into C57BL/6J mice using CRISPR-Cas9. Investigating nearby deafness-associated genes revealed a pronounced downregulation of Cdh23 mRNA and a complete absence of full-length CDH23 protein, which is critical for the development and maintenance of inner ear hair cells, in whole head extracts from AWG neonates. Manual inspection of WGS read depth plots of the Cdh23 locus revealed a putative 10.4 kb genomic deletion of exons 11 and 12 that was validated by PCR and Sanger sequencing. This study underscores the imperative to refine variant detection strategies to permit identification of pathogenic CNVs easily missed by conventional variant calling to enhance diagnostic precision and ultimately improve clinical outcomes for individuals with genetically heterogenous disorders such as hearing loss.

Copy number profiling implicates thin high-grade squamous intraepithelial lesions as a true precursor of cervical human papillomavirus-induced squamous cell cancer.

Full-thickness high-grade squamous intraepithelial lesions (HSIL) are precursors of invasive cervical squamous cell carcinoma (SCC). The World Health Organization (WHO) and Lower Anogenital Squamous Terminology (LAST) Standardization Project for human papilloma virus (HPV)-associated lesions divide full-thickness HSIL of the cervix into thin HSIL with one to nine cell layers thickness and the typical full-thickness HSIL of more than ten cell layers. Although HPV oncogene transcripts and p16ink4a overexpression, as markers of transforming HPV infection, are detectable in thin HSIL, the biological significance of thin HSIL in cervical carcinogenesis remains poorly understood. To further characterize thin HSIL, we performed a comparative study of chromosomal copy number variations (CNV), analysis of dysregulated genes present in the segments with CNV, and a generalized genetic complexity calculation for 31 thin HSIL, 31 thick HSIL, 24 microinvasive SCC (pT1a SCC), and 22 highly invasive SCC samples. Thin HSIL share various CNV and specific dysregulated gene pathways with thick HSIL and invasive SCC. Thin HSIL exhibited an average CNV of 11.6%, compared with 14.1% for thick HSIL, 15.5% for pT1a SCC, and 26.6% for highly invasive SCC. The CNV included gains at 1q and 3q (40 and 43%, respectively), partial loss of 3p, and loss of chromosomes 11 (18%), 16 (50%), 20 (35%), and 22 (40%). Pathways affected solely in thin HSIL were those enhancing immune evasion and primarily involved the interleukin (IL)6, IL21, and IL23 genes. ILs are transiently upregulated in response to infection and play a crucial role in mounting antitumor T-cell activity. Deregulation reflects an attempt by the HPV to evade the initial immune response of the host. The primary pathways shared by thick HSIL and invasive SCC were interactions between lymphoid and non-lymphoid cells, Notch2 signaling, tight junction (TJ) interactions (primarily of the claudin family), and FGR2 alternative splicing. Our results show that thin HSIL carry similar genetic changes as thick HSIL and SCC, indicating that thin HSIL are true precursor lesions that can progress to thick HSIL and SCC.

Genome-Wide Detection of Copy Number Variations and Their Potential Association with Carcass and Meat Quality Traits in Pingliang Red Cattle.

Copy number variation (CNV) serves as a significant source of genetic diversity in mammals and exerts substantial effects on various complex traits. Pingliang red cattle, an outstanding indigenous resource in China, possess remarkable breeding value attributed to their tender meat and superior marbling quality. However, the genetic mechanisms influencing carcass and meat quality traits in Pingliang red cattle are not well understood. We generated a comprehensive genome-wide CNV map for Pingliang red cattle using the GGP Bovine 100K SNP chip. A total of 755 copy number variable regions (CNVRs) spanning 81.03 Mb were identified, accounting for approximately 3.24% of the bovine autosomal genome. Among these, we discovered 270 potentially breed-specific CNVRs in Pingliang red cattle, including 143 gains, 73 losses, and 54 mixed events. Functional annotation analysis revealed significant associations between these specific CNVRs and important traits such as carcass and meat quality, reproduction, exterior traits, growth traits, and health traits. Additionally, our network and transcriptome analysis highlighted CACNA2D1, CYLD, UBXN2B, TG, NADK, and ITGA9 as promising candidate genes associated with carcass weight and intramuscular fat deposition. The current study presents a genome-wide CNV map in Pingliang red cattle, highlighting breed-specific CNVRs, and transcriptome findings provide valuable insights into the underlying genetic characteristics of Pingliang red cattle. These results offer potential avenues for enhancing meat quality through a targeted breeding program.

Evaluation of somatic copy number variation detection by NGS technologies and bioinformatics tools on a hyper-diploid cancer genome.

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.

Copy number variation and population-specific immune genes in the model vertebrate zebrafish.

Copy number variation in large gene families is well characterized for plant resistance genes, but similar studies are rare in animals. The zebrafish (Danio rerio) has hundreds of NLR immune genes, making this species ideal for studying this phenomenon. By sequencing 93 zebrafish from multiple wild and laboratory populations, we identified a total of 1513 NLRs, many more than the previously known 400. Approximately half of those are present in all wild populations, but only 4% were found in 80% or more of the individual fish. Wild fish have up to two times as many NLRs per individual and up to four times as many NLRs per population than laboratory strains. In contrast to the massive variability of gene copies, nucleotide diversity in zebrafish NLR genes is very low: around half of the copies are monomorphic and the remaining ones have very few polymorphisms, likely a signature of purifying selection.

Humans and other animals have immune systems that protect them from bacteria, viruses and other potentially harmful microbes. Members of a family of genes known as the NLR family play various roles in helping to recognize and destroy these microbes. Different species have varying numbers of NLR genes, for example, humans have 22 NLRs, but fish can have hundreds. 400 have been found in the small tropical zebrafish, also known as zebra danios. Zebrafish are commonly used as model animals in research studies because they reproduce quickly and are easy to keep in fish tanks. Much of what we know about fish biology comes from studying strains of those laboratory zebrafish, including the 400 NLRs found in a specific laboratory strain. Many NLRs in zebrafish are extremely similar, suggesting that they have only evolved fairly recently through gene duplication. It remains unclear why laboratory zebrafish have so many almost identical NLRs, or if wild zebrafish also have lots of these genes. To find out more, Schäfer et al. sequenced the DNA of NLRs from almost 100 zebrafish from multiple wild and laboratory populations. The approach identified over 1,500 different NLR genes, most of which, were previously unknown. Computational modelling suggested that each wild population of zebrafish may harbour up to around 2,000 NLR genes, but laboratory strains had much fewer NLRs. The numbers of NLR genes in individual zebrafish varied greatly ­ only 4% of the genes were present in 80% or more of the fish. Many genes were only found in specific populations or single individuals. Together, these findings suggest that the NLR family has expanded in zebrafish as part of an ongoing evolutionary process that benefits the immune system of the fish. Similar trends have also been observed in the NLR genes of plants, indicating there may be an evolutionary strategy across all living things to continuously diversify large families of genes. Additionally, this work highlights the lack of diversity in the genes of laboratory animals compared with those of their wild relatives, which may impact how results from laboratory studies are used to inform conservation efforts or are interpreted in the context of human health.

Human Genetics of Ventricular Septal Defect.

Ventricular septal defects (VSDs) are recognized as one of the commonest congenital heart diseases (CHD), accounting for up to 40% of all cardiac malformations, and occur as isolated CHDs as well as together with other cardiac and extracardiac congenital malformations in individual patients and families. The genetic etiology of VSD is complex and extraordinarily heterogeneous. Chromosomal abnormalities such as aneuploidy and structural variations as well as rare point mutations in various genes have been reported to be associated with this cardiac defect. This includes both well-defined syndromes with known genetic cause (e.g., DiGeorge syndrome and Holt-Oram syndrome) and so far undefined syndromic forms characterized by unspecific symptoms. Mutations in genes encoding cardiac transcription factors (e.g., NKX2-5 and GATA4) and signaling molecules (e.g., CFC1) have been most frequently found in VSD cases. Moreover, new high-resolution methods such as comparative genomic hybridization enabled the discovery of a high number of different copy number variations, leading to gain or loss of chromosomal regions often containing multiple genes, in patients with VSD. In this chapter, we will describe the broad genetic heterogeneity observed in VSD patients considering recent advances in this field.

Correlation between types of ventricular septal defect and chromosomal abnormalities in low-risk non-invasive prenatal testing.

PURPOSE: The aim of this study was to examine whether there is a correlation between different types of ventricular septal defects (VSD) and chromosomal abnormalities in the low-risk setting of non-invasive prenatal testing (NIPT) and to evaluate the prognosis of fetuses with varying types of VSD. METHODS: Cases of pregnant women who underwent amniocentesis due to fetal VSD were collected by Tianjin Central Hospital of Obstetrics and Gynecology from May 2017 to May 2022. Exclusions were made for those without NIPT, with high-risk NIPT results, genetic disorders, and those lost to follow-up. Data collected included ultrasound classification of VSD, prenatal NIPT results, copy-number variations (CNVs) results, and neonatal outcomes. RESULTS: The prevalence of pathogenic CNVs was investigated in 74 cases of VSDs. Of these cases, 45 were isolated VSDs (9 muscular and 36 non-muscular) and 29 were non-isolated VSDs (10 with intracardiac and 19 with extra-cardiac structural anomalies). The results revealed that the incidence of pathogenic CNVs was lower in isolated VSDs compared to non-isolated VSDs in a low-risk NIPT condition (χ2 = 9.344, P = 0.002). There was no significant difference in the prevalence of pathogenic CNVs between VSDs with intracardiac and extra-cardiac structural anomalies (P = 0.541). Moreover, VSDs associated with intracardiac structural anomalies had the highest rate of surgical intervention. CONCLUSION: When NIPT is low-risk and VSD is isolated, the likelihood of fetal chromosomal defects is not increased. However, if there are intra- or extra-cardiac structural abnormalities present alongside VSD, the possibility of pathogenic CNV is considerably greater, necessitating invasive prenatal diagnosis. Isolated muscular VSDs usually do not require surgery, which can be used as a basis for prenatal counseling regarding fetal VSD.

Single-nucleus multi-omic profiling of polyploid heart nuclei identifies fusion-derived cardiomyocytes in the human heart.

Understanding the mechanisms of polyploidization in cardiomyocytes is crucial for advancing strategies to stimulate myocardial regeneration. Although endoreplication has long been considered the primary source of polyploid human cardiomyocytes, recent animal work suggests the potential for cardiomyocyte fusion. Moreover, the effects of polyploidization on the genomic-transcriptomic repertoire of human cardiomyocytes have not been studied previously. We applied single-nuclei whole genome sequencing, single nuclei RNA sequencing, and multiome ATAC + gene expression (from the same nuclei) techniques to nuclei isolated from 11 healthy hearts. Utilizing post-zygotic non-inherited somatic mutations occurring during development as "endogenous barcodes," to reconstruct lineage relationships of polyploid cardiomyocytes. Of 482 cardiomyocytes from multiple healthy donor hearts 75.7% can be sorted into several developmental clades marked by one or more somatic single-nucleotide variants (SNVs). At least ~10% of tetraploid cardiomyocytes contain cells from distinct clades, indicating fusion of lineally distinct cells, whereas 60% of higher-ploidy cardiomyocytes contain fused cells from distinct clades. Combined snRNA-seq and snATAC-seq revealed transcriptome and chromatin landscapes of polyploid cardiomyocytes distinct from diploid cardiomyocytes, and show some higher-ploidy cardiomyocytes with transcriptional signatures suggesting fusion between cardiomyocytes and endothelial and fibroblast cells. These observations provide the first evidence for cell and nuclear fusion of human cardiomyocytes, raising the possibility that cell fusion may contribute to developing or maintaining polyploid cardiomyocytes in the human heart.

Does the esv3587290 Copy Number Variation in the VANGL1 Gene Differ as a Genetic Factor for Developing Nephritis in Mexican Childhood-Onset Systemic Lupus Erythematosus Patients?

A ~3-kb deletion-type DNA copy number variation (CNV, esv3587290) located at intron 7 of the VANGL1 gene (1p13.1, MIM*610132) has been proposed as a genetic factor in lupus nephritis (LN) development in adult systemic lupus erythematosus (SLE) patients across European-descent populations, but its replication in other ethnicities has been inconsistent and its association with LN in childhood-onset SLE (cSLE) remains unknown. Here, we performed an exploratory association study in a sample of 66 unrelated cSLE Mexican patients (11 males, 55 females; ages 7.8 to 18.6 years). Two stratified groups were compared: cSLE patients with (N = 39) or without (N = 27) LN, as diagnosed by renal biopsy (N = 17), proteinuria (N = 33), urinary protein-creatinine ratio > 0.2 (N = 34), and erythrocyturia and/or granular casts in urinary sediment (N = 16). For esv3587290 CNV genotyping, we performed an end-point PCR assay with breakpoint confirmation using Sanger sequencing. We also determined the allelic frequencies of the esv3587290 CNV in 181 deidentified ethnically matched individuals (reference group). The obtained genotypes were tested for Hardy-Weinberg equilibrium using the χ2 test. Associations between LN and esv3587290 CNV were tested by calculating the odds ratio (OR) and using Pearson's χ2 tests, with a 95% confidence interval and p ≤ 0.05. The esv3587290 CNV allele (OR 0.108, 95% CI 0.034-0.33, p = 0.0003) and the heterozygous genotype (OR 0.04, 95% CI 0.119-0.9811, p = 0.002) showed a significant protective effect against LN development. Finally, we characterized the precise breakpoint of the esv3587290 CNV to be NG_016548.1(NM_138959.3):c.1314+1339_1315-897del in our population. This report supports the notion that a broad genetic heterogeneity underlies the susceptibility for developing LN.

Identifying Genetic Etiology in Patients with Intellectual Disability: An Experience in Public Health Services in Northeastern Brazil.

Intellectual disability (ID) is considered a common neuropsychiatric disorder that affects up to 3% of the population. The etiologic origin of ID may be genetic, environmental, and multifactorial. Chromosomopathies are relatively common among the genetic causes of ID, especially in the most severe cases and those associated with dysmorphic features. Currently, the application of new molecular cytogenetics technologies has increasingly allowed the identification of microdeletions, microduplications, and unbalanced translocations as causes of ID. The objective of this study was to investigate the etiology of ID in patients admitted to a public hospital in Northeastern Brazil. In total, 119 patients with ID who had normal karyotypes and fragile X exams participated in this study. The patients were initially physically examined for microdeletion syndromes and then tested using fluorescence in situ hybridization (FISH), multiplex ligation-dependent probe amplification (MLPA), methylation-sensitive polymerase chain reaction (MS-PCR), and chromosome microarray analysis (CMA), according to clinical suspicion. Patients with no diagnoses after FISH, MLPA, and/or MS-PCR evaluations were subsequently tested by CMA. The rate of etiologic diagnoses of ID in the current study was 28%. FISH diagnosed 25 out of 79 tested (31%), MLPA diagnosed 26 out of 79 tested (32%), MS-PCR diagnosed 7 out of 20 tested (35%), and the single nucleotide polymorphism array diagnosed 6 out of 27 tested (22%). Although the CMA is the most complete and recommended tool for the diagnosis of microdeletions, microduplications, and unbalance translocations in patients with ID, FISH, MLPA, and MS-PCR testing can be used as the first tests for specific syndromes, as long as the patients are first physically screened clinically, especially in the public health networks system in Brazil, where resources are scarce.

Detection and characterization of copy number variation in three differentially-selected Nellore cattle populations.

Introduction: Nellore cattle (Bos taurus indicus) is the main beef cattle breed raised in Brazil. This breed is well adapted to tropical conditions and, more recently, has experienced intensive genetic selection for multiple performance traits. Over the past 43 years, an experimental breeding program has been developed in the Institute of Animal Science (IZ, Sertaozinho, SP, Brazil), which resulted in three differentially-selected lines known as Nellore Control (NeC), Nellore Selection (NeS), and Nellore Traditional (NeT). The primary goal of this selection experiment was to determine the response to selection for yearling weight (YW) and residual feed intake (RFI) on Nellore cattle. The main objectives of this study were to: 1) identify copy number variation (CNVs) in Nellore cattle from three selection lines; 2) identify and characterize CNV regions (CNVR) on these three lines; and 3) perform functional enrichment analyses of the CNVR identified. Results: A total of 14,914 unique CNVs and 1,884 CNVRs were identified when considering all lines as a single population. The CNVRs were non-uniformly distributed across the chromosomes of the three selection lines included in the study. The NeT line had the highest number of CNVRs (n = 1,493), followed by the NeS (n = 823) and NeC (n = 482) lines. The CNVRs covered 23,449,890 bp (0.94%), 40,175,556 bp (1.61%), and 63,212,273 bp (2.54%) of the genome of the NeC, NeS, and NeT lines, respectively. Two CNVRs were commonly identified between the three lines, and six, two, and four exclusive regions were identified for NeC, NeS, and NeT, respectively. All the exclusive regions overlap with important genes, such as SMARCD3, SLC15A1, and MAPK1. Key biological processes associated with the candidate genes were identified, including pathways related to growth and metabolism. Conclusion: This study revealed large variability in CNVs and CNVRs across three Nellore lines differentially selected for YW and RFI. Gene annotation and gene ontology analyses of the exclusive CNVRs to each line revealed specific genes and biological processes involved in the expression of growth and feed efficiency traits. These findings contribute to the understanding of the genetic mechanisms underlying the phenotypic differences among the three Nellore selection lines.

High prevalence of copy number variations in the Japanese participants with suspected MODY.

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.

Fetal agenesis of the corpus callosum: Clinical and genetic analysis in a series of 40 patients.

OBJECTIVES: This study aimed to explore the genetic causes of agenesis of the corpus callosum (ACC) and assess the utility of karyotype analysis, copy number variation sequencing (CNV-seq), and whole-exome sequencing (WES) to genetically diagnose fetal ACC. METHODS: We retrospectively examined 40 fetuses diagnosed with ACC who underwent prenatal ultrasonography or magnetic resonance imaging between January 2019 and October 2023. Genetic tests were conducted on the fetuses using karyotype analysis or CNV-seq as the first-line diagnosis. WES was performed if aneuploid and pathogenic CNVs were excluded. RESULTS: Among the 40 fetuses, 29 (72 %) had non-isolated ACC and 11 (28 %) had isolated ACC. Cerebellar dysplasia and hydrocephalus were the most common abnormal developments in the central nervous system. Twenty-eight patients underwent karyotype analysis, with a detection rate of 14 % (4/28). Twenty-six patients underwent CNV-seq; three patients were found to have pathogenic CNVs, with a detection rate of 12 % (3/26). Thirty-three fetuses with no findings of karyotype analysis or CNV-seq were subsequently tested using WES, with a detection rate of 36 % (12/33). Overall, the total diagnostic yield was 48 % (19/40), and monogenic etiology accounted for 30 % (12/40). The genetic detection rate of fetal non-isolated ACC (62 %, 18/29) was higher than that of isolated ACC (9 %, 1/11). CONCLUSION: Prenatal genetic analysis of fetuses with ACC is clinically significant, with monogenic disorders being the main cause. WES may enhance the detection rate of fetuses with ACC with negative karyotype analysis or CNV-seq results.

Genome-wide CNV 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 variant (CNVs) spanning 133.29 Mb of the buffalo genome, resulting in 2,100 copy number variant regions (CNVRs), with 1,993 shared CNVRs being found within the studied buffalo types. Analyzing CNVRs 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 CNVRs that exhibited significant divergence between different buffalo breeds, including genes linked to milk production traits. Second, expression quantitative loci (eQTL) analysis revealed differential expression of CNVR-driven genes (DECGs) associated with milk production traits. Notably, known milk production-related genes were among these DECGs, validating their relevance. Last, a genome-wide association study (GWAS) identified 3 CNVRs 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.

Efficiency of copy number variation sequencing combined with karyotyping in fetuses with congenital heart disease and the following outcomes.

BACKGROUND: Both copy number variant-sequencing (CNV-seq) and karyotype analysis have been used as powerful tools in the genetic aetiology of fetuses with congenital heart diseases (CHD). However, CNV-seq brings clinicians more confusions to interpret the detection results related to CHD with or without extracardiac abnormalities. Hence, we conducted this study to investigate the clinical value of CNV-seq in fetuses with CHD. RESULTS: A total of 167 patients with fetal CHD including 36 single CHD (sCHD), 41 compound CHD (cCHD) and 90 non-isolated CHD (niCHD) were recruited into the study. 28 cases (16.77%, 28/167) were revealed with chromosomal abnormalities at the level of karyotype. The pathogenic detection rate (DR) of CNV-seq (23.17%, 19/82) was higher than that of karyotyping (15.85%, 13/82) in 82 cases by CNV-seq and karyotyping simultaneously. The DR of pathogenic copy number variations (PCNVs) (31.43%) was higher in niCHD subgroup than that in sCHD and cCHD (9.52% and 23.08%). Conotruncal defect (CTD) was one of the most common heart malformations with the highest DR of PCNVs (50%) in 7 categories of CHD. In terms of all the pregnancy outcomes, 67 (40.12%) cases were terminated and 100 (59.88%) cases were live neonates. Only two among 34 cases with a pathogenic genetic result chose to continue the pregnancy. CONCLUSIONS: CNV-seq combined with karyotyping is a reliable and accurate prenatal technique for identifying pathogenic chromosomal abnormalities associated with fetal CHD with or without extracardiac abnormalities, which can assist clinicians to perform detailed genetic counselling with regard to the etiology and related outcomes of CHD.

Identification of the hybrid gene LILRB5-3 by long-read sequencing and implication of its novel signaling function.

Leukocyte immunoglobulin (Ig)-like receptors (LILRs) on human chromosome 19q13.4 encode 11 immunoglobulin superfamily receptors, exhibiting genetic diversity within and between human populations. Among the LILR genes, the genomic region surrounding LILRB3 and LILRA6 has yet to be fully characterized due to their significant sequence homology, which makes it difficult to differentiate between them. To examine the LILRB3 and LILRA6 genomic region, a tool named JoGo-LILR CN Caller, which can call copy number from short-read whole genome sequencing (srWGS) data, was applied to an extensive international srWGS dataset comprising 2,504 samples. During this process, a previously unreported loss of both LILRB3 and LILRA6 was detected in three samples. Using long-read sequencing of these samples, we have discovered a novel large deletion (33,692 bp) in the LILRB3 and LILRA6 genomic regions in the Japanese population. This deletion spanned three genes, LILRB3, LILRA6, and LILRB5, resulting in LILRB3 exons 12-13 being located immediately downstream of LILRB5 exons 1-12 with the loss of LILRA6, suggesting the potential expression of a hybrid gene between LILRB5 and LILRB3 (LILRB5-3). Transcription and subsequent translation of the LILRB5-3 hybrid gene were also verified. The hybrid junction was located within the intracellular domain, resulting in an LILRB5 extracellular domain fused to a partial LILRB3 intracellular domain with three immunoreceptor tyrosine-based inhibitory motifs (ITIMs), suggesting that LILRB5-3 acquired a novel signaling function. Further application of the JoGo-LILR tool to srWGS samples suggested the presence of the LILRB5-3 hybrid gene in the CEU population. Our findings provide insight into the genetic and functional diversity of the LILR family.