Open Access
Open Peer Review

This article has Open Peer Review reports available.

How does Open Peer Review work?

Novel VPS13B Mutations in Three Large Pakistani Cohen Syndrome Families Suggests a Baloch Variant with Autistic-Like Features

  • Muhammad Arshad Rafiq1, 20,
  • Claire S Leblond2, 21,
  • Muhammad Arif Nadeem Saqib3,
  • Akshita K. Vincent1,
  • Amirthagowri Ambalavanan2, 21,
  • Falak Sher Khan3,
  • Muhammad Ayaz4,
  • Naseema Shaheen5,
  • Dan Spiegelman2, 21,
  • Ghazanfar Ali6,
  • Muhammad Amin-ud-din7,
  • Sandra Laurent2, 21,
  • Huda Mahmood1,
  • Mehtab Christian2,
  • Nadir Ali3,
  • Alanna Fennell1,
  • Zohair Nanjiani8,
  • Gerald Egger1, 9,
  • Chantal Caron10, 11,
  • Ahmed Waqas3,
  • Muhammad Ayub4, 12,
  • Saima Rasheed13,
  • Baudouin Forgeot d’Arc10, 11, 14,
  • Amelie Johnson11, 14,
  • Joyce So15, 16, 17,
  • Muhammad Qasim Brohi18,
  • Laurent Mottron10, 11, 14,
  • Muhammad Ansar3,
  • John B Vincent1, 19Email author and
  • Lan Xiong2, 11, 14Email author
Contributed equally
BMC Medical Genetics201516:41

https://doi.org/10.1186/s12881-015-0183-0

Received: 3 December 2014

Accepted: 29 May 2015

Published: 25 June 2015

Abstract

Background

Cohen Syndrome (COH1) is a rare autosomal recessive disorder, principally identified by ocular, neural and muscular deficits. We identified three large consanguineous Pakistani families with intellectual disability and in some cases with autistic traits.

Methods

Clinical assessments were performed in order to allow comparison of clinical features with other VPS13B mutations. Homozygosity mapping followed by whole exome sequencing and Sanger sequencing strategies were used to identify disease-related mutations.

Results

We identified two novel homozygous deletion mutations in VPS13B, firstly a 1 bp deletion, NM_017890.4:c.6879delT; p.Phe2293Leufs*24, and secondly a deletion of exons 37-40, which co-segregate with affected status. In addition to COH1-related traits, autistic features were reported in a number of family members, contrasting with the “friendly” demeanour often associated with COH1. The c.6879delT mutation is present in two families from different regions of the country, but both from the Baloch sub-ethnic group, and with a shared haplotype, indicating a founder effect among the Baloch population.

Conclusion

We suspect that the c.6879delT mutation may be a common cause of COH1 and similar phenotypes among the Baloch population. Additionally, most of the individuals with the c.6879delT mutation in these two families also present with autistic like traits, and suggests that this variant may lead to a distinct autistic-like COH1 subgroup.

Background

Intellectual disability (ID) affects ~1 % of the population [1], however for up to 80 % of ID patients the underlying etiology is unknown. Chromosomal aberrations are thought to be the most common cause [2]. In a cohort of 1170 unselected patients referred for developmental delay or ID, Cohen syndrome (COH1) [MIM# 216550] was found among several common genetic diagnoses for ID; and with a frequency of 0.7 %, COH1 is ranked right after Fragile-X syndrome (1.2 %) among patients with clear diagnosis of ID [2]. COH1 is an autosomal recessively inherited ID disorder, first described by Cohen and colleagues in 1973 [3]. Kivitie-Kallio and Norio [4] proposed that the essential features for a COH1 diagnosis include non-progressive psychomotor retardation, motor clumsiness, microcephaly, distinctive facial features, childhood hypotonia, joint laxity, chorioretinal dystrophy, myopia and granulocytopenia [4]. After much debate on these diagnostic criteria in ethnically diverse groups, newly amended diagnostic guidelines were given in which a child with significant learning difficulties must show at least two of the three features: typical facial gestalt, pigmentary retinopathy or neutropenia [5, 6].

To date, COH1 has been mainly attributed to mutations in the VPS13B gene (MIM# 607817) among patients from many different ethnic backgrounds [715].

For approximately 70 % of COH1 patients both mutated alleles of VPS13B were detected, whereas in approximately 20-30 % of the cases a single heterozygous mutation was found [8, 10]. A total of 159 VPS13B gene mutations (HGMD, 2013.4) have been identified which have since led to further modification of the clinical diagnostic criteria for differentiation between Cohen and Cohen-like syndrome [8, 14, 1618]. In the present study, we describe the ascertainment and clinical assessment of three large consanguineous families from Pakistan segregating ID in an autosomal recessive pattern. We also describe the mapping by autozygosity of the disease locus for each family to 8q22, followed by screening and mutation detection in VPS13B through either candidate gene sequencing or whole exome sequencing.

Methods

Patients and families

Two of the Pakistani families (RQMR10 and ANMR51) were initially identified through referral to the study on genetics of autosomal recessive ID through collaborations between the Lahore Institute of Research & Development (LIRD), Pakistan, and the Centre for Addiction & Mental Health (CAMH), Toronto, Canada, and between Quaid-i-Azam University, Islamabad and CAMH. These two families are from Dera Ghazi Khan District, Punjab Province, and from Naushahro Feroze district in Sindh Province, respectively. The third Pakistani family (ATM02) was identified and recruited as a recessive form of pervasive developmental disorder (PDD) from Dadu District, Sindh Province, through collaboration between the Sir Cowasji Jehangir Institute of Psychiatry (SCJIP), Hyderabad, Pakistan and CHUM Research Center (CRCHUM), University of Montreal, Canada. Research Ethics Board approval was obtained from the involved institutions (LIRD, Quaid-i-Azam University, CAMH, SCJIP and CRCHUM), and written informed consent was obtained from the parents/guardians and other family members for publication of this article and any accompanying images, in compliance with the Helsinki Declaration. The pedigrees (Fig. 1) were constructed by interviewing different family members. 293 healthy Pakistani controls were also recruited, after written informed consent was given, so that variants could be checked for frequency in an ethnically matched population.
Fig. 1

Pedigrees of the three Pakistani families. For family RQMR10, although the degree of relatedness between the parents of V-1, V-2, and V-3 could not be established, marriages within this rural community are strictly within the clan system, and thus they are almost certainly related. ATM02 is part of a much larger pedigree segregating ID, with 18 affected individuals, however only members of the portion indicated here was assessed in this study. This full pedigree is shown in Additional files 1, 2, 3 and 4

Clinical investigation of families RQMR10 and ANMR51

All affected family members were examined clinically, with special attention to neurological, psychological, ophthalmological, skeletal, behavioral and dental anomalies. One affected member from RQMR10 family (VI-2) was subjected to fundoscopic examination, in order to investigate the ophthalmological condition associated with the disease in that family. Complete biochemical testing of blood of affected as well as several unaffected individuals was performed for both families, in particular to confirm or exclude neutropenia — one of the characteristic features of COH1.

Clinical investigation of ATM02 family

The three affected brothers were examined by local clinicians, including a psychiatrist, a clinical psychologist and a neurologist, and evaluated retrospectively using the Childhood Autism Rating Scale (CARS) [19]. One of the affected individuals was also available for additional clinical tests: including routine blood/urine test, biochemistry test (serum TSH, B12 and folic acid, CPK, and random blood glucose test) and brain MRI. Their clinical synopsis, CARS ratings, lab test results and diagnosis have been further reviewed and validated by three pediatric psychiatrists (C.C., B.F.A., L.M.) from Montreal.

DNA extraction, genotyping, homozygosity mapping and sequencing

Genomic DNA was extracted from either whole blood or from saliva of affected and unaffected family members, by following standard protocols. For RQMR10, DNA samples from three affected (V-2, VI-2, VI-3) and two unaffected individuals (V-3, V-4), were genotyped using the Affymetrix Genome-Wide Human SNP 6.0 microarray at The Centre for Applied Genomics (TCAG), Toronto. DNA samples from eight affected individuals (IV-3, IV-4, IV-5, V-1, V-3, V-6, V-7, V-8) and one unaffected individual (V-10) from ANMR51 were genotyped using the AffymetrixGeneChip Mapping 500 K NspI array, at the London Regional Genomics Centre (LRGC, University of Western Ontario). For ATM02, the three affected individuals (V-1, V-3, V-5) and four unaffected individuals (IV-7, IV-8, V-4 and V-6) were genotyped using the Illumina HumanOmniExpress BeadChip at the Genome Quebec and McGill University Innovation Centre.

For RQMR10 and ANMR51 homozygosity-by-descent (HBD) or autozygosity mapping and copy number variation (CNV) analysis was performed using dChip software [20, 21]. For ATM02, KaryoStudio program from Illumina Inc. was used to exclude cytogenetic abnormalities, and HBD mapping was performed by using the HomozygosityMapper [22].

For RQMR10 and ANMR51, direct Sanger sequencing of VPS13B was performed, based on homozygosity mapping results, for all 62 exons and splice junctions of VPS13B. PCR primers were designed using Primer3 (version 0.4.0). PCR amplicons were sequenced at TCAG. Sequencing data were aligned to the annotated human genome sequence (hg19) using the BLAT tool of the UCSC Genome Browser, and variants confirmed by checking the sequencing chromatograms. Once identified, the inheritance pattern of the variant was checked across the pedigree for cosegregation with the phenotype. Potential candidate variants were further validated in Pakistani population controls through Sanger sequencing.

For ATM02 family, whole exome capturing and sequencing in two affected individuals (V-1 and V-5) were performed using Agilent SureSelect™ 50 Mb Human All Exon Kit and sequenced by HiSeq2000 at the Genome Quebec and McGill University Innovation Centre. Sequence data was processed and aligned using a bioinformatic pipeline including FASTQC, BWA [23], GATK [24] and ANNOVAR programs [25]. An in-house script was created to generate familial variant segregation files. Potential candidate variants were further validated in additional family members, and Pakistani population controls through Sanger sequencing.

Deletion break-point mapping

A PCR-based approach was used in order to map the deletion breakpoints in family ANMR51. Sets of primers were designed using Primer 3 (Primer3 software, version 0.4.0) for regions flanking the putative deleted exons (exons 37–40). PCR amplification was performed using a forward primer from the closest proximal amplicon and a reverse primer from the closest distal set that that amplified in affected individuals, in order to identify the junction fragment. The junction fragment was then sequenced, and the exact physical co-ordinates were mapped using the BLAT tool of the UCSC Genome Browser.

Results

Clinical findings

Families RQMR10 and ANMR51 (Figure 1): All affected, as well as several unaffected family members were examined clinically, with special attention to presence or absence of microcephaly, short stature, speech development, facial dysmorphisms, myopia, pigmentary retinopathy, pes planus, pes cavus, neutropenia and neuropsychiatric presentation, except for RQMR10 VI-1, who died of liver failure after joining the study, results shown in Table 1. Analysis of photographs (Fig. 2) by a dysmorphologist (J.S.) indicated for RQMR10 affected individuals having a boxy shaped face, straight eyebrows, mild synophrys, hypoplastic nares, prominent lips, and a small chin with mild prognathia (V-2, VI-3, VI-2), left eye exotropia (V-2, VI-3), ptosis (VI-2), slender nose with hypoplastic nares and a hanging columella (V-2). VI-2 also shows patchy skin hyperpigmentation over the temples, cheeks, nose, chin, and lateral aspects of the face and neck. For family ANMR51 individuals V-6 and V-7 have straight, bushy eyebrows with synophrys, deep-set eyes, short philtrum, and thin lips. V-6 also has a large chin, prominent nasal bridge with hypoplastic nares, pointed nasal tip and hanging columella. V-7 shows bilateral ptosis, exotropia, a broad nasal root and bridge and nasal tip with hypoplastic nares and hanging columella. Fundoscopy images for RQMR10 individual VI:2 were examined (E.H.), and showed abnormality compatible with a retinal dystrophy and with Cohen syndrome. The optic nerves were pale and the macular looked thinned and the blood vessels (arterioles) are attenuated (see Additional file 1: Figure S1). Night blindness and progressive vision loss were also observed in the affected individuals.
Table 1

Clinical features of affected members of families ANMR51, RQMR10 and ATM02. n.a. = not assessed

Family

ANMR51

RQMR10

ATM02

Patient ID

IV-5

IV-4

IV-3

V-3

V-1

V-6

V-7

V-8

VI-2

VI-3

VI-1a

V-2

V-1

V-3

V-5

Sex

M

M

F

M

M

M

M

F

M

F

M

F

M

M

M

Age (Yrs)

25

27

17

20

15

>50

>50

>50

17

7

8

20

23

18

15

Neurodevelopment

Developmental delay (motor, cognitive, speech)

Yes

Yes

Yes

Yes

Yes

Yes

Yes

Yes

Yes

Yes

Yes

Yes

Yes

Yes

Yes

Intellectual Disability

Severe

Severe

Severe

Severe

Severe

Severe

Severe

Severe

Severe

Severe

Severe

Severe

Moderate

Moderate

Moderate

Language

Single words

Single words

Single words

Single words

Single words

Single words

Single words

Single words

Single words

No speech

Single words

Single words

Small sentences

Small sentences

Small sentences

Autistic behaviors

Nb

Nb

Nb

Nb

Nb

Nb

Nb

Nb

Nbc

Nb

Nb

Y

Y

Y

Y

Age 1st walked

<2

<2

<2

<2

<2

<2

<2

<2

2

4

2

<2

2.5

2.4

2.0

Hypotonia?

n.a.

n.a.

n.a.

n.a.

n.a.

n.a.

n.a.

n.a.

n.a.

n.a.

n.a.

n.a.

yes

yes

yes

Coordination

poor

poor

poor

poor

poor

poor

poor

poor

poor

poor

poor

poor

poor

poor

poor

Facies

OFC (cm)

49

49

47

n.a.

n.a.

47.5

48

48

n.a.

n.a.

n.a.

n.a.

51

49

49

Thick eyebrows

n.a.

n.a.

n.a.

n.a.

n.a.

No

Yes

n.a.

Yes

No

n.a.

No

Yes

Yes

Yes

Prominent nasal bridge

n.a.

n.a.

n.a.

n.a.

n.a.

No

Yes

n.a.

No

No

n.a.

No

Yes

Yes

Yes

short philtrum

n.a.

n.a.

n.a.

n.a.

n.a.

Yes

Yes

n.a.

No

No

n.a.

No

Yes

Yes

Yes

wide or wave shaped palpebral fissures

n.a.

n.a.

n.a.

n.a.

n.a.

Yes

No

n.a.

No

No

n.a.

No

Yes

Yes

Yes

Prominent incisors

n.a.

n.a.

n.a.

n.a.

n.a.

n.d.

n.d.

n.a.

Yes

Yes

n.a.

Yes

Yes

Yes

Yes

Mild synophrys

n.a.

n.a.

n.a.

n.a.

n.a.

Yes

Yes

n.a.

Yes

Yes

n.a.

Yes

Yes

Yes

Yes

Eyes

Myopia

+

+

+

+

+

Blind

Blind

+

+

+

+

+

+

+

+

Retinopathy

n.a.

n.a.

n.a.

n.a.

n.a.

n.a.

n.a.

n.a.

+ retinal dystrophy

n.a.

n.a.

n.a.

n.a.

n.a.

n.a.

Neutropenia

n.a.

n.a.

n.a.

n.a.

n.a.

n.a.

n.a.

n.a.

-

Pan-cytopenia

n.a.

n.a.

+

n.a.

n.a.

Miscellaneous

  

Total loss of teeth

 

-

Prominent chin

-

-

Ptosis; hyperactivity

  

Pes planus & pes cavus; Truncal obesity

Pes planus hypoactive

Pes planus hyperactive recurrent skin infections

Pes planus; pica

aNow deceased

bthese individuals were aloof, isolated, not sociable; however there were no repetitive behaviors

cHyperactivity also noted

Fig. 2

Photos of affected individuals from the three families

Some autistic-like traits were observed in affected members of RQMR10, who were reported to be aloof, isolated and not sociable; however in the context of their ID it was not possible to assess for autism spectrum disorder. No repetitive behaviors were noted. Based on clinical interview, the level of ID was estimated as moderate to severe in the affected individuals in RQMR10. In regards to social functioning, one female (V-2) was initially considered aloof, and later developed severe depressive disorder. RQMR10 VI-2 is reported as being hyperactive, having no sense of danger, and would frequently wander away from home.

Family ATM02 (Figure 1) is part of a much larger pedigree segregating with ID across many branches, with 18 affected individuals (see Additional file 2: Figure S2). In this study we were able to access just the portion of the family shown in Fig. 1. Three affected brothers presented with Cohen syndrome features (Table 1). In addition, individual V-1 had normal blood and urine tests, except for blood neutrophil at 2.0 × 10−9/L, and his brain MRI indicates non-specific mild to moderate brain atrophy with enlarged ventricles and subarachnoid spaces (see Additional file 3: Figure S3). However, due to lack of information on neutropenia and retinopathy, the diagnosis of COH1 was not made prior to the genetic investigation. They were all initially diagnosed with pervasive developmental disorder not otherwise specified (PDD-NOS). Ophthalmological examination later confirmed that the three affected brothers all present with pigmentary retinopathy.

Head circumference measurements (occipito-frontal circumference) were available for ANMR51 and ATM02, and show measures consistently > 3SD below the mean, suggesting that microcephaly is also a common feature in these families (Table 1).

Mutation Detection in RQMR10, ANMR51 and ATM02 families

RQMR10 Family

We identified only one ~25.430 Mb HBD region shared by the affected individuals, but not unaffected family members, on 8q22, flanked by SNPs rs1994035 and rs3098237 in family RQMR10. Sanger sequence analysis of all PCR amplicons revealed a deletion of 1 bp, T, (Fig. 3) at position chr8:100,732,719 (hg19) on chromosome 8q22.2. This deletion, NM_017890.4:c.6879delT, introduces a framseshift, resulting in premature protein truncation (p.Phe2293Leufs*24). Exon38 of VPS13B was sequenced in family members and in 293 normal unrelated healthy individuals from Pakistan. The homozygous c.6879delT mutation cosegregates with the clinical phenotype correctly within the family, and was not present in 293 healthy Pakistani controls. The variant is present as heterozygote in the Exome Aggregation Consortium Browser (ExAC; Cambridge, MA; http://exac.broadinstitute.org; Oct 2014), which includes over 8,000 exomes from South Asian individuals, but only one out of 122494 alleles total, and one out of 16,592 South Asian alleles.
Fig. 3

a Ideogram of chromosome 8, indicating location of HBD at 8q22. b Genomic organization of VPS13B/COH1, indicating location of the NM_017890.4:c.6879delT mutation (RQMR10 and ATM02) in exon 38, as well as the deleted region (exons 37–39) in ANMR51. c Electropherograms indicating the wild type homozygous (top), heterozygous (middle) and mutant homozygous (bottom) sequence across NM_017890.4:c.6879 (RQMR10 and ATM02). d Electropherogram from mRNA sequencing from an affected individual from ANMR51, indicating the (aberrant) splicing of exons 36 and 40

ANMR51 Family

Analysis of genotypes for family ANMR51 revealed only one ~4.413 Mb HBD region shared by the affected individuals, but not unaffected family members, on 8q22 flanked by SNPs rs11783149 and rs2679749. Four consecutive exons (37–40) failed to amplify by PCR, yet amplified in unaffected family members, suggesting a possible deletion. The deletion fragment was confirmed and mapped by PCR amplification and sequencing of the junction fragment (forward primer-VPS13B-JF-F-5′-TCCTAAAATGTGCCA TTGGTT-3′; reverse primer-VPS13B-JF-R: 5′-TGCTGAAGAATAGATTGCTGAA-3′) (Fig. 3). Sequence analysis of the junction fragment showed the deletion happens between chr8:100,728,884 and 100,779,911 (hg19), spanning a region of 51,027 bp (hg19:g.chr8:100728884_chr8:100779911 del). Assuming utilization of the same acceptor and donor splice sites as for the wild type for the remaining part of the gene, the deletion is predicted to lead to a truncated protein, i.e., p.Gly2177Alafs*16. The deletion was co-segregating as homozygous mutation within the family. The deletion was not present in 293 unrelated healthy Pakistani individuals.

ATM02 family

Homozygosity mapping identified four HBD regions shared in three affected individuals: 1p36.32-p36.31 (rs10799181 to rs9308476, 1.87 Mb); 1p13.2-q21.3 (rs11804649 to rs6655975, 37.19 Mb), 8q21.11-q24.21 (rs7824974 to rs4407842, 52.49 Mb), and 16p13.3-p12.1 (rs2541696 to rs8047869, 24.27 Mb). Exome sequencing and analysis of two affected individuals (V1 and V5) led to the identification of the same 1 bp homozygous deletion at c.6879 T found in RQMR10. In the meantime, we have excluded any additional loss of function variant(s) (both homozygote and heterozygote) in genes known to carry pathogenic mutations for ID and ASD (SFARI gene list), in the two affected individuals, separately and jointly. The deletion at c.6879 T was confirmed by Sanger sequencing as homozygous in three affected individuals, but absent in unaffected family members and 293 controls, as mentioned for RQMR10.

Haplotype sharing between RQMR10 and ATM02

Interestingly, RQMR10 and ATM02 were found to carry the same c.6879delT mutation independently. Subsequent joint analysis of shared HBD regions showed that among 592 common SNPs shared between Affymetrix and Illumina SNP chips within the overlapping HBD region on chr 8q22.2 in both families, 208 SNPs from a 2.14 Mb region showed the same haplotype (see in Fig. 4, also Additional file 4: Table S1), indicating a founder effect of this mutation within the population of Pakistan. Neither family is aware of any common ancestry, and has lived in regions separated by a distance of 534 km, although both are from the Baloch ethnic group. Based on the size of shared haplotype, the two families may share common ancestors approximately 50 generations ago (mean distance between the nearest flanking crossovers of the disease locus is 2/N morgans, where N is the number of meioses) [26].
Fig. 4

Homozygosity by descent of the Pakistani families and shared haplotype region for RQMR10 and ATM02. The homozygous-by-descent regions segregating with the phenotype within families RQMR10, ATM02 and ANMR51 are indicated with a dark blue, a green and a purple bar, respectively. The red bar indicates the region of shared haplotype for families RQMR10 and ATM02 and spans the entire VPS13B

Discussion

We present a study of three large consanguineous Pakistani families with autosomal recessive ID and in some cases autistic behaviors, segregating with deletion mutations within the VPS13B gene. This large gene, spanning a region of 864 Kb with 62 exons, maps to chromosome 8q22.2 [7, 8]. The gene encodes a protein of 4,022 amino acids (plus a number of shorter isoforms) with several short regions showing homology with yeast vacuolar protein sorting-associated protein 13 (Vps13b), suggesting a possible role in vesicle-mediated protein sorting (VPS) and intracellular protein transportation [13, 14, 27, 28]. Functional studies show evidence for VPS13B co-localizes strongly with the cis-Golgi matrix protein GM130 [29]. The study also showed that VPS13B/COH1 is required for the maintenance of Golgi stacks as well as for reassembly of Golgi cisternae from ministacks, and can regulate the formation of membrane tubules from the Golgi. These functions suggest similarity with Golgins (Golgi autoantigens), except that VPS13B/COH1 does not contain coiled-coil domains. Absence or depletion of VPS13B/COH1 could lead to disruption in the Golgi complex integrity and functionality [29]. More recent work further demonstrated that VPS13B mutations, responsible for COH1, are associated with strong protein glycosylation defects, highlighting an important role for VPS13B in Golgi glycosylation and morphology, as well as in lysosomal–endosomal pathway maintenance [30].

The first mutation we report, a 51 Kb homozygous deletion encompassing four coding exons, assuming correct splicing of the flanking exons 36 to 41 during RNA maturation and processing, 793 bp of coding sequence would be deleted, leading to a shift in reading frame, with the substitution glycine to alanine at amino acid 2177, followed by a premature stop codon after a further 14 residues, leading to a 2191 amino acid residue protein (compared to 4022 for wild type, NM_017890). This premature stop codon would either lead to nonsense-mediated mRNA decay (NMD) or production of a truncated protein lacking a number of functional domains. This potential shortened protein (p.Gly2177Alafs*16) would lack domains DUF1162 and ATG_C (using SMART). The genomic region of the deletion contains no known segmental duplications, large repeats or other sequences that may be likely to predispose towards genomic instability (UCSC Genome Browser, hg19), and no CNVs have been identified at this region, according to the Database of Genomic Variants.

The 1 bp deletion at c.6879 T [NM_017890.4 (VPS13B_v001):c.6879delT], leads to a shift in frame and a premature stop codon (PTC) after amino acid 2315 (p.Phe2293Leufs*24). This premature stop codon would also either lead to NMD or production of nonfunctional protein. The truncated protein would lack many of the main predicted functional domains of the VPS13B protein, including Pfam predicted domains DUF1162 and ATG_C.

Since the two mutations would have a very similar functional effect, with truncation at similar points along the protein (2191 residues for RQMR10 and ATM02 mutation; 2315 residues for ANMR51 mutation), one might expect the resulting phenotype to be very similar. However, the affected individuals from these three families showed remarkable phenotypic variations, which could be due to other factors, including the degree of NMD that may occur for the two mutations, which in itself may depend on the nucleotide context of the premature stop codon and the assembly of the NMD complex [31]. Interestingly, most affected members of the two families with the c.6879delT mutation were either evaluated as having autism spectrum disorder (ASD) or of possessing autistic-like traits. Mutations in VPS13B have recently been found in additional cases of ASD or ID individuals with autistic features. Yu et al. in a study of 163 consanguineous and/or multiplex families with ASD through exome sequencing, found a homozygous frameshift mutation (Ser3945Glnfs*22) in a proband with ASD and mild dysmorphic features, as well as a rare missense mutation (p.S824A) in a proband with autism and mild facial dysmorphism and joint laxity [32]. In the same study, screening VPS13B gene in 612 families from the Simons Simplex Collection, one affected male child with ASD was found to be compound heterozygous for two different mutations in VPS13B (p.W963*/p.G2704R), and a second, unrelated male child affected with autism was compound heterozygous for two rare point mutations in VPS13B, p.S3303R and p.A3691T, both altering highly conserved residues of VPS13B [32]. Re-analyzing exome sequencing data of 488 ASD cases and 372 controls from the ARRA Autism Sequencing Collaboration (AASC) project, Ionita-Laza I et al. found two loss-of-function variants (one nonsense: p.Ser3383* and one splice site: c:2650 + 2 T > G), and a rare, homozygous probably-damaging variant (p.Arg3198Trp) in VPS13B that could contribute to autism risk [33]. Although only a single previous report of VPS13B mutations and COH1 in individuals of Pakistani origin has been published [34], identification of three unrelated large families from Pakistan carrying homozygous mutations in VPS13B indicates that this is an important gene, not only for COH1, but also for undiagnosed or misdiagnosed ID and ASD patients, in the Pakistani population. We believe that Cohen syndrome individuals with this “Baloch” variant may be enriched for ASD or autistic-like behaviors, which contrasts with the “friendly disposition” frequently associated with COH1. Early genetic screen for the extended families or clans, and local communities would be extremely useful in genetic counseling, and clinical diagnosis.

Conclusion

We suspect that the c.6879delT mutation may be a common cause of COH1 and similar phenotypes among the Baloch population. Additionally, most of the individuals with the c.6879delT mutation in these two families also present with autistic like traits, and suggests that this variant may lead to a distinct autistic-like COH1 subgroup.

Bioinformatic resources

  1. 1.

    Primer3 software (version 0.4.0): http://frodo.wi.mit.edu/primer3/

     
  2. 2.

    UCSC Genome Browser: http://genome.ucsc.edu/

     
  3. 3.

    Exome Aggregation Consortium browser (ExAC; Cambridge, MA): http://exac.broadinstitute.org

     
  4. 4.
     
  5. 5.
     
  6. 6.
     
  7. 7.
     

Notes

Declarations

Acknowledgements

This research was supported by grants to JBV (#MOP-102758) and LX (#MOP-259333) from the Canadian Institutes of Health Research. We thank all the patients and their families, without their participation this project would not have been possible. We also thank Dr. Elise Héon of the Department of Ophthalmology and Vision Sciences, The Hospital for Sick Children, University of Toronto, Toronto, Canada, for assistance interpreting fundoscopy images.

Authors’ Affiliations

(1)
Molecular Neuropsychiatry & Development Laboratory, Campbell Family Mental Health Research Institute, Centre for Addiction and Mental Health
(2)
CHUM Research Center - Notre Dame Hospital
(3)
Department of Biochemistry, Quaid-I-Azam University, and Pakistan Medical Research Council
(4)
The Lahore Institute for Research and Development
(5)
University of Education, Township Campus
(6)
Department of Biotechnology, University of Azad Jammu and Kashmir
(7)
Dept: zoology, University of Education, Lahore
(8)
Ma Ayesha Memorial Centre
(9)
Institute of Human Genetics, Medical University of Graz
(10)
Hôpital Rivière-des-Prairies
(11)
Département de Psychiatrie, Université de Montréal
(12)
Division of Developmental Disabilities, Department of Psychiatry, Queen’s University
(13)
Autism Institute
(14)
Research Centre, Montreal Mental Health University Institute
(15)
The Fred A. Litwin Family Centre in Genetic Medicine, University Health Network and Mount Sinai Hospital
(16)
The Centre for Addiction and Mental Health
(17)
Department of Laboratory Medicine and Pathobiology, University of Toronto
(18)
Sir Cowasji Jehangir Institute of Psychiatry
(19)
Department of Psychiatry, University of Toronto
(20)
Currently at: Department of Physiology and Experimental Medicine (PEM), Hospital for Sick Children
(21)
Currently at: Montreal Neurological Institute, McGill University

References

  1. Maulik PK, Mascarenhas MN, Mathers CD, Dua T, Saxena S. Prevalence of intellectual disability: a meta-analysis of population-based studies. Res Dev Disabil. 2011;32(2):419–36.PubMedGoogle Scholar
  2. Rauch A, Hoyer J, Guth S, et al. Diagnostic yield of various genetic approaches in patients with unexplained developmental delay or mental retardation. Am J Med Genet. 2006;140A(19):2063–74.Google Scholar
  3. Cohen Jr MM, Hall BD, Smith DW, Graham CB, Lampert KJ. A new syndrome with hypotonia, obesity, mental deficiency, and facial, oral, ocular and limb anomalies. J Pediat. 1973;83(2):280–4.PubMedGoogle Scholar
  4. Kivitie-Kallio S, Norio R. Cohen Syndrome: Essential features, natural history, and heterogeneity. Am J Med Genet. 2001;102(2):125–35.PubMedGoogle Scholar
  5. Chandler KE, Kidd A, Al-Gazali L, et al. Diagnostic criteria, clinical characteristics, and natural history of Cohen syndrome. J Med Genet. 2003;40(4):233–41.PubMedPubMed CentralGoogle Scholar
  6. Douzgou S, Peterson MB. Clinical variability of genetic isolates of Cohen syndrome. Clin Genet. 2011;79(6):501–6.PubMedGoogle Scholar
  7. Kolehmainen J, Black GC, Saarinen A, et al. Cohen syndrome is caused by mutations in a novel gene, COH1, encoding a transmembrane protein with a presumed role in vesicle mediated sorting and intracellular protein transport. Am J Hum Genet. 2003;72(6):1359–69.PubMedPubMed CentralGoogle Scholar
  8. Kolehmainen J, Wilkinson R, Lehesjoki AE, et al. Delineation of Cohen syndrome following a large-scale genotype-phenotype screen. Am J Hum Genet. 2004;75(1):122–7.PubMedPubMed CentralGoogle Scholar
  9. Hennies HC, Rauch A, Seifert W, et al. Allelic heterogeneity in the COH1 gene explains clinical variability in Cohen syndrome. Am J Hum Genet. 2004;75(1):138–45.PubMedPubMed CentralGoogle Scholar
  10. Seifert W, Holder-Espinasse M, Spranger S, et al. Mutational spectrum of COH1 and clinical heterogeneity in Cohen syndrome. J Med Genet. 2006;43(5):e22.PubMedPubMed CentralGoogle Scholar
  11. Katzaki E, Pescucci C, Uliana V, et al. Clinical and molecular characterization of Italian patients affected by Cohen syndrome. J Hum Genet. 2007;52(12):1011–7.PubMedGoogle Scholar
  12. Bugiani M, Gyftodimou Y, Tsimpouka P, et al. Cohen syndrome resulting from a novel large intragenic COH1 deletion segregating in an isolated Greek island population. Am J MedGenet A. 2008;146A(17):2221–6.Google Scholar
  13. Balikova I, Lehesjoki AE, de Ravel TJ, et al. Deletions in the VPS13B (COH1) gene as a cause of Cohen syndrome. Hum Mutat. 2009;30(9):E845–54.Google Scholar
  14. El Chehadeh S, Aral B, Gigot N, et al. Search for the best indicators for the presence of a VPS13B gene mutation and confirmation of diagnostic criteria in a series of 34 patients genotyped for suspected Cohen syndrome. J Med Genet. 2010;47(8):549–53.PubMedGoogle Scholar
  15. Rivera-Brugu’es N, Albrecht B, Wieczorek D, et al. Cohen syndrome diagnosis using whole genome arrays. J Med Genet. 2011;48(2):136–40.Google Scholar
  16. Parri V, Katzaki E, Uliana V, et al. High frequency of COH1intragenic deletions and duplications detected by MLPA inpatients with Cohen syndrome. Eur J Hum Genet. 2010;18(10):1133–40.PubMedPubMed CentralGoogle Scholar
  17. El Chehadeh-Djebbar S, Blair E, Holder-Espinasse M, et al. Changing facial phenotype in Cohen syndrome: towards clues for an earlier diagnosis. Eur J Hum Genet. 2013;21(7):736–42.PubMedGoogle Scholar
  18. Gueneau L, Duplomb L, Sarda P, et al. Congenital neutropenia with retinopathy, a new phenotype without intellectual deficiency or obesity secondary to VPS13B mutations. Am J Med Genet A. 2014;164A(2):522–7.PubMedGoogle Scholar
  19. Schopler E, Reichler RJ, DeVellis RF, Daly K. Toward objective classification of childhood autism: Childhood Autism Rating Scale (CARS). J Autism Dev Disord. 1980;10(1):91–103.PubMedGoogle Scholar
  20. Zhao X, Li C, Paez JG, Chin K, Janne PA, Chen TH, et al. An integrated view of copy number and allelic alterations in the cancer genome using single nucleotide polymorphism arrays. Cancer Res. 2004;64(9):3060–71.PubMedGoogle Scholar
  21. Zhao X, Weir BA, LaFramboise T, Lin M, Beroukhim R, Garraway L, et al. Homozygous deletions and chromosome amplifications in human lung carcinomas revealed by single nucleotide polymorphism array analysis. Cancer Res. 2005;65(13):5561–70.PubMedGoogle Scholar
  22. Seelow D, Schuelke M, Hildebrandt F, Nürnberg P. HomozygosityMapper--an interactive approach to homozygosity mapping. Nucleic Acids Res. 2009;37(Web Server issue):W593–599.PubMedPubMed CentralGoogle Scholar
  23. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler Transform. Bioinformatics. 2009;25(14):1754–60.PubMedPubMed CentralGoogle Scholar
  24. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.PubMedPubMed CentralGoogle Scholar
  25. Wang K, Li M, Hakonarson H. ANNOVAR: Functional annotation of genetic variants from next-generation sequencing data. Nucleic Acids Res. 2010;38(16):e164.PubMedPubMed CentralGoogle Scholar
  26. Boehnke M. Limits of resolution of genetic linkage studies: implications for the positional cloning of human disease genes. Am J Hum Genet. 1994;55(2):379–90.PubMedPubMed CentralGoogle Scholar
  27. Redding K, Brickner JH, Marschall LG, et al. Allele-specific suppression of a defective trans-Golgi network (TGN) localization signal in Kex2p identifies three genes involved in localization of TGN transmembrane proteins. Mol Cell Biol. 1996;16(11):6208–17.PubMedPubMed CentralGoogle Scholar
  28. Seifert W, Holder-Espinasse M, Kühnisch J, et al. Expanded mutational spectrum in Cohen syndrome, tissue expression, and transcript variants of COH1. Hum Mutat. 2009;30(2):E404–20.PubMedGoogle Scholar
  29. Seifert W, Kühnisch J, Maritzen T, et al. Cohen syndrome-associated protein, COH1, is a novel, giant Golgi matrix protein required for Golgi integrity. J Biol Chem. 2011;286(43):37665–75.PubMedPubMed CentralGoogle Scholar
  30. Duplomb L, Duvet S, Picot D, et al. Cohen syndrome is associated with major glycosylation defects. Hum Mol Genet. 2014;23(9):2391–9.PubMedGoogle Scholar
  31. Kervestin S, Jacobson A. NMD: a multifaceted response to premature translational termination. Nat Rev Mol Cell Biol. 2012;13(11):700–12.PubMedPubMed CentralGoogle Scholar
  32. Yu TW, Chahrour MH, Coulter ME, et al. Using whole exome sequencing to identify inherited causes of autism. Neuron. 2013;77(2):259–73.PubMedPubMed CentralGoogle Scholar
  33. Ionita-Laza I, Capanu M, De Rubeis S, McCallum K, Buxbaum JD. Identification of Rare Causal Variants in Sequence-Based Studies: Methods and Applications to VPS13B, a Gene Involved in Cohen Syndrome and Autism. Zeggini E, ed. PLoS Genetics. 2014;10(12):e1004729. doi:10.1371/journal.pgen.1004729.PubMedPubMed CentralGoogle Scholar
  34. Waite A, Somer M, O’Driscoll M, Millen K, Manson FD, Chandler KE. Cerebellar hypoplasia and Cohen syndrome: a confirmed association. Am J Med Genet A. 2010;152A(9):2390–3.PubMedGoogle Scholar

Copyright

© Rafiq et al. 2015

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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Advertisement