Tribolium Wnts: evidence for a larger repertoire in insects with overlapping expression patterns that suggest multiple redundant functions in embryogenesis

Wingless (wg)/Wnt family genes encode secreted glycoproteins that function as signalling molecules in the development of vertebrates as well as invertebrates. In a survey of Wnt family genes in the newly sequenced Tribolium genome, we found a total of nine Wnt genes. In addition to wg or Wnt1, Tribolium contains orthologs of the vertebrate Wnt5-7 and Wnt9-11 genes. As in Drosophila, Wnt1, Wnt6 and Wnt10 are clustered in the genome. Comparative genomics indicates that Wnt9 is also a conserved member of this cluster in several insects for which genome sequence is available. One of the Tribolium Wnt genes appears to be a member of the WntA family, members of which have been identified in Anopheles and other invertebrates but not in Drosophila or vertebrates. Careful phylogenetic examination suggests an Apis Wnt gene, previously identified as a Wnt4 homolog, is also a member of the WntA family. The ninth Tribolium Wnt gene is related to the diverged Drosophila WntD gene, both of which phylogenetically group with Wnt8 genes. Some of the Tribolium Wnt genes display multiple overlapping expression patterns, suggesting that they may be functionally redundant in segmentation, brain, appendage and hindgut development. In contrast, the unique expression patterns of Wnt5, Wnt7 and Wnt11 in developing appendages likely indicate novel functions.     

Gene Arrangement within the Unique Long Genome Region of Infectious Laryngotracheitis Virus Is Distinct from That of Other Alphaherpesviruses

The genome of the avian alphaherpesvirus infectious laryngotracheitis virus (ILTV) comprises ca. 155 kbp of which ca. one-third have been sequenced so far. To gain additional sequence information we analyzed two stretches of 15.5 and 1.9 kbp of the ILTV unique long (U_L) genome region. The larger fragment contains homologs of the herpes simplex virus (HSV) UL23 (thymidine kinase) and UL22 (glycoprotein H) genes followed by five open reading frames (ORF) encoding putative proteins of 334 to 410 amino acids which exhibit no homology to any known herpesvirus protein. RNA analyses showed that these unique ILTV genes are indeed expressed. An origin of replication separates this cluster of unique genes from a conserved gene cluster consisting of the UL45, UL46, UL48, UL49, UL49.5, and UL50 homologs. The absence of UL47 from this position coincides with the localization of a UL47-homologous ORF within the unique short (U_S) region of the ILTV genome (M. Wild, S. Cook, and M. Cochran, Virus Genes 12:107–116, 1996). Within the second analyzed region the ILTV UL21 homolog was found adjacent to the UL44 gene. We thus identified five novel herpesvirus genes in ILTV and present evidence for a large internal inversion in the ILTV U_L region, in contrast to the collinear genomes of other alphaherpesviruses. Interestingly, a similar inversion is also present in the porcine alphaherpesvirus pseudorabies virus.     

A novel RNA helicase gene tightly linked to the Triplo-lethal locus of Drosophila.

The Triplo-lethal (Tpl) locus of Drosophila is the only known locus which is lethal when present in three copies rather than the normal two. After recovering a hybrid-dysgenesis-induced mutation of Tpl we used a rapid combination of transposon tagging, chromosome microdissection and PCR to clone the P element that had transposed into the Tpl region. That P element is located within the gene for a new and unique member of the RNA helicase family. This new helicase differs from all others known by having glycine-rich repeats at both the amino and carboxyl termini. Curiously, genetic analysis shows that the P element inserted into this gene is not responsible for the Tpl mutant phenotype. We present possible explanations for these findings.     

Identification of the Anopheles gambiae ATP-binding cassette transporter superfamily genes

The Anopheles gambiae genome sequence has been analyzed to find ATP-binding cassette protein genes based on deduced protein similarity to known family members. A nonredundant collection of 44 putative genes was identified including five genes not detected by the original Anopheles genome project machine annotation. These genes encode at least one member of all the human and Drosophila melanogaster ATP-binding protein subgroups. Like D. melanogaster, A. gambiae has subgroup ABCH genes encoding proteins different from the ABC proteins found in other complex organisms. The largest Anopheles subgroup is the ABCC genes which includes one member that can potentially encode ten different isoforms of the protein by differential splicing. As with Drosophila, the second largest Anopheles group is the ABCG subgroup with 12 genes compared to 15 genes in D. melanogaster, but only 5 genes in the human genome. In contrast, fewer ABCA and ABCB genes were identified in the mosquito genome than in the human or Drosophila genomes. Gene duplication is very evident in the Anopheles ABC genes with two groups of four genes, one group with three genes and three groups with two head to tail duplicated genes. These characteristics argue that the A. gambiae is actively using gene duplication as a mechanism to drive genetic variation in this important gene group.     

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.     

Phylogenetic analysis of the rapidly evolving SuUR gene in insects

Different genome regions differ in replication timing during the S phase. Late-replicating sequences are often underreplicated in the Drosophila salivary-gland polytene chromosomes. The SuUR gene, whose mutation changes the replication time of late-replicating regions in salivary-gland cells, has been identified in Drosophila melanogaster. The SUUR protein lacks homologs by a BLAST search, and only moderate homology is observed between its N-terminal end and chromatin-remodeling proteins of the SWI2/SNF2 family. The gene and the protein were analyzed in insects. Orthologs of the SuUR gene were found in all annotated Drosophila species. The number of amino acid substitutions in the SUUR protein proved to be extremely high, corresponding to that of rapidly evolving genes. Orthologs with low homology were found in mosquitoes Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus. No orthologs of the SuUR gene were detected beyond Diptera.     

The Gr family of candidate gustatory and olfactory receptors in the yellow-fever mosquito Aedes aegypti

The gustatory receptor (Gr) protein family contains most of the diversity in the insect chemoreceptor superfamily, including within it not only taste receptors but select olfactory receptors as well. Manual annotation of the Gr family in the genome sequence of the yellow-fever mosquito, Aedes aegypti, yielded a total of 114 potential proteins encoded by 79 genes. In the sequenced genome, 23 of these genes and protein isoforms are pseudogenic, leaving 91 putatively functional Grs. Comparison with our previously published set of 76 Grs encoded by 52 genes in the distantly related Anopheles gambiae mosquito revealed 13 new AgGrs encoded by 8 genes. Phylogenetic analysis reveals the conservation of carbon dioxide, sugar, and several orphan receptors in these 2 mosquitoes and Drosophila flies. On the other hand, most of these Grs are unique to mosquitoes and many are specific to the Aedes or Anopheles lineages, indicating their involvement in mosquito-specific aspects of both gustatory and olfactory perception. In particular, most instances of alternative splicing in orthologous loci appear to have evolved after the culicine-anopheline split +/-150 million years ago.     

The Triplo-Lethal Locus of Drosophila: Reexamination of Mutants and Discovery of a Second-Site Suppressor

In the genome of Drosophila melanogaster there is a single locus, Triplo-lethal (Tpl), that causes lethality when present in either one or three copies in an otherwise diploid animal. Previous attempts to mutagenize Tpl produced alleles that were viable over a chromosome bearing a duplication of Tpl, but were not lethal in combination with a wild-type chromosome, as deficiencies for Tpl are. These mutations were interpreted as hypomorphic alleles of Tpl. In this work, we show that these alleles are not mutations at Tpl; rather, they are dominant mutations in a tightly linked, but cytologically distant, locus that we have named Suppressor-of-Tpl (Su(Tpl)). Su(Tpl) mutations suppress the lethality associated with three copies of the Triplo-lethal locus and are recessive lethal. We have mapped Su(Tpl) to the approximate map position 3-46.5, within the cytological region 76B-76D.     

The Unusual Spectrum of Mutations Induced by Hybrid Dysgenesis at the Triplo-Lethal Locus of Drosophila Melanogaster

The Triplo-lethal locus (Tpl) is unique in its dosage sensitivity; no other locus in Drosophila has been identified that is lethal when present in three doses. Tpl is also haplo-lethal, and its function is still a mystery. Previous workers have found it nearly impossible to mutationally inactive Tpl other than by completely deleting the chromosomal region in which Tpl resides (83DE). We have utilized P-M hybrid dysgenesis in an effort to obtain new mutations of Tpl. We recovered 19 new duplications of Tpl, 15 hypomorphic mutations of Tpl (a previously rare class of mutation), and no null mutations. Surprisingly, 14 of the 15 hypomorphic alleles have no detectable P element sequences at the locus. The difficulty in recovering null mutations in Tpl suggests that it may be a complex locus, perhaps consisting of several genes with redundant functions. The relative ease with which we recovered hypomorphic alleles is in sharp contrast to previous attempts by others to mutagenize Tpl. A higher mutation rate with hybrid dysgenesis than with radiation or chemicals also suggests a peculiar genetic organization for the locus.     

Phospholipase C-gamma contains introns shared by src homology 2 domains in many unrelated proteins

Many proteins with novel functions were created by exon shuffling around the time of the metazoan radiation. Phospholipase C-gamma (PLC-gamma) is typical of proteins that appeared at this time, containing several different modules that probably originated elsewhere. To gain insight into both PLC-gamma evolution and structure-function relationships within the Drosophila PLC-gamma encoded by small wing (sl), we cloned and sequenced the PLC-gamma homologs from Drosophila pseudoobscura and D. virilis and compared their gene structure and predicted amino acid sequences with PLC-gamma homologs in other animals. PLC-gamma has been well conserved throughout, although structural differences suggest that the role of tyrosine phosphorylation in enzyme activation differs between vertebrates and invertebrates. Comparison of intron positions demonstrates that extensive intron loss has occurred during invertebrate evolution and also reveals the presence of conserved introns in both the N- and C-terminal PLC-gamma SH2 domains that are present in SH2 domains in many other genes. These and other conserved SH2 introns suggest that the SH2 domains in PLC-gamma are derived from an ancestral domain that was shuffled not only into PLC-gamma, but also into many other unrelated genes during animal evolution.     

Comparative analysis of BAC and whole genome shotgun sequences from an Anopheles gambiae region related to Plasmodium encapsulation

The only natural mechanism of malaria transmission in sub-Saharan Africa is the mosquito, generally Anopheles gambiae. Blocking malaria parasite transmission by stopping the development of Plasmodium in the insect vector would provide a useful alternative to the current methods of malaria control. Toward this end, it is important to understand the molecular basis of the malaria parasite refractory phenotype in An. gambiae mosquito strains. We have selected and sequenced six bacterial artificial chromosome (BAC) clones from the Pen-1 region that is the major quantitative trait locus involved in Plasmodium encapsulation. The sequence and the annotation of five overlapping BAC clones plus one adjacent, but not contiguous clone, totaling 585kb of genomic sequence from the centromeric end of the Pen-1 region of the PEST strain were compared to that of the genome sequence of the same strain produced by the whole genome shotgun technique. This project identified 23 putative mosquito genes plus putative copies of the retrotransposable elements BEL12 and TRANSIBN1_AG in the six BAC clones. Nineteen of the predicted genes are most similar to their Drosophila melanogaster homologs while one is more closely related to vertebrate genes. Comparison of these new BAC sequences plus previously published BAC sequences to the cognate region of the assembled genome sequence identified three retrotransposons present in one sequence version but not the other. One of these elements, Indy, has not been previously described. These observations provide evidence for the recent active transposition of these elements and demonstrate the plasticity of the Anopheles genome. The BAC sequences strongly support the public whole genome shotgun assembly and automatic annotation while also demonstrating the benefit of complementary genome sequences and of human curation. Importantly, the data demonstrate the differences in the genome sequence of an individual mosquito compared to that of a hypothetical, average genome sequence generated by whole genome shotgun assembly.     

New selenoproteins identified in silico from the genome of Anopheles gambiae

Selenoprotein is biosynthesized by the incorporation of selenocysteine into proteins,where the TGA codon in the open reading frame does not act as a stop signal but is translated into selenocysteine.The dual functions of TGA result in mis-annotation or lack of selenoproteins in the sequenced genomes of many species.Available computational tools fail to correctly predict selenoproteins.Thus,we devel-oped a new method to identify selenoproteins from the genome of Anopheles gambiae computationally.Based on released genomic information,several programs were edited with PERL language to identify selenocysteine insertion sequence(SECIS)element,the coding potential of TGA codons,and cys-teine-containing homologs of selenoprotein genes.Our results showed that 11365 genes were termi-nated with TGA codons,918 of which contained SECIS elements.Similarity search revealed that 58 genes contained Sec/Cys pairs and similar flanking regions around in-frame TGA codons.Finally,7 genes were found to fully meet requirements for selenoproteins,although they have not been anno-tated as selenoproteins in NCBI databases.Deduced from their basic properties,the newly found se-lenoproteins in the genome of Anopheles gambiae are possibly related to in vivo oxidation tolerance and protein regulation in order to interfere with anopheles' vectorial capacity of Plasmodium.This study may also provide theoretical bases for the prevention of malaria from anopheles transmission.     

New selenoproteins identified in silico from the genome of Anopheles gambiae

Selenoprotein is biosynthesized by the incorporation of selenocysteine into proteins, where the TGA codon in the open reading frame does not act as a stop signal but is translated into selenocysteine. The dual functions of TGA result in mis-annotation or lack of selenoproteins in the sequenced genomes of many species. Available computational tools fail to correctly predict selenoproteins. Thus, we developed a new method to identify selenoproteins from the genome of Anopheles gambiae computationally.Based on released genomic information, several programs were edited with PERL language to identify selenocysteine insertion sequence (SECIS) element, the coding potential of TGA codons, and cysteine-containing homologs of selenoprotein genes. Our results showed that 11365 genes were terminated with TGA codons, 918 of which contained SECIS elements. Similarity search revealed that 58genes contained Sec/Cys pairs and similar flanking regions around in-frame TGA codons. Finally, 7genes were found to fully meet requirements for selenoproteins, although they have not been annotated as selenoproteins in NCBI databases. Deduced from their basic properties, the newly found selenoproteins in the genome of Anopheles gambiae are possibly related to in vivo oxidation tolerance and protein regulation in order to interfere with anopheles' vectorial capacity of Plasmodium. This study may also provide theoretical bases for the prevention of malaria from anopheles transmission.     

A peritrophin-like protein expressed in the embryonic tracheae of Drosophila melanogaster.

We have cloned and sequenced a cDNA from Drosophila melanogaster that encodes a protein homologous to the peritrophins, a family of chitin-binding proteins from the peritrophic matrix of insects. Unexpectedly, the gene, Gasp, is expressed in the embryonic tracheae. We suggest that this family of proteins may be present in other tissues than the peritrophic matrix, particularly where nutrient or gas exchange are important, and/or where invasion by parasites or viruses is possible. We have also mapped two similar genes that had been sequenced by the Berkeley Drosophila Genome Project, and find that these three very similar genes are not clustered, but are located on three different chromosomes.     

Complete DNA sequence of the rat cytomegalovirus genome

We have determined the complete genome sequence of the Maastricht strain of rat cytomegalovirus (RCMV). The RCMV genome has a length of 229,896 bp and is arranged as a single unique sequence flanked by 504-bp terminal direct repeats. RCMV was found to have counterparts of all but one of the open reading frames (ORFs) that are conserved between murine CMV (MCMV) and human CMV (HCMV). Like HCMV, RCMV lacks homologs of the genes belonging to the MCMV m02 glycoprotein gene family. However, RCMV contains 15 ORFs with homology to members of the MCMV m145 glycoprotein gene family. Four ORFs are predicted to encode homologs of host proteins; R33 and R78 both putatively encode G protein-coupled receptors, whereas r144 and r131 encode homologs of major histocompatibility class I heavy chains and CC chemokines, respectively. An intriguing feature of the RCMV genome is the presence of an ORF, r127, with similarity to the rep gene of parvoviruses as well as ORF U94 of human herpesvirus 6A (HHV-6A) and HHV-6B. Counterparts of these ORFs have not been found in the other sequenced herpesviruses.     

Computational identification of novel chitinase-like proteins in the Drosophila melanogaster genome

MOTIVATION: Multiple chitinases as well as lectins closely related to them have been characterized previously from many insect species and the corresponding genes/cDNAs have been cloned. However, the identification of the entire assortment of genes for chitinase family proteins and their differences in biochemical properties have not been carried out in any individual insect species. The completion of the entire DNA sequence of Drosophila melanogaster (fruit fly) genome and identification of open reading frames presents an opportunity to study the structures and functions of chitinase-like proteins, and also to identify new members of this family in DROSOPHILA: We are, therefore, interested in studying the functional genomics of chitinase-like gene families in insects. METHODS: We searched the Drosophila protein sequences database using fully characterized insect chitinase sequences and BLASTP software, identified all the putative chitinase-like proteins encoded in Drosophila genome, and predicted their structures using domain analysis tools. A phylogenetic analysis of the chitinase-like proteins from Drosophila and several other insect species was carried out. The structures of these chitinases were modeled using homology modeling software. RESULTS: Our analysis revealed the presence of 18 chitinase-like proteins in the Drosophila protein database. Among these are seven novel chitinase-like proteins that contain four signature amino acid sequences of chitinases belonging to family 18 glycosylhydrolases, including both acidic and hydrophobic amino acid residues critical for enzyme activity. All the proteins contain at least one catalytic domain with one having four catalytic domains. Phylogenetic analysis of chitinase-like proteins from Drosophila and other insects revealed an evolutionary relationship among all these proteins, which indicated gene duplication and domain shuffling to generate the observed diversity in the encoded proteins. Homology modeling showed that all the Drosophila chitinase-like proteins contain one or more catalytic domains with a (alpha/beta)8 barrel-like structure. Our results suggest that insects utilize multiple family 18 chitinolytic enzymes and also non-enzymatic chitinase-like proteins for degrading/remodeling/binding to chitin in different insect anatomical extracellular structures, such as the cuticle, peritrophic membrane, trachea and mouth parts during insect development, and possibly for other roles including chitin synthesis. AVAILABILITY: Perl program and supplementary material are available at http://www.ksu.edu/bioinformatics/supplementary.htm     

A Recombinational Hotspot at the Triplo-Lethal Locus of Drosophila Melanogaster

In the genome of Drosophila melanogaster there is only one locus, Tpl, that is triplo-lethal; it is also haplo-lethal. Previous work has identified 3 hypomorphic alleles of Tpl which rescue animals carrying a duplication of Tpl, but which are not dominant lethals as null mutations or deficiencies would be. We have found that all three hypomorphic alleles act as site-specific hotspots for recombination when heterozygous with a wild-type homolog. Recombination between the flanking markers ri and Ki is increased 6.5-10.5-fold in the presence of Tpl hypomorphic alleles. The increased recombination was found to occur between Tpl and Ki, while recombination in other adjacent regions is unchanged. The use of isogenic Tpl+ controls, and the use of flanking intervals in the mutant chromosomes allows us to rule out the interchromosomal effect as a cause. We have also observed premeiotic recombination occurring at the Tpl hypomorphic alleles in male heterozygotes. We hypothesize that transposons are responsible for both the hypomorphic phenotype and the high frequency of recombination.     

Acyl-CoA Dehydrogenases: Dynamic History of Protein Family Evolution

The acyl-CoA dehydrogenases (ACADs) are enzymes that catalyze the α,β-dehydrogenation of acyl-CoA esters in fatty acid and amino acid catabolism. Eleven ACADs are now recognized in the sequenced human genome, and several homologs have been reported from bacteria, fungi, plants, and nematodes. We performed a systematic comparative genomic study, integrating homology searches with methods of phylogenetic reconstruction, to investigate the evolutionary history of this family. Sequence analyses indicate origin of the family in the common ancestor of Archaea, Bacteria, and Eukaryota, illustrating its essential role in the metabolism of early life. At least three ACADs were already present at that time: ancestral glutaryl-CoA dehydrogenase (GCD), isovaleryl-CoA dehydrogenase (IVD), and ACAD10/11. Two gene duplications were unique to the eukaryotic domain: one resulted in the VLCAD and ACAD9 paralogs and another in the ACAD10 and ACAD11 paralogs. The overall patchy distribution of specific ACADs across the tree of life is the result of dynamic evolution that includes numerous rounds of gene duplication and secondary losses, interdomain lateral gene transfer events, alteration of cellular localization, and evolution of novel proteins by domain acquisition. Our finding that eukaryotic ACAD species are more closely related to bacterial ACADs is consistent with endosymbiotic origin of ACADs in eukaryotes and further supported by the localization of all nine previously studied ACADs in mitochondria. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.