Dynamic evolution of the innate immune system in Drosophila

The availability of complete genome sequence from 12 Drosophila species presents the opportunity to examine how natural selection has affected patterns of gene family evolution and sequence divergence among different components of the innate immune system. We have identified orthologs and paralogs of 245 Drosophila melanogaster immune-related genes in these recently sequenced genomes. Genes encoding effector proteins, and to a lesser extent genes encoding recognition proteins, are much more likely to vary in copy number across species than genes encoding signaling proteins. Furthermore, we can trace the apparent recent origination of several evolutionarily novel immune-related genes and gene families. Using codon-based likelihood methods, we show that immune-system genes, and especially those encoding recognition proteins, evolve under positive darwinian selection. Positively selected sites within recognition proteins cluster in domains involved in recognition of microorganisms, suggesting that molecular interactions between hosts and pathogens may drive adaptive evolution in the Drosophila immune system.     

Distribution of fitness effects among beneficial mutations before selection in experimental populations of bacteria

The extent to which a population diverges from its ancestor through adaptive evolution depends on variation supplied by novel beneficial mutations. Extending earlier work, recent theory makes two predictions that seem to be robust to biological details: the distribution of fitness effects among beneficial mutations before selection should be (i) exponential and (ii) invariant, meaning it is always exponential regardless of the fitness rank of the wild-type allele. Here we test these predictions by assaying the fitness of 665 independently derived single-step mutations in the bacterium Pseudomonas fluorescens across a range of environments. We show that the distribution of fitness effects among beneficial mutations is indistinguishable from an exponential despite marked variation in the fitness rank of the wild type across environments. These results suggest that the initial step in adaptive evolution--the production of novel beneficial mutants from which selection sorts--is very general, being characterized by an approximately exponential distribution with many mutations of small effect and few of large effect. We also document substantial variation in the pleiotropic costs of antibiotic resistance, a result that may have implications for strategies aimed at eliminating resistant pathogens in animal and human populations.     

Loss of BBS proteins causes anosmia in humans and defects in olfactory cilia structure and function in the mouse

Defects in cilia are associated with several human disorders, including Kartagener syndrome, polycystic kidney disease, nephronophthisis and hydrocephalus. We proposed that the pleiotropic phenotype of Bardet-Biedl syndrome (BBS), which encompasses retinal degeneration, truncal obesity, renal and limb malformations and developmental delay, is due to dysfunction of basal bodies and cilia. Here we show that individuals with BBS have partial or complete anosmia. To test whether this phenotype is caused by ciliary defects of olfactory sensory neurons, we examined mice with deletions of Bbs1 or Bbs4. Loss of function of either BBS protein affected the olfactory, but not the respiratory, epithelium, causing severe reduction of the ciliated border, disorganization of the dendritic microtubule network and trapping of olfactory ciliary proteins in dendrites and cell bodies. Our data indicate that BBS proteins have a role in the microtubule organization of mammalian ciliated cells and that anosmia might be a useful determinant of other pleiotropic disorders with a suspected ciliary involvement.     

Mutation of the matrix metalloproteinase 2 gene (MMP2) causes a multicentric osteolysis and arthritis syndrome.

The inherited osteolyses or 'vanishing bone' syndromes are a group of rare disorders of unknown etiology characterized by destruction and resorption of affected bones. The multicentric osteolyses are notable for interphalangeal joint erosions that mimic severe juvenile rheumatoid arthritis (OMIMs 166300, 259600, 259610 and 277950). We recently described an autosomal recessive form of multicentric osteolysis with carpal and tarsal resorption, crippling arthritic changes, marked osteoporosis, palmar and plantar subcutaneous nodules and distinctive facies in a number of consanguineous Saudi Arabian families. We localized the disease gene to 16q12-21 by using members of these families for a genome-wide search for homozygous-by-descent microsatellite markers. Haplotype analysis narrowed the critical region to a 1.2-cM region that spans the gene encoding MMP-2 (gelatinase A, collagenase type IV; (ref. 3). We detected no MMP2 enzymatic activity in the serum or fibroblasts of affected family members. We identified two family-specific homoallelic MMP2 mutations: R101H and Y244X. The nonsense mutation effects a deletion of the substrate-binding and catalytic sites and the fibronectin type II-like and hemopexin/TIMP2 binding domains. Based on molecular modeling, the missense mutation disrupts hydrogen bond formation within the highly conserved prodomain adjacent to the catalytic zinc ion.     

Parallel adaptive origins of digestive RNases in Asian and African leaf monkeys

Similar morphological or physiological changes occurring in multiple evolutionary lineages are not uncommon. Such parallel changes are believed to be adaptive, because a complex character is unlikely to originate more than once by chance. However, the occurrence of adaptive parallel amino acid substitutions is debated. Here I propose four requirements for establishing adaptive parallel evolution at the protein sequence level and use these criteria to demonstrate such a case. I report that the gene encoding pancreatic ribonuclease was duplicated independently in Asian and African leaf-eating monkeys. Statistical analyses of DNA sequences, functional assays of reconstructed ancestral proteins and site-directed mutagenesis show that the new genes acquired enhanced digestive efficiencies through parallel amino acid replacements driven by darwinian selection. They also lost a non-digestive function independently, under a relaxed selective constraint. These results demonstrate that despite the overall stochasticity, even molecular evolution has a certain degree of repeatability and predictability under the pressures of natural selection.     

An empirical test of the mutational landscape model of adaptation using a single-stranded DNA virus

The primary impediment to formulating a general theory for adaptive evolution has been the unknown distribution of fitness effects for new beneficial mutations. By applying extreme value theory, Gillespie circumvented this issue in his mutational landscape model for the adaptation of DNA sequences, and Orr recently extended Gillespie's model, generating testable predictions regarding the course of adaptive evolution. Here we provide the first empirical examination of this model, using a single-stranded DNA bacteriophage related to phiX174, and find that our data are consistent with Orr's predictions, provided that the model is adjusted to incorporate mutation bias. Orr's work suggests that there may be generalities in adaptive molecular evolution that transcend the biological details of a system, but we show that for the model to be useful as a predictive or inferential tool, some adjustments for the biology of the system will be necessary.     

Mechanisms and convergence of compensatory evolution in mammalian mitochondrial tRNAs

The function of protein and RNA molecules depends on complex epistatic interactions between sites. Therefore, the deleterious effect of a mutation can be suppressed by a compensatory second-site substitution. In relating a list of 86 pathogenic mutations in human tRNAs encoded by mitochondrial genes to the sequences of their mammalian orthologs, we noted that 52 pathogenic mutations were present in normal tRNAs of one or several nonhuman mammals. We found at least five mechanisms of compensation for 32 pathogenic mutations that destroyed a Watson-Crick pair in one of the four tRNA stems: restoration of the affected Watson-Crick interaction (25 cases), strengthening of another pair (4 cases), creation of a new pair (8 cases), changes of multiple interactions in the affected stem (11 cases) and changes involving the interaction between the loop and stem structures (3 cases). A pathogenic mutation and its compensating substitution are fixed in a lineage in rapid succession, and often a compensatory interaction evolves convergently in different clades. At least 10%, and perhaps as many as 50%, of all nucleotide substitutions in evolving mammalian tRNAs participate in such interactions, indicating that the evolution of tRNAs proceeds along highly epistatic fitness ridges.     

Parallel bacterial evolution within multiple patients identifies candidate pathogenicity genes

Bacterial pathogens evolve during the infection of their human host(1-8), but separating adaptive and neutral mutations remains challenging(9-11). Here we identify bacterial genes under adaptive evolution by tracking recurrent patterns of mutations in the same pathogenic strain during the infection of multiple individuals. We conducted a retrospective study of a Burkholderia dolosa outbreak among subjects with cystic fibrosis, sequencing the genomes of 112 isolates collected from 14 individuals over 16 years. We find that 17 bacterial genes acquired nonsynonymous mutations in multiple individuals, which indicates parallel adaptive evolution. Mutations in these genes affect important pathogenic phenotypes, including antibiotic resistance and bacterial membrane composition and implicate oxygen-dependent regulation as paramount in lung infections. Several genes have not previously been implicated in pathogenesis and may represent new therapeutic targets. The identification of parallel molecular evolution as a pathogen spreads among multiple individuals points to the key selection forces it experiences within human hosts.     

Mutations in a member of the Ras superfamily of small GTP-binding proteins causes Bardet-Biedl syndrome

RAB, ADP-ribosylation factors (ARFs) and ARF-like (ARL) proteins belong to the Ras superfamily of small GTP-binding proteins and are essential for various membrane-associated intracellular trafficking processes. None of the approximately 50 known members of this family are linked to human disease. Using a bioinformatic screen for ciliary genes in combination with mutational analyses, we identified ARL6 as the gene underlying Bardet-Biedl syndrome type 3, a multisystemic disorder characterized by obesity, blindness, polydactyly, renal abnormalities and cognitive impairment. We uncovered four different homozygous substitutions in ARL6 in four unrelated families affected with Bardet-Biedl syndrome, two of which disrupt a threonine residue important for GTP binding and function of several related small GTP-binding proteins. Analysis of the Caenorhabditis elegans ARL6 homolog indicates that it is specifically expressed in ciliated cells, and that, in addition to the postulated cytoplasmic functions of ARL proteins, it undergoes intraflagellar transport. These findings implicate a small GTP-binding protein in ciliary transport and the pathogenesis of a pleiotropic disorder.     

A genome-wide RNA interference screen in Drosophila melanogaster cells for new components of the Hh signaling pathway

Members of the Hedgehog (Hh) family of signaling proteins are powerful regulators of developmental processes in many organisms and have been implicated in many human disease states. Here we report the results of a genome-wide RNA interference screen in Drosophila melanogaster cells for new components of the Hh signaling pathway. The screen identified hundreds of potential new regulators of Hh signaling, including many large protein complexes with pleiotropic effects, such as the coat protein complex I (COPI) complex, the ribosome and the proteasome. We identified the multimeric protein phosphatase 2A (PP2A) and two new kinases, the D. melanogaster orthologs of the vertebrate PITSLRE and cyclin-dependent kinase-9 (CDK9) kinases, as Hh regulators. We also identified a large group of constitutive and alternative splicing factors, two nucleoporins involved in mRNA export and several RNA-regulatory proteins as potent regulators of Hh signal transduction, indicating that splicing regulation and mRNA transport have a previously unrecognized role in Hh signaling. Finally, we showed that several of these genes have conserved roles in mammalian Hh signaling.     

The gene defective in leukocyte adhesion deficiency II encodes a putative GDP-fucose transporter

Leukocyte adhesion deficiency II (LAD II) is characterized by the lack of fucosylated glycoconjugates, including selectin ligands, causing immunodeficiency and severe mental and growth retardation. No deficiency in fucosyltransferase activities or in the activities of enzymes involved in GDP-fucose biosynthesis has been found. Instead, the transport of GDP-fucose into isolated Golgi vesicles of LAD II cells appeared to be reduced. To identify the gene mutated in LAD II, we cloned 12 cDNAs from Caenorhabditis elegans, encoding multi-spanning transmembrane proteins with homology to known nucleotide sugar transporters, and transfected them into fibroblasts from an LAD II patient. One of these clones re-established expression of fucosylated glycoconjugates with high efficiency and allowed us to identify a human homolog with 55% identity, which also directed re-expression of fucosylated glycoconjugates. Both proteins were localized to the Golgi. The corresponding endogenous protein in LAD II cells had an R147C amino acid change in the conserved fourth transmembrane region. Overexpression of this mutant protein in cells from a patient with LAD II did not rescue fucosylation, demonstrating that the point mutation affected the activity of the protein. Thus, we have identified the first putative GDP-fucose transporter, which has been highly conserved throughout evolution. A point mutation in its gene is responsible for the disease in this patient with LAD II.     

Mutant chromatin remodeling protein SMARCAL1 causes Schimke immuno-osseous dysplasia.

Schimke immuno-osseous dysplasia (SIOD, MIM 242900) is an autosomal-recessive pleiotropic disorder with the diagnostic features of spondyloepiphyseal dysplasia, renal dysfunction and T-cell immunodeficiency. Using genome-wide linkage mapping and a positional candidate approach, we determined that mutations in SMARCAL1 (SWI/SNF2-related, matrix-associated, actin-dependent regulator of chromatin, subfamily a-like 1), are responsible for SIOD. Through analysis of data from persons with SIOD in 26 unrelated families, we observed that affected individuals from 13 of 23 families with severe disease had two alleles with nonsense, frameshift or splicing mutations, whereas affected individuals from 3 of 3 families with milder disease had a missense mutation on each allele. These observations indicate that some missense mutations allow retention of partial SMARCAL1 function and thus cause milder disease.     

Mutations in NYX, encoding the leucine-rich proteoglycan nyctalopin, cause X-linked complete congenital stationary night blindness

During development, visual photoreceptors, bipolar cells and other neurons establish connections within the retina enabling the eye to process visual images over approximately 7 log units of illumination. Within the retina, cells that respond to light increment and light decrement are separated into ON- and OFF-pathways. Hereditary diseases are known to disturb these retinal pathways, causing either progressive degeneration or stationary deficits. Congenital stationary night blindness (CSNB) is a group of stable retinal disorders that are characterized by abnormal night vision. Genetic subtypes of CSNB have been defined and different disease actions have been postulated. The molecular bases have been elucidated in several subtypes, providing a better understanding of the disease mechanisms and developmental retinal neurobiology. Here we have studied 22 families with 'complete' X-linked CSNB (CSNB1; MIM 310500; ref. 4) in which affected males have night blindness, some photopic vision loss and a defect of the ON-pathway. We have found 14 different mutations, including 1 founder mutation in 7 families from the United States, in a novel candidate gene, NYX. NYX, which encodes a glycosylphosphatidyl (GPI)-anchored protein called nyctalopin, is a new and unique member of the small leucine-rich proteoglycan (SLRP) family. The role of other SLRP proteins suggests that mutant nyctalopin disrupts developing retinal interconnections involving the ON-bipolar cells, leading to the visual losses seen in patients with complete CSNB.     

Heterozygous mutations in the gene encoding noggin affect human joint morphogenesis

The secreted polypeptide noggin (encoded by the Nog gene) binds and inactivates members of the transforming growth factor beta superfamily of signalling proteins (TGFbeta-FMs), such as BMP4 (ref. 1). By diffusing through extracellular matrices more efficiently than TGFbeta-FMs, noggin may have a principal role in creating morphogenic gradients. During mouse embryogenesis, Nog is expressed at multiple sites, including developing bones. Nog-/- mice die at birth from multiple defects that include bony fusion of the appendicular skeleton. We have identified five dominant human NOG mutations in unrelated families segregating proximal symphalangism (SYM1; OMIM 185800) and a de novo mutation in a patient with unaffected parents. We also found a dominant NOG mutation in a family segregating multiple synostoses syndrome (SYNS1; OMIM 186500); both SYM1 and SYNS1 have multiple joint fusion as their principal feature. All seven NOG mutations alter evolutionarily conserved amino acid residues. The findings reported here confirm that NOG is essential for joint formation and suggest that NOG requirements during skeletogenesis differ between species and between specific skeletal elements within species.     

Loss-of-function mutations in LEMD3 result in osteopoikilosis, Buschke-Ollendorff syndrome and melorheostosis

Osteopoikilosis, Buschke-Ollendorff syndrome (BOS) and melorheostosis are disorders characterized by increased bone density. The occurrence of one or more of these phenotypes in the same individual or family suggests that these entities might be allelic. We collected data from three families in which affected individuals had osteopoikilosis with or without manifestations of BOS or melorheostosis. A genome-wide linkage analysis in these families, followed by the identification of a microdeletion in an unrelated individual with these diseases, allowed us to map the gene that is mutated in osteopoikilosis. All the affected individuals that we investigated were heterozygous with respect to a loss-of-function mutation in LEMD3 (also called MAN1), which encodes an inner nuclear membrane protein. A somatic mutation in the second allele of LEMD3 could not be identified in fibroblasts from affected skin of an individual with BOS and an individual with melorheostosis. XMAN1, the Xenopus laevis ortholog, antagonizes BMP signaling during embryogenesis. In this study, LEMD3 interacted with BMP and activin-TGFbeta receptor-activated Smads and antagonized both signaling pathways in human cells.     

Mutations in the skeletal muscle alpha-actin gene in patients with actin myopathy and nemaline myopathy.

Muscle contraction results from the force generated between the thin filament protein actin and the thick filament protein myosin, which causes the thick and thin muscle filaments to slide past each other. There are skeletal muscle, cardiac muscle, smooth muscle and non-muscle isoforms of both actin and myosin. Inherited diseases in humans have been associated with defects in cardiac actin (dilated cardiomyopathy and hypertrophic cardiomyopathy), cardiac myosin (hypertrophic cardiomyopathy) and non-muscle myosin (deafness). Here we report that mutations in the human skeletal muscle alpha-actin gene (ACTA1) are associated with two different muscle diseases, 'congenital myopathy with excess of thin myofilaments' (actin myopathy) and nemaline myopathy. Both diseases are characterized by structural abnormalities of the muscle fibres and variable degrees of muscle weakness. We have identified 15 different missense mutations resulting in 14 different amino acid changes. The missense mutations in ACTA1 are distributed throughout all six coding exons, and some involve known functional domains of actin. Approximately half of the patients died within their first year, but two female patients have survived into their thirties and have children. We identified dominant mutations in all but 1 of 14 families, with the missense mutations being single and heterozygous. The only family showing dominant inheritance comprised a 33-year-old affected mother and her two affected and two unaffected children. In another family, the clinically unaffected father is a somatic mosaic for the mutation seen in both of his affected children. We identified recessive mutations in one family in which the two affected siblings had heterozygous mutations in two different exons, one paternally and the other maternally inherited. We also identified de novo mutations in seven sporadic probands for which it was possible to analyse parental DNA.