Clonal multi-omics reveals Bcor as a negative regulator of emergency dendritic cell development.

Despite advances in single-cell multi-omics, a single stem or progenitor cell can only be tested once. We developed clonal multi-omics, in which daughters of a clone act as surrogates of the founder, thereby allowing multiple independent assays per clone. With SIS-seq, clonal siblings in parallel "sister" assays are examined either for gene expression by RNA sequencing (RNA-seq) or for fate in culture. We identified, and then validated using CRISPR, genes that controlled fate bias for different dendritic cell (DC) subtypes. This included Bcor as a suppressor of plasmacytoid DC (pDC) and conventional DC type 2 (cDC2) numbers during Flt3 ligand-mediated emergency DC development. We then developed SIS-skew to examine development of wild-type and Bcor-deficient siblings of the same clone in parallel. We found Bcor restricted clonal expansion, especially for cDC2s, and suppressed clonal fate potential, especially for pDCs. Therefore, SIS-seq and SIS-skew can reveal the molecular and cellular mechanisms governing clonal fate.

An intronic GAA repeat expansion in FGF14 causes the autosomal-dominant adult-onset ataxia SCA50/ATX-FGF14.

Adult-onset cerebellar ataxias are a group of neurodegenerative conditions that challenge both genetic discovery and molecular diagnosis. In this study, we identified an intronic (GAA) repeat expansion in fibroblast growth factor 14 (FGF14). Genetic analysis of 95 Australian individuals with adult-onset ataxia identified four (4.2%) with (GAA)>300 and a further nine individuals with (GAA)>250. PCR and long-read sequence analysis revealed these were pure (GAA) repeats. In comparison, no control subjects had (GAA)>300 and only 2/311 control individuals (0.6%) had a pure (GAA)>250. In a German validation cohort, 9/104 (8.7%) of affected individuals had (GAA)>335 and a further six had (GAA)>250, whereas 10/190 (5.3%) control subjects had (GAA)>250 but none were (GAA)>335. The combined data suggest (GAA)>335 are disease causing and fully penetrant (p = 6.0 × 10-8, OR = 72 [95% CI = 4.3-1,227]), while (GAA)>250 is likely pathogenic with reduced penetrance. Affected individuals had an adult-onset, slowly progressive cerebellar ataxia with variable features including vestibular impairment, hyper-reflexia, and autonomic dysfunction. A negative correlation between age at onset and repeat length was observed (R2 = 0.44, p = 0.00045, slope = -0.12) and identification of a shared haplotype in a minority of individuals suggests that the expansion can be inherited or generated de novo during meiotic division. This study demonstrates the power of genome sequencing and advanced bioinformatic tools to identify novel repeat expansions via model-free, genome-wide analysis and identifies SCA50/ATX-FGF14 as a frequent cause of adult-onset ataxia.

Genetics of reticular pseudodrusen in age-related macular degeneration.

Reticular pseudodrusen (RPD) are subretinal deposits that, when observed with age-related macular degeneration (AMD), form a distinct phenotype, often associated with late-stage disease. To date, RPD genetic risk associations overlap six well-established AMD-risk regions. Determining RPD-specific underlying genetic causes by using adequate imaging methods should improve our understanding of the pathophysiology of RPD.

A founder event causing a dominant childhood epilepsy survives 800 years through weak selective pressure.

Genetic epilepsy with febrile seizures plus (GEFS+) is an autosomal dominant familial epilepsy syndrome characterized by distinctive phenotypic heterogeneity within families. The SCN1B c.363C>G (p.Cys121Trp) variant has been identified in independent, multi-generational families with GEFS+. Although the variant is present in population databases (at very low frequency), there is strong clinical, genetic, and functional evidence to support pathogenicity. Recurrent variants may be due to a founder event in which the variant has been inherited from a common ancestor. Here, we report evidence of a single founder event giving rise to the SCN1B c.363C>G variant in 14 independent families with epilepsy. A common haplotype was observed in all families, and the age of the most recent common ancestor was estimated to be approximately 800 years ago. Analysis of UK Biobank whole-exome-sequencing data identified 74 individuals with the same variant. All individuals carried haplotypes matching the epilepsy-affected families, suggesting all instances of the variant derive from a single mutational event. This unusual finding of a variant causing an autosomal dominant, early-onset disease in an outbred population that has persisted over many generations can be attributed to the relatively mild phenotype in most carriers and incomplete penetrance. Founder events are well established in autosomal recessive and late-onset disorders but are rarely observed in early-onset, autosomal dominant diseases. These findings suggest variants present in the population at low frequencies should be considered potentially pathogenic in mild phenotypes with incomplete penetrance and may be more important contributors to the genetic landscape than previously thought.

Iron homeostasis governs erythroid phenotype in polycythemia vera.

Polycythemia vera (PV) is a myeloproliferative neoplasm driven by activating mutations in JAK2 that result in unrestrained erythrocyte production, increasing patients' hematocrit and hemoglobin concentrations, placing them at risk of life-threatening thrombotic events. Our genome-wide association study of 440 PV cases and 403 351 controls using UK Biobank data showed that single nucleotide polymorphisms in HFE known to cause hemochromatosis are highly associated with PV diagnosis, linking iron regulation to PV. Analysis of the FinnGen dataset independently confirmed overrepresentation of homozygous HFE variants in patients with PV. HFE influences the expression of hepcidin, the master regulator of systemic iron homeostasis. Through genetic dissection of mouse models of PV, we show that the PV erythroid phenotype is directly linked to hepcidin expression: endogenous hepcidin upregulation alleviates erythroid disease whereas hepcidin ablation worsens it. Furthermore, we demonstrate that in PV, hepcidin is not regulated by expanded erythropoiesis but is likely governed by inflammatory cytokines signaling via GP130-coupled receptors. These findings have important implications for understanding the pathophysiology of PV and offer new therapeutic strategies for this disease.

Population-level genome-wide STR discovery and validation for population structure and genetic diversity assessment of Plasmodium species.

Short tandem repeats (STRs) are highly informative genetic markers that have been used extensively in population genetics analysis. They are an important source of genetic diversity and can also have functional impact. Despite the availability of bioinformatic methods that permit large-scale genome-wide genotyping of STRs from whole genome sequencing data, they have not previously been applied to sequencing data from large collections of malaria parasite field samples. Here, we have genotyped STRs using HipSTR in more than 3,000 Plasmodium falciparum and 174 Plasmodium vivax published whole-genome sequence data from samples collected across the globe. High levels of noise and variability in the resultant callset necessitated the development of a novel method for quality control of STR genotype calls. A set of high-quality STR loci (6,768 from P. falciparum and 3,496 from P. vivax) were used to study Plasmodium genetic diversity, population structures and genomic signatures of selection and these were compared to genome-wide single nucleotide polymorphism (SNP) genotyping data. In addition, the genome-wide information about genetic variation and other characteristics of STRs in P. falciparum and P. vivax have been available in an interactive web-based R Shiny application PlasmoSTR (https://github.com/bahlolab/PlasmoSTR).

Sporadic hypothalamic hamartoma is a ciliopathy with somatic and bi-allelic contributions.

Hypothalamic hamartoma with gelastic seizures is a well-established cause of drug-resistant epilepsy in early life. The development of novel surgical techniques has permitted the genomic interrogation of hypothalamic hamartoma tissue. This has revealed causative mosaic variants within GLI3, OFD1 and other key regulators of the sonic-hedgehog pathway in a minority of cases. Sonic-hedgehog signalling proteins localize to the cellular organelle primary cilia. We therefore explored the hypothesis that cilia gene variants may underlie hitherto unsolved cases of sporadic hypothalamic hamartoma. We performed high-depth exome sequencing and chromosomal microarray on surgically resected hypothalamic hamartoma tissue and paired leukocyte-derived DNA from 27 patients. We searched for both germline and somatic variants under both dominant and bi-allelic genetic models. In hamartoma-derived DNA of seven patients we identified bi-allelic (one germline, one somatic) variants within one of four cilia genes-DYNC2I1, DYNC2H1, IFT140 or SMO. In eight patients, we identified single somatic variants in the previously established hypothalamic hamartoma disease genes GLI3 or OFD1. Overall, we established a plausible molecular cause for 15/27 (56%) patients. Here, we expand the genetic architecture beyond single variants within dominant disease genes that cause sporadic hypothalamic hamartoma to bi-allelic (one germline/one somatic) variants, implicate three novel cilia genes and reconceptualize the disorder as a ciliopathy.

Progressive myoclonus epilepsies-Residual unsolved cases have marked genetic heterogeneity including dolichol-dependent protein glycosylation pathway genes.

Progressive myoclonus epilepsies (PMEs) comprise a group of clinically and genetically heterogeneous rare diseases. Over 70% of PME cases can now be molecularly solved. Known PME genes encode a variety of proteins, many involved in lysosomal and endosomal function. We performed whole-exome sequencing (WES) in 84 (78 unrelated) unsolved PME-affected individuals, with or without additional family members, to discover novel causes. We identified likely disease-causing variants in 24 out of 78 (31%) unrelated individuals, despite previous genetic analyses. The diagnostic yield was significantly higher for individuals studied as trios or families (14/28) versus singletons (10/50) (OR = 3.9, p value = 0.01, Fisher's exact test). The 24 likely solved cases of PME involved 18 genes. First, we found and functionally validated five heterozygous variants in NUS1 and DHDDS and a homozygous variant in ALG10, with no previous disease associations. All three genes are involved in dolichol-dependent protein glycosylation, a pathway not previously implicated in PME. Second, we independently validate SEMA6B as a dominant PME gene in two unrelated individuals. Third, in five families, we identified variants in established PME genes; three with intronic or copy-number changes (CLN6, GBA, NEU1) and two very rare causes (ASAH1, CERS1). Fourth, we found a group of genes usually associated with developmental and epileptic encephalopathies, but here, remarkably, presenting as PME, with or without prior developmental delay. Our systematic analysis of these cases suggests that the small residuum of unsolved cases will most likely be a collection of very rare, genetically heterogeneous etiologies.

Familial Mesial Temporal Lobe Epilepsy: Clinical Spectrum and Genetic Evidence for a Polygenic Architecture.

OBJECTIVE: Familial mesial temporal lobe epilepsy (FMTLE) is an important focal epilepsy syndrome; its molecular genetic basis is unknown. Clinical descriptions of FMTLE vary between a mild syndrome with prominent déjà vu to a more severe phenotype with febrile seizures and hippocampal sclerosis. We aimed to refine the phenotype of FMTLE by analyzing a large cohort of patients and asked whether common risk variants for focal epilepsy and/or febrile seizures, measured by polygenic risk scores (PRS), are enriched in individuals with FMTLE. METHODS: We studied 134 families with ≥ 2 first or second-degree relatives with temporal lobe epilepsy, with clear mesial ictal semiology required in at least one individual. PRS were calculated for 227 FMTLE cases, 124 unaffected relatives, and 16,077 population controls. RESULTS: The age of patients with FMTLE onset ranged from 2.5 to 70 years (median = 18, interquartile range = 13-28 years). The most common focal seizure symptom was déjà vu (62% of cases), followed by epigastric rising sensation (34%), and fear or anxiety (22%). The clinical spectrum included rare cases with drug-resistance and/or hippocampal sclerosis. FMTLE cases had a higher mean focal epilepsy PRS than population controls (odds ratio = 1.24, 95% confidence interval = 1.06, 1.46, p = 0.007); in contrast, no enrichment for the febrile seizure PRS was observed. INTERPRETATION: FMTLE is a generally mild drug-responsive syndrome with déjà vu being the commonest symptom. In contrast to dominant monogenic focal epilepsy syndromes, our molecular data support a polygenic basis for FMTLE. Furthermore, the PRS data suggest that sub-genome-wide significant focal epilepsy genome-wide association study single nucleotide polymorphisms are important risk variants for FMTLE. ANN NEUROL 2023;94:825-835.

Genetic aetiologies for childhood speech disorder: novel pathways co-expressed during brain development.

Childhood apraxia of speech (CAS), the prototypic severe childhood speech disorder, is characterized by motor programming and planning deficits. Genetic factors make substantive contributions to CAS aetiology, with a monogenic pathogenic variant identified in a third of cases, implicating around 20 single genes to date. Here we aimed to identify molecular causation in 70 unrelated probands ascertained with CAS. We performed trio genome sequencing. Our bioinformatic analysis examined single nucleotide, indel, copy number, structural and short tandem repeat variants. We prioritised appropriate variants arising de novo or inherited that were expected to be damaging based on in silico predictions. We identified high confidence variants in 18/70 (26%) probands, almost doubling the current number of candidate genes for CAS. Three of the 18 variants affected SETBP1, SETD1A and DDX3X, thus confirming their roles in CAS, while the remaining 15 occurred in genes not previously associated with this disorder. Fifteen variants arose de novo and three were inherited. We provide further novel insights into the biology of child speech disorder, highlighting the roles of chromatin organization and gene regulation in CAS, and confirm that genes involved in CAS are co-expressed during brain development. Our findings confirm a diagnostic yield comparable to, or even higher, than other neurodevelopmental disorders with substantial de novo variant burden. Data also support the increasingly recognised overlaps between genes conferring risk for a range of neurodevelopmental disorders. Understanding the aetiological basis of CAS is critical to end the diagnostic odyssey and ensure affected individuals are poised for precision medicine trials.

Stuttering associated with a pathogenic variant in the chaperone protein cyclophilin 40.

Stuttering is a common speech disorder that interrupts speech fluency and tends to cluster in families. Typically, stuttering is characterized by speech sounds, words or syllables which may be repeated or prolonged and speech that may be further interrupted by hesitations or 'blocks'. Rare variants in a small number of genes encoding lysosomal pathway proteins have been linked to stuttering. We studied a large four-generation family in which persistent stuttering was inherited in an autosomal dominant manner with disruption of the cortico-basal-ganglia-thalamo-cortical network found on imaging. Exome sequencing of three affected family members revealed the PPID c.808C>T (p.Pro270Ser) variant that segregated with stuttering in the family. We generated a Ppid p.Pro270Ser knock-in mouse model and performed ex vivo imaging to assess for brain changes. Diffusion-weighted MRI in the mouse revealed significant microstructural changes in the left corticospinal tract, as previously implicated in stuttering. Quantitative susceptibility mapping also detected changes in cortico-striatal-thalamo-cortical loop tissue composition, consistent with findings in affected family members. This is the first report to implicate a chaperone protein in the pathogenesis of stuttering. The humanized Ppid murine model recapitulates network findings observed in affected family members.

The Genetics of Epilepsy.

Epilepsy encompasses a group of heterogeneous brain diseases that affect more than 50 million people worldwide. Epilepsy may have discernible structural, infectious, metabolic, and immune etiologies; however, in most people with epilepsy, no obvious cause is identifiable. Based initially on family studies and later on advances in gene sequencing technologies and computational approaches, as well as the establishment of large collaborative initiatives, we now know that genetics plays a much greater role in epilepsy than was previously appreciated. Here, we review the progress in the field of epilepsy genetics and highlight molecular discoveries in the most important epilepsy groups, including those that have been long considered to have a nongenetic cause. We discuss where the field of epilepsy genetics is moving as it enters a new era in which the genetic architecture of common epilepsies is starting to be unraveled.

Germline and Mosaic Variants in PRKACA and PRKACB Cause a Multiple Congenital Malformation Syndrome.

PRKACA and PRKACB code for two catalytic subunits (Cα and Cß) of cAMP-dependent protein kinase (PKA), a pleiotropic holoenzyme that regulates numerous fundamental biological processes such as metabolism, development, memory, and immune response. We report seven unrelated individuals presenting with a multiple congenital malformation syndrome in whom we identified heterozygous germline or mosaic missense variants in PRKACA or PRKACB. Three affected individuals were found with the same PRKACA variant, and the other four had different PRKACB mutations. In most cases, the mutations arose de novo, and two individuals had offspring with the same condition. Nearly all affected individuals and their affected offspring shared an atrioventricular septal defect or a common atrium along with postaxial polydactyly. Additional features included skeletal abnormalities and ectodermal defects of variable severity in five individuals, cognitive deficit in two individuals, and various unusual tumors in one individual. We investigated the structural and functional consequences of the variants identified in PRKACA and PRKACB through the use of several computational and experimental approaches, and we found that they lead to PKA holoenzymes which are more sensitive to activation by cAMP than are the wild-type proteins. Furthermore, expression of PRKACA or PRKACB variants detected in the affected individuals inhibited hedgehog signaling in NIH 3T3 fibroblasts, thereby providing an underlying mechanism for the developmental defects observed in these cases. Our findings highlight the importance of both Cα and Cß subunits of PKA during human development.

A survey of RNA editing at single-cell resolution links interneurons to schizophrenia and autism.

Conversion of adenosine to inosine in RNA by ADAR enzymes, termed "RNA editing," is essential for healthy brain development. Editing is dysregulated in neuropsychiatric diseases, but has not yet been investigated at scale at the level of individual neurons. We quantified RNA editing sites in nuclear transcriptomes of 3055 neurons from six cortical regions of a neurotypical female donor, and found 41,930 sites present in at least ten nuclei. Most sites were located within Alu repeats in introns or 3' UTRs, and approximately 80% were cataloged in public RNA editing databases. We identified 9285 putative novel editing sites, 29% of which were also detectable in unrelated donors. Intersection with results from bulk RNA-seq studies provided cell-type and spatial context for 1730 sites that are differentially edited in schizophrenic brain donors, and 910 such sites in autistic donors. Autism-related genes were also enriched with editing sites predicted to modify RNA structure. Inhibitory neurons showed higher overall transcriptome editing than excitatory neurons, and the highest editing rates were observed in the frontal cortex. We used generalized linear models to identify differentially edited sites and genes between cell types. Twenty nine genes were preferentially edited in excitatory neurons, and 43 genes were edited more heavily in inhibitory neurons, including RBFOX1, its target genes, and genes in the autism-associated Prader-Willi locus (15q11). The abundance of SNORD115/116 genes from locus 15q11 was positively associated with editing activity across the transcriptome. We contend that insufficient editing of autism-related genes in inhibitory neurons may contribute to the specific perturbation of those cells in autism.

Heterozygous PNPT1 Variants Cause Spinocerebellar Ataxia Type 25.

OBJECTIVE: Dominant spinocerebellar ataxias (SCA) are characterized by genetic heterogeneity. Some mapped and named loci remain without a causal gene identified. Here we applied next generation sequencing (NGS) to uncover the genetic etiology of the SCA25 locus. METHODS: Whole-exome and whole-genome sequencing were performed in families linked to SCA25, including the French family in which the SCA25 locus was originally mapped. Whole exome sequence data were interrogated in a cohort of 796 ataxia patients of unknown etiology. RESULTS: The SCA25 phenotype spans a slowly evolving sensory and cerebellar ataxia, in most cases attributed to ganglionopathy. A pathogenic variant causing exon skipping was identified in the gene encoding Polyribonucleotide Nucleotidyltransferase PNPase 1 (PNPT1) located in the SCA25 linkage interval. A second splice variant in PNPT1 was detected in a large Australian family with a dominant ataxia also mapping to SCA25. An additional nonsense variant was detected in an unrelated individual with ataxia. Both nonsense and splice heterozygous variants result in premature stop codons, all located in the S1-domain of PNPase. In addition, an elevated type I interferon response was observed in blood from all affected heterozygous carriers tested. PNPase notably prevents the abnormal accumulation of double-stranded mtRNAs in the mitochondria and leakage into the cytoplasm, associated with triggering a type I interferon response. INTERPRETATION: This study identifies PNPT1 as a new SCA gene, responsible for SCA25, and highlights biological links between alterations of mtRNA trafficking, interferonopathies and ataxia. ANN NEUROL 2022;92:122-137.

Spatial distribution of metabolites in the retina and its relevance to studies of metabolic retinal disorders.

INTRODUCTION: The primate retina has evolved regional specialisations for specific visual functions. The macula is specialised towards high acuity vision and is an area that contains an increased density of cone photoreceptors and signal processing neurons. Different regions in the retina display unique susceptibility to pathology, with many retinal diseases primarily affecting the macula. OBJECTIVES: To better understand the properties of different retinal areas we studied the differential distribution of metabolites across the retina. METHODS: We conducted an untargeted metabolomics analysis on full-thickness punches from three different regions (macula, temporal peri-macula and periphery) of healthy primate retina. RESULTS: Nearly half of all metabolites identified showed differential abundance in at least one comparison between the three regions. Furthermore, mapping metabolomics results from macula-specific eye diseases onto our region-specific metabolite distributions revealed differential abundance defining systemic metabolic dysregulations that were region specific. CONCLUSIONS: The unique metabolic phenotype of different retinal regions is likely due to the differential distribution of different cell types in these regions reflecting the specific metabolic requirements of each cell type. Our results may help to better understand the pathobiology of retinal diseases with region specificity.

Genes4Epilepsy: An epilepsy gene resource.

OBJECTIVE: "How many epilepsy genes are there?" is a frequently asked question. We sought to (1) provide a curated list of genes that cause monogenic epilepsies, and (2) compare and contrast epilepsy gene panels from multiple sources. METHODS: We compared genes included on the epilepsy panels (as of July 29, 2022) of four clinical diagnostic providers: Invitae, GeneDx, Fulgent Genetics, and Blueprint Genetics; and two research resources: PanelApp Australia and ClinGen. A master list of all unique genes was supplemented by additional genes identified via PubMed searches up until August 15, 2022, using the search terms "genetics" AND/OR "epilepsy" AND/OR "seizures". Evidence supporting a monogenic role for all genes was manually reviewed; those with limited or disputed evidence were excluded. All genes were annotated according to inheritance pattern and broad epilepsy phenotype. RESULTS: The comparison of genes included on epilepsy clinical panels revealed high heterogeneity in both number of genes (range: 144-511) and content. Just 111 genes (15.5%) were included on all four clinical panels. Subsequent manual curation of all "epilepsy genes" identified >900 monogenic etiologies. Almost 90% of genes were associated with developmental and epileptic encephalopathies. By comparison only 5% of genes were associated with monogenic causes of "common epilepsies" (i.e., generalized and focal epilepsy syndromes). Autosomal recessive genes were most frequent (56% of genes); however, this varied according to the associated epilepsy phenotype(s). Genes associated with common epilepsy syndromes were more likely to be dominantly inherited and associated with multiple epilepsy types. SIGNIFICANCE: Our curated list of monogenic epilepsy genes is publicly available: github.com/bahlolab/genes4epilepsy and will be regularly updated. This gene resource can be utilized to target genes beyond those included on clinical gene panels, for gene enrichment methods and candidate gene prioritization. We invite ongoing feedback and contributions from the scientific community via genes4-epilepsy@unimelb.edu.au.

WWOX developmental and epileptic encephalopathy: Understanding the epileptology and the mortality risk.

OBJECTIVE: WWOX is an autosomal recessive cause of early infantile developmental and epileptic encephalopathy (WWOX-DEE), also known as WOREE (WWOX-related epileptic encephalopathy). We analyzed the epileptology and imaging features of WWOX-DEE, and investigated genotype-phenotype correlations, particularly with regard to survival. METHODS: We studied 13 patients from 12 families with WWOX-DEE. Information regarding seizure semiology, comorbidities, facial dysmorphisms, and disease outcome were collected. Electroencephalographic (EEG) and brain magnetic resonance imaging (MRI) data were analyzed. Pathogenic WWOX variants from our cohort and the literature were coded as either null or missense, allowing individuals to be classified into one of three genotype classes: (1) null/null, (2) null/missense, (3) missense/missense. Differences in survival outcome were estimated using the Kaplan-Meier method. RESULTS: All patients experienced multiple seizure types (median onset = 5 weeks, range = 1 day-10 months), the most frequent being focal (85%), epileptic spasms (77%), and tonic seizures (69%). Ictal EEG recordings in six of 13 patients showed tonic (n = 5), myoclonic (n = 2), epileptic spasms (n = 2), focal (n = 1), and migrating focal (n = 1) seizures. Interictal EEGs demonstrated slow background activity with multifocal discharges, predominantly over frontal or temporo-occipital regions. Eleven of 13 patients had a movement disorder, most frequently dystonia. Brain MRIs revealed severe frontotemporal, hippocampal, and optic atrophy, thin corpus callosum, and white matter signal abnormalities. Pathogenic variants were located throughout WWOX and comprised both missense and null changes including five copy number variants (four deletions, one duplication). Survival analyses showed that patients with two null variants are at higher mortality risk (p-value = .0085, log-rank test). SIGNIFICANCE: Biallelic WWOX pathogenic variants cause an early infantile developmental and epileptic encephalopathy syndrome. The most common seizure types are focal seizures and epileptic spasms. Mortality risk is associated with mutation type; patients with biallelic null WWOX pathogenic variants have significantly lower survival probability compared to those carrying at least one presumed hypomorphic missense pathogenic variant.

Global diversity and balancing selection of 23 leading Plasmodium falciparum candidate vaccine antigens.

Investigation of the diversity of malaria parasite antigens can help prioritize and validate them as vaccine candidates and identify the most common variants for inclusion in vaccine formulations. Studies of vaccine candidates of the most virulent human malaria parasite, Plasmodium falciparum, have focused on a handful of well-known antigens, while several others have never been studied. Here we examine the global diversity and population structure of leading vaccine candidate antigens of P. falciparum using the MalariaGEN Pf3K (version 5.1) resource, comprising more than 2600 genomes from 15 malaria endemic countries. A stringent variant calling pipeline was used to extract high quality antigen gene 'haplotypes' from the global dataset and a new R-package named VaxPack was used to streamline population genetic analyses. In addition, a newly developed algorithm that enables spatial averaging of selection pressure on 3D protein structures was applied to the dataset. We analysed the genes encoding 23 leading and novel candidate malaria vaccine antigens including csp, trap, eba175, ama1, rh5, and CelTOS. Our analysis shows that current malaria vaccine formulations are based on rare haplotypes and thus may have limited efficacy against natural parasite populations. High levels of diversity with evidence of balancing selection was detected for most of the erythrocytic and pre-erythrocytic antigens. Measures of natural selection were then mapped to 3D protein structures to predict targets of functional antibodies. For some antigens, geographical variation in the intensity and distribution of these signals on the 3D structure suggests adaptation to different human host or mosquito vector populations. This study provides an essential framework for the diversity of P. falciparum antigens to be considered in the design of the next generation of malaria vaccines.

Atypical development of Broca's area in a large family with inherited stuttering.

Developmental stuttering is a condition of speech dysfluency, characterized by pauses, blocks, prolongations and sound or syllable repetitions. It affects around 1% of the population, with potential detrimental effects on mental health and long-term employment. Accumulating evidence points to a genetic aetiology, yet gene-brain associations remain poorly understood due to a lack of MRI studies in affected families. Here we report the first neuroimaging study of developmental stuttering in a family with autosomal dominant inheritance of persistent stuttering. We studied a four-generation family, 16 family members were included in genotyping analysis. T1-weighted and diffusion-weighted MRI scans were conducted on seven family members (six male; aged 9-63 years) with two age and sex matched controls without stuttering (n = 14). Using Freesurfer, we analysed cortical morphology (cortical thickness, surface area and local gyrification index) and basal ganglia volumes. White matter integrity in key speech and language tracts (i.e. frontal aslant tract and arcuate fasciculus) was also analysed using MRtrix and probabilistic tractography. We identified a significant age by group interaction effect for cortical thickness in the left hemisphere pars opercularis (Broca's area). In affected family members this region failed to follow the typical trajectory of age-related thinning observed in controls. Surface area analysis revealed the middle frontal gyrus region was reduced bilaterally in the family (all cortical morphometry significance levels set at a vertex-wise threshold of P < 0.01, corrected for multiple comparisons). Both the left and right globus pallidus were larger in the family than in the control group (left P = 0.017; right P = 0.037), and a larger right globus pallidus was associated with more severe stuttering (rho = 0.86, P = 0.01). No white matter differences were identified. Genotyping identified novel loci on chromosomes 1 and 4 that map with the stuttering phenotype. Our findings denote disruption within the cortico-basal ganglia-thalamo-cortical network. The lack of typical development of these structures reflects the anatomical basis of the abnormal inhibitory control network between Broca's area and the striatum underpinning stuttering in these individuals. This is the first evidence of a neural phenotype in a family with an autosomal dominantly inherited stuttering.