Candidate Luminal B Breast Cancer Genes Identified by Genome,Gene Expression and DNA Methylation Profiling

Breast cancers (BCs) of the luminal B subtype are estrogen receptor-positive (ER+), highly proliferative, resistant to standard therapies and have a poor prognosis. To better understand this subtype we compared DNA copy number aberrations (CNAs), DNA promoter methylation, gene expression profiles, and somatic mutations in nine selected genes, in 32 luminal B tumors with those observed in 156 BCs of the other molecular subtypes. Frequent CNAs included 8p11-p12 and 11q13.1-q13.2 amplifications, 7q11.22-q34, 8q21.12-q24.23, 12p12.3-p13.1, 12q13.11-q24.11, 14q21.1-q23.1, 17q11.1-q25.1, 20q11.23-q13.33 gains and 6q14.1-q24.2, 9p21.3-p24,3, 9q21.2, 18p11.31-p11.32 losses. A total of 237 and 101 luminal B-specific candidate oncogenes and tumor suppressor genes (TSGs) presented a deregulated expression in relation with their CNAs, including 11 genes previously reported associated with endocrine resistance. Interestingly, 88% of the potential TSGs are located within chromosome arm 6q, and seven candidate oncogenes are potential therapeutic targets. A total of 100 candidate oncogenes were validated in a public series of 5,765 BCs and the overexpression of 67 of these was associated with poor survival in luminal tumors. Twenty-four genes presented a deregulated expression in relation with a high DNA methylation level. FOXO3, PIK3CA and TP53 were the most frequent mutated genes among the nine tested. In a meta-analysis of next-generation sequencing data in 875 BCs, KCNB2 mutations were associated with luminal B cases while candidate TSGs MDN1 (6q15) and UTRN (6q24), were mutated in this subtype. In conclusion, we have reported luminal B candidate genes that may play a role in the development and/or hormone resistance of this aggressive subtype.     

Copy Number Gains at 8q24 and 20q11-q13 in Gastric Cancer Are More Common in Intestinal-Type than Diffuse-Type

The present study was aimed at discovering DNA copy number alterations (CNAs) involved in the carcinogenesis of stomach and at understanding their clinicopathological significances in the Korean population. DNA copy numbers were analyzed using Agilent 244K or 400K array comparative genomic hybridization (aCGH) in fresh-frozen tumor and matched normal tissues from 40 gastric cancer patients. Some of the detected CNA regions were validated using multiplex ligation-dependent probe amplification (MLPA) in six of the 40 patients and customized Agilent 60K aCGH in an independent set of 48 gastric cancers. The mRNA levels of genes at common CNA regions were analyzed using quantitative real-time PCR. Copy number gains were more common than losses across the entire genome in tumor tissues compared to matched normal tissues. The mean number of alterations per case was 64 for gains and 40 for losses, and the median aberration length was 44016 bp for gains and 4732 bp for losses. Copy number gains were frequently detected at 7p22.1 (20%), 8q24.21 (27%–30%), 8q24.3 (22%–48%), 13q34 (20%–31%), and 20q11-q13 (25%–30%), and losses at 3p14.2 (43%), 4q35.2 (27%), 6q26 (23%), and 17p13.3 (20%–23%). CNAs at 7p22.1, 13q34, and 17p13.3 have not been reported in other populations. Most of the copy number losses were associated with down-regulation of mRNA levels, but the correlation between copy number gains and mRNA expression levels varied in a gene-dependent manner. In addition, copy number gains tended to occur more commonly in intestinal-type cancers than in diffuse-type cancers. In conclusion, the present study suggests that copy number gains at 8q24 and 20q11-q13 and losses at 3p14.2 may be common events in gastric cancer but CNAs at 7p22.1, 13q34, and 17p13.3 may be Korean-specific.     

A comprehensive characterization of genome-wide copy number aberrations in colorectal cancer reveals novel oncogenes and patterns of alterations

To develop a comprehensive overview of copy number aberrations (CNAs) in stage-II/III colorectal cancer (CRC), we characterized 302 tumors from the PETACC-3 clinical trial. Microsatellite-stable (MSS) samples (n = 269) had 66 minimal common CNA regions, with frequent gains on 20 q (72.5%), 7 (41.8%), 8 q (33.1%) and 13 q (51.0%) and losses on 18 (58.6%), 4 q (26%) and 21 q (21.6%). MSS tumors have significantly more CNAs than microsatellite-instable (MSI) tumors: within the MSI tumors a novel deletion of the tumor suppressor WWOX at 16 q23.1 was identified (p<0.01). Focal aberrations identified by the GISTIC method confirmed amplifications of oncogenes including EGFR, ERBB2, CCND1, MET, and MYC, and deletions of tumor suppressors including TP53, APC, and SMAD4, and gene expression was highly concordant with copy number aberration for these genes. Novel amplicons included putative oncogenes such as WNK1 and HNF4A, which also showed high concordance between copy number and expression. Survival analysis associated a specific patient segment featured by chromosome 20 q gains to an improved overall survival, which might be due to higher expression of genes such as EEF1B2 and PTK6. The CNA clustering also grouped tumors characterized by a poor prognosis BRAF-mutant-like signature derived from mRNA data from this cohort. We further revealed non-random correlation between CNAs among unlinked loci, including positive correlation between 20 q gain and 8 q gain, and 20 q gain and chromosome 18 loss, consistent with co-selection of these CNAs. These results reinforce the non-random nature of somatic CNAs in stage-II/III CRC and highlight loci and genes that may play an important role in driving the development and outcome of this disease.     

Detection of DNA copy number abnormality by microarray expression analysis

Gene copy-number abnormalities (CNAs) are characteristic of solid tumors and are found in association with developmental abnormalities and/or mental retardation. The ultimate impact of CNAs is exerted by the altered expression of encoded genes. We have utilized high-density oligonucleotide arrays from Affymetrix to identify DNA CNAs via their impact on mRNA expression levels. In these studies, we have used three different trisomic cell lines (trisomy 9, trisomy 18, trisomy 21) as models of CNAs and have compared mRNA expression in those trisomic cells with that observed in diploid cell lines of matched tissue origin. Our data clearly show that genes from CNA chromosome regions are substantially over-represented (P<0.000001 by chi-square analysis) in the differentially expressed subset from comparisons of all three trisomic cell lines with normal matching cells. In addition, we have been able to detect the origin of the duplication by a statistical scan for over-expressed genes. These data show that microarray detection of differential mRNA expression can be used to identify significant DNA CNAs.     

Impact of Gene Dosage on Gene Expression,Biological Processes and Survival in Cervical Cancer: A Genome-Wide Follow-Up Study

We investigated the role of tumor copy number (CN)–altered genome (CN-AG) in the carcinogenesis of cervical cancer (CC), especially its effect on gene expression, biological processes, and patient survival. Fifty-nine human papillomavirus 16 (HPV16)-positive CCs were investigated with microarrays–31 for mapping CN-AG and 55 for global gene expression, with 27 CCs in common. Five-year survival was investigated in 55 patients. Deletions and amplifications >2.5 Mb were defined as CN alterations. The %CN-AG varied from 0 to 32.2% (mean = 8.1±8.9). Tumors were classified as low (mean = 0.5±0.6, n = 11), medium (mean = 5.4±2.4, n = 10), or high (mean = 19.2±6.6, n = 10) CN. The highest %CN-AG was found in 3q, which contributed an average of 55% of all CN alterations. Genome-wide, only 5.3% of CN-altered genes were deregulated directly by gene dosage. In contrast, the rate in fully duplicated 3q was twice as high. Amplification of 3q explained 23.2% of deregulated genes in whole tumors (r² = 0.232, p = 0.006; analysis of variance), including genes located in 3q and other chromosomes. A total of 862 genes were deregulated exclusively in high-CN tumors, but only 22.9% were CN altered. This suggests that the remaining genes are not deregulated directly by gene dosage, but by mechanisms induced in trans by CN-altered genes. Anaphase-promoting complex/cyclosome (APC/C)-dependent proteasome proteolysis, glycolysis, and apoptosis were upregulated, whereas cell adhesion and angiogenesis were downregulated exclusively in high-CN tumors. The high %CN-AG and upregulated gene expression profile of APC/C-dependent proteasome proteolysis were associated with poor patient survival (p<0.05, log-rank test). Along with glycolysis, they were linearly associated with FIGO stage (r>0.38, p<0.01, Spearman test). Therefore, inhibition of APC/C-dependent proteasome proteolysis and glycolysis could be useful for CC treatment. However, whether they are indispensable for tumor growth remains to be demonstrated.     

Flexible and accurate detection of genomic copy-number changes from aCGH

Genomic DNA copy-number alterations (CNAs) are associated with complex diseases, including cancer: CNAs are indeed related to tumoral grade, metastasis, and patient survival. CNAs discovered from array-based comparative genomic hybridization (aCGH) data have been instrumental in identifying disease-related genes and potential therapeutic targets. To be immediately useful in both clinical and basic research scenarios, aCGH data analysis requires accurate methods that do not impose unrealistic biological assumptions and that provide direct answers to the key question, "What is the probability that this gene/region has CNAs?" Current approaches fail, however, to meet these requirements. Here, we introduce reversible jump aCGH (RJaCGH), a new method for identifying CNAs from aCGH; we use a nonhomogeneous hidden Markov model fitted via reversible jump Markov chain Monte Carlo; and we incorporate model uncertainty through Bayesian model averaging. RJaCGH provides an estimate of the probability that a gene/region has CNAs while incorporating interprobe distance and the capability to analyze data on a chromosome or genome-wide basis. RJaCGH outperforms alternative methods, and the performance difference is even larger with noisy data and highly variable interprobe distance, both commonly found features in aCGH data. Furthermore, our probabilistic method allows us to identify minimal common regions of CNAs among samples and can be extended to incorporate expression data. In summary, we provide a rigorous statistical framework for locating genes and chromosomal regions with CNAs with potential applications to cancer and other complex human diseases.     

Methylation and expression analyses of Pallister-Killian syndrome reveal partial dosage compensation of tetrasomy 12p and hypomethylation of gene-poor regions on 12p

To ascertain the epigenomic features, i.e., the methylation, non-coding RNA, and gene expression patterns, associated with gain of i(12p) in Pallister-Killian syndrome (PKS), we investigated single cell clones, harboring either disomy 12 or tetrasomy 12p, from a patient with PKS. The i(12p)-positive cells displayed a characteristic expression and methylation signature. Of all the genes on 12p, 13% were overexpressed, including the ATN1, COPS7A, and NECAP1 genes in 12p13.31, a region previously implicated in PKS. However, the median expression fold change (1.3) on 12p was lower than expected by tetrasomy 12p. Thus, partial dosage compensation occurs in cells with i(12p). The majority (89%) of the significantly deregulated genes were not situated on 12p, indicating that global perturbation of gene expression is a key pathogenetic event in PKS. Three genes—ATP6V1G1 in 9q32, GMPS in 3q25.31, and TBX5 in 12q24.21—exhibited concomitant hypermethylation and decreased expression. The i(12p)-positive cells displayed global hypomethylation of gene-poor regions on 12p, a footprint previously associated with constitutional and acquired gains of whole chromosomes as well as with X-chromosome inactivation in females. We hypothesize that this non-genic hypomethylation is associated with chromatin processing that facilitates cellular adaptation to excess genetic material.     

利用cDNA微阵列-CGH筛选支气管上皮细胞恶性转化中的扩增基因

目的:基因组不稳定是导致肺癌发生与发展的重要分子机理之一。本研究旨在筛选支气管上皮细胞恶性转化过程中拷贝数扩增的基因。方法:利用业已建立的支气管上皮细胞体外恶性转化模型,通过cDNA微阵列-CGH技术对支气管上皮来源的永生化细胞和癌变细胞的基因拷贝数进行了检测,并对部分结果进行了实时PCR验证。结果:永生化BEP2D细胞染色体中的某些区域存在不同程度的扩增,包括5q31.3、9q32-33.1、14q22.2-23.1、19p13.12-13.13、20q13.12-13.31;恶性转化BERP35T2细胞染色体中的扩增区域集中在1p12-13.1、5q33.1、5q31.3、9q32、19p13.12-13.13;5q31.3、9q32、19q13.12-13.13是以上2种细胞系中的共同扩增区域。共检测到201个基因的拷贝数发生扩增,其中PCNA、TP53及GADD45A基因的异常扩增已经实时PCR进一步验证。结论:在支气管上皮细胞恶性转化过程中,病毒与低剂量辐射的双重作用使得某些重要基因的拷贝数发生扩增,因基因剂量增加而导致某些癌基因高表达可能是细胞恶性转化的重要机制之一。     

Identification of Novel Genomic Aberrations in AML-M5 in a Level of Array CGH

To assess the possible existence of unbalanced chromosomal abnormalities and delineate the characterization of copy number alterations (CNAs) of acute myeloid leukemia-M5 (AML-M5), R-banding karyotype, oligonucelotide array CGH and FISH were performed in 24 patients with AML-M5. A total of 117 CNAs with size ranging from 0.004 to 146.263 Mb was recognized in 12 of 24 cases, involving all chromosomes other than chromosome 1, 4, X and Y. Cryptic CNAs with size less than 5 Mb accounted for 59.8% of all the CNAs. 12 recurrent chromosomal alterations were mapped. Seven out of them were described in the previous AML studies and five were new candidate AML-M5 associated CNAs, including gains of 3q26.2-qter and 13q31.3 as well as losses of 2q24.2, 8p12 and 14q32. Amplication of 3q26.2-qter was the sole large recurrent chromosomal anomaly and the pathogenic mechanism in AML-M5 was possibly different from the classical recurrent 3q21q26 abnormality in AML. As a tumor suppressor gene, FOXN3, was singled out from the small recurrent CNA of 14q32, however, it is proved that deletion of FOXN3 is a common marker of myeloid leukemia rather than a specific marker for AML-M5 subtype. Moreover, the concurrent amplication of MLL and deletion of CDKN2A were noted and it might be associated with AML-M5. The number of CNA did not show a significant association with clinico-biological parameters and CR number of the 22 patients received chemotherapy. This study provided the evidence that array CGH served as a complementary platform for routine cytogenetic analysis to identify those cryptic alterations in the patients with AML-M5. As a subtype of AML, AML-M5 carries both common recurrent CNAs and unique CNAs, which may harbor novel oncogenes or tumor suppressor genes. Clarifying the role of these genes will contribute to the understanding of leukemogenic network of AML-M5.     

Integrated analysis of gene copy number, copy neutral LOH, and microRNA profiles in adult acute lymphoblastic leukemia

We adopted an integrated analysis of gene copy number alterations (CNAs), copy number neutral loss of heterozygosity (CNN LOH), and microRNA (miRNA) profiling in 21 adult acute lymphoblastic leukemia (ALL) patients. This study revealed the most frequent CNAs to be at chromosomes 9p, 7, and 17 and recurrent CNN LOH at 5p, 9p, and Xq. As for the most differentially expressed miRNAs, they included 8 upregulated and 14 downregulated miRNAs, of which miR-148a at 7p15.2, miR-22 at 17p13.3, miR-223 at Xq12, as well as miR-101-2 at 9p24.1 exhibited recurrent CNAs or CNN LOH. miR-101-2 was recurrently downregulated, and although the related CNN LOH was detected only in BCR-ABL1 negative cases (2/14), deletions of miR-101-2 were observed solely in BCR-ABL1 positive cases (4/7). Finally, BCR-ABL1 positive cases, in contrast to negative ones, were characterized by slightly, but still significantly, higher expression levels of miR-29b.     

Deletions in chromosome 4 differentially associated with the development of cervical cancer: evidence of slit2 as a candidate tumor suppressor gene

The aim of this study was to locate the candidate tumor suppressor genes (TSGs) loci in the chromosomal 4p15-16, 4q22-23 and 4q34-35 regions associated with the development of uterine cervical carcinoma (CA-CX). Deletion mapping of the regions by microsatellite markers identified six discrete areas with high frequency of deletions, viz. 4p16.2 (D1: 40%), 4p15.31 (D2: 35–38%), 4p15.2 (D3: 37–40%), 4q22.2 (D4: 34%), 4q34.2-34.3 (D5: 37–59%) and 4q35.1 (D6: 40–50%). Significant correlation was noted among the deleted regions D1, D2 and D3. The deletions in D1, D2, D5 and D6 regions are suggested to be associated with the cervical intraepithelial neoplasia (CIN), and deletions in the D2, D3, D5 and D6 regions seems to be associated with progression of CA-CX. The deletions in the D2 and D6 regions showed significant prognostic implications (P = 0.001; 0.02). The expression of the candidate TSG SLIT2 mapped to D2 region gradually reduced from normal cervix uteri →CIN → CA-CX. SLIT2 promoter hypermethylation was seen in 28% CIN samples and significantly increased with tumor progression (P = 0.04). Significant correlation was seen between SLIT2 deletion and its promoter methylation (P = 0.001), indicating that both these phenomena could occur simultaneously to inactivate this gene. Immunohistochemical analysis showed reduced expression of SLIT2 in cervical lesions and CA-CX cell lines. Although no mutation was detected in the SLIT2 promoter region (−432 to + 55 bp), CC and AA haplotypes were seen in −227 and −195 positions, respectively. Thus, it indicates that inactivation of SLIT2-ROBO1 signaling pathway may have an important role in CA-CX development.     

Integration of DNA copy number alterations and transcriptional expression analysis in human gastric cancer

Background

Genomic instability with frequent DNA copy number alterations is one of the key hallmarks of carcinogenesis. The chromosomal regions with frequent DNA copy number gain and loss in human gastric cancer are still poorly defined. It remains unknown how the DNA copy number variations contributes to the changes of gene expression profiles, especially on the global level.

Principal Findings

We analyzed DNA copy number alterations in 64 human gastric cancer samples and 8 gastric cancer cell lines using bacterial artificial chromosome (BAC) arrays based comparative genomic hybridization (aCGH). Statistical analysis was applied to correlate previously published gene expression data obtained from cDNA microarrays with corresponding DNA copy number variation data to identify candidate oncogenes and tumor suppressor genes. We found that gastric cancer samples showed recurrent DNA copy number variations, including gains at 5p, 8q, 20p, 20q, and losses at 4q, 9p, 18q, 21q. The most frequent regions of amplification were 20q12 (7/72), 20q12–20q13.1 (12/72), 20q13.1–20q13.2 (11/72) and 20q13.2–20q13.3 (6/72). The most frequent deleted region was 9p21 (8/72). Correlating gene expression array data with aCGH identified 321 candidate oncogenes, which were overexpressed and showed frequent DNA copy number gains; and 12 candidate tumor suppressor genes which were down-regulated and showed frequent DNA copy number losses in human gastric cancers. Three networks of significantly expressed genes in gastric cancer samples were identified by ingenuity pathway analysis.

Conclusions

This study provides insight into DNA copy number variations and their contribution to altered gene expression profiles during human gastric cancer development. It provides novel candidate driver oncogenes or tumor suppressor genes for human gastric cancer, useful pathway maps for the future understanding of the molecular pathogenesis of this malignancy, and the construction of new therapeutic targets.     

A transcriptomic computational analysis of mastic oil-treated Lewis lung carcinomas reveals molecular mechanisms targeting tumor cell growth and survival

Background

DNA copy number alterations are frequently observed in ovarian cancer, but it remains a challenge to identify the most relevant alterations and the specific causal genes in those regions.

Methods

We obtained high-resolution 500K SNP array data for 52 ovarian tumors and identified the most statistically significant minimal genomic regions with the most prevalent and highest-level copy number alterations (recurrent CNAs). Within a region of recurrent CNA, comparison of expression levels in tumors with a given CNA to tumors lacking that CNA and to whole normal ovary samples was used to select genes with CNA-specific expression patterns. A public expression array data set of laser capture micro-dissected (LCM) non-malignant fallopian tube epithelia and LCM ovarian serous adenocarcinoma was used to evaluate the effect of cell-type mixture biases.

Results

Fourteen recurrent deletions were detected on chromosomes 4, 6, 9, 12, 13, 15, 16, 17, 18, 22 and most prevalently on X and 8. Copy number and expression data suggest several apoptosis mediators as candidate drivers of the 8p deletions. Sixteen recurrent gains were identified on chromosomes 1, 2, 3, 5, 8, 10, 12, 15, 17, 19, and 20, with the most prevalent gains localized to 8q and 3q. Within the 8q amplicon, PVT1, but not MYC, was strongly over-expressed relative to tumors lacking this CNA and showed over-expression relative to normal ovary. Likewise, the cell polarity regulators PRKCI and ECT2 were identified as putative drivers of two distinct amplicons on 3q. Co-occurrence analyses suggested potential synergistic or antagonistic relationships between recurrent CNAs. Genes within regions of recurrent CNA showed an enrichment of Cancer Census genes, particularly when filtered for CNA-specific expression.

Conclusion

These analyses provide detailed views of ovarian cancer genomic changes and highlight the benefits of using multiple reference sample types for the evaluation of CNA-specific expression changes.     

Identification of DNA copy-number aberrations by array-comparative genomic hybridization in patients with schizophrenia

Chromosomal abnormalities are implicated as important markers for the pathogenesis in patients with schizophrenia. In this study, with using bacterial artificial chromosome (BAC) array-based comparative genomic hybridization (CGH), we analyzed DNA copy-number changes among 30 patients with schizophrenia. The most frequent changes were partial gain of Xq23 (52%) and loss of 3q13.12 (32%). Other frequent gains were found in: 1p, 6q, 10p, 11p, 11q, 14p, and 15q regions, and frequent losses were found in: 2p, 9q, 10q, 14q, 20q, and 22q regions. The set of abnormal regions was confirmed by real-time PCR (9q12, 9q34.2, 11p15.4, 14q32.33, 15q15.1, 22q11.21, and Xq23). All real-time PCR results were consistent with the array-CGH results. Therefore, it is suggested that array-CGH and real-time PCR analysis could be used as powerful tools in screening for schizophrenia-related genes. Our results might be useful for further exploration of candidate genomic regions in the pathogenesis of schizophrenia.     

Polymorphisms near a chromosome 6q QTL area are associated with modulation of fetal hemoglobin levels in sickle cell anemia.

In patients with sickle cell anemia, fetal hemoglobin (HbF) concentrations vary by 2 orders of magnitude. This variance may be a result of heterogeneity in gene regulatory elements; accordingly, we searched for single nucleotide polymorphisms (SNPs) that might identify this variation. More than 180 SNPs were studied in 38 genes in 280 sickle cell anemia patients. The strongest association with HbF was found with SNPs near a QTL previously localized on chromosome 6q22.3-q23.2. Initially, two SNPs were identified in intergenic portions of this QTL and were associated with about a 20% difference in percent HbF. Subsequently, we genotyped 44 additional SNPs in the genomic region between 136.1 Mb and 137.5 Mb on chromosome 6q. Twelve SNPs, associated with a 20%-30% difference in HbF concentrations, were located in the introns of four genes, PDE7B, MAP7, MAP3K5 and PEX7. In K562 cells, the p38-MAPK pathway has been associated with the activation of gamma-globin gene expression by histone deacetylase inhibitors. Haplotypes C-T-T-T in MAP7 and T-C-C in PEX7 were significantly associated with increases in concentration of HbF, both showing strong dominance. Genetic elements abutting the 6q22.3-q23.2 QTL, may harbor trans-acting elements that help modulate baseline HbF level in sickle cell anemia.     

Detailed Analysis of Focal Chromosome Arm 1q and 6p Amplifications in Urothelial Carcinoma Reveals Complex Genomic Events on 1q,and SOX4 as a Possible Auxiliary Target on 6p

Background

Urothelial carcinoma shows frequent amplifications at 6p22 and 1q21–24. The main target gene at 6p22 is believed to be E2F3, frequently co-amplified with CDKAL1 and SOX4. There are however reports on 6p22 amplifications that do not include E2F3. Previous analyses have identified frequent aberrations occurring at 1q21–24. However, due to complex rearrangements it has been difficult to identify specific 1q21–24 target regions and genes.

Methods

We selected 29 cases with 6p and 37 cases with 1q focal genomic amplifications from 261 cases of urothelial carcinoma analyzed by array-CGH for high resolution zoom-in oligonucleotide array analyses. Genomic analyses were combined with gene expression data and genomic sequence analyses to characterize and fine map 6p22 and 1q21–24 amplifications.

Results

We show that the most frequently amplified gene at 6p22 is SOX4 and that SOX4 can be amplified and overexpressed without the E2F3 or CDKAL1 genes being included in the amplicon. Hence, our data point to SOX4 as an auxiliary amplification target at 6p22. We further show that at least three amplified regions are observed at 1q21–24. Copy number data, combined with gene expression data, highlighted BCL9 and CHD1L as possible targets in the most proximal region and MCL1, SETDB1, and HIF1B as putative targets in the middle region, whereas no obvious targets could be determined in the most distal amplicon. We highlight enrichment of G4 quadruplex sequence motifs and a high number of intraregional sequence duplications, both known to contribute to genomic instability, as prominent features of the 1q21–24 region.

Conclusions

Our detailed analyses of the 6p22 amplicon suggest SOX4 as an auxiliary target gene for amplification. We further demonstrate three separate target regions for amplification at 1q21–24 and identified BCL9, CHD1L, and MCL1, SETDB1, and HIF1B as putative target genes within these regions.     

Dense‐map genome scan for dyslexia supports loci at 4q13, 16p12, 17q22; suggests novel locus at 7q36

Analysis of genetic linkage to dyslexia was performed using 133,165 array‐based SNPs genotyped in 718 persons from 101 dyslexia‐affected families. Results showed five linkage peaks with lod scores >2.3 (4q13.1, 7q36.1‐q36.2, 7q36.3, 16p12.1, and 17q22). Of these five regions, three have been previously implicated in dyslexia (4q13.1, 16p12.1, and 17q22), three have been implicated in attention‐deficit hyperactivity disorder (ADHD, which highly co‐occurs with dyslexia; 4q13.1, 7q36.3, 16p12.1) and four have been implicated in autism (a condition characterized by language deficits; 7q36.1‐q36.2, 7q36.3, 16p12.1, and 17q22). These results highlight the reproducibility of dyslexia linkage signals, even without formally significant lod scores, and suggest dyslexia predisposing genes with relatively major effects and locus heterogeneity. The largest lod score (2.80) occurred at 17q22 within the MSI2 gene, involved in neuronal stem cell lineage proliferation. Interestingly, the 4q13.1 linkage peak (lod 2.34) occurred immediately upstream of the LPHN3 gene, recently reported both linked and associated with ADHD. Separate analyses of larger pedigrees revealed lods >2.3 at 1–3 regions per family; one family showed strong linkage (lod 2.9) to a known dyslexia locus (18p11) not detected in our overall data, demonstrating the value of analyzing single large pedigrees. Association analysis identified no SNPs with genome‐wide significance, although a borderline significant SNP (P = 6 × 10^–7) occurred at 5q35.1 near FGF18, involved in laminar positioning of cortical neurons during development. We conclude that dyslexia genes with relatively major effects exist, are detectable by linkage analysis despite genetic heterogeneity, and show substantial overlapping predisposition with ADHD and autism.