Presence Of Copy Number Variants Research Articles

Presence of Copy Number Variants Associated With Esotropia in Patients With Exotropia

Strabismus is a common ocular disorder of childhood. There is a clear genetic component to strabismus, but it is not known if esotropia and exotropia share genetic risk factors. To determine whether genetic duplications associated with esotropia are also associated with exotropia. This was a cross-sectional study conducted from November 2005 to December 2023. Individuals with constant or intermittent exotropia of any magnitude or a history of surgery for exotropia were recruited from pediatric ophthalmic practices. Data were analyzed from March to December 2023. Genetic duplication. Presence of genetic duplications at 2p11.2, 4p15.2, and 10q11.22 assessed by digital droplet polymerase chain reaction. Orthoptic measurements and history of strabismus surgery were performed. A total of 234 individuals (mean [SD] age, 19.5 [19.0] years; 127 female [54.3%]) were included in this study. The chromosome 2 duplication was present in 1.7% of patients with exotropia (4 of 234; P = .40), a similar proportion to the 1.4% of patients with esotropia (23 of 1614) in whom it was previously reported and higher than the 0.1% of controls (4 of 3922) previously reported (difference, 1.6%; 95% CI, 0%-3.3%; P < .001). The chromosome 4 duplication was present in 3.0% of patients with exotropia (7 of 234; P = .10), a similar proportion to the 1.7% of patients with esotropia (27 of 1614) and higher than the 0.2% of controls (6 of 3922) in whom it was previously reported (difference, 2.8%; 95% CI, 0.6%-5.0%; P < .001). The chromosome 10 duplication was present in 6.0% of patients with exotropia (14 of 234; P = .08), a similar proportion to the 4% of patients with esotropia (64 of 1614) and higher than the 0.4% of controls (18 of 3922) in whom it was previously reported (difference, 5.6%; 95% CI, 2.5%-8.6%; P < .001). Individuals with a duplication had higher mean (SD) magnitude of deviation (31 [13] vs 22 [14] prism diopters [PD]; difference, 9 PD; 95% CI, 1-16 PD; P = .03), were more likely to have constant (vs intermittent) exotropia (70% vs 29%; difference, 41%; 95% CI, 20.8%-61.2%; P < .001), and had a higher rate of exotropia surgery than those without a duplication (58% vs 34%; difference, 24%; 95% CI, 3%-44%; P = .02). In this cross-sectional study, results suggest that the genetic duplications on chromosomes 2, 4, and 10 were risk factors for exotropia as well as esotropia. These findings support the possibility that esotropia and exotropia have shared genetic risk factors. Whether esotropia or exotropia develops in the presence of these duplications may be influenced by other shared or independent genetic variants or by environmental factors.

Contribution of copy number variants to the genetic etiology of dilated cardiomyopathy

Abstract Background/Introduction Dilated cardiomyopathy (DCM) is one of the most common cardiomyopathies, with a prevalence of 1 in 250. Single nucleotide variants (missense, nonsense, or frameshift) in definitive DCM genes are the most frequently identified in positive genetic studies. Despite major advances in genetic diagnostic technologies in recent decades, and more than 100 genes potentially associated with DCM identified to date, the diagnostic yield is still lower than 40%. Copy number variants (CNVs) have been identified as cause of the disease in several studies, especially in genes intolerant to haploinsufficiency. However, the contribution of CNVs to the etiology of DCM has not been evaluated in large cohorts of DCM patients and remains unknown. Purpose To determine the percentage of cases in which CNVs have been identified as the cause of disease in a large cohort of DCM patients evaluated genetically. Methods Retrospective study in which the presence of CNVs was evaluated in a cohort of patients diagnosed with DCM and studied by next-generation-sequencing (NGS). CNVs were detected using a read depth approach and confirmed by an orthogonal molecular technique (Sanger sequencing, MLPA, or dPCR). Results 6,512 DCM patients referred for genetic testing were sequenced in our center by NGS, and the median age at the time of genetic study was 51.1 years (range 0-93). A pathogenic (P) or likely-pathogenic (LP) disease-causing variant was identified in 1,495 of the patients (diagnostic yield of 22.9%). In 67 (1.02%) of the cases, a CNV classified as P/LP was identified. DMD was the gene with the highest number of CNVs identified (n=40), followed by TTN (n=7) and LMNA (n=5). Of the 67 cases, 45 (67%) corresponded to deletions, with a lesser extent to duplications (n=9; 13%), complex alterations (n=7; 11%), and deletion/insertion (n=6; 9%). Only two of the patients harboring disease-causing CNVs had other variants potentially associated with DCM: a likely-pathogenic variant in MYH7 and a homozygous pathogenic variant in the ALMS1 gene. Conclusion In our cohort of DCM patients, 1.02% harbored CNVs considered P/LP, being DMD the most frequently identified gene. These results show that systematic CNVs analysis in NGS studies improves the yield of genetic testing in DCM and highlights the importance of a complete genetic study that includes structural variant detection in the diagnosis of the disease.

A comprehensive framework for detecting copy number variants from single nucleotide polymorphism data: 'rCNV', a versatile r package for paralogue and CNV detection.

Recent studies have highlighted the significant role of copy number variants (CNVs) in phenotypic diversity, environmental adaptation and species divergence across eukaryotes. The presence of CNVs also has the potential to introduce genotyping biases, which can pose challenges to accurate population and quantitative genetic analyses. However, detecting CNVs in genomes, particularly in non-model organisms, presents a formidable challenge. To address this issue, we have developed a statistical framework and an accompanying r software package that leverage allelic-read depth from single nucleotide polymorphism (SNP) data for accurate CNV detection. Our framework capitalises on two key principles. First, it exploits the distribution of allelic-read depth ratios in heterozygotes for individual SNPs by comparing it against an expected distribution based on binomial sampling. Second, it identifies SNPs exhibiting an apparent excess of heterozygotes under Hardy-Weinberg equilibrium. By employing multiple statistical tests, our method not only enhances sensitivity to sampling effects but also effectively addresses reference biases, resulting in optimised SNP classification. Our framework is compatible with various NGS technologies (e.g. RADseq, Exome-capture). This versatility enables CNV calling from genomes of diverse complexities. To streamline the analysis process, we have implemented our framework in the user-friendly r package 'rCNV', which automates the entire workflow seamlessly. We trained our models using simulated data and validated their performance on four datasets derived from different sequencing technologies, including RADseq (Chinook salmon-Oncorhynchus tshawytscha), Rapture (American lobster-Homarus americanus), Exome-capture (Norway spruce-Picea abies) and WGS (Malaria mosquito-Anopheles gambiae).

Extensive MHC class IIβ diversity across multiple loci in the small-spotted catshark (Scyliorhinus canicula)

The major histocompatibility complex (MHC) is a multigene family responsible for pathogen detection, and initiation of adaptive immune responses. Duplication, natural selection, recombination, and their resulting high functional genetic diversity spread across several duplicated loci are the main hallmarks of the MHC. Although these features were described in several jawed vertebrate lineages, a detailed MHC IIβ characterization at the population level is still lacking for chondrichthyans (chimaeras, rays and sharks), i.e. the most basal lineage to possess an MHC-based adaptive immune system. We used the small-spotted catshark (Scyliorhinus canicula, Carcharhiniformes) as a case-study species to characterize MHC IIβ diversity using complementary molecular tools, including publicly available genome and transcriptome datasets, and a newly developed high-throughput Illumina sequencing protocol. We identified three MHC IIβ loci within the same genomic region, all of which are expressed in different tissues. Genetic screening of the exon 2 in 41 individuals of S. canicula from a single population revealed high levels of sequence diversity, evidence for positive selection, and footprints of recombination. Moreover, the results also suggest the presence of copy number variation in MHC IIβ genes. Thus, the small-spotted catshark exhibits characteristics of functional MHC IIβ genes typically observed in other jawed vertebrates.

An association between copy number variation of enhancer involved in craniofacial development and biogeographic ancestry.

Human facial morphology is a combination of many complex traits and is determined by a large number of genes and enhancers. Here, we report a Copy Number Variation (CNV) study of enhancer hs1431 in populations of Central European and South Siberian ancestry. Central European samples included 97 Poles, while South Siberian samples included 78 Buryats and 27 Tuvinians. CNVs were detected by real-time PCR, using ViiA™ 7 Real-Time PCR System (Applied Biosystems). We revealed significant differences in CNV of hs1431 enhancer between Polish and Buryat population (p=0.0378), but not between Central European and South Siberian population (p=0.1225). Our results suggest that an increase in copy number variation of hs1431 enhancer is associated with biogeographic ancestry. However, this result needs extending and replicating in larger cohorts. This is the first study revealing the presence of copy number variation of enhancer hs1431 in humans.

GATA 4 Deletions Associated with Congenital Heart Diseases in South Brazil.

The normal development of the heart comprises a highly regulated machinery of genetic events, involving transcriptional factors. Congenital heart disease (CHD), have been associated with chromosomal abnormalities and copy number variants (CNVs). Our goal was to investigate through the multiplex ligation-dependent probe amplification (MLPA) technique, the presence of CNVs in reference genes for normal cardiac development in patients with CHD. GATA4 , NKX2-5 , TBX5 , BMP4 , and CRELD1 genes and 22q11.2 chromosome region were analyzed in 207 children with CHD admitted for the first time in a cardiac intensive care unit from a pediatric hospital. CNVs were detected in seven patients (3.4%): four had a 22q11.2 deletion (22q11DS) (1.9%), two had a GATA4 deletion (1%) and one had a 22q11.2 duplication (0.5%). No patients with CNVs in the NKX2-5 , TBX5 , BMP4 , and CRELD1 genes were identified. GATA4 deletions appear to be present in a significant number of CHD patients, especially those with septal defects, persistent left superior vena cava, pulmonary artery abnormalities, and extracardiac findings. GATA4 screening seems to be more effective when directed to these CHDs. The investigation of CNVs in GATA4 and 22q11 chromosome region in patients with CHD is important to anticipating the diagnosis, and to contributing to family planning.

Copy number variation burden does not predict severity of neurodevelopmental phenotype in children with a sex chromosome trisomy.

Sex chromosome trisomies (SCTs) (XXX, XXY, and XYY karyotypes) are associated with an elevated risk of neurodevelopmental disorders. The range of severity of the phenotype is substantial. We considered whether this variable outcome was related to the presence of copy number variants (CNVs)-stretches of duplicated or deleted DNA. A sample of 125 children with an SCT were compared with 181 children of normal karyotype who had been given the same assessments. First, we compared the groups on measures of overall CNV burden: number of CNVs, total span of CNVs, and likely functional impact (probability of loss-of-function intolerance, pLI, summed over CNVs). Differences between groups were small relative to within-group variance and not statistically significant on overall test. Next, we considered whether a measure of general neurodevelopmental impairment was predicted by pLI summed score, SCT versus comparison group, or the interaction between them. There was a substantial effect of SCT/comparison status but the pLI score was not predictive of outcomes in either group. We conclude that variable presence of CNVs is not a likely explanation for the wide phenotypic variation in children with SCTs. We discuss methodological challenges of testing whether CNVs are implicated in causing neurodevelopmental problems.

Copy Number Variation of the PIGY Gene in Sheep and Its Association Analysis with Growth Traits.

Simple SummaryThe PIGY (phosphatidylinositol glycan anchor biosynthesis class Y) gene is a member of the PIG gene family and encodes the glycosylphosphatidylinositol-N-acetylglucosaminyltransferase (GPI-GnT) complex. It initiates the biosynthesis of GPI and plays an important role in cell–cell interactions. Sequencing has revealed a 3600 bp copy number variation (CNV) in exon 2 of the PIGY gene in sheep, potentially altering a functional part of the protein. The CNV overlaps 28 quantitative trait loci that are relevant to some economic traits like muscle density and carcass weight. We screened for this CNV of the PIGY gene in 569 individuals, namely, 240 Chaka sheep (CKS), 168 Hu sheep (HS), and 161 small-tailed Han sheep (STHS), and analyzed the association between the presence of this CNV and sheep body size traits. The results showed that the loss-type CNV was more prevalent than other types in these three breeds, and there were significant effects of the PIGY gene CNV on body weight, chest circumference, and circumference of cannon bone of sheep. The results showed that sheep with gain-type CNV had better growth traits than those with other types. The findings reveal the relationship between the CNV of the PIGY gene and growth traits of sheep, suggesting that CNV could be utilized for improved molecular breeding of sheep.Copy number variation (CNV) is a type of genomic variation with an important effect on animal phenotype. We found that the PIGY gene contains a 3600 bp copy number variation (CNV) region located in chromosome 6 of sheep (Oar_v4.0 36,121,601–36,125,200 bp). This region overlaps with multiple quantitative trait loci related to phenotypes like muscle density and carcass weight. Therefore, in this study, the copy number variation of the PIGY gene was screened in three Chinese sheep breeds, namely, Chaka sheep (CKS, May of 2018, Wulan County, Qinghai Province, China), Hu sheep (HS, May of 2015, Mengjin County, Henan Province, China), and small-tailed Han sheep (STHS, May of 2016, Yongjing, Gansu Province, China). Association analyses were performed on the presence of CNV and sheep body size traits. We used real-time quantitative PCR (qPCR) to detect the CNV for association analysis. According to the results, the loss-type CNV was more common than other types in the three breeds (global average: loss = 61.5%, normal = 17.5%, and gain = 21.0%). The association analysis also showed significant effects of the PIGY gene CNV on body weight, chest circumference, and circumference of the cannon bone of sheep. Sheep with gain-type CNV had better growth traits than those with other types. The results indicate a clear relationship between the PIGY gene CNV and growth traits of sheep, suggesting the use of CNV as a new molecular breeding marker.

Copy number variation in human genomes from three major ethno-linguistic groups in Africa

BackgroundCopy number variation is an important class of genomic variation that has been reported in 75% of the human genome. However, it is underreported in African populations. Copy number variants (CNVs) could have important impacts on disease susceptibility and environmental adaptation. To describe CNVs and their possible impacts in Africans, we sequenced genomes of 232 individuals from three major African ethno-linguistic groups: (1) Niger Congo A from Guinea and Côte d’Ivoire, (2) Niger Congo B from Uganda and the Democratic Republic of Congo and (3) Nilo-Saharans from Uganda. We used GenomeSTRiP and cn.MOPS to identify copy number variant regions (CNVRs).ResultsWe detected 7608 CNVRs, of which 2172 were only deletions, 2384 were only insertions and 3052 had both. We detected 224 previously un-described CNVRs. The majority of novel CNVRs were present at low frequency and were not shared between populations. We tested for evidence of selection associated with CNVs and also for population structure. Signatures of selection identified previously, using SNPs from the same populations, were overrepresented in CNVRs. When CNVs were tagged with SNP haplotypes to identify SNPs that could predict the presence of CNVs, we identified haplotypes tagging 3096 CNVRs, 372 CNVRs had SNPs with evidence of selection (iHS > 3) and 222 CNVRs had both. This was more than expected (p < 0.0001) and included loci where CNVs have previously been associated with HIV, Rhesus D and preeclampsia. When integrated with 1000 Genomes CNV data, we replicated their observation of population stratification by continent but no clustering by populations within Africa, despite inclusion of Nilo-Saharans and Niger-Congo populations within our dataset.ConclusionsNovel CNVRs in the current study increase representation of African diversity in the database of genomic variants. Over-representation of CNVRs in SNP signatures of selection and an excess of SNPs that both tag CNVs and are subject to selection show that CNVs may be the actual targets of selection at some loci. However, unlike SNPs, CNVs alone do not resolve African ethno-linguistic groups. Tag haplotypes for CNVs identified may be useful in predicting African CNVs in future studies where only SNP data is available.

ENVIRONMENTAL AND GENETIC EVALUATION IN 46,XY DISORDERS OF SEX DEVELOPMENT

The phenotype of Disorders of Sex Development (DSD)patients depends on many factors including the presence of Copy Number Variation (CNVs) of different genes. The unbalanced rearrangement and the presence of deletions or duplications affect dramatically the phenotypic sex. The aim of this study is to correlate genotypic abnormalities with clinical phenotype in 46,XY DSD patients by Multiplex Ligation dependant Probe Amplification (MLPA) technique for accurate diagnosis in these patients and study possible paternal and maternal exposure to environmental risk factors in these cases. This study reported on forty patients with variable presentations of disorders of sex development (DSD) presenting with ambiguous genitalia, hypospadias, micropenis or with female phenotype with primary amenorrhea or short stature, with exclusion of cases with disorders in androgen synthesis. A complete personal, family history and clinical examination were done. Parents were asked to fill a questionnaire about frequency of dealing with some environmental factors including smoking, caffeine and using some materials having estrogenic effect as plastics, insecticides, and others. Conventional cytogenetics studies and fluorescence in situ hybridization (FISH) on peripheral blood as well as DNA extraction and Multiplex Ligation-dependent Probe Amplification (MLPA) were done for all cases for some genes; DMRT1, CYP17A1, SRD5A2, HSD17B3, DAX1, CXor21, SOX9, SRY, ZFY, WNT4 and SF. Results showed 46,XY was shown in all patients. MLPA analysis showed Copy Number Variation (CNVs) in 15% of cases. The study showed duplication of DMRT1 in 5% patients, deletion of SRY 2.5% patients, deletion of SOX9 in 5% patients, duplication of DAX1 in 12.5% patient, duplication of CYP17A41 in 5% patients, deletion of DMRT1 in 12.5% patient, duplication of SRD5A in 2.5% patient and duplication of HSD17B3 in 2.5% patient. Maternal exposure during pregnancy and paternal exposure to Endocrine Disrupting Compounds (EDCs) were insignificantly associated with DSD compared to control cases. This study demonstrates the importance of proper detection of genetic mutation in DSD patients showing a discrepancy between their karyotype and gonadal phenotype. It was concluded that using MLPA is recommended for better understanding of the phenotype with other recent techniques as for better diagnosis and follow up.

The clinical benefit of array-based comparative genomic hybridization for detection of copy number variants in Czech children with intellectual disability and developmental delay

BackgroundChromosomal microarray analysis has been shown to be a valuable and cost effective assay for elucidating copy number variants (CNVs) in children with intellectual disability and developmental delay (ID/DD).MethodsIn our study, we performed array-based comparative genomic hybridization (array-CGH) analysis using oligonucleotide-based platforms in 542 Czech patients with ID/DD, autism spectrum disorders and multiple congenital abnormalities. Prior to the array-CGH analysis, all the patients were first examined karyotypically using G-banding. The presence of CNVs and their putative derivation was confirmed using fluorescence in situ hybridization (FISH), multiplex ligation-dependent probe amplification (MLPA) and predominantly relative quantitative polymerase chain reaction (qPCR).ResultsIn total, 5.9% (32/542) patients were positive for karyotypic abnormalities. Pathogenic/likely pathogenic CNVs were identified in 17.7% of them (96/542), variants of uncertain significance (VOUS) were detected in 4.8% (26/542) and likely benign CNVs in 9.2% of cases (50/542). We identified 6.6% (36/542) patients with known recurrent microdeletion (24 cases) and microduplication (12 cases) syndromes, as well as 4.8% (26/542) patients with non-recurrent rare microdeletions (21 cases) and microduplications (5 cases). In the group of patients with submicroscopic pathogenic/likely pathogenic CNVs (13.3%; 68/510) we identified 91.2% (62/68) patients with one CNV, 5.9% (4/68) patients with two likely independent CNVs and 2.9% (2/68) patients with two CNVs resulting from cryptic unbalanced translocations. Of all detected CNVs, 21% (31/147) had a de novo origin, 51% (75/147) were inherited and 28% (41/147) of unknown origin.In our cohort pathogenic/likely pathogenic microdeletions were more frequent than microduplications (69%; 51/74 vs. 31%; 23/74) ranging in size from 0.395 Mb to 10.676 Mb (microdeletions) and 0.544 Mb to 8.156 Mb (microduplications), but their sizes were not significantly different (P = 0.83). The pathogenic/likely pathogenic CNVs (median 2.663 Mb) were significantly larger than benign CNVs (median 0.394 Mb) (P < 0.00001) and likewise the pathogenic/likely pathogenic CNVs (median 2.663 Mb) were significantly larger in size than VOUS (median 0.469 Mb) (P < 0.00001).ConclusionsOur results confirm the benefit of array-CGH in the current clinical genetic diagnostics leading to identification of the genetic cause of ID/DD in affected children.

Evaluation of SHOX defects in the era of next-generation sequencing.

Short stature homeobox (SHOX) haploinsufficiency is a frequent cause of short stature. Despite advances in sequencing technologies, the identification of SHOX mutations continues to be performed using standard methods, including multiplex ligation-dependent probe amplification (MLPA) followed by Sanger sequencing. We designed a targeted panel of genes associated with growth impairment, including SHOX genomic and enhancer regions, to improve the resolution of next-generation sequencing for SHOX analysis. We used two software packages, CONTRA and Nexus Copy Number, in addition to visual analysis to investigate the presence of copy number variants (CNVs). We evaluated 15 patients with previously known SHOX defects, including point mutations, deletions and a duplication, and 77 patients with idiopathic short stature (ISS). The panel was able to confirm all known defects in the validation analysis. During the prospective evaluation, we identified two new partial SHOX deletions (one detected only by visual analysis), including an intragenic deletion not detected by MLPA. Additionally, we were able to determine the breakpoints in four cases. Our results show that the designed panel can be used for the molecular investigation of patients with ISS, and it may even detect CNVs in SHOX and its enhancers, which may be present in a significant fraction of patients.

Molecular diagnosis and functional study of a pedigree affected with Lubs X-linked mental retardation syndrome

To explore the genetic basis for a pedigree affected with X-linked mental retardation. The proband was subjected to chromosomal karyotyping, FMR1 mutation testing and copy number variation analysis with a single nucleotide polymorphism microarray (SNP array). His family members were subjected to multiplex ligation-dependent probe amplification (MLPA) assaying. Expression of genes within the repeated region were analyzed. The proband had a normal chromosomal karyotype and normal number of CGG repeats within the FMR1 gene. SNP array identified a 370 kb duplication in Xq28 (ChrX: 153 027 633-153 398 515), which encompassed 14 genes including MECP2. The patient was diagnosed as Lubs X-linked mental retardation syndrome (MRXSL). MLPA confirmed the presence of copy number variation, its co-segregation with the disease, in addition with the carrier status of females. Genes from the duplicated region showed higher levels of expression (1.79 to 5.38 folds) within peripheral blood nucleated cells of the proband. The patients were diagnosed with MRXSL. The expression of affected genes was up-regulated due to the duplication. Genetic counseling and prenatal diagnosis may be provided based on the results.

Copy number variation in the MSRB3 gene enlarges porcine ear size through a mechanism involving miR-584-5p

BackgroundThe size and type of ears are important conformation characteristics that distinguish pig breeds. A significant quantitative trait locus (QTL) for ear size has been identified on SSC5 (SSC for Sus scrofa chromosome) but the underlying causative gene and mutation remain unknown. Thus, our aim was to identify the gene responsible for enlarged ears in pig.ResultsFirst, we narrowed down the QTL region on SSC5 to a 137.85-kb interval that harbors only the methionine sulfoxide reductase B3 (MSRB3) gene. Then, we identified a 38.7-kb copy number variation (CNV) that affects the last two exons of MSRB3 and could be the candidate causative mutation for this QTL. This CNV showed complete concordance with genotype at the QTL of the founder animals in a white Duroc × Erhualian F2 intercross and was found only in pigs from six Chinese indigenous breeds with large ears and from the Landrace breed with half-floppy ears. Moreover, it accounted for the significant association with ear size on SSC5 across the five pig populations tested. eQTL mapping revealed that this CNV was significantly associated with the expression of the microRNA (miRNA) miR-584-5p, which interacts with MSRB3, one of its target genes. In vivo and in vitro experiments confirmed that miR-584-5p inhibits the translation of MSRB3 mRNA. Taken together, these results led us to conclude that presence of the 38.7-kb CNV in the genome of some pig breeds affects ear size by altering the expression of miR-584-5p, which consequently hinders the expression of one of its target genes (e.g. MSRB3).ConclusionsOur findings shed insight into the underlying mechanism of development of external ears in mammals and contribute to a better understanding of how the presence of CNV can regulate gene expression.

Copy number variants in hypoplastic right heart syndrome.

Hypoplastic right heart syndrome (HRHS) is a rare congenital defect characterized by underdeveloped and malformed structures of the right heart. Familial recurrence of HRHS indicates genetic factors contribute to its etiology. Our study investigates the presence of copy number variants (CNVs) in HRHS cases. We genotyped 42 HRHS cases identified from live births throughout California (2003-2010) using the Illumina HumanOmni2.5-8 array. We identified 14 candidate CNVs in 14 HRHS cases (33%) based on the genes included in the CNVs and their functions. Duplications overlapping part of ERBB4 were identified in two unrelated cases. ERBB4 is a neuregulin receptor with a pivotal role in cardiomyocyte differentiation and heart development. We also described a 7.5 Mb duplication at 16q11-12. Multiple genes in the duplicated region have previously been linked to heart defects and cardiac development, including RPGRIP1L, RBL2, SALL1, and MYLK3. Of the 14 validated CNVs, we identified four CNVs in close proximity to genes linked to the Wnt signaling pathway. This study expands on our previous work supporting the role of genetics in HRHS. We identified CNVs affecting crucial genes and signaling pathways involved in right heart development. ERBB4 and duplication of the 16q11-12 region are important areas for future investigation.

High Incidence of Copy Number Variants in Adults with Intellectual Disability and Co-morbid Psychiatric Disorders

A genetic analysis of unexplained mild-moderate intellectual disability and co-morbid psychiatric or behavioural disorders is not systematically conducted in adults. A cohort of 100 adult patients affected by both phenotypes were analysed in order to identify the presence of copy number variants (CNVs) responsible for their condition identifying a yield of 12.8% of pathogenic CNVs (19% when including clinically recognizable microdeletion syndromes). Moreover, there is a detailed clinical description of an additional 11% of the patients harbouring possible pathogenic CNVs—including a 7q31 deletion (IMMP2L) in two unrelated patients and duplications in 3q29, 9p24.2p24.1 and 15q14q15.1—providing new evidence of its contribution to the phenotype. This study adds further proof of including chromosomal microarray analysis (CMA) as a mandatory test to improve the diagnosis in the adult patients in psychiatric services.

Polymorphisms in the SOX9 region and testicular disorder of sex development (38,XX; SRY-negative) in pigs

Testicular XX disorder of sex development (XX DSD, SRY-negative), causing sterility, is quite frequently diagnosed in pig populations, however, its molecular background is still unknown. This disorder may affect carcass quality (boar taint) due to the presence of testicular tissue in animals with ambiguous female external genitalia. A chromosome fragment encompassing the SOX9 gene was previously considered as a candidate region in functional and association studies. We analyzed distribution of polymorphic variants in 5’UTR and in 3’-flanking regions of the SOX9 gene in a cohort of 12 XX DSD pigs and a panel of 116 normal females. Altogether, 8 previously known polymorphic sites were analyzed and no significant association under Bonferroni correction (P>0.0063) between DSD phenotype and identified SNPs was found. The region harboring the SOX9 gene was also searched for the presence of copy number variation (CNV) by fluorescence in situ hybridization technique (FISH), using a set of 7 BAC probes. A potential CNV region, manifested by size variation (large and small) of fluorescence signals, produced by a single BAC probe hybridizing downstream (approx. 500kb) of the SOX9 gene, was observed in 4 DSD cases. We conclude, that a new CNV polymorphism in the region harboring SOX9 gene can be considered as a promising marker for the DSD phenotype.

Emergence of a Viral RNA Polymerase Variant during Gene Copy Number Amplification Promotes Rapid Evolution of Vaccinia Virus.

Viruses can evolve quickly to defeat host immune functions. For poxviruses, little is known about how multiple adaptive mutations emerge in populations at the same time. In this study, we uncovered a means of vaccinia virus adaptation involving the accumulation of distinct genetic variants within a single population. We identified adaptive point mutations in the viral RNA polymerase gene A24R and, surprisingly, found that one of these mutations activates the nucleic acid sensing factor PKR. We also found that gene copy number variation (CNV) can provide dual benefits to evolving virus populations, including evidence that CNV facilitates the accumulation of a point mutation distant from the expanded locus. Our data suggest that transient CNV can accelerate the fixation of mutations conferring modest benefits, or even fitness trade-offs, and highlight how structural variation might aid poxvirus adaptation through both direct and indirect actions.

Family-assisted inference of the genetic architecture of major histocompatibility complex variation.

With their direct link to individual fitness, genes of the major histocompatibility complex (MHC) are a popular system to study the evolution of adaptive genetic diversity. However, owing to the highly dynamic evolution of the MHC region, the isolation, characterization and genotyping of MHC genes remain a major challenge. While high-throughput sequencing technologies now provide unprecedented resolution of the high allelic diversity observed at the MHC, in many species, it remains unclear (i) how alleles are distributed among MHC loci, (ii) whether MHC loci are linked or segregate independently and (iii) how much copy number variation (CNV) can be observed for MHC genes in natural populations. Here, we show that the study of allele segregation patterns within families can provide significant insights in this context. We sequenced two MHC class I (MHC-I) loci in 1267 European barn owls (Tyto alba), including 590 offspring from 130 families using Illumina MiSeq technology. Coupled with a high per-individual sequencing coverage (~3000×), the study of allele segregation patterns within families provided information on three aspects of the architecture of MHC-I variation in barn owls: (i) extensive sharing of alleles among loci, (ii) strong linkage of MHC-I loci indicating tandem architecture and (iii) the presence of CNV in the barn owl MHC-I. We conclude that the additional information that can be gained from high-coverage amplicon sequencing by investigating allele segregation patterns in families not only helps improving the accuracy of MHC genotyping, but also contributes towards enhanced analyses in the context of MHC evolutionary ecology.

The Intolerance of Regulatory Sequence to Genetic Variation Predicts Gene Dosage Sensitivity.

Noncoding sequence contains pathogenic mutations. Yet, compared with mutations in protein-coding sequence, pathogenic regulatory mutations are notoriously difficult to recognize. Most fundamentally, we are not yet adept at recognizing the sequence stretches in the human genome that are most important in regulating the expression of genes. For this reason, it is difficult to apply to the regulatory regions the same kinds of analytical paradigms that are being successfully applied to identify mutations among protein-coding regions that influence risk. To determine whether dosage sensitive genes have distinct patterns among their noncoding sequence, we present two primary approaches that focus solely on a gene’s proximal noncoding regulatory sequence. The first approach is a regulatory sequence analogue of the recently introduced residual variation intolerance score (RVIS), termed noncoding RVIS, or ncRVIS. The ncRVIS compares observed and predicted levels of standing variation in the regulatory sequence of human genes. The second approach, termed ncGERP, reflects the phylogenetic conservation of a gene’s regulatory sequence using GERP++. We assess how well these two approaches correlate with four gene lists that use different ways to identify genes known or likely to cause disease through changes in expression: 1) genes that are known to cause disease through haploinsufficiency, 2) genes curated as dosage sensitive in ClinGen’s Genome Dosage Map, 3) genes judged likely to be under purifying selection for mutations that change expression levels because they are statistically depleted of loss-of-function variants in the general population, and 4) genes judged unlikely to cause disease based on the presence of copy number variants in the general population. We find that both noncoding scores are highly predictive of dosage sensitivity using any of these criteria. In a similar way to ncGERP, we assess two ensemble-based predictors of regional noncoding importance, ncCADD and ncGWAVA, and find both scores are significantly predictive of human dosage sensitive genes and appear to carry information beyond conservation, as assessed by ncGERP. These results highlight that the intolerance of noncoding sequence stretches in the human genome can provide a critical complementary tool to other genome annotation approaches to help identify the parts of the human genome increasingly likely to harbor mutations that influence risk of disease.