Skip to main content
Download PDF

Copy number loss of (src homology 2 domain containing)-transforming protein 2 (SHC2) gene: discordant loss in monozygotic twins and frequent loss in patients with multiple system atrophy

Molecular Brain volume 4, Article number: 24 (2011) Cite this article

Abstract

Background

Multiple system atrophy (MSA) is a sporadic disease. Its pathogenesis may involve multiple genetic and nongenetic factors, but its etiology remains largely unknown. We hypothesized that the genome of a patient with MSA would demonstrate copy number variations (CNVs) in the genes or genomic regions of interest. To identify genomic alterations increasing the risk for MSA, we examined a pair of monozygotic (MZ) twins discordant for the MSA phenotype and 32 patients with MSA.

Results

By whole-genome CNV analysis using a combination of CNV beadchip and comparative genomic hybridization (CGH)-based CNV microarrays followed by region-targeting, high-density, custom-made oligonucleotide tiling microarray analysis, we identified disease-specific copy number loss of the (Src homology 2 domain containing)-transforming protein 2 (SHC2) gene in the distal 350-kb subtelomeric region of 19p13.3 in the affected MZ twin and 10 of the 31 patients with MSA but not in 2 independent control populations (p = 1.04 × 10-8, odds ratio = 89.8, Pearson's chi-square test).

Conclusions

Copy number loss of SHC2 strongly indicates a causal link to MSA. CNV analysis of phenotypically discordant MZ twins is a powerful tool for identifying disease-predisposing loci. Our results would enable the identification of novel diagnostic measure, therapeutic targets and better understanding of the etiology of MSA.

Background

Multiple system atrophy (MSA; MIM146500) is a progressive neurodegenerative disease clinically characterized by a variable combination of cerebellar ataxia, autonomic disturbance, and parkinsonism with a poor response to levodopa. X-ray computed tomography or magnetic resonance imaging (MRI) studies of the brain usually detect atrophy of the cerebellum and brain stem; MRI also reveals abnormal signal intensity in the white matter of these structures and in the putamen. The neuropathologic features of MSA are neuronal loss, astrogliosis, and argyrophilic glial cytoplasmic inclusions (GCIs) in oligodendrocytes [1]. GCIs involve the aggregation of insoluble fibrillar α-synuclein [2], and the SYNA locus has been associated with MSA in some genetic studies [3, 4], but not in others [5]. The neuroimaging findings, clinical features, and neuropathologic features constitute the current diagnostic criteria for MSA [6]. MSA essentially is a sporadic disorder with onset in adulthood. The involvement of environmental factors and epigenetic mechanisms in its pathogenesis has been postulated; however, its etiology remains largely unknown.

Recently, copy number variations (CNVs) have been recognized as important interindividual structural variations often located in the complex repetitive regions of the human genome. They account for more nucleotide variations between individuals than single-nucleotide polymorphisms (SNPs). Recent studies have indicated that CNVs are considerable contributors to genomic diseases and disease susceptibility in humans [7, 8]. In addition, researchers have developed a molecular strategy that takes advantage of the unique genetic characteristics of monozygotic (MZ) twins: they predicted that analysis of MZ twins discordant for a given disorder (i.e., one is affected and the other is not) will provide important clues for identifying the underlying genetic mechanism of that disorder [9, 10]. Because a genetic difference between MZ twins is an extreme form of somatic mosaicism, de novo CNVs observed in the affected co-twin is a potential characteristic of the CNV region harboring the disease-susceptibility gene(s).

Accordingly, we hypothesized that the genome of a patient with MSA would demonstrate CNVs in the genes or genomic regions of interest. To test this hypothesis, we examined whether a pair of MZ twins discordant for the MSA phenotype display differences in genomic structure as de novo CNVs and searched for CNVs in the genomes of 33 patients with MSA, including MZ twins. We conducted whole-genome CNV analysis by using a genome-wide oligonucleotide CNV microarray and CNV beadchips, both of which focus on the CNV-rich region of the human genome. To validate the results of the whole-genome analysis, we characterized the genetic alteration by using a region-targeting, high-density, custom-made oligonucleotide tiling microarray. Our research strategy is illustrated in Figure 1. Throughout these procedures, we searched for CNVs rare in normal individuals but frequent in patients with MSA. A gene whose transcription is affected by such CNVs and whose expression is mostly confined to the brain must be a crucial candidate for the pathogenesis of MSA.

Figure 1
figure 1

Strategy for identifying genes susceptible to MSA.

Full size image

Methods

Subjects

The study comprised 33 unrelated patients with MSA including a pair of MZ twins discordant for the MSA phenotype (Table 1). All the patients were neurologically evaluated at the Department of Neurology, Hokkaido University Hospital, and examined by brain MRI. The mean (SD) age of MSA onset was 58.1 (8.2) years (age range = 37-72 years). According to the current consensus criteria, 31 patients had probable MSA (10 with MSA-P and 21 with MSA-C) and 2 had possible MSA (one with MSA-P and one with MSA-C) [6]. Of the latter 2 patients, the affected MZ twin (aged 67 years) subsequently satisfied the criteria for probable MSA-C, with obvious atrophy of the brainstem and cerebellum on brain MRI. His disorder began at age 59, but his twin was asymptomatic during this study (Table 2).

Table 1 Profile of the subjects
Full size table
Table 2 Clinical features and phenotypes of the 33 patients
Full size table

One hundred control subjects were randomly selected from community-dwelling elderly individuals with no neurologic diseases (first set of healthy controls; Table 1) [11]. Fisher's exact test showed no significant differences in the mean age or male-to-female ratio between the MSA-affected patients and the first set of control subjects. In addition, chi-square test showed the absence of significant deviation in distribution between the groups. We extracted DNA from peripheral blood leukocytes but not from lymphoblastoid cell lines.

All the subjects gave written informed consent for genetic analysis. The study was approved by the Medical Ethics Committees of Hokkaido University Graduate School of Medicine and Yamagata University Faculty of Medicine.

Whole-genome CNV microarray analysis of the MZ twins

After labeling with Cy5 (test) or Cy3 (reference) dye, the DNAs from the MZ twins discordant for the MSA phenotype were competitively hybridized in a human CNV microarray (SurePrint G3 Human CNV 400K Microarray, Agilent Technologies, Santa Clara, CA) and were washed and scanned according to previously described procedures [12]. Further, the DNAs were separately assayed against a single reference sample (NA19000, a HapMap Japanese male). All the hybridizations were performed in duplicate, and the dyes were interchanged between the DNAs in the second hybridization to eliminate a possible dye-specific bias. Moreover, to validate monozygosity with multiple SNP probes, the DNAs were hybridized with CytoSNP-12 beadchips (Illumina, Inc., San Diego, CA) and scanned for the resulting 300,000 SNP calls.

Whole-genome CNV beadchip analysis of the patients with MSA

We used whole-genome CNV beadchips (57K, i-select format, Illumina Infinium system; deCODE Genetics, Inc., Reykjavik, Iceland) with CNV probes to target the CNV-rich region of the whole genome as previously described [11]. The probe content comprised 15,559 CNV segments covering 190 Mb or 6% of the human genome. The platform has been tested among 4,000 Icelandic and HapMap samples. The test showed that 7,880 of the 42,800 univariant probes called were CNVs with minor allele frequency of > 5%. Over 3,742 CNV segments of size range 5-60 kb were defined; these segments are 18-fold enriched for CNVs compared with the genome-wide average density [13].

High-density custom-made oligonucleotide tiling microarray analysis

We fabricated a custom-made microarray comprising 60-mer probes (Agilent Technologies) targeting a 350-kb genomic region in the distal subtelomeric region of 19p13.3 (Chr. 19: 250,000-600,000 [NCBI Build 36.1, hg18]) and used the Agilent website [12] to select and design the custom tiling array for DNA analysis of the patients with MSA. The tiling microarray experiments were performed as described previously [11]. For the twin analysis, the DNAs from the MZ twins were competitively hybridized with custom tiling microarrays as described for the whole-genome CNV microarray analysis (CNV 400K microarray).

Statistical analysis

Data analysis of the microarray experiments was conducted by using the Aberration Detection Method-2 statistical algorithm (Agilent Technologies) on the basis of the combined log2 ratios at a threshold of 5 as was done in our previous study[14]. The data were centralized, and calls with average log2 ratios of < 0.15 were filtered to exclude false positives.

Data analysis of the CNV beadchip experiments was conducted by using DosageMiner software (deCODE Genetics). The loss/gain analysis consisted of the following 4 steps: (1) intensity normalization and GC content correction, (2) removal of batch effects by using principal component analysis, (3) calling of clusters by using a Gaussian mixture model, and (4) determination of CNV type by using graphical constraints. In brief, CNVs were identified when CNV events were conspicuous among the data, because all sample intensities for CNV probes should increase or decrease relative to those for the neighboring probes (not in the CNV region). To determine deviations in signal intensity, we first normalized the intensities. The normalized intensities for each color channel were determined by using an equation and fit formula (deCODE Genetics). A stretch with more than one marker showing abnormality in copy number compared with a consecutive stretch in the genome is considered more likely evidence of deletion or gain [13]. Statistical analysis of the clinical data was performed with chi-square test or Fisher's exact test in statistical program R [15].

Results

Whole-genome CNV array analysis of MZ twins

According to a previous molecular strategy, we investigated a pair of MZ twins discordant for the MSA phenotype [9, 10]. We confirmed their monozygosity on the basis of > 99% concordant genotypes by SNP beadchip analysis containing more than 300,000 SNPs.

CGH-based whole-genome CNV oligonucleotide microarray analysis for initial screening of genome alterations in the twins revealed putative variations in each twin. DNA assay against a normal Japanese reference sample (HapMap NA1900) confirmed the true genotype of the twins. We identified the following 3 regions with specific copy number loss in the MSA-affected twin (HK33): a CNV region on 2p25.3 (genomic positions 3,642,547-3,643,266, with 2 probes [16]), a CNV region on 4q35.2 (genomic positions 187,590,335-187,594,679, with 32 probes [16]), and a CNV region on 19p13.3 (genomic positions 249,367-252,260, with 35 probes [16]). The patterns of the 3 CNV regions are displayed in Figure S1 in Additional file 1, Figure S2 in Additional file 2, and Figure S3 in Additional file 3 respectively.

The last region on 19p13.3 is located within the 19p subtelomere (250 kb from the 19p telomere). Because a subtelomere is associated with abundant multicopy repeats or CNVs known as rearrangement hotspots, predisposing individuals to deletion or duplication events and possible pathogenic alterations, we hypothesized that this copy number loss on 19p13.3 predisposes individuals to MSA and thus focused on this region for subsequent analysis.

Whole-genome CNV beadchip analysis of patients with MSA

We screened 33 patients with sporadic MSA, including the MZ twins, and a set of 100 controls by whole-genome CNV beadchip analysis (Table 1). Of the markers that showed CNVs, 3 genomic regions (i.e., CNV regions on 4p16.3, 20q13.3, and 19p13.3) displayed frequent copy number loss in 10 patients with MSA but not the controls (p < 1.0 × 10-8).

Among the 3 CNV regions with structural differences between the MZ twins, frequent copy number loss was observed on 19p13.3, but not on 2p25.3 (Figure S4 in Additional file 4) and 4q35.2 (Figure S5 in Additional file 5), in 10 patients with MSA compared with the first set of controls (p = 1.04 × 10-8, odds ratio = 89.8, Pearson's chi-square test; Figure 2). Because the 19p13.3 region (comprising the distal 350-kb region and the proximal 1-Mb region in Figure 2) showed copy number loss in both the affected MZ twin and 10 patients with MSA, we focused on this region for further analysis.

Figure 2
figure 2

Genomic region harboring copy number loss on 19p13.3 in patients with MSA. Data measured by CNV 57K beadchip analysis were analyzed by using the Hidden Malcov Model. The genomic structures of 100 normal control subjects (top) and 33 patients with MSA (bottom) are horizontally aligned from position 000,000 (left) to position 5,000,000 (right). Each blue square represents copy number loss at each CNV probe site whereas each red square represents copy number gain. The horizontal blue bars at the bottom show the distal 350-kb region and proximal 1-Mb region with frequent copy number loss. The asterisks to the bottom left indicate cases of SHC2 deletion. The bottom map shows the positions of putative genes in the region [16].

Full size image

High-density tiling microarray analysis of the MZ twins

To validate the results of the CNV screening in both the affected MZ twin and 10 patients with MSA, we conducted high-density, custom-made oligonucleotide tiling microarray analysis. The DNAs from the twins were compared by competitive hybridization (Figure 3a), the hybridization results were verified by a dye-swap experiment (Figure 3b), and the DNAs were separately hybridized against a normal reference sample (HapMap NA1900; Figure 3c and 3d). We found heterozygous copy number loss (deletion) in the distal 350-kb region on the 19p subtelomere in the MSA-affected twin (HK33). The deleted region encompasses 4 genes including the (Src homology 2 domain containing)-transforming protein 2 (SHC2) gene on 19p13.3 (genomic positions 250-400).

Figure 3
figure 3

High-density custom-made tiling microarray analysis of the MZ twins discordant for the MSA phenotype. (a) Competitive hybridization of genomic DNA from the MSA-affected twin (HK33) versus that from his twin (HK34). (b) Dye-swap experiment of the normal twin (HK34) versus his affected twin (HK33). (c) Competitive hybridization of genomic DNA from the MSA-affected twin (HK33) versus that from a reference Japanese male (HapMap NA1900). (d) Competitive hybridization of genomic DNA from the MSA-affected twin (HK34) versus that from the reference Japanese male (HapMap NA19000). Each blue line represents a moving average ratio of log2 (Cy5/Cy3). The blue regions indicate deletions; these loci were defined by dye-swap experiments (red region). The top map shows the positions of putative genes in the region [16].

Full size image

High-density tiling microarray analysis of the patients with MSA

After confirming the absence of statistical differences in age and gender distribution, we further analyzed the distal 350-kb region showing copy number loss on the 19p subtelomere in the 33 patients with MSA against an additional set of 25 control subjects (second set of controls, Table 1) by using the same custom-made tiling array as for the MZ twins. Again, we found frequent heterozygous copy number loss in this region in 10 of the 33 patients but not in any of the 25 control subjects. Figure 4 demonstrates the moving average pattern of the distal 350-kb region with copy number loss on 19p13.3 in 6 of the 10 patients, and Figure 5 shows a deletion map of this region in the 10 patients.

Figure 4
figure 4

Structure of copy number loss in the 350-kb subtelomeric region on 19p13.3 resolved by high-density tiling microarray. The moving average log2 ratio (y-axis) is plotted against the genomic position along the chromosome (x-axis). HK06, HK12, HK15, HK19, HK23, and HK28 represent 6 patients with MSA and HK60c and HK78c represent controls. The dark lines indicate the copy number loss. The light lines and dotted lines indicate the normal range and median of the average log2 ratios for probes among normal individuals (n = 25), respectively.

Full size image
Figure 5
figure 5

Extent of copy number loss in the 350-kb subtelomeric region on 19p13.3 in the 10 patients with MSA. The horizontal bars represent the length of the copy number loss region between genomic positions 250,000 (left) and 600,000 (right) in each patient. The extent of deletion observed in the affected MZ twin (HK33) is also shown. The top map shows the positions of putative genes in the region [16].

Full size image

Several genes are located in the subtelomeric region on 19p13.3, which were frequently deleted in the patients with MSA and the affected MZ twin (Figure 5). Of these genes, SHC2, which is expressed in the nervous system, is a prime candidate for MSA predisposition, because SHC2 impairment reportedly caused neurologic defects in a mouse model [17, 18]. Therefore, heterozygous copy number loss of SHC2 would predispose individuals to MSA. Other genes in the region, such as hyperpolarization-activated cyclic nucleotide-gated potassium channel 2 (HCN2), mucosal vascular addressin cell adhesion molecule 1 (MADCAM1), and fibroblast growth factor 22 (FGF22) genes, may also have an etiologic link to MSA, because they exert their functions in the nervous system.

Correlation between copy number loss of SHC2 and the MSA phenotype

Finally, we estimated the copy number loss of SHC2 in each patient and tested it for correlation with the phenotype at sampling (MSA-P vs. MSA-C) and at the onset of MSA (parkinsonism vs. cerebellar ataxia with/without autonomic failure) by chi-square test or Fisher's exact probability test. However, no significant correlations were noted (data not shown).

Discussion

We employed 2 strategies for investigating MSA-specific CNVs: examination of de novo CNVs in a pair of MZ twins discordant for the MSA phenotype and CNV analysis of 33 patients with MSA, including the MZ twins. The CNV beadchip analysis enabled direct measurement of the fluorescent intensity of each marker, which is suitable for screening sporadic cases and control populations. The comparative genomic hybridization (CGH)-based whole-genome CNV oligonucleotide microarray enabled comparison of 2 similar genomes, which is ideal for detecting subtle differences between MZ twins. The combination of CGH-based genome-wide oligonucleotide CNV microarray and CNV beadchip analysis followed by region-targeting, high-density, custom-made oligonucleotide tiling microarray analysis led us to identify copy number loss of SHC2 and its neighboring genes in the 19p13.3 subtelomeric region of the affected MZ twin and 10 (30%) unrelated patients with MSA. Despite the limited number of subjects, this frequency is quite high when compared with those of other disorders: in schizophrenia, for instance, less than 1% in the patients has pathogenic CNV markers [13]. The copy number loss of SHC2 in the MSA-affected twin, which was validated by the population analysis, strongly suggests a causal link between genomic alteration and MSA phenotype.

Our results support the idea that CNV analysis of phenotypically discordant MZ twins is a powerful tool for identifying disease-predisposing loci. They also support the opinion that molecular analysis of phenotypically discordant MZ twins is an excellent focus for studying disease, because genotypic differences between twins derived from the same zygote highlight somatic variation [9]. Instances of somatic mosaicism by mutation in specific genes or chromosomal aberrations with links to disease have been described previously [19, 20]. Our results also suggest that somatic mosaicism by pathogenic mutations affecting disease-susceptibility genes is observable frequently, rather than as an exception [9]. Further studies of a larger cohort of MZ twins discordant for MSA will be particularly useful for more detailed characterization of the genetic factors predisposing individuals to MSA.

We identified frequent heterozygous copy number loss of SHC2 and the surrounding area on the 19p13.3 subtelomeric region. This region contains a substantially larger number of low copy repeats and segmental duplications in a small area than that in the whole human genome. This complicated repetitive structure has prevented accurate mapping and sequencing analysis of the region [21]. The presence of numerous low copy repeats results in instability and can trigger frequent segmental loss by unequal crossover or end-joining events; therefore, such unstable CNV-rich subtelomeric regions are predisposed to deletion or duplication events [22]. Noteworthily, locus-specific mutation rates for CNVs or structural rearrangements are between 10-6 and 10-4, at least 2 to 4 orders of magnitude (100- to 10,000-fold) greater than those for point mutations or SNPs [23]. Considering the causes of sporadic disease such as MSA, the new high mutation rate of CNVs is fascinating, and should be studied while searching for gene(s) related to such disorders [8]. CNVs in the distal 350-kb region of the 19p13.3 subtelomere have not been described for other neurologic diseases. For instance, our preliminary CNV beadchip analysis did not detect specific copy number loss in patients with multiple sclerosis (unpublished data).

SHC2 was deleted in the affected MZ twin and frequently deleted in the patients with MSA. SHC2 mRNA is expressed in various human adult tissues, including the nervous system. None of the CNVs around SHC2 described in the Database of Genomic Variants were associated with human neurologic disease [16]. In adult mice, Shc2 expression is limited to the nervous system, and in rat and mouse embryos, the expression of Shc2 mRNA is the highest in the dorsal root ganglia and superior cervical ganglia [17, 18]. Shc proteins act as molecular switches in neuronal cell development from proliferation to survival and/or differentiation.

Several other genes frequently deleted in patients with MSA, such as HCN2, MADCAM1, and FGF22, are also expressed in nervous tissue. HCN2 contributes to spontaneous rhythmic activity in the external segment of the globus pallidus [24, 25]. Augmented currents through the channel of an HCN2 variant was described in patients with febrile seizure syndrome [26]. Human MADCAM-1 mRNA transcripts are expressed in the brain and abnormal elevation of MADCAM-1 has been described in patients with multiple sclerosis [27]. FGF22 plays an essential role in nervous system development and is expressed by cerebellar granule cells [28]. Copy number loss of these genes may alter their expression or produce aberrant or unstable mRNA or protein products. Furthermore, CNVs influence the expression of not only the genes harboring them but also the genes in their vicinity extending up to half a megabase [29]. The mechanism of the influence of CNVs on the phenotype over multiple gene expressions is yet to be studied.

Specific copy number loss indicates the possibility that MSA is a genomic disorder derived from genetic hotspots harboring unstable structures predisposed to CNVs. In this regard, rare clustering of MSA within families could provide information on inherited genetic variants related to the development of this disorder [30].

Conclusions

Our study provides compelling evidence of heterozygous copy number loss in the SHC2 region in the affected MZ twin and one-third of the patients with MSA. Further studies on the function of SHC2 in the nervous system help to elucidate the pathogenesis of MSA and identify novel therapeutic targets for the disease. Studies of independent populations with MSA and different ethnicities are warranted to fully understand the etiology of MSA, including yet unidentified genetic and/or nongenetic factors.

Abbreviations

BSG:

basigin

CDC34:

cell division cycle 34 homolog (S. cerevisiae)

CNV:

copy number variation

CGH:

comparative genomic hybridization

C2CD4C:

C2 calcium-dependent domain containing 4C

FGF22:

fibroblast growth factor 22

GCI:

glial cytoplasmic inclusion

GZMM:

granzyme M

HCN2:

hyperpolarization-activated cyclic nucleotide-gated potassium channel 2

MADCAM1:

mucosal vascular addressin cell adhesion molecule 1

MIER2:

mesoderm induction early response 1, family member 2

MRI:

magnetic resonance imaging

MSA:

multiple system atrophy

MZ:

monozygotic

POLRMT:

polymerase (RNA) mitochondrial (DNA directed)

SHC2:

(Src homology 2 domain containing)-transforming protein 2

SNP:

single-nucleotide polymorphism

THEG:

testicular haploid expressed gene homolog (mouse)

References

  1. Papp MI, Lantos PL: The distribution of oligodendroglial inclusions in multiple system atrophy and its relevance to clinical symptomatology. Brain. 1994, 117: 235-243. 10.1093/brain/117.2.235.

    Article  PubMed  Google Scholar 

  2. Tu P, Galvin JE, Baba M, Giasson B, Tomita T, Leight S, Nakajo S, Iwatsubo T, Trojanowski JQ, Lee VMY: Glial cytoplasmic inclusions in white matter oligodendrocytes of multiple system atrophy brains contain insoluble α-synuclein. Ann Neurol. 1998, 44: 415-422. 10.1002/ana.410440324.

    Article  CAS  PubMed  Google Scholar 

  3. Scholz SW, Houlden H, Schulte C, Sharma M, Li A, Berg D, Melchers A, Paudel R, Gibbs JR, Simon-Sanchez J, Paisan-Ruiz C, Bras J, Ding J, Chen H, Traynor BJ, Arepalli S, Zonozi RR, Revesz T, Holton J, Wood N, Lees A, Oertel W, Wüllner U, Goldwurm S, Pellecchia MT, Illig T, Riess O, Fernandez HH, Rodriguez RL, Okun MS, Poewe W, Wenning GK, Hardy JA, Singleton AB, Del Sorbo F, Schneider S, Bhatia KP, Gasser T: SNCA variants are associated with increased risk for multiple system atrophy. Ann Neurol. 2009, 65: 610-614. 10.1002/ana.21685.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Al-Chalabi A, Dürr A, Wood NW, Parkinson MH, Camuzat A, Hulot JS, Morrison KE, Renton A, Sussmuth SD, Landwehrmeyer BG, Ludolph A, Agid Y, Brice A, Leigh PN, Bensimon G, NNIPPS Genetic Study Group: Genetic variants of the α-synuclein gene SNCA are associated with multiple system atrophy. PLoS One. 2009, 4: e7114-10.1371/journal.pone.0007114.

    Article  PubMed Central  PubMed  Google Scholar 

  5. Morris HR, Vaughan JR, Datta SR, Bandopadhyay R, Rohan De Silva HA, Schrag A, Cairns NJ, Burn D, Nath U, Lantos PL, Daniel S, Lees AJ, Quinn NP, Wood NW: Multiple system atrophy/progressive supranuclear palsy: α-synuclein, synphilin, tau, and APOE. Neurology. 2000, 55: 1918-1920.

    Article  CAS  PubMed  Google Scholar 

  6. Gilman S, Wenning GK, Low PA, Brooks DJ, Mathias CJ, Trojanowski JQ, Wood NW, Colosimo C, Dürr A, Fowler CJ, Kaufmann H, Klockgether T, Lees A, Poewe W, Quinn N, Revesz T, Robertson D, Sandroni P, Seppi K, Vidailhet M: Second consensus statement on the diagnosis of multiple system atrophy. Neurology. 2008, 71: 670-676. 10.1212/01.wnl.0000324625.00404.15.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Conrad DF, Pinto D, Redon R, Feuk L, Gokcumen O, Zhang Y, Aerts J, Andrews TD, Barnes C, Campbell P, Fitzgerald T, Hu M, Ihm CH, Kristiansson K, Macarthur DG, Macdonald JR, Onyiah I, Pang AW, Robson S, Stirrups K, Valsesia A, Walter K, Wei J, Wellcome Trust Case Control Consortium, Tyler-Smith C, Carter NP, Lee C, Scherer SW, Hurles ME: Origins and functional impact of copy number variation in the human genome. Nature. 2010, 464: 704-712. 10.1038/nature08516.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Eichler EE, Flint J, Gibson G, Kong A, Leal SM, Moore JH, Nadeau JH: Missing heritability and strategies for finding the underlying causes of complex disease. Nat Rev Genet. 2010, 11: 446-450. 10.1038/nrg2809.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Bruder CEG, Piotrowski A, Gijsbers AACJ, Andersson R, Erickson S, Diaz de Ståhl T, Menzel U, Sandgren J, von Tell D, Poplawski A, Crowley M, Crasto C, Partridge EC, Tiwari H, Allison DB, Komorowski J, van Ommen G-JB, Boomsma DI, Pedersen NL, den Dunnen JT, Wirdefeldt K, Dumanski JP: Phenotypically concordant and discordant monozygotic twins display different DNA copy-number-variation profiles. Am J Hum Genet. 2008, 82: 763-771. 10.1016/j.ajhg.2007.12.011.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Baranzini SE, Mudge J, van Velkinburgh JC, Khankhanian P, Khrebtukova I, Miller NA, Zhang L, Farmer AD, Bell CJ, Kim RW, May GD, Woodward JE, Caillier SJ, McElroy JP, Gomez R, Pando MJ, Clendenen LE, Ganusova EE, Schilkey FD, Ramaraj T, Khan OA, Huntley JJ, Luo S, Kwok P, Wu TD, Schroth GP, Oksenberg JR, Hauser SL, Kingsmore SF: Genome, epigenome and RNA sequences of monozygotic twins discordant for multiple sclerosis. Nature. 2010, 464: 1351-1356. 10.1038/nature08990.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Kato T, Emi M, Sato H, Arawaka S, Wada M, Kawanami T, Katagiri T, Tsuburaya K, Toyoshima I, Tanaka F, Sobue G, Matsubara K: Segmental copy-number gain within the region of isopentenyl diphosphate isomerase genes in sporadic amyotrophic lateral sclerosis. Biochem Biophys Res Commun. 2010, 402: 438-442. 10.1016/j.bbrc.2010.10.056.

    Article  CAS  PubMed  Google Scholar 

  12. Agilent Technologies eArray. [http://earray.chem.agilent.com/earray/]

  13. Stefansson H, Rujescu D, Cichon S, Pietiläinen OP, Ingason A, Steinberg S, Fossdal R, Sigurdsson E, Sigmundsson T, Buizer-Voskamp JE, Hansen T, Jakobsen KD, Muglia P, Francks C, Matthews PM, Gylfason A, Halldorsson BV, Gudbjartsson D, Thorgeirsson TE, Sigurdsson A, Jonasdottir A, Jonasdottir A, Bjornsson A, Mattiasdottir S, Blondal T, Haraldsson M, Magnusdottir BB, Giegling I, Möller HJ, Hartmann A, et al: Large recurrent microdeletions associated with schizophrenia. Nature. 2008, 455: 232-236. 10.1038/nature07229.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Nakayama M, Nozu K, Goto Y, Kamei K, Ito S, Sato H, Emi M, Nakanishi K, Tsuchiya S, Iijima K: HNF1B alterations associated with congenital anomalies of the kidney and urinary tract. Pediatr Nephrol. 2010, 25: 1073-1079. 10.1007/s00467-010-1454-9.

    Article  PubMed  Google Scholar 

  15. R Development Core Team: The R Project for Statistical Computing. [http://www.R-project.org/]

  16. Database of Genomic Variants. [http://projects.tcag.ca/variation/]

  17. Nakamura T, Muraoka S, Sanokawa R, Mori N: N-Shc and Sck, two neuronally expressed Shc adapter homologs. Their differential regional expression in the brain and roles in neurotrophin and Src signaling. J Biol Chem. 1998, 273: 6960-6967. 10.1074/jbc.273.12.6960.

    Article  CAS  PubMed  Google Scholar 

  18. Sakai R, Henderson JT, O'Bryan JP, Elia AJ, Saxton TM, Pawson T: The mammalian ShcB and ShcC phosphotyrosine docking proteins function in the maturation of sensory and sympathetic neurons. Neuron. 2000, 28: 819-833. 10.1016/S0896-6273(00)00156-2.

    Article  CAS  PubMed  Google Scholar 

  19. Youssoufian H, Pyeritz RE: Mechanisms and consequences of somatic mosaicism in humans. Nat Rev Genet. 2002, 3: 748-758.

    Article  CAS  PubMed  Google Scholar 

  20. Erickson RP: Somatic gene mutation and human disease other than cancer. Mutat Res. 2003, 543: 125-136. 10.1016/S1383-5742(03)00010-3.

    Article  CAS  PubMed  Google Scholar 

  21. Itsara A, Cooper GM, Baker C, Girirajan S, Li J, Absher D, Krauss RM, Myers RM, Ridker PM, Chasman DI, Mefford H, Ying P, Nickerson DA, Eichler EE: Population analysis of large copy number variants and hotspots of human genetic disease. Am J Hum Genet. 2009, 84: 148-161. 10.1016/j.ajhg.2008.12.014.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Sharp AJ, Locke DP, McGrath SD, Cheng Z, Bailey JA, Vallente RU, Pertz LM, Clark RA, Schwartz S, Segraves R, Oseroff VV, Albertson DG, Pinkel D, Eichler EE: Segmental duplications and copy-number variation in the human genome. Am J Hum Genet. 2005, 77: 78-88. 10.1086/431652.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Lupski JR: Genomic rearrangements and sporadic disease. Nat Genet. 2007, 39: S43-47. 10.1038/ng2084.

    Article  CAS  PubMed  Google Scholar 

  24. Santoro B, Grant SGN, Bartsch D, Kandel ER: Interactive cloning with the SH3 domain of N-src identifies a new brain specific ion channel protein, with homology to Eag and cyclic nucleotide-gated channels. Proc Natl Acad Sci USA. 1997, 94: 14815-14820. 10.1073/pnas.94.26.14815.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Ludwig A, Zong X, Stieber J, Hullin R, Hofmann F, Biel M: Two pacemaker channels from human heart with profoundly different activation kinetics. EMBO J. 1999, 18: 2323-2329. 10.1093/emboj/18.9.2323.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Dibbens LM, Reid CA, Hodgson B, Thomas EA, Phillips AM, Gazina E, Cromer BA, Clarke AL, Baram TZ, Scheffer IE, Berkovic SF, Petrou S: Augmented currents of an HCN2 variant in patients with febrile seizure syndromes. Ann Neurol. 2010, 67: 542-546.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Allavena R, Noy S, Andrews M, Pullen N: CNS elevation of vascular and not mucosal addressin cell adhesion molecules in patients with multiple sclerosis. Am J Pathol. 2010, 176: 556-562. 10.2353/ajpath.2010.090437.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Umemori H, Linhoff MW, Ornitz DM, Sanes JR: FGF22 and its close relatives are presynaptic organizing molecules in the mammalian brain. Cell. 2004, 118: 257-270. 10.1016/j.cell.2004.06.025.

    Article  CAS  PubMed  Google Scholar 

  29. Henrichsen CN, Vinckenbosch N, Zöllner S, Chaignat E, Pradervand S, Schütz F, Ruedi M, Kaessmann H, Reymond A: Segmental copy number variation shapes tissue transcriptomes. Nat Genet. 2009, 41: 424-429. 10.1038/ng.345.

    Article  CAS  PubMed  Google Scholar 

  30. Hara K, Momose Y, Tokiguchi S, Shimohata M, Terajima K, Onodera O, Kakita A, Yamada M, Takahashi H, Hirasawa M, Mizuno Y, Ogata K, Goto J, Kanazawa I, Nishizawa M, Tsuji S: Multiplex families with multiple system atrophy. Arch Neurol. 2007, 64: 545-551. 10.1001/archneur.64.4.545.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors thank Y. Hama, S. Shimizu, and other staff of the Department of Neurology, Graduate School of Medicine, Hokkaido University, for their support in the clinical research and clerical assistance as well as the patients and their families, who enabled this research.

The study was supported in part by a Grant-in-Aid for Scientific Research from the Japanese Ministry of Education, Culture, Sports, Science and Technology (to H. Sasaki); a Grant-in-Aid from the Research Committee for Ataxic Diseases, Japanese Ministry of Health, Labour and Welfare (to H. Sasaki); a Grant-in-Aid from the Matching Program for Innovations in Future Drug Discovery and Medical Care of Japan (to H. Sasaki); and a Grant-in-Aid from the Global COE Program (F03) of the Japan Society for the Promotion of Science (to TK).

Author information

Author notes

    Authors and Affiliations

    1. Department of Neurology, Graduate School of Medicine, Hokkaido University, North 15, West 7, Kita-ku, Sapporo, 060-8638, Japan

      Hidenao Sasaki & Ichiro Yabe

    2. CNV Laboratory, DNA Chip Research Institute, 1-1-43 Suehirocho, Tsurumi-ku, Yokohama, Kanagawa, 230-0045, Japan

      Mitsuru Emi, Hiroshi Iijima, Noriko Ito, Hidenori Sato & Kenichi Matsubara

    3. Department of Neurology, Haematology, Metabolism, Endocrinology, and Diabetology, Yamagata University Faculty of Medicine, 2-2-2 Iida-Nishi, Yamagata, Yamagata, 990-9585, Japan

      Takeo Kato

    4. Creative Research Institution, Hokkaido University, North 21, West 10, Kita-ku, Sapporo, 001-0021, Japan

      Jun Utsumi

    5. Graduate School of Pharmaceutical Sciences, Kyoto University, 46-29 Yoshida-Shimo-Adachi-cho, Sakyo-ku, Kyoto, 606-8501, Japan

      Jun Utsumi

    Authors
    1. Hidenao Sasaki
      View author publications

      You can also search for this author in PubMed Google Scholar

    2. Mitsuru Emi
      View author publications

      You can also search for this author in PubMed Google Scholar

    3. Hiroshi Iijima
      View author publications

      You can also search for this author in PubMed Google Scholar

    4. Noriko Ito
      View author publications

      You can also search for this author in PubMed Google Scholar

    5. Hidenori Sato
      View author publications

      You can also search for this author in PubMed Google Scholar

    6. Ichiro Yabe
      View author publications

      You can also search for this author in PubMed Google Scholar

    7. Takeo Kato
      View author publications

      You can also search for this author in PubMed Google Scholar

    8. Jun Utsumi
      View author publications

      You can also search for this author in PubMed Google Scholar

    9. Kenichi Matsubara
      View author publications

      You can also search for this author in PubMed Google Scholar

    Corresponding author

    Correspondence to Hidenao Sasaki.

    Additional information

    Competing interests

    The authors declare that they have no competing interests.

    Authors' contributions

    H. Sasaki and ME conceived the study, conducted the molecular genetic studies and data interpretation, had full access to the data, accept responsibility for the integrity of the work and accuracy of the data analysis, and drafted the manuscript as well as provided administrative, technical, and material support. HI conducted the molecular genetic studies and data interpretation, participated in designing the study, and performed the statistical analysis. NI and H. Sato conducted the molecular genetic studies and data interpretation. IY participated in designing the study and performed the statistical analysis as well as provided administrative, technical, and material support. TK conducted the molecular genetic studies and data interpretation as well as provided administrative, technical, and material support. KM and JU participated in supervising the study and critically revising the manuscript content. All authors read and approved the final manuscript.

    Hidenao Sasaki, Mitsuru Emi contributed equally to this work.

    Electronic supplementary material

    13041_2011_128_MOESM1_ESM.PDF

    Additional file 1: Figure S1 - Pattern of the CNV region on 2p25.3 in the MZ twins discordant for the MSA phenotype by CGH-based whole-genome CNV microarray analysis. (Top panel) Competitive hybridization of genomic DNA from the MSA-affected twin (HK33) versus that from his twin (HK34). (Bottom panel) Dye-swap experiment of the normal twin (HK34) versus his affected twin (HK33). Each blue line represents a moving average ratio of log2 (Cy5/Cy3). The blue region indicates deletion. These loci were defined by dye-swap experiments (red region). (PDF 43 KB)

    13041_2011_128_MOESM2_ESM.PDF

    Additional file 2: Figure S2 - Pattern of the CNV region on 4q35.2 in the MZ twins discordant for the MSA phenotype by CGH-based whole-genome CNV microarray analysis. (Top panel) Competitive hybridization of genomic DNA from the MSA-affected twin (HK33) versus that from his twin (HK34). (Bottom panel) Dye-swap experiment of the normal twin (HK34) versus his affected twin (HK33). Each blue line represents a moving average ratio of log2 (Cy5/Cy3). The blue region indicates deletion. These loci were defined by dye-swap experiments (red region). (PDF 36 KB)

    13041_2011_128_MOESM3_ESM.PDF

    Additional file 3: Figure S3 - Pattern of the CNV region on 19p13.3 in the MZ twins discordant for the MSA phenotype by CGH-based whole-genome CNV microarray analysis. (Top panel) Competitive hybridization of genomic DNA from the MSA-affected twin (HK33) versus that from his twin (HK34). (Bottom panel) Dye-swap experiment of the normal twin (HK34) versus his affected twin (HK33). Each blue line represents a moving average ratio of log2 (Cy5/Cy3). The blue region indicates deletion. These loci were defined by dye-swap experiments (red region). (PDF 47 KB)

    13041_2011_128_MOESM4_ESM.PDF

    Additional file 4: Figure S4 - Pattern of the CNV region on 2p25.3 in the patients with MSA and controls. Data measured by CNV 57K beadchip analysis were analyzed by the Hidden Malcov Model. The genomic structures of 100 normal control subjects (top) and 33 patients with MSA (bottom) are horizontally aligned from position 000,000 (left) to position 5,000,000 (right). Each blue square represents copy number loss at each CNV probe site whereas each red square represents copy number gain. (PDF 458 KB)

    13041_2011_128_MOESM5_ESM.PDF

    Additional file 5: Figure S5 - Pattern of the CNV region on 4q35.2 in the patients with MSA and controls. Data measured by CNV 57K beadchip analysis were analyzed by the Hidden Malcov Model. The genomic structures of 100 normal control subjects (top) and 33 patients with MSA (bottom) are horizontally aligned from position 186,000,000 (left) to position 191,000,000 (right). Each blue square represents copy number loss at each CNV probe site whereas each red square represents copy number gain. (PDF 472 KB)

    Authors’ original submitted files for images

    Rights and permissions

    This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

    Reprints and Permissions

    About this article

    Cite this article

    Sasaki, H., Emi, M., Iijima, H. et al. Copy number loss of (src homology 2 domain containing)-transforming protein 2 (SHC2) gene: discordant loss in monozygotic twins and frequent loss in patients with multiple system atrophy. Mol Brain 4, 24 (2011). https://doi.org/10.1186/1756-6606-4-24

    Download citation

    • Received:

    • Accepted:

    • Published:

    • DOI: https://doi.org/10.1186/1756-6606-4-24

    Keywords

    • Multiple system atrophy
    • copy number variation
    • phenotypically discordant monozygotic twins
    • (Src homology 2 domain containing)-transforming protein 2
    • subtelomere
    • ataxia
    • parkinsonism
    • disease-susceptibility gene

    Molecular Brain

    ISSN: 1756-6606