Copy number variations and cognitive phenotypes in unselected populations

IMPORTANCE: The association of copy number variations (CNVs), differing numbers of copies of genetic sequence at locations in the genome, with phenotypes such as intellectual disability has been almost exclusively evaluated using clinically ascertained cohorts. The contribution of these genetic variants to cognitive phenotypes in the general population remains unclear. OBJECTIVE: To investigate the clinical features conferred by CNVs associated with known syndromes in adult carriers without clinical preselection and to assess the genome-wide consequences of rare CNVs (frequency ≤0.05%; size ≥250 kilobase pairs [kb]) on carriers' educational attainment and intellectual disability prevalence in the general population. DESIGN, SETTING, AND PARTICIPANTS: The population biobank of Estonia contains 52,000 participants enrolled from 2002 through 2010. General practitioners examined participants and filled out a questionnaire of health- and lifestyle-related questions, as well as reported diagnoses. Copy number variant analysis was conducted on a random sample of 7877 individuals and genotype-phenotype associations with education and disease traits were evaluated. Our results were replicated on a high-functioning group of 993 Estonians and 3 geographically distinct populations in the United Kingdom, the United States, and Italy. MAIN OUTCOMES AND MEASURES: Phenotypes of genomic disorders in the general population, prevalence of autosomal CNVs, and association of these variants with educational attainment (from less than primary school through scientific degree) and prevalence of intellectual disability. RESULTS: Of the 7877 in the Estonian cohort, we identified 56 carriers of CNVs associated with known syndromes. Their phenotypes, including cognitive and psychiatric problems, epilepsy, neuropathies, obesity, and congenital malformations are similar to those described for carriers of identical rearrangements ascertained in clinical cohorts. A genome-wide evaluation of rare autosomal CNVs (frequency, ≤0.05%; ≥250 kb) identified 831 carriers (10.5%) of the screened general population. Eleven of 216 (5.1%) carriers of a deletion of at least 250 kb (odds ratio [OR], 3.16; 95% CI, 1.51-5.98; P = 1.5e-03) and 6 of 102 (5.9%) carriers of a duplication of at least 1 Mb (OR, 3.67; 95% CI, 1.29-8.54; P = .008) had an intellectual disability compared with 114 of 6819 (1.7%) in the Estonian cohort. The mean education attainment was 3.81 (P = 1.06e-04) among 248 (≥250 kb) deletion carriers and 3.69 (P = 5.024e-05) among 115 duplication carriers (≥1 Mb). Of the deletion carriers, 33.5% did not graduate from high school (OR, 1.48; 95% CI, 1.12-1.95; P = .005) and 39.1% of duplication carriers did not graduate high school (OR, 1.89; 95% CI, 1.27-2.8; P = 1.6e-03). Evidence for an association between rare CNVs and lower educational attainment was supported by analyses of cohorts of adults from Italy and the United States and adolescents from the United Kingdom. CONCLUSIONS AND RELEVANCE: Known pathogenic CNVs in unselected, but assumed to be healthy, adult populations may be associated with unrecognized clinical sequelae. Additionally, individually rare but collectively common intermediate-size CNVs may be negatively associated with educational attainment. Replication of these findings in additional population groups is warranted given the potential implications of this observation for genomics research, clinical care, and public health.

Genome-wide DNA profiling of marginal zone lymphomas identifies subtype-specific lesions with an impact on the clinical outcome

Marginal zone B-cell lymphomas (MZLs) have been divided into 3 distinct subtypes (extranodal MZLs of mucosa-associated lymphoid tissue [MALT] type, nodal MZLs, and splenic MZLs). Nevertheless, the relationship between the subtypes is still unclear. We performed a comprehensive analysis of genomic DNA copy number changes in a very large series of MZL cases with the aim of addressing this question. Samples from 218 MZL patients (25 nodal, 57 MALT, 134 splenic, and 2 not better specified MZLs) were analyzed with the Affymetrix Human Mapping 250K SNP arrays, and the data combined with matched gene expression in 33 of 218 cases. MALT lymphoma presented significantly more frequently gains at 3p, 6p, 18p, and del(6q23) (TNFAIP3/A20), whereas splenic MZLs was associated with del(7q31), del(8p). Nodal MZLs did not show statistically significant differences compared with MALT lymphoma while lacking the splenic MZLs-related 7q losses. Gains of 3q and 18q were common to all 3 subtypes. del(8p) was often present together with del(17p) (TP53). Although del(17p) did not determine a worse outcome and del(8p) was only of borderline significance, the presence of both deletions had a highly significant negative impact on the outcome of splenic MZLs.

USING THE R PACKAGE crlmm FOR GENOTYPING AND COPY NUMBER ESTIMATION

Genotyping platforms such as Affymetrix can be used to assess genotype-phenotype as well as copy number-phenotype associations at millions of markers. While genotyping algorithms are largely concordant when assessed on HapMap samples, tools to assess copy number changes are more variable and often discordant. One explanation for the discordance is that copy number estimates are susceptible to systematic differences between groups of samples that were processed at different times or by different labs. Analysis algorithms that do not adjust for batch effects are prone to spurious measures of association. The R package crlmm implements a multilevel model that adjusts for batch effects and provides allele-specific estimates of copy number. This paper illustrates a workflow for the estimation of allele-specific copy number, develops markerand study-level summaries of batch effects, and demonstrates how the marker-level estimates can be integrated with complimentary Bioconductor software for inferring regions of copy number gain or loss. All analyses are performed in the statistical environment R. A compendium for reproducing the analysis is available from the author’s website (http://www.biostat.jhsph.edu/~rscharpf/crlmmCompendium/index.html).

A compendium of myeloma-associated chromosomal copy number abnormalities and their prognostic value.

To obtain a comprehensive genomic profile of presenting multiple myeloma cases we performed high-resolution single nucleotide polymorphism mapping array analysis in 114 samples alongside 258 samples analyzed by U133 Plus 2.0 expression array (Affymetrix). We examined DNA copy number alterations and loss of heterozygosity (LOH) to define the spectrum of minimally deleted regions in which relevant genes of interest can be found. The most frequent deletions are located at 1p (30%), 6q (33%), 8p (25%), 12p (15%), 13q (59%), 14q (39%), 16q (35%), 17p (7%), 20 (12%), and 22 (18%). In addition, copy number-neutral LOH, or uniparental disomy, was also prevalent on 1q (8%), 16q (9%), and X (20%), and was associated with regions of gain and loss. Based on fluorescence in situ hybridization and expression quartile analysis, genes of prognostic importance were found to be located at 1p (FAF1, CDKN2C), 1q (ANP32E), and 17p (TP53). In addition, we identified common homozygously deleted genes that have functions relevant to myeloma biology. Taken together, these analyses indicate that the crucial pathways in myeloma pathogenesis include the nuclear factor-κB pathway, apoptosis, cell-cycle regulation, Wnt signaling, and histone modifications. This study was registered at http://isrctn.org as ISRCTN68454111.

Genomic imbalances associated with mullerian aplasia

Background: Aplasia of the mullerian ducts leads to absence of the uterine corpus, uterine cervix, and upper (superior) vagina. Patients with mullerian aplasia (MA) often exhibit additional clinical features such as renal, vertebral and cardiac defects. A number of different syndromes have been associated with MA, and in most cases its aetiology remains poorly understood. Objective and methods: 14 syndromic patients with MA and 46, XX G-banded karyotype were screened for DNA copy number changes by similar to 1 Mb whole genome bacterial artificial chromosome (BAC) array based comparative genomic hybridisation (CGH). The detected alterations were validated by an independent method and further mapped by high resolution oligo-arrays. Results: Submicroscopic genomic imbalances affecting the 1q21.1, 17q12, 22q11.21, and Xq21.31 chromosome regions were detected in four probands. Presence of the alterations in the normal mother of one patient suggests incomplete penetrance and/or variable expressivity. Conclusion: 4 of the 14 patients (29%) were found to have cryptic genomic alterations. The imbalances on 22q11.21 support recent findings by us and others that alterations in this chromosome region may result in impairment of mullerian duct development. The remaining imbalances indicate involvement of previously unknown chromosome regions in MA, and point specifically to LHX1 and KLHL4 as candidate genes.

Li-Fraumeni-like syndrome associated with a large BRCA1 intragenic deletion

Background: Li-Fraumeni (LFS) and Li-Fraumeni-like (LFL) syndromes are associated to germline TP53 mutations, and are characterized by the development of central nervous system tumors, sarcomas, adrenocortical carcinomas, and other early-onset tumors. Due to the high frequency of breast cancer in LFS/LFL families, these syndromes clinically overlap with hereditary breast cancer (HBC). Germline point mutations in BRCA1, BRCA2, and TP53 genes are associated with high risk of breast cancer. Large rearrangements involving these genes are also implicated in the HBC phenotype. Methods: We have screened DNA copy number changes by MLPA on BRCA1, BRCA2, and TP53 genes in 23 breast cancer patients with a clinical diagnosis consistent with LFS/LFL; most of these families also met the clinical criteria for other HBC syndromes. Results: We found no DNA copy number alterations in the BRCA2 and TP53 genes, but we detected in one patient a 36.4 Kb BRCA1 microdeletion, confirmed and further mapped by array-CGH, encompassing exons 9-19. Breakpoints sequencing analysis suggests that this rearrangement was mediated by flanking Alu sequences. Conclusion: This is the first description of a germline intragenic BRCA1 deletion in a breast cancer patient with a family history consistent with both LFL and HBC syndromes. Our results show that large rearrangements in these known cancer predisposition genes occur, but are not a frequent cause of cancer susceptibility.

The molecular genetics of breast cancer: The contribution of comparative genomic hybridization

Comparative genomic hybridization (CGH) has been the technique of choice over the last 10 years for mapping DNA copy number changes in human tumors. Here we review the literature to demonstrate how CGH has contributed to the comprehension of molecular aspects of breast tumorigenesis. At least two distinct molecular pathways of breast cancer have been characterized that show a strong correlation with histological grade. It seems that grade I invasive ductal carcinomas (IDCs) arise from well-differentiated ductal carcinoma in situ (DCIS), whereas grade III IDCs come from poorly differentiated DCIS. In addition, dedifferentiation from a low- to a high-grade breast cancer has proven an unlikely phenomenon. CGH has been instrumental in dissecting distinct molecular pathways toward breast malignancy and in establishing a direct relationship between genotype and clinical pathological features. (C) 2005 Elsevier GrnbH. All rights reserved.

Microarray-based comparative genomic hybridisation of breast cancer patients receiving neoadjuvant chemotherapy

We analysed the molecular genetic profiles of breast cancer samples before and after neoadjuvant chemotherapy with combination doxorubicin and cyclophosphamide (AC). DNA was obtained from microdissected frozen breast core biopsies from 44 patients before chemotherapy. Additional samples were obtained before the second course of chemotherapy (D21) and after the completion of the treatment (surgical specimens) in 17 and 21 patients, respectively. Microarray-based comparative genome hybridisation was performed using a platform containing approx5800 bacterial artificial chromosome clones (genome-wide resolution: 0.9 Mb). Analysis of the 44 pretreatment biopsies revealed that losses of 4p, 4q, 5q, 12q13.11–12q13.12, 17p11.2 and 17q11.2; and gains of 1p, 2p, 7q, 9p, 11q, 19p and 19q were significantly associated with oestrogen receptor negativity. 16q21–q22.1 losses were associated with lobular and 8q24 gains with ductal types. Losses of 5q33.3–q4 and 18p11.31 and gains of 6p25.1–p25.2 and Xp11.4 were associated with HER2 amplification. No correlations between DNA copy number changes and clinical response to AC were found. Microarray-based comparative genome hybridisation analysis of matched pretreatment and D21 biopsies failed to identify statistically significant differences, whereas a comparison between matched pretreatment and surgical samples revealed a statistically significant acquired copy number gain on 11p15.2–11p15.5. The modest chemotherapy-driven genomic changes, despite profound loss of cell numbers, suggest that there is little therapeutic selection of resistant non-modal cell lineages.

Characterising the TP53-deleted subgroup of chronic lymphocytic leukemia: an analysis of additional cytogenetic abnormalities detected by interphase fluorescence in situ hybridisation and array-based comparative genomic hybridisation.

Deletion of the TP53 gene on chromosome 17p13.1 is the prognostic factor associated with the shortest survival in CLL. We used array-based comparative genomic hybridisation (arrayCGH) to identify additional DNA copy number changes in peripheral blood samples from 74 LRF CLL4 trial patients, 37 with >or=5% and 37 without TP53-deleted cells. ArrayCGH reliably detected deletions on 17p, including the TP53 locus, in cases with >or=50%TP53-deleted cells detected by fluorescence in situ hybridisation, plus seven additional cases with deleted regions on 17p excluding TP53. Losses on chromosomal regions 18p and/or 20p were found exclusively in cases with >or=5%TP53-deleted cells (por=5%TP53-deleted cases (p=0.02). In particular, amplification of 2p and deletion of 6q were both more frequent. Cases with >20%TP53-deleted cells had the worst prognosis in the LRF CLL4 trial.

Investigation of Copy Number Variation in Children with Conotruncal Heart Defects

Background:Congenital heart defects (CHD) are the most prevalent group of structural abnormalities at birth and one of the main causes of infant morbidity and mortality. Studies have shown a contribution of the copy number variation in the genesis of cardiac malformations.Objectives:Investigate gene copy number variation (CNV) in children with conotruncal heart defect.Methods:Multiplex ligation-dependent probe amplification (MLPA) was performed in 39 patients with conotruncal heart defect. Clinical and laboratory assessments were conducted in all patients. The parents of the probands who presented abnormal findings were also investigated.Results:Gene copy number variation was detected in 7/39 patients: 22q11.2 deletion, 22q11.2 duplication, 15q11.2 duplication, 20p12.2 duplication, 19p deletion, 15q and 8p23.2 duplication with 10p12.31 duplication. The clinical characteristics were consistent with those reported in the literature associated with the encountered microdeletion/microduplication. None of these changes was inherited from the parents.Conclusions:Our results demonstrate that the technique of MLPA is useful in the investigation of microdeletions and microduplications in conotruncal congenital heart defects. Early diagnosis of the copy number variation in patients with congenital heart defect assists in the prevention of morbidity and decreased mortality in these patients.

Phenotypic consequences of copy number variation: insights from Smith-Magenis and Potocki-Lupski syndrome mouse models.

A large fraction of genome variation between individuals is comprised of submicroscopic copy number variation of genomic DNA segments. We assessed the relative contribution of structural changes and gene dosage alterations on phenotypic outcomes with mouse models of Smith-Magenis and Potocki-Lupski syndromes. We phenotyped mice with 1n (Deletion/+), 2n (+/+), 3n (Duplication/+), and balanced 2n compound heterozygous (Deletion/Duplication) copies of the same region. Parallel to the observations made in humans, such variation in gene copy number was sufficient to generate phenotypic consequences: in a number of cases diametrically opposing phenotypes were associated with gain versus loss of gene content. Surprisingly, some neurobehavioral traits were not rescued by restoration of the normal gene copy number. Transcriptome profiling showed that a highly significant propensity of transcriptional changes map to the engineered interval in the five assessed tissues. A statistically significant overrepresentation of the genes mapping to the entire length of the engineered chromosome was also found in the top-ranked differentially expressed genes in the mice containing rearranged chromosomes, regardless of the nature of the rearrangement, an observation robust across different cell lineages of the central nervous system. Our data indicate that a structural change at a given position of the human genome may affect not only locus and adjacent gene expression but also "genome regulation." Furthermore, structural change can cause the same perturbation in particular pathways regardless of gene dosage. Thus, the presence of a genomic structural change, as well as gene dosage imbalance, contributes to the ultimate phenotype.

DNA Topoisomerase I Gene Copy Number and mRNA Expression Assessed as Predictive Biomarkers for Adjuvant Irinotecan in Stage II/III Colon Cancer.

PURPOSE: Prospective-retrospective assessment of theTOP1gene copy number andTOP1mRNA expression as predictive biomarkers for adjuvant irinotecan in stage II/III colon cancer. EXPERIMENTAL DESIGN: Formalin-fixed, paraffin-embedded tissue microarrays were obtained from an adjuvant colon cancer trial (PETACC3) where patients were randomized to 5-fluorouracil/folinic acid with or without additional irinotecan.TOP1copy number status was analyzed by fluorescencein situhybridization (FISH) using aTOP1/CEN20 dual-probe combination.TOP1mRNA data were available from previous analyses. RESULTS: TOP1FISH and follow-up data were obtained from 534 patients.TOP1gain was identified in 27% using a single-probe enumeration strategy (≥4TOP1signals per cell) and in 31% when defined by aTOP1/CEN20 ratio ≥ 1.5. The effect of additional irinotecan was not dependent onTOP1FISH status.TOP1mRNA data were available from 580 patients with stage III disease. Benefit of irinotecan was restricted to patients characterized byTOP1mRNA expression ≥ third quartile (RFS: HRadjusted, 0.59;P= 0.09; OS: HRadjusted, 0.44;P= 0.03). The treatment byTOP1mRNA interaction was not statistically significant, but in exploratory multivariable fractional polynomial interaction analysis, increasingTOP1mRNA values appeared to be associated with increasing benefit of irinotecan. CONCLUSIONS: In contrast to theTOP1copy number, a trend was demonstrated for a predictive property ofTOP1mRNA expression. On the basis ofTOP1mRNA, it might be possible to identify a subgroup of patients where an irinotecan doublet is a clinically relevant option in the adjuvant setting of colon cancer.Clin Cancer Res; 22(7); 1621-31. ©2015 AACR.

Sequence and Copy Number Analyses of HEXB Gene in Patients Affected by Sandhoff Disease: Functional Characterization of 9 Novel Sequence Variants

Sandhoff disease (SD) is a lysosomal disorder caused by mutations in the HEXB gene. To date, 43 mutations of HEXB have been described, including 3 large deletions. Here, we have characterized 14 unrelated SD patients and developed a Multiplex Ligation-dependent Probe Amplification (MLPA) assay to investigate the presence of large HEXB deletions. Overall, we identified 16 alleles, 9 of which were novel, including 4 sequence variation leading to aminoacid changes [c.626C>T (p.T209I), c.634C>A (p.H212N), c.926G>T (p.C309F), c.1451G>A (p.G484E)] 3 intronic mutations (c.1082+5G>A, c.1242+1G>A, c.1169+5G>A), 1 nonsense mutation c.146C>A (p.S49X) and 1 small in-frame deletion c.1260_1265delAGTTGA (p.V421_E422del). Using the new MLPA assay, 2 previously described deletions were identified. In vitro expression studies showed that proteins bearing aminoacid changes p.T209I and p.G484E presented a very low or absent activity, while proteins bearing the p.H212N and p.C309F changes retained a significant residual activity. The detrimental effect of the 3 novel intronic mutations on the HEXB mRNA processing was demonstrated using a minigene assay. Unprecedentedly, minigene studies revealed the presence of a novel alternative spliced HEXB mRNA variant also present in normal cells. In conclusion, we provided new insights into the molecular basis of SD and validated an MLPA assay for detecting large HEXB deletions.

Quantifying single gene copy number by measuring fluorescent probe lengths on combed genomic DNA

An approach was developed for the quantification of subtle gains and losses of genomic DNA. The approach relies on a process called molecular combing. Molecular combing consists of the extension and alignment of purified molecules of genomic DNA on a glass coverslip. It has the advantage that a large number of genomes can be combed per coverslip, which allows for a statistically adequate number of measurements to be made on the combed DNA. Consequently, a high-resolution approach to mapping and quantifying genomic alterations is possible. The approach consists of applying fluorescence hybridization to the combed DNA by using probes to identify the amplified region. Measurements then are made on the linear hybridization signals to ascertain the region's exact size. The reliability of the approach first was tested for low copy number amplifications by determining the copy number of chromosome 21 in a normal and trisomy 21 cell line. It then was tested for high copy number amplifications by quantifying the copy number of an oncogene amplified in the tumor cell line GTL-16. These results demonstrate that a wide range of amplifications can be accurately and reliably quantified. The sensitivity and resolution of the approach likewise was assessed by determining the copy number of a single allele (160 kb) alteration.