Towards a molecular characterization of autism spectrum disorders: an exome sequencing and systems approach.

The hypothetical 'AXAS' gene network model that profiles functional patterns of heterogeneous DNA variants overrepresented in autism spectrum disorder (ASD), X-linked intellectual disability, attention deficit and hyperactivity disorder and schizophrenia was used in this current study to analyze whole exome sequencing data from an Australian ASD cohort. An optimized DNA variant filtering pipeline was used to identify loss-of-function DNA variations. Inherited variants from parents with a broader autism phenotype and de novo variants were found to be significantly associated with ASD. Gene ontology analysis revealed that putative rare causal variants cluster in key neurobiological processes and are overrepresented in functions involving neuronal development, signal transduction and synapse development including the neurexin trans-synaptic complex. We also show how a complex gene network model can be used to fine map combinations of inherited and de novo variations in families with ASD that converge in the L1CAM pathway. Our results provide an important step forward in the molecular characterization of ASD with potential for developing a tool to analyze the pathogenesis of individual affected families.

Neurodevelopmental and neuropsychiatric disorders represent an interconnected molecular system.

Many putative genetic factors that confer risk to neurodevelopmental disorders such as autism spectrum disorders (ASDs) and X-linked intellectual disability (XLID), and to neuropsychiatric disorders including attention deficit hyperactivity disorder (ADHD) and schizophrenia (SZ) have been identified in individuals from diverse human populations. Although there is significant aetiological heterogeneity within and between these conditions, recent data show that genetic factors contribute to their comorbidity. Many studies have identified candidate gene associations for these mental health disorders, albeit this is often done in a piecemeal fashion with little regard to the inherent molecular complexity. Here, we sought to abstract relationships from our knowledge of systems level biology to help understand the unique and common genetic drivers of these conditions. We undertook a global and systematic approach to build and integrate available data in gene networks associated with ASDs, XLID, ADHD and SZ. Complex network concepts and computational methods were used to investigate whether candidate genes associated with these conditions were related through mechanisms of gene regulation, functional protein-protein interactions, transcription factor (TF) and microRNA (miRNA) binding sites. Although our analyses show that genetic variations associated with the four disorders can occur in the same molecular pathways and functional domains, including synaptic transmission, there are patterns of variation that define significant differences between disorders. Of particular interest is DNA variations located in intergenic regions that comprise regulatory sites for TFs or miRNA. Our approach provides a hypothetical framework, which will help discovery and analysis of candidate genes associated with neurodevelopmental and neuropsychiatric disorders.

Big ideas for small brains: what can psychiatry learn from worms, flies, bees and fish?

While the research community has accepted the value of rodent models as informative research platforms, there is less awareness of the utility of other small vertebrate and invertebrate animal models. Neuroscience is increasingly turning to smaller, non-rodent models to understand mechanisms related to neuropsychiatric disorders. Although they can never replace clinical research, there is much to be learnt from 'small brains'. In particular, these species can offer flexible genetic 'tool kits' that can be used to explore the expression and function of candidate genes in different brain regions. Very small animals also offer efficiencies with respect to high-throughput screening programs. This review provides a concise overview of the utility of models based on worm, fruit fly, honeybee and zebrafish. Although these species may have small brains, they offer the neuropsychiatric research community opportunities to explore some of the most important research questions in our field.

Metabolic enzymes associated with xenobiotic and chemosensory responses in Nasonia vitripennis.

The numbers of glutathione S-transferase, cytochrome P450 and esterase genes in the genome of the hymenopteran parasitoid Nasonia vitripennis are about twice those found in the genome of another hymenopteran, the honeybee Apis mellifera. Some of the difference is associated with clades of these families implicated in xenobiotic resistance in other insects and some is in clades implicated in hormone and pheromone metabolism. The data support the hypothesis that the eusocial behaviour of the honeybee and the concomitant homeostasis of the nest environment may obviate the need for as many gene/enzyme systems associated with xenobiotic metabolism as are found in other species, including N. vitripennis, that are thought to encounter a wider range of potentially toxic xenobiotics in their diet and habitat.

Jhe in Gryllus assimilis: cloning, sequence-activity associations and phylogeny.

The 458 amino acid sequence of a mature JHE protein from the cricket Gryllus assimilis was identified after isolating the partial cDNA sequence encoding this protein from a fat body and midgut cDNA library. This hemimetabolan JHE sequence shows over 40% amino acid similarity to the known JHE sequences of several holometabolous insects. It also includes previously determined peptide sequences for G. assimilis JHE as well as two other motifs associated with JHE enzymes in holometabolous insects. The predicted molecular weight of the protein agrees with that of the JHE previously purified from G. assimilis. Partial genomic sequence encoding the Jhe contains two large (1330 and 2918bp) introns. No coding DNA sequence variation was observed over a 1293bp region between selected lines differing six to eight-fold in hemolymph JHE activity. However, a 19bp indel was found in one of the introns; the insertion was strongly associated with elevated hemolymph activity, both in the selected lines and in the F(2) progeny of crosses between them. Phylogenetic analyses localised the G. assimilis JHE to a clade containing dipteran and coleopteran JHEs, with lepidopteran JHEs occurring in a separate clade.

A deficit of detoxification enzymes: pesticide sensitivity and environmental response in the honeybee.

The honeybee genome has substantially fewer protein coding genes ( approximately 11 000 genes) than Drosophila melanogaster ( approximately 13 500) and Anopheles gambiae ( approximately 14 000). Some of the most marked differences occur in three superfamilies encoding xenobiotic detoxifying enzymes. Specifically there are only about half as many glutathione-S-transferases (GSTs), cytochrome P450 monooxygenases (P450s) and carboxyl/cholinesterases (CCEs) in the honeybee. This includes 10-fold or greater shortfalls in the numbers of Delta and Epsilon GSTs and CYP4 P450s, members of which clades have been recurrently associated with insecticide resistance in other species. These shortfalls may contribute to the sensitivity of the honeybee to insecticides. On the other hand there are some recent radiations in CYP6, CYP9 and certain CCE clades in A. mellifera that could be associated with the evolution of the hormonal and chemosensory processes underpinning its highly organized eusociality.

The flightless I protein and the gelsolin family in nuclear hormone receptor-mediated signalling.

The Drosophila melanogaster flightless I protein and its homologues in higher eukaryotes (FliI) are conserved members of the gelsolin family of actin-binding proteins. Members of the gelsolin family generally contain three or six copies of a 125-amino-acid residue gelsolin-related repeating unit, and may contain additional domains including the C-terminal villin-related 'headpiece' or N-terminal extensions such as the leucine-rich repeat of the FliI protein. Numerous studies including work done with mouse knockouts for gelsolin, villin and CapG support a role for the family in cytoskeletal actin dynamics. In both fruitfly and mouse, the FliI protein is also essential for early development. Recent studies indicate that supervillin, gelsolin and FliI are involved in intracellular signalling via nuclear hormone receptors including the androgen, oestrogen and thyroid hormone receptors. This unexpected role in signalling has opened a new area in research on the gelsolin family and is providing important new insights into the mechanisms of gene regulation via nuclear receptors.

Knockout of the rodent malaria parasite chitinase pbCHT1 reduces infectivity to mosquitoes.

During mosquito transmission, malaria ookinetes must cross a chitin-containing structure known as the peritrophic matrix (PM), which surrounds the infected blood meal in the mosquito midgut. In turn, ookinetes produce multiple chitinase activities presumably aimed at disrupting this physical barrier to allow ookinete invasion of the midgut epithelium. Plasmodium chitinase activities are demonstrated targets for human and avian malaria transmission blockade with the chitinase inhibitor allosamidin. Here, we identify and characterize the first chitinase gene of a rodent malaria parasite, Plasmodium berghei. We show that the gene, named PbCHT1, is a structural ortholog of PgCHT1 of the avian malaria parasite Plasmodium gallinaceum and a paralog of PfCHT1 of the human malaria parasite Plasmodium falciparum. Targeted disruption of PbCHT1 reduced parasite infectivity in Anopheles stephensi mosquitoes by up to 90%. Reductions in infectivity were also observed in ookinete feeds-an artificial situation where midgut invasion occurs before PM formation-suggesting that PbCHT1 plays a role other than PM disruption. PbCHT1 null mutants had no residual ookinete-derived chitinase activity in vitro, suggesting that P. berghei ookinetes express only one chitinase gene. Moreover, PbCHT1 activity appeared insensitive to allosamidin inhibition, an observation that raises questions about the use of allosamidin and components like it as potential malaria transmission-blocking drugs. Taken together, these findings suggest a fundamental divergence among rodent, avian, and human malaria parasite chitinases, with implications for the evolution of Plasmodium-mosquito interactions.

Identification of a juvenile hormone esterase gene by matching its peptide mass fingerprint with a sequence from the Drosophila genome project.

Juvenile hormone esterase (JHE, EC 3.1.1.1) from whole Drosophila melanogaster prepupae has previously been purified by selective precipitations, isoelectric focussing and two column chromatography steps. JHE bands from dried silver-stained SDS-PAGE gels of that material were digested with trypsin. The masses of the tryptic digest peptides were determined by MALDI-TOF mass spectrometry. Only one predicted gene product (CG8425) from the D. melanogaster genome matches the JHE tryptic fingerprint with high confidence. This predicted JHE sequence includes features that are conserved among all active members of the serine carboxylesterase multigene family as well as features peculiar to JHEs from other species. Also we show that this JHE can be purified by an alternative method using anion exchange chromotography followed by trifluoromethylketone affinity chromatography. A cDNA encoding this JHE was isolated using 3' and 5' RACE. This sequence is in agreement with the Drosophila genome project's prediction except that the sixth predicted intron is not removed; instead there is a stop codon followed by a polyadenylation signal and a polyA tail.

Overproduction of a P450 that metabolizes diazinon is linked to a loss-of-function in the chromosome 2 ali-esterase (MdalphaE7) gene in resistant house flies.

Up-regulation of detoxifying enzymes in insecticide-resistant strains of the house fly is a common mechanism for metabolic resistance. However, the molecular basis of this increased insecticide metabolism is not well understood. In the multiresistant Rutgers strain, several cytochromes P450 and glutathione S-transferases are constitutively overexpressed at the transcriptional level. Overexpression is the result of trans-regulation, and a regulatory gene has been located on chromosome 2. A Gly137 to Asp point mutation in alphaE7 esterase gene, leading to the loss of carboxylesterase activity, has been associated with organophosphate resistance in the house fly and the sheep blowfly. We show here that purified recombinant CYP6A1 is able to detoxify diazinon with a high efficiency. We also show that either the Gly137 to Asp point mutation in alphaE7 esterase gene or a deletion at this locus confer resistance and overproduction of the CYP6A1 protein. Based on these findings, we propose it is the absence of the wild-type Gly137 allele of the alphaE7 gene that releases the transcriptional repression of genes coding for detoxification enzymes such as CYP6A1, thereby leading to metabolic resistance to diazinon.

Human and mouse homologues of the Drosophila melanogaster tweety (tty) gene: a novel gene family encoding predicted transmembrane proteins.

We have cloned cDNA for TTYH1, a human homologue of the Drosophila melanogaster tweety (tty) gene. The 450-residue predicted protein shows 27% amino acid sequence identity (51% similarity) to the Drosophila protein, which contains an additional C-terminal repetitive region. A second Drosophila homologue exhibits 42% identity (65% similarity) to the tty protein. Mouse (Ttyh1), macaque, and Caenorhabditis elegans homologues were also identified, and the complete coding sequence for the mouse gene was determined. The mouse protein is 91% identical to the human protein. Hydrophobicity analysis of the tty-related proteins indicates that they represent a new family of membrane proteins with five potential membrane-spanning regions. The yeast FTR1 and FTH1 iron transporter proteins and the mammalian neurotensin receptors 1 and 2 have a similar hydrophobicity profile, although there is no detectable sequence homology to the tty-related proteins. This suggests that the tweety-related proteins could be involved in transport of iron or other divalent cations or alternatively that they may be membrane-bound receptors. TTYH1 was mapped to chromosome 19q13.4 by FISH and by radiation hybrid mapping using the Stanford G3 panel.

Reconstructing the diversification of alpha-esterases: comparing the gene clusters of Drosophila buzzatii and D. melanogaster.

A cluster composed of 10 active alpha-esterase genes and a pseudogene is distributed over 60 kb in the Drosophila melanogaster genome. This paper describes the corresponding cluster in Drosophila buzzatii, whose lineage diverged from that of D. melanogaster when the subgenera Drosophila and Sophophora diverged about 50 Mya. With three exceptions we find that the composition of the cluster is conserved in the two lineages. The location of alpha E1 in D. melanogaster differs from that of its nearest relative in D. buzzatii, and alpha E4 has duplicated independently in the two lineages. The nature of these differences indicates that a mechanism exists whereby copies of genes can be placed in opposite orientation and nonadjacent positions within a gene cluster, although this does not seem to be a feature of earlier events in the cluster's evolution. The rates of amino acid change are not significantly different between orthologs, but the rates differ sevenfold among paralogs, indicating that very different selective forces are acting on the genes of the cluster. Mapping of sequence differences onto a model of the tertiary structure of the enzymes indicates that motifs contributing to substrate binding and catalysis have changed radically in the alphaE4s and suggest that this subgroup of alpha-esterases may be evolving into a substantially different functional niche.

An episode of accelerated amino acid change in Drosophila esterase-6 associated with a change in physiological function.

In most lineages of the subgenus Sophophora esterase-6 is a homodimeric haemolymph protein. In the melanogaster subgroup of species it has become a monomer which is mainly expressed in the male sperm ejaculatory duct. Our analyses of esterase-6 sequences from three melanogaster subgroup species and two close relatives reveal a brief period of accelerated amino acid sequence change during the transition between the ancestral and derived states. In this period of 2-6Myr the ratio of replacement to silent site substitutions (0.51) is about three times higher than the values in other lineages of the phylogeny. There are about 50 more replacements in this period than would be predicted from the ratios of replacement to silent site substitutions found elsewhere in the phylogeny. Modelling on the known structure of a related acetylcholinesterase suggests that an unusually high proportion of the replacements in the transitional branch are non-conservative changes on the protein surface. Up to half the accelerated replacement rate can be accounted for by clusters of changes to the face of the molecule containing the opening of the active site gorge. This includes changes in and around regions homologous to peripheral substrate binding sites in acetylcholinesterase. There are also three changes in glycosylation status. One region predicted to lie on the protein surface which becomes markedly more hydrophilic is proposed to be the ancestral dimerisation site that is lost in the transitional branch.

Carboxyl/cholinesterases: a case study of the evolution of a successful multigene family.

The evolution of organismal diversity among the Metazoa is dependent on the proliferation of genes and diversification of functions in multigene families. Here we analyse these processes for one highly successful family, the carboxyl/cholinesterases. One key to the expansion of the functional niche of this group of enzymes is associated with versatile substrate binding and catalytic machinery. Qualitatively new functions can be obtained by substitution of one or a very few amino acids. This crudely adapted new functionality is then refined rapidly by a pulse of change elsewhere in the molecule; in one case about 13% amino acid divergence occurred in 5-10 million years. Furthermore, we postulate that the versatility of the substrate binding motifs underpins the recruitment of several family members to additional noncatalytic signal transduction functions.

The same amino acid substitution in orthologous esterases confers organophosphate resistance on the house fly and a blowfly.

Organophosphate (OP) insecticide resistance in certain strains of Musca domestica is associated with reduction in the carboxylesterase activity of a particular esterase isozyme. This has been attributed to a 'mutant ali-esterase hypothesis', which invokes a structural mutation to an ali-esterase resulting in the loss of its carboxylesterase activity but acquisition of OP hydrolase activity. It has been shown that the mutation in Lucilia cuprina is a Gly137-->Asp substitution in the active site of an esterase encoded by the Lc alpha E7 gene (Newcomb, R.D., Campbell, P.M., Ollis, D.L., Cheah, E., Russell, R.J., Oakeshott, J.G., 1997. A single amino acid substitution converts a carboxylesterase to an organophosphate hydrolase and confers insecticide resistance on a blowfly. Proc. Natl. Acad. Sci. USA 94, 7464-7468). We now report the cloning and characterisation of the orthologous M. domestica Md alpha E7 gene, including the sequencing of cDNAs from the OP resistant Rutgers and OP susceptible sbo and WHO strains. The Md alpha E7 gene has the same intron structure as Lc alpha E7 and encodes a protein with 76% amino acid identity to Lc alpha E7. Comparisons between susceptible and resistance alleles show resistance in M. domestica is associated with the same Gly137-->Asp mutation as in L. cuprina. Bacterial expression of the Rutgers allele shows its product has OP hydrolase activity. The data indicate identical catalytic mechanisms have evolved in orthologous Md alpha E7 and Lc alpha E7 molecules to endow diazinon-type resistance on the two species of higher Diptera.

Genomic structure, evolution, and expression of human FLII, a gelsolin and leucine-rich-repeat family member: overlap with LLGL.

The Drosophila melanogaster flightless-I gene is involved in cellularization processes in early embryogenesis and in the structural organization of indirect flight muscle. The encoded protein contains a gelsolin-like actin binding domain and an N-terminal leucine-rich repeat protein-protein interaction domain. The homologous human FLII gene encodes a 1269-residue protein with 58% amino acid sequence identity and is deleted in Smith-Magenis syndrome. We have cloned the FLII gene and determined its nucleotide sequence (14.1 kb). FLII has 29 introns, compared with 13 in Caenorhabditis elegans and 3 in D. melanogaster. The positions of several introns are conserved in FLII-related genes and in the domains and subdomains of the gelsolin-like regions giving indications of gelsolin gene family evolution. In keeping with its function in indirect flight muscle in Drosophila, the human FLII gene was most highly expressed in muscle. The FLII gene lies adjacent to LLGL, the human homologue of the D. melanogaster tumor suppressor gene lethal(2) giant larvae. The 3' end of the FLII transcript overlaps the 3' end of the LLGL transcript, and the corresponding mouse genes Fliih and Llglh also overlap. The overlap region contains poly(A) signals for both genes and is strongly conserved between human and mouse.

Biochemistry of esterases associated with organophosphate resistance in Lucilia cuprina with comparisons to putative orthologues in other Diptera.

Esterase activities associated with organophosphate insecticide resistance in the Australian sheep blowfly, Lucilia cuprina, are compared with similar activities in other Diptera. The enzymes making the major contribution to methyl butyrate hydrolysis ("ali-esterase") in L. cuprina, M. domestica, and D. melanogaster comigrate during electrophoresis. The enzymes in L. cuprina and D. melanogaster correspond to the naphthyl acetate hydrolyzing E3 and EST23 isozymes of those species. These and previously published data suggest that the ali-esterases of all three species are orthologous. Strains of L. cuprina fall into four groups on the basis of quantitative determinations of their ali-estesterase, OP hydrolase, and malathion carboxylesterase activities and these groups correspond to their status with respect to two types of OP resistance. Strains susceptible to OP's have high ali-esterase, low OP hydrolase, and intermediate MCE activities; those resistant to malathion but not diazinon have low ali-esterase, intermediate OP hydrolase, and high MCE activities; those resistant to diazinon but not malathion have low ali-esterase, high OP hydrolase, and low MCE activities; those resistant to both OPs have low ali-esterase, high OP hydrolase, and high MCE activities. The correlated changes among the three biochemical and two resistance phenotypes suggest that they are all properties of one gene/enzyme system; three major allelic variants of that system explain OP susceptibility and the two types of OP resistance. Models are proposed to explain the joint contribution of OP hydrolase and MCE activities to malathion resistance and the invariant association of low ali-esterase and elevated OP hydrolase activities in either type of resistance.

The novel flightless-I gene brings together two gene families, actin-binding proteins related to gelsolin and leucine-rich-repeat proteins involved in Ras signal transduction.

The Drosophila melanogaster gene flightless-I, involved in gastrulation and muscle degeneration, has Caenorhabditis elegans and human homologues. In these highly conserved genes, two previously known gene families have been brought together, families encoding the actin-binding proteins related to gelsolin and the leucine-rich-repeat (LRR) group of proteins involved in protein-protein interactions. Both these gene families exhibit characteristics of molecular changes involving replication slippage and exon shuffling. Phylogenetic analyses of 19 amino acid sequences of 6 related protein types indicate that actin-associated proteins related to gelsolin are monophyletic to a common ancestor and include flightless proteins. Conversely, comparison of 24 amino acid sequences of LRR proteins including the flightless proteins indicates that flightless proteins are members of a structurally related subgroup. Included in the flightless cluster are human and mouse rsp-1 proteins involved in suppressing v-Ras transformation of cells and the membrane-associated yeast (Saccharomyces cerevisae) adenylate cyclase whose analogous LRRs are required for interaction with Ras proteins. There is a strong possibility that ligands for this group could be related and that flightless may have a similar role in Ras signal transduction. It is hypothesized that an ancestral monomeric gelsolin precursor protein has undergone at least four independent gene reorganization events to account for the structural diversity of the extant family of gelsolin-related proteins and that gene duplication and exon shuffling events occurred prior to or at the beginning of multicellular life, resulting in the evolution of some members of the family soon after the appearance of actin-type proteins.

The Drosophila melanogaster flightless-I gene involved in gastrulation and muscle degeneration encodes gelsolin-like and leucine-rich repeat domains and is conserved in Caenorhabditis elegans and humans.

Mutations at the flightless-I locus (fliI) of Drosophila melanogaster cause flightlessness or, when severe, incomplete cellularization during early embryogenesis, with subsequent abnormalities in mesoderm invagination and in gastrulation. After chromosome walking, deficiency mapping, and transgenic analysis, we have isolated and characterized flightless-I cDNAs, enabling prediction of the complete amino acid sequence of the 1256-residue protein. Data base searches revealed a homologous gene in Caenorhabditis elegans, and we have isolated and characterized corresponding cDNAs. By using the polymerase chain reaction with nested sets of degenerate oligonucleotide primers based on conserved regions of the C. elegans and D. melanogaster proteins, we have cloned a homologous human cDNA. The predicted C. elegans and human proteins are, respectively, 49% and 58% identical to the D. melanogaster protein. The predicted proteins have significant sequence similarity to the actin-binding protein gelsolin and related proteins and, in addition, have an N-terminal domain consisting of a repetitive amphipathic leucine-rich motif. This repeat is found in D. melanogaster, Saccharomyces cerevisiae, and mammalian proteins known to be involved in cell adhesion and in binding to other proteins. The structure of the maternally expressed flightless-I protein suggests that it may play a key role in embryonic cellularization by interacting with both the cytoskeleton and other cellular components. The presence of a highly conserved homologue in nematodes, flies, and humans is indicative of a fundamental role for this protein in many metazoans.

Inositol 1,4,5-trisphosphate phosphatase deficiency and malignant hyperpyrexia in swine.

The sarcoplasmic reticulum from muscle of swine which are susceptible to malignant hyperpyrexia is deficient in inositol 1,4,5-trisphosphate phosphatase (InsP35-ase) activity, which leads to high intracellular concentrations of inositol 1,4,5-trisphosphate (InsP3) and of calcium ions. Halothane inhibits InsP35-ase and further increases myoplasmic InsP3 and calcium ion concentrations, and produces the clinical features of malignant hyperpyrexia.