Mode and tempo of molecular evolution in the nematode caenorhabditis: cytochrome oxidase II and calmodulin sequences

Eimeria zuernii (Eimeriidae: Coccidia): mitochondrial genome and genetic diversity in the Chinese yak

BACKGROUND

Coccidiosis caused by Eimeria zuernii (Eimeriidae: Coccidia) represents a significant economic threat to the bovine industry. Understanding the evolutionary and genetic biology of E. zuernii can assist in new interaction developments for the prevention and control of this protozoosis.

METHODS

We defined the evolutionary and genetic characteristics of E. zuernii by sequencing the complete mitogenome and analyzing the genetic diversity and population structure of 51 isolates collected from eight yak breeding parks in China.

RESULTS

The 6176-bp mitogenome of E. zuernii was linear and encoded typical mitochondrial contents of apicomplexan parasites, including three protein-coding genes [PCGs; cytochrome c oxidase subunits I and III (cox1 and cox3), and cytochrome b (cytb)], seven fragmented small subunit (SSU) and 12 fragmented large subunit (LSU) rRNAs. Genome-wide comparative and evolutionary analyses showed cytb and cox3 to be the most and least conserved Eimeria PCGs, respectively, and placed E. zuernii more closely related to Eimeria mephitidis than other Eimeria species. Furthermore, cox1-based genetic structure defined 24 haplotypes of E. zuernii with high haplotype diversities and low nucleotide diversities across eight geographic populations, supporting a low genetic structure and rapid evolutionary rate as well as a previous expansion event among E. zuernii populations.

CONCLUSIONS

To our knowledge, this is the first study presenting the phylogeny, genetic diversity, and population structure of the yak E. zuernii, and such information, together with its mitogenomic data, should contribute to a better understanding of the genetic and evolutionary biological studies of apicomplexan parasites in bovines.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans

Human mitochondria contain a high copy number, maternally transmitted genome (mtDNA) that encodes 13 proteins required for oxidative phosphorylation. Heteroplasmy arises when multiple mtDNA variants co-exist in an individual and can exhibit complex dynamics in disease and in aging. As all proteins involved in mtDNA replication and maintenance are nuclear-encoded, heteroplasmy levels can, in principle, be under nuclear genetic control, however this has never been shown in humans. Here, we develop algorithms to quantify mtDNA copy number (mtCN) and heteroplasmy levels using blood-derived whole genome sequences from 274,832 individuals of diverse ancestry and perform GWAS to identify nuclear loci controlling these traits. After careful correction for blood cell composition, we observe that mtCN declines linearly with age and is associated with 92 independent nuclear genetic loci. We find that nearly every individual carries heteroplasmic variants that obey two key patterns: (1) heteroplasmic single nucleotide variants are somatic mutations that accumulate sharply after age 70, while (2) heteroplasmic indels are maternally transmitted as mtDNA mixtures with resulting levels influenced by 42 independent nuclear loci involved in mtDNA replication, maintenance, and novel pathways. These nuclear loci do not appear to act by mtDNA mutagenesis, but rather, likely act by conferring a replicative advantage to specific mtDNA molecules. As an illustrative example, the most common heteroplasmy we identify is a length variant carried by >50% of humans at position m.302 within a G-quadruplex known to serve as a replication switch. We find that this heteroplasmic variant exerts cis -acting genetic control over mtDNA abundance and is itself under trans -acting genetic control of nuclear loci encoding protein components of this regulatory switch. Our study showcases how nuclear haplotype can privilege the replication of specific mtDNA molecules to shape mtCN and heteroplasmy dynamics in the human population.

High Mitochondrial Genome Diversity and Intricate Population Structure of Bursaphelenchus xylophilus in Kyushu, Japan

Mitogenomic diversity and genetic population structure of the pinewood nematode (PWN) Bursaphelenchus xylophilus inhabiting Kyushu, Japan were analyzed. A method for performing long PCR using single nematodes and sequencing nematode mitochondrial genomes individually is presented here. About 8 kb (∼55%) of the complete mitochondrial genome was successfully obtained from 285 individuals collected from 12 populations. The 158 single nucleotide polymorphisms detected corresponded to 30 haplotypes, clearly classified into two clades. Haplotype diversity was 0.83, evidencing a remarkable high diversity within Kyushu. The high genetic differentiation among the 12 populations (0.331) might be due to past invasion and expansion routes of PWN in northeastern and southeastern Kyushu. The distinct genetic composition of populations within the northwestern, central western, and southwestern Kyushu seems to be mostly related to the extinction of pine forests and long-range migration of PWN due to human activity. Overall, direct long PCR and sequencing of single nematode individuals are effective methods for investigating mitochondrial polymorphisms, and these are effective tools for PWN population genetics and other intraspecific studies.

Genetic characterisation and phylogenetic status of whipworms (Trichuris spp.) from captive non-human primates in China, determined by nuclear and mitochondrial sequencing

Background

Whipworms (Nematoda: Trichuridae), among the most common soil-transmitted helminths (STHs), can cause the socioeconomically important disease trichuriasis in various mammalian hosts including humans and non-human primates. For many years, Trichuris from non-human primates has been assigned to the same species as the one infecting humans Trichuris trichiura. More recently, several molecular reports challenged this assumption following recognition of a Trichuris species complex observed in humans and non-human primates. A refined concept for species limits within Trichuris contributes to an understanding of diversity and the potential (zoonotic) transmission among humans and non-human primates. In this study, we expanded previous investigations by exploring the diversity of Trichuris among eight primates including three Asian autochthonous species (i.e. Rhinopithecus roxellana, Rhinopithecus bieti and Nomascus leucogenys). Species-level identification, whether novel or assignable to known lineages of Trichuris, was based on analyses of nuclear internal transcribed spacers (ITS) and mitochondrial cytochrome c oxidase subunit 1 (cox1) genes.

Results

In total, seven genetically distinct subgroups of whipworms were determined to be present among the primates sampled. Most Trichuris lineages, including Subgroups 1, 1’, 3, 5 and 6, showed a broad host range and were not restricted to particular primate species; in addition to T. trichiura, a complex of Trichuris species was shown infecting primates. Furthermore, it was assumed that Trichuris spp. from either N. leucogenys and P. hamadryas or R. roxellana and R. bieti, respectively, were conspecific. Each pair was indicated to be a discrete lineage of Trichuris, designated, respectively, as Subgroups 1 or 1’ and 2, based on integrated genetic and phylogenetic evidence.

Conclusion

These results emphasise that the taxonomy and genetic variations of Trichuris are more complicated than previously acknowledged. These cumulative molecular and phylogenetic data provide a better understanding of the taxonomy, genetics and evolutionary biology of the whipworms.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3100-5) contains supplementary material, which is available to authorized users.

Large Variation in the Ratio of Mitochondrial to Nuclear Mutation Rate across Animals: Implications for Genetic Diversity and the Use of Mitochondrial DNA as a Molecular Marker

It is commonly assumed that mitochondrial DNA (mtDNA) evolves at a faster rate than nuclear DNA (nuDNA) in animals. This has contributed to the popularity of mtDNA as a molecular marker in evolutionary studies. Analyzing 121 multilocus data sets and four phylogenomic data sets encompassing 4,676 species of animals, we demonstrate that the ratio of mitochondrial over nuclear mutation rate is highly variable among animal taxa. In nonvertebrates, such as insects and arachnids, the ratio of mtDNA over nuDNA mutation rate varies between 2 and 6, whereas it is above 20, on average, in vertebrates such as scaled reptiles and birds. Interestingly, this variation is sufficient to explain the previous report of a similar level of mitochondrial polymorphism, on average, between vertebrates and nonvertebrates, which was originally interpreted as reflecting the effect of pervasive positive selection. Our analysis rather indicates that the among-phyla homogeneity in within-species mtDNA diversity is due to a negative correlation between mtDNA per-generation mutation rate and effective population size, irrespective of the action of natural selection. Finally, we explore the variation in the absolute per-year mutation rate of both mtDNA and nuDNA using a reduced data set for which fossil calibration is available, and discuss the potential determinants of mutation rate variation across genomes and taxa. This study has important implications regarding DNA-based identification methods in predicting that mtDNA barcoding should be less reliable in nonvertebrates than in vertebrates.

A Model for Evolutionary Ecology of Disease: The Case for Caenorhabditis Nematodes and Their Natural Parasites

Many of the outstanding questions in disease ecology and evolution call for combining observation of natural host-parasite populations with experimental dissection of interactions in the field and the laboratory. The "rewilding" of model systems holds great promise for this endeavor. Here, we highlight the potential for development of the nematode Caenorhabditis elegans and its close relatives as a model for the study of disease ecology and evolution. This powerful laboratory model was disassociated from its natural habitat in the 1960s. Today, studies are uncovering that lost natural history, with several natural parasites described since 2008. Studies of these natural Caenorhabditis-parasite interactions can reap the benefits of the vast array of experimental and genetic tools developed for this laboratory model. In this review, we introduce the natural parasites of C. elegans characterized thus far and discuss resources available to study them, including experimental (co)evolution, cryopreservation, behavioral assays, and genomic tools. Throughout, we present avenues of research that are interesting and feasible to address with caenorhabditid nematodes and their natural parasites, ranging from the maintenance of outcrossing to the community dynamics of host-associated microbes. In combining natural relevance with the experimental power of a laboratory supermodel, these fledgling host-parasite systems can take on fundamental questions in evolutionary ecology of disease.

Expression of calmodulin and myosin light chain kinase during larval settlement of the Barnacle Balanus amphitrite

Barnacles are one of the most common organisms in intertidal areas. Their life cycle includes seven free-swimming larval stages and sessile juvenile and adult stages. The transition from the swimming to the sessile stages, referred to as larval settlement, is crucial for their survivor success and subsequent population distribution. In this study, we focused on the involvement of calmodulin (CaM) and its binding proteins in the larval settlement of the barnacle, Balanus ( = Amphibalanus) amphitrite. The full length of CaM gene was cloned from stage II nauplii of B. amphitrite (referred to as Ba-CaM), encoding 149 amino acid residues that share a high similarity with published CaMs in other organisms. Quantitative real-time PCR showed that Ba-CaM was highly expressed in cyprids, the stage at which swimming larvae are competent to attach and undergo metamorphosis. In situ hybridization revealed that the expressed Ba-CaM gene was localized in compound eyes, posterior ganglion and cement glands, all of which may have essential functions during larval settlement. Larval settlement assays showed that both the CaM inhibitor compound 48/80 and the CaM-dependent myosin light chain kinase (MLCK) inhibitor ML-7 effectively blocked barnacle larval settlement, whereas Ca(2+)/CaM-dependent kinase II (CaMKII) inhibitors did not show any clear effects. The subsequent real-time PCR assay showed a higher expression level of Ba-MLCK gene in larval stages than in adults, suggesting an important role of Ba-MLCK gene in larval development and competency. Overall, the results suggest that CaM and CaM-dependent MLCK function during larval settlement of B. amphitrite.

Influence of finite-sites mutation, population subdivision and sampling schemes on patterns of nucleotide polymorphism for species with molecular hyperdiversity

Molecular hyperdiversity has been documented in viruses, prokaryotes and eukaryotes. Such organisms undermine the assumptions of the infinite-sites mutational model, because multiple mutational events at a site comprise a non-negligible portion of polymorphisms. Moreover, different sampling schemes of individuals from species with subdivided populations can profoundly influence resulting patterns and interpretations of molecular variation. Inspired by molecular hyperdiversity in the nematode Caenorhabditis sp. 5, which exhibits average pairwise differences among synonymous sites of >5% as well as modest population structure, we investigated via coalescent simulation the joint effects of a finite-sites mutation (FSM) process and population subdivision on the variant frequency spectrum. From many demes interconnected through a stepping-stone migration model, we constructed local samples from a single deme, pooled samples from several demes and scattered samples of a single individual from numerous demes. Compared with a single panmictic population at equilibrium, we find that high population mutation rates induce a deficit of rare variants (positive Tajima's D) under a FSM model. Population structure also induces such a skew for local samples when migration is high and for pooled samples when migration is low. Contrasts of sampling schemes for C. sp. 5 imply high mutational input coupled with high migration. We propose that joint analysis of local, pooled and scattered samples for species with subdivided populations provides a means of improving inference of demographic history, by virtue of the partially distinct patterns of polymorphism that manifest when sequences are analyzed according to differing sampling schemes.

Shifting patterns of natural variation in the nuclear genome of caenorhabditis elegans

Background

Genome wide analysis of variation within a species can reveal the evolution of fundamental biological processes such as mutation, recombination, and natural selection. We compare genome wide sequence differences between two independent isolates of the nematode Caenorhabditis elegans (CB4856 and CB4858) and the reference genome (N2).

Results

The base substitution pattern when comparing N2 against CB4858 reveals a transition over transversion bias (1.32:1) that is not present in CB4856. In CB4856, there is a significant bias in the direction of base substitution. The frequency of A or T bases in N2 that are G or C bases in CB4856 outnumber the opposite frequencies for transitions as well as transversions. These differences were not observed in the N2/CB4858 comparison. Similarly, we observed a strong bias for deletions over insertions in CB4856 (1.44: 1) that is not present in CB4858. In both CB4856 and CB4858, there is a significant correlation between SNP rate and recombination rate on the autosomes but not on the X chromosome. Furthermore, we identified numerous significant hotspots of variation in the CB4856-N2 comparison.

In both CB4856 and CB4858, based on a measure of the strength of selection (k_a/k_s), all the chromosomes are under negative selection and in CB4856, there is no difference in the strength of natural selection in either the autosomes versus X or between any of the chromosomes. By contrast, in CB4858, k_a/k_s values are smaller in the autosomes than in the X chromosome. In addition, in CB4858, k_a/k_s values differ between chromosomes.

Conclusions

The clear bias of deletions over insertions in CB4856 suggests that either the CB4856 genome is becoming smaller or the N2 genome is getting larger. We hypothesize the hotspots found represent alleles that are shared between CB4856 and CB4858 but not N2. Because the k_a/k_s ratio in the X chromosome is higher than the autosomes on average in CB4858, purifying selection is reduced on the X chromosome.

Complete mitochondrial genomes of Baylisascaris schroederi, Baylisascaris ailuri and Baylisascaris transfuga from giant panda, red panda and polar bear

Roundworms of the genus Baylisascaris are the most common parasitic nematodes of the intestinal tracts of wild mammals, and most of them have significant impacts in veterinary and public health. Mitochondrial (mt) genomes provide a foundation for studying epidemiology and ecology of these parasites and therefore may be used to assist in the control of Baylisascariasis. Here, we determined the complete sequences of mtDNAs for Baylisascaris schroederi, Baylisascaris ailuri and Baylisascaris transfuga, with 14,778 bp, 14,657 bp and 14,898 bp in size, respectively. Each mtDNA encodes 12 protein-coding genes, 22 transfer RNAs and 2 ribosomal RNAs, typical for other chromadorean nematodes. The gene arrangements for the three Baylisascaris species are the same as those of the Ascaridata species, but radically different from those of the Spirurida species. Phylogenetic analysis based on concatenated amino acid sequences of 12 protein-coding genes from nine nematode species indicated that the three Baylisascaris species are more closely related to Ascaris suum than to the three Toxocara species (Toxocara canis, Toxocara cati and Toxocara malaysiensis) and Anisakis simplex, and that B. ailuri is more closely related to B. transfuga than to B. schroeder. The determination of the complete mt genome sequences for these three Baylisascaris species (the first members of the genus Baylisascaris ever sequenced) is of importance in refining the phylogenetic relationships within the order Ascaridida, and provides new molecular data for population genetic, systematic, epidemiological and ecological studies of parasitic nematodes of socio-economic importance in wildlife.

Extremely high molecular diversity within the East Asian nematode Caenorhabditis sp. 5

Most relatives of the self-fertilizing hermaphroditic nematode model organism Caenorhabditis elegans reproduce via obligate outbreeding between males and females, which also represents the ancestral mode of reproduction within the genus. However, little is known about the scope of genetic diversity and differentiation within such gonochoristic species, especially those found outside of temperate Europe and North America. It is critical to understand the evolutionary processes operating in these species to provide a framework for deciphering the evolution of hermaphroditism and a baseline for the application of outcrossing Caenorhabditis to problems in evolutionary genetics. Here, we investigate for the first time molecular sequence variation for Caenorhabditis sp. 5, a species found commonly in eastern Asia. We identify enormous levels of standing genetic variation that approach the levels observed in the marine broadcast-spawning sea squirt, Ciona savignyi. Although we document significant isolation by distance, we demonstrate that the high polymorphism within C. sp. 5 is not because of strong differentiation among populations or to the presence of cryptic species. These findings illustrate that molecular population genetic approaches to studying obligately outbreeding species of Caenorhabditis will prove powerful in identifying and characterizing functionally and evolutionarily important features of the genome.

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.

Recombinational landscape and population genomics of Caenorhabditis elegans

Recombination rate and linkage disequilibrium, the latter a function of population genomic processes, are the critical parameters for mapping by linkage and association, and their patterns in Caenorhabditis elegans are poorly understood. We performed high-density SNP genotyping on a large panel of recombinant inbred advanced intercross lines (RIAILs) of C. elegans to characterize the landscape of recombination and, on a panel of wild strains, to characterize population genomic patterns. We confirmed that C. elegans autosomes exhibit discrete domains of nearly constant recombination rate, and we show, for the first time, that the pattern holds for the X chromosome as well. The terminal domains of each chromosome, spanning about 7% of the genome, exhibit effectively no recombination. The RIAILs exhibit a 5.3-fold expansion of the genetic map. With median marker spacing of 61 kb, they are a powerful resource for mapping quantitative trait loci in C. elegans. Among 125 wild isolates, we identified only 41 distinct haplotypes. The patterns of genotypic similarity suggest that some presumed wild strains are laboratory contaminants. The Hawaiian strain, CB4856, exhibits genetic isolation from the remainder of the global population, whose members exhibit ample evidence of intercrossing and recombining. The population effective recombination rate, estimated from the pattern of linkage disequilibrium, is correlated with the estimated meiotic recombination rate, but its magnitude implies that the effective rate of outcrossing is extremely low, corroborating reports of selection against recombinant genotypes. Despite the low population, effective recombination rate and extensive linkage disequilibrium among chromosomes, which are techniques that account for background levels of genomic similarity, permit association mapping in wild C. elegans strains.

C. elegans is a model system for diverse fields of biology, but its ability to serve as a model for quantitative trait gene mapping depends on its recombination rate in the laboratory and in nature. The latter is a function of how worms mate and migrate in the wild. We examined the patterns of recombination in a population that we put through thousands of meioses in the laboratory and in a collection of strains isolated from nature. The data suggest that meiotic recombination rate is highly regular in worms, with discrete domains whose boundaries we identify. The pattern in natural strains suggests that population structure, population size, outcrossing rate, and selection combine to suppress the overall effects of recombination. Moreover, some “wild” strains appear to be laboratory contaminants. Nevertheless, the history of recombination in wild worms is sufficient to permit correlations between genotype and phenotype to pinpoint the loci responsible for phenotypic variation.

Distinct patterns of genetic variation in Pristionchus pacificus and Caenorhabditis elegans, two partially selfing nematodes with cosmopolitan distribution

Hermaphroditism has evolved several times independently in nematodes. The model organism Caenorhabditis elegans and Pristionchus pacificus are self-fertile hermaphrodites with rare facultative males. Both species are members of different families: C. elegans belongs to the Rhabditidae and P. pacificus to the Diplogastridae. Also, both species differ in their ecology: C. elegans is a soil-dwelling nematode that is often found in compost heaps. In contrast, field studies in Europe and North America indicate that Pristionchus nematodes are closely associated with scarab beetles. In C. elegans, several recent studies have found low genetic diversity and rare out-crossing events. Little is known about diversity levels and population structure in free-living hermaphroditic nematodes outside the genus Caenorhabditis. Taking a comparative approach, we analyse patterns of molecular diversity and linkage disequilibrium in 18 strains of P. pacificus from eight countries and four continents. Mitochondrial sequence data of P. pacificus isolates reveal a substantially higher genetic diversity on a global scale when compared to C. elegans. A mitochondrial-derived hermaphrodite phylogeny shows little geographic structuring, indicating several worldwide dispersal events. Amplified fragment length polymorphism and single strand conformation polymorphism analyses demonstrate a high degree of genome-wide linkage disequilibrium, which also extends to the mitochondrial genome. Together, these findings indicate distinct patterns of genetic variation of the two species. The low level of genetic diversity observed in C. elegans might reflect a recent human-associated dispersal, whereas the P. pacificus diversity might reflect a long-lasting and ongoing insect association. Thus, despite similar lifestyle characteristics in the laboratory, the reproductive mode of hermaphroditism with rare facultative males can result in distinct genetic variability patterns in different ecological settings.

Population genetics of anthelmintic resistance in parasitic nematodes

A key aim of anthelmintic resistance research is to identify molecular markers that could form the basis of sensitive and accurate diagnostic tests. These would provide powerful tools to study the origin and spread of anthelmintic resistance in the field and to monitor strategies aimed at preventing and managing resistance. Molecular markers could also form the basis of routine diagnostic tests for use in surveillance and clinical veterinary practice. Much of the research conducted to date has focused on the investigation of possible associations of particular candidate genes with the resistance phenotype. In the future, as full parasite genome sequences become available, there will be an opportunity to apply genome-wide approaches to identify the genetic loci that underlie anthelmintic resistance. Both the interpretation of candidate gene studies and the application of genome-wide approaches require a good understanding of the genetics and population biology of the relevant parasites as well as knowledge of how resistance mutations arise and are selected in populations. Unfortunately, much of this information is lacking for parasitic nematodes. This review deals with a number of aspects of genetics and population biology that are pertinent to these issues. We discuss the possible origins of resistance mutations and the likely effects of subsequent selection on the genetic variation at the resistance-conferring locus. We also review some of the experimental approaches that have been used to test associations between candidate genes and anthelmintic resistance phenotypes and highlight implications for future genome-wide studies.

The frataxin-encoding operon of Caenorhabditis elegans shows complex structure and regulation

The genome of the nematode Caenorhabditis elegans is unusual among eukaryotes, in that it contains operons. Approximately 15% of genes in the worm are clustered into groups of between two and eight genes, which are under the control of shared regulatory sequences. Polycistronic transcripts from such operons are trans-spliced, during transcription, to produce mature monocistronic messengers. The C. elegans frataxin gene, frh-1, is encoded in the operon CEOP2232. This is one of the largest operons identified thus far in the C. elegans genome. Here we describe in detail the structure of all of the coding units within this operon. The operon is composed of eight genes of a diverse nature, organized in a complex structure. We have produced transgenic strains carrying fusions between gfp and a number of genes from the operon. These constructs show complex differential expression patterns that suggest the presence of internal promoters and regulatory sequences in the operon. This organization would permit both coordinated expression and differential expression of the components of the CEOP2232 operon. The heterogeneity of the genes, and their complex expression patterns, suggests that the clustering of CEOP2232 is not due to a need for synchronized expression of genes involved in the same physiological pathway.

An improved molecular phylogeny of the Nematoda with special emphasis on marine taxa

Phylogenetic reconstructions of relations within the phylum Nematoda are inherently difficult but have been advanced with the introduction of large-scale molecular-based techniques. However, the most recent revisions were heavily biased towards terrestrial and parasitic species and greater representation of clades containing marine species (e.g. Araeolaimida, Chromadorida, Desmodorida, Desmoscolecida, Enoplida, and Monhysterida) is needed for accurate coverage of known taxonomic diversity. We now add small subunit ribosomal DNA (SSU rDNA) sequences for 100 previously un-sequenced species of nematodes, including 46 marine taxa. SSU rDNA sequences for >200 taxa have been analysed based on Bayesian inference and LogDet-transformed distances. The resulting phylogenies provide support for (i) the re-classification of the Secernentea as the order Rhabditida that derived from a common ancestor of chromadorean orders Araeolaimida, Chromadorida, Desmodorida, Desmoscolecida, and Monhysterida and (ii) the position of Bunonema close to the Diplogasteroidea in the Rhabditina. Other, previously controversial relationships can now be resolved more clearly: (a) Alaimus, Campydora, and Trischistoma belong in the Enoplida, (b) Isolaimium is placed basally to a big clade containing the Axonolaimidae, Plectidae, and Rhabditida, (c) Xyzzors belongs in the Desmodoridae, (d) Comesomatidae and Cyartonema belongs in the Monhysterida, (e) Globodera belongs in the Hoplolaimidae and (f) Paratylenchus dianeae belongs in the Criconematoidea. However, the SSU gene did not provide significant support for the class Chromadoria or clear evidence for the relationship between the three classes, Enoplia, Dorylaimia, and Chromadoria. Furthermore, across the whole phylum, the phylogenetically informative characters of the SSU gene are not informative in a parsimony analysis, highlighting the short-comings of the parsimony method for large-scale phylogenetic modelling.

Patterns of nucleotide polymorphism distinguish temperate and tropical wild isolates of Caenorhabditis briggsae

Caenorhabditis briggsae provides a natural comparison species for the model nematode C. elegans, given their similar morphology, life history, and hermaphroditic mode of reproduction. Despite C. briggsae boasting a published genome sequence and establishing Caenorhabditis as a model genus for genetics and development, little is known about genetic variation across the geographic range of this species. In this study, we greatly expand the collection of natural isolates and characterize patterns of nucleotide variation for six loci in 63 strains from three continents. The pattern of polymorphisms reveals differentiation between C. briggsae strains found in temperate localities in the northern hemisphere from those sampled near the Tropic of Cancer, with diversity within the tropical region comparable to what is found for C. elegans in Europe. As in C. elegans, linkage disequilibrium is pervasive, although recombination is evident among some variant sites, indicating that outcrossing has occurred at a low rate in the history of the sample. In contrast to C. elegans, temperate regions harbor extremely little variation, perhaps reflecting colonization and recent expansion of C. briggsae into northern latitudes. We discuss these findings in relation to their implications for selection, demographic history, and the persistence of self-fertilization.

Population genetics of the bovine/cattle lungworm (Dictyocaulus viviparus) based on mtDNA and AFLP marker techniques

Mitochondrial DNA (mtDNA) sequence data and amplified fragment length polymorphism (AFLP) patterns were compared for the lungworm Dictyocaulus viviparus, a nematode parasite of cattle. Eight individual D. viviparus samples from each of 8 herds in Sweden and 1 laboratory isolate were analysed, with the aim of describing the diversity and genetic structure in populations using different genetic markers on exactly the same DNA samples. There was qualitative agreement between the whole-genome AFLP data and the mtDNA sequence data, both indicating relatively strong genetic differentiation among the Swedish farms. However, the AFLP data detected much more genetic variation than did the mtDNA data, even after allowing for the different inheritance patterns of the markers, and indicated that there was much less differentiation among the populations. The mtDNA data therefore seemed to be more informative about the most recent history of the parasite populations, as the general patterns were less obscured by detailed inter-relationships among individual worms. The 4 mtDNA genes sequenced (1542 bp) showed consistent patterns, although there was more genetic variation in the protein-coding genes than in the structural RNA genes. Furthermore, there appeared to be at least 3 distinct genetic groups of D. viviparus infecting Swedish cattle, 1 of which was predominant and showed considerable differentiation between farms, but not necessarily within farms. Second, the 2 smaller genetic groups occurred on farms where the predominant group also occurred, suggesting that these farms have had multiple introductions of D. viviparus.

Nucleotide polymorphism and linkage disequilibrium in wild populations of the partial selfer Caenorhabditis elegans

An understanding of the relative contributions of different evolutionary forces on an organism's genome requires an accurate description of the patterns of genetic variation within and between natural populations. To this end, I report a survey of nucleotide polymorphism in six loci from 118 strains of the nematode Caenorhabditis elegans. These strains derive from wild populations of several regions within France, Germany, and new localities in Scotland, in addition to stock center isolates. Overall levels of silent-site diversity are low within and between populations of this self-fertile species, averaging 0.2% in European samples and 0.3% worldwide. Population structure is present despite a lack of association of sequences with geography, and migration appears to occur at all geographic scales. Linkage disequilibrium is extensive in the C. elegans genome, extending even between chromosomes. Nevertheless, recombination is clearly present in the pattern of polymorphisms, indicating that outcrossing is an infrequent, but important, feature in this species ancestry. The range of outcrossing rates consistent with the data is inferred from linkage disequilibrium, using "scattered" samples representing the collecting phase of the coalescent process in a subdivided population. I propose that genetic variation in this species is shaped largely by population subdivision due to self-fertilization coupled with long- and short-range migration between subpopulations.

Evolution of base-substitution gradients in primate mitochondrial genomes

Inferences of phylogenies and dates of divergence rely on accurate modeling of evolutionary processes; they may be confounded by variation in substitution rates among sites and changes in evolutionary processes over time. In vertebrate mitochondrial genomes, substitution rates are affected by a gradient along the genome of the time spent being single-stranded during replication, and different types of substitutions respond differently to this gradient. The gradient is controlled by biological factors including the rate of replication and functionality of repair mechanisms; little is known, however, about the consistency of the gradient over evolutionary time, or about how evolution of this gradient might affect phylogenetic analysis. Here, we evaluate the evolution of response to this gradient in complete primate mitochondrial genomes, focusing particularly on A-->G substitutions, which increase linearly with the gradient. We developed a methodology to evaluate the posterior probability densities of the response parameter space, and used likelihood ratio tests and mixture models with different numbers of classes to determine whether groups of genomes have evolved in a similar fashion. Substitution gradients usually evolve slowly in primates, but there have been at least two large evolutionary jumps: on the lineage leading to the great apes, and a convergent change on the lineage leading to baboons (Papio). There have also been possible convergences at deeper taxonomic levels, and different types of substitutions appear to evolve independently. The placements of the tarsier and the tree shrew within and in relation to primates may be incorrect because of convergence in these factors.

Host isolation and patterns of genetic variability in three populations of Teladorsagia from sheep

We have used a mitochondrial marker to explore the population genetics of an economically important parasite of sheep, Teladorsagia. We examined diversity within and between parasites from three very different host populations, as well as within and between individual hosts. One of our study populations, the Soay sheep on Hirta, St Kilda, is unusually isolated with no sheep having been introduced to the island since 1932. Worm haplotypes from Hirta were compared with those from two other host populations. Remarkably, despite its historical isolation the Hirta population shows similar levels of within-population diversity to the other study populations. No divergence between the three Teladorsagia populations was found, consistent with gene flow between the populations. The high diversity within Teladorsagia populations provides compelling evidence that this variability is a general feature of parasitic nematode populations. Such diversity may be caused by high effective population size, coupled with an increased mutation rate for mtDNA, which has important implications for the spread of anthelmintic resistance in nematode populations.

Macroevolutionary relationships of species of Drosophila melanogaster group based on mtDNA sequences

The phylogenetic relationships among the Drosophila melanogaster group species were analyzed using approximately 1700 nucleotide-long sequences of the mitochondrial DNA. Phylogenetic analysis was performed using this region consisting of a part of the cytochrome b (cytb) coding gene, the entire coding sequences of tRNA-Leu, tRNA-Ser and the first subunit of NADH dehydrogenase (NADH1), and a part of the 16S-rRNA gene. The study of these sequences showed that this region of mtDNA is very invariable, as regards with the type of the genes that it contains, as well as the order that they are located on it. The resulting phylogenetic trees reveal a topology that separates the species into three main ancestral lines, leading to the following subgroups: (a) ananassae subgroup, (b) montium subgroup, and (c) melanogaster and Oriental subgroups. The inferred topology complements and generally agrees with previously proposed classifications based on morphological and molecular data.

Rates of deleterious mutation and the evolution of sex in Caenorhabditis

A variety of models propose that the accumulation of deleterious mutations plays an important role in the evolution of breeding systems. These models make predictions regarding the relative rates of protein evolution and deleterious mutation in taxa with contrasting modes of reproduction. Here we compare available coding sequences from one obligately outcrossing and two primarily selfing species of Caenorhabditis to explore the potential for mutational models to explain the evolution of breeding system in this clade. If deleterious mutations interact synergistically, the mutational deterministic hypothesis predicts that a high genomic deleterious mutation rate (U) will offset the reproductive disadvantage of outcrossing relative to asexual or selfing reproduction. Therefore, C. elegans and C. briggsae (both largely selfing) should both exhibit lower rates of deleterious mutation than the obligately outcrossing relative C. remanei. Using a comparative approach, we estimate U to be equivalent (and < 1) among all three related species. Stochastic mutational models, Muller's ratchet and Hill-Robertson interference, are expected to cause reductions in the effective population size in species that rarely outcross, thereby allowing deleterious mutations to accumulate at an elevated rate. We find only limited support for more rapid molecular evolution in selfing lineages. Overall, our analyses indicate that the evolution of breeding system in this group is unlikely to be explained solely by available mutational models.

Molecular evolution and quantitative variation for chemosensory behaviour in the nematode genus Caenorhabditis

Caenorhabditis elegans is a model organism in biology, yet despite the tremendous information generated from genetic, genomic and functional analyses, C. elegans has rarely been used to address questions in ecological genetics. Here, we analyse genetic variation for chemosensory behaviour, an ecologically important trait that is also genetically well characterized, at both the phenotypic and molecular levels within three species of the genus Caenorhabditis. We show that the G-protein ODR-3 plays an important role in chemosensory avoidance behaviour and identify orthologues of odr-3 in C. briggsae and C. remanei. Both quantitative genetic analysis of chemosensory behaviour and molecular population genetic analysis of odr-3 show that there is little genetic variation among a worldwide collection of isolates of the primarily selfing C. elegans, whereas there is substantially more variation within a single population of the outcrossing C. remanei. Although there are a large number of substitutions at silent sites within odr-3 among the three species, molecular evolution at the protein level is extremely conserved, suggesting that odr-3 plays an important role in cell signalling during chemosensation and/or neuronal cilia development in C. remanei and in C. briggsae as it does in C. elegans. Our results suggest that C. remanei may be a more suitable subject for ecological and evolutionary genetic studies than C. elegans.

Sex-determination gene and pathway evolution in nematodes

The pathway that controls sexual fate in the nematode Caenorhabditis elegans has been well characterized at the molecular level. By identifying differences between the sex-determination mechanisms in C. elegans and other nematode species, it should be possible to understand how complex sex-determining pathways evolve. Towards this goal, orthologues of many of the C. elegans sex regulators have been isolated from other members of the genus Caenorhabditis. Rapid sequence evolution is observed in every case, but several of the orthologues appear to have conserved sex-determining roles. Thus extensive sequence divergence does not necessarily coincide with changes in pathway structure, although the same forces may contribute to both. This review summarizes recent findings and, with reference to results from other animals, offers explanations for why sex-determining genes and pathways appear to be evolving rapidly. Experimental strategies that hold promise for illuminating pathway differences between nematodes are also discussed.

Molecular identification of natural hybrids between Trichinella nativa and Trichinella T6 provides evidence of gene flow and ongoing genetic divergence

To date, there are no data available on the population genetics of Trichinella due to the lack of genetic markers and the difficulty of working with such small parasites. In the Arctic region of North America and along the Rocky Mountains, there exist two genotypes of Trichinella, Trichinella nativa and Trichinella T6, respectively, which are well differentiated by biochemical and molecular characters. However, both are resistant to freezing, show other common biological characters (e.g. low or no infectivity to rodents and swine) and produce fertile F1 offspring upon interbreeding. To data, these two genotypes have been considered allopatric. In this study, we detected both genotypes in wolves of the same wolf packs in Alaska, suggesting sympatry. A single GTT trinucleotide present in the ITS-2 sequence of T. nativa but not in Trichinella T6 was used as a genetic marker to study gene flow for this character in both a murine infection model and in larvae from naturally-infected Alaskan wolves. Only F1 larvae originating from a cross between T. nativa male and Trichinella T6 female were able to produce F2 offspring. Larvae (F1) originating from a cross between Trichinella T6 male and T. nativa female were not reproductively viable. As expected, all F1 larvae showed a heterozygote pattern for the GTT character upon heteroduplex analysis; however, within the F2 population, the number of observed heterozygotes (n=52) was substantially higher than expected (n=39.08), as supported by the F(is) index, and was not in the Hardy-Weinberg equilibrium. Larvae from two of the 16 Trichinella positive Alaskan wolves, showed the Trichinella T6 pattern or the T. nativa/Trichinella T6 hybrid pattern. Our data demonstrate that T. nativa and Trichinella T6 live in sympatry at least in Alaskan wolves, where T. nativa occurs more frequently (69%) than Trichinella T6 (31%). One explanation for this phenomenon is that glacial periods may have caused a geographical relocation, colonisation and independent evolution of T. nativa within the Rocky Mountains, resulting in a bifurcation of the freeze-resistant genotype. Additional studies will be required to test this hypothesis.

Levels of DNA polymorphism vary with mating system in the nematode genus caenorhabditis

Self-fertilizing species often harbor less genetic variation than cross-fertilizing species, and at least four different models have been proposed to explain this trend. To investigate further the relationship between mating system and genetic variation, levels of DNA sequence polymorphism were compared among three closely related species in the genus Caenorhabditis: two self-fertilizing species, Caenorhabditis elegans and C. briggsae, and one cross-fertilizing species, C. remanei. As expected, estimates of silent site nucleotide diversity were lower in the two self-fertilizing species. For the mitochondrial genome, diversity in the selfing species averaged 42% of diversity in C. remanei. Interestingly, the reduction in genetic variation was much greater for the nuclear than for the mitochondrial genome. For two nuclear genes, diversity in the selfing species averaged 6 and 13% of diversity in C. remanei. We argue that either population bottlenecks or the repeated action of natural selection, coupled with high levels of selfing, are likely to explain the observed reductions in species-wide genetic diversity.

Conservation of long-range synteny and microsynteny between the genomes of two distantly related nematodes

BACKGROUND

Comparisons between the genomes of the closely related nematodes Caenorhabditis elegans and Caenorhabditis briggsae reveal high rates of rearrangement, with a bias towards within-chromosome events. To assess whether this pattern is true of nematodes in general, we have used genome sequence to compare two nematode species that last shared a common ancestor approximately 300 million years ago: the model C. elegans and the filarial parasite Brugia malayi.

RESULTS

An 83 kb region flanking the gene for Bm-mif-1 (macrophage migration inhibitory factor, a B. malayi homolog of a human cytokine) was sequenced. When compared to the complete genome of C. elegans, evidence for conservation of long-range synteny and microsynteny was found. Potential C. elegans orthologs for II of the 12 protein-coding genes predicted in the B. malayi sequence were identified. Ten of these orthologs were located on chromosome I, with eight clustered in a 2.3 Mb region. While several, relatively local, intrachromosomal rearrangements have occurred, the order, composition, and configuration of two gene clusters, each containing three genes, was conserved. Comparison of B. malayi BAC-end genome survey sequence to C. elegans also revealed a bias towards intrachromosome rearrangements.

CONCLUSIONS

We suggest that intrachromosomal rearrangement is a major force driving chromosomal organization in nematodes, but is constrained by the interdigitation of functional elements of neighboring genes.

Species-specific identification of human hookworms by PCR of the mitochondrial cytochrome oxidase I gene

Significant differences in the life histories of the human hookworms Ancylostoma duodenale and Necator americanus necessitate their differentiation for epidemiological studies and the design of control programs. Current methods of identification require time-consuming, labor-intensive techniques. A polymerase chain reaction (PCR)-based method that enables rapid species identification is described. The mitochondrial cytochrome oxidase I genes of both species were sequenced, and species-specific primer sets were designed. The primers were used in PCR to amplify 585-bp fragments of the cytochrome oxidase gene from individual hookworm eggs, larvae, and adults. The technique was also able to identify mixed infections containing equal amounts of eggs from each species. The technique is rapid, technically simple, and sensitive and will permit the accurate identification of human hookworms in epidemiological field studies.

High direct estimate of the mutation rate in the mitochondrial genome of Caenorhabditis elegans

Mutations in the mitochondrial genome have been implicated in numerous human genetic disorders and offer important data for phylogenetic, forensic, and population genetic studies. Using a long-term series of Caenorhabditis elegans mutation accumulation lines, we performed a wide-scale screen for mutations in the mitochondrial genome that revealed a mutation rate that is two orders of magnitude higher than previous indirect estimates, a highly biased mutational spectrum, multiple mutations affecting coding function, as well as mutational hotspots at homopolymeric nucleotide stretches.

Regulatory elements required for development of caenorhabditis elegans hermaphrodites are conserved in the tra-2 homologue of C. remanei, a male/female sister species

The Caenorhabditis elegans hermaphrodite is essentially a female that produces sperm. In C. elegans, tra-2 promotes female fates and must be repressed to achieve hermaphrodite spermatogenesis. In an effort to learn how mating systems evolve, we have cloned tra-2 from C. remanei, the closest gonochoristic relative of C. elegans. We found its structure to be similar to that of Ce-tra-2 but its sequence to be divergent. RNA interference demonstrates that Cr-tra-2 promotes female fates. Two sites of tra-2 regulation are required for the onset of hermaphrodite spermatogenesis in C. elegans. One, the MX region of TRA-2, is as well conserved in C. remanei as it is in C. briggsae (another male/hermaphrodite species), suggesting that this control is not unique to hermaphrodites. Another, the DRE/TGE element of the tra-2 3' UTR, was not detected by sequence analysis. However, gel-shift assays demonstrate that a factor in C. remanei can bind specifically to the Cr-tra-2 3' UTR, suggesting that this translational control is also conserved. We propose that both controls are general and do not constitute a novel "switch" that enables sexual mosaicism in hermaphrodites. However, subtle quantitative or qualitative differences in their employment may underlie differences in mating system seen in Caenorhabditis.

Secretion of a novel class of iFABPs in nematodes: coordinate use of the Ascaris/Caenorhabditis model systems

A novel fatty acid binding protein, As-p18, is secreted into both the perivitelline and perienteric fluids of the parasitic nematode, Ascaris suum, and at least eight potential homologues of As-p18 have been identified in the Caenorhabditis elegans genome. The products of the three most closely related homologues are fatty acid binding proteins (LBP-1, LBP-2 and LBP-3) which contain putative secretory signals. Phylogenetic analysis revealed that these secreted fatty acid binding proteins comprise a distinct gene class within the fatty acid binding protein family and are possibly unique to nematodes. To examine the potential sites of As-p18 secretion, the expression of the putative promoters of the C. elegans homologues was examined with GFP reporter constructs. The developmental expression of lbp-1 was identical to that of As-p18 and consistent with the secretion of LBP-1 from the hypodermis to the perivitelline fluid. The expression patterns of lbp-2 and lbp-3 were consistent with the secretion of LBP-2 and LBP-3 from muscle into the perienteric fluid later in development. These studies demonstrate that at least some perivitelline fluid proteins appear to be secreted from the hypodermis prior to the formation of the cuticle and, perhaps more importantly, that this coordinate C. elegans/A. suum approach may be potentially useful for examining a number of key physiological processes in parasitic nematodes.

Onchocerca volvulus: limited heterogeneity in the nuclear and mitochondrial genomes

West African populations of Onchocerca volvulus endemic to the rain forest and savanna bioclimes of West Africa differ in their ability to induce ocular disease in infected individuals. In recent years, both clinical- and animal-model-based studies have implicated particular parasite antigens in the development of ocular onchocerciasis. To test the hypothesis that the difference in pathogenic potential of blinding and nonblinding parasites might be reflected in qualitative differences in antigens that have been implicated in the development of ocular onchocerciasis, we compared the sequences of two parasite antigens implicated in the development of ocular disease in blinding- and nonblinding-strain parasites. The results demonstrated a high level of homogeneity between the parasite strains in these genes. The study was extended to include additional nuclear genes encoding antigens that are commonly recognized by individuals infected with O. volvulus and to the mitochondrial genome of the parasite. The results demonstrate a high degree of homogeneity in both the nuclear and the mitochondrial genomes among O. volvulus isolates collected from several different sites in Africa and in the Americas. This high degree of genetic homogeneity may reflect the passage of the parasite through a recent genetic bottleneck.

Mitochondrial DNA diversity in nematodes

The relatively small literature on mitochondrial DNA (mtDNA) diversity in nematode species is summarized here. Nematodes show a wide range of overall genetic diversities and population genetic structures. Species-wide levels of diversity correlate strongly with the breeding system and other life cycle features that control effective population size. Obligate outcrossers that parasitize mobile vertebrate hosts are the most diverse, species having hermaphroditic stages are less so, and species having asexual reproductive stages appear even less diverse. Nevertheless, these conclusions are preliminary because there exist so few data on DNA diversity in nematodes. What is needed are more comparative studies using similar sampling designs and the same DNA markers, including nuclear loci and further work with mtDNA.

An evolutionary framework for the study of developmental evolution in a set of nematodes related to Caenorhabditis elegans

Nematodes are known to be a useful system for studies of comparative development. Here we perform a molecular phylogenetic analysis to allow for the independent interpretation of the developmental and morphological changes observed among a selected set of nematode species. Our molecular phylogenetic analysis is based on coding regions of the genes for RNA polymerase II, the small subunit rRNA and an expansion segment of the large subunit rRNA. Sequences were compared from five species in the family (Rhabditidae) that includes the developmental model organism Caenorhabditis elegans and from an outgroup taxon Aduncospiculum halicti (Diplogasterina). The phylogenetic analysis does not support the monophyly of the subfamily Mesorhabditinae and identifies the unnamed strain PS1010 as a sister taxon of C. elegans despite its morphologically divergent buccal capsule. On the basis of the inferred framework, we can begin to interpret the evolution of vulval development and of morphological differences among these nematode species.

Quantitative trait loci controlling halothane sensitivity in Caenorhabditis elegans

Genetic analysis is an essential tool for defining the molecular mechanisms whereby volatile anesthetics (VA) disrupt nervous system function. However, the degree of natural variation of the genetic determinants of VA sensitivity has not been determined nor have mutagenesis approaches been very successful at isolating significantly resistant mutant strains. Thus, a quantitative genetic approach was taken toward these goals. Recombinant-inbred strains derived from two evolutionarily distinct lineages of the nematode Caenorhabditis elegans were tested for sensitivity to clinically relevant concentrations (0.3-0.5 mM) of the VA halothane. The halothane sensitivities of coordinated movement and male mating behavior were highly variant among the recombinant-inbred strains with a range of EC50 values of 13- and 4-fold, respectively. Both traits were highly heritable (H2 = 0.82, 0.87, respectively). Several strains were found to be significantly resistant to halothane when compared with the wild-type strain N2. A major locus or loci mapping to the middle of chromosome V accounted for more than 40% of the phenotypic variance for both traits. Five weaker loci, four of which interact, explained most of the remaining variance. None of the halothane-sensitivity quantitative trait loci significantly affected behavior in the absence of halothane or halothane's potency for C. elegans immobilization, which requires 5-fold higher drug concentrations. Thus, the quantitative trait loci are unlikely to result from differences in halothane-independent (native) behavior or differences in halothane metabolism or permeability. Rather, these loci may code for targets and/or downstream effectors of halothane in the C. elegans nervous system or for modifiers of such gene products.

Natural variation and copulatory plug formation in Caenorhabditis elegans

Most of the available natural isolates of the nematode Caenorhabditis elegans have been examined and compared with the standard laboratory wild type (Bristol N2). Molecular markers, in particular transposon restriction fragment length polymorphisms, were used to assign these isolates to 22 different races, for which brood size and spontaneous male frequency were determined. Several distinctive traits were observed in some of these races. One example is mab-23, in a race from Vancouver, which leads to severe distortion of male genitalia and prevents male mating. Another is gro-1, segregating in a Californian race, which is associated with slow growth, heat resistance and longevity. Many races differ from N2 in carrying a dominant allele at the plg-1 locus, causing copulatory plug formation by males. Properties and possible advantages of the plugging trait have been investigated. The dominant plg-1 allele does not lead to increased male mating efficiency, but males from a Stanford race (CB4855), in which the plugging trait was first observed, are much more virile than N2 males. Crosses between N2 and CB4855 indicate that the higher virility is due to multiple factors. Size differences between N2 and CB4855 are associated with factors mapping to LGV and LGX.

Evolution of sex determination in caenorhabditis: unusually high divergence of tra-1 and its functional consequences

The tra-1 gene is a terminal regulator of somatic sex in Caenorhabditis elegans: high tra-1 activity elicits female development, low tra-1 activity elicits male development. To investigate the function and evolution of tra-1, we examined the tra-1 gene from the closely related nematode C. briggsae. Ce-tra-1 and Cb-tra-1 are unusually divergent. Each gene generates two transcripts, but only one of these is present in both species. This common transcript encodes TRA-1A, which shows only 44% amino acid identity between the species, a figure much lower than that for previously compared genes. A Cb-tra-1 transgene rescues many tissues of tra-1(null) mutants of C. elegans but not the somatic gonad or germ line. This transgene also causes nongonadal feminization of XO animals, indicating incorrect sexual regulation. Alignment of Ce-TRA-1A and Cb-TRA-1A defines several conserved regions likely to be important for tra-1 function. The phenotypic differences between Ce-tra-1(null) mutants rescued by Cb-tra-1 transgenes and wild-type C. elegans indicate significant divergence of regulatory regions. These molecular and functional studies suggest that evolution of sex determination in nematodes is rapid and genetically complex.

Caenorhabditis globin genes: rapid intronic divergence contrasts with conservation of silent exonic sites

Globin genes from the Caenorhabditis species briggsae and remanei were identified and compared with a previously described C. elegans globin gene. The encoded globins share between 86% and 93% amino acid identity, with most of the changes in or just before the putative B helix. C. remanei was found to have two globin alleles, Crg1-1 and Crg1-2. The coding sequence for each is interrupted by a single intron in the same position. The exons of the two genes are only 1% divergent at the nucleotide level and encode identical polypeptides. In contrast, intron sequence divergence is 16% and numerous insertions and deletions have significantly altered the size and content of both introns. Genetic crosses show that Crg1-1 and Crg1-2 segregate as alleles. Homozygous lines for each allele were constructed and northern analysis confirmed the expression of both alleles. These data reveal an unusual situation wherein two alleles encoding identical proteins have diverged much more rapidly in their introns than the silent sites of their coding sequences, suggesting multiple gene conversion events.

On transition bias in mitochondrial genes of pocket gophers

The relative contribution of mutation and purifying selection to transition bias has not been quantitatively assessed in mitochondrial protein genes. The observed transition/transversion (s/v) ratio is (micros Ps)/(microv Pv), where micros and microv denote mutation rate of transitions and transversions, respectively, and Ps and Pv denote fixation probabilities of transitions and transversions, respectively. Because selection against synonymous transitions can be assumed to be roughly equal to that against synonymous transversions, Ps/Pv approximately 1 at fourfold degenerate sites, so that the s/v ratio at fourfold degenerate sites is approximately micros/microv, which is a measure of mutational contribution to transition bias. Similarly, the s/v ratio at nondegenerate sites is also an estimate of micros/microv if we assume that selection against nonsynonymous transitions is roughly equal to that against nonsynonymous transversions. In two mitochondrial genes, cytochrome oxidase subunit I (COI) and cytochrome b (cyt-b) in pocket gophers, the s/v ratio is about two at nondegenerate and fourfold degenerate sites for both the COI and the cyt-b genes. This implies that mutation contribution to transition bias is relatively small. In contrast, the s/v ratio is much greater at twofold degenerate sites, being 48 for COI and 40 for cyt-b. Given that the micros/microv ratio is about 2, the Ps/Pv ratio at twofold degenerate sites must be on the order of 20 or greater. This suggests a great effect of purifying selection on transition bias in mitochondrial protein genes because transitions are synonymous and transversions are nonsynonymous at twofold degenerate sites in mammalian mitochondrial genes. We also found that nonsynonymous mutations at twofold degenerate sites are more neutral than nonsynonymous mutations at nondegenerate sites, and that the COI gene is subject to stronger purifying selection than is the cyt-b gene. A model is presented to integrate the effect of purifying selection, codon bias, DNA repair and GC content on s/v ratio of protein-coding genes.

Molecular genetic approaches to parasite identification: their value in diagnostic parasitology and systematics

A wide range of approaches is available to parasitologists to aid in specific parasite identification and to formulate phylogenetic relationships. This review emphasises the usefulness of molecular genetic techniques, especially DNA-based procedures, in addressing problems of identification, characterisation and phylogeny of parasites. It should be stressed that an understanding of the various DNA approaches, techniques and target genes most likely to be effective in addressing key issues in diagnostic parasitology and systematics is still developing. Nevertheless, DNA methods clearly have great potential with regard to specificity and sensitivity, and applications will increase further with technological advance. Indeed, because of the minimal requirements for material, PCR-based methods especially should prove of immense value in future studies with parasites.