The large srh family of chemoreceptor genes in Caenorhabditis nematodes reveals processes of genome evolution involving large duplications and deletions and intron gains and losses

G protein-coupled receptor-based thermosensation determines temperature acclimatization of Caenorhabditis elegans

Animals must sense and acclimatize to environmental temperatures for survival, yet their thermosensing mechanisms other than transient receptor potential (TRP) channels remain poorly understood. We identify a trimeric G protein-coupled receptor (GPCR), SRH-40, which confers thermosensitivity in sensory neurons regulating temperature acclimatization in Caenorhabditis elegans. Systematic knockdown of 1000 GPCRs by RNAi reveals GPCRs involved in temperature acclimatization, among which srh-40 is highly expressed in the ADL sensory neuron, a temperature-responsive chemosensory neuron, where TRP channels act as accessorial thermoreceptors. In vivo Ca2+ imaging demonstrates that an srh-40 mutation reduced the temperature sensitivity of ADL, resulting in supranormal temperature acclimatization. Ectopically expressing SRH-40 in a non-warmth-sensing gustatory neuron confers temperature responses. Moreover, temperature-dependent SRH-40 activation is reconstituted in Drosophila S2R+ cells. Overall, SRH-40 may be involved in thermosensory signaling underlying temperature acclimatization. We propose a dual thermosensing machinery through a GPCR and TRP channels in a single sensory neuron.

Dietary vitamin B12 regulates chemosensory receptor gene expression via the MEF2 transcription factor in Caenorhabditis elegans

Dynamic changes in chemoreceptor gene expression levels in sensory neurons are one strategy that an animal can use to modify their responses to dietary changes. However, the mechanisms underlying diet-dependent modulation of chemosensory gene expression are unclear. Here, we show that the expression of the srh-234 chemoreceptor gene localized in a single ADL sensory neuron type of Caenorhabditis elegans is downregulated when animals are fed a Comamonas aquatica bacterial diet, but not on an Escherichia coli diet. Remarkably, this diet-modulated effect on srh-234 expression is dependent on the micronutrient vitamin B12 endogenously produced by Comamonas aq. bacteria. Excess propionate and genetic perturbations in the canonical and shunt propionate breakdown pathways are able to override the repressive effects of vitamin B12 on srh-234 expression. The vitamin B12-mediated regulation of srh-234 expression levels in ADL requires the MEF-2 MADS domain transcription factor, providing a potential mechanism by which dietary vitamin B12 may transcriptionally tune individual chemoreceptor genes in a single sensory neuron type, which in turn may change animal responses to biologically relevant chemicals in their diet.

Cytogenomic analysis unveils mixed molecular evolution and recurrent chromosomal rearrangements shaping the multigene families on Schistocerca grasshopper genomes

Multigene families are essential components of eukaryotic genomes and play key roles either structurally and functionally. Their modes of evolution remain elusive even in the era of genomics, because multiple multigene family sequences coexist in genomes, particularly in large repetitive genomes. Here, we investigate how the multigene families 18S rDNA, U2 snDNA, and H3 histone evolved in 10 species of Schistocerca grasshoppers with very large and repeat-enriched genomes. Using sequenced genomes and fluorescence in situ hybridization mapping, we find substantial differences between species, including the number of chromosomal clusters, changes in sequence abundance and nucleotide composition, pseudogenization, and association with transposable elements (TEs). The intragenomic analysis of Schistocerca gregaria using long-read sequencing and genome assembly unveils conservation for H3 histone and recurrent pseudogenization for 18S rDNA and U2 snDNA, likely promoted by association with TEs and sequence truncation. Remarkably, TEs were frequently associated with truncated copies, were also among the most abundant in the genome, and revealed signatures of recent activity. Our findings suggest a combined effect of concerted and birth-and-death models driving the evolution of multigene families in Schistocerca over the last 8 million years, and the occurrence of intra- and interchromosomal rearrangements shaping their chromosomal distribution. Despite the conserved karyotype in Schistocerca, our analysis highlights the extensive reorganization of repetitive DNAs in Schistocerca, contributing to the advance of comparative genomics for this important grasshopper genus.

Chemosensory signal transduction in Caenorhabditis elegans

Chemosensory neurons translate perception of external chemical cues, including odorants, tastants, and pheromones, into information that drives attraction or avoidance motor programs. In the laboratory, robust behavioral assays, coupled with powerful genetic, molecular and optical tools, have made Caenorhabditis elegans an ideal experimental system in which to dissect the contributions of individual genes and neurons to ethologically relevant chemosensory behaviors. Here, we review current knowledge of the neurons, signal transduction molecules and regulatory mechanisms that underlie the response of C. elegans to chemicals, including pheromones. The majority of identified molecules and pathways share remarkable homology with sensory mechanisms in other organisms. With the development of new tools and technologies, we anticipate that continued study of chemosensory signal transduction and processing in C. elegans will yield additional new insights into the mechanisms by which this animal is able to detect and discriminate among thousands of chemical cues with a limited sensory neuron repertoire.

Large genetic diversity and strong positive selection in F-box and GPCR genes among the wild isolates of Caenorhabditis elegans

The F-box and chemosensory GPCR (csGPCR) gene families are greatly expanded in nematodes, including the model organism Caenorhabditis elegans, compared with insects and vertebrates. However, the intraspecific evolution of these two gene families in nematodes remain unexamined. In this study, we analyzed the genomic sequences of 330 recently sequenced wild isolates of C. elegans using a range of population genetics approaches. We found that F-box and csGPCR genes, especially the Srw family csGPCRs, showed much more diversity than other gene families. Population structure analysis and phylogenetic analysis divided the wild strains into eight non-Hawaiian and three Hawaiian subpopulations. Some Hawaiian strains appeared to be more ancestral than all other strains. F-box and csGPCR genes maintained a great amount of the ancestral variants in the Hawaiian subpopulation and their divergence among the non-Hawaiian subpopulations contributed significantly to population structure. F-box genes are mostly located at the chromosomal arms and high recombination rate correlates with their large polymorphism. Moreover, using both neutrality tests and extended haplotype homozygosity analysis, we identified signatures of strong positive selection in the F-box and csGPCR genes among the wild isolates, especially in the non-Hawaiian population. Accumulation of high-frequency-derived alleles in these genes was found in non-Hawaiian population, leading to divergence from the ancestral genotype. In summary, we found that F-box and csGPCR genes harbor a large pool of natural variants, which may be subjected to positive selection. These variants are mostly mapped to the substrate-recognition domains of F-box proteins and the extracellular and intracellular regions of csGPCRs, possibly resulting in advantages during adaptation by affecting protein degradation and the sensing of environmental cues, respectively.

Chemogenomic approach to identifying nematode chemoreceptor drug targets in the entomopathogenic nematode Heterorhabditis bacteriophora

Parasitic nematodes constitute one of the major threats to human health, causing diseases of major socioeconomic importance worldwide. Recent estimates indicate that more than 1 billion people are infected with parasitic nematodes around the world. Current measures to combat parasitic nematode infections include anthelmintic drugs. However, heavy exposure to anthelmintics has selected populations of livestock parasitic nematodes that are no longer susceptible to the drugs, rendering several anthelmintics useless for parasitic nematode control in many areas of the world. The rapidity with which anthelmintic resistance developed in response to these drugs suggests that increasing the selective pressure on human parasitic nematodes will also rapidly generate resistant worm populations. Therefore, development of new anthelmintics is of major importance before resistance becomes widespread in human parasitic nematode populations. G-Protein Coupled Receptors (GPCRs) represent an important target for many pharmacological interventions due to their ubiquitous expression in various cell types. GPCRs contribute to numerous physiological processes, and their ligand binding sites located on cell surfaces make them accessible targets and attractive substrates in terms of druggability. In fact, ∼35 % of Food and Drug Administration (FDA) and European Medicines Agency (EMA) approved drugs target GPCRs and their associated proteins, with over 300 additional drugs targeting GPCRs at the clinical trial stage. Nematode Chemosensory GPCRs (NemChRs) are unique to nematodes, and therefore represent ideal substrates for target-based drug discovery. Here we set out to identify NemChRs that are transcriptionally active inside the host, and to use these NemChRs in a reverse pharmacological screen to impede parasitic development. Our data identified several NemChRs, and we focused on one that was expressed in neuronal cells and exhibited the highest fold change in transcription after host activation. Next, we performed homology modelling and molecular dynamics simulations of this NemChR in order to conduct a virtual screening campaign to identify candidate drug targets which were ranked and selected for experimental testing in bioassays. Taken together, our results identify and characterize a candidate NemChR drug target, and provide a chemogenomic pipeline for identifying nematicide substrates.

Comprehensive analyses of ZFP gene family and characterization of expression profiles during plant hormone response in cotton

BACKGROUND

Zinc finger proteins (ZFPs) containing only a single zinc finger domain play important roles in the regulation of plant growth and development, as well as in biotic and abiotic stress responses. To date, the evolutionary history and functions of the ZFP gene family have not been identified in cotton.

RESULTS

In this paper, we identified 29 ZFP genes in Gossypium hirsutum. This gene family was divided into seven subfamilies, 22 of which were distributed over 17 chromosomes. Bioinformatic analysis revealed that 20 GhZFP genes originated from whole genome duplications and two originated from dispersed duplication events, indicating that whole genome duplication is the main force in the expansion of the GhZFP gene family. Most GhZFP8 subfamily genes, except for GhZFP8-3, were highly expressed during fiber cell growth, and were induced by brassinosteroids in vitro. Furthermore, we found that a large number of GhZFP genes contained gibberellic acid responsive elements, auxin responsive elements, and E-box elements in their promoter regions. Exogenous application of these hormones significantly stimulated the expression of these genes.

CONCLUSIONS

Our findings reveal that GhZFP8 genes are involved in cotton fiber development and widely induced by auxin, gibberellin and BR, which provides a foundation for the identification of more downstream genes with potential roles in phytohormone stimuli, and a basis for breeding better cotton varieties in the future.

Genome-Wide Identification and Analysis of Class III Peroxidases in Allotetraploid Cotton (Gossypium hirsutum L.) and their Responses to PK Deficiency

Class III peroxidases (PODs), commonly known as secretable class III plant peroxidases, are plant-specific enzymes that play critical roles in not only plant growth and development but also the responses to biotic and abiotic stress. In this study, we identified 198 nonredundant POD genes, designated GhPODs, with 180 PODs being predicted to secrete into apoplast. These POD genes were divided into 10 sub-groups based on their phylogenetic relationships. We performed systematic bioinformatic analysis of the POD genes, including analysis of gene structures, phylogenetic relationships, and gene expression profiles. The GhPODs are unevenly distributed on both upland cotton sub-genome A and D chromosomes. Additionally, these genes have undergone 15 segmental and 12 tandem duplication events, indicating that both segmental and tandem duplication contributed to the expansion of the POD gene family in upland cotton. Ka/Ks analysis suggested that most duplicated GhPODs experienced negative selection, with limited functional divergence during the duplication events. High-throughput RNA-seq data indicated that most highly expressed genes might play significant roles in root, stem, leaf, and fiber development. Under K or P deficiency conditions, PODs showed different expression patterns in cotton root and leaf. This study provides useful information for further functional analysis of the POD gene family in upland cotton.

Complex Evolutionary Dynamics of Massively Expanded Chemosensory Receptor Families in an Extreme Generalist Chelicerate Herbivore

While mechanisms to detoxify plant produced, anti-herbivore compounds have been associated with plant host use by herbivores, less is known about the role of chemosensory perception in their life histories. This is especially true for generalists, including chelicerate herbivores that evolved herbivory independently from the more studied insect lineages. To shed light on chemosensory perception in a generalist herbivore, we characterized the chemosensory receptors (CRs) of the chelicerate two-spotted spider mite, Tetranychus urticae, an extreme generalist. Strikingly, T. urticae has more CRs than reported in any other arthropod to date. Including pseudogenes, 689 gustatory receptors were identified, as were 136 degenerin/Epithelial Na+ Channels (ENaCs) that have also been implicated as CRs in insects. The genomic distribution of T. urticae gustatory receptors indicates recurring bursts of lineage-specific proliferations, with the extent of receptor clusters reminiscent of those observed in the CR-rich genomes of vertebrates or C. elegans Although pseudogenization of many gustatory receptors within clusters suggests relaxed selection, a subset of receptors is expressed. Consistent with functions as CRs, the genomic distribution and expression of ENaCs in lineage-specific T. urticae expansions mirrors that observed for gustatory receptors. The expansion of ENaCs in T. urticae to > 3-fold that reported in other animals was unexpected, raising the possibility that ENaCs in T. urticae have been co-opted to fulfill a major role performed by unrelated CRs in other animals. More broadly, our findings suggest an elaborate role for chemosensory perception in generalist herbivores that are of key ecological and agricultural importance.

Characterization of the arginine kinase isoforms in Caenorhabditis elegans

Phosphagen kinases (PKs) are well-studied enzymes involved in energy homeostasis in a wide range of animal, protozoan, and even some bacterial species. Recent genome efforts have allowed comparative work on the PKs to extend beyond the biochemistry of individual proteins to the comparative cellular physiology and examining of the role of all PK family members in an organism. The sequencing of the Caenorhabditis elegans genome and availability of sophisticated genetic tools within that system affords the opportunity to conduct a detailed physiological analysis of the PKs from a well known invertebrate for comparison with the extensive work conducted on vertebrate systems. As a first step in this effort we have carried out a detailed molecular genetic and biochemical characterization of the PKs in C. elegans. Our results reveal that C. elegans has five PK genes encoding arginine kinases that range in catalytic efficiency (kcat/KM(Arg)) from (3.1±0.6)×10(4) to (9±4)×10(5) M(-1) s(-1). This range is generally within the range seen for arginine kinases from a variety of species. Our molecular genetic and phylogenetic analysis reveals that the gene family has undergone extensive intron loss and gain within the suborder Rhabditina. In addition, within C. elegans we find evidence of gene duplication and loss. The analysis described here for the C. elegans AKs represents one of the most complete biochemical and molecular genetic analysis of a PK family within a genetically tractable invertebrate system and opens up the possibility of conducting detailed physiological comparisons with vertebrate systems using the sophisticated tools available with this model invertebrate system.

Inventory of the cichlid olfactory receptor gene repertoires: identification of olfactory genes with more than one coding exon

BACKGROUND

To help understand the molecular mechanisms underlying the remarkable phenotypic diversity displayed by cichlids, the genome sequences of O. niloticus, P. nyererei, H. burtoni, N. brichardi and M. zebra were recently determined. Here, we present the contents of the olfactory receptor (OR) repertoires in the genomes of these five fishes.

RESULTS

We performed an exhaustive TBLASTN search of the five cichlid genomes to identify their OR repertoires as completely as possible. We used as bait a set of ORs described in the literature. The cichlid repertoires thereby extracted contained large numbers of complete genes (O. niloticus 158; H. burtoni 90; M. zebra 102; N. brichardi 69; P. nyererei 88), a small numbers of pseudogenes and many "edge genes" corresponding to incomplete genes located at the ends of contigs. A phylogenetic tree was constructed and showed these repertoires include a large number of families and subfamilies. It also allowed the identification of a large number of OR analogues between cichlids with very high amino-acid identity (≥ 99%). Nearly 9% of the full-length cichlid OR genes are composed of several coding exons. This is very unusual for vertebrate OR genes. Nevertheless, the evidence is strong, and includes the donor and acceptor splice junction sequences; also, the positions of these genes in the phylogenetic tree indicate that they constitute subfamilies well apart from non-OR G protein-coupled receptor families.

CONCLUSIONS

Cichlid OR repertoires are made up of a larger number of genes and fewer pseudogenes than those in other teleosts except zebrafish. These ORs share all identified properties common to all fish ORs; however, the large number of families and subfamilies, each containing few ORs implies that they have evolved more rapidly. This high level of OR diversity is consistent with the substantial phenotypic diversity that characterizes cichlids.

Multigenic natural variation underlies Caenorhabditis elegans olfactory preference for the bacterial pathogen Serratia marcescens

The nematode Caenorhabditis elegans can use olfaction to discriminate among different kinds of bacteria, its major food source. We asked how natural genetic variation contributes to choice behavior, focusing on differences in olfactory preference behavior between two wild-type C. elegans strains. The laboratory strain N2 strongly prefers the odor of Serratia marcescens, a soil bacterium that is pathogenic to C. elegans, to the odor of Escherichia coli, a commonly used laboratory food source. The divergent Hawaiian strain CB4856 has a weaker attraction to Serratia than the N2 strain, and this behavioral difference has a complex genetic basis. At least three quantitative trait loci (QTLs) from the CB4856 Hawaii strain (HW) with large effect sizes lead to reduced Serratia preference when introgressed into an N2 genetic background. These loci interact and have epistatic interactions with at least two antagonistic QTLs from HW that increase Serratia preference. The complex genetic architecture of this C. elegans trait is reminiscent of the architecture of mammalian metabolic and behavioral traits.

Neuropeptide gene families in Caenorhabditis elegans

Neuropeptides are short sequences ofamino acids that function in all multicellular organisms to communicate information between cells. The first sequence ofa neuropeptide was reported in 1970' and the number of identified neuropeptides remained relatively small until the 1990s when the DNA sequence of multiple genomes revealed treasure troves ofinformation. Byblasting away at the genome, gene families, the sizes ofwhich were previously unknown, could now be determined. This information has led to an exponential increase in the number of putative neuropeptides and their respective gene families. The molecular biology age greatly benefited the neuropeptide field in the nematode Caenorhabditis elegans. Its genome was among the first to be sequenced and this allowed us the opportunity to screen the genome for neuropeptide genes. Initially, the screeningwas slow, as the Genefinder and BLAST programs had difficulty identifying small genes and peptides. However, as the bioinformatics programs improved, the extent of the neuropeptide gene families in C. elegans gradually emerged.

RNAi screen of DAF-16/FOXO target genes in C. elegans links pathogenesis and dauer formation

The DAF-16/FOXO transcription factor is the major downstream output of the insulin/IGF1R signaling pathway controlling C. elegans dauer larva development and aging. To identify novel downstream genes affecting dauer formation, we used RNAi to screen candidate genes previously identified to be regulated by DAF-16. We used a sensitized genetic background [eri-1(mg366); sdf-9(m708)], which enhances both RNAi efficiency and constitutive dauer formation (Daf-c). Among 513 RNAi clones screened, 21 displayed a synthetic Daf-c (SynDaf) phenotype with sdf-9. One of these genes, srh-100, was previously identified to be SynDaf, but twenty have not previously been associated with dauer formation. Two of the latter genes, lys-1 and cpr-1, are known to participate in innate immunity and six more are predicted to do so, suggesting that the immune response may contribute to the dauer decision. Indeed, we show that two of these genes, lys-1 and clc-1, are required for normal resistance to Staphylococcus aureus. clc-1 is predicted to function in epithelial cohesion. Dauer formation exhibited by daf-8(m85), sdf-9(m708), and the wild-type N2 (at 27°C) were all enhanced by exposure to pathogenic bacteria, while not enhanced in a daf-22(m130) background. We conclude that knockdown of the genes required for proper pathogen resistance increases pathogenic infection, leading to increased dauer formation in our screen. We propose that dauer larva formation is a behavioral response to pathogens mediated by increased dauer pheromone production.

Genome-wide comparative gene family classification

Correct classification of genes into gene families is important for understanding gene function and evolution. Although gene families of many species have been resolved both computationally and experimentally with high accuracy, gene family classification in most newly sequenced genomes has not been done with the same high standard. This project has been designed to develop a strategy to effectively and accurately classify gene families across genomes. We first examine and compare the performance of computer programs developed for automated gene family classification. We demonstrate that some programs, including the hierarchical average-linkage clustering algorithm MC-UPGMA and the popular Markov clustering algorithm TRIBE-MCL, can reconstruct manual curation of gene families accurately. However, their performance is highly sensitive to parameter setting, i.e. different gene families require different program parameters for correct resolution. To circumvent the problem of parameterization, we have developed a comparative strategy for gene family classification. This strategy takes advantage of existing curated gene families of reference species to find suitable parameters for classifying genes in related genomes. To demonstrate the effectiveness of this novel strategy, we use TRIBE-MCL to classify chemosensory and ABC transporter gene families in C. elegans and its four sister species. We conclude that fully automated programs can establish biologically accurate gene families if parameterized accordingly. Comparative gene family classification finds optimal parameters automatically, thus allowing rapid insights into gene families of newly sequenced species.

Orthology confers intron position conservation

Background

With the wealth of genomic data available it has become increasingly important to assign putative protein function through functional transfer between orthologs. Therefore, correct elucidation of the evolutionary relationships among genes is a critical task, and attempts should be made to further improve the phylogenetic inference by adding relevant discriminating features. It has been shown that introns can maintain their position over long evolutionary timescales. For this reason, it could be possible to use conservation of intron positions as a discriminating factor when assigning orthology. Therefore, we wanted to investigate whether orthologs have a higher degree of intron position conservation (IPC) compared to non-orthologous sequences that are equally similar in sequence.

Results

To this end, we developed a new score for IPC and applied it to ortholog groups between human and six other species. For comparison, we also gathered the closest non-orthologs, meaning sequences close in sequence space, yet falling just outside the ortholog cluster. We found that ortholog-ortholog gene pairs on average have a significantly higher degree of IPC compared to ortholog-closest non-ortholog pairs. Also pairs of inparalogs were found to have a higher IPC score than inparalog-closest non-inparalog pairs. We verified that these differences can not simply be attributed to the generally higher sequence identity of the ortholog-ortholog and the inparalog-inparalog pairs.

Furthermore, we analyzed the agreement between IPC score and the ortholog score assigned by the InParanoid algorithm, and found that it was consistently high for all species comparisons. In a minority of cases, the IPC and InParanoid score ranked inparalogs differently. These represent cases where sequence and intron position divergence are discordant. We further analyzed the discordant clusters to identify any possible preference for protein functions by looking for enriched GO terms and Pfam protein domains. They were enriched for functions important for multicellularity, which implies a connection between shifts in intronic structure and the origin of multicellularity.

Conclusions

We conclude that orthologous genes tend to have more conserved intron positions compared to non-orthologous genes. As a consequence, our IPC score is useful as an additional discriminating factor when assigning orthology.

Copy number variation in the genomes of twelve natural isolates of Caenorhabditis elegans

Background

Copy number variation is an important component of genetic variation in higher eukaryotes. The extent of natural copy number variation in C. elegans is unknown outside of 2 highly divergent wild isolates and the canonical N2 Bristol strain.

Results

We have used array comparative genomic hybridization (aCGH) to detect copy number variation in the genomes of 12 natural isolates of Caenorhabditis elegans. Deletions relative to the canonical N2 strain are more common in these isolates than duplications, and indels are enriched in multigene families on the autosome arms. Among the strains in our study, the Hawaiian and Madeiran strains (CB4856 and JU258) carry the largest number of deletions, followed by the Vancouver strain (KR314). Overall we detected 510 different deletions affecting 1136 genes, or over 5% of the genes in the canonical N2 genome. The indels we identified had a median length of 2.7 kb. Since many deletions are found in multiple isolates, deletion loci were used as markers to derive an unrooted tree to estimate genetic relatedness among the strains.

Conclusion

Copy number variation is extensive in C. elegans, affecting over 5% of the genes in the genome. The deletions we have detected in natural isolates of C. elegans contribute significantly to the number of deletion alleles available to researchers. The relationships between strains are complex and different regions of the genome possess different genealogies due to recombination throughout the natural history of the species, which may not be apparent in studies utilizing smaller numbers of genetic markers.

Evolutionary mirages: selection on binding site composition creates the illusion of conserved grammars in Drosophila enhancers

The clustering of transcription factor binding sites in developmental enhancers and the apparent preferential conservation of clustered sites have been widely interpreted as proof that spatially constrained physical interactions between transcription factors are required for regulatory function. However, we show here that selection on the composition of enhancers alone, and not their internal structure, leads to the accumulation of clustered sites with evolutionary dynamics that suggest they are preferentially conserved. We simulated the evolution of idealized enhancers from Drosophila melanogaster constrained to contain only a minimum number of binding sites for one or more factors. Under this constraint, mutations that destroy an existing binding site are tolerated only if a compensating site has emerged elsewhere in the enhancer. Overlapping sites, such as those frequently observed for the activator Bicoid and repressor Krüppel, had significantly longer evolutionary half-lives than isolated sites for the same factors. This leads to a substantially higher density of overlapping sites than expected by chance and the appearance that such sites are preferentially conserved. Because D. melanogaster (like many other species) has a bias for deletions over insertions, sites tended to become closer together over time, leading to an overall clustering of sites in the absence of any selection for clustered sites. Since this effect is strongest for the oldest sites, clustered sites also incorrectly appear to be preferentially conserved. Following speciation, sites tend to be closer together in all descendent species than in their common ancestors, violating the common assumption that shared features of species' genomes reflect their ancestral state. Finally, we show that selection on binding site composition alone recapitulates the observed number of overlapping and closely neighboring sites in real D. melanogaster enhancers. Thus, this study calls into question the common practice of inferring "cis-regulatory grammars" from the organization and evolutionary dynamics of developmental enhancers.

Identification of novel aspartic proteases from Strongyloides ratti and characterisation of their evolutionary relationships, stage-specific expression and molecular structure

BACKGROUND

Aspartic proteases are known to play an important role in the biology of nematode parasitism. This role is best characterised in blood-feeding nematodes, where they digest haemoglobin, but they are also likely to play important roles in the biology of nematode parasites that do not feed on blood. In the present work, we investigate the evolution and expression of aspartic proteases in Strongyloides ratti, which permits a unique comparison between parasitic and free-living adult forms within its life-cycle.

RESULTS

We identified eight transcribed aspartic protease sequences and a further two genomic sequences and compared these to homologues in Caenorhabditis elegans and other nematode species. Phylogenetic analysis demonstrated a complex pattern of gene evolution, such that some S. ratti sequences had a one-to-one correspondence with orthologues of C. elegans but that lineage-specific expansions have occurred for other aspartic proteases in these two nematodes. These gene duplication events may have contributed to the adaptation of the two species to their different lifestyles. Among the set of S. ratti aspartic proteases were two closely-related isoforms that showed differential expression during different life stages: ASP-2A is highly expressed in parasitic females while ASP-2B is predominantly found in free-living adults. Molecular modelling of the ASP-2 isoforms reveals that their substrate specificities are likely to be very similar, but that ASP-2B is more electrostatically negative over its entire molecular surface than ASP-2A. This characteristic may be related to different pH values of the environments in which these two isoforms operate.

CONCLUSIONS

We have demonstrated that S. ratti provides a powerful model to explore the genetic adaptations associated with parasitic versus free-living life-styles. We have discovered gene duplication of aspartic protease genes in Strongyloides and identified a pair of paralogues differentially expressed in either the parasitic or the free-living phase of the nematode life-cycle, consistent with an adaptive role for aspartic proteases in the evolution of nematode parasitism.

Polymorphic segmental duplication in the nematode Caenorhabditis elegans

BACKGROUND

The nematode Caenorhabditis elegans was the first multicellular organism to have its genome fully sequenced. Over the last 10 years since the original publication in 1998, the C. elegans genome has been scrutinized and the last gaps were filled in November 2002, which present a unique opportunity for examining genome-wide segmental duplications.

RESULTS

Here, we performed analysis of the C. elegans genome in search for segmental duplications using a new tool -- OrthoCluster -- we have recently developed. We detected 3,484 duplicated segments -- duplicons -- ranging in size from 234 bp to 108 Kb. The largest pair of duplicons, 108 kb in length located on the left arm of Chromosome V, was further characterized. They are nearly identical at the DNA level (99.7% identity) and each duplicon contains 26 putative protein coding genes. Genotyping of 76 wild-type strains obtained from different labs in the C. elegans community revealed that not all strains contain this duplication. In fact, only 29 strains carry this large segmental duplication, suggesting a very recent duplication event in the C. elegans genome.

CONCLUSION

This report represents the first demonstration that the C. elegans laboratory wild-type N2 strains has acquired large-scale differences.

Power-laws in the genomic distribution of coding segments in several organisms: an evolutionary trace of segmental duplications, possible paleopolyploidy and gene loss

Large-scale features of the spatial arrangement of protein-coding segments (PCS) are investigated by means of the inter-PCS spacers' size distributions, which have been found to follow power-laws. Linearity in double-logarithmic scale extends to several orders of magnitude in the genomes of organisms as disparate as mammals, insects and plants. This feature is also present in the most compact eukaryotic genomes and in half of the examined bacteria, despite their very limited non-coding space. We have tried to determine the sequence of events in the course of genomes' evolution which may account for the formation of the observed size distributions. The proposed mechanism essentially includes two types of events: (i) segmental duplications (and possibly paleopolyploidy), and (ii) the subsequent loss of most of the duplicated genes. It is shown by computer simulations that the formulated scenario generates power-law-like inter-PCS spacers' size distributions, which remain robust for a variety of parameter choices, even if insertion of external sequences, such as viruses or proliferating retroelements is included. Moreover, power-laws are preserved after most of the non-coding DNA has been removed, thus explaining the finding of this pattern in genomes as compact as that of Takifugu rubripes.

Deletional bias across the three domains of life

Elevated levels of genetic drift are hypothesized to be a dominant factor that influences genome size evolution across all life-forms. However, increased levels of drift appear to be correlated with genome expansion in eukaryotes but with genome contraction in bacteria, suggesting that these two groups of organisms experience vastly different mutational inputs and selective constraints. To determine the contribution of small insertion and deletion events to the differences in genome organization between eukaryotes and prokaryotes, we systematically surveyed 17 taxonomic groups across the three domains of life. Based on over 5,000 indel events in noncoding regions, we found that deletional events outnumbered insertions in all groups examined. The extent of deletional bias, when measured by the total length of insertions to deletions, revealed a marked disparity between eukaryotes and prokaryotes, whereas the ratio was close to one in the three eukaryotic groups examined, deletions outweighed insertions by at least a factor of 10 in most prokaryotes. Moreover, the strength of deletional bias is associated with the proportion of coding regions in prokaryotic genomes. Considering that genetic drift is a stochastic process and does not discriminate the exact nature of mutations, the degree of bias toward deletions provides an explanation to the differential responses of eukaryotes and prokaryotes to elevated levels of drift. Furthermore, deletional bias, rather than natural selection, is the primary mechanism by which the compact gene packing within most prokaryotic genomes is maintained.

Evolution of the Caenorhabditis elegans genome

A fundamental problem in genome biology is to elucidate the evolutionary forces responsible for generating nonrandom patterns of genome organization. As the first metazoan to benefit from full-genome sequencing, Caenorhabditis elegans has been at the forefront of research in this area. Studies of genomic patterns, and their evolutionary underpinnings, continue to be augmented by the recent push to obtain additional full-genome sequences of related Caenorhabditis taxa. In the near future, we expect to see major advances with the onset of whole-genome resequencing of multiple wild individuals of the same species. In this review, we synthesize many of the important insights to date in our understanding of genome organization and function that derive from the evolutionary principles made explicit by theoretical population genetics and molecular evolution and highlight fertile areas for future research on unanswered questions in C. elegans genome evolution. We call attention to the need for C. elegans researchers to generate and critically assess nonadaptive hypotheses for genomic and developmental patterns, in addition to adaptive scenarios. We also emphasize the potential importance of evolution in the gonochoristic (female and male) ancestors of the androdioecious (hermaphrodite and male) C. elegans as the source for many of its genomic and developmental patterns.

Evolution of the gene lineage encoding the carbon dioxide receptor in insects

A heterodimer of the insect chemoreceptors Gr21a and Gr63a has been shown to be the carbon dioxide receptor in Drosophila melanogaster (Meigen) (Diptera: Drosophilidae). Comparison of the genes encoding these two proteins across the 12 available drosophilid fly genomes allows refined definition of their N-termini. These genes are highly conserved, along with a paralog of Gr21a, in the Anopheles gambiae, Aedes aegypti, and Culex pipiens mosquitoes, as well as in the silk moth Bombyx mori and the red flour beetle Tribolium castaneum. In the latter four species we name these three proteins Gr1, Gr2, and Gr3. Intron evolution within this distinctive three gene lineage is considerable, with at least 13 inferred gains and 39 losses. Surprisingly, this entire ancient gene lineage is absent from all other available more basal insect and related arthropod genomes, specifically the honey bee, parasitoid wasp, human louse, pea aphid, waterflea, and blacklegged tick genomes. At least two of these species can detect carbon dioxide, suggesting that they evolved other means to do so.

High nucleotide divergence in developmental regulatory genes contrasts with the structural elements of olfactory pathways in caenorhabditis

Almost all organismal function is controlled by pathways composed of interacting genetic components. The relationship between pathway structure and the evolution of individual pathway components is not completely understood. For the nematode Caenorhabditis elegans, chemosensory pathways regulate critical aspects of an individual's life history and development. To help understand how olfaction evolves in Caenorhabditis and to examine patterns of gene evolution within transduction pathways in general, we analyzed nucleotide variation within and between species across two well-characterized olfactory pathways, including regulatory genes controlling the fate of the cells in which the pathways are expressed. In agreement with previous studies, we found much higher levels of polymorphism within C. remanei than within the related species C. elegans and C. briggsae. There are significant differences in the rates of nucleotide evolution for genes across the two pathways but no particular association between evolutionary rate and gene position, suggesting that the evolution of functional pathways must be considered within the context of broader gene network structure. However, developmental regulatory genes show both higher levels of divergence and polymorphism than the structural genes of the pathway. These results show that, contrary to the emerging paradigm in the evolution of development, important structural changes can accumulate in transcription factors.

The Caenorhabditis globin gene family reveals extensive nematode-specific radiation and diversification

Background

Globin isoforms with variant properties and functions have been found in the pseudocoel, body wall and cuticle of various nematode species and even in the eyespots of the insect-parasite Mermis nigrescens. In fact, much higher levels of complexity exist, as shown by recent whole genome analysis studies. In silico analysis of the genome of Caenorhabditis elegans revealed an unexpectedly high number of globin genes featuring a remarkable diversity in gene structure, amino acid sequence and expression profiles.

Results

In the present study we have analyzed whole genomic data from C. briggsae, C. remanei, Pristionchus pacificus and Brugia malayi and EST data from several other nematode species to study the evolutionary history of the nematode globin gene family. We find a high level of conservation of the C. elegans globin complement, with even distantly related nematodes harboring orthologs to many Caenorhabditis globins. Bayesian phylogenetic analysis resolves all nematode globins into two distinct globin classes. Analysis of the globin intron-exon structures suggests extensive loss of ancestral introns and gain of new positions in deep nematode ancestors, and mainly loss in the Caenorhabditis lineage. We also show that the Caenorhabditis globin genes are expressed in distinct, mostly non-overlapping, sets of cells and that they are all under strong purifying selection.

Conclusion

Our results enable reconstruction of the evolutionary history of the globin gene family in the nematode phylum. A duplication of an ancestral globin gene occurred before the divergence of the Platyhelminthes and the Nematoda and one of the duplicated genes radiated further in the nematode phylum before the split of the Spirurina and Rhabditina and was followed by further radiation in the lineage leading to Caenorhabditis. The resulting globin genes were subject to processes of subfunctionalization and diversification leading to cell-specific expression patterns. Strong purifying selection subsequently dampened further evolution and facilitated fixation of the duplicated genes in the genome.

The Caenorhabditis chemoreceptor gene families

BACKGROUND

Chemoreceptor proteins mediate the first step in the transduction of environmental chemical stimuli, defining the breadth of detection and conferring stimulus specificity. Animal genomes contain families of genes encoding chemoreceptors that mediate taste, olfaction, and pheromone responses. The size and diversity of these families reflect the biology of chemoperception in specific species.

RESULTS

Based on manual curation and sequence comparisons among putative G-protein-coupled chemoreceptor genes in the nematode Caenorhabditis elegans, we identified approximately 1300 genes and 400 pseudogenes in the 19 largest gene families, most of which fall into larger superfamilies. In the related species C. briggsae and C. remanei, we identified most or all genes in each of the 19 families. For most families, C. elegans has the largest number of genes and C. briggsae the smallest number, suggesting changes in the importance of chemoperception among the species. Protein trees reveal family-specific and species-specific patterns of gene duplication and gene loss. The frequency of strict orthologs varies among the families, from just over 50% in two families to less than 5% in three families. Several families include large species-specific expansions, mostly in C. elegans and C. remanei.

CONCLUSION

Chemoreceptor gene families in Caenorhabditis species are large and evolutionarily dynamic as a result of gene duplication and gene loss. These dynamics shape the chemoreceptor gene complements in Caenorhabditis species and define the receptor space available for chemosensory responses. To explain these patterns, we propose the gray pawn hypothesis: individual genes are of little significance, but the aggregate of a large number of diverse genes is required to cover a large phenotype space.

Use of shotgun proteomics for the identification, confirmation, and correction of C. elegans gene annotations

We describe a general mass spectrometry-based approach for gene annotation of any organism and demonstrate its effectiveness using the nematode Caenorhabditis elegans. We detected 6779 C. elegans proteins (67,047 peptides), including 384 that, although annotated in WormBase WS150, lacked cDNA or other prior experimental support. We also identified 429 new coding sequences that were unannotated in WS150. Nearly half (192/429) of the new coding sequences were confirmed with RT-PCR data. Thirty-three (approximately 8%) of the new coding sequences had been predicted to be pseudogenes, 151 (approximately 35%) reveal apparent errors in gene models, and 245 (57%) appear to be novel genes. In addition, we verified 6010 exon-exon splice junctions within existing WormBase gene models. Our work confirms that mass spectrometry is a powerful experimental tool for annotating sequenced genomes. In addition, the collection of identified peptides should facilitate future proteomics experiments targeted at specific proteins of interest.

Genome evolution in Caenorhabditis

Since the completion of the Caenorhabditis elegans genome sequence 10 years ago, efforts of the large community of C. elegans geneticists have resulted in a high-quality annotation of the structures and sequence relatedness of nearly all the protein encoding and RNA genes. Based on increasingly accurate gene counts in other species, it now appears that C. elegans has more functional genes than most insects and approximately the same number as most mammals. In the last few years, draft genome sequences for several other nematodes have been published (C. briggsae and Brugia malayi) or publicly released (C. remanei, C. brenneri, C. japonica, Pristionchus pacificus, Trichinella spiralis and Haemonchus contortus). Comparisons of gene content within the phylum and to other phyla reveal complex patterns of genome evolution. These patterns include substantial numbers of genes conserved across all the major metazoan phyla (core metazoan genes) and many nematode-specific genes and gene families. Nematode-specific genes are located predominantly on autosomal arms, which also have higher recombination rates. It appears that evolutionary innovations occur mostly in these regions, probably facilitated by higher recombination. Few of these genes have gross phenotypes when knocked down by RNAi, suggesting that many of them function in specific aspects of nematode biology that were not tested, including chemosensation, pathogen response and xenobiotic detoxification.

Searching for the ligands of odorant receptors

Through the sense of smell mammals can detect and discriminate between a large variety of odorants present in the surrounding environment. Odorants bind to a large repertoire of odorant receptors located in the cilia of olfactory sensory neurons of the nose. Each olfactory neuron expresses one single type of odorant receptor, and neurons expressing the same type of receptor project their axons to one or a few glomeruli in the olfactory bulb, creating a map of odorant receptor inputs. The information is then passed on to other regions of the brain, leading to odorant perception. To understand how the olfactory system discriminates between odorants, it is necessary to determine the odorant specificities of individual odorant receptors. These studies are complicated by the extremely large size of the odorant receptor family and by the poor functional expression of these receptors in heterologous cells. This article provides an overview of the methods that are currently being used to investigate odorant receptor-ligand interactions.

Generation and modulation of chemosensory behaviors in C. elegans

C. elegans recognizes and discriminates among hundreds of chemical cues using a relatively compact chemosensory nervous system. Chemosensory behaviors are also modulated by prior experience and contextual cues. Because of the facile genetics and genomics possible in this organism, C. elegans provides an excellent system in which to explore the generation of chemosensory behaviors from the level of a single gene to the motor output. This review summarizes the current knowledge on the molecular and neuronal substrates of chemosensory behaviors and chemosensory behavioral plasticity in C. elegans.

Mechanisms and evolution of environmental responses in Caenorhabditis elegans

We review mechanistic and evolutionary aspects of interactions between the model organism Caenorhabditis elegans and its environment. In particular, we focus on environmental effects affecting developmental mechanisms. We describe natural and laboratory environments of C. elegans and provide an overview of the different environmental responses of this organism. We then show how two developmental processes respond to changes in the environment. First, we discuss the development of alternative juvenile stages, the dauer and non-dauer larva. This example illustrates how development responds to variation in the environment to generate complex phenotypic variation. Second, we discuss the development of the C. elegans vulva. This example illustrates how development responds to variation in the environment while generating an invariant final phenotype.

Origins and Evolution of Spliceosomal Introns

Research into the origins of introns is at a critical juncture in the resolution of theories on the evolution of early life (which came first, RNA or DNA?), the identity of LUCA (the last universal common ancestor, was it prokaryotic- or eukaryotic-like?), and the significance of noncoding nucleotide variation. One early notion was that introns would have evolved as a component of an efficient mechanism for the origin of genes. But alternative theories emerged as well. From the debate between the "introns-early" and "introns-late" theories came the proposal that introns arose before the origin of genetically encoded proteins and DNA, and the more recent "introns-first" theory, which postulates the presence of introns at that early evolutionary stage from a reconstruction of the "RNA world." Here we review seminal and recent ideas about intron origins. Recent discoveries about the patterns and causes of intron evolution make this one of the most hotly debated and exciting topics in molecular evolutionary biology today.

Noncell- and cell-autonomous G-protein-signaling converges with Ca2+/mitogen-activated protein kinase signaling to regulate str-2 receptor gene expression in Caenorhabditis elegans

In the sensory system of C. elegans, the candidate odorant receptor gene str-2 is strongly expressed in one of the two AWC neurons and weakly in both ASI neurons. Asymmetric AWC expression results from suppression of str-2 expression by a Ca2+/MAPK signaling pathway in one of the AWC neurons early in development. Here we show that the same Ca2+/MAPK pathway promotes str-2 expression in the AWC and ASI neurons together with multiple cell-autonomous and noncell-autonomous G-protein-signaling pathways. In first-stage larvae and adult animals, signals mediated by the Galpha subunits ODR-3, GPA-2, GPA-5, and GPA-6 and a Ca2+/MAPK pathway involving the Ca2+ channel subunit UNC-36, the CaMKII UNC-43, and the MAPKK kinase NSY-1 induce strong str-2 expression. Cell-specific rescue experiments suggest that ODR-3 and the Ca2+/MAPK genes function in the AWC neurons, but that GPA-5 and GPA-6 function in the AWA and ADL neurons, respectively. In Dauer larvae, the same network of genes promotes strong str-2 expression in the ASI neurons, but ODR-3 functions in AWB and ASH and GPA-6 in AWB. Our results reveal a complex signaling network, encompassing signals from multiple cells, that controls the level of receptor gene expression at different developmental stages.

The chemoreceptor superfamily in the honey bee, Apis mellifera: expansion of the odorant, but not gustatory, receptor family

The honey bee genome sequence reveals a remarkable expansion of the insect odorant receptor (Or) family relative to the repertoires of the flies Drosophila melanogaster and Anopheles gambiae, which have 62 and 79 Ors respectively. A total of 170 Or genes were annotated in the bee, of which seven are pseudogenes. These constitute five bee-specific subfamilies in an insect Or family tree, one of which has expanded to a total of 157 genes encoding proteins with 15%-99% amino acid identity. Most of the Or genes are in tandem arrays, including one with 60 genes. This bee-specific expansion of the Or repertoire presumably underlies their remarkable olfactory abilities, including perception of several pheromone blends, kin recognition signals, and diverse floral odors. The number of Apis mellifera Ors is approximately equal to the number of glomeruli in the bee antennal lobe (160-170), consistent with a general one-receptor/one-neuron/one-glomerulus relationship. The bee genome encodes just 10 gustatory receptors (Grs) compared with the D. melanogaster and A. gambiae repertoires of 68 and 76 Grs, respectively. A lack of Gr gene family expansion primarily accounts for this difference. A nurturing hive environment and a mutualistic relationship with plants may explain the lack of Gr family expansion. The Or family is the most dramatic example of gene family expansion in the bee genome, and characterizing their caste- and sex-specific gene expression may provide clues to their specific roles in detection of pheromone, kin, and floral odors.

Ribosomal DNA in the grasshopper Podisma pedestris: escape from concerted evolution

Eukaryote nuclear ribosomal DNA (rDNA) typically exhibits strong concerted evolution: a pattern in which several hundred rDNA sequences within any one species show little or no genetic diversity, whereas the sequences of different species diverge. We report a markedly different pattern in the genome of the grasshopper Podisma pedestris. Single individuals contain several highly divergent ribosomal DNA groups. Analysis of the magnitude of divergence indicates that these groups have coexisted in the Podisma lineage for at least 11 million years. There are two putatively functional groups, each estimated to be at least 4 million years old, and several pseudogene groups, many of which are transcribed. Southern hybridization and real-time PCR experiments show that only one of the putatively functional types occurs at high copy number. However, this group is scarcely amplified under standard PCR conditions, which means that phylogenetic inference on the basis of standard PCR would be severely distorted. The analysis suggests that concerted evolution has been remarkably ineffective in P. pedestris. We propose that this outcome may be related to the species' exceptionally large genome and the associated low rate of deletion per base pair, which may allow pseudogenes to persist.

Evolutionary genomics of the recently duplicated amphioxus Hairy genes

Amphioxus Hairy genes have gone through a number of lineage-specific duplications, resulting in eight members, some of which are differentially expressed in the embryo. In order to gain insights into the evolution and function of this gene family we have compared their genomic structure and searched for conserved non-coding sequence elements. We have found that introns have been lost independently from these genes at least twice and after the duplication events. By carrying out phylogenetic footprinting between paralogues expressed in the embryo, we have found a differential distribution of conserved elements that could explain the limited overlap in expression patterns of Hairy genes in the amphioxus embryo. Furthermore, clustering of RBP-Jk binding sites in these conserved elements suggests that amphioxus Hairy genes are downstream targets of the Notch signaling pathway, as occurs in vertebrates. All of this evidence suggests that amphioxus Hairy genes have gone through a process of subfunctionalization shortly after their duplication, representing an extreme and rapid case of the duplication-degeneration-complementation model.

Characterization of human gene encoding SLA/LP autoantigen and its conserved homologs in mouse, fish, fly, and worm

AIM: To approach the elusive function of the SLA/LP molecule, we have characterized genomic organization and conservation of the major antigenic and functional properties of the SLA/LP molecule in various species.

METHODS: By means of computational biology, we have characterized the complete SLA/LP gene, mRNA and deduced protein sequences in man, mouse, zebrafish, fly, and worm.

RESULTS: The human SLA/LP gene sequence of approximately 39 kb, which maps to chromosome 4p15.2, is organized in 11 exons, of which 10 or 11 are translated, depending on the splice variant. Homologous molecules were identified in several biological model organisms. The various homologous protein sequences showed a high degree of similarity or homology, notably at those residues that are of functional importance. The only domain of the human protein sequence that lacks significant homology with homologous sequences is the major antigenic epitope recognized by autoantibodies from autoimmune hepatitis (AIH) patients.

CONCLUSION: The SLA/LP molecule and its functionally relevant residues have been highly conserved throughout the evolution, suggesting an indispensable function of the molecule. The finding that the only non-conserved domain is the dominant antigenic epitope of the human SLA/LP sequence, suggests that SLA/LP autoimmunity is autoantigen-driven rather than being driven by molecular mimicry.

A general model of G protein-coupled receptor sequences and its application to detect remote homologs

G protein-coupled receptors (GPCRs) constitute a large superfamily involved in various types of signal transduction pathways triggered by hormones, odorants, peptides, proteins, and other types of ligands. The superfamily is so diverse that many members lack sequence similarity, although they all span the cell membrane seven times with an extracellular N and a cytosolic C terminus. We analyzed a divergent set of GPCRs and found distinct loop length patterns and differences in amino acid composition between cytosolic loops, extracellular loops, and membrane regions. We configured GPCRHMM, a hidden Markov model, to fit those features and trained it on a large dataset representing the entire superfamily. GPCRHMM was benchmarked to profile HMMs and generic transmembrane detectors on sets of known GPCRs and non-GPCRs. In a cross-validation procedure, profile HMMs produced an error rate nearly twice as high as GPCRHMM. In a sensitivity-selectivity test, GPCRHMM's sensitivity was about 15% higher than that of the best transmembrane predictors, at comparable false positive rates. We used GPCRHMM to search for novel members of the GPCR superfamily in five proteomes. All in all we detected 120 sequences that lacked annotation and are potentially novel GPCRs. Out of those 102 were found in Caenorhabditis elegans, four in human, and seven in mouse. Many predictions (65) belonged to Pfam domains of unknown function. GPCRHMM strongly rejected a family of arthropod-specific odorant receptors believed to be GPCRs. A detailed analysis showed that these sequences are indeed very different from other GPCRs. GPCRHMM is available at http://gpcrhmm.cgb.ki.se.

Concerted evolution of two novel protein families in Caenorhabditis species

Among a large number of homologous gene clusters in C. elegans, two gene families that appear to undergo concerted evolution were studied in detail. Both gene families are nematode specific and encode small secreted proteins of unknown function. For both families in three Caenorhabditis species, concerted groups of genes are characterized by close genomic proximity and by genes in inverted orientation. The rate of protein evolution in one of the two families could be calibrated by comparison with a closely related nonconcerted singleton gene with one-to-one orthologs in all three species. This comparison suggests that protein evolution in concerted gene clusters is two- to sevenfold accelerated. A broader survey of clustered gene families, focused on adjacent inverted gene pairs, identified an additional seven families in which concerted evolution probably occurs. All nine identified families encode relatively small proteins, eight of them encode putative secreted proteins, and most of these have very unusual amino acid composition or sequence. I speculate that these genes encode rapidly evolving antimicrobial peptides.

The excess of 5' introns in eukaryotic genomes

In this work, 21 completely sequenced eukaryotic genomes were analyzed using an intragene comparison approach. We found that all of these genomes show a significant 5'-biased distribution of introns of protein-coding genes. Our findings are different from previous studies based on the intergene method, where introns are biased towards the 5' end of genes only in intron-poor genomes, but are evenly distributed in intron-rich genomes. In addition, by analyzing the patterns of intron distribution of a set of well-compiled housekeeping genes from human and their respective orthologs identified by a bidirectional best BLAST hit method from the other genomes, we found that the trend of 5'-biased intron positions of the set of housekeeping genes for each genome is much more skewed than that of all genes of the same genome, and rarely if any of the housekeeping genes examined have an extremely 3'-biased position distribution in which all introns of a gene are located only at the 3' portion of the gene. The most parsimonious explanation for our findings may be the model in which intron loss is caused by homologous recombination between the genomic copy of a gene and a reverse transcriptase product of a spliced mRNA.

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

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

Explosive lineage-specific expansion of the orphan nuclear receptor HNF4 in nematodes

The nuclear receptor superfamily expanded in at least two episodes: one early in metazoan evolution, the second within the vertebrate lineage. An exception to this pattern is the genome of the nematode Caenorhabditis elegans, which encodes more than 270 nuclear receptors, most of them highly divergent. We generated 128 cDNA sequences for 76 C. elegans nuclear receptors, confirming that these are active genes. Among these numerous receptors are 13 orthologues of nuclear receptors found in arthropods and/or vertebrates. We show that the supplementary nuclear receptors (supnrs) originated from an explosive burst of duplications of a unique orphan receptor, HNF4. This origin has specific implications for the role of ligand binding in the function and evolution of the nematode supplementary nuclear receptors. Moreover, the supplementary nuclear receptors include a group of very rapidly evolving genes found primarily on chromosome V. We propose a model of lineage-specific duplications from a chromosome on which duplication and substitution rates are highly increased. Our results provide a framework to study nuclear receptors in nematodes, as well as to consider the functional and evolutionary consequences of lineage-specific duplications.

Characterization of novel GPCR gene coding locus in amphioxus genome: gene structure, expression, and phylogenetic analysis with implications for its involvement in chemoreception

Chemosensation is the primary sensory modality in almost all metazoans. The vertebrate olfactory receptor genes exist as tandem clusters in the genome, so that identifying their evolutionary origin would be useful for understanding the expansion of the sensory world in relation to a large-scale genomic duplication event in a lineage leading to the vertebrates. In this study, I characterized a novel GPCR (G-protein-coupled receptor) gene-coding locus from the amphioxus genome. The genomic DNA contains an intronless ORF whose deduced amino acid sequence encodes a seven-transmembrane protein with some amino acid residues characteristic of vertebrate olfactory receptors (ORs). Surveying counterparts in the Ciona intestinalis (Asidiacea, Urochordata) genome by querying BLAST programs against the Ciona genomic DNA sequence database resulted in the identification of a remotely related gene. In situ hybridization analysis labeled primary sensory neurons in the rostral epithelium of amphioxus adults. Based on these findings, together with comparison of the developmental gene expression between amphioxus and vertebrates, I postulate that chemoreceptive primary sensory neurons in the rostrum are an ancient cell population traceable at least as far back in phylogeny as the common ancestor of amphioxus and vertebrates.

Functional constraint and divergence in the G protein family in Caenorhabditis elegans and Caenorhabditis briggsae

Part of the challenge of the post-genomic world is to identify functional elements within the wide array of information generated by genome sequencing. Although cross-species comparisons and investigation of rates of sequence divergence are an efficient approach, the relationship between sequence divergence and functional conservation is not clear. Here, we use a comparative approach to examine questions of evolutionary rates and conserved function within the guanine nucleotide-binding protein (G protein) gene family in nematodes of the genus Caenorhabditis. In particular, we show that, in cases where the Caenorhabditis elegans ortholog shows a loss-of-function phenotype, G protein genes of C. elegans and Caenorhabditis briggsae diverge on average three times more slowly than G protein genes that do not exhibit any phenotype when mutated in C. elegans, suggesting that genes with loss of function phenotypes are subject to stronger selective constraints in relation to their function in both species. Our results also indicate that selection is as strong on G proteins involved in environmental perception as it is on those controlling other important processes. Finally, using phylogenetic footprinting, we identify a conserved non-coding motif present in multiple copies in the genomes of four species of Caenorhabditis. The presence of this motif in the same intron in the gpa-1 genes of C. elegans, C. briggsae and Caenorhabditis remanei suggests that it plays a role in the regulation of gpa-1, as well as other loci.

Rates of intron loss and gain: implications for early eukaryotic evolution

We study the intron-exon structures of 684 groups of orthologs from seven diverse eukaryotic genomes and provide maximum likelihood estimates for rates and numbers of intron losses and gains in these same genes for a variety of lineages. Rates of intron loss vary from approximately 2 x 10(-9) to 2 x 10(-10) per year. Rates of gain vary from 6 x 10(-13) to 4 x 10(-12) per possible intron insertion site per year. There is an inverse correspondence between rates of intron loss and gain, leading to a 20-fold variation among lineages in the ratio of the rates of the two processes. The observed rates of intron gain are insufficient to explain the large number of introns estimated to have been present in the plant-animal ancestor, suggesting that introns present in early eukaryotes may have been created by a fundamentally different process than more recently gained introns.

Adaptive evolution in the SRZ chemoreceptor families of Caenorhabditis elegans and Caenorhabditis briggsae

We investigated the possibility of positive selection acting on members of the putative seven-pass chemoreceptor superfamily in Caenorhabditis elegans, which comprises approximately 1,300 genes encoding seven-pass G protein-coupled receptors (GPCRs). Using a maximum-likelihood approach, we conducted statistical tests for evidence of codon sites where the ratio of nonsynonymous substitutions per nonsynonymous site to synonymous substitutions per synonymous site (d(N)/d(S)) was >1. Evidence for positive selection was found only for the srz family, about which virtually nothing specific is known. We extended the annotation of the srz gene family, establishing gene models for 60 srz genes in C. elegans and 28 srz genes in Caenorhabditis briggsae. d(N)/d(S) ratios varied dramatically in different regions of the SRZ proteins, peaking in predicted extracellular regions. These regions included 23 sites where evidence of positive selection was highly significant, corresponding remarkably well with regions implicated in ligand binding in other GPCR family members. We interpret these results as indicating that the srz family is under positive selection, probably driven by ligand binding.

Complex early genes

We use the pattern of intron conservation in 684 groups of orthologs from seven fully sequenced eukaryotic genomes to provide maximum likelihood estimates of the number of introns present in the same orthologs in various eukaryotic ancestors. We find: (i) intron density in the plant-animal ancestor was high, perhaps two-thirds that of humans and three times that of Drosophila; and (ii) intron density in the ancestral bilateran was also high, equaling that of humans and four times that of Drosophila. We further find that modern introns are generally very old, with two-thirds of modern bilateran introns dating to the ancestral bilateran and two-fifths of modern plant, animal, and fungus introns dating to the plant-animal ancestor. Intron losses outnumber gains over a large range of eukaryotic lineages. These results show that early eukaryotic gene structures were very complex, and that simplification, not embellishment, has dominated subsequent evolution.