Molecular Evolution of Drosophila Cuticular Protein Genes

Several multigene families have been described that together encode scores of structural cuticular proteins in Drosophila, although the functional significance of this diversity remains to be explored. Here I investigate the evolutionary histories of several multigene families (CPR, Tweedle, CPLCG, and CPF/CPFL) that vary in age, size, and sequence complexity, using sequenced Drosophila genomes and mosquito outgroups. My objective is to describe the rates and mechanisms of ‘cuticle-ome’ divergence, in order to identify conserved and rapidly evolving elements. I also investigate potential examples of interlocus gene conversion and concerted evolution within these families during Drosophila evolution. The absolute rate of change in gene number (per million years) is an order of magnitude lower for cuticular protein families within Drosophila than it is among Drosophila and the two mosquito taxa, implying that major transitions in the cuticle proteome have occurred at higher taxonomic levels. Several hotspots of intergenic conversion and/or gene turnover were identified, e.g. some gene pairs have independently undergone intergenic conversion within different lineages. Some gene conversion hotspots were characterized by conversion tracts initiating near nucleotide repeats within coding regions, and similar repeats were found within concertedly evolving cuticular protein genes in Anopheles gambiae. Rates of amino-acid substitution were generally severalfold higher along the branch connecting the Sophophora and Drosophila species groups, and 13 genes have Ka/Ks significantly greater than one along this branch, indicating adaptive divergence. Insect cuticular proteins appear to be a source of adaptive evolution within genera and, at higher taxonomic levels, subject to periods of gene-family expansion and contraction followed by quiescence. However, this relative stasis is belied by hotspots of molecular evolution, particularly concerted evolution, during the diversification of Drosophila. The prominent association between interlocus gene conversion and repeats within the coding sequence of interacting genes suggests that the latter promote strand exchange.     

Genetics of lipopolysaccharide biosynthesis in enteric bacteria.

From a historical perspective, the study of both the biochemistry and the genetics of lipopolysaccharide (LPS) synthesis began with the enteric bacteria. These organisms have again come to the forefront as the blocks of genes involved in LPS synthesis have been sequenced and analyzed. A number of new and unanticipated genes were found in these clusters, indicating a complexity of the biochemical pathways which was not predicted from the older studies. One of the most dramatic areas of LPS research has been the elucidation of the lipid A biosynthetic pathway. Four of the genes in this pathway have now been identified and sequenced, and three of them are located in a complex operon which also contains genes involved in DNA and phospholipid synthesis. The rfa gene cluster, which contains many of the genes for LPS core synthesis, includes at least 17 genes. One of the remarkable findings in this cluster is a group of several genes which appear to be involved in the synthesis of alternate rough core species which are modified so that they cannot be acceptors for O-specific polysaccharides. The rfb gene clusters which encode O-antigen synthesis have been sequenced from a number of serotypes and exhibit the genetic polymorphism anticipated on the basis of the chemical complexity of the O antigens. These clusters appear to have originated by the exchange of blocks of genes among ancestral organisms. Among the large number of LPS genes which have now been sequenced from these rfa and rfb clusters, there are none which encode proteins that appear to be secreted across the cytoplasmic membrane and surprisingly few which encode integral membrane proteins or proteins with extensive hydrophobic domains. These data, together with sequence comparison and complementation experiments across strain and species lines, suggest that the LPS biosynthetic enzymes may be organized into clusters on the inner surface of the cytoplasmic membrane which are organized around a few key membrane proteins.     

Concerted evolution of sea anemone neurotoxin genes is revealed through analysis of the Nematostella vectensis genome

Gene families, which encode toxins, are found in many poisonous animals, yet there is limited understanding of their evolution at the nucleotide level. The release of the genome draft sequence for the sea anemone Nematostella vectensis enabled a comprehensive study of a gene family whose neurotoxin products affect voltage-gated sodium channels. All gene family members are clustered in a highly repetitive approximately 30-kb genomic region and encode a single toxin, Nv1. These genes exhibit extreme conservation at the nucleotide level which cannot be explained by purifying selection. This conservation greatly differs from the toxin gene families of other animals (e.g., snakes, scorpions, and cone snails), whose evolution was driven by diversifying selection, thereby generating a high degree of genetic diversity. The low nucleotide diversity at the Nv1 genes is reminiscent of that reported for DNA encoding ribosomal RNA (rDNA) and 2 hsp70 genes from Drosophila, which have evolved via concerted evolution. This evolutionary pattern was experimentally demonstrated in yeast rDNA and was shown to involve unequal crossing-over. Through sequence analysis of toxin genes from multiple N. vectensis populations and 2 other anemone species, Anemonia viridis and Actinia equina, we observed that the toxin genes for each sea anemone species are more similar to one another than to those of other species, suggesting they evolved by manner of concerted evolution. Furthermore, in 2 of the species (A. viridis and A. equina) we found genes that evolved under diversifying selection, suggesting that concerted evolution and accelerated evolution may occur simultaneously.     

Analysis of homologous gene clusters in Caenorhabditis elegans reveals striking regional cluster domains

An algorithm for detecting local clusters of homologous genes was applied to the genome of Caenorhabditis elegans. Clusters of two or more homologous genes are abundant, totaling 1391 clusters containing 4607 genes, over one-fifth of all genes in C. elegans. Cluster genes are distributed unevenly in the genome, with the large majority located on autosomal chromosome arms, regions characterized by higher genetic recombination and more repeat sequences than autosomal centers and the X chromosome. Cluster genes are transcribed at much lower levels than average and very few have gross phenotypes as assayed by RNAi-mediated reduction of function. The molecular identity of cluster genes is unusual, with a preponderance of nematode-specific gene families that encode putative secreted and transmembrane proteins, and enrichment for genes implicated in xenobiotic detoxification and innate immunity. Gene clustering in Drosophila melanogaster is also substantial and the molecular identity of clustered genes follows a similar pattern. I hypothesize that autosomal chromosome arms in C. elegans undergo frequent local gene duplication and that these duplications support gene diversification and rapid evolution in response to environmental challenges. Although specific gene clusters have been documented in C. elegans, their abundance, genomic distribution, and unusual molecular identities were previously unrecognized.     

A novel polypeptide secreted by activated human T lymphocytes

We have identified two cDNA clones, I-309 and G-26, which define genes expressed abundantly in activated human PBMC, but at low or undetectable levels in resting PBMC. Based upon nucleotide sequence analysis, both clones are predicted to encode small, structurally related polypeptides, each containing a hydrophobic leader sequence characteristic of secreted proteins and a motif of four conserved cysteine residues. Further, I-309 and G-26 are structurally related to a growing family of genes that apparently encode small polypeptides whose secretion is induced upon cell activation. I-309 represents a previously undescribed human gene. We have generated an anti-peptide antiserum to the I-309 gene product which recognizes proteins in culture supernatants of an activated T cell clone and of COS cells transfected with the I-309 cDNA, supporting the idea that I-309 encodes a secreted protein. Because I-309 encodes a small protein secreted by activated T cells that displays structural features similar to other cytokines, we believe that it defines a novel cytokine with as yet unknown function.     

Purifying selection and birth-and-death evolution in the histone H4 gene family

Histones are small basic proteins encoded by a multigene family and are responsible for the nucleosomal organization of chromatin in eukaryotes. Because of the high degree of protein sequence conservation, it is generally believed that histone genes are subject to concerted evolution. However, purifying selection can also generate a high degree of sequence homogeneity. In this study, we examined the long-term evolution of histone H4 genes to determine whether concerted evolution or purifying selection was the major factor for maintaining sequence homogeneity. We analyzed the proportion (p(S)) of synonymous nucleotide differences between the H4 genes from 59 species of fungi, plants, animals, and protists and found that p(S) is generally very high and often close to the saturation level (p(S) ranging from 0.3 to 0.6) even though protein sequences are virtually identical for all H4 genes. A small proportion of genes showed a low level of p(S) values, but this appeared to be caused by recent gene duplication. Our findings suggest that the members of this gene family evolve according to the birth-and-death model of evolution under strong purifying selection. Using histone-like genes in archaebacteria as outgroups, we also showed that H1, H2A, H2B, H3, and H4 histone genes in eukaryotes form separate clusters and that these classes of genes diverged nearly at the same time, before the eukaryotic kingdoms diverged.     

Comparative genome analysis of filamentous fungi reveals gene family expansions associated with fungal pathogenesis

Fungi and oomycetes are the causal agents of many of the most serious diseases of plants. Here we report a detailed comparative analysis of the genome sequences of thirty-six species of fungi and oomycetes, including seven plant pathogenic species, that aims to explore the common genetic features associated with plant disease-causing species. The predicted translational products of each genome have been clustered into groups of potential orthologues using Markov Chain Clustering and the data integrated into the e-Fungi object-oriented data warehouse (http://www.e-fungi.org.uk/). Analysis of the species distribution of members of these clusters has identified proteins that are specific to filamentous fungal species and a group of proteins found only in plant pathogens. By comparing the gene inventories of filamentous, ascomycetous phytopathogenic and free-living species of fungi, we have identified a set of gene families that appear to have expanded during the evolution of phytopathogens and may therefore serve important roles in plant disease. We have also characterised the predicted set of secreted proteins encoded by each genome and identified a set of protein families which are significantly over-represented in the secretomes of plant pathogenic fungi, including putative effector proteins that might perturb host cell biology during plant infection. The results demonstrate the potential of comparative genome analysis for exploring the evolution of eukaryotic microbial pathogenesis.     

Identification of two novel clusters of ultrahigh-sulfur keratin-associated protein genes on human chromosome 11

We analyzed two novel clusters of keratin-associated protein (KAP) genes on human chromosome 11 (11p15.5 and 11q13.5) in which we identified two known human KRTAP5 genes, KerA (=KRN1) and KerB, and nine novel KRTAP5 family genes. RT-PCR analysis of these KAP genes showed preferential expression in human hair root, suggesting these gene products are required for hair formation. Based on the deduced amino acid sequences, all these KAP proteins were classified into an ultrahigh-sulfur (UHS) type KAP with high cysteine content (> 30 mol%). These KAPs also showed high glycine and serine contents (average 24.30 and 21.13 mol%, respectively), distinguishing from other UHS/HS KAP families located on human chromosomes 17 and 21. Dot-matrix analysis revealed a significant similarity between these two KAP gene clusters. We postulated a mechanism by which these two KAP gene clusters are generated via genomic duplication of a primordial gene cluster followed by genetic modification during evolution.     

Evolutionary history of x-tox genes in three lepidopteran species: Origin,evolution of primary and secondary structure and alternative splicing,generating a repertoire of immune-related proteins

The proteins of the X-tox family have imperfectly conserved tandem repeats of several defensin-like motifs known as cysteine-stabilized αβ (CS-αβ) motifs. These immune-related proteins are inducible and expressed principally in hemocytes, but they have lost the antimicrobial properties of the ancestral defensins from which they evolved. We compared x-tox gene structure and expression in three lepidopteran species (Spodoptera frugiperda, Helicoverpa armigera and Bombyx mori). Synteny and phylogenetic analyses showed that the x-tox exons encoding CS-αβ motifs were phylogenetically closely related to defensin genes mapping to chromosomal positions close to the x-tox genes. We were able to define two groups of paralogous x-tox exons (three in Noctuids) that each followed the expected species tree. These results suggest that the ancestor of the three species already possessed an x-tox gene with at least two proto-domains, and an additional duplication/fusion should have occurred in the ancestor of the two noctuid species. An expansion of the number of exons subsequently occurred in each lineage. Alternatively, the proto x-tox gene possessed more copy and each group of x-tox domains might undergo concerted evolution through gene conversion. Accelerated protein evolution was detected in x-tox domains when compared to related defensins, concomitantly to multiplication of exons and/or the possible activation of concerted evolution. The x-tox genes of the three species have similar structural organizations, with repeat motifs composed of CS-αβ-encoding exons flanked by introns in phase 1. Diverse mechanisms underlie this organization: (i) the acquisition of new repeat motifs, (ii) the duplication of preexisting repeat motifs and (iii) the duplication of modules. A comparison of gDNA and cDNA structures showed that alternative splicing results in the production of multiple X-tox protein isoforms from the x-tox genes. Differences in the number and sequence of CS-αβ motifs in these isoforms were found between species, but also between individuals of the same species. Thus, our analysis of the genetic organization and expression of x-tox genes in three lepidopteran species suggests a rapid evolution of the organization of these genes.     

Concerted evolution within a trypsin gene cluster in Drosophila.

A cluster of four trypsin genes has previously been localized to cytological position 47D-F of the Drosophila melanogaster genome. One of these genes had been sequenced, and the presence of the other three genes was identified by cross-hybridization. Here, we present the DNA sequence of the entire genomic region encoding these four trypsin genes. In addition to the four previously inferred genes, we have identified a fifth trypsin-coding sequence located within this gene cluster. This new gene shows a high degree of sequence divergence (more than 30%) from the other four genes, although it retains all of the functional motifs that are characteristic of trypsin-coding sequences. In order to trace the molecular evolution of this gene cluster, we isolated and sequenced the homologous 7-kb region from the closely related species Drosophila erecta. A comparison of the DNA sequences between the two species provides strong evidence for the concerted evolution of some members of this gene family. Two genes within the cluster are evolving in concert, while a third gene appears to be evolving independently. The remaining two genes show an intermediate pattern of evolution. We propose a simple model, involving chromosome looping and gene conversion, to explain the relatively complex patterns of molecular evolution within this gene cluster.     

Identification of Xanthomonas campestris pv. campestris galactose utilization genes from transcriptome data

A 70 mer oligonucleotide microarray was constructed to analyze genome-wide expression profiles of Xanthomonas campestris pv. campestris B100, a plant-pathogenic bacterium that is industrially employed to produce the exopolysaccharide xanthan gum which has many applications as a stabilizing, thickening, gelling, and emulsifying agent in food, pharmaceutical, and cosmetic industries. As an application example, global changes of gene expression were monitored during growth of X. campestris pv. campestris B100 on two different carbon sources. Exponential growing bacterial cultures were incubated either for 1h or permanently in minimal medium supplemented with 1% galactose in comparison to growth in minimal medium supplemented with 1% glucose. Six genes were identified that were significantly increased in gene expression under both growth conditions. These genes were located in three distinguished chromosomal regions in operon-like gene clusters. Genes from these clusters encode secreted glycosidases, which were predicted to be specific for galactose-containing carbohydrates, as well as transport proteins probably located in the outer and inner cell membrane. Finally genes from one cluster code for cytoplasmic enzymes of a metabolic pathway specific for the breakdown of galactose to intermediates of glycolysis.     

Recent duplications dominate NBS-encoding gene expansion in two woody species

Most disease resistance genes in plants encode NBS-LRR proteins. However, in woody species, little is known about the evolutionary history of these genes. Here, we identified 459 and 330 respective NBS-LRRs in grapevine and poplar genomes. We subsequently investigated protein motif composition, phylogenetic relationships and physical locations. We found significant excesses of recent duplications in perennial species, compared with those of annuals, represented by rice and Arabidopsis. Consequently, we observed higher nucleotide identity among paralogs and a higher percentage of NBS-encoding genes positioned in numerous clusters in the grapevine and poplar. These results suggested that recent tandem duplication played a major role in NBS-encoding gene expansion in perennial species. These duplication events, together with a higher probability of recombination revealed in this study, could compensate for the longer generation time in woody perennial species e.g. duplication and recombination could serve to generate novel resistance specificities. In addition, we observed extensive species-specific expansion in TIR-NBS-encoding genes. Non-TIR-NBS-encoding genes were poly- or paraphyletic, i.e. genes from three or more plant species were nested in different clades, suggesting different evolutionary patterns between these two gene types.     

The mitochondrial genome of Saprolegnia ferax: organization, gene content and nucleotide sequence

The mitochondrial genome of the peronosporomycete water mold Saprolegnia ferax has been characterized as a 46 930 bp circle containing an 8618 bp large inverted repeat (LIR). Eighteen reading frames encode identified subunits of respiratory complexes I, III, IV and V; 16 encode polypeptides of small and large mitoribosome subunits; and one encodes a subunit of the sec-independent protein translocation pathway. Of four additional putative reading frames three are homologues of those found in the related Phytophthora infestans genome. Protein encoding loci in the tightly compacted genome typically are arranged in operon-like clusters including three abutting and two overlapping pairs of reading frames. Translational RNAs include the mitochondrial small and large subunit rRNAs and 25 tRNA species. No tRNAs are encoded to enable translation of any threonine or the arginine CGR codons. The LIR separates the molecule into 19 274 bp large and 10 420 bp small single copy regions, and it encodes intact duplicate copies of four reading frames encoding known proteins, both rRNAs, and five tRNAs. Partial 3' sequences of three additional reading frames are duplicated at single copy sequence junctions. Active recombination between LIR elements generates two distinctive gene orders and uses the duplicated 3' sequences to maintain intact copies of the partially duplicated loci.     

Extensive gene amplification and concerted evolution within the CPR family of cuticular proteins in mosquitoes

Annotation of the Anopheles gambiae genome has revealed a large increase in the number of genes encoding cuticular proteins with the Rebers and Riddiford Consensus (the CPR gene family) relative to Drosophila melanogaster. This increase reflects an expansion of the RR-2 group of CPR genes, particularly the amplification of sets of highly similar paralogs. Patterns of nucleotide variation indicate that extensive concerted evolution is occurring within these clusters. The pattern of concerted evolution is complex, however, as sequence similarity within clusters is uncorrelated with gene order and orientation, and no comparable clusters occur within similarly compact arrays of the RR-1 group in mosquitoes or in either group in D. melanogaster. The dearth of pseudogenes suggests that sequence clusters are maintained by selection for high gene-copy number, perhaps due to selection for high expression rates. This hypothesis is consistent with the apparently parallel evolution of compact gene architectures within sequence clusters relative to single-copy genes. We show that RR-2 proteins from sequence-cluster genes have complex repeats and extreme amino-acid compositions relative to single-copy CPR proteins in An. gambiae, and that the amino-acid composition of the N-terminal and C-terminal sequence flanking the chitin-binding consensus region evolves in a correlated fashion.     

Molecular evolution of the nontandemly repeated genes of the histone 3 multigene family.

In some species, histone gene clusters consist of tandem arrays of each type of histone gene, whereas in other species the genes may be clustered but not arranged in tandem. In certain species, however, histone genes are found scattered across several different chromosomes. This study examines the evolution of histone 3 (H3) genes that are not arranged in large clusters of tandem repeats. Although H3 amino acid sequences are highly conserved both within and between species, we found that the nucleotide sequence divergence at synonymous sites is high, indicating that purifying selection is the major force for maintaining H3 amino acid sequence homogeneity over long-term evolution. In cases where synonymous-site divergence was low, recent gene duplication appeared to be a better explanation than gene conversion. These results, and other observations on gene inactivation, organization, and phylogeny, indicated that these H3 genes evolve according to a birth-and-death process under strong purifying selection. Thus, we found little evidence to support previous claims that all H3 proteins, regardless of their genome organization, undergo concerted evolution. Further analyses of the structure of H3 proteins revealed that the histones of higher eukaryotes might have evolved from a replication-independent-like H3 gene.     

Exploring the plant transcriptome through phylogenetic profiling

Publicly available protein sequences represent only a small fraction of the full catalog of genes encoded by the genomes of different plants, such as green algae, mosses, gymnosperms, and angiosperms. By contrast, an enormous amount of expressed sequence tags (ESTs) exists for a wide variety of plant species, representing a substantial part of all transcribed plant genes. Integrating protein and EST sequences in comparative and evolutionary analyses is not straightforward because of the heterogeneous nature of both types of sequence data. By combining information from publicly available EST and protein sequences for 32 different plant species, we identified more than 250,000 plant proteins organized in more than 12,000 gene families. Approximately 60% of the proteins are absent from current sequence databases but provide important new information about plant gene families. Analysis of the distribution of gene families over different plant species through phylogenetic profiling reveals interesting insights into plant gene evolution, and identifies species- and lineage-specific gene families, orphan genes, and conserved core genes across the green plant lineage. We counted a similar number of approximately 9,500 gene families in monocotyledonous and eudicotyledonous plants and found strong evidence for the existence of at least 33,700 genes in rice (Oryza sativa). Interestingly, the larger number of genes in rice compared to Arabidopsis (Arabidopsis thaliana) can partially be explained by a larger amount of species-specific single-copy genes and species-specific gene families. In addition, a majority of large gene families, typically containing more than 50 genes, are bigger in rice than Arabidopsis, whereas the opposite seems true for small gene families.     

Genomic clustering of cyanogenic glucoside biosynthetic genes aids their identification in Lotus japonicus and suggests the repeated evolution of this chemical defence pathway

Cyanogenic glucosides are amino acid-derived defence compounds found in a large number of vascular plants. Their hydrolysis by specific β-glucosidases following tissue damage results in the release of hydrogen cyanide. The cyanogenesis deficient1 (cyd1) mutant of Lotus japonicus carries a partial deletion of the CYP79D3 gene, which encodes a cytochrome P450 enzyme that is responsible for the first step in cyanogenic glucoside biosynthesis. The genomic region surrounding CYP79D3 contains genes encoding the CYP736A2 protein and the UDP-glycosyltransferase UGT85K3. In combination with CYP79D3, these genes encode the enzymes that constitute the entire pathway for cyanogenic glucoside biosynthesis. The biosynthetic genes for cyanogenic glucoside biosynthesis are also co-localized in cassava (Manihot esculenta) and sorghum (Sorghum bicolor), but the three gene clusters show no other similarities. Although the individual enzymes encoded by the biosynthetic genes in these three plant species are related, they are not necessarily orthologous. The independent evolution of cyanogenic glucoside biosynthesis in several higher plant lineages by the repeated recruitment of members from similar gene families, such as the CYP79s, is a likely scenario.     

Two large Arabidopsis thaliana gene families are homologous to the Brassica gene superfamily that encodes pollen coat proteins and the male component of the self-incompatibility response

The male component of the self-incompatibility response in Brassica has recently been shown to be encoded by the S locus cysteine-rich gene (SCR). SCR is related, at the sequence level, to the pollen coat protein (PCP) gene family whose members encode small, cysteine-rich proteins located in the proteo-lipidic surface layer (tryphine) of Brassica pollen grains. Here we show that the Arabidopsis genome includes two large gene families with homology to SCR and to the PCP gene family, respectively. These genes are poorly predicted by gene-identification algorithms and, with few exceptions, have been missed in previous annotations. Based on sequence comparison and an analysis of the expression patterns of several members of each family, we discuss the possible functions of these genes. In particular, we consider the possibility that SCR-related genes in Arabidopsis may encode ligands for the S gene family of receptor-like kinases in this species.