Association of <Emphasis Type="Italic">TSHR</Emphasis> Gene Copy Number Variation with TSH Abnormalities

Thyroid-stimulating hormone (TSH) is secreted by the pituitary gland and promotes thyroid growth and function, with increased TSH levels typically associated with hypothyroidism. By consulting the literature, we found that the TSHR, PAX8, and PDE4B genes are associated with thyroid function. Recently, copy number variations (CNVs) have been used as genetic markers to investigate inter-individual variation. Therefore, we investigated the relationship between the TSHR, PAX8, and PDE4B gene CNVs and TSH abnormalities, by calculating variations in gene copy number. Four hundred and eighty-one participants, 232 healthy controls and 249 patients with TSH abnormalities, were selected from three distinct areas in China with different iodine statuses. RT-PCR was used to detect CNVs. Urinary iodine concentrations (UIC) were measured by As³⁺–Ce⁴⁺ catalytic spectrophotometry. There was an association between a CNV at the TSHR gene and TSH abnormalities (p?=?0.002). The distribution of PAX8 and PDE4B gene CNVs between patients with TSH abnormalities and healthy controls was not significantly different. UIC >?200 μg/l (OR?=?1.49, 95% CI?=?1.01–2.22) and the TSHR gene (OR?=?6.01, 95% CI?=?1.96–18.41) were found to be risk factors for TSH abnormalities. PAX8 and PDE4B gene CNVs were not significantly associated with TSH abnormalities. There was no significant interaction between UIC and any of the examined CNVs. In conclusion, the TSHR gene CNV was associated with the development of TSH abnormalities. No significant associations were revealed between urinary iodine levels and candidate gene CNVs.     

Integrating CNVs into meta-QTL identified <Emphasis Type="Italic">GBP4</Emphasis> as positional candidate for adult cattle stature

Copy number variation (CNV) of DNA sequences, functionally significant but yet fully ascertained, is believed to confer considerable increments in unexplained heritability of quantitative traits. Identification of phenotype-associated CNVs (paCNVs) therefore is a pressing need in CNV studies to speed up their exploitation in cattle breeding programs. Here, we provided a new avenue to achieve this goal that is to project the published CNV data onto meta-quantitative trait loci (meta-QTL) map which connects causal genes with phenotypes. Any CNVs overlapping meta-QTL therefore will be potential paCNVs. This study reported potential paCNVs in Bos taurus autosome 3 (BTA3). Notably, overview indexes and CNVs both highlighted a narrower region (BTA3 54,500,000–55,000,000 bp, named BTA3_INQTL_6) within one constructed meta-QTL. Then, we ascertained guanylate-binding protein 4 (GBP4) among the nine positional candidate genes was significantly associated with adult cattle stature, including body weight (BW, P?<?0.05) and withers height (WHT, P?<?0.05), fitting GBP4 CNV either with three levels or with six levels in the model. Although higher copy number downregulated the mRNA levels of GBP2 (P?<?0.05) and GBP4 (P?<?0.05) in 1-Mb window (54.0–55.0 Mb) in muscle and adipose, additional analyses will be needed to clarify the causality behind the ascertained association.     

Number of rare germline CNVs and <Emphasis Type="Italic">TP53</Emphasis> mutation types

Background

The Li-Fraumeni syndrome (LFS), an inherited rare cancer predisposition syndrome characterized by a variety of early-onset tumors, is caused by different highly penetrant germline mutations in the TP53 gene; each separate mutation has dissimilar functional and phenotypic effects, which partially clarifies the reported heterogeneity between LFS families. Increases in copy number variation (CNV) have been reported in TP53 mutated individuals, and are also postulated to contribute to LFS phenotypic variability. The Brazilian p.R337H TP53 mutation has particular functional and regulatory properties that differ from most other common LFS TP53 mutations, by conferring a strikingly milder phenotype.

Methods

We compared the CNV profiles of controls, and LFS individuals carrying either p.R337H or DNA binding domain (DBD) TP53 mutations by high resolution array-CGH.

Results

Although we did not find any significant difference in the frequency of CNVs between LFS patients and controls, our data indicated an increased proportion of rare CNVs per genome in patients carrying DBD mutations compared to both controls (p=0.0002***) and p.R337H (0.0156*) mutants.

Conclusions

The larger accumulation of rare CNVs in DBD mutants may contribute to the reported anticipation and severity of the syndrome; likewise the fact that p.R337H individuals do not present the same magnitude of rare CNV accumulation may also explain the maintenance of this mutation at relatively high frequency in some populations.

     

Diagnostic efficacy and new variants in isolated and complex autism spectrum disorder using molecular karyotyping

Autism spectrum disorder (ASD) is a group of the neurodevelopment disorders presenting as an isolated ASD or more complex forms, where a broader clinical phenotype comprised of developmental delay and intellectual disability is present. Both the isolated and complex forms have a significant causal genetic component and submicroscopic genomic copy number variations (CNV) are the most common identifiable genetic factor in these patients. The data on microarray testing in ASD cohorts are still accumulating and novel loci are often identified; therefore, we aimed to evaluate the diagnostic efficacy of the method and the relevance of implementing it into routine genetic testing in ASD patients. A genome-wide CNV analysis using the Agilent microarrays was performed in a group of 150 individuals with an isolated or complex ASD. Altogether, 11 (7.3%) pathogenic CNVs and 15 (10.0%) variants of unknown significance (VOUS) were identified, with the highest proportion of pathogenic CNVs in the subgroup of the complex ASD patients (14.3%). An interesting case of previously unreported partial UPF3B gene deletion was identified among the pathogenic CNVs. Among the CNVs with unknown significance, four VOUS involved genes with possible correlation to ASD, namely genes SNTG2, PARK2, CADPS2 and NLGN4X. The diagnostic efficacy of aCGH in our cohort was comparable with those of the previously reported and identified an important proportion of genetic ASD cases. Despite the continuum of published studies on the CNV testing in ASD cohorts, a considerable number of VOUS CNVs is still being identified, namely 10.0% in our study.     

Analyses of copy number variation reveal putative susceptibility loci in MTX‐induced mouse neural tube defects

Copy number variations (CNVs) are thought to act as an important genetic mechanism underlying phenotypic heterogeneity. Impaired folate metabolism can result in neural tube defects (NTDs). However, the precise nature of the relationship between low folate status and NTDs remains unclear. Using an array‐comparative genomic hybridization (aCGH) assay, we investigated whether CNVs could be detected in the NTD embryonic neural tissues of methotrexate (MTX)‐induced folate dysmetabolism pregnant C57BL/6 mice and confirmed the findings with quantitative real‐time PCR (qPCR). The CNVs were then comprehensively investigated using bioinformatics methods to prioritize candidate genes. We measured dihydrofolate reductase (DHFR) activity and concentrations of folate and relevant metabolites in maternal serum using enzymologic method and liquid chromatography/tandem mass spectrometry (LC/MS/MS). Three high confidence CNVs on XqA1.1, XqA1.1‐qA2, and XqE3 were found in the NTD embryonic neural tissues. Twelve putative genes and three microRNAs were identified as potential susceptibility candidates in MTX‐induced NTDs and possible roles in NTD pathogenesis. DHFR activity and 5‐methyltetrahydrofolate (5‐MeTHF), 5‐formyltetrahydrofolate (5‐FoTHF), and S‐adenosylmethionine (SAM) concentrations of maternal serum decreased significantly after MTX injection. These findings suggest that CNVs caused by defects in folate metabolism lead to NTD, and further support the hypothesis that folate dysmetabolism is a direct cause for CNVs in MTX‐induced NTDs. © 2014 Wiley Periodicals, Inc. Develop Neurobiol 74: 877–893, 2014     

Association between copy-number variation on metabolic phenotypes and HDL-C levels in patients with polycystic ovary syndrome

Polygenic diseases with a broad phenotypic spectrum, such as polycystic ovary syndrome (PCOS), present a particular challenge in terms of identifying the underlying genetic mechanisms, nevertheless genetic variants have impact on the individual phenotype. We aimed to determine if next to genetic variations like SNPs further mechanisms might play a role in the pathogenesis of PCOS. We examined the effect of copy-number variations (CNVs) on metabolic phenotypes in PCOS. The intragenic rs1244979, rs2815752 in NEGR1 gene, and rs780094 in GCKR gene were genotyped and CNVs were determined by droplet digital polymerase chain reaction (ddPCR) in PCOS patients (n?=?153) and controls without metabolic syndrome (n?=?142). The study indicated that SNPs are not associated with the pathogenesis of PCOS but affect metabolic phenotypes. The CNVs investigated show a lower variability in PCOS than in CON. Furthermore, we provided direct evidence that the copy number, but not the genotype of the CNV in the genomic regions of rs780094(GCKR) is associated with low level of high-density lipoprotein cholesterol in PCOS. This study supports the hypothesis that not only genetic variants, but also CNVs in metabolically relevant genes, have an effect on metabolic phenotypes in our group of PCOS patients.     

Y-chromosome polymorphisms of the domestic Bactrian camel in China

Single-nucleotide polymorphisms (SNPs), microsatellites and copy number variation (CNV) were studied on the Y chromosome to understand the paternal origin and phylogenetic relationships for resource protection, rational development and utilization of the domestic Bactrian camel in China. Our sample set consisted of 94 Chinese domestic Bactrian camels from four regions (Inner Mongolia, Gansu, Qinghai and Xinjiang), we screened 29 Y-chromosome-specific loci for SNPs, analysed 40 bovine-derived microsatellite loci and measured CNVs of HSFY and SRY through Sanger sequencing, automated fluorescence-based microsatellite analysis and quantitative real-time PCR, respectively. A multicopy gene, SRY, was first found, and sequence variation was only detected in SRY in a screen of 29 loci in 13 DNA pools of individual camels. In addition, a TG repeat in the USP9Y gene was identified as the first polymorphic microsatellite in the camel Y chromosome, whereas microsatellite based on bovine sequences were not detected. The frequency of each allele varied among different populations. For the Nanjiang, Hexi and Alashan populations, a 243-bp allele was found. For the Sunite population, 241-bp, 243-bp and 247-bp alleles were detected, and the frequencies of these alleles were \(22.2\%\), \(44.5\%\) and \(33.3\%\), respectively; 241-bp and 243-bp alleles were found in other populations. Finally, CNVs in two Y-chromosomal genes were detected; CNV for HSFY and SRY ranged from 1 to 3 and from 1 to 9, respectively.     

Copy Number Variations in Tilapia Genomes

Discovering the nature and pattern of genome variation is fundamental in understanding phenotypic diversity among populations. Although several millions of single nucleotide polymorphisms (SNPs) have been discovered in tilapia, the genome-wide characterization of larger structural variants, such as copy number variation (CNV) regions has not been carried out yet. We conducted a genome-wide scan for CNVs in 47 individuals from three tilapia populations. Based on 254 Gb of high-quality paired-end sequencing reads, we identified 4642 distinct high-confidence CNVs. These CNVs account for 1.9% (12.411 Mb) of the used Nile tilapia reference genome. A total of 1100 predicted CNVs were found overlapping with exon regions of protein genes. Further association analysis based on linear model regression found 85 CNVs ranging between 300 and 27,000 base pairs significantly associated to population types (R ² > 0.9 and P > 0.001). Our study sheds first insights on genome-wide CNVs in tilapia. These CNVs among and within tilapia populations may have functional effects on phenotypes and specific adaptation to particular environments.     

The Role of Copy Number Variation in Susceptibility to Amyotrophic Lateral Sclerosis: Genome-Wide Association Study and Comparison with Published Loci

Background

The genetic contribution to sporadic amyotrophic lateral sclerosis (ALS) has not been fully elucidated. There are increasing efforts to characterise the role of copy number variants (CNVs) in human diseases; two previous studies concluded that CNVs may influence risk of sporadic ALS, with multiple rare CNVs more important than common CNVs. A little-explored issue surrounding genome-wide CNV association studies is that of post-calling filtering and merging of raw CNV calls. We undertook simulations to define filter thresholds and considered optimal ways of merging overlapping CNV calls for association testing, taking into consideration possibly overlapping or nested, but distinct, CNVs and boundary estimation uncertainty.

Methodology and Principal Findings

In this study we screened Illumina 300K SNP genotyping data from 730 ALS cases and 789 controls for copy number variation. Following quality control filters using thresholds defined by simulation, a total of 11321 CNV calls were made across 575 cases and 621 controls. Using region-based and gene-based association analyses, we identified several loci showing nominally significant association. However, the choice of criteria for combining calls for association testing has an impact on the ranking of the results by their significance. Several loci which were previously reported as being associated with ALS were identified here. However, of another 15 genes previously reported as exhibiting ALS-specific copy number variation, only four exhibited copy number variation in this study. Potentially interesting novel loci, including EEF1D, a translation elongation factor involved in the delivery of aminoacyl tRNAs to the ribosome (a process which has previously been implicated in genetic studies of spinal muscular atrophy) were identified but must be treated with caution due to concerns surrounding genomic location and platform suitability.

Conclusions and Significance

Interpretation of CNV association findings must take into account the effects of filtering and combining CNV calls when based on early genome-wide genotyping platforms and modest study sizes.     

A novel Alu-mediated microdeletion in the <Emphasis Type="BoldItalic">RUNX2</Emphasis> gene in a Chinese patient with cleidocranial dysplasia

Cleidocranial dysplasia (CCD; OMIM: 119600) is a rare autosomal dominant skeletal dysplasia caused by RUNX2 gene mutations. The present study described a sporadic case with CCD. The clinical data of the proband with CCD was reported and genetic analysis was performed. The proband presented with typical CCD features including supernumerary impacted teeth, bilateral clavicle dysplasia, delayed closure of cranial sutures, and short stature; while his hands were normal. Sequencing analysis of the entire coding region of the RUNX2 gene revealed no pathogenic changes; however, copy-number analysis with the Affymetrix HD array found \(\sim \)500 kb genomic microdeletion. Real-time quantitative PCR validated this microdeletion in the 1–4 exons of the RUNX2 gene. The junction point of the breaking DNA was located in the directly oriented AluSz6 and AluSx repetitive elements, indicating that this microdeletion might be generated through an Alu–Alu mediated mechanism. In addition, this microdeletion existed in 21.8% of the asymptomatic mother’s peripheral blood cells, demonstrating that the mosaicism was not associated with CCD phenotypes. In summary, a pathogenic microdeletion in the RUNX2 gene located on chromosome 6 was responsible for CCD.     

A genome‐wide association study of copy number variations with umbilical hernia in swine

Umbilical hernia (UH) is one of the most common congenital defects in pigs, leading to considerable economic loss and serious animal welfare problems. To test whether copy number variations (CNVs) contribute to pig UH, we performed a case–control genome‐wide CNV association study on 905 pigs from the Duroc, Landrace and Yorkshire breeds using the Porcine SNP60 BeadChip and penncnv algorithm. We first constructed a genomic map comprising 6193 CNVs that pertain to 737 CNV regions. Then, we identified eight CNVs significantly associated with the risk for UH in the three pig breeds. Six of seven significantly associated CNVs were validated using quantitative real‐time PCR. Notably, a rare CNV (CNV14:13030843–13059455) encompassing the NUGGC gene was strongly associated with UH (permutation‐corrected P = 0.0015) in Duroc pigs. This CNV occurred exclusively in seven Duroc UH‐affected individuals. SNPs surrounding the CNV did not show association signals, indicating that rare CNVs may play an important role in complex pig diseases such as UH. The NUGGC gene has been implicated in human omphalocele and inguinal hernia. Our finding supports that CNVs, including the NUGGC CNV, contribute to the pathogenesis of pig UH.     

Analysis of chromosomal structural variation in patients with congenital left‐sided cardiac lesions

Background: We sought to characterize the landscape of structural variation associated with the subset of congenital cardiac defects characterized by left‐sided obstruction. Methods: Cases with left‐sided cardiac defects (LSCD) and pediatric controls were uniformly genotyped and assessed for copy number variant (CNV) calls. Significance testing was performed to ascertain differences in overall CNV incidence, and for CNV enrichment of specific genes and gene functions in LSCD cases relative to controls. Results: A total of 257 cases of European descent and 962 ethnically matched, disease‐free pediatric controls were included. Although there was no difference in CNV rate between cases and controls, a significant enrichment in rare LSCD CNVs was detected overall (p = 7.30 × 10^?3, case/control ratio = 1.26) and when restricted either to deletions (p = 7.58 × 10^?3, case/control ratio = 1.20) or duplications (3.02 × 10^?3, case/control ratio = 1.43). Neither gene‐based, functional nor knowledge‐based analyses identified genes, loci or pathways that were significantly enriched in cases as compared to controls when appropriate corrections for multiple tests were applied. However, several genes of interest were identified by virtue of their association with cardiac development, known human conditions, or reported disruption by CNVs in other patient cohorts. Conclusion: This study examines the largest cohort to date with LSCD for structural variation. These data suggest that CNVs play a role in disease risk and identify numerous genes disrupted by CNVs of potential disease relevance. These findings further highlight the genetic heterogeneity and complexity of these disorders. Birth Defects Research (Part A) 100:951–964, 2014. © 2014 Wiley Periodicals, Inc.     

Clinical analysis of germline copy number variation in DMD using a non-conjugate hierarchical Bayesian model

Background

Detection of copy number variants (CNVs) is an important aspect of clinical testing for several disorders, including Duchenne muscular dystrophy, and is often performed using multiplex ligation-dependent probe amplification (MLPA). However, since many genetic carrier screens depend instead on next-generation sequencing (NGS) for wider discovery of small variants, they often do not include CNV analysis. Moreover, most computational techniques developed to detect CNVs from exome sequencing data are not suitable for carrier screening, as they require matched normals, very large cohorts, or extensive gene panels.

Methods

We present a computational software package, geneCNV (http://github.com/vkozareva/geneCNV), which can identify exon-level CNVs using exome sequencing data from only a few genes. The tool relies on a hierarchical parametric model trained on a small cohort of reference samples.

Results

Using geneCNV, we accurately inferred heterozygous CNVs in the DMD gene across a cohort of 15 test subjects. These results were validated against MLPA, the current standard for clinical CNV analysis in DMD. We also benchmarked the tool’s performance against other computational techniques and found comparable or improved CNV detection in DMD using data from panels ranging from 4,000 genes to as few as 8 genes.

Conclusions

geneCNV allows for the creation of cost-effective screening panels by allowing NGS sequencing approaches to generate results equivalent to bespoke genotyping assays like MLPA. By using a parametric model to detect CNVs, it also fulfills regulatory requirements to define a reference range for a genetic test. It is freely available and can be incorporated into any Illumina sequencing pipeline to create clinical assays for detection of exon duplications and deletions.

     

Precision-mapping and statistical validation of quantitative trait loci by machine learning

Background

DNA sequence diversity within the human genome may be more greatly affected by copy number variations (CNVs) than single nucleotide polymorphisms (SNPs). Although the importance of CNVs in genome wide association studies (GWAS) is becoming widely accepted, the optimal methods for identifying these variants are still under evaluation. We have previously reported a comprehensive view of CNVs in the HapMap DNA collection using high density 500 K EA (Early Access) SNP genotyping arrays which revealed greater than 1,000 CNVs ranging in size from 1 kb to over 3 Mb. Although the arrays used most commonly for GWAS predominantly interrogate SNPs, CNV identification and detection does not necessarily require the use of DNA probes centered on polymorphic nucleotides and may even be hindered by the dependence on a successful SNP genotyping assay.

Results

In this study, we have designed and evaluated a high density array predicated on the use of non-polymorphic oligonucleotide probes for CNV detection. This approach effectively uncouples copy number detection from SNP genotyping and thus has the potential to significantly improve probe coverage for genome-wide CNV identification. This array, in conjunction with PCR-based, complexity-reduced DNA target, queries over 1.3 M independent NspI restriction enzyme fragments in the 200 bp to 1100 bp size range, which is a several fold increase in marker density as compared to the 500 K EA array. In addition, a novel algorithm was developed and validated to extract CNV regions and boundaries.

Conclusion

Using a well-characterized pair of DNA samples, close to 200 CNVs were identified, of which nearly 50% appear novel yet were independently validated using quantitative PCR. The results indicate that non-polymorphic probes provide a robust approach for CNV identification, and the increasing precision of CNV boundary delineation should allow a more complete analysis of their genomic organization.     

Copy-number variant analysis of classic heterotaxy highlights the importance of body patterning pathways

Classic heterotaxy consists of congenital heart defects with abnormally positioned thoracic and abdominal organs. We aimed to uncover novel, genomic copy-number variants (CNVs) in classic heterotaxy cases. A microarray containing 2.5 million single-nucleotide polymorphisms (SNPs) was used to genotype 69 infants (cases) with classic heterotaxy identified from California live births from 1998 to 2009. CNVs were identified using the PennCNV software. We identified 56 rare CNVs encompassing genes in the NODAL (NIPBL, TBX6), BMP (PPP4C), and WNT (FZD3) signaling pathways, not previously linked to classic heterotaxy. We also identified a CNV involving FGF12, a gene previously noted in a classic heterotaxy case. CNVs involving RBFOX1 and near MIR302F were detected in multiple cases. Our findings illustrate the importance of body patterning pathways for cardiac development and left/right axes determination. FGF12, RBFOX1, and MIR302F could be important in human heterotaxy, because they were noted in multiple cases. Further investigation into genes involved in the NODAL, BMP, and WNT body patterning pathways and into the dosage effects of FGF12, RBFOX1, and MIR302F is warranted.     

Genome-Wide Analysis of Copy Number Variation Identifies Candidate Gene Loci Associated with the Progression of Non-Alcoholic Fatty Liver Disease

Between 10 and 25% of individuals with non-alcoholic fatty liver disease (NAFLD) develop hepatic fibrosis leading to cirrhosis and hepatocellular carcinoma (HCC). To investigate the molecular basis of disease progression, we performed a genome-wide analysis of copy number variation (CNV) in a total of 49 patients with NAFLD [10 simple steatosis and 39 non-alcoholic steatohepatitis (NASH)] and 49 matched controls using high-density comparative genomic hybridization (CGH) microarrays. A total of 11 CNVs were found to be unique to individuals with simple steatosis, whilst 22 were common between simple steatosis and NASH, and 224 were unique to NASH. We postulated that these CNVs could be involved in the pathogenesis of NAFLD progression. After stringent filtering, we identified four rare and/or novel CNVs that may influence the pathogenesis of NASH. Two of these CNVs, located at 13q12.11 and 12q13.2 respectively, harbour the exportin 4 (XPO4) and phosphodiesterase 1B (PDE1B) genes which are already known to be involved in the etiology of liver cirrhosis and HCC. Cross-comparison of the genes located at these four CNV loci with genes already known to be associated with NAFLD yielded a set of genes associated with shared biological processes including cell death, the key process involved in ‘second hit’ hepatic injury. To our knowledge, this pilot study is the first to provide CNV information of potential relevance to the NAFLD spectrum. These data could prove invaluable in predicting patients at risk of developing NAFLD and more importantly, those who will subsequently progress to NASH.     

Pathogenic copy number variations are associated with foetal short femur length in a tertiary referral centre study

Shortened foetal femur length (FL) is a common abnormal phenotype that often causes anxiety in pregnant women, and standard clinical treatments remain unavailable. We investigated the clinical characteristics, genetic aetiology and obstetric pregnancy outcomes of foetuses with short FL and provided a reference for perinatal management of such cases. Chromosomal microarray analysis was used to analyse the copy number variations (CNV) in short FL foetuses. Of the 218 foetuses with short FL, 33 foetuses exhibited abnormal CNVs, including 19 with pathogenic CNVs and 14 with variations of uncertain clinical significance. Of the 19 foetuses with pathogenic CNVs, four had aneuploidy, 14 had deletions/duplications, and one had pathogenic uniparental diploidy. The 7q11.23 microdeletion was detected in three foetuses. The severity of short FL was not associated with the rate of pathogenic CNVs. The duration of short FL for the intrauterine ultrasound phenotype in foetuses carrying a pathogenic CNV was independent of the gestational age. Further, maternal age was not associated with the incidence of foetal pathogenic CNVs. Adverse pregnancy outcomes occurred in 77 cases, including termination of pregnancy in 63 cases, postnatal dwarfed foetuses with intellectual disability in 11 cases, and three deaths within 3 months of birth. Pathogenic CNVs closely related to foetal short FL were identified, among which the 7q11.23 microdeletion was highly associated with short FL development. This study provides a reference for the perinatal management of foetuses with short FL.