Comparison of clustering tools in R for medium-sized 10x Genomics single-cell RNA-sequencing data [version 1; referees: 1 approved, 2 approved with reservations]

Background: The commercially available 10x Genomics protocol to generate droplet-based single-cell RNA-seq (scRNA-seq) data is enjoying growing popularity among researchers. Fundamental to the analysis of such scRNA-seq data is the ability to cluster similar or same cells into non-overlapping groups. Many competing methods have been proposed for this task, but there is currently little guidance with regards to which method to use. Methods: Here we use one gold standard 10x Genomics dataset, generated from the mixture of three cell lines, as well as three silver standard 10x Genomics datasets generated from peripheral blood mononuclear cells to examine not only the accuracy but also robustness of a dozen methods. Results: We found that some methods, including Seurat and Cell Ranger, outperform other methods, although performance seems to be dependent on the complexity of the studied system. Furthermore, we found that solutions produced by different methods have little in common with each other. Conclusions: In light of this, we conclude that the choice of clustering tool crucially determines interpretation of scRNA-seq data generated by 10x Genomics. Hence practitioners and consumers should remain vigilant about the outcome of 10x Genomics scRNA-seq analysis

Analysis of genome-wide association study data using the protein knowledge base

Article 9

Neuropathology of childhood‐onset basal ganglia degeneration caused by mutation of VAC14

ObjectiveTo characterize the clinical features and neuropathology associated with recessive VAC14 mutations.MethodsWhole‐exome sequencing was used to identify the genetic etiology of a rapidly progressive neurological disease presenting in early childhood in two deceased siblings with distinct neuropathological features on post mortem examination.ResultsWe identified compound heterozygous variants in VAC14 in two deceased siblings with early childhood onset of severe, progressive dystonia, and neurodegeneration. Their clinical phenotype is consistent with the VAC14–related childhood‐onset, striatonigral degeneration recently described in two unrelated children. Post mortem examination demonstrated prominent vacuolation associated with degenerating neurons in the caudate nucleus, putamen, and globus pallidus, similar to previously reported ex vivo vacuoles seen in the late‐endosome/lysosome of VAC14‐deficient neurons. We identified upregulation of ubiquitinated granules within the cell cytoplasm and lysosomal‐associated membrane protein (LAMP2) around the vacuole edge to suggest a process of vacuolation of lysosomal structures associated with active autophagocytic‐associated neuronal degeneration.InterpretationOur findings reveal a distinct clinicopathological phenotype associated with recessive VAC14 mutations.Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/142276/1/acn3487_am.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/142276/2/acn3487.pd

Genetic disruption of serine biosynthesis is a key driver of macular telangiectasia type 2 aetiology and progression.

BACKGROUND: Macular telangiectasia type 2 (MacTel) is a rare, heritable and largely untreatable retinal disorder, often comorbid with diabetes. Genetic risk loci subtend retinal vascular calibre and glycine/serine/threonine metabolism genes. Serine deficiency may contribute to MacTel via neurotoxic deoxysphingolipid production; however, an independent vascular contribution is also suspected. Here, we use statistical genetics to dissect the causal mechanisms underpinning this complex disease. METHODS: We integrated genetic markers for MacTel, vascular and metabolic traits, and applied Mendelian randomisation and conditional and interaction genome-wide association analyses to discover the causal contributors to both disease and spatial retinal imaging sub-phenotypes. RESULTS: Genetically induced serine deficiency is the primary causal metabolic driver of disease occurrence and progression, with a lesser, but significant, causal contribution of type 2 diabetes genetic risk. Conversely, glycine, threonine and retinal vascular traits are unlikely to be causal for MacTel. Conditional regression analysis identified three novel disease loci independent of endogenous serine biosynthetic capacity. By aggregating spatial retinal phenotypes into endophenotypes, we demonstrate that SNPs constituting independent risk loci act via related endophenotypes. CONCLUSIONS: Follow-up studies after GWAS integrating publicly available data with deep phenotyping are still rare. Here, we describe such analysis, where we integrated retinal imaging data with MacTel and other traits genomics data to identify biochemical mechanisms likely causing this disorder. Our findings will aid in early diagnosis and accurate prognosis of MacTel and improve prospects for effective therapeutic intervention. Our integrative genetics approach also serves as a useful template for post-GWAS analyses in other disorders

Whole exome sequencing combined with linkage analysis identifies a novel 3 bp deletion in NR5A1

Disorders of sex development (DSDs) encompass a broad spectrum of conditions affecting the development of the gonads and genitalia. The underlying causes for DSDs include gain or loss of function variants in genes responsible for gonad development or steroidogenesis. Most patients with DSD have an unknown genetic etiology and cannot be given an

Comparing genotyping algorithms for Illumina's Infinium whole-genome SNP BeadChips

Background: Illumina's Infinium SNP BeadChips are extensively used in both small and large-scale genetic studies. A fundamental step in any analysis is the processing of raw allele A and allele B intensities from each SNP into genotype calls (AA, AB, BB). Various algorithms which make use of different statistical models are available for this task. We compare four methods (GenCall, Illuminus, GenoSNP and CRLMM) on data where the true genotypes are known in advance and data from a recently published genome-wide association study.Results: In general, differences in accuracy are relatively small between the methods evaluated, although CRLMM and GenoSNP were found to consistently outperform GenCall. The performance of Illuminus is heavily dependent on sample size, with lower no call rates and improved accuracy as the number of samples available increases. For X chromosome SNPs, methods with sex-dependent models (Illuminus, CRLMM) perform better than methods which ignore gender information (GenCall, GenoSNP). We observe that CRLMM and GenoSNP are more accurate at calling SNPs with low minor allele frequency than GenCall or Illuminus. The sample quality metrics from each of the four methods were found to have a high level of agreement at flagging samples with unusual signal characteristics.Conclusions: CRLMM, GenoSNP and GenCall can be applied with confidence in studies of any size, as their performance was shown to be invariant to the number of samples available. Illuminus on the other hand requires a larger number of samples to achieve comparable levels of accuracy and its use in smaller studies (50 or fewer individuals) is not recommended

A Mouse Model of Harlequin Ichthyosis Delineates a Key Role for Abca12 in Lipid Homeostasis

Harlequin Ichthyosis (HI) is a severe and often lethal hyperkeratotic skin disease caused by mutations in the ABCA12 transport protein. In keratinocytes, ABCA12 is thought to regulate the transfer of lipids into small intracellular trafficking vesicles known as lamellar bodies. However, the nature and scope of this regulation remains unclear. As part of an original recessive mouse ENU mutagenesis screen, we have identified and characterised an animal model of HI and showed that it displays many of the hallmarks of the disease including hyperkeratosis, loss of barrier function, and defects in lipid homeostasis. We have used this model to follow disease progression in utero and present evidence that loss of Abca12 function leads to premature differentiation of basal keratinocytes. A comprehensive analysis of lipid levels in mutant epidermis demonstrated profound defects in lipid homeostasis, illustrating for the first time the extent to which Abca12 plays a pivotal role in maintaining lipid balance in the skin. To further investigate the scope of Abca12's activity, we have utilised cells from the mutant mouse to ascribe direct transport functions to the protein and, in doing so, we demonstrate activities independent of its role in lamellar body function. These cells have severely impaired lipid efflux leading to intracellular accumulation of neutral lipids. Furthermore, we identify Abca12 as a mediator of Abca1-regulated cellular cholesterol efflux, a finding that may have significant implications for other diseases of lipid metabolism and homeostasis, including atherosclerosis

Progressive Myoclonus Epilepsy Caused by a Homozygous Splicing Variant of SLC7A6OS

Exome sequencing was performed in 2 unrelated families with progressive myoclonus epilepsy. Affected individuals from both families shared a rare, homozygous c.191A > G variant affecting a splice site in SLC7A6OS. Analysis of cDNA from lymphoblastoid cells demonstrated partial splice site abolition and the creation of an abnormal isoform. Quantitative reverse transcriptase polymerase chain reaction and Western blot showed a marked reduction of protein expression. Haplotype analysis identified a similar to 0.85cM shared genomic region on chromosome 16q encompassing the c.191A > G variant, consistent with a distant ancestor common to both families. Our results suggest that biallelic loss-of-function variants in SLC7A6OS are a novel genetic cause of progressive myoclonus epilepsy. ANN NEUROL 2020Peer reviewe

Increasingly inbred and fragmented populations of Plasmodium vivax associated with the eastward decline in malaria transmission across the Southwest Pacific

The human malaria parasite Plasmodium vivax is more resistant to malaria control strategies than Plasmodium falciparum, and maintains high genetic diversity even when transmission is low. To investigate whether declining P. vivax transmission leads to increasing population structure that would facilitate elimination, we genotyped samples from across the Southwest Pacific region, which experiences an eastward decline in malaria transmission, as well as samples from two time points at one site (Tetere, Solomon Islands) during intensified malaria control. Analysis of 887 P. vivax microsatellite haplotypes from hyperendemic Papua New Guinea (PNG, n = 443), meso-hyperendemic Solomon Islands (n = 420), and hypoendemic Vanuatu (n = 24) revealed increasing population structure and multilocus linkage disequilibrium yet a modest decline in diversity as transmission decreases over space and time. In Solomon Islands, which has had sustained control efforts for 20 years, and Vanuatu, which has experienced sustained low transmission for many years, significant population structure was observed at different spatial scales. We conclude that control efforts will eventually impact P. vivax population structure and with sustained pressure, populations may eventually fragment into a limited number of clustered foci that could be targeted for elimination