Expanding the genetic and phenotypic spectrum of congenital myasthenic syndrome: new homozygous VAMP1 splicing variants in 2 novel individuals.

We report the cases of two Spanish pediatric patients with hypotonia, muscle weakness and feeding difficulties at birth. Whole-exome sequencing (WES) uncovered two new homozygous VAMP1 (Vesicle Associated Membrane Protein 1) splicing variants, NM_014231.5:c.129+5 G > A in the boy patient (P1) and c.341-24_341-16delinsAGAAAA in the girl patient (P2). This gene encodes the vesicle-associated membrane protein 1 (VAMP1) that is a component of a protein complex involved in the fusion of synaptic vesicles with the presynaptic membrane. VAMP1 has a highly variable C-terminus generated by alternative splicing that gives rise to three main isoforms (A, B and D), being VAMP1A the only isoform expressed in the nervous system. In order to assess the pathogenicity of these variants, expression experiments of RNA for VAMP1 were carried out. The c.129+5 G > A and c.341-24_341-16delinsAGAAAA variants induced aberrant splicing events resulting in the deletion of exon 2 (r.5_131del; p.Ser2TrpfsTer7) in the three isoforms in the first case, and the retention of the last 14 nucleotides of the 3' of intron 4 (r.340_341ins341-14_341-1; p.Ile114AsnfsTer77) in the VAMP1A isoform in the second case. Pathogenic VAMP1 variants have been associated with autosomal dominant spastic ataxia 1 (SPAX1) and with autosomal recessive presynaptic congenital myasthenic syndrome (CMS). Our patients share the clinical manifestations of CMS patients with two important differences: they do not show the typical electrophysiological pattern that suggests pathology of pre-synaptic neuromuscular junction, and their muscular biopsies present hypertrophic fibers type 1. In conclusion, our data expand both genetic and phenotypic spectrum associated with VAMP1 variants.

Clinical description and genetic analysis of a novel familial skeletal dysplasia characterized by high bone mass and lucent bone lesions.

High bone mass (HBM) disorders are a clinically and genetically heterogeneous subgroup of rare skeletal dysplasias. Here we present a case of a previously unreported familial skeletal dysplasia characterized by HBM and lucent bone lesions that we aimed to clinically characterize and genetically investigate. For phenotyping, we reviewed past clinical records and imaging tests, and performed physical examination (PE), bone densitometry, and mineral panels in affected individuals, including a male proband, his son and daughter, in addition to unaffected controls, including the proband's wife and brother. Affected individuals also underwent impact microindentation (IMI). In an effort to elucidate the disorder's molecular etiology, whole exome sequencing (WES) was performed in all individuals to filter for rare variants present only in affected ones. The cases displayed a unique skeletal phenotype with a mix of sclerotic features and lucent bone lesions, and high IMI values. Bone mineral density was very elevated in the proband and his daughter. The proband's daughter also exhibited idiopathic scoliosis (IS), in addition to mild thrombocytopenia and mild structural thyroid abnormalities, which were the only extra-skeletal abnormalities identified. WES analysis yielded 5 rare putative pathogenic variants in affected members in genes that are associated with bone metabolism including: SEM4AD, TBX18, PTCH1, PTK7, and ADGRE5. The PTK7 variant appeared as possibly implicated in the development of IS while the TBX18 and SEMA4D variants stood out as the strongest candidates for the lucent bone lesions and HBM, respectively, given their high predicted pathogenicity and putative role in bone biology. Variant functionality should be addressed in the future to assess their implication in skeletal metabolism as it is the first time that mutations in TBX18 and SEMA4D have been associated to bone developmental lesions and mineral metabolism in a clinical setting.

Five new cases of syndromic intellectual disability due to KAT6A mutations: widening the molecular and clinical spectrum.

BACKGROUND: Pathogenic variants of the lysine acetyltransferase 6A or KAT6A gene are associated with a newly identified neurodevelopmental disorder characterized mainly by intellectual disability of variable severity and speech delay, hypotonia, and heart and eye malformations. Although loss of function (LoF) mutations were initially reported as causing this disorder, missense mutations, to date always involving serine residues, have recently been associated with a form of the disorder without cardiac involvement. RESULTS: In this study we present five new patients, four with truncating mutations and one with a missense change and the only one not presenting with cardiac anomalies. The missense change [p.(Gly359Ser)], also predicted to affect splicing by in silico tools, was functionally tested in the patient's lymphocyte RNA revealing a splicing effect for this allele that would lead to a frameshift and premature truncation. CONCLUSIONS: An extensive revision of the clinical features of these five patients revealed high concordance with the 80 cases previously reported, including developmental delay with speech delay, feeding difficulties, hypotonia, a high bulbous nose, and recurrent infections. Other features present in some of these five patients, such as cryptorchidism in males, syndactyly, and trigonocephaly, expand the clinical spectrum of this syndrome.

Genetic factors contributing to autism spectrum disorder in Williams-Beuren syndrome.

BACKGROUND: The hallmark of the neurobehavioural phenotype of Williams-Beuren syndrome (WBS) is increased sociability and relatively preserved language skills, often described as opposite to autism spectrum disorders (ASD). However, the prevalence of ASD in WBS is 6-10 times higher than in the general population. We have investigated the genetic factors that could contribute to the ASD phenotype in individuals with WBS. METHODS: We studied four males and four females with WBS and a confirmed diagnosis of ASD by the Autism Diagnostic Interview-Revised. We performed a detailed molecular characterisation of the deletion and searched for genomic variants using exome sequencing. RESULTS: A de novo deletion of 1.55 Mb (6 cases) or 1.83 Mb (2 cases) at 7q11.23 was detected, being in 7/8 patients of paternal origin. No common breakpoint, deletion mechanism or size was found. Two cases were hemizygous for the rare T allele at rs12539160 in MLXIPL, previously associated with ASD. Inherited rare variants in ASD-related or functionally constrained genes and a de novo nonsense mutation in the UBR5 gene were identified in six cases, with higher burden in females compared with males (p=0.016). CONCLUSIONS: The increased susceptibility to ASD in patients with WBS might be due to additive effects of the common WBS deletion, inherited and de novo rare sequence variants in ASD-related genes elsewhere in the genome, with higher burden of deleterious mutations required for females, and possible hypomorphic variants in the hemizygous allele or cis-acting mechanisms on imprinting.

bigSCale: an analytical framework for big-scale single-cell data.

Single-cell RNA sequencing (scRNA-seq) has significantly deepened our insights into complex tissues, with the latest techniques capable of processing tens of thousands of cells simultaneously. Analyzing increasing numbers of cells, however, generates extremely large data sets, extending processing time and challenging computing resources. Current scRNA-seq analysis tools are not designed to interrogate large data sets and often lack sensitivity to identify marker genes. With bigSCale, we provide a scalable analytical framework to analyze millions of cells, which addresses the challenges associated with large data sets. To handle the noise and sparsity of scRNA-seq data, bigSCale uses large sample sizes to estimate an accurate numerical model of noise. The framework further includes modules for differential expression analysis, cell clustering, and marker identification. A directed convolution strategy allows processing of extremely large data sets, while preserving transcript information from individual cells. We evaluated the performance of bigSCale using both a biological model of aberrant gene expression in patient-derived neuronal progenitor cells and simulated data sets, which underlines the speed and accuracy in differential expression analysis. To test its applicability for large data sets, we applied bigSCale to assess 1.3 million cells from the mouse developing forebrain. Its directed down-sampling strategy accumulates information from single cells into index cell transcriptomes, thereby defining cellular clusters with improved resolution. Accordingly, index cell clusters identified rare populations, such as reelin (Reln)-positive Cajal-Retzius neurons, for which we report previously unrecognized heterogeneity associated with distinct differentiation stages, spatial organization, and cellular function. Together, bigSCale presents a solution to address future challenges of large single-cell data sets.

Detectable clonal mosaicism in blood as a biomarker of cancer risk in Fanconi anemia.

Detectable clonal mosaicism for large chromosomal events has been associated with aging and an increased risk of hematological and some solid cancers. We hypothesized that genetic cancer predisposition disorders, such as Fanconi anemia (FA), could manifest a high rate of chromosomal mosaic events (CMEs) in peripheral blood, which could be used as early biomarkers of cancer risk. We studied the prevalence of CMEs by single-nucleotide polymorphism (SNP) array in 130 FA patients' blood DNA and their impact on cancer risk. We detected 51 CMEs (4.4-159 Mb in size) in 16 out of 130 patients (12.3%), of which 9 had multiple CMEs. The most frequent events were gains at 3q (n = 6) and 1q (n = 5), both previously associated with leukemia, as well as rearrangements with breakpoint clustering within the major histocompatibility complex locus (P = 7.3 × 10-9). Compared with 15 743 age-matched population controls, FA patients had a 126 to 140 times higher risk of detectable CMEs in blood (P < 2.2 × 10-16). Prevalent and incident hematologic and solid cancers were more common in CME carriers (odds ratio [OR] = 11.6, 95% confidence interval [CI] = 3.4-39.3, P = 2.8 × 10-5), leading to poorer prognosis. The age-adjusted hazard risk (HR) of having cancer was almost 5 times higher in FA individuals with CMEs than in those without CMEs. Regarding survival, the HR of dying was 4 times higher in FA individuals having CMEs (HR = 4.0, 95% CI = 2.0-7.9, P = 5.7 × 10-5). Therefore, our data suggest that molecular karyotyping with SNP arrays in easy-to-obtain blood samples could be used for better monitoring of bone marrow clonal events, cancer risk, and overall survival of FA patients.

Targeting tryptophan and tyrosine metabolism by liquid chromatography tandem mass spectrometry.

An imbalance in tryptophan (Trp) and tyrosine (Tyr) metabolites is associated with neurological and inflammatory disorders. The accurate and precise measurement of these compounds in biological specimens is a powerful tool to understand the biochemical state in several diseases. In this study, a rapid, accurate and sensitive method based on liquid chromatography-tandem mass spectrometry (LC-MS/MS) for the targeted analysis of the metabolism of Trp and Tyr has been developed and validated. The method allows for the adequate quantification of Trp, Tyr and, eight Trp metabolites, three Tyr metabolites, together with four competitive large neutral amino acids. Serotonin, 5-hydroxyindoleacetic acid, kynurenine, kynurenic acid, dopamine, and homovanilic acid were among the targeted compounds. Sample preparation, chromatographic separation and mass spectrometric detection were optimized in human urine, human plasma and mice prefrontal cortex extracts. The method was shown to be linear (r>0.98) in the range of endogenous concentrations for all studied metabolites. In general, the limits of detection were suitable for the detection of the endogenous levels. Intra- and inter-assay precisions below 25% and accuracies ranging from 80 to 120% were found for most of the analytes. The use of labeled internal standards corrected the moderate matrix effect observed for some compounds. The applicability of the method was confirmed by analyzing urine samples collected from 13 healthy volunteers and comparing the results with previously established normal ranges. In addition, urine samples from two patients and a heterozygous carrier of a family with disturbed monoamine metabolism due to a loss of function mutation in the MAOA gene (X-linked) were analyzed and compared with samples from controls. All data together show the potential of the developed approach for targeted metabolomic studies.

Integrated analysis of whole-exome sequencing and transcriptome profiling in males with autism spectrum disorders.

BACKGROUND: Autism spectrum disorders (ASD) are a group of neurodevelopmental disorders with high heritability. Recent findings support a highly heterogeneous and complex genetic etiology including rare de novo and inherited mutations or chromosomal rearrangements as well as double or multiple hits. METHODS: We performed whole-exome sequencing (WES) and blood cell transcriptome by RNAseq in a subset of male patients with idiopathic ASD (n = 36) in order to identify causative genes, transcriptomic alterations, and susceptibility variants. RESULTS: We detected likely monogenic causes in seven cases: five de novo (SCN2A, MED13L, KCNV1, CUL3, and PTEN) and two inherited X-linked variants (MAOA and CDKL5). Transcriptomic analyses allowed the identification of intronic causative mutations missed by the usual filtering of WES and revealed functional consequences of some rare mutations. These included aberrant transcripts (PTEN, POLR3C), deregulated expression in 1.7% of mutated genes (that is, SEMA6B, MECP2, ANK3, CREBBP), allele-specific expression (FUS, MTOR, TAF1C), and non-sense-mediated decay (RIT1, ALG9). The analysis of rare inherited variants showed enrichment in relevant pathways such as the PI3K-Akt signaling and the axon guidance. CONCLUSIONS: Integrative analysis of WES and blood RNAseq data has proven to be an efficient strategy to identify likely monogenic forms of ASD (19% in our cohort), as well as additional rare inherited mutations that can contribute to ASD risk in a multifactorial manner. Blood transcriptomic data, besides validating 88% of expressed variants, allowed the identification of missed intronic mutations and revealed functional correlations of genetic variants, including changes in splicing, expression levels, and allelic expression.

Common and rare variants of microRNA genes in autism spectrum disorders.

OBJECTIVES: MicroRNAs (miRNAs) are post-transcriptional regulators that have been shown to be involved in disease susceptibility. Here we explore the possible contribution of common and rare variants in miRNA genes in autism spectrum disorders (ASD). METHODS: A total of 350 tag SNPs from 163 miRNA genes were genotyped in 636 ASD cases and 673 controls. A replication study was performed in a sample of 449 ASD cases and 415 controls. Additionally, rare variants in 701 miRNA genes of 41 ASD patients were examined using whole-exome sequencing. RESULTS: The most significant association in the discovery sample was obtained for the miR-133b/miR-206 cluster (rs16882131, P = 0.00037). The replication study did not reach significance. However, the pooled analysis (1,085 cases and 1,088 controls) showed association with two miRNA clusters: miR-133b/miR-206 (rs16882131, permP = 0.037) and miR-17/miR-18a/miR-19a/miR-20a/miR-19b-1/miR92a-1 (rs6492538, permP = 0.019). Both miR-133b and miR-206 regulate the MET gene, previously associated with ASD. Rare variant analysis identified mutations in several miRNA genes, among them miR-541, a brain-specific miRNA that regulates SYN1, found mutated in ASD. CONCLUSIONS: Although our results do not establish a clear role for miRNAs in ASD, we pinpointed a few candidate genes. Further exome and GWAS studies are warranted to get more insight into their potential contribution to the disorder.

Metabolic abnormalities in Williams-Beuren syndrome.

BACKGROUND: Williams-Beuren syndrome (WBS, OMIM-194050) is a neurodevelopmental disorder with multisystemic manifestations caused by a 1.55-1.83 Mb deletion at 7q11.23 including 26-28 genes. Reported endocrine and metabolic abnormalities include transient hypercalcaemia of infancy, subclinical hypothyroidism in ∼ 30% of children and impaired glucose tolerance in ∼ 75% of adult individuals. The purpose of this study was to further study metabolic alterations in patients with WBS, as well as in several mouse models, to establish potential candidate genes. METHODS: We analysed several metabolic parameters in a cohort of 154 individuals with WBS (data available from 69 to 151 cases per parameter), as well as in several mouse models with complete and partial deletions of the orthologous WBS locus, and searched for causative genes and potential modifiers. RESULTS: Triglyceride plasma levels were significantly decreased in individuals with WBS while cholesterol levels were slightly decreased compared with controls. Hyperbilirubinemia, mostly unconjugated, was found in 18.3% of WBS cases and correlated with subclinical hypothyroidism and hypotriglyceridemia, suggesting common pathogenic mechanisms. Haploinsufficiency at MLXIPL and increased penetrance for hypomorphic alleles at the UGT1A1 gene promoter might underlie the lipid and bilirubin alterations. Other disturbances included increased protein and iron levels, as well as the known subclinical hypothyroidism and glucose intolerance. CONCLUSIONS: Our results show that several unreported biochemical alterations, related to haploinsufficiency for specific genes at 7q11.23, are relatively common in WBS. The early diagnosis, follow-up and management of these metabolic disturbances could prevent long-term complications in this disorder.

Neurotransmitter systems and neurotrophic factors in autism: association study of 37 genes suggests involvement of DDC.

OBJECTIVES: Neurotransmitter systems and neurotrophic factors can be considered strong candidates for autism spectrum disorder (ASD). The serotoninergic and dopaminergic systems are involved in neurotransmission, brain maturation and cortical organization, while neurotrophic factors (NTFs) participate in neurodevelopment, neuronal survival and synapses formation. We aimed to test the contribution of these candidate pathways to autism through a case-control association study of genes selected both for their role in central nervous system functions and for pathophysiological evidences. METHODS: The study sample consisted of 326 unrelated autistic patients and 350 gender-matched controls from Spain. We genotyped 369 tagSNPs to perform a case-control association study of 37 candidate genes. RESULTS: A significant association was obtained between the DDC gene and autism in the single-marker analysis (rs6592961, P = 0.00047). Haplotype-based analysis pinpointed a four-marker combination in this gene associated with the disorder (rs2329340C-rs2044859T-rs6592961A-rs11761683T, P = 4.988e-05). No significant results were obtained for the remaining genes after applying multiple testing corrections. However, the rs167771 marker in DRD3, associated with ASD in a previous study, displayed a nominal association in our analysis (P = 0.023). CONCLUSIONS: Our data suggest that common allelic variants in the DDC gene may be involved in autism susceptibility.

Clinical utility of chromosomal microarray analysis in invasive prenatal diagnosis.

Novel methodologies for detection of chromosomal abnormalities have been made available in the recent years but their clinical utility in prenatal settings is still unknown. We have conducted a comparative study of currently available methodologies for detection of chromosomal abnormalities after invasive prenatal sampling.A multicentric collection of a 1-year series of fetal samples with indication for prenatal invasive sampling was simultaneously evaluated using three screening methodologies: (1) karyotype and quantitative fluorescent polymerase chain reaction (QF-PCR), (2) two panels of multiplex ligation-dependent probe amplification (MLPA), and (3) chromosomal microarray-based analysis (CMA) with a targeted BAC microarray. A total of 900 pregnant women provided informed consent to participate (94% acceptance rate). Technical performance was excellent for karyotype, QF-PCR, and CMA (~1% failure rate), but relatively poor for MLPA (10% failure). Mean turn-around time (TAT) was 7 days for CMA or MLPA, 25 for karyotype, and two for QF-PCR, with similar combined costs for the different approaches. A total of 57 clinically significant chromosomal aberrations were found (6.3%), with CMA yielding the highest detection rate (32% above other methods). The identification of variants of uncertain clinical significance by CMA (17, 1.9%) tripled that of karyotype and MLPA, but most alterations could be classified as likely benign after proving they all were inherited. High acceptability, significantly higher detection rate and lower TAT, could justify the higher cost of CMA and favor targeted CMA as the best method for detection of chromosomal abnormalities in at-risk pregnancies after invasive prenatal sampling.

Reassessing the role of mitochondrial DNA mutations in autism spectrum disorder.

BACKGROUND: There is increasing evidence that impairment of mitochondrial energy metabolism plays an important role in the pathophysiology of autism spectrum disorders (ASD; OMIM number: 209850). A significant proportion of ASD cases display biochemical alterations suggestive of mitochondrial dysfunction and several studies have reported that mutations in the mitochondrial DNA (mtDNA) molecule could be involved in the disease phenotype. METHODS: We analysed a cohort of 148 patients with idiopathic ASD for a number of mutations proposed in the literature as pathogenic in ASD. We also carried out a case control association study for the most common European haplogroups (hgs) and their diagnostic single nucleotide polymorphisms (SNPs) by comparing cases with 753 healthy and ethnically matched controls. RESULTS: We did not find statistical support for an association between mtDNA mutations or polymorphisms and ASD. CONCLUSIONS: Our results are compatible with the idea that mtDNA mutations are not a relevant cause of ASD and the frequent observation of concomitant mitochondrial dysfunction and ASD could be due to nuclear factors influencing mitochondrion functions or to a more complex interplay between the nucleus and the mitochondrion/mtDNA.

Intermediate FMR1 alleles and cognitive and/or behavioural phenotypes.

During the last few years, several studies have reported an excess of intermediate FMR1 alleles in patients with cognitive and/or behavioural phenotypes. Here, we report the frequency of intermediate alleles (IAs) in three pathologies, intellectual disabilities (IDs), attention-deficit/hyperactivity disorder and autism, from different Spanish regions. We found 142 IAs among 9015 patients with ID (1.6%), 4 among the 415 ADHD patients (0.96%) and 4 among the 300 autistic patients (1.3%), similar to the frequency reported in our control population. No evidence was found of an excess of IA at the FRAXA locus in any of the study populations, although geographical variability was detected. Moreover, the analysis of 100 transmissions of IAs showed that 95% of these alleles were stable. Only 3% expanded within the same range and 2% expanded to a full mutation in two generations. No evidence of an association between IAs and behavioural or cognitive phenotypes was found, suggesting that IAs are not clearly implicated in these pathologies.

Evolutionary mechanisms shaping the genomic structure of the Williams-Beuren syndrome chromosomal region at human 7q11.23.

About 5% of the human genome consists of segmental duplications or low-copy repeats, which are large, highly homologous (>95%) fragments of sequence. It has been estimated that these segmental duplications emerged during the past approximately 35 million years (Myr) of human evolution and that they correlate with chromosomal rearrangements. Williams-Beuren syndrome (WBS) is a segmental aneusomy syndrome that is the result of a frequent de novo deletion at 7q11.23, mediated by large (approximately 400-kb) region-specific complex segmental duplications composed of different blocks. We have precisely defined the structure of the segmental duplications on human 7q11.23 and characterized the copy number and structure of the orthologous regions in other primates (macaque, orangutan, gorilla, and chimpanzee). Our data indicate a recent origin and rapid evolution of the 7q11.23 segmental duplications, starting before the diversification of hominoids (approximately 12-16 million years ago [Mya]), with species-specific duplications and intrachromosomal rearrangements that lead to significant differences among those genomes. Alu sequences are located at most edges of the large hominoid-specific segmental duplications, suggesting that they might have facilitated evolutionary rearrangements. We propose a mechanistic model based on Alu-mediated duplicated transposition along with nonallelic homologous recombination for the generation and local expansion of the segmental duplications. The extraordinary rate of evolutionary turnover of this region, rich in segmental duplications, results in important genomic variation among hominoid species, which could be of functional relevance and predispose to disease.