Diversity and phylogenetic relationships of haemosporidian and hemogregarine parasites in Australian lizards.

Apicomplexa is a large monophyletic phylum of unicellular, parasitic organisms. Reptiles are hosts to both haemosporidian (Haemosporida) and hemogregarine (Eucoccidiorida) apicomplexan blood parasites. Within reptiles our understanding of their diversity remains limited, with a paucity of information from Australia, despite a high diversity of squamates (snakes and lizards). We provide a preliminary assessment of haemosporidian and hemogregarine diversity occurring in lizards across northern tropical Australia, building on existing data with results from a microscopy and genetic assessment. We screened total of 233 blood slides using microscopy and detected hemogregarines in 25 geckos, 2 skinks and 1 agamid, while haemosporidians were detected in 13 geckos. DNA sequencing of 28 samples of the hemogregarine 18S rRNA (∼900 bp) nuclear gene revealed five lineages of Australian lizard hemogregarines within heteroxenous adeleids. We sequenced 10 samples of Haemosporida mtDNA (cytb & coI: ∼1313 bp) and phylogenetic analysis with 30 previously published sequences revealed that the Australian Haemosporida grouped within the Haemoproteidae but were not supported as a monophyletic clade. Our results demonstrate that there is significant undocumented evolutionary diversity in Australian lizard haemosporidian and hemogregarine parasites, with preliminary evidence of significantly higher infection rates in geckos.

Neuronal ceroid lipofuscinosis type 2: an Australian case series.

AIM: Late infantile neuronal ceroid lipofuscinosis type 2 (CLN2) disease is a rare neurodegenerative disorder presenting in children aged 2-4 years with seizures and loss of motor and language skills, followed by blindness and death in late childhood. Initial presenting features are similar to a range of common epilepsies. We aim to highlight typical clinical and radiological features that may prompt diagnosis of CLN2 disease in early disease stages. METHODS: We present a series of 13 Australian patients with CLN2 disease, describing clinical features, disease evolution, neuroimaging, electroencephalogram, biochemical and genetic results. Expert neuroradiological magnetic resonance imaging (MRI) analysis was retrospectively performed on 10 cases. RESULTS: Twelve patients presented with seizures, with initial seizures being focal (n = 4), generalised tonic-clonic (n = 3), absence (n = 3) and febrile (n = 2). Eleven patients (85%) had a language delay before the onset of seizures. Cerebellar or cerebral atrophy was noted in all patients on centralised MRI review, with abnormalities of the brain-stem, ventricles, corpus callosum and hippocampi. CONCLUSIONS: Early language delay with the onset of seizures at 2-4 years of age is the hallmark of CLN2 disease. MRI findings of early subtle atrophy in the cerebellum or posterior cortical regions should hasten testing for CLN2 disease to enable early initiation of enzyme replacement therapy.

BRAT1 encephalopathy: a recessive cause of epilepsy of infancy with migrating focal seizures.

Epilepsy of infancy with migrating focal seizures (EIMFS), one of the most severe developmental and epileptic encephalopathy syndromes, is characterized by seizures that migrate from one hemisphere to the other. EIMFS is genetically heterogeneous with 33 genes. We report five patients with EIMFS caused by recessive BRAT1 variants, identified via next generation sequencing. Recessive pathogenic variants in BRAT1 cause the rigidity and multifocal seizure syndrome, lethal neonatal with hypertonia, microcephaly, and intractable multifocal seizures. The epileptology of BRAT1 encephalopathy has not been well described. All five patients were profoundly impaired with seizure onset in the first week of life and focal seizure migration between hemispheres. We show that BRAT1 is an important recessive cause of EIMFS with onset in the first week of life, profound impairment, and early death. Early recognition of this genetic aetiology will inform management and reproductive counselling.

The Genetic Landscape of Epilepsy of Infancy with Migrating Focal Seizures.

OBJECTIVE: Epilepsy of infancy with migrating focal seizures (EIMFS) is one of the most severe developmental and epileptic encephalopathies. We delineate the genetic causes and genotype-phenotype correlations of a large EIMFS cohort. METHODS: Phenotypic and molecular data were analyzed on patients recruited through an international collaborative study. RESULTS: We ascertained 135 patients from 128 unrelated families. Ninety-three of 135 (69%) had causative variants (42/55 previously reported) across 23 genes, including 9 novel EIMFS genes: de novo dominant GABRA1, GABRB1, ATP1A3; X-linked CDKL5, PIGA; and recessive ITPA, AIMP1, KARS, WWOX. The most frequently implicated genes were KCNT1 (36/135, 27%) and SCN2A (10/135, 7%). Mosaicism occurred in 2 probands (SCN2A, GABRB3) and 3 unaffected mothers (KCNT1). Median age at seizure onset was 4 weeks, with earlier onset in the SCN2A, KCNQ2, and BRAT1 groups. Epileptic spasms occurred in 22% patients. A total of 127 patients had severe to profound developmental impairment. All but 7 patients had ongoing seizures. Additional features included microcephaly, movement disorders, spasticity, and scoliosis. Mortality occurred in 33% at median age 2 years 7 months. INTERPRETATION: We identified a genetic cause in 69% of patients with EIMFS. We highlight the genetic heterogeneity of EIMFS with 9 newly implicated genes, bringing the total number to 33. Mosaicism was observed in probands and parents, carrying critical implications for recurrence risk. EIMFS pathophysiology involves diverse molecular processes from gene and protein regulation to ion channel function and solute trafficking. ANN NEUROL 2019;86:821-831.

AMPA receptor GluA2 subunit defects are a cause of neurodevelopmental disorders.

AMPA receptors (AMPARs) are tetrameric ligand-gated channels made up of combinations of GluA1-4 subunits encoded by GRIA1-4 genes. GluA2 has an especially important role because, following post-transcriptional editing at the Q607 site, it renders heteromultimeric AMPARs Ca2+-impermeable, with a linear relationship between current and trans-membrane voltage. Here, we report heterozygous de novo GRIA2 mutations in 28 unrelated patients with intellectual disability (ID) and neurodevelopmental abnormalities including autism spectrum disorder (ASD), Rett syndrome-like features, and seizures or developmental epileptic encephalopathy (DEE). In functional expression studies, mutations lead to a decrease in agonist-evoked current mediated by mutant subunits compared to wild-type channels. When GluA2 subunits are co-expressed with GluA1, most GRIA2 mutations cause a decreased current amplitude and some also affect voltage rectification. Our results show that de-novo variants in GRIA2 can cause neurodevelopmental disorders, complementing evidence that other genetic causes of ID, ASD and DEE also disrupt glutamatergic synaptic transmission.

De Novo Pathogenic Variants in CACNA1E Cause Developmental and Epileptic Encephalopathy with Contractures, Macrocephaly, and Dyskinesias.

Developmental and epileptic encephalopathies (DEEs) are severe neurodevelopmental disorders often beginning in infancy or early childhood that are characterized by intractable seizures, abundant epileptiform activity on EEG, and developmental impairment or regression. CACNA1E is highly expressed in the central nervous system and encodes the α1-subunit of the voltage-gated CaV2.3 channel, which conducts high voltage-activated R-type calcium currents that initiate synaptic transmission. Using next-generation sequencing techniques, we identified de novo CACNA1E variants in 30 individuals with DEE, characterized by refractory infantile-onset seizures, severe hypotonia, and profound developmental impairment, often with congenital contractures, macrocephaly, hyperkinetic movement disorders, and early death. Most of the 14, partially recurring, variants cluster within the cytoplasmic ends of all four S6 segments, which form the presumed CaV2.3 channel activation gate. Functional analysis of several S6 variants revealed consistent gain-of-function effects comprising facilitated voltage-dependent activation and slowed inactivation. Another variant located in the domain II S4-S5 linker results in facilitated activation and increased current density. Five participants achieved seizure freedom on the anti-epileptic drug topiramate, which blocks R-type calcium channels. We establish pathogenic variants in CACNA1E as a cause of DEEs and suggest facilitated R-type calcium currents as a disease mechanism for human epilepsy and developmental disorders.

Identifying hybridization and admixture using SNPs: application of the DArTseq platform in phylogeographic research on vertebrates.

Next-generation sequencing (NGS) approaches are increasingly being used to generate multi-locus data for phylogeographic and evolutionary genetics research. We detail the applicability of a restriction enzyme-mediated genome complexity reduction approach with subsequent NGS (DArTseq) in vertebrate study systems at different evolutionary and geographical scales. We present two case studies using SNP data from the DArTseq molecular marker platform. First, we used DArTseq in a large phylogeographic study of the agamid lizard Ctenophorus caudicinctus, including 91 individuals and spanning the geographical range of this species across arid Australia. A low-density DArTseq assay resulted in 28 960 SNPs, with low density referring to a comparably reduced set of identified and sequenced markers as a cost-effective approach. Second, we applied this approach to an evolutionary genetics study of a classic frog hybrid zone (Litoria ewingii-Litoria paraewingi) across 93 individuals, which resulted in 48 117 and 67 060 SNPs for a low- and high-density assay, respectively. We provide a docker-based workflow to facilitate data preparation and analysis, then analyse SNP data using multiple methods including Bayesian model-based clustering and conditional likelihood approaches. Based on comparison of results from the DArTseq platform and traditional molecular approaches, we conclude that DArTseq can be used successfully in vertebrates and will be of particular interest to researchers working at the interface between population genetics and phylogenetics, exploring species boundaries, gene exchange and hybridization.

Expression of a hindlimb-determining factor Pitx1 in the forelimb of the lizard Pogona vitticeps during morphogenesis.

With over 9000 species, squamates, which include lizards and snakes, are the largest group of reptiles and second-largest order of vertebrates, spanning a vast array of appendicular skeletal morphology. As such, they provide a promising system for examining developmental and molecular processes underlying limb morphology. Using the central bearded dragon (Pogona vitticeps) as the primary study model, we examined limb morphometry throughout embryonic development and characterized the expression of three known developmental genes (GHR, Pitx1 and Shh) from early embryonic stage through to hatchling stage via reverse transcription quantitative polymerase chain reaction (RT-qPCR) and immunohistochemistry (IHC). In this study, all genes were found to be transcribed in both the forelimbs and hindlimbs of P. vitticeps. While the highest level of GHR expression occurred at the hatchling stage, Pitx1 and Shh expression was greatest earlier during embryogenesis, which coincides with the onset of the differentiation between forelimb and hindlimb length. We compared our finding of Pitx1 expression-a hindlimb-determining gene-in the forelimbs of P. vitticeps to that in a closely related Australian agamid lizard, Ctenophorus pictus, where we found Pitx1 expression to be more highly expressed in the hindlimb compared with the forelimb during early and late morphogenesis-a result consistent with that found across other tetrapods. Expression of Pitx1 in forelimbs has only rarely been documented, including via in situ hybridization in a chicken and a frog. Our findings from both RT-qPCR and IHC indicate that further research across a wider range of tetrapods is needed to more fully understand evolutionary variation in molecular processes underlying limb morphology.

Inhibitor of cysteine proteases is critical for motility and infectivity of Plasmodium sporozoites.

UNLABELLED: Malaria is transmitted when motile sporozoites are injected into the dermis by an infected female Anopheles mosquito. Inside the mosquito vector, sporozoites egress from midgut-associated oocysts and eventually penetrate the acinar cells of salivary glands. Parasite-encoded factors with exclusive vital roles in the insect vector can be studied by classical reverse genetics. Here, we characterized the in vivo roles of Plasmodium berghei falstatin/ICP (inhibitor of cysteine proteases). This protein was previously suggested to act as a protease inhibitor during erythrocyte invasion. We show by targeted gene disruption that loss of ICP function does not affect growth inside the mammalian host but causes a complete defect in sporozoite transmission. Sporogony occurred normally in icp(-) parasites, but hemocoel sporozoites showed a defect in continuous gliding motility and infectivity for salivary glands, which are prerequisites for sporozoite transmission to the mammalian host. Absence of ICP correlates with enhanced cleavage of circumsporozoite protein, in agreement with a role as a protease regulator. We conclude that ICP is essential for only the final stages of sporozoite maturation inside the mosquito vector. This study is the first genetic evidence that an ICP is necessary for the productive motility of a eukaryotic parasitic cell. IMPORTANCE: Cysteine proteases and their inhibitors are considered ideal drug targets for the treatment of a wide range of diseases, including cancer and parasitic infections. In protozoan parasites, including Leishmania, Trypanosoma, and Plasmodium, cysteine proteases play important roles in life cycle progression. A mouse malaria model provides an unprecedented opportunity to study the roles of a parasite-encoded inhibitor of cysteine proteases (ICP) over the entire parasite life cycle. By precise gene deletion, we found no evidence that ICP influences disease progression or parasite virulence. Instead, we discovered that this factor is necessary for parasite movement and malaria transmission from mosquitoes to mammals. This finding in a fast-moving unicellular protozoan has important implications for malaria intervention strategies and the roles of ICPs in the regulation of eukaryotic cell migration.

Arrested oocyst maturation in Plasmodium parasites lacking type II NADH:ubiquinone dehydrogenase.

The Plasmodium mitochondrial electron transport chain has received considerable attention as a potential target for new antimalarial drugs. Atovaquone, a potent inhibitor of Plasmodium cytochrome bc(1), in combination with proguanil is recommended for chemoprophylaxis and treatment of malaria. The type II NADH:ubiquinone oxidoreductase (NDH2) is considered an attractive drug target, as its inhibition is thought to lead to the arrest of the mitochondrial electron transport chain and, as a consequence, pyrimidine biosynthesis, an essential pathway for the parasite. Using the rodent malaria parasite Plasmodium berghei as an in vivo infection model, we studied the role of NDH2 during Plasmodium life cycle progression. NDH2 can be deleted by targeted gene disruption and, thus, is dispensable for the pathogenic asexual blood stages, disproving the candidacy for an anti-malarial drug target. After transmission to the insect vector, NDH2-deficient ookinetes display an intact mitochondrial membrane potential. However, ndh2(-) parasites fail to develop into mature oocysts in the mosquito midgut. We propose that Plasmodium blood stage parasites rely on glycolysis as the main ATP generating process, whereas in the invertebrate vector, a glucose-deprived environment, the malaria parasite is dependent on an intact mitochondrial respiratory chain.

Expression of P. falciparum var genes involves exchange of the histone variant H2A.Z at the promoter.

Plasmodium falciparum employs antigenic variation to evade the human immune response by switching the expression of different variant surface antigens encoded by the var gene family. Epigenetic mechanisms including histone modifications and sub-nuclear compartmentalization contribute to transcriptional regulation in the malaria parasite, in particular to control antigenic variation. Another mechanism of epigenetic control is the exchange of canonical histones with alternative variants to generate functionally specialized chromatin domains. Here we demonstrate that the alternative histone PfH2A.Z is associated with the epigenetic regulation of var genes. In many eukaryotic organisms the histone variant H2A.Z mediates an open chromatin structure at promoters and facilitates diverse levels of regulation, including transcriptional activation. Throughout the asexual, intraerythrocytic lifecycle of P. falciparum we found that the P. falciparum ortholog of H2A.Z (PfH2A.Z) colocalizes with histone modifications that are characteristic of transcriptionally-permissive euchromatin, but not with markers of heterochromatin. Consistent with this finding, antibodies to PfH2A.Z co-precipitate the permissive modification H3K4me3. By chromatin-immunoprecipitation we show that PfH2A.Z is enriched in nucleosomes around the transcription start site (TSS) in both transcriptionally active and silent stage-specific genes. In var genes, however, PfH2A.Z is enriched at the TSS only during active transcription in ring stage parasites. Thus, in contrast to other genes, temporal var gene regulation involves histone variant exchange at promoter nucleosomes. Sir2 histone deacetylases are important for var gene silencing and their yeast ortholog antagonises H2A.Z function in subtelomeric yeast genes. In immature P. falciparum parasites lacking Sir2A or Sir2B high var transcription levels correlate with enrichment of PfH2A.Z at the TSS. As Sir2A knock out parasites mature the var genes are silenced, but PfH2A.Z remains enriched at the TSS of var genes; in contrast, PfH2A.Z is lost from the TSS of de-repressed var genes in mature Sir2B knock out parasites. This result indicates that PfH2A.Z occupancy at the active var promoter is antagonized by PfSir2A during the intraerythrocytic life cycle. We conclude that PfH2A.Z contributes to the nucleosome architecture at promoters and is regulated dynamically in active var genes.

Transcribed var genes associated with placental malaria in Malawian women.

Determining the diversity of PfEMP1 sequences expressed by Plasmodium falciparum-infected erythrocytes isolated from placentas is important for attempts to develop a pregnancy-specific malaria vaccine. The DBLgamma and var2csa DBL3x domains of PfEMP1 molecules are believed to mediate placental sequestration of infected erythrocytes, so the sequences encoding these domains were amplified from the cDNAs of placental parasites by using degenerate oligonucleotides. The levels of specific var cDNAs were then determined by quantitative reverse transcription-PCR. Homologues of var2csa DBL3x were the predominant sequences amplified from the cDNAs of most placental but not most children's parasites. There was 56% identity between all placental var2csa sequences. Many different DBLgamma domains were amplified from the cDNAs of placental and children's isolates. var2csa transcripts were the most abundant var transcripts of those tested in 11 of 12 placental isolates and 1 of 6 children's isolates. Gravidity did not affect the levels of var2csa transcripts. We concluded that placental malaria is frequently associated with transcription of var2csa but that other var genes are also expressed, and parasites expressing high levels of var2csa are not restricted to pregnant women. The diversity of var2csa sequences may be important for understanding immunity and for the development of vaccines for malaria during pregnancy.