Large-scale taxonomic profiling of eukaryotic model organisms: a comparison of orthologous proteins encoded by the human, fly, nematode, and yeast genomes

Dopaminergic modulation by quercetin: In silico and in vivo evidence using Caenorhabditis elegans as a model

Quercetin is a flavonol widely distributed in plants and has various described biological functions. Several studies have reported on its ability to restore neuronal function in a wide variety of disease models, including animal models of neurodegenerative disorders such as Parkinson's disease. Quercetin per se can act as a neuroprotector/neuromodulator, especially in diseases related to impaired dopaminergic neurotransmission. However, little is known about how quercetin interacts with the dopaminergic machinery. Here we employed the nematode Caenorhabditis elegans to study this putative interaction. After observing behavioral modulation, mutant analysis and gene expression in C. elegans upon exposure to quercetin at a concentration that does not protect against MPTP, we constructed a homology-based dopamine transporter protein model to conduct a docking study. This led to suggestive evidence on how quercetin may act as a dopaminergic modulator by interacting with C. elegans' dopamine transporter and alter the nematode's exploratory behavior. Consistent with this model, quercetin controls C. elegans behavior in a way dependent on the presence of both the dopamine transporter (dat-1), which is up-regulated upon quercetin exposure, and the dopamine receptor 2 (dop-2), which appears to be mandatory for dat-1 up-regulation. Our data propose an interaction with the dopaminergic machinery that may help to establish the effects of quercetin as a neuromodulator.

OrthoVenn3: an integrated platform for exploring and visualizing orthologous data across genomes

Advancements in comparative genomics research have led to a growing interest in studying species evolution and genetic diversity. To facilitate this research, OrthoVenn3 has been developed as a powerful, web-based tool that enables users to efficiently identify and annotate orthologous clusters and infer phylogenetic relationships across a range of species. The latest upgrade of OrthoVenn includes several important new features, including enhanced orthologous cluster identification accuracy, improved visualization capabilities for numerous sets of data, and wrapped phylogenetic analysis. Furthermore, OrthoVenn3 now provides gene family contraction and expansion analysis to support researchers better understanding the evolutionary history of gene families, as well as collinearity analysis to detect conserved and variable genomic structures. With its intuitive user interface and robust functionality, OrthoVenn3 is a valuable resource for comparative genomics research. The tool is freely accessible at https://orthovenn3.bioinfotoolkits.net.

S-adenosyl-L-homocysteine hydrolase FgSah1 is required for fungal development and virulence in Fusarium graminearum

The S-adenosyl-L-homocysteine hydrolase (Sah1) plays a crucial role in methylation and lipid metabolism in yeast and mammals, yet its function remains elusive in filamentous fungi. In this study, we characterized Sah1 in the phytopathogenic fungus F. graminearum by generating knockout and knockout-complemented strains of FgSAH1. We found that the FgSah1-GFP fusion protein was localized to the cytoplasm, and that deletion of FgSAH1 resulted in defects in vegetative growth, asexual and sexual reproduction, stress responses, virulence, lipid metabolism, and tolerance against fungicides. Moreover, the accumulations of S-adenosyl-L-homocysteine (AdoHcy) and S-adenosyl-L-methionine (AdoMet) (the methyl group donor in most methyl transfer reactions) in ΔFgSah1 were seven- and ninefold higher than those in the wild-type strain, respectively. All of these defective phenotypes in ΔFgSah1 mutants were rescued by target gene complementation. Taken together, these results demonstrate that FgSah1 plays essential roles in methylation metabolism, fungal development, full virulence, multiple stress responses, lipid metabolism, and fungicide sensitivity in F. graminearum. To our knowledge, this is the first report on the systematic functional characterization of Sah1 in F. graminearum.

Transcriptional Profiling of Three Pseudomonas syringae pv. actinidiae Biovars Reveals Different Responses to Apoplast-Like Conditions Related to Strain Virulence on the Host

Pseudomonas syringae pv. actinidiae is a phytopathogen that causes devastating bacterial canker in kiwifruit. Among five biovars defined by genetic, biochemical, and virulence traits, P. syringae pv. actinidiae biovar 3 (Psa3) is the most aggressive and is responsible for the most recent reported outbreaks; however, the molecular basis of its heightened virulence is unclear. Therefore, we designed the first P. syringae multistrain whole-genome microarray, encompassing biovars Psa1, Psa2, and Psa3 and the well-established model P. syringae pv. tomato, and analyzed early bacterial responses to an apoplast-like minimal medium. Transcriptomic profiling revealed i) the strong activation in Psa3 of all hypersensitive reaction and pathogenicity (hrp) and hrp conserved (hrc) cluster genes, encoding components of the type III secretion system required for bacterial pathogenicity and involved in responses to environmental signals; ii) potential repression of the hrp/hrc cluster in Psa2; and iii) activation of flagellum-dependent cell motility and chemotaxis genes in Psa1. The detailed investigation of three gene families encoding upstream regulatory proteins (histidine kinases, their cognate response regulators, and proteins with diguanylate cyclase or phosphodiesterase domains) indicated that cyclic di-GMP may be a key regulator of virulence in P. syringae pv. actinidiae biovars. The gene expression data were supported by the quantification of biofilm formation. Our findings suggest that diverse early responses to the host apoplast, even among bacteria belonging to the same pathovar, can lead to different virulence strategies and may explain the differing outcomes of infections. Based on our detailed structural analysis of hrp operons, we also propose a revision of hrp cluster organization and operon regulation in P. syringae.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Functional and Pathological Roles of AHCY

Adenosylhomocysteinase (AHCY) is a unique enzyme and one of the most conserved proteins in living organisms. AHCY catalyzes the reversible break of S-adenosylhomocysteine (SAH), the by-product and a potent inhibitor of methyltransferases activity. In mammals, AHCY is the only enzyme capable of performing this reaction. Controlled subcellular localization of AHCY is believed to facilitate local transmethylation reactions, by removing excess of SAH. Accordingly, AHCY is recruited to chromatin during replication and active transcription, correlating with increasing demands for DNA, RNA, and histone methylation. AHCY deletion is embryonic lethal in many organisms (from plants to mammals). In humans, AHCY deficiency is associated with an incurable rare recessive disorder in methionine metabolism. In this review, we focus on the AHCY protein from an evolutionary, biochemical, and functional point of view, and we discuss the most recent, relevant, and controversial contributions to the study of this enzyme.

Variation across Species and Levels: Implications for Model Species Research

The selection of model species tends to involve two typically unstated assumptions, namely: (1) that the similarity between species decreases steadily with phylogenetic distance, and (2) that similarities are greater at lower levels of biological organization. The first assumption holds on average, but species similarities tend to decrease with the square root of divergence time, rather than linearly, and lineages with short generation times (which includes most model species) tend to diverge faster than average, making the decrease in similarity non-monotonic. The second assumption is more difficult to test. Comparative molecular research has traditionally emphasized species similarities over differences, whereas comparative research at higher levels of organization frequently highlights the species differences. However, advances in comparative genomics have brought to light a great variety of species differences, not just in gene regulation but also in protein coding genes. Particularly relevant are cases in which homologous high-level characters are based on non-homologous genes. This phenomenon of non-orthologous gene displacement, or "deep non-homology," indicates that species differences at the molecular level can be surprisingly large. Given these observations, it is not surprising that some findings obtained in model species do not generalize across species as well as researchers had hoped, even if the research is molecular.

Homocysteine regulates fatty acid and lipid metabolism in yeast

S-Adenosyl-l-homocysteine hydrolase (AdoHcy hydrolase; Sah1 in yeast/AHCY in mammals) degrades AdoHcy, a by-product and strong product inhibitor of S-adenosyl-l-methionine (AdoMet)-dependent methylation reactions, to adenosine and homocysteine (Hcy). This reaction is reversible, so any elevation of Hcy levels, such as in hyperhomocysteinemia (HHcy), drives the formation of AdoHcy, with detrimental consequences for cellular methylation reactions. HHcy, a pathological condition linked to cardiovascular and neurological disorders, as well as fatty liver among others, is associated with a deregulation of lipid metabolism. Here, we developed a yeast model of HHcy to identify mechanisms that dysregulate lipid metabolism. Hcy supplementation to wildtype cells up-regulated cellular fatty acid and triacylglycerol content and induced a shift in fatty acid composition, similar to changes observed in mutants lacking Sah1. Expression of the irreversible bacterial pathway for AdoHcy degradation in yeast allowed us to dissect the impact of AdoHcy accumulation on lipid metabolism from the impact of elevated Hcy. Expression of this pathway fully suppressed the growth deficit of sah1 mutants as well as the deregulation of lipid metabolism in both the sah1 mutant and Hcy-exposed wildtype, showing that AdoHcy accumulation mediates the deregulation of lipid metabolism in response to elevated Hcy in yeast. Furthermore, Hcy supplementation in yeast led to increased resistance to cerulenin, an inhibitor of fatty acid synthase, as well as to a concomitant decline of condensing enzymes involved in very long-chain fatty acid synthesis, in line with the observed shift in fatty acid content and composition.

Comparative Analysis of the Flavobacterium columnare Genomovar I and II Genomes

Columnaris disease caused by Gram-negative rod Flavobacterium columnare is one of the most common diseases of catfish. F. columnare is also a common problem in other cultured fish species worldwide. F. columnare has three major genomovars; we have sequenced a representative strain from genomovar I (ATCC 49512, which is avirulent in catfish) and genomovar II (94-081, which is highly pathogenic in catfish). Here, we present a comparative analysis of the two genomes. Interestingly, F. columnare ATCC 49512 and 94-081 meet criteria to be considered different species based on the Average Nucleotide Identity (90.71% similar) and DNA–DNA Hybridization (42.6% similar). Genome alignment indicated the two genomes have a large number of rearrangements. However, function-based comparative genomics analysis indicated that the two strains have similar functional capabilities with 2,263 conserved orthologous clusters; strain ATCC 49512 has 290 unique orthologous clusters while strain 94-081 has 391. Both strains carry type I secretion system, type VI secretion system, and type IX secretion system. The two genomes also have similar CRISPR capacities. The F. columnare strain ATCC 49512 genome contains a higher number of insertion sequence families and phage regions, while the F. columnare strain 94-081 genome has more genomic islands and more regulatory gene capacity. Transposon mutagenesis using Tn4351 in pathogenic strain 94-081 yielded six mutants, and experimental infections of fish showed hemolysin and glycine cleavage protein mutants had 15 and 10% mortalities, respectively, while the wild-type strain caused 100% mortalities. Our comparative and mutational analysis yielded important information on classification of genomovars I and II F. columnare as well as potential virulence genes in F. columnare strain 94-081.

Cellular identity at the single-cell level

A single cell creates surprising heterogeneity in a multicellular organism. While every organismal cell shares almost an identical genome, molecular interactions in cells alter the use of DNA sequences to modulate the gene of interest for specialization of cellular functions. Each cell gains a unique identity through molecular coding across the DNA, RNA, and protein conversions. On the other hand, loss of cellular identity leads to critical diseases such as cancer. Most cell identity dissection studies are based on bulk molecular assays that mask differences in individual cells. To probe cell-to-cell variability in a population, we discuss single cell approaches to decode the genetic, epigenetic, transcriptional, and translational mechanisms for cell identity formation. In combination with molecular instructions, the physical principles behind cell identity determination are examined. Deciphering and reprogramming cellular types impact biology and medicine.

OrthoVenn: a web server for genome wide comparison and annotation of orthologous clusters across multiple species

Genome wide analysis of orthologous clusters is an important component of comparative genomics studies. Identifying the overlap among orthologous clusters can enable us to elucidate the function and evolution of proteins across multiple species. Here, we report a web platform named OrthoVenn that is useful for genome wide comparisons and visualization of orthologous clusters. OrthoVenn provides coverage of vertebrates, metazoa, protists, fungi, plants and bacteria for the comparison of orthologous clusters and also supports uploading of customized protein sequences from user-defined species. An interactive Venn diagram, summary counts, and functional summaries of the disjunction and intersection of clusters shared between species are displayed as part of the OrthoVenn result. OrthoVenn also includes in-depth views of the clusters using various sequence analysis tools. Furthermore, OrthoVenn identifies orthologous clusters of single copy genes and allows for a customized search of clusters of specific genes through key words or BLAST. OrthoVenn is an efficient and user-friendly web server freely accessible at http://probes.pw.usda.gov/OrthoVenn or http://aegilops.wheat.ucdavis.edu/OrthoVenn.

Allelic variation within the S-adenosyl-L-homocysteine hydrolase gene family is associated with wood properties in Chinese white poplar (Populus tomentosa)

Background

S-adenosyl-l-homocysteine hydrolase (SAHH) is the only eukaryotic enzyme capable of S-adenosyl-l-homocysteine (SAH) catabolism for the maintenance of cellular transmethylation potential. Recently, biochemical and genetic studies in herbaceous species have obtained important discoveries in the function of SAHH, and an extensive characterization of SAHH family in even one tree species is essential, but currently lacking.

Results

Here, we first identified the SAHH family from Populus tomentosa using molecular cloning method. Phylogenetic analyses of 28 SAHH proteins from dicotyledons, monocotyledons, and lower plants revealed that the sequences formed two monophyletic groups: the PtrSAHHA with PtoSAHHA and PtrSAHHB with PtoSAHHB. Examination of tissue-specific expression profiles of the PtoSAHH family revealed similar expression patterns; high levels of expression in xylem were found. Nucleotide diversity and linkage disequilibrium (LD) in the PtoSAHH family, sampled from P. tomentosa natural distribution, revealed that PtoSAHH harbors high single-nucleotide polymorphism (SNP) diversity (π=0.01059±0.00122 and 0.00930±0.00079,respectively) and low LD (r² > 0.1, within 800 bp and 2,200 bp, respectively). Using an LD-linkage analysis approach, two noncoding SNPs (PtoSAHHB_1065 and PtoSAHHA_2203) and the corresponding haplotypes were found to significantly associate with α-cellulose content, and a nonsynonymous SNP (PtoSAHHB_410) within the SAHH signature motifs showed significant association with fiber length, with an average of 3.14% of the phenotypic variance explained.

Conclusions

The present study demonstrates that PtoSAHHs were split off prior to the divergence of interspecies in Populus, and SAHHs may play a key role promoting transmethylation reactions in the secondary cell walls biosynthesis in trees. Hence, our findings provide insights into SAHH function and evolution in woody species and also offer a theoretical basis for marker-aided selection breeding to improve the wood quality of Populus.

Pleiotropic constraints, expression level, and the evolution of miRNA sequences

Post-transcriptional gene regulation mediated by microRNAs (miRNAs) plays critical roles during development by modulating gene expression and conferring robustness to stochastic errors. Phylogenetic analyses suggest that miRNA acquisition could play a role in phenotypic innovation. Moreover, miRNA-induced regulation strongly impacts genome evolution, increasing selective constraints on 3'UTRs, protein sequences, and expression level divergence. Thus, it is essential to understand the factors governing sequence evolution for this important class of regulatory molecules. Investigation of the patterns of molecular evolution at miRNA loci have been limited in Caenorhabditis elegans because of the lack of a close outgroup. Instead, I used Caenorhabditis briggsae as the focus point of this study because of its close relationship to Caenorhabditis sp. 9. I also corroborated the patterns of sequence evolution in Caenorhabditis using published orthologous relationships among miRNAs in Drosophila. In nematodes and in flies, miRNA sequence divergence is not influenced by the genomic neighborhood (i.e., intronic or intergenic) but is nevertheless affected by the genomic context because X-linked miRNAs evolve faster than autosomal miRNAs. However, this effect of chromosomal linkage can be explained by differential expression levels rather than a fast-X effect. The results presented here support a universal negative relationship between rates of molecular evolution and expression level, and suggest that mutations in highly expressed miRNAs are more likely to be deleterious because they potentially affect a larger number of target genes. Finally, I show that many single family member miRNAs evolve faster than miRNAs from multigene families and have limited functional scope, suggesting that they are not strongly integrated in gene regulatory networks.

Probing the boundaries of orthology: the unanticipated rapid evolution of Drosophila centrosomin

The rapid evolution of essential developmental genes and their protein products is both intriguing and problematic. The rapid evolution of gene products with simple protein folds and a lack of well-characterized functional domains typically result in a low discovery rate of orthologous genes. Additionally, in the absence of orthologs it is difficult to study the processes and mechanisms underlying rapid evolution. In this study, we have investigated the rapid evolution of centrosomin (cnn), an essential gene encoding centrosomal protein isoforms required during syncytial development in Drosophila melanogaster. Until recently the rapid divergence of cnn made identification of orthologs difficult and questionable because Cnn violates many of the assumptions underlying models for protein evolution. To overcome these limitations, we have identified a group of insect orthologs and present conserved features likely to be required for the functions attributed to cnn in D. melanogaster. We also show that the rapid divergence of Cnn isoforms is apparently due to frequent coding sequence indels and an accelerated rate of intronic additions and eliminations. These changes appear to be buffered by multi-exon and multi-reading frame maximum potential ORFs, simple protein folds, and the splicing machinery. These buffering features also occur in other genes in Drosophila and may help prevent potentially deleterious mutations due to indels in genes with large coding exons and exon-dense regions separated by small introns. This work promises to be useful for future investigations of cnn and potentially other rapidly evolving genes and proteins.

Comparative description of ten transcriptomes of newly sequenced invertebrates and efficiency estimation of genomic sampling in non-model taxa

INTRODUCTION

Traditionally, genomic or transcriptomic data have been restricted to a few model or emerging model organisms, and to a handful of species of medical and/or environmental importance. Next-generation sequencing techniques have the capability of yielding massive amounts of gene sequence data for virtually any species at a modest cost. Here we provide a comparative analysis of de novo assembled transcriptomic data for ten non-model species of previously understudied animal taxa.

RESULTS

cDNA libraries of ten species belonging to five animal phyla (2 Annelida [including Sipuncula], 2 Arthropoda, 2 Mollusca, 2 Nemertea, and 2 Porifera) were sequenced in different batches with an Illumina Genome Analyzer II (read length 100 or 150 bp), rendering between ca. 25 and 52 million reads per species. Read thinning, trimming, and de novo assembly were performed under different parameters to optimize output. Between 67,423 and 207,559 contigs were obtained across the ten species, post-optimization. Of those, 9,069 to 25,681 contigs retrieved blast hits against the NCBI non-redundant database, and approximately 50% of these were assigned with Gene Ontology terms, covering all major categories, and with similar percentages in all species. Local blasts against our datasets, using selected genes from major signaling pathways and housekeeping genes, revealed high efficiency in gene recovery compared to available genomes of closely related species. Intriguingly, our transcriptomic datasets detected multiple paralogues in all phyla and in nearly all gene pathways, including housekeeping genes that are traditionally used in phylogenetic applications for their purported single-copy nature.

CONCLUSIONS

We generated the first study of comparative transcriptomics across multiple animal phyla (comparing two species per phylum in most cases), established the first Illumina-based transcriptomic datasets for sponge, nemertean, and sipunculan species, and generated a tractable catalogue of annotated genes (or gene fragments) and protein families for ten newly sequenced non-model organisms, some of commercial importance (i.e., Octopus vulgaris). These comprehensive sets of genes can be readily used for phylogenetic analysis, gene expression profiling, developmental analysis, and can also be a powerful resource for gene discovery. The characterization of the transcriptomes of such a diverse array of animal species permitted the comparison of sequencing depth, functional annotation, and efficiency of genomic sampling using the same pipelines, which proved to be similar for all considered species. In addition, the datasets revealed their potential as a resource for paralogue detection, a recurrent concern in various aspects of biological inquiry, including phylogenetics, molecular evolution, development, and cellular biochemistry.

S-adenosyl-L-homocysteine hydrolase and methylation disorders: yeast as a model system

S-adenosyl-L-methionine (AdoMet)-dependent methylation is central to the regulation of many biological processes: more than 50 AdoMet-dependent methyltransferases methylate a broad spectrum of cellular compounds including nucleic acids, proteins and lipids. Common to all AdoMet-dependent methyltransferase reactions is the release of the strong product inhibitor S-adenosyl-L-homocysteine (AdoHcy), as a by-product of the reaction. S-adenosyl-L-homocysteine hydrolase is the only eukaryotic enzyme capable of reversible AdoHcy hydrolysis to adenosine and homocysteine and, thus, relief from AdoHcy inhibition. Impaired S-adenosyl-L-homocysteine hydrolase activity in humans results in AdoHcy accumulation and severe pathological consequences. Hyperhomocysteinemia, which is characterized by elevated levels of homocysteine in blood, also exhibits a similar phenotype of AdoHcy accumulation due to the reversal of the direction of the S-adenosyl-L-homocysteine hydrolase reaction. Inhibition of S-adenosyl-L-homocysteine hydrolase is also linked to antiviral effects. In this review the advantages of yeast as an experimental system to understand pathologies associated with AdoHcy accumulation will be discussed.

Evolutionary patterns of RNA-based gene duplicates in Caenorhabditis nematodes coincide with their genomic features

BACKGROUND

RNA-based gene duplicates (retrocopies) played pivotal roles in many physiological processes. Nowadays, functional retrocopies have been systematically identified in several mammals, fruit flies, plants, zebrafish and other chordates, etc. However, studies about this kind of duplication in Caenorhabditis nematodes have not been reported.

FINDINGS

We identified 43, 48, 43, 9, and 42 retrocopies, of which 6, 15, 18, 3, and 13 formed chimeric genes in C. brenneri, C. briggsae, C. elegans, C. japonica, and C. remanei, respectively. At least 5 chimeric types exist in Caenorhabditis species, of which retrocopy recruiting both N and C terminus is the commonest one. Evidences from different analyses demonstrate many retrocopies and almost all chimeric genes may be functional in these species. About half of retrocopies in each species has coordinates in other species, and we suggest that retrocopies in closely related species may be helpful in identifying retrocopies for one certain species.

CONCLUSIONS

A number of retrocopies and chimeric genes exist in Caenorhabditis genomes, and some of them may be functional. The evolutionary patterns of these genes may correlate with their genomic features, such as the activity of retroelements, the high rate of mutation and deletion rate, and a large proportion of genes subject to trans-splicing.

Genome-wide mapping of histone modifications and mass spectrometry reveal H4 acetylation bias and H3K36 methylation at gene promoters in fission yeast

AIMS

To map histone modifications with unprecedented resolution both globally and locus-specifically, and to link modification patterns to gene expression.

MATERIALS & METHODS

Using correlations between quantitative mass spectrometry and chromatin immunoprecipitation/microarray analyses, we have mapped histone post-translational modifications in fission yeast (Schizosaccharomyces pombe).

RESULTS

Acetylations at lysine 9, 18 and 27 of histone H3 give the best positive correlations with gene expression in this organism. Using clustering analysis and gene ontology search tools, we identified promoter histone modification patterns that characterize several classes of gene function. For example, gene promoters of genes involved in cytokinesis have high H3K36me2 and low H3K4me2, whereas the converse pattern is found ar promoters of gene involved in positive regulation of the cell cycle. We detected acetylation of H4 preferentially at lysine 16 followed by lysine 12, 8 and 5. Our analysis shows that this H4 acetylation bias in the coding regions is dependent upon gene length and linked to gene expression. Our analysis also reveals a role for H3K36 methylation at gene promoters where it functions in a crosstalk between the histone methyltransferase Set2(KMT3) and the histone deacetylase Clr6, which removes H3K27ac leading to repression of transcription.

CONCLUSION

Histone modification patterns could be linked to gene expression in fission yeast.

The draft genome of the parasitic nematode Trichinella spiralis

Genome-based studies of metazoan evolution are most informative when phylogenetically diverse species are incorporated in the analysis. As such, evolutionary trends within and outside the phylum Nematoda have been less revealing by focusing only on comparisons involving Caenorhabditis elegans. Herein, we present a draft of the 64 megabase nuclear genome of Trichinella spiralis, containing 15,808 protein coding genes. This parasitic nematode is an extant member of a clade that diverged early in the evolution of the phylum enabling identification of archetypical genes and molecular signatures exclusive to nematodes. Comparative analyses support intrachromosomal rearrangements across the phylum, disproportionate numbers of protein family deaths over births in parasitic vs. a non-parasitic nematode, and a preponderance of gene loss and gain events in nematodes relative to Drosophila melanogaster. This sequence and the panphylum characteristics identified herein will advance evolutionary studies and strategies to combat global parasites of humans, food animals and crops.

Characterization of the deleted in autism 1 protein family: implications for studying cognitive disorders

Autism spectrum disorders (ASDs) are a group of commonly occurring, highly-heritable developmental disabilities. Human genes c3orf58 or Deleted In Autism-1 (DIA1) and cXorf36 or Deleted in Autism-1 Related (DIA1R) are implicated in ASD and mental retardation. Both gene products encode signal peptides for targeting to the secretory pathway. As evolutionary medicine has emerged as a key tool for understanding increasing numbers of human diseases, we have used an evolutionary approach to study DIA1 and DIA1R. We found DIA1 conserved from cnidarians to humans, indicating DIA1 evolution coincided with the development of the first primitive synapses. Nematodes lack a DIA1 homologue, indicating Caenorhabditis elegans is not suitable for studying all aspects of ASD etiology, while zebrafish encode two DIA1 paralogues. By contrast to DIA1, DIA1R was found exclusively in vertebrates, with an origin coinciding with the whole-genome duplication events occurring early in the vertebrate lineage, and the evolution of the more complex vertebrate nervous system. Strikingly, DIA1R was present in schooling fish but absent in fish that have adopted a more solitary lifestyle. An additional DIA1-related gene we named DIA1-Like (DIA1L), lacks a signal peptide and is restricted to the genomes of the echinoderm Strongylocentrotus purpuratus and cephalochordate Branchiostoma floridae. Evidence for remarkable DIA1L gene expansion was found in B. floridae. Amino acid alignments of DIA1 family gene products revealed a potential Golgi-retention motif and a number of conserved motifs with unknown function. Furthermore, a glycine and three cysteine residues were absolutely conserved in all DIA1-family proteins, indicating a critical role in protein structure and/or function. We have therefore identified a new metazoan protein family, the DIA1-family, and understanding the biological roles of DIA1-family members will have implications for our understanding of autism and mental retardation.

Comparing biological networks via graph compression

BACKGROUND

Comparison of various kinds of biological data is one of the main problems in bioinformatics and systems biology. Data compression methods have been applied to comparison of large sequence data and protein structure data. Since it is still difficult to compare global structures of large biological networks, it is reasonable to try to apply data compression methods to comparison of biological networks. In existing compression methods, the uniqueness of compression results is not guaranteed because there is some ambiguity in selection of overlapping edges.

RESULTS

This paper proposes novel efficient methods, CompressEdge and CompressVertices, for comparing large biological networks. In the proposed methods, an original network structure is compressed by iteratively contracting identical edges and sets of connected edges. Then, the similarity of two networks is measured by a compression ratio of the concatenated networks. The proposed methods are applied to comparison of metabolic networks of several organisms, H. sapiens, M. musculus, A. thaliana, D. melanogaster, C. elegans, E. coli, S. cerevisiae, and B. subtilis, and are compared with an existing method. These results suggest that our methods can efficiently measure the similarities between metabolic networks.

CONCLUSIONS

Our proposed algorithms, which compress node-labeled networks, are useful for measuring the similarity of large biological networks.

Characterizing Ancylostoma caninum transcriptome and exploring nematode parasitic adaptation

BACKGROUND

Hookworm infection is one of the most important neglected diseases in developing countries, with approximately 1 billion people infected worldwide. To better understand hookworm biology and nematode parasitism, the present study generated a near complete transcriptome of the canine hookworm Ancylostoma caninum to a very high coverage using high throughput technology, and compared it to those of the free-living nematode Caenorhabditis elegans and the parasite Brugia malayi.

RESULTS

The generated transcripts from four developmental stages, infective L3, serum stimulated L3, adult male and adult female, covered 93% of the A. caninum transcriptome. The broad diversity among nematode transcriptomes was confirmed, and an impact of parasitic adaptation on transcriptome diversity was inferred. Intra-population analysis showed that A. caninum has higher coding sequence diversity than humans. Examining the developmental expression profiles of A. caninum revealed major transitions in gene expression from larval stages to adult. Adult males expressed the highest number of selectively expressed genes, but adult female expressed the highest number of selective parasitism-related genes. Genes related to parasitism adaptation and A. caninum specific genes exhibited more expression selectivity while those conserved in nematodes tend to be consistently expressed. Parasitism related genes were expressed more selectively in adult male and female worms. The comprehensive analysis of digital expression profiles along with transcriptome comparisons enabled identification of a set of parasitism genes encoding secretory proteins in animal parasitic nematode.

CONCLUSIONS

This study validated the usage of deep sequencing for gene expression profiling. Parasitic adaptation of the canine hookworm is related to its diversity and developmental dynamics. This comprehensive comparative genomic and expression study substantially improves our understanding of the basic biology and parasitism of hookworms and, is expected, in the long run, to accelerate research toward development of vaccines and novel anthelmintics.

The nematode story: Hox gene loss and rapid evolution

The loss in some taxa of conserved developmental control genes that are present in the vast majority of animal lineages is an understudied phenomenon. It is likely that in those lineages in which loss has occurred it may be a strong signal of the mode, tempo and direction of developmental evolution and thus identify ways of generating morphological novelties. Intuitively we might expect these novelties to be particularly those associated with morphological simplifications. One striking example of this has occurred within the nematodes. It appears that over half the ancestral bilaterian Hox cluster has been lost from the model organism Caenorhabditis elegans and its closest related species. Studying the Hox gene complement of nematodes across the phylum has shown that many, if not all these losses occurred within the phylum. Other nematode clades only distantly related to C. elegans have additional Hox genes orthologous to those present in the ancestral bilaterian but absent from the model nematode. In some of these cases rapid sequence evolution of the homeodomain itself obscures orthology assignment until comparison is made with sequences from multiple nematode clades with slower evolving Hox genes. Across the phylum the homeodomains of the Hox genes that are present are evolving very rapidly. In one particular case the genomic arrangement of two homeodomains suggests a mechanism for gene loss. Studying the function in nematodes of the Hox genes absent from C. elegans awaits further research and the establishment of new nematode models. However, what we do know about Hox gene functions suggests that the genetic circuits within which Hox genes act have changed significantly within C. elegans and its close relatives.

Hazardous effect of tannery solid waste leachates on development and reproduction in Drosophila melanogaster: 70kDa heat shock protein as a marker of cellular damage

Rapid industrialization has increased the burden of chemicals in the environment. These chemicals may be harmful to development and reproduction of any organism. We therefore analyzed the adverse effects of leachates from a tannery solid waste on development and reproduction using Drosophila. We show a significant delay in mean emergence of flies observed at the higher concentrations of the leachates, indicating their effect on the organism's development. Significant leachate-induced effect on reproduction of the organism was also observed. Sub-organismal analyses revealed Hsp70 expression and tissue damage in a sex-specific manner. Refractoriness of Hsp70 expression in accessory glands of male flies and ovaries of females was concurrent with tissue damage. Genes encoding certain seminal proteins (Acp70A and Acp36DE) from accessory glands were significantly down-regulated at higher concentrations of the leachates. The study suggests that (i) sub-organismal adverse responses are reflected at organismal level, (ii) tannery waste leachates cause adverse effects on the expression of genes encoding seminal proteins that facilitate normal reproduction and (iii) Hsp70 may be used as a marker of cellular damage for reproductive organs.

OrthoSelect: a protocol for selecting orthologous groups in phylogenomics

BACKGROUND

Phylogenetic studies using expressed sequence tags (EST) are becoming a standard approach to answer evolutionary questions. Such studies are usually based on large sets of newly generated, unannotated, and error-prone EST sequences from different species. A first crucial step in EST-based phylogeny reconstruction is to identify groups of orthologous sequences. From these data sets, appropriate target genes are selected, and redundant sequences are eliminated to obtain suitable sequence sets as input data for tree-reconstruction software. Generating such data sets manually can be very time consuming. Thus, software tools are needed that carry out these steps automatically.

RESULTS

We developed a flexible and user-friendly software pipeline, running on desktop machines or computer clusters, that constructs data sets for phylogenomic analyses. It automatically searches assembled EST sequences against databases of orthologous groups (OG), assigns ESTs to these predefined OGs, translates the sequences into proteins, eliminates redundant sequences assigned to the same OG, creates multiple sequence alignments of identified orthologous sequences and offers the possibility to further process this alignment in a last step by excluding potentially homoplastic sites and selecting sufficiently conserved parts. Our software pipeline can be used as it is, but it can also be adapted by integrating additional external programs. This makes the pipeline useful for non-bioinformaticians as well as to bioinformatic experts. The software pipeline is especially designed for ESTs, but it can also handle protein sequences.

CONCLUSION

OrthoSelect is a tool that produces orthologous gene alignments from assembled ESTs. Our tests show that OrthoSelect detects orthologs in EST libraries with high accuracy. In the absence of a gold standard for orthology prediction, we compared predictions by OrthoSelect to a manually created and published phylogenomic data set. Our tool was not only able to rebuild the data set with a specificity of 98%, but it detected four percent more orthologous sequences. Furthermore, the results OrthoSelect produces are in absolut agreement with the results of other programs, but our tool offers a significant speedup and additional functionality, e.g. handling of ESTs, computing sequence alignments, and refining them. To our knowledge, there is currently no fully automated and freely available tool for this purpose. Thus, OrthoSelect is a valuable tool for researchers in the field of phylogenomics who deal with large quantities of EST sequences. OrthoSelect is written in Perl and runs on Linux/Mac OS X. The tool can be downloaded at (http://gobics.de/fabian/orthoselect.php).

Bilaterian phylogeny: a broad sampling of 13 nuclear genes provides a new Lophotrochozoa phylogeny and supports a paraphyletic basal acoelomorpha

During the past decade, great progress has been made in clarifying the relationships among bilaterian animals. Studies based on a limited number of markers established new hypotheses such as the existence of three superclades (Deuterostomia, Ecdysozoa, and Lophotrochozoa) but left major questions unresolved. The data sets used to the present either bear few characters for many taxa (i.e., the ribosomal genes) or present many characters but lack many phyla (such as recent phylogenomic approaches) failing to provide definitive answers for all the regions of the bilaterian tree. We performed phylogenetic analyses using a molecular matrix with a high number of characters and bilaterian phyla. This data set is built from 13 genes (8,880 bp) belonging to 90 taxa from 27 bilaterian phyla. Probabilistic analyses robustly support the three superclades, the monophyly of Chordata, a spiralian clade including Brachiozoa, the basal position of a paraphyletic Acoelomorpha, and point to an ecdysozoan affiliation for Chaetognatha. This new phylogeny not only agrees with most classical molecular results but also provides new insights into the relationships between lophotrochozoans and challenges the results obtained using high-throughput strategies, highlighting the problems associated with the current trend to increase gene number rather than taxa.

Phylogenetic-signal dissection of nuclear housekeeping genes supports the paraphyly of sponges and the monophyly of Eumetazoa

The relationships at the base of the metazoan tree have been difficult to robustly resolve, and there are several different hypotheses regarding the interrelationships among sponges, cnidarians, ctenophores, placozoans, and bilaterians, with each hypothesis having different implications for the body plan of the last common ancestor of animals and the paleoecology of the late Precambrian. We have sequenced seven nuclear housekeeping genes from 17 new sponges, bringing the total to 29 species analyzed, including multiple representatives of the Demospongiae, Calcarea, Hexactinellida, and Homoscleromorpha, and analyzed a data set also including six nonmetazoan outgroups and 36 eumetazoans using a variety of phylogenetic methods and evolutionary models. We used leaf stability to identify rogue taxa and investigate their effect on the support of the nodes in our trees, and we identified clades most likely to represent phylogenetic artifacts through the comparison of trees derived using different methods (and models) and through site-stripping analyses. Further, we investigated compositional heterogeneity and tested whether amino acid composition bias affected our results. Finally, we used Bayes factors to compare our results against previously published phylogenies. All our maximum likelihood (ML) and Bayesian analyses find sponges to be paraphyletic, with all analyses finding three extant paraphyletic sponge lineages, Demospongiae plus Hexactinellida, Calcarea, and Homoscleromorpha. All but one of our ML and Bayesian analyses support the monophyly of Eumetazoa (here Cnidaria + Bilateria) and a sister group relationship between Placozoa (here Trichoplax adhaerens) and Eumetazoa. Bayes factors invariably provide decisive support in favor of poriferan paraphyly when compared against either a sister group relationship between Porifera and Cnidaria or with a monophyletic Porifera with respect to a monophyletic Eumetazoa. Although we were able to recover sponge monophyly using our data set, this was only possible under unrealistic evolutionary models, if poorly performing phylogenetic methods were used, or in situations where the potential for the generation of tree reconstruction artifacts was artificially exacerbated. Everything considered, our data set does not provide any support for a monophyletic Diploblastica (here Placozoa + Cnidaria + Porifera) and suggests that a monophyletic Porifera may be better seen as a phylogenetic artifact.

What was the set of ubiquitin and ubiquitin-like conjugating enzymes in the eukaryote common ancestor?

Ubiquitin (Ub)-conjugating enzymes (E2) are key enzymes in ubiquitination or Ub-like modifications of proteins. We searched for all proteins belonging to the E2 enzyme super-family in seven species (Homo sapiens, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans, Schizosaccharomyces pombe, Saccharomyces cerevisiae, and Arabidopsis thaliana) to identify families and to reconstruct each family’s phylogeny. Our phylogenetic analysis of 207 genes led us to define 17 E2 families, with 37 E2 genes, in the human genome. The subdivision of E2 into four classes did not correspond to the phylogenetic tree. The sequence signature HPN (histidine–proline–asparagine), followed by a tryptophan residue at 16 (up to 29) amino acids, was highly conserved. When present, the active cysteine was found 7 to 8 amino acids from the C-terminal end of HPN. The secondary structures were characterized by a canonical alpha/beta fold. Only family 10 deviated from the common organization because the proteins were devoid of enzymatic activity. Family 7 had an insertion between beta strands 1 and 2; families 3, 5 and 14 had an insertion between the active cysteine and the conserved tryptophan. The three-dimensional data of these proteins highlight a strong structural conservation of the core domain. Our analysis shows that the primitive eukaryote ancestor possessed a diversified set of E2 enzymes, thus emphasizing the importance of the Ub pathway. This comprehensive overview of E2 enzymes emphasizes the diversity and evolution of this superfamily and helps clarify the nomenclature and true orthologies. A better understanding of the functions of these enzymes is necessary to decipher several human diseases.

Specific retention of the protostome-specific PsGEF may parallel with the evolution of mushroom bodies in insect and lophotrochozoan brains

Background

Gene gain and subsequent retention or loss during evolution may be one of the underlying mechanisms involved in generating the diversity of metazoan nervous systems. However, the causal relationships acting therein have not been studied extensively.

Results

We identified the gene PsGEF (protostome-specific GEF), which is present in all the sequenced genomes of insects and limpet but absent in those of sea anemones, deuterostomes, and nematodes. In Drosophila melanogaster, PsGEF encodes a short version of a protein with the C2 and PDZ domains, as well as a long version with the C2, PDZ, and RhoGEF domains through alternative splicing. Intriguingly, the exons encoding the RhoGEF domain are specifically deleted in the Daphnia pulex genome, suggesting that Daphnia PsGEF contains only the C2 and PDZ domains. Thus, the distribution of PsGEF containing the C2, PDZ, and RhoGEF domains among metazoans appears to coincide with the presence of mushroom bodies. Mushroom bodies are prominent neuropils involved in the processing of multiple sensory inputs as well as associative learning in the insect, platyhelminth, and annelid brains. In the adult Drosophila brain, PsGEF is expressed in mushroom bodies, antennal lobe, and optic lobe, where it is necessary for the correct axon branch formation of alpha/beta neurons in mushroom bodies. PsGEF genetically interacts with Rac1 but not other Rho family members, and the RhoGEF domain of PsGEF induces actin polymerization in the membrane, thus resulting in the membrane ruffling that is observed in cultured cells with activated forms of Rac.

Conclusion

The specific acquisition of PsGEF by the last common ancestor of protostomes followed by its retention or loss in specific animal species during evolution demonstrates that there are some structural and/or functional features common between insect and lophotrochozoan nervous systems (for example, mushroom bodies), which are absent in all deuterostomes and cnidarians. PsGEF is therefore one of genes associated with the diversity of metazoan nervous systems.

A novel hybrid yeast-human network analysis reveals an essential role for FNBP1L in antibacterial autophagy

Autophagy is a conserved cellular process required for the removal of defective organelles, protein aggregates, and intracellular pathogens. We used a network analysis strategy to identify novel human autophagy components based upon the yeast interactome centered on the core yeast autophagy proteins. This revealed the potential involvement of 14 novel mammalian genes in autophagy, several of which have known or predicted roles in membrane organization or dynamics. We selected one of these membrane interactors, FNBP1L (formin binding protein 1-like), an F-BAR-containing protein (also termed Toca-1), for further study based upon a predicted interaction with ATG3. We confirmed the FNBP1L/ATG3 interaction biochemically and mapped the FNBP1L domains responsible. Using a functional RNA interference approach, we determined that FNBP1L is essential for autophagy of the intracellular pathogen Salmonella enterica serovar Typhimurium and show that the autophagy process serves to restrict the growth of intracellular bacteria. However, FNBP1L appears dispensable for other forms of autophagy induced by serum starvation or rapamycin. We present a model where FNBP1L is essential for autophagy of intracellular pathogens and identify FNBP1L as a differentially used molecule in specific autophagic contexts. By using network biology to derive functional biological information, we demonstrate the utility of integrated genomics to novel molecule discovery in autophagy.

Expression profiles of urbilaterian genes uniquely shared between honey bee and vertebrates

BACKGROUND

Large-scale comparison of metazoan genomes has revealed that a significant fraction of genes of the last common ancestor of Bilateria (Urbilateria) is lost in each animal lineage. This event could be one of the underlying mechanisms involved in generating metazoan diversity. However, the present functions of these ancient genes have not been addressed extensively. To understand the functions and evolutionary mechanisms of such ancient Urbilaterian genes, we carried out comprehensive expression profile analysis of genes shared between vertebrates and honey bees but not with the other sequenced ecdysozoan genomes (honey bee-vertebrate specific, HVS genes) as a model.

RESULTS

We identified 30 honey bee and 55 mouse HVS genes. Many HVS genes exhibited tissue-selective expression patterns; intriguingly, the expression of 60% of honey bee HVS genes was found to be brain enriched, and 24% of mouse HVS genes were highly expressed in either or both the brain and testis. Moreover, a minimum of 38% of mouse HVS genes demonstrated neuron-enriched expression patterns, and 62% of them exhibited expression in selective brain areas, particularly the forebrain and cerebellum. Furthermore, gene ontology (GO) analysis of HVS genes predicted that 35% of genes are associated with DNA transcription and RNA processing.

CONCLUSION

These results suggest that HVS genes include genes that are biased towards expression in the brain and gonads. They also demonstrate that at least some of Urbilaterian genes retained in the specific animal lineage may be selectively maintained to support the species-specific phenotypes.

S-adenosyl-L-homocysteine hydrolase, key enzyme of methylation metabolism, regulates phosphatidylcholine synthesis and triacylglycerol homeostasis in yeast: implications for homocysteine as a risk factor of atherosclerosis

In eukaryotes, S-adenosyl-L-homocysteine hydrolase (Sah1) offers a single way for degradation of S-adenosyl-L-homocysteine, a product and potent competitive inhibitor of S-adenosyl-L-methionine (AdoMet)-dependent methyltransferases. De novo phosphatidylcholine (PC) synthesis requires three AdoMet-dependent methylation steps. Here we show that down-regulation of SAH1 expression in yeast leads to accumulation of S-adenosyl-L-homocysteine and decreased de novo PC synthesis in vivo. This decrease is accompanied by an increase in triacylglycerol (TG) levels, demonstrating that Sah1-regulated methylation has a major impact on cellular lipid homeostasis. TG accumulation is also observed in cho2 and opi3 mutants defective in methylation of phosphatidylethanolamine to PC, confirming that PC de novo synthesis and TG synthesis are metabolically coupled through the efficiency of the phospholipid methylation reaction. Indeed, because both types of lipids share phosphatidic acid as a precursor, we find in cells with down-regulated Sah1 activity major alterations in the expression of the INO1 gene as well as in the localization of Opi1, a negative regulatory factor of phospholipid synthesis, which binds and is retained in the endoplasmic reticulum membrane by phosphatidic acid in conjunction with VAMP/synaptobrevin-associated protein, Scs2. The addition of homocysteine, by the reversal of the Sah1-catalyzed reaction, also leads to TG accumulation in yeast, providing an attractive model for the role of homocysteine as a risk factor of atherosclerosis in humans.

Clustering orthologous proteins across phylogenetically distant species

The quality of orthologous protein clusters (OPCs) is largely dependent on the results of the reciprocal BLAST (basic local alignment search tool) hits among genomes. The BLAST algorithm is very efficient and fast, but it is very difficult to get optimal solution among phylogenetically distant species because the genomes with large evolutionary distance typically have low similarity in their protein sequences. To reduce the false positives in the OPCs, thresholding is often employed on the BLAST scores. However, the thresholding also eliminates large numbers of true positives as the orthologs from distant species likely have low BLAST scores. To rectify this problem, we introduce a new hybrid method combining the Recursive and the Markov CLuster (MCL) algorithms without using the BLAST thresholding. In the first step, we use InParanoid to produce n(n-1)/2 ortholog tables from n genomes. After combining all the tables into one, our clustering algorithm clusters ortholog pairs recursively in the table. Then, our method employs MCL algorithm to compute the clusters and refines the clusters by adjusting the inflation factor. We tested our method using six different genomes and evaluated the results by comparing against Kegg Orthology (KO) OPCs, which are generated from manually curated pathways. To quantify the accuracy of the results, we introduced a new intuitive similarity measure based on our Least-move algorithm that computes the consistency between two OPCs. We compared the resulting OPCs with the KO OPCs using this measure. We also evaluated the performance of our method using InParanoid as the baseline approach. The experimental results show that, at the inflation factor 1.3, we produced 54% more orthologs than InParanoid sacrificing a little less accuracy (1.7% less) than InParanoid, and at the factor 1.4, produced not only 15% more orthologs than InParanoid but also a higher accuracy (1.4% more) than InParanoid.

Ancient origin of the Parkinson disease gene LRRK2

Dominant mutations in the LRRK2 gene, a member of the Roco family, cause both familial and sporadic Parkinson disease. LRRK genes had so far been detected only in bilaterian animals. In deuterostomes, including humans, two LRRK genes (LRRK1 and LRRK2) exist, while in protostomes a single LRRK gene has been found. In this study, I combine structural and phylogenetic analyses to show that the cnidarian Nematostella vectensis has four LRRK genes. One of them is a bona fide orthologue of the human LRRK2 gene, demonstrating that this gene has an ancient origin. Two others are, respectively, orthologues of the deuterostome LRRK1 and the protostome LRRK genes. The fourth gene is probably cnidarian-specific. This precise characterization of the early evolution of LRRK genes in animals has important implications, because it indicates that the Drosophila and Caenorhabditis LRRK genes, which are studied to gain an understanding of LRRK2 function, are not true orthologues of the human Parkinson disease gene. Novel functional insights are also gained by comparison of the structures of LRRK2 genes in distantly related species.

Homoplasy in genome-wide analysis of rare amino acid replacements: the molecular-evolutionary basis for Vavilov's law of homologous series

BACKGROUND

Rare genomic changes (RGCs) that are thought to comprise derived shared characters of individual clades are becoming an increasingly important class of markers in genome-wide phylogenetic studies. Recently, we proposed a new type of RGCs designated RGC_CAMs (after Conserved Amino acids-Multiple substitutions) that were inferred using genome-wide identification of amino acid replacements that were: i) located in unambiguously aligned regions of orthologous genes, ii) shared by two or more taxa in positions that contain a different, conserved amino acid in a much broader range of taxa, and iii) require two or three nucleotide substitutions. When applied to animal phylogeny, the RGC_CAM approach supported the coelomate clade that unites deuterostomes with arthropods as opposed to the ecdysozoan (molting animals) clade. However, a non-negligible level of homoplasy was detected.

RESULTS

We provide a direct estimate of the level of homoplasy caused by parallel changes and reversals among the RGC_CAMs using 462 alignments of orthologous genes from 19 eukaryotic species. It is shown that the impact of parallel changes and reversals on the results of phylogenetic inference using RGC_CAMs cannot explain the observed support for the Coelomata clade. In contrast, the evidence in support of the Ecdysozoa clade, in large part, can be attributed to parallel changes. It is demonstrated that parallel changes are significantly more common in internal branches of different subtrees that are separated from the respective common ancestor by relatively short times than in terminal branches separated by longer time intervals. A similar but much weaker trend was detected for reversals. The observed evolutionary trend of parallel changes is explained in terms of the covarion model of molecular evolution. As the overlap between the covarion sets in orthologous genes from different lineages decreases with time after divergence, the likelihood of parallel changes decreases as well.

CONCLUSION

The level of homoplasy observed here appears to be low enough to justify the utility of RGC_CAMs and other types of RGCs for resolution of hard problems in phylogeny. Parallel changes, one of the major classes of events leading to homoplasy, occur much more often in relatively recently diverged lineages than in those separated from their last common ancestor by longer time intervals of time. This pattern seems to provide the molecular-evolutionary underpinning of Vavilov's law of homologous series and is readily interpreted within the framework of the covarion model of molecular evolution.

Comparative analysis of in vivo interactions between Rev1 protein and other Y-family DNA polymerases in animals and yeasts

Eukaryotes are endowed with multiple specialized DNA polymerases, some (if not all) of which are believed to play important roles in the tolerance of base damage during DNA replication. Among these DNA polymerases, Rev1 protein (a deoxycytidyl transferase) from vertebrates interacts with several other specialized polymerases via a highly conserved C-terminal region. The present studies assessed whether these interactions are retained in more experimentally tractable model systems, including yeasts, flies, and the nematode C. elegans. We observed a physical interaction between Rev1 protein and other Y-family polymerases in the fruit fly Drosophila melanogaster. However, despite the fact that the C-terminal region of Drosophila and yeast Rev1 are conserved from vertebrates to a similar extent, such interactions were not observed in Saccharomyces cerevisiae or Schizosaccharomyces pombe. With respect to regions in specialized DNA polymerases that are required for interaction with Rev1, we find predicted disorder to be an underlying structural commonality. The results of this study suggest that special consideration should be exercised when making mechanistic extrapolations regarding translesion DNA synthesis from one eukaryotic system to another.

The tardigrade Hypsibius dujardini, a new model for studying the evolution of development

Studying development in diverse taxa can address a central issue in evolutionary biology: how morphological diversity arises through the evolution of developmental mechanisms. Two of the best-studied developmental model organisms, the arthropod Drosophila and the nematode Caenorhabditis elegans, have been found to belong to a single protostome superclade, the Ecdysozoa. This finding suggests that a closely related ecdysozoan phylum could serve as a valuable model for studying how developmental mechanisms evolve in ways that can produce diverse body plans. Tardigrades, also called water bears, make up a phylum of microscopic ecdysozoan animals. Tardigrades share many characteristics with C. elegans and Drosophila that could make them useful laboratory models, but long-term culturing of tardigrades historically has been a challenge, and there have been few studies of tardigrade development. Here, we show that the tardigrade Hypsibius dujardini can be cultured continuously for decades and can be cryopreserved. We report that H. dujardini has a compact genome, a little smaller than that of C. elegans or Drosophila, and that sequence evolution has occurred at a typical rate. H. dujardini has a short generation time, 13-14 days at room temperature. We have found that the embryos of H. dujardini have a stereotyped cleavage pattern with asymmetric cell divisions, nuclear migrations, and cell migrations occurring in reproducible patterns. We present a cell lineage of the early embryo and an embryonic staging series. We expect that these data can serve as a platform for using H. dujardini as a model for studying the evolution of developmental mechanisms.

Support for the Coelomata clade of animals from a rigorous analysis of the pattern of intron conservation

Many intron positions are conserved in varying subsets of eukaryotic genomes and, consequently, comprise a potentially informative class of phylogenetic characters. Roy and Gilbert developed a method of phylogenetic reconstruction using the patterns of intron presence-absence in eukaryotic genes and, applying this method to the analysis of animal phylogeny, obtained support for an Ecdysozoa clade (Roy SW, Gilbert W. 2005. Resolution of a deep animal divergence by the pattern of intron conservation. Proc Natl Acad Sci USA. 102:4403-4408). The critical assumption in the method was the independence of intron loss in different branches of the phylogenetic tree. Here, this assumption is refuted by showing that the branch-specific intron loss rates are strongly correlated. We show that different tree topologies are obtained, in each case with a significant statistical support, when different subsets of intron positions are analyzed. The analysis of the conserved intron positions supports the Coelomata topology, that is, a clade comprised of arthropods and chordates, whereas the analysis of more variable intron positions favors the Ecdysozoa topology, that is, a clade of arthropods and nematodes. We show, however, that the support for Ecdysozoa is fully explained by parallel loss of introns in nematodes and arthropods, a factor that does not contribute to the analysis of the conserved introns. The developed procedure for the identification and analysis of conserved introns and other characters with minimal or no homoplasy is expected to be useful for resolving many hard phylogenetic problems.

Dating the early evolution of plants: detection and molecular clock analyses of orthologs

Orthologs generally are under selective pressure against loss of function, while paralogs usually accumulate mutations and finally die or deviate in terms of function or regulation. Most ortholog detection methods contaminate the resulting datasets with a substantial amount of paralogs. Therefore we aimed to implement a straightforward method that allows the detection of ortholog clusters with a reduced amount of paralogs from completely sequenced genomes. The described cross-species expansion of the reciprocal best BLAST hit method is a time-effective method for ortholog detection, which results in 68% truly orthologous clusters and the procedure specifically enriches single-copy orthologs. The detection of true orthologs can provide a phylogenetic toolkit to better understand evolutionary processes. In a study across six photosynthetic eukaryotes, nuclear genes of putative mitochondrial origin were shown to be over-represented among single copy orthologs. These orthologs are involved in fundamental biological processes like amino acid metabolism or translation. Molecular clock analyses based on this dataset yielded divergence time estimates for the red/green algae (1,142 MYA), green algae/land plant (725 MYA), mosses/seed plant (496 MYA), gymno-/angiosperm (385 MYA) and monocotyledons/core eudicotyledons (301 MYA) divergence times.

Evolution of secretin family GPCR members in the metazoa

Background

Comparative approaches using protostome and deuterostome data have greatly contributed to understanding gene function and organismal complexity. The family 2 G-protein coupled receptors (GPCRs) are one of the largest and best studied hormone and neuropeptide receptor families. They are suggested to have arisen from a single ancestral gene via duplication events. Despite the recent identification of receptor members in protostome and early deuterostome genomes, relatively little is known about their function or origin during metazoan divergence. In this study a comprehensive description of family 2 GPCR evolution is given based on in silico and expression analyses of the invertebrate receptor genes.

Results

Family 2 GPCR members were identified in the invertebrate genomes of the nematodes C. elegans and C. briggsae, the arthropods D. melanogaster and A. gambiae (mosquito) and in the tunicate C. intestinalis. This suggests that they are of ancient origin and have evolved through gene/genome duplication events. Sequence comparisons and phylogenetic analyses have demonstrated that the immediate gene environment, with regard to gene content, is conserved between the protostome and deuterostome receptor genomic regions. Also that the protostome genes are more like the deuterostome Corticotrophin Releasing Factor (CRF) and Calcitonin/Calcitonin Gene-Related Peptide (CAL/CGRP) receptors members than the other family 2 GPCR members. The evolution of family 2 GPCRs in deuterostomes is characterised by acquisition of new family members, with SCT (Secretin) receptors only present in tetrapods. Gene structure is characterised by an increase in intron number with organismal complexity with the exception of the vertebrate CAL/CGRP receptors.

Conclusion

The family 2 GPCR members provide a good example of gene duplication events occurring in tandem with increasing organismal complexity during metazoan evolution. The putative ancestral receptors are proposed to be more like the deuterostome CAL/CGRP and CRF receptors and this may be associated with their fundamental role in calcium regulation and the stress response, both of which are essential for survival.

Mitogenomics and phylogenomics reveal priapulid worms as extant models of the ancestral Ecdysozoan

Research into arthropod evolution is hampered by the derived nature and rapid evolution of the best-studied out-group: the nematodes. We consider priapulids as an alternative out-group. Priapulids are a small phylum of bottom-dwelling marine worms; their tubular body with spiny proboscis or introvert has changed little over 520 million years and recognizable priapulids are common among exceptionally preserved Cambrian fossils. Using the complete mitochondrial genome and 42 nuclear genes from Priapulus caudatus, we show that priapulids are slowly evolving ecdysozoans; almost all these priapulid genes have evolved more slowly than nematode orthologs and the priapulid mitochondrial gene order may be unchanged since the Cambrian. Considering their primitive bodyplan and embryology and the great conservation of both nuclear and mitochondrial genomes, priapulids may deserve the popular epithet of "living fossil." Their study is likely to yield significant new insights into the early evolution of the Ecdysozoa and the origins of the arthropods and their kin as well as aiding inference of the morphology of ancestral Ecdysozoa and Bilateria and their genomes.

Further use of nearly complete 28S and 18S rRNA genes to classify Ecdysozoa: 37 more arthropods and a kinorhynch

This work expands on a study from 2004 by Mallatt, Garey, and Shultz [Mallatt, J.M., Garey, J.R., Shultz, J.W., 2004. Ecdysozoan phylogeny and Bayesian inference: first use of nearly complete 28S and 18S rRNA gene sequences to classify the arthropods and their kin. Mol. Phylogenet. Evol. 31, 178-191] that evaluated the phylogenetic relationships in Ecdysozoa (molting animals), especially arthropods. Here, the number of rRNA gene-sequences was effectively doubled for each major group of arthropods, and sequences from the phylum Kinorhyncha (mud dragons) were also included, bringing the number of ecdysozoan taxa to over 80. The methods emphasized maximum likelihood, Bayesian inference and statistical testing with parametric bootstrapping, but also included parsimony and minimum evolution. Prominent findings from our combined analysis of both genes are as follows. The fundamental subdivisions of Hexapoda (insects and relatives) are Insecta and Entognatha, with the latter consisting of collembolans (springtails) and a clade of proturans plus diplurans. Our rRNA-gene data provide the strongest evidence to date that the sister group of Hexapoda is Branchiopoda (fairy shrimps, tadpole shrimps, etc.), not Malacostraca. The large, Pancrustacea clade (hexapods within a paraphyletic Crustacea) divided into a few basic subclades: hexapods plus branchiopods; cirripedes (barnacles) plus malacostracans (lobsters, crabs, true shrimps, isopods, etc.); and the basally located clades of (a) ostracods (seed shrimps) and (b) branchiurans (fish lice) plus the bizarre pentastomids (tongue worms). These findings about Pancrustacea agree with a recent study by Regier, Shultz, and Kambic that used entirely different genes [Regier, J.C., Shultz, J.W., Kambic, R.E., 2005a. Pancrustacean phylogeny: hexapods are terrestrial crustaceans and maxillopods are not monophyletic. Proc. R. Soc. B 272, 395-401]. In Malacostraca, the stomatopod (mantis shrimp) was not at the base of the eumalacostracans, as is widely claimed, but grouped instead with an euphausiacean (krill). Within centipedes, Craterostigmus was the sister to all other pleurostigmophorans, contrary to the consensus view. Our trees also united myriapods (millipedes and centipedes) with chelicerates (horseshoe crabs, spiders, scorpions, and relatives) and united pycnogonids (sea spiders) with chelicerates, but with much less support than in the previous rRNA-gene study. Finally, kinorhynchs joined priapulans (penis worms) at the base of Ecdysozoa.

One hundred years of high-throughput Drosophila research

From the beginning, Drosophila was a high-throughput model organism. Unbiased and genome-wide efforts ranging from Morgan's search for spontaneous mutations and subsequent saturating loss-of-function and gain-of-function screens up to more recent techniques such as microarrays, proteomics and cellular assays have been and will continue to be the backbone of Drosophila research. Integrating these large datasets is one of the next challenges. However, once achieved, a plethora of information far exceeding the information content of the singular experiments will be revealed. Several high-throughput techniques and experimental strategies highlighting the unbiased and integrative nature of Drosophila research during the last century will be discussed.

Poriferan mtDNA and animal phylogeny based on mitochondrial gene arrangements

Phylogenetic relationships among the metazoan phyla are the subject of an ongoing controversy. Analysis of mitochondrial gene arrangements is a powerful tool to investigate these relationships; however, its previous application outside of individual animal phyla has been hampered by the lack of informative out-group data. To address this shortcoming, we determined complete mitochondrial DNA sequences for the demosponges Geodia neptuni and Tethya actinia, two representatives of the most basal animal phylum, the Porifera. With sponges as an outgroup, we investigated phylogenetic relationships of nine bilaterian phyla using both breakpoint analysis of global mitochondrial gene arrangements and maximum parsimony analysis of mitochondrial gene adjacencies. Our results provide strong support for a group that includes protostome (but not deuterostome) coelomate, pseudocoelomate, and acoelomate animals, thus clearly rejecting the Coelomata hypothesis. Two other groups of bilaterian animals, Lophotrochozoa and Ambulacraria, are also supported by our analyses. However, due to the remarkable stability of mitochondrial gene arrangements in Deuterostomia and the Ecdysozoa, conclusions on their evolutionary history cannot be drawn.

Conservation of novel Mahya genes shows the existence of neural functions common between Hymenoptera and Deuterostome

Honeybees have been shown to exhibit cognitive performances that were thought to be specific to some vertebrates. However, the molecular and cellular mechanisms of such cognitive abilities of the bees have not been understood. We have identified a novel gene, Mahya, expressed in the brain of the honeybee, Apis mellifera, and other Hymenoptera. Mahya orthologues are present in Deuterostomes but are absent or highly diverged in nematodes and, intriguingly, in two dipteran insects (fruit fly and mosquito) and Lepidoptera (silk moth). Mahya genes encode novel secretory proteins with a follistatin-like domain (Kazal-type serine/threonine protease inhibitor domain and EF-hand calcium-binding domain), two immunoglobulin domains, and a C-terminal novel domain. Honeybee Mahya is expressed in the mushroom bodies and antennal lobes of the brain. Zebra fish Mahya orthologues are expressed in the olfactory bulb, telencephalon, habenula, optic tectum, and cerebellum of the brain. Mouse Mahya orthologues are expressed in the olfactory bulb, hippocampus, and cerebellum of the brain. These results suggest that Mahya may be involved in learning and memory and in processing of sensory information in Hymenoptera and vertebrates. Furthermore, the limited existence of Mahya in the genomes of Hymenoptera and Deuterostomes supports the hypothesis that the genes typically represented by Mahya were lost or highly diverged during the evolution of the central nervous system of specific Bilaterian branches under the specific selection and subsequent adaptation associated with different ecologies and life histories.