Showing papers by "Stylianos E. Antonarakis published in 2018"

Somatic Activating KRAS Mutations in Arteriovenous Malformations of the Brain

Abstract: Background Sporadic arteriovenous malformations of the brain, which are morphologically abnormal connections between arteries and veins in the brain vasculature, are a leading cause of hemorrhagic stroke in young adults and children. The genetic cause of this rare focal disorder is unknown. Methods We analyzed tissue and blood samples from patients with arteriovenous malformations of the brain to detect somatic mutations. We performed exome DNA sequencing of tissue samples of arteriovenous malformations of the brain from 26 patients in the main study group and of paired blood samples from 17 of those patients. To confirm our findings, we performed droplet digital polymerase-chain-reaction (PCR) analysis of tissue samples from 39 patients in the main study group (21 with matching blood samples) and from 33 patients in an independent validation group. We interrogated the downstream signaling pathways, changes in gene expression, and cellular phenotype that were induced by activating KRAS mutations,...

Extensive cellular heterogeneity of X inactivation revealed by single-cell allele-specific expression in human fibroblasts.

Abstract: X-chromosome inactivation (XCI) provides a dosage compensation mechanism where, in each female cell, one of the two X chromosomes is randomly silenced. However, some genes on the inactive X chromosome and outside the pseudoautosomal regions escape from XCI and are expressed from both alleles (escapees). We investigated XCI at single-cell resolution combining deep single-cell RNA sequencing with whole-genome sequencing to examine allelic-specific expression in 935 primary fibroblast and 48 lymphoblastoid single cells from five female individuals. In this framework we integrated an original method to identify and exclude doublets of cells. In fibroblast cells, we have identified 55 genes as escapees including five undescribed escapee genes. Moreover, we observed that all genes exhibit a variable propensity to escape XCI in each cell and cell type and that each cell displays a distinct expression profile of the escapee genes. A metric, the Inactivation Score-defined as the mean of the allelic expression profiles of the escapees per cell-enables us to discover a heterogeneous and continuous degree of cellular XCI with extremes represented by "inactive" cells, i.e., cells exclusively expressing the escaping genes from the active X chromosome and "escaping" cells expressing the escapees from both alleles. We found that this effect is associated with cell-cycle phases and, independently, with the XIST expression level, which is higher in the quiescent phase (G0). Single-cell allele-specific expression is a powerful tool to identify novel escapees in different tissues and provide evidence of an unexpected cellular heterogeneity of XCI.

Biallelic variants in KIF14 cause intellectual disability with microcephaly.

Abstract: Kinesin proteins are critical for various cellular functions such as intracellular transport and cell division, and many members of the family have been linked to monogenic disorders and cancer. We report eight individuals with intellectual disability and microcephaly from four unrelated families with parental consanguinity. In the affected individuals of each family, homozygosity for likely pathogenic variants in KIF14 were detected; two loss-of-function (p.Asn83Ilefs*3 and p.Ser1478fs), and two missense substitutions (p.Ser841Phe and p.Gly459Arg). KIF14 is a mitotic motor protein that is required for spindle localization of the mitotic citron rho-interacting kinase, CIT, also mutated in microcephaly. Our results demonstrate the involvement of KIF14 in development and reveal a wide phenotypic variability ranging from fetal lethality to moderate developmental delay and microcephaly.

Slightly deleterious genomic variants and transcriptome perturbations in Down syndrome embryonic selection

Abstract: The majority of aneuploid fetuses are spontaneously miscarried. Nevertheless, some aneuploid individuals survive despite the strong genetic insult. Here, we investigate if the survival probability of aneuploid fetuses is affected by the genome-wide burden of slightly deleterious variants. We analyzed two cohorts of live-born Down syndrome individuals (388 genotyped samples and 16 fibroblast transcriptomes) and observed a deficit of slightly deleterious variants on Chromosome 21 and decreased transcriptome-wide variation in the expression level of highly constrained genes. We interpret these results as signatures of embryonic selection, and propose a genetic handicap model whereby an individual bearing an extremely severe deleterious variant (such as aneuploidy) could escape embryonic lethality if the genome-wide burden of slightly deleterious variants is sufficiently low. This approach can be used to study the composition and effect of the numerous slightly deleterious variants in humans and model organisms.

Bi-allelic Loss-of-Function Variants in DNMBP Cause Infantile Cataracts.

Abstract: Infantile and childhood-onset cataracts form a heterogeneous group of disorders; among the many genetic causes, numerous pathogenic variants in additional genes associated with autosomal-recessive infantile cataracts remain to be discovered. We identified three consanguineous families affected by bilateral infantile cataracts. Using exome sequencing, we found homozygous loss-of-function variants in DNMBP: nonsense variant c.811C>T (p.Arg271∗) in large family F385 (nine affected individuals; LOD score = 5.18 at θ = 0), frameshift deletion c.2947_2948del (p.Asp983∗) in family F372 (two affected individuals), and frameshift variant c.2852_2855del (p.Thr951Metfs∗41) in family F3 (one affected individual). The phenotypes of all affected individuals include infantile-onset cataracts. RNAi-mediated knockdown of the Drosophila ortholog still life (sif), enriched in lens-secreting cells, affects the development of these cells as well as the localization of E-cadherin, alters the distribution of septate junctions in adjacent cone cells, and leads to a ∼50% reduction in electroretinography amplitudes in young flies. DNMBP regulates the shape of tight junctions, which correspond to the septate junctions in invertebrates, as well as the assembly pattern of E-cadherin in human epithelial cells. E-cadherin has an important role in lens vesicle separation and lens epithelial cell survival in humans. We therefore conclude that DNMBP loss-of-function variants cause infantile-onset cataracts in humans.

Biallelic loss of function variants in PPP1R21 cause a neurodevelopmental syndrome with impaired endocytic function

Abstract: Next-generation sequencing (NGS) has been instrumental in solving the genetic basis of rare inherited diseases, especially neurodevelopmental syndromes. However, functional workup is essential for precise phenotype definition and to understand the underlying disease mechanisms. Using whole exome (WES) and whole genome sequencing (WGS) in four independent families with hypotonia, neurodevelopmental delay, facial dysmorphism, loss of white matter, and thinning of the corpus callosum, we identified four previously unreported homozygous truncating PPP1R21 alleles: c.347delT p.(Ile116Lysfs*25), c.2170_2171insGGTA p.(Ile724Argfs*8), c.1607dupT p.(Leu536Phefs*7), c.2063delA p.(Lys688Serfs*26) and found that PPP1R21 was absent in fibroblasts of an affected individual, supporting the allele's loss of function effect. PPP1R21 function had not been studied except that a large scale affinity proteomics approach suggested an interaction with PIBF1 defective in Joubert syndrome. Our co-immunoprecipitation studies did not confirm this but in contrast defined the localization of PPP1R21 to the early endosome. Consistent with the subcellular expression pattern and the clinical phenotype exhibiting features of storage diseases, we found patient fibroblasts exhibited a delay in clearance of transferrin-488 while uptake was normal. In summary, we delineate a novel neurodevelopmental syndrome caused by biallelic PPP1R21 loss of function variants, and suggest a role of PPP1R21 within the endosomal sorting process or endosome maturation pathway.

Visual impairment and progressive phthisis bulbi caused by recessive pathogenic variant in MARK3.

Abstract: Developmental eye defects often severely reduce vision. Despite extensive efforts, for a substantial fraction of these cases the molecular causes are unknown. Recessive eye disorders are frequent in consanguineous populations and such large families with multiple affected individuals provide an opportunity to identify recessive causative genes. We studied a Pakistani consanguineous family with three affected individuals with congenital vision loss and progressive eye degeneration. The family was analyzed by exome sequencing of one affected individual and genotyping of all family members. We have identified a non-synonymous homozygous variant (NM_001128918.2: c.1708C > G: p.Arg570Gly) in the MARK3 gene as the likely cause of the phenotype. Given that MARK3 is highly conserved in flies (I: 55%; S: 67%) we knocked down the MARK3 homologue, par-1, in the eye during development. This leads to a significant reduction in eye size, a severe loss of photoreceptors and loss of vision based on electroretinogram (ERG) recordings. Expression of the par-1 p.Arg792Gly mutation (equivalent to the MARK3 variant found in patients) in developing fly eyes also induces loss of eye tissue and reduces the ERG signals. The data in flies and human indicate that the MARK3 variant corresponds to a loss of function. We conclude that the identified mutation in MARK3 establishes a new gene-disease link, since it likely causes structural abnormalities during eye development and visual impairment in humans, and that the function of MARK3/par-1 is evolutionarily conserved in eye development.

Genomic, Proteomic and Phenotypic Heterogeneity in HeLa Cells across Laboratories: Implications for Reproducibility of Research Results

Abstract: The independent reproduction of research results is a cornerstone of experimental research, yet it is beset by numerous challenges, including the quality and veracity of reagents and materials. Much of life science research depends on life materials, including human tissue culture cells. In this study we aimed at determining the degree of variability in the molecular makeup and the ensuing phenotypic consequences in commonly used human tissue culture cells. We collected 14 stock HeLa aliquots from 13 different laboratories across the globe, cultured them in uniform conditions and profiled the genome-wide copy numbers, mRNAs, proteins and protein turnover rates via genomic techniques and SWATH mass spectrometry, respectively. We also phenotyped each cell line with respect to the ability of transfected Let7 mimics to modulate Salmonella infection. We discovered significant heterogeneity between HeLa variants, especially between lines of the CCL2 and Kyoto variety. We also observed progressive divergence within a specific cell line over 50 successive passages. From the aggregate multi-omic datasets we quantified the response of the cells to genomic variability across the transcriptome and proteome. We discovered organelle-specific proteome remodeling and buffering of protein abundance by protein complex stoichiometry, mediated by the adaptation of protein turnover rates. By associating quantitative proteotype and phenotype measurements we identified protein patterns that explained the varying response of the different cell lines to Salmonella infection. Altogether the results indicate a striking degree of genomic variability, the rapid evolution of genomic variability in culture and its complex translation into distinctive expressed molecular and phenotypic patterns. The results have broad implications for the interpretation and reproducibility of research results obtained from HeLa cells and provide important basis for a general discussion of the value and requirements for communicating research results obtained from human tissue culture cells.

Novel NEXMIF pathogenic variant in a boy with severe autistic features, intellectual disability, and epilepsy, and his mildly affected mother.

Abstract: Intellectual disability (ID) and autism spectrum disorders are complex neurodevelopmental disorders occurring among all ethnic and socioeconomic groups. Pathogenic variants in the neurite extension and migration factor (NEXMIF) gene (formerly named KIAA2022) on the X chromosome are responsible for ID, autistic behavior, epilepsy, or dysmorphic features in males. Most affected females described had a milder phenotype or were asymptomatic obligate carriers. We report here for the first time mother-to-son transmission of a novel NEXMIF truncating variant without X-inactivation skewing in the blood. Truncating gene variant leads to symptomatic mother to severely affected son transmission. Our findings emphasize that NEXMIF sequencing should be strongly considered in patients with unexplained autism spectrum disorder, ID, and epilepsy, irrespective of gender. Such testing could increase our knowledge of the pathogenicity of NEXMIF variants and improve genetic counseling.

Time and space dimensions of gene dosage imbalance of aneuploidies revealed by single cell transcriptomes.

Abstract: The mechanisms underlying cellular and organismal phenotypes due to copy number alterations (CNA) are not fully understood Aneuploidy is a major source of gene dosage imbalance due to CNA and viable human trisomies are model disorders of altered gene expression To understand the cellular impact of gene dosage imbalance, we studied gene and allele specific expression (ASE) of 9668 single-cell fibroblasts in trisomies T21, T18, T13 and T8 To limit the bias of interindividual noise, all comparisons between euploid and trisomic single-cells were performed on an isogenic setting for all trisomies studied Initially we examined 928 single cells with deep RNA-Seq For T21 we used fibroblasts from one pair of monozygotic twins discordant for T21 and from mosaic T21 For T18, T13 and T8 we analyzed single cells from mosaic individuals Single-cell analyses revealed inconsistencies concerning the overexpression of some genes observed in differential trisomic vs euploid bulk RNAseq while this imbalance was not detectable in trisomic vs euploid single cells Moreover, ASE profiling of all single cells uncovered a substantial monoallelic pattern of expression in the trisomic fraction of the genome By classifying genes according to the level of mono and bi-allelic transcription, we have observed that, for genes with monoallelic and low-to-average expression, the altered gene dosage is mainly due to the higher fraction of cells simultaneously expressing these genes in the trisomic samples These results were confirmed in a further experiment of 8740 single fibroblasts from the monozygotic twins discordant for T21 samples We conclude that gene dosage imbalance is of bidimensional nature: over time (simultaneous expression of all alleles resulting in increased accumulation of RNA of copy altered genes in each single cell) as previously stated, and over space (increased fraction of cells simultaneously expressing copy altered genes) These results strongly suggest that each class of genes contributes to the phenotypic variability of trisomies according to its temporal and spatial behavior and propose an improved model to understand the effects of copy number alterations

Extensive cellular heterogeneity of X inactivation revealed by single-cell allele-specific expression in human fibroblasts.

Abstract: In eutherian mammals, X chromosome inactivation (XCI) provides a dosage compensation mechanism where in each female cell one of the two X chromosomes is randomly silenced. However, some genes on the inactive X chromosome and outside the pseudoautosomal regions escape from XCI and are expressed from both alleles (escapees). Given the relevance of the escapees in biology and medicine, we investigated XCI at an unprecedented single-cell resolution. We combined deep single-cell RNA sequencing with whole genome sequencing to examine allelic specific expression (ASE) in 935 primary fibroblast and 48 lymphoblastoid single cells from five female individuals. In this framework we integrated an original method to identify and exclude doublets of cells. We have identified 55 genes as escapees including 5 novel escapee genes. Moreover, we observed that all genes exhibit a variable propensity to escape XCI in each cell and cell type, and that each cell displays a distinct expression profile of the escapee genes. We devised a novel metric, the Inactivation Score (IS), defined as the mean of the allelic expression profiles of the escapees per cell, and discovered a heterogeneous and continuous degree of cellular XCI with extremes represented by "inactive" cells, i.e., exclusively expressing the escaping genes from the active X chromosome, and "escaping" cells, expressing the escapees from both alleles. Intriguingly we found that XIST is the major genetic determinant of IS, and that XIST expression, higher in G0 phase, is negatively correlated with the expression of escapees, inactivated and pseudoautosomal genes. In this study we use single-cell allele specific expression to identify novel escapees in different tissues and provide evidence of an unexpected cellular heterogeneity of XCI driven by a possible regulatory activity of XIST.