Effects Of Copy Number Variation Research Articles

Effects of Copy Number Variations in the Plectin (PLEC) Gene on the Growth Traits and Meat Quality of Leizhou Black Goats

The plectin (PLEC) gene is crucial in regulating muscle development and maintaining the cytoskeleton. An abnormal expression of PLEC can lead to muscle atrophy and muscular dystrophy. In a previous study, we found that Leizhou black goats exhibit abundant structural variations in the PLEC gene. However, the genetic effects of these variations on growth traits and meat quality in goats are not fully understood. In this study, three PLEC copy number variations (CNVs) were identified in a population of 417 Leizhou black goats, using quantitative polymerase chain reaction (qPCR) technology. Population distribution analysis revealed a high abundance of various types of these three CNVs. PLEC mRNA was found to be highly expressed in muscle tissue and remained consistently high from 1 month to 24 months after birth. Specifically, the gain type of CNV-1 (chr14: 81056401-81064800) showed a significant association with PLEC mRNA expression in muscle tissue (p < 0.01). The sequence of CNV-1 in PLEC shares similarities with three domain superfamilies associated with muscle development and skin disease. Furthermore, there were significant differences in chest circumference, body weight, carcass weight, the cross-sectional area of the longissimus dorsi lumbar muscle, and shear stress between different types of CNV-1 (p < 0.05). Notably, goats with the CNV-1 gain type demonstrated better phenotypic values compared to those with loss and normal types. These findings suggest that PLEC CNV-1 could play a crucial role in the growth and muscle development of Leizhou black goats, making it a potential marker for assisted selection in goat breeding.

Effects of copy number variations on longevity in late-onset Alzheimer's disease patients: insights from a causality network analysis.

Alzheimer's disease (AD), particularly late-onset Alzheimer's disease (LOAD), is a prevalent form of dementia that significantly affects patients' cognitive and behavioral capacities and longevity. Although approximately 70 genetic risk factors linked with AD have been identified, their influence on patient longevity remains unclear. Further, recent studies have associated copy number variations (CNVs) with the longevity of healthy individuals and immune-related pathways in AD patients. This study aims to investigate the role of CNVs on the longevity of AD patients by integrating the Whole Genome Sequencing (WGS) and transcriptomics data from the Religious Orders Study/Memory and Aging Project (ROSMAP) cohort through causality network inference. Our comprehensive analysis led to the construction of a CNV-Gene-Age of Death (AOD) causality network. We successfully identified three key CNVs (DEL5006, mCNV14192, and DUP42180) and seven AD-longevity causal genes (PLGRKT, TLR1, PLAU, CALB2, SYTL2, OTOF, and NT5DC1) impacting AD patient longevity, independent of disease severity. This outcome emphasizes the potential role of plasminogen activation and chemotaxis in longevity. We propose several hypotheses regarding the role of identified CNVs and the plasminogen system on patient longevity. However, experimental validation is required to further corroborate these findings and uncover precise mechanisms. Despite these limitations, our study offers promising insights into the genetic influence on AD patient longevity and contributes to paving the way for potential therapeutic interventions.

Recurrent copy number variations in the human fungal pathogen Candida parapsilosis.

Candida parapsilosis is an opportunistic fungal pathogen with increasing incidence in hospital settings worldwide; however, we lack a comprehensive understanding of the mechanisms promoting its virulence and drug resistance. Bergin et al. systematically quantify the frequency and effect of copy number variation (CNV) across 170 diverse clinical and environmental isolates of C. parapsilosis (Bergin SA, Zhao F, Ryan AP, Müller CA, Nieduszynski CA, Zhai B, Rolling T, Hohl TM, Morio F, Scully J, Wolfe KH, Butler G, 2022, mBio, https://doi.org/10.1128/mbio.01777-22). Using a combination of both short- and long-read whole genome sequencing techniques, they determine the structure and copy number of two CNVs that arose recurrently throughout the evolution of these isolates. Each CNV predominantly amplifies one coding sequence (ARR3 or RTA3); however, the amplitude and recombination breakpoints are variable across the isolates. Amplification of RTA3 correlates with drug resistance and deletion causes drug susceptibility. This study highlights the need for further research into the mechanisms and dynamics of CNV formation and the impact of these CNVs on virulence and drug resistance across diverse fungal pathogens.

17.3 The Effect of Copy Number Variants on Cognitive Ability and ASD Is Determined by the Spatial Temporal and Cell-Type Gene Expression Patterns

17.3 The Effect of Copy Number Variants on Cognitive Ability and ASD Is Determined by the Spatial Temporal and Cell-Type Gene Expression Patterns

Combination of expert guidelines-based and machine learning-based approaches leads to superior accuracy of automated prediction of clinical effect of copy number variations

Clinical interpretation of copy number variants (CNVs) is a complex process that requires skilled clinical professionals. General recommendations have been recently released to guide the CNV interpretation based on predefined criteria to uniform the decision process. Several semiautomatic computational methods have been proposed to recommend appropriate choices, relieving clinicians of tedious searching in vast genomic databases. We have developed and evaluated such a tool called MarCNV and tested it on CNV records collected from the ClinVar database. Alternatively, the emerging machine learning-based tools, such as the recently published ISV (Interpretation of Structural Variants), showed promising ways of even fully automated predictions using broader characterization of affected genomic elements. Such tools utilize features additional to ACMG criteria, thus providing supporting evidence and the potential to improve CNV classification. Since both approaches contribute to evaluation of CNVs clinical impact, we propose a combined solution in the form of a decision support tool based on automated ACMG guidelines (MarCNV) supplemented by a machine learning-based pathogenicity prediction (ISV) for the classification of CNVs. We provide evidence that such a combined approach is able to reduce the number of uncertain classifications and reveal potentially incorrect classifications using automated guidelines. CNV interpretation using MarCNV, ISV, and combined approach is available for non-commercial use at https://predict.genovisio.com/.

Population history modulates the fitness effects of Copy Number Variation in the Roma

We provide the first whole genome Copy Number Variant (CNV) study addressing Roma, along with reference populations from South Asia, the Middle East and Europe. Using CNV calling software for short-read sequence data, we identified 3171 deletions and 489 duplications. Taking into account the known population history of the Roma, as inferred from whole genome nucleotide variation, we could discern how this history has shaped CNV variation. As expected, patterns of deletion variation, but not duplication, in the Roma followed those obtained from single nucleotide polymorphisms (SNPs). Reduced effective population size resulting in slightly relaxed natural selection may explain our observation of an increase in intronic (but not exonic) deletions within Loss of Function (LoF)-intolerant genes. Over-representation analysis for LoF-intolerant gene sets hosting intronic deletions highlights a substantial accumulation of shared biological processes in Roma, intriguingly related to signaling, nervous system and development features, which may be related to the known profile of private disease in the population. Finally, we show the link between deletions and known trait-related SNPs reported in the genome-wide association study (GWAS) catalog, which exhibited even frequency distributions among the studied populations. This suggests that, in general human populations, the strong association between deletions and SNPs associated to biomedical conditions and traits could be widespread across continental populations, reflecting a common background of potentially disease/trait-related CNVs.

Transcriptome and pan-cancer system analysis identify PM2.5-induced stanniocalcin 2 as a potential prognostic and immunological biomarker for cancers.

Epidemiological studies have shown that air pollution and particulate matter (PM) are closely related to the occurrence of cancer. However, the potential prognostic and immunological biomarkers for air pollution related cancers are lacking. In this study, we proved PM2.5 exposure was correlated with lung cancer through transcriptome analysis. Importantly, we identified STC2 as a key gene regulated by PM2.5, whose expression in epithelial cells was significantly increased after PM2.5 treatment and validated by using RT-qPCR and immunofluorescence. Kaplan-Meier OS curves suggested that high STC2 expression positively correlated with a poor prognosis in lung cancer. Furthermore, we discovered that STC2 was associated with multiple cancers and pathways in cancer. Next, Pan-Cancer Expression Landscape of STC2 showed that STC2 exhibited inconsistent expression across 26 types of human cancer, lower in KIRP in cancer versus adjacent normal tissues, and significantly higher in another cancers. Cox regression results suggested that STC2 expression was positively or negatively associated with prognosis in different cancers. Moreover, STC2 expression was associated with clinical phenotypes including age, gender, stage and grade. Mutation features of STC2 were also analyzed, in which the highest alteration frequency of STC2 was presented in KIRC with amplification. Meanwhile, the effects of copy number variation (CNV) on STC2 expression were investigated across various tumor types, suggesting that STC2 expression was significantly correlated with CNV in tumors. Additionally, STC2 was closely related to tumor heterogeneity, tumor stemness and tumor immune microenvironment like immune cell infiltration. In the meantime, we analyzed methylation modifications and immunological correlation of STC2. The results demonstrated that STC2 expression positively correlated with most RNA methylation genes and immunomodulators across tumors. Taken together, the findings revealed that PM2.5-induced STC2 might be a potential prognostic and immunological biomarker for cancers related to air pollution.

Genetic factors associated with serum amylase in a Japanese population: combined analysis of copy-number and single-nucleotide variants

Amylase activity and levels in humans are heritable quantitative traits. Although many studies exist on the effects of copy-number variants (CNVs) in amylase genes (AMY) on human phenotypes, such as body mass index (BMI), the genetic factors controlling interindividual variation in amylase levels remain poorly understood. Here, we conducted a genome-wide association study (GWAS) of serum amylase levels (SAL) in 814 Japanese individuals to identify associated single-nucleotide variants (SNVs), after adjusting for non-genetic factors. Diploid copy numbers (CN) of AMY (AMY1, AMY2A, and AMY2B) were measured using droplet digital PCR to examine the association between each diploid CN and SAL. We further assessed the relative contribution of the GWAS-lead SNV and AMY CNVs to SAL. GWAS identified 14 significant SNVs (p < 5 × 10−8) within a linkage disequilibrium block near the AMY cluster on chromosome 1. The association analyses of AMY CNVs and SAL showed a significant association between AMY1 diploid CN and SAL (p = 1.89 × 10−19), while no significant association with SAL was found for AMY2A CN (p = 0.54) or AMY2B CN (p = 0.15). In a joint association analysis with SAL using the GWAS-lead SNV and AMY1 diploid CN, AMY1 CN remained significant (p = 5.4 ×10−13), while the association of the lead SNV was marginal (p = 0.08). We also found no association between AMY1 diploid CN and BMI (p = 0.14). Our results indicate that AMY1 CNV is the major genetic factor for Japanese SAL, with no significant association with BMI.

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

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

Cancer cells modulate the effects of copy number variation on protein abundance

Cancer cells modulate the effects of copy number variation on protein abundance

Evaluation of the Effect of Gene Duplication by Genome Editing on Drug Resistance in Plasmodium falciparum.

The emergence and spread of drug-resistant Plasmodium falciparum have compromised antimalarial efficacy and threatened the global malaria elimination campaign using artemisinin combination therapies. The impacts of amino acid substitutions in antimalarial drug resistance-associated genes on drug susceptibility have been investigated; however, the effects of amplification of those genes remain unexplored due to the lack of robust genetic approaches. Here, we generated transgenic P. falciparum parasites with an additional copy of a drug resistance-associated gene using the highly efficient CRISPR/Cas9 system and investigated their drug response. Insertion of a drug resistance-associated gene expression cassette in the genome resulted in a roughly twofold increase of mRNA levels of the target gene mdr1, which encodes multidrug resistance protein 1. The gene duplication event contributed to resistance to mefloquine, lumefantrine, and dihydroartemisinin, while the duplication of a genomic region encoding plasmepsin 2 and plasmepsin 3 did not affect resistance to antimalarial drugs, including piperaquine. We further demonstrated that mdr1 mRNA expression levels are strongly associated with mefloquine resistance in several field-derived P. falciparum lines with various genetic backgrounds. This study provides compelling evidence that mdr1 could serve as a molecular marker for the surveillance of mefloquine-resistant parasites. Long DNA integration into parasite genomes using the CRISPR/Cas9 system provides a useful tool for the evaluation of the effect of copy number variation on drug response.

Prevalence and Phenotypic Effects of Copy Number Variants in Isolated Hypogonadotropic Hypogonadism.

The genetic architecture of isolated hypogonadotropic hypogonadism (IHH) has not been completely defined. To determine the role of copy number variants (CNVs) in IHH pathogenicity and define their phenotypic spectrum. Exome sequencing (ES) data in IHH probands (n = 1394) (Kallmann syndrome [IHH with anosmia; KS], n = 706; normosmic IHH [nIHH], n = 688) and family members (n = 1092) at the Reproductive Endocrine Unit and the Center for Genomic Medicine of Massachusetts General Hospital were analyzed for CNVs and single nucleotide variants (SNVs)/indels in 62 known IHH genes. IHH subjects without SNVs/indels in known genes were considered "unsolved." Phenotypes associated with CNVs were evaluated through review of patient medical records. A total of 29 CNVs in 13 genes were detected (overall IHH cohort prevalence: ~2%). Almost all (28/29) CNVs occurred in unsolved IHH cases. While some genes (eg, ANOS1 and FGFR1) frequently harbor both CNVs and SNVs/indels, the mutational spectrum of others (eg, CHD7) was restricted to SNVs/indels. Syndromic phenotypes were seen in 83% and 63% of IHH subjects with multigenic and single gene CNVs, respectively. CNVs in known genes contribute to ~2% of IHH pathogenesis. Predictably, multigenic contiguous CNVs resulted in syndromic phenotypes. Syndromic phenotypes resulting from single gene CNVs validate pleiotropy of some IHH genes. Genome sequencing approaches are now needed to identify novel genes and/or other elusive variants (eg, noncoding/complex structural variants) that may explain the remaining missing etiology of IHH.

Effect of copy number variation of PLA2G2A gene to growth traits in Chinese cattle.

As a member of genetic polymorphism, copy number variation has been a commonly used method in the world for investigating effect of genetic polymorphism on gene expression. Effect of genetic polymorphism made on livestock development has been more and more important in beef cattle molecular breeding. The characteristics of Chinese cattle are excellent meat quality, tolerant to rough feeding, good environmental adaptability and so on. But there are some obvious weaknesses still exist in the process of cattle growth and development, such as weak hindquarters and growth slowly. To improve the growth performance and market competitiveness of Chinese cattle, a lot of studies have been made about finding and investigating effective molecular marker. In this study, Q-PCR and data association analysis were used for PLA2G2A gene copy number variation detection and related effect analysis in Chinese cattle. Results showed that PLA2G2A gene has a significant effect on two breeds of Chinese cattle on growth traits, which could be a basic materials and effective information of cattle molecular markers breeding. PLA2G2A, member of secreted phospholipases A2 (sPLA2) in superfamily of phospholipase A2, could catalyze the process of glycerophospholipids hydrolysis from position of sn-2 acyl with the release of free fatty acids and lysophospholipids. Researches about PLA2G2A gene are mostly focus on disease, including tumors and diabetes, the number of study occurred on animal breeding is weak. In this study, blood samples were collected from five breeds of Chinese cattle (Qingchuan cattle, Xianan cattle, Yunling cattle, Pinan cattle and Guyuan cattle) for PLA2G2A gene CNV type detection. SPSS 20.0 software and method of ANOVA were used to analyzed the association between types of CNV and growth traits. Results reveal that the distribution of different copy number types in different cattle breeds is different. In QC, XN and GY cattle, the frequencies of Deletion and Duplication are about 40%; in YL cattle, the frequency of Deletion type exceeds 60%; in PN cattle, the frequency of Duplication is closed to 80%. Association analysis indicate that CNV of PLA2G2A gene showed a positive effect in cattle growth: in QC cattle, Chest depth with Normal type copy number possess a increased trend (P<0.05); individuals with Deletion type copy number have better performance on Height at sacrum, Heart girth and Body height in GY cattle (P<0.05). The functional role and molecular mechanism of PLA2G2A gene in animal growth and development are still unclear, and it is necessary for processing a further research. This research aims to provide basic materials for molecular breeding of Chinese cattle.

Associations between brain abnormalities and common genetic variants for schizophrenia: a narrative review of structural and functional neuroimaging findings.

Schizophrenia is a common psychotic disease worldwide, and several genetic variants have been suggested to influence susceptibility to schizophrenia. However, the exact genetic and neural mechanisms underlying such relationships remain unclear. Neuroimaging endophenotypes have been considered to be desirable for finding genes influencing brain structure and function. Because of its high penetrance of genetic variants and wide applications in clinic practice, neuroimaging could be a great option in schizophrenia study. In this narrative review, we aim to provide a summary of current progress in neuroimaging findings to detect genetic variants that may influence the brain's structure and function, and thus improve understanding of how this interaction affects the onset of schizophrenia. Multiple common genetic variants for schizophrenia including the ZNF804A, DTNBP1, DAOA, AKT1, NTRK3, and ERBB4 genes are reviewed. In summary, the current investigations have reported effects of these genetic variants on both structural and functional brain abnormalities (e.g., altered grey matter volumes, white matter integrity and functional connectivity patterns) in schizophrenia within multiple brain regions, especially the frontal, temporal and hippocampal areas. These findings facilitate our understanding of schizophrenia and helping us with early diagnosis, prognosis prediction and medication. However, a deep insight into this field necessitate further investigations and there are a number of challenges to be settled, such as enlarging the sample size; solving the statistical, methodological, and technological limitations; performing not only cross-sectional but also longitudinal observations; focusing on the effect of copy number variants; figuring out the complexity of epistasis and gene-environment interaction; achieving clinical translations; and depict the effect of a polygenic risk score beyond those of single genetic markers.

Effects of copy number variations on brain structure and risk for psychiatric illness: Large-scale studies from the ENIGMA working groups on CNVs.

The Enhancing NeuroImaging Genetics through Meta‐Analysis copy number variant (ENIGMA‐CNV) and 22q11.2 Deletion Syndrome Working Groups (22q‐ENIGMA WGs) were created to gain insight into the involvement of genetic factors in human brain development and related cognitive, psychiatric and behavioral manifestations. To that end, the ENIGMA‐CNV WG has collated CNV and magnetic resonance imaging (MRI) data from ~49,000 individuals across 38 global research sites, yielding one of the largest studies to date on the effects of CNVs on brain structures in the general population. The 22q‐ENIGMA WG includes 12 international research centers that assessed over 533 individuals with a confirmed 22q11.2 deletion syndrome, 40 with 22q11.2 duplications, and 333 typically developing controls, creating the largest‐ever 22q11.2 CNV neuroimaging data set. In this review, we outline the ENIGMA infrastructure and procedures for multi‐site analysis of CNVs and MRI data. So far, ENIGMA has identified effects of the 22q11.2, 16p11.2 distal, 15q11.2, and 1q21.1 distal CNVs on subcortical and cortical brain structures. Each CNV is associated with differences in cognitive, neurodevelopmental and neuropsychiatric traits, with characteristic patterns of brain structural abnormalities. Evidence of gene‐dosage effects on distinct brain regions also emerged, providing further insight into genotype–phenotype relationships. Taken together, these results offer a more comprehensive picture of molecular mechanisms involved in typical and atypical brain development. This “genotype‐first” approach also contributes to our understanding of the etiopathogenesis of brain disorders. Finally, we outline future directions to better understand effects of CNVs on brain structure and behavior.

Piecing together the Pediatric Brain Tumor Puzzle

A recent study by Petralia et al. of 218 pediatric brain tumors across seven different entities applied an integrated approach incorporating proteomics, phosphoproteomics, whole-genome sequencing, and RNA sequencing. This elegant study unveiled new signaling pathways, the composition of tumor microenvironments, and functional effects of copy number variants and somatic mutations.

The nine ADAMs family members serve as potential biomarkers for immune infiltration in pancreatic adenocarcinoma.

BackgroundThe functional significance of ADAMs family members in the immune infiltration of pancreatic adenocarcinoma (PAAD) awaits elucidation.MethodsADAMs family members with significant expression were identified among differentially expressed genes of PAAD based on The Cancer Genome Atlas (TCGA) database followed by a verification based on the Oncomine database. The correlation of ADAMs in PAAD was estimated with the Spearman’s rho value. The pathway enrichment of ADAMs was performed by STRING and GSEALite, respectively. The protein–protein interaction and Gene Ontology analyses of ADAMs and their similar genes were exanimated in STRING and visualized by Cytoscape. Subsequently, the Box-Whisker plot was used to show a correlation between ADAMs and different tumor grade 1/2/3/4 with Student’s t-test. TIMER was applied to estimate a correlation of ADAMs expressions with immune infiltrates and immune checkpoint blockade (ICB) immunotherapy-related molecules. Furthermore, the effect of copy number variation (CNV) of ADAMs genes was assessed on the immune infiltration levels.ResultADAM8/9/10/12/15/19/28/TS2/TS12 were over-expressed in PAAD. Most of the nine ADAMs had a significant correlation. ADAM8/12/15/19 expression was remarkably increased in the comparison between grade 1 and grade 2/3 of PAAD. ADAM8/9/10/12/19/28/TS2/TS12 had a positive correlation with almost five immune infiltrates. ADAM12/19/TS2/TS12 dramatically related with ICB immunotherapy-related molecules. CNV of ADAMs genes potentially influenced the immune infiltration levels.ConclusionKnowledge of the expression level of the ADAMs family could provide a reasonable strategy for improved immunotherapies to PAAD.

Copy Number Variations and Expression Levels of Guanylate-Binding Protein 6 Gene Associated with Growth Traits of Chinese Cattle

Simple SummaryOur study identified the copy number variations (CNVs) and transcriptional expressions of guanylate-binding protein 6 (GBP6) in Chinese cattle, as well as analyzed the association of CNVs and expressions with the growth traits of Chinese cattle. The results showed Xianan cattle with gain types (Log22−ΔΔCt > 0.5) of CNVs of GBP6 had highest relative gene expression levels in the muscle tissues and displayed superior phenotypic values of body weight, cannon circumference and chest circumference. Our study suggested that CNV gain types of GBP6 could be used as the candidate markers in the cattle-breeding program for growth traits.Association studies have indicated profound effects of copy number variations (CNVs) on various phenotypes in different species. In this study, we identified the CNV distributions and expression levels of guanylate-binding protein 6 (GBP6) associated with the growth traits of Chinese cattle. The results showed that the phenotypic values of body size and weight of Xianan (XN) cattle were higher than those of Nanyang (NY) cattle. The medium CNV types were mostly identified in the XN and NY breeds, but their CNV distributions were significantly different (adjusted p < 0.05). The association analysis revealed that the body weight, cannon circumference and chest circumference of XN cattle had significantly different values in different CNV types (p < 0.05), with CNV gain types (Log22−ΔΔCt > 0.5) displaying superior phenotypic values. We also found that transcription levels varied in different tissues (p < 0.001) and the CNV gain types showed the highest relative gene expression levels in the muscle tissue, consistent with the highest phenotypic values of body weight and cannon circumference among the three CNV types. Consequently, our results suggested that CNV gain types of GBP6 could be used as the candidate markers in the cattle-breeding program for growth traits.

Differential Allele-Specific Expression Uncovers Breast Cancer Genes Dysregulated by Cis Noncoding Mutations.

Identifying cancer-relevant mutations in noncoding regions is challenging due to the large numbers of such mutations, their low levels of recurrence, and difficulties in interpreting their functional impact. To uncover genes that are dysregulated due to somatic mutations in cis, we build upon the concept of differential allele-specific expression (ASE) and introduce methods to identify genes within an individual's cancer whose ASE differs from what is found in matched normal tissue. When applied to breast cancer tumor samples, our methods detect the known allele-specific effects of copy number variation and nonsense-mediated decay. Further, genes that are found to recurrently exhibit differential ASE across samples are cancer relevant. Genes with cis mutations are enriched for differential ASE, and we find 147 potentially functional noncoding mutations cis to genes that exhibit significant differential ASE. We conclude that differential ASE is a promising means for discovering gene dysregulation due to cis noncoding mutations.

Intragenic CNTN4 copy number variants associated with a spectrum of neurobehavioral phenotypes

Deletions and duplications involving the CNTN4 gene, which encodes for the contactin 4 protein, have been reported in children with autism spectrum disorder (ASD) and other neurodevelopmental phenotypes. In this study, we performed clinical and genetic characterization of three individuals from unrelated families with copy number variants (CNV) (one deletion and two duplications) within CNTN4. The patients exhibited cognitive delay (3/3), growth restriction (3/3), motor delay (2/3), and febrile seizure/epilepsy (2/3). In contrast to previous reports, all probands presented with speech apraxia or delay with no diagnosis of ASD. Parental studies for the proband with the deletion and one of the 2 probands with the duplication revealed paternal origin of the CNTN4 CNV. Interestingly, previously documented CNV involving this gene were mostly inherited from unaffected fathers, raising questions regarding reduced penetrance and potential parent-of-origin effect. Our findings are compared with previously reported patients and patients in the DECIPHER database. The speech impairment in the three probands suggests a role for CNTN4 in language development. We discuss potential factors contributing to phenotypic heterogeneity and reduced penetrance and attempt to find possible genotype-phenotype correlation. Larger cohorts are needed for comprehensive and unbiased phenotyping and molecular characterization that may lead to better understanding of the underlying mechanisms of reduced penetrance, variable expressivity, and potential parent-of-origin effect of copy number variants encompassing CNTN4.