Ticks from the genus Rhipicephalus have enormous global economic impact as ectoparasites of cattle. Rhipicephalus microplus and Rhipicephalus annulatus are known to harbor infectious pathogens such as Babesia bovis, Babesia bigemina, and Anaplasma marginale. Having reference quality genomes of these ticks would advance research to identify druggable targets for chemical entities with acaricidal activity and refine anti-tick vaccine approaches. We sequenced and assembled the genomes of R. microplus and R. annulatus, using Pacific Biosciences and HiSeq 4000 technologies on very high molecular weight genomic DNA. We used 22 and 29 SMRT cells on the Pacific Biosciences Sequel for R. microplus and R. annulatus, respectively, and 3 lanes of the Illumina HiSeq 4000 platform for each tick. The PacBio sequence yields for R. microplus and R. annulatus were 21.0 and 27.9 million subreads, respectively, which were assembled with Canu v. 1.7. The final Canu assemblies consisted of 92,167 and 57,796 contigs with an average contig length of 39,249 and 69,055 bp for R. microplus and R. annulatus, respectively. Annotated genome quality was assessed by BUSCO analysis to provide quantitative measures for each assembled genome. Over 82% and 92% of the 1066 member BUSCO gene set was found in the assembled genomes of R. microplus and R. annulatus, respectively. For R. microplus, only 189 of the 1066 BUSCO genes were missing and only 140 were present in a fragmented condition. For R. annulatus, only 75 of the BUSCO genes were missing and only 109 were present in a fragmented condition. The raw sequencing reads and the assembled contigs/scaffolds are archived at the National Center for Biotechnology Information.

Angelman syndrome (AS) is a neurodevelopmental disorder with unique behavioral phenotypes, seizures, and distinctive electroencephalographic (EEG) patterns. Recent studies identified motor, social communication, and learning and memory deficits in a CRISPR engineered rat model with a complete maternal deletion of the Ube3a gene. It is unknown whether this model recapitulates other aspects of the clinical disorder. We report here the effect of Ube3a maternal deletion in the rat on epileptiform activity, seizure threshold, and quantitative EEG. Using video-synchronized EEG (vEEG) monitoring, we assessed spectral power and epileptiform activity early postnatally through adulthood. While EEG power was similar to wild-type (WT) at 1.5 weeks postnatally, at all other ages analyzed, our findings were similar to the AS phenotype in mice and humans with significantly increased δ power. Analysis of epileptiform activity in juvenile and adult rats showed increased time spent in epileptiform activity in AS compared with WT rats. We evaluated seizure threshold using pentylenetetrazol (PTZ), audiogenic stimulus, and hyperthermia to provoke febrile seizures (FSs). Behavioral seizure scoring following PTZ induction revealed no difference in seizure threshold in AS rats, however behavioral recovery from the PTZ-induced seizure was longer in the adult group with significantly increased hippocampal epileptiform activity during this phase. When exposed to hyperthermia, AS rat pups showed a significantly lower temperature threshold to first seizure than WT. Our findings highlight an age-dependence for the EEG and epileptiform phenotypes in a preclinical model of AS, and support the use of quantitative EEG and increased δ power as a potential biomarker of AS.

Focal epilepsies are the largest epilepsy subtype and associated with significant morbidity. Somatic variation is a newly recognized genetic mechanism underlying a subset of focal epilepsies, but little is known about the processes through which somatic mosaicism causes seizures, the cell types carrying the pathogenic variants, or their developmental origin. Meanwhile, the inception of single cell biology has completely revolutionized the study of neurological diseases and has the potential to answer some of these key questions. Focusing on single cell genomics, transcriptomics, and epigenomics in focal epilepsy research, circumvents the averaging artifact associated with studying bulk brain tissue and offers the kind of granularity that is needed for investigating the consequences of somatic mosaicism. Here we have provided a brief overview of some of the most developed single cell techniques and the major considerations around applying them to focal epilepsy research.

The cognitive abilities of humans are distinctive among primates, but their molecular and cellular substrates are poorly understood. We used comparative single-nucleus transcriptomics to analyze samples of the middle temporal gyrus (MTG) from adult humans, chimpanzees, gorillas, rhesus macaques, and common marmosets to understand human-specific features of the neocortex. Human, chimpanzee, and gorilla MTG showed highly similar cell-type composition and laminar organization as well as a large shift in proportions of deep-layer intratelencephalic-projecting neurons compared with macaque and marmoset MTG. Microglia, astrocytes, and oligodendrocytes had more-divergent expression across species compared with neurons or oligodendrocyte precursor cells, and neuronal expression diverged more rapidly on the human lineage. Only a few hundred genes showed human-specific patterning, suggesting that relatively few cellular and molecular changes distinctively define adult human cortical structure.

While germline copy-number variants (CNVs) contribute to schizophrenia (SCZ) risk, the contribution of somatic CNVs (sCNVs)-present in some but not all cells-remains unknown. We identified sCNVs using blood-derived genotype arrays from 12,834 SCZ cases and 11,648 controls, filtering sCNVs at loci recurrently mutated in clonal blood disorders. Likely early-developmental sCNVs were more common in cases (0.91%) than controls (0.51%, p = 2.68e-4), with recurrent somatic deletions of exons 1-5 of the NRXN1 gene in five SCZ cases. Hi-C maps revealed ectopic, allele-specific loops forming between a potential cryptic promoter and non-coding cis-regulatory elements upon 5' deletions in NRXN1. We also observed recurrent intragenic deletions of ABCB11, encoding a transporter implicated in anti-psychotic response, in five treatment-resistant SCZ cases and showed that ABCB11 is specifically enriched in neurons forming mesocortical and mesolimbic dopaminergic projections. Our results indicate potential roles of sCNVs in SCZ risk.

Alzheimer's disease (AD) is an age-associated neurodegenerative disorder characterized by progressive neuronal loss and pathological accumulation of the misfolded proteins amyloid-β and tau1,2. Neuroinflammation mediated by microglia and brain-resident macrophages plays a crucial role in AD pathogenesis1-5, though the mechanisms by which age, genes, and other risk factors interact remain largely unknown. Somatic mutations accumulate with age and lead to clonal expansion of many cell types, contributing to cancer and many non-cancer diseases6,7. Here we studied somatic mutation in normal aged and AD brains by three orthogonal methods and in three independent AD cohorts. Analysis of bulk RNA sequencing data from 866 samples from different brain regions revealed significantly higher (~two-fold) overall burdens of somatic single-nucleotide variants (sSNVs) in AD brains compared to age-matched controls. Molecular-barcoded deep (>1000X) gene panel sequencing of 311 prefrontal cortex samples showed enrichment of sSNVs and somatic insertions and deletions (sIndels) in cancer driver genes in AD brain compared to control, with recurrent, and often multiple, mutations in genes implicated in clonal hematopoiesis (CH)8,9. Pathogenic sSNVs were enriched in CSF1R+ microglia of AD brains, and the high proportion of microglia (up to 40%) carrying some sSNVs in cancer driver genes suggests mutation-driven microglial clonal expansion (MiCE). Analysis of single-nucleus RNA sequencing (snRNAseq) from temporal neocortex of 62 additional AD cases and controls exhibited nominally increased mosaic chromosomal alterations (mCAs) associated with CH10,11. Microglia carrying mCA showed upregulated pro-inflammatory genes, resembling the transcriptomic features of disease-associated microglia (DAM) in AD. Our results suggest that somatic driver mutations in microglia are common with normal aging but further enriched in AD brain, driving MiCE with inflammatory and DAM signatures. Our findings provide the first insights into microglial clonal dynamics in AD and identify potential new approaches to AD diagnosis and therapy.

Mutations in several genes encoding centrosomal proteins dramatically decrease the size of the human brain. We show that Aspm (abnormal spindle-like, microcephaly-associated) and Wdr62 (WD repeat-containing protein 62) interact genetically to control brain size, with mice lacking Wdr62, Aspm, or both showing gene dose-related centriole duplication defects that parallel the severity of the microcephaly and increased ectopic basal progenitors, suggesting premature delamination from the ventricular zone. Wdr62 and Aspm localize to the proximal end of the mother centriole and interact physically, with Wdr62 required for Aspm localization, and both proteins, as well as microcephaly protein Cep63, required to localize CENPJ/CPAP/Sas-4, a final common target. Unexpectedly, Aspm and Wdr62 are required for normal apical complex localization and apical epithelial structure, providing a plausible unifying mechanism for the premature delamination and precocious differentiation of progenitors. Together, our results reveal links among centrioles, apical proteins, and cell fate, and illuminate how alterations in these interactions can dynamically regulate brain size.

Comparative analyses have identified genomic regions potentially involved in human evolution but do not directly assess function. Human accelerated regions (HARs) represent conserved genomic loci with elevated divergence in humans. If some HARs regulate human-specific social and behavioral traits, then mutations would likely impact cognitive and social disorders. Strikingly, rare biallelic point mutations-identified by whole-genome and targeted "HAR-ome" sequencing-showed a significant excess in individuals with ASD whose parents share common ancestry compared to familial controls, suggesting a contribution in 5% of consanguineous ASD cases. Using chromatin interaction sequencing, massively parallel reporter assays (MPRA), and transgenic mice, we identified disease-linked, biallelic HAR mutations in active enhancers for CUX1, PTBP2, GPC4, CDKL5, and other genes implicated in neural function, ASD, or both. Our data provide genetic evidence that specific HARs are essential for normal development, consistent with suggestions that their evolutionary changes may have altered social and/or cognitive behavior. PAPERCLIP.

We sought to explore the genetic and molecular causes of Troyer syndrome, one of several complicated hereditary spastic paraplegias (HSPs). Troyer syndrome had been thought to be restricted to the Amish; however, we identified 2 Omani families with HSP, short stature, dysarthria and developmental delay-core features of Troyer syndrome-and a novel mutation in the SPG20 gene, which is also mutated in the Amish. In addition, we analyzed SPG20 expression throughout development to infer how disruption of this gene might generate the constellation of developmental and degenerative Troyer syndrome phenotypes.

The dorsal thalamus acts as a gateway and modulator for information going to and from the cerebral cortex. This activity requires the formation of reciprocal topographic axon connections between thalamus and cortex. The axons grow along a complex multistep pathway, making sharp turns, crossing expression boundaries, and encountering intermediate targets. However, the cellular and molecular components mediating these steps remain poorly understood.

Hypogonadism, alopecia, diabetes mellitus, mental retardation, and extrapyramidal syndrome (also referenced as Woodhouse-Sakati syndrome) is a rare autosomal recessive multisystemic disorder. We have identified a founder mutation consisting of a single base-pair deletion in C2orf37 in eight families of Saudi origin. Three other loss-of-function mutations were subsequently discovered in patients of different ethnicities. The gene encodes a nucleolar protein of unknown function, and the cellular phenotype observed in patient lymphoblasts implicates a role for the nucleolus in the pathogenesis of this disease. Our findings expand the list of human disorders linked to the nucleolus and further highlight the developmental and/or maintenance functions of this organelle.

The growing body of DNA microarray data has the potential to advance our understanding of the molecular basis of disease. However annotating microarray datasets with clinically useful information is not always possible, as this often requires access to detailed patient records. In this study we introduce GLAD, a new Semi-Supervised Learning (SSL) method for combining independent annotated datasets and unannotated datasets with the aim of identifying more robust sample classifiers. In our method, independent models are developed using subsets of genes for the annotated and unannotated datasets. These models are evaluated according to a scoring function that incorporates terms for classification accuracy on annotated data, and relative cluster separation in unannotated data. Improved models are iteratively generated using a genetic algorithm feature selection technique. Our results show that the addition of unannotated data into training, significantly improves classifier robustness.

Charged multivesicular body protein 1A (CHMP1A; also known as chromatin-modifying protein 1A) is a member of the ESCRT-III (endosomal sorting complex required for transport-III) complex but is also suggested to localize to the nuclear matrix and regulate chromatin structure. Here, we show that loss-of-function mutations in human CHMP1A cause reduced cerebellar size (pontocerebellar hypoplasia) and reduced cerebral cortical size (microcephaly). CHMP1A-mutant cells show impaired proliferation, with increased expression of INK4A, a negative regulator of stem cell proliferation. Chromatin immunoprecipitation suggests loss of the normal INK4A repression by BMI in these cells. Morpholino-based knockdown of zebrafish chmp1a resulted in brain defects resembling those seen after bmi1a and bmi1b knockdown, which were partially rescued by INK4A ortholog knockdown, further supporting links between CHMP1A and BMI1-mediated regulation of INK4A. Our results suggest that CHMP1A serves as a critical link between cytoplasmic signals and BMI1-mediated chromatin modifications that regulate proliferation of central nervous system progenitor cells.

Mutations in several known or putative glycosyltransferases cause glycosylation defects in α-dystroglycan (α-DG), an integral component of the dystrophin glycoprotein complex. The hypoglycosylation reduces the ability of α-DG to bind laminin and other extracellular matrix ligands and is responsible for the pathogenesis of an inherited subset of muscular dystrophies known as the dystroglycanopathies. By exome and Sanger sequencing we identified two individuals affected by a dystroglycanopathy with mutations in β-1,3-N-acetylgalactosaminyltransferase 2 (B3GALNT2). B3GALNT2 transfers N-acetyl galactosamine (GalNAc) in a β-1,3 linkage to N-acetyl glucosamine (GlcNAc). A subsequent study of a separate cohort of individuals identified recessive mutations in four additional cases that were all affected by dystroglycanopathy with structural brain involvement. We show that functional dystroglycan glycosylation was reduced in the fibroblasts and muscle (when available) of these individuals via flow cytometry, immunoblotting, and immunocytochemistry. B3GALNT2 localized to the endoplasmic reticulum, and this localization was perturbed by some of the missense mutations identified. Moreover, knockdown of b3galnt2 in zebrafish recapitulated the human congenital muscular dystrophy phenotype with reduced motility, brain abnormalities, and disordered muscle fibers with evidence of damage to both the myosepta and the sarcolemma. Functional dystroglycan glycosylation was also reduced in the b3galnt2 knockdown zebrafish embryos. Together these results demonstrate a role for B3GALNT2 in the glycosylation of α-DG and show that B3GALNT2 mutations can cause dystroglycanopathy with muscle and brain involvement.

In mammals, centrioles participate in brain development, and human mutations affecting centriole duplication cause microcephaly. Here, we identify a role for the mammalian homologue of yeast SFI1, involved in the duplication of the yeast spindle pole body, as a critical regulator of centriole duplication in mammalian cells. Mammalian SFI1 interacts with USP9X, a deubiquitylase associated with human syndromic mental retardation. SFI1 localizes USP9X to the centrosome during S phase to deubiquitylate STIL, a critical regulator of centriole duplication. USP9X-mediated deubiquitylation protects STIL from degradation. Consistent with a role for USP9X in stabilizing STIL, cells from patients with USP9X loss-of-function mutations have reduced STIL levels. Together, these results demonstrate that SFI1 is a centrosomal protein that localizes USP9X to the centrosome to stabilize STIL and promote centriole duplication. We propose that the USP9X protection of STIL to facilitate centriole duplication underlies roles of both proteins in human neurodevelopment.

Focal cortical dysplasia (FCD) and hemimegalencephaly (HME) are epileptogenic neurodevelopmental malformations caused by mutations in mTOR pathway genes. Deep sequencing of these genes in FCD/HME brain tissue identified an etiology in 27 of 66 cases (41%). Radiographically indistinguishable lesions are caused by somatic activating mutations in AKT3, MTOR, and PIK3CA and germline loss-of-function mutations in DEPDC5, NPRL2, and TSC1/2, including TSC2 mutations in isolated HME demonstrating a "two-hit" model. Mutations in the same gene cause a disease continuum from FCD to HME to bilateral brain overgrowth, reflecting the progenitor cell and developmental time when the mutation occurred. Single-cell sequencing demonstrated mTOR activation in neurons in all lesions. Conditional Pik3ca activation in the mouse cortex showed that mTOR activation in excitatory neurons and glia, but not interneurons, is sufficient for abnormal cortical overgrowth. These data suggest that mTOR activation in dorsal telencephalic progenitors, in some cases specifically the excitatory neuron lineage, causes cortical dysplasia.

Genes disrupted in human microcephaly (meaning "small brain") define key regulators of neural progenitor proliferation and cell-fate specification. In comparison, genes mutated in human lissencephaly (lissos means smooth and cephalos means brain) highlight critical regulators of neuronal migration. Here, we report two families with extreme microcephaly and grossly simplified cortical gyral structure, a condition referred to as microlissencephaly, and show that they carry homozygous frameshift mutations in NDE1, which encodes a multidomain protein that localizes to the centrosome and mitotic spindle poles. Both human mutations in NDE1 truncate the C-terminal NDE1domains, which are essential for interactions with cytoplasmic dynein and thus for regulation of cytoskeletal dynamics in mitosis and for cell-cycle-dependent phosphorylation of NDE1 by Cdk1. We show that the patient NDE1 proteins are unstable, cannot bind cytoplasmic dynein, and do not localize properly to the centrosome. Additionally, we show that CDK1 phosphorylation at T246, which is within the C-terminal region disrupted by the mutations, is required for cell-cycle progression from the G2 to the M phase. The role of NDE1 in cell-cycle progression probably contributes to the profound neuronal proliferation defects evident in Nde1-null mice and patients with NDE1 mutations, demonstrating the essential role of NDE1 in human cerebral cortical neurogenesis.

Progressive microcephaly is a heterogeneous condition with causes including mutations in genes encoding regulators of neuronal survival. Here, we report the identification of mutations in QARS (encoding glutaminyl-tRNA synthetase [QARS]) as the causative variants in two unrelated families affected by progressive microcephaly, severe seizures in infancy, atrophy of the cerebral cortex and cerebellar vermis, and mild atrophy of the cerebellar hemispheres. Whole-exome sequencing of individuals from each family independently identified compound-heterozygous mutations in QARS as the only candidate causative variants. QARS was highly expressed in the developing fetal human cerebral cortex in many cell types. The four QARS mutations altered highly conserved amino acids, and the aminoacylation activity of QARS was significantly impaired in mutant cell lines. Variants p.Gly45Val and p.Tyr57His were located in the N-terminal domain required for QARS interaction with proteins in the multisynthetase complex and potentially with glutamine tRNA, and recombinant QARS proteins bearing either substitution showed an over 10-fold reduction in aminoacylation activity. Conversely, variants p.Arg403Trp and p.Arg515Trp, each occurring in a different family, were located in the catalytic core and completely disrupted QARS aminoacylation activity in vitro. Furthermore, p.Arg403Trp and p.Arg515Trp rendered QARS less soluble, and p.Arg403Trp disrupted QARS-RARS (arginyl-tRNA synthetase 1) interaction. In zebrafish, homozygous qars loss of function caused decreased brain and eye size and extensive cell death in the brain. Our results highlight the importance of QARS during brain development and that epilepsy due to impairment of QARS activity is unusually severe in comparison to other aminoacyl-tRNA synthetase disorders.

Channelopathies are disorders caused by abnormal ion channel function in differentiated excitable tissues. We discovered a unique neurodevelopmental channelopathy resulting from pathogenic variants in SCN3A, a gene encoding the voltage-gated sodium channel NaV1.3. Pathogenic NaV1.3 channels showed altered biophysical properties including increased persistent current. Remarkably, affected individuals showed disrupted folding (polymicrogyria) of the perisylvian cortex of the brain but did not typically exhibit epilepsy; they presented with prominent speech and oral motor dysfunction, implicating SCN3A in prenatal development of human cortical language areas. The development of this disorder parallels SCN3A expression, which we observed to be highest early in fetal cortical development in progenitor cells of the outer subventricular zone and cortical plate neurons and decreased postnatally, when SCN1A (NaV1.1) expression increased. Disrupted cerebral cortical folding and neuronal migration were recapitulated in ferrets expressing the mutant channel, underscoring the unexpected role of SCN3A in progenitor cells and migrating neurons.

Single cell whole-genome sequencing (scWGS) is providing novel insights into the nature of genetic heterogeneity in normal and diseased cells. However, the whole-genome amplification process required for scWGS introduces biases into the resulting sequencing that can confound downstream analysis. Here, we present a statistical method, with an accompanying package PaSD-qc (Power Spectral Density-qc), that evaluates the properties and quality of single cell libraries. It uses a modified power spectral density to assess amplification uniformity, amplicon size distribution, autocovariance and inter-sample consistency as well as to identify chromosomes with aberrant read-density profiles due either to copy alterations or poor amplification. These metrics provide a standard way to compare the quality of single cell samples as well as yield information necessary to improve variant calling strategies. We demonstrate the usefulness of this tool in comparing the properties of scWGS protocols, identifying potential chromosomal copy number variation, determining chromosomal and subchromosomal regions of poor amplification, and selecting high-quality libraries from low-coverage data for deep sequencing. The software is available free and open-source at https://github.com/parklab/PaSDqc.