New copy number variations in schizophrenia

Genome-wide screenings for copy number variations (CNVs) in patients with schizophrenia have demonstrated the presence of several CNVs that increase the risk of developing the disease and a growing number of large rare CNVs; the contribution of these rare CNVs to schizophrenia remains unknown. Using Affymetrix 6.0 arrays, we undertook a systematic search for CNVs in 172 patients with schizophrenia and 160 healthy controls, all of Italian origin, with the aim of confirming previously identified loci and identifying novel schizophrenia susceptibility genes. We found five patients with a CNV occurring in one of the regions most convincingly implicated as risk factors for schizophrenia: NRXN1 and the 16p13.1 regions were found to be deleted in single patients and 15q11.2 in 2 patients, whereas the 15q13.3 region was duplicated in one patient. Furthermore, we found three distinct patients with CNVs in 2q12.2, 3q29 and 17p12 loci, respectively. These loci were previously reported to be deleted or duplicated in patients with schizophrenia but were never formally associated with the disease. We found 5 large CNVs (>900 kb) in 4q32, 5q14.3, 8q23.3, 11q25 and 17q12 in five different patients that could include some new candidate schizophrenia susceptibility genes. In conclusion, the identification of previously reported CNVs and of new, rare, large CNVs further supports a model of schizophrenia that includes the effect of multiple, rare, highly penetrant variants.     

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.     

Investigation of 15q11-q13, 16p11.2 and 22q13 CNVs in Autism Spectrum Disorder Brazilian Individuals with and without Epilepsy

Copy number variations (CNVs) are an important cause of ASD and those located at 15q11-q13, 16p11.2 and 22q13 have been reported as the most frequent. These CNVs exhibit variable clinical expressivity and those at 15q11-q13 and 16p11.2 also show incomplete penetrance. In the present work, through multiplex ligation-dependent probe amplification (MLPA) analysis of 531 ethnically admixed ASD-affected Brazilian individuals, we found that the combined prevalence of the 15q11-q13, 16p11.2 and 22q13 CNVs is 2.1% (11/531). Parental origin could be determined in 8 of the affected individuals, and revealed that 4 of the CNVs represent de novo events. Based on CNV prediction analysis from genome-wide SNP arrays, the size of those CNVs ranged from 206 kb to 2.27 Mb and those at 15q11-q13 were limited to the 15q13.3 region. In addition, this analysis also revealed 6 additional CNVs in 5 out of 11 affected individuals. Finally, we observed that the combined prevalence of CNVs at 15q13.3 and 22q13 in ASD-affected individuals with epilepsy (6.4%) was higher than that in ASD-affected individuals without epilepsy (1.3%; p<0.014). Therefore, our data show that the prevalence of CNVs at 15q13.3, 16p11.2 and 22q13 in Brazilian ASD-affected individuals is comparable to that estimated for ASD-affected individuals of pure or predominant European ancestry. Also, it suggests that the likelihood of a greater number of positive MLPA results might be found for the 15q13.3 and 22q13 regions by prioritizing ASD-affected individuals with epilepsy.     

Genome-wide profiling of gene amplification and deletion in cancer

Accumulations of genetic changes in somatic cells induce phenotypic transformations leading to cancer. Among these genetic changes, gene amplification and deletion are most frequently observed in several kinds of cancers. Amplification of oncogene and/or deletion of tumor suppressor gene, together with dysfunction of the gene by point mutation, are the main causes of cancer. Genome-wide analysis of amplification and deletion of genes in cancers is basic to resolving the mechanisms of carcinogenesis. Comparative genomic hybridization (CGH) developed in 1992 has been utilized to identify DNA copy number abnormalities in various kind of cancers and several reports have shown its usefulness in screening of the genes involved in carcinogenesis, and also in the identification of prognostic factors in cancer. We have shown that 1q23 gain is associated with neuroblastomas that are resistant to aggressive treatment, and have poor prognosis, and 1q and 13q gains are possibly related to drug resistance in ovarian cancers. Recently, the "rough draft" of the human genome was reported and we are ready to utilize the vast information on genomic sequences in cancer research. Moreover, microarray technology enables us to analyze more than ten thousand genes at a time and revealed genetic abnormalities in cancers at a genome-wide level. By combination of microarray and CGH, a powerful screening method for oncogenes and tumor suppressor genes in cancers, called array-CGH, has been developed by several groups. In this article, we overview these genome-wide analytical methods, CGH and array-CGH, and discuss their potential in molecular characterization of cancers.     

A molecular cytogenetic approach to studying platinum resistance

The technique of comparative genomic hybridisation (CGH) has until recently been used to screen for common genomic abnormalities in fresh tumour material; it has identified previously unrecognised regions of amplification associated with poor prognosis subtypes of breast cancer and lymphoma. Our group has applied this technique to resistant cell lines and their sensitive counterparts in order to define chromosomal abnormalities associated with acquired drug resistance. We have demonstrated the applicability of this technique to the study of drug resistance using cell lines with known mechanisms of resistance. The ability to detect novel genomic alterations in cell lines with novel mechanisms of resistance was also demonstrated. We subsequently examined the CGH profiles of seven different cell lines made resistant to three platinum analogues and showed the most consistent abnormalities to involve over-representation of regions 4q and 6q. More recently, we have applied the CGH technique to a series of testicular germ cell tumours (TGCTs) collected as formalin-fixed paraffin-embedded biopsy specimens from patients, both pre- and post-therapy using a platinum-based regimen (POMB/ACE). Previous reports have shown over-representation of X, 7q, 8q and 12p and loss of 13q to occur in 25% of primary TGCTs. Over-representation of 12p was confirmed in the majority of these biopsy samples; deletion of 13q was noted in the initial biopsies of several patients. We also demonstrated alterations of 4p, 4q, 5q and 6q in this series of patients. Newly acquired deletions of 2q and 18q and amplifications of 8q were frequently observed in post-chemotherapy samples from resistant tumours. The CGH studies on these patients with TGCT will not only enable us to correlate our observations on clinical material with those from long-term cell lines, but should also identify sites of key genes involved in clinical platinum resistance.     

Copy number variants at Williams–Beuren syndrome 7q11.23 region

Copy number variants (CNVs) of the Williams–Beuren syndrome (WBS) 7q11.23 region are responsible for neurodevelopmental disorders with multi-system involvement and variable expressivity. Typical features of WBS microdeletion comprise a recognizable pattern of facial dysmorphisms, supravalvular aortic stenosis, connective tissue abnormalities, hypercalcemia, and a distinctive neurobehavioral phenotype. Conversely, the phenotype of patients carrying the 7q11.23 reciprocal duplications includes less distinctive facial dysmorphisms and prominent speech delay. The common deletion/duplication ranges in size from 1.5 to 1.8 Mb and encompasses approximately 28 genes. This region is flanked by low copy repeats (LCRs) with greater than ~97% identity, which can mediate non-allelic homologous recombination resulting from misalignment of LCRs during meiosis. A clear genotype–phenotype correlation has been established in WBS only for the elastin gene, which is responsible for the vascular and connective tissue abnormalities. The molecular substrates underlying the other clinical features of 7q11.23 CNVs, including the neurocognitive phenotypes, are still debated. Recent studies suggest that besides the role of the genes in the deleted/duplicated interval, multiple factors such as regulatory sequences, epigenetic mechanisms, parental origin of the CNV, and nucleotide variations in the non-deleted/duplicated allele may be important in determining the variable expressivity of 7q11.23 CNV phenotypes. Here, we review the clinical and molecular findings and the recent insights on genomic disorders associated with CNVs involving the 7q11.23 region.     

From amplification to gene in thyroid cancer: a high-resolution mapped bacterial-artificial-chromosome resource for cancer chromosome aberrations guides gene discovery after comparative genome hybridization.

Chromosome rearrangements associated with neoplasms provide a rich resource for definition of the pathways of tumorigenesis. The power of comparative genome hybridization (CGH) to identify novel genes depends on the existence of suitable markers, which are lacking throughout most of the genome. We now report a general approach that translates CGH data into higher-resolution genomic-clone data that are then used to define the genes located in aneuploid regions. We used CGH to study 33 thyroid-tumor DNAs and two tumor-cell-line DNAs. The results revealed amplifications of chromosome band 2p21, with less-intense amplification on 2p13, 19q13.1, and 1p36 and with least-intense amplification on 1p34, 1q42, 5q31, 5q33-34, 9q32-34, and 14q32. To define the 2p21 region amplified, a dense array of 373 FISH-mapped chromosome 2 bacterial artificial chromosomes (BACs) was constructed, and 87 of these were hybridized to a tumor-cell line. Four BACs carried genomic DNA that was amplified in these cells. The maximum amplified region was narrowed to 3-6 Mb by multicolor FISH with the flanking BACs, and the minimum amplicon size was defined by a contig of 420 kb. Sequence analysis of the amplified BAC 1D9 revealed a fragment of the gene, encoding protein kinase C epsilon (PKCepsilon), that was then shown to be amplified and rearranged in tumor cells. In summary, CGH combined with a dense mapped resource of BACs and large-scale sequencing has led directly to the definition of PKCepsilon as a previously unmapped candidate gene involved in thyroid tumorigenesis.     

Next Generation Sequencing of Acute Myeloid Leukemia: Influencing Prognosis

Background

Epilepsy is genetically complex neurological disorder affecting millions of people of different age groups varying in its type and severity. Copy number variants (CNVs) are key players in the genetic etiology of numerous neurodevelopmental disorders and prior findings also revealed that chromosomal aberrations are more susceptible against the pathogenesis of epilepsy. Novel technologies, such as array comparative genomic hybridization (array-CGH), may help to uncover the pathogenic CNVs in patients with epilepsy.

Results

This study was carried out by high density whole genome array-CGH analysis with blood DNA samples from a cohort of 22 epilepsy patients to search for CNVs associated with epilepsy. Pathogenic rearrangements which include 6p12.1 microduplications in 5 patients covering a total region of 99.9kb and 7q32.3 microdeletions in 3 patients covering a total region of 63.9kb were detected. Two genes BMP5 and PODXL were located in the predicted duplicated and deleted regions respectively. Furthermore, these CNV findings were confirmed by qPCR.

Conclusion

We have described, for the first time, several novel CNVs/genes implicated in epilepsy in the Saudi population. These findings enable us to better describe the genetic variations in epilepsy, and could provide a foundation for understanding the critical regions of the genome which might be involved in the development of epilepsy.

     

Gene amplification profiling of esophageal squamous cell carcinomas by DNA array CGH

Gene amplification is one of the basic mechanisms that lead to overexpression of oncogenes. DNA array comparative genomic hybridization (CGH) has great potential for comprehensive analysis of both a relative gene-copy number and altered chromosomal regions in cancers, which enables us to identify new amplified genes and unstable chromosomal loci. We examined the amplification status in 32 esophageal squamous cell carcinomas (ESCCs) and 13 ESCC cell lines on 51 frequently amplified loci in a variety of cancers by both DNA array CGH and Southern blot analyses. The 1p34 locus containing MYCL1, 2p24 (MYCN), 7p12 (EGFR), and 12q14 (MDM2) were amplified in one of the 32 cases (3%), and the 17q12 locus (ERBB2) and 8p11 (FGFR1) in two of the 32 cases (6%), while only the 11q13 locus (Cyclin D1, FGF4, and EMS1) was frequently amplified (28%, 9/32), demonstrating this locus to be a major target in ESCCs. One locus, 8q24 (c-MYC) was found to be amplified only in the cell lines. Eight out of 51 loci (15.7%) were found to be amplified in at least one of the 32 primary ESCCs or the 13 ESCC cell lines, suggesting that chromosomal loci frequently amplified in a type of human cancer may also be amplified in other types of cancers. This paper is the first report of an application of DNA array CGH to ESCCs.     

Constitutional promoter methylation and risk of familial melanoma

Constitutional epigenetic changes detected in blood or non-disease involving tissues have been associated with disease susceptibility. We measured promoter methylation of CDKN2A (p16 and p14ARF) and 13 melanoma-related genes using bisulfite pyrosequencing of blood DNA from 114 cases and 122 controls in 64 melanoma-prone families (26 segregating CDKN2A germline mutations). We also obtained gene expression data for these genes using microarrays from the same blood samples. We observed that CDKN2A epimutation is rare in melanoma families, and therefore is unlikely to cause major susceptibility in families without CDKN2A mutations. Although methylation levels for most gene promoters were very low (<5%), we observed a significantly reduced promoter methylation (odds ratio = 0.63, 95% confidence interval = 0.50, 0.80, P < 0.001) and increased expression (fold change = 1.27, P = 0.048) for TNFRSF10C in melanoma cases. Future research in large prospective studies using both normal and melanoma tissues is required to assess the significance of TNFRSF10C methylation and expression changes in melanoma susceptibility.     

BAGE genes generated by juxtacentromeric reshuffling in the Hominidae lineage are under selective pressure

In this paper, we show that the BAGE (B melanoma antigen) gene family was generated by chromosome rearrangements that occurred during the evolution of hominoids. An 84-kb DNA fragment derived from the phylogenetic 7q36 region was duplicated in the juxtacentromeric region of either chromosome 13 or chromosome 21. The duplicated region contained a fragment of the MLL3 gene, which, after juxtacentromeric reshuffling, generated the ancestral BAGE gene. Then, this ancestral gene gave rise to several independent genes through successive rounds of inter- and intrachromosome duplications. Comparison of synonymous and nonsynonymous mutations in putative coding regions shows that BAGE genes, but not the BAGE gene fragments, are under selective pressure. Our data strongly suggest that BAGE proteins have a function and that juxtacentromeric regions, whose plasticity is now largely proved, are not a simple junkyard of gene fragments, but may be the birth site of novel genes.     

NIPA2 located in 15q11.2 is mutated in patients with childhood absence epilepsy

While pathogenic copy number variations (CNVs) in 15q11.2 were recently identified in Caucasian patients with idiopathic generalized epilepsies (IGEs), the epilepsy-associated gene(s) in this region is/are still unknown. Our study investigated whether the CNVs in 15q11.2 are associated with childhood absence epilepsy (CAE) in Chinese patients and whether the selective magnesium transporter NIPA2 gene affected by 15q11.2 microdeletions is a susceptive gene for CAE. We assessed IGE-related CNVs by Affymetrix SNP 5.0 microarrays in 198 patients with CAE and 198 controls from northern China, and verified the identified CNVs by high-density oligonucleotide-based CGH microarrays. The coding region and exon-intron boundaries of NIPA2 were sequenced in all 380 patients with CAE and 400 controls. 15q11.2 microdeletions were detected in 3 of 198 (1.5%) patients and in no controls. Furthermore, we identified point mutations or indel in a heterozygous state of the NIPA2 gene in 3 out of 380 patients, whereas they were absent in 700 controls (P = 0.043). These mutations included two novel missense mutations (c.532A>T, p.I178F; c.731A>G, p.N244S) and one small novel insertion (c.1002_1003insGAT, p.N334_335EinsD). No NIPA2 mutation was found in 400 normal controls. We first identified that NIPA2, encoding a selective magnesium transporter, is a susceptible gene of CAE, and 15q11.2 microdeletions are important pathogenic CNVs for CAE with higher frequency in Chinese populations than that previously reported in Caucasians. The haploinsufficiency of NIPA2 may be a mechanism underlying the neurological phenotypes of 15q11.2 microdeletions.     

Constitutional promoter methylation and risk of familial melanoma

Constitutional epigenetic changes detected in blood or non-disease involving tissues have been associated with disease susceptibility. We measured promoter methylation of CDKN2A (p16 and p14ARF) and 13 melanoma-related genes using bisulfite pyrosequencing of blood DNA from 114 cases and 122 controls in 64 melanoma-prone families (26 segregating CDKN2A germline mutations). We also obtained gene expression data for these genes using microarrays from the same blood samples. We observed that CDKN2A epimutation is rare in melanoma families, and therefore is unlikely to cause major susceptibility in families without CDKN2A mutations. Although methylation levels for most gene promoters were very low (<5%), we observed a significantly reduced promoter methylation (odds ratio = 0.63, 95% confidence interval = 0.50, 0.80, P < 0.001) and increased expression (fold change = 1.27, P = 0.048) for TNFRSF10C in melanoma cases. Future research in large prospective studies using both normal and melanoma tissues is required to assess the significance of TNFRSF10C methylation and expression changes in melanoma susceptibility.     

Multiple recurrent de novo CNVs, including duplications of the 7q11.23 Williams syndrome region, are strongly associated with autism

We have undertaken a genome-wide analysis of rare copy-number variation (CNV) in 1124 autism spectrum disorder (ASD) families, each comprised of a single proband, unaffected parents, and, in most kindreds, an unaffected sibling. We find significant association of ASD with de novo duplications of 7q11.23, where the reciprocal deletion causes Williams-Beuren syndrome, characterized by a highly social personality. We identify rare recurrent de novo CNVs at five additional regions, including 16p13.2 (encompassing genes USP7 and C16orf72) and Cadherin 13, and implement a rigorous approach to evaluating the statistical significance of these observations. Overall, large de novo CNVs, particularly those encompassing multiple genes, confer substantial risks (OR = 5.6; CI = 2.6-12.0, p = 2.4 × 10(-7)). We estimate there are 130-234 ASD-related CNV regions in the human genome and present compelling evidence, based on cumulative data, for association of rare de novo events at 7q11.23, 15q11.2-13.1, 16p11.2, and Neurexin 1.     

Genome-wide mapping of copy number variation in humans: comparative analysis of high resolution array platforms

Accurate and efficient genome-wide detection of copy number variants (CNVs) is essential for understanding human genomic variation, genome-wide CNV association type studies, cytogenetics research and diagnostics, and independent validation of CNVs identified from sequencing based technologies. Numerous, array-based platforms for CNV detection exist utilizing array Comparative Genome Hybridization (aCGH), Single Nucleotide Polymorphism (SNP) genotyping or both. We have quantitatively assessed the abilities of twelve leading genome-wide CNV detection platforms to accurately detect Gold Standard sets of CNVs in the genome of HapMap CEU sample NA12878, and found significant differences in performance. The technologies analyzed were the NimbleGen 4.2 M, 2.1 M and 3×720 K Whole Genome and CNV focused arrays, the Agilent 1×1 M CGH and High Resolution and 2×400 K CNV and SNP+CGH arrays, the Illumina Human Omni1Quad array and the Affymetrix SNP 6.0 array. The Gold Standards used were a 1000 Genomes Project sequencing-based set of 3997 validated CNVs and an ultra high-resolution aCGH-based set of 756 validated CNVs. We found that sensitivity, total number, size range and breakpoint resolution of CNV calls were highest for CNV focused arrays. Our results are important for cost effective CNV detection and validation for both basic and clinical applications.     

Reference-unbiased copy number variant analysis using CGH microarrays

Comparative genomic hybridization (CGH) microarrays have been used to determine copy number variations (CNVs) and their effects on complex diseases. Detection of absolute CNVs independent of genomic variants of an arbitrary reference sample has been a critical issue in CGH array experiments. Whole genome analysis using massively parallel sequencing with multiple ultra-high resolution CGH arrays provides an opportunity to catalog highly accurate genomic variants of the reference DNA (NA10851). Using information on variants, we developed a new method, the CGH array reference-free algorithm (CARA), which can determine reference-unbiased absolute CNVs from any CGH array platform. The algorithm enables the removal and rescue of false positive and false negative CNVs, respectively, which appear due to the effects of genomic variants of the reference sample in raw CGH array experiments. We found that the CARA remarkably enhanced the accuracy of CGH array in determining absolute CNVs. Our method thus provides a new approach to interpret CGH array data for personalized medicine.     

Cutaneous melanoma: fishing with chips

DNA microarray technology is a versatile platform that allows rapid genetic analysis to take place on a genome-wide scale and has revolutionized the way cancers are studied. This platform has enabled researchers to characterize mechanisms central to tumorigenesis and understand important molecular events in the multi-step tumor progression model of cutaneous melanoma and other cancers. In melanoma, multiple global gene expression profiling studies using various DNA microarray platforms and various experimental designs have been performed. Each study has been able to capture and characterize either the involvement of a novel pathway or a novel cause-effect-relationship. The use of microarrays to define subclasses, to identify differentially regulated genes within a mutational context to analyze epigenetically regulated genes has resulted in an unprecedented understanding of the biology of cutaneous melanoma that may lead to more accurate diagnosis, more comprehensive prognosis, prediction and more effective therapeutic interventions. Related DNA microarray platforms like array-comparative genomic hybridization (CGH) have also been instrumental to identify many non-random chromosomal alterations; however, studies identifying validated targets as a result of CGH are limited. Thus, there exists significant opportunity to discover novel melanoma genes and translate such discoveries into meaningful clinical endpoints. In this review, we focus on various DNA microarray-based studies performed in cutaneous melanoma and summarize our current understanding of the genetics and biology of melanoma progression derived from accumulating genomic information.     

The functional myosin light chain kinase (MYLK) gene localizes with marker D3S3552 on human chromosome 3q21 in a >5-Mb yeast artificial chromosome region and is not linked to olfactory receptor genes

The myosin light chain kinase (MYLK) gene is duplicated on human chromosome 3 (3q13-->q21; 3p13), two sites known to contain olfactory receptor (OR) genes. The 3p13 site contains a MYLK pseudogene (MYLKP) associated with a cluster of OR pseudogenes and therefore could have arisen from the duplication of a large region in 3q13-->q21. Here, we present the localization of the MYLK gene in a >5-Mb region of the chromosome 3q21 integrated map. MYLK colocalizes with marker D3S3552. OR genes are absent from this region, suggesting that the 3p13 duplicated region incurred further rearrangements during evolution.