Sequencing and annotation of mitochondrial genomes from individual parasitic helminths

Mitochondrial (mt) genomics has significant implications in a range of fundamental areas of parasitology, including evolution, systematics, and population genetics as well as explorations of mt biochemistry, physiology, and function. Mt genomes also provide a rich source of markers to aid molecular epidemiological and ecological studies of key parasites. However, there is still a paucity of information on mt genomes for many metazoan organisms, particularly parasitic helminths, which has often related to challenges linked to sequencing from tiny amounts of material. The advent of next-generation sequencing (NGS) technologies has paved the way for low cost, high-throughput mt genomic research, but there have been obstacles, particularly in relation to post-sequencing assembly and analyses of large datasets. In this chapter, we describe protocols for the efficient amplification and sequencing of mt genomes from small portions of individual helminths, and highlight the utility of NGS platforms to expedite mt genomics. In addition, we recommend approaches for manual or semi-automated bioinformatic annotation and analyses to overcome the bioinformatic "bottleneck" to research in this area. Taken together, these approaches have demonstrated applicability to a range of parasites and provide prospects for using complete mt genomic sequence datasets for large-scale molecular systematic and epidemiological studies. In addition, these methods have broader utility and might be readily adapted to a range of other medium-sized molecular regions (i.e., 10-100 kb), including large genomic operons, and other organellar (e.g., plastid) and viral genomes.

A functional genomic screen for evolutionarily conserved genes required for lifespan and immunity in germline-deficient C. elegans

The reproductive system regulates lifespan in insects, nematodes and vertebrates. In Caenorhabditis elegans removal of germline increases lifespan by 60% which is dependent upon insulin signaling, nuclear hormone signaling, autophagy and fat metabolism and their microRNA-regulators. Germline-deficient C. elegans are also more resistant to various bacterial pathogens but the underlying molecular mechanisms are largely unknown. Firstly, we demonstrate that previously identified genes that regulate the extended lifespan of germline-deficient C. elegans (daf-2, daf-16, daf-12, tcer-1, mir-7.1 and nhr-80) are also essential for resistance to the pathogenic bacterium Xenorhabdus nematophila. We then use a novel unbiased approach combining laser cell ablation, whole genome microarrays, RNAi screening and exposure to X. nematophila to generate a comprehensive genome-wide catalog of genes potentially required for increased lifespan and innate immunity in germline-deficient C. elegans. We find 3,440 genes to be upregulated in C. elegans germline-deficient animals in a gonad dependent manner, which are significantly enriched for genes involved in insulin signaling, fatty acid desaturation, translation elongation and proteasome complex function. Using RNAi against a subset of 150 candidate genes selected from the microarray results, we show that the upregulated genes such as transcription factor DAF-16/FOXO, the PTEN homolog lipid phosphatase DAF-18 and several components of the proteasome complex (rpn-6.1, rpn-7, rpn-9, rpn-10, rpt-6, pbs-3 and pbs-6) are essential for both lifespan and immunity of germline deficient animals. We also identify a novel role for genes including par-5 and T12G3.6 in both lifespan-extension and increased survival on X. nematophila. From an evolutionary perspective, most of the genes differentially expressed in germline deficient C. elegans also show a conserved expression pattern in germline deficient Pristionchus pacificus, a nematode species that diverged from C. elegans 250-400 MYA.

Operons are a conserved feature of nematode genomes

The organization of genes into operons, clusters of genes that are co-transcribed to produce polycistronic pre-mRNAs, is a trait found in a wide range of eukaryotic groups, including multiple animal phyla. Operons are present in the class Chromadorea, one of the two main nematode classes, but their distribution in the other class, the Enoplea, is not known. We have surveyed the genomes of Trichinella spiralis, Trichuris muris, and Romanomermis culicivorax and identified the first putative operons in members of the Enoplea. Consistent with the mechanism of polycistronic RNA resolution in other nematodes, the mRNAs produced by genes downstream of the first gene in the T. spiralis and T. muris operons are trans-spliced to spliced leader RNAs, and we are able to detect polycistronic RNAs derived from these operons. Importantly, a putative intercistronic region from one of these potential enoplean operons confers polycistronic processing activity when expressed as part of a chimeric operon in Caenorhabditis elegans. We find that T. spiralis genes located in operons have an increased likelihood of having operonic C. elegans homologs. However, operon structure in terms of synteny and gene content is not tightly conserved between the two taxa, consistent with models of operon evolution. We have nevertheless identified putative operons conserved between Enoplea and Chromadorea. Our data suggest that operons and "spliced leader" (SL) trans-splicing predate the radiation of the nematode phylum, an inference which is supported by the phylogenetic profile of proteins known to be involved in nematode SL trans-splicing.

Polymerase theta-mediated end joining of replication-associated DNA breaks in C. elegans

DNA lesions that block replication fork progression are drivers of cancer-associated genome alterations, but the error-prone DNA repair mechanisms acting on collapsed replication are incompletely understood, and their contribution to genome evolution largely unexplored. Here, through whole-genome sequencing of animal populations that were clonally propagated for more than 50 generations, we identify a distinct class of deletions that spontaneously accumulate in C. elegans strains lacking translesion synthesis (TLS) polymerases. Emerging DNA double-strand breaks are repaired via an error-prone mechanism in which the outermost nucleotide of one end serves to prime DNA synthesis on the other end. This pathway critically depends on the A-family polymerase theta, which protects the genome against gross chromosomal rearrangements. By comparing the genomes of isolates of C. elegans from different geographical regions, we found that in fact most spontaneously evolving structural variations match the signature of polymerase theta-mediated end joining (TMEJ), illustrating that this pathway is an important source of genetic diversification.

A novel flavivirus in the soybean cyst nematode

Heterodera glycines, the soybean cyst nematode (SCN), is a subterranean root pathogen that causes the most damaging disease of soybean in the USA. A novel nematode virus genome, soybean cyst nematode virus 5 (SbCNV-5), was identified in RNA sequencing data from SCN eggs and second-stage juveniles. The SbCNV-5 RNA-dependent RNA polymerase and RNA helicase domains had homology to pestiviruses in the family Flaviviridae, suggesting that SbCNV-5 is a positive-polarity ssRNA virus. SbCNV-5 RNA was present in all nematode developmental stages, indicating a transovarial mode of transmission, but is also potentially sexually transmitted via the male. SbCNV-5 was common in SCN laboratory cultures and in nematode populations isolated from the field. Transmission electron microscopy of sections from a female SCN showed virus particles budding from the endoplasmic reticulum and in endosomes. The size of the viral genome was 19 191 nt, which makes it much larger than other known pestiviruses. Additionally, the presence of a methyltransferase in the SbCNV-5 genome is atypical for a pestivirus. When cDNA sequences were mapped to the genome of SbCNV-5, a disproportionate number aligned to the 3' NTR, suggesting that SbCNV-5 produces a subgenomic RNA, which was confirmed by RNA blot analysis. As subgenomic RNAs and methyltransferases do not occur in pestiviruses, we conclude that SbCNV-5 is a new flavivirus infecting SCNs.

Diversity and expression of microRNAs in the filarial parasite, Brugia malayi

Human filarial parasites infect an estimated 120 million people in 80 countries worldwide causing blindness and the gross disfigurement of limbs and genitals. An understanding of RNA-mediated regulatory pathways in these parasites may open new avenues for treatment. Toward this goal, small RNAs from Brugia malayi adult females, males and microfilariae were cloned for deep-sequencing. From ∼30 million sequencing reads, 145 miRNAs were identified in the B. malayi genome. Some microRNAs were validated using the p19 RNA binding protein and qPCR. B. malayi miRNAs segregate into 99 families each defined by a unique seed sequence. Sixty-one of the miRNA families are highly conserved with homologues in arthropods, vertebrates and helminths. Of those miRNAs not highly conserved, homologues of 20 B. malayi miRNA families were found in vertebrates. Nine B. malayi miRNA families appear to be filarial-specific as orthologues were not found in other organisms. The miR-2 family is the largest in B. malayi with 11 members. Analysis of the sequences shows that six members result from a recent expansion of the family. Library comparisons found that 1/3 of the B. malayi miRNAs are differentially expressed. For example, miR-71 is 5–7X more highly expressed in microfilariae than adults. Studies suggest that in C.elegans, miR-71 may enhance longevity by targeting the DAF-2 pathway. Characterization of B. malayi miRNAs and their targets will enhance our understanding of their regulatory pathways in filariads and aid in the search for novel therapeutics.

The genome and life-stage specific transcriptomes of Globodera pallida elucidate key aspects of plant parasitism by a cyst nematode

BACKGROUND

Globodera pallida is a devastating pathogen of potato crops, making it one of the most economically important plant parasitic nematodes. It is also an important model for the biology of cyst nematodes. Cyst nematodes and root-knot nematodes are the two most important plant parasitic nematode groups and together represent a global threat to food security.

RESULTS

We present the complete genome sequence of G. pallida, together with transcriptomic data from most of the nematode life cycle, particularly focusing on the life cycle stages involved in root invasion and establishment of the biotrophic feeding site. Despite the relatively close phylogenetic relationship with root-knot nematodes, we describe a very different gene family content between the two groups and in particular extensive differences in the repertoire of effectors, including an enormous expansion of the SPRY domain protein family in G. pallida, which includes the SPRYSEC family of effectors. This highlights the distinct biology of cyst nematodes compared to the root-knot nematodes that were, until now, the only sedentary plant parasitic nematodes for which genome information was available. We also present in-depth descriptions of the repertoires of other genes likely to be important in understanding the unique biology of cyst nematodes and of potential drug targets and other targets for their control.

CONCLUSIONS

The data and analyses we present will be central in exploiting post-genomic approaches in the development of much-needed novel strategies for the control of G. pallida and related pathogens.

Genome of the human hookworm Necator americanus

The hookworm Necator americanus is the predominant soil-transmitted human parasite. Adult worms feed on blood in the small intestine, causing iron-deficiency anemia, malnutrition, growth and development stunting in children, and severe morbidity and mortality during pregnancy in women. We report sequencing and assembly of the N. americanus genome (244 Mb, 19,151 genes). Characterization of this first hookworm genome sequence identified genes orchestrating the hookworm's invasion of the human host, genes involved in blood feeding and development, and genes encoding proteins that represent new potential drug targets against hookworms. N. americanus has undergone a considerable and unique expansion of immunomodulator proteins, some of which we highlight as potential treatments against inflammatory diseases. We also used a protein microarray to demonstrate a postgenomic application of the hookworm genome sequence. This genome provides an invaluable resource to boost ongoing efforts toward fundamental and applied postgenomic research, including the development of new methods to control hookworm and human immunological diseases.

The genome of Romanomermis culicivorax: revealing fundamental changes in the core developmental genetic toolkit in Nematoda

BACKGROUND

The genetics of development in the nematode Caenorhabditis elegans has been described in exquisite detail. The phylum Nematoda has two classes: Chromadorea (which includes C. elegans) and the Enoplea. While the development of many chromadorean species resembles closely that of C. elegans, enoplean nematodes show markedly different patterns of early cell division and cell fate assignment. Embryogenesis of the enoplean Romanomermis culicivorax has been studied in detail, but the genetic circuitry underpinning development in this species has not been explored.

RESULTS

We generated a draft genome for R. culicivorax and compared its gene content with that of C. elegans, a second enoplean, the vertebrate parasite Trichinella spiralis, and a representative arthropod, Tribolium castaneum. This comparison revealed that R. culicivorax has retained components of the conserved ecdysozoan developmental gene toolkit lost in C. elegans. T. spiralis has independently lost even more of this toolkit than has C. elegans. However, the C. elegans toolkit is not simply depauperate, as many novel genes essential for embryogenesis in C. elegans are not found in, or have only extremely divergent homologues in R. culicivorax and T. spiralis. Our data imply fundamental differences in the genetic programmes not only for early cell specification but also others such as vulva formation and sex determination.

CONCLUSIONS

Despite the apparent morphological conservatism, major differences in the molecular logic of development have evolved within the phylum Nematoda. R. culicivorax serves as a tractable system to contrast C. elegans and understand how divergent genomic and thus regulatory backgrounds nevertheless generate a conserved phenotype. The R. culicivorax draft genome will promote use of this species as a research model.

A rapid protocol for integrating extrachromosomal arrays with high transmission rate into the C. elegans genome

Microinjecting DNA into the cytoplasm of the syncytial gonad of Caenorhabditis elegans is the main technique used to establish transgenic lines that exhibit partial and variable transmission rates of extrachromosomal arrays to the next generation. In addition, transgenic animals are mosaic and express the transgene in a variable number of cells. Extrachromosomal arrays can be integrated into the C. elegans genome using UV irradiation to establish nonmosaic transgenic strains with 100% transmission rate of the transgene. To that extent, F1 progenies of UV irradiated transgenic animals are screened for animals carrying a heterozygous integration of the transgene, which leads to a 75% Mendelian transmission rate to the F2 progeny. One of the challenges of this method is to distinguish between the percentage of transgene transmission in a population before (X% transgenic animals) and after integration (≥75% transgenic F2 animals). Thus, this method requires choosing a nonintegrated transgenic line with a percentage of transgenic animals that is significantly lower than the Mendelian segregation of 75%. Consequently, nonintegrated transgenic lines with an extrachromosomal array transmission rate to the next generation ≤60% are usually preferred for integration, and transgene integration in highly transmitting strains is difficult. Here we show that the efficiency of extrachromosomal arrays integration into the genome is increased when using highly transmitting transgenic lines (≥80%). The described protocol allows for easy selection of several independent lines with homozygous transgene integration into the genome after UV irradiation of transgenic worms exhibiting a high rate of extrachromosomal array transmission. Furthermore, this method is quite fast and low material consuming. The possibility of rapidly generating different lines that express a particular integrated transgene is of great interest for studies focusing on gene expression pattern and regulation, protein localization, and overexpression, as well as for the development of subcellular markers.

Extensively duplicated and transcriptionally active recent lateral gene transfer from a bacterial Wolbachia endosymbiont to its host filarial nematode Brugia malayi

BACKGROUND

Lymphatic filariasis is a neglected tropical disease afflicting more than 120 million people, while another 1.3 billion people are at risk of infection. The nematode worm Brugia malayi is one of the causative agents of the disease and exists in a mutualistic symbiosis with Wolbachia bacteria. Since extensive lateral gene transfer occurs frequently between Wolbachia and its hosts, we sought to measure the extent of such LGT in B. malayi by whole genome sequencing of Wolbachia-depleted worms.

RESULTS

A considerable fraction (at least 115.4-kbp, or 10.6%) of the 1.08-Mbp Wolbachia wBm genome has been transferred to its nematode host and retains high levels of similarity, including 227 wBm genes and gene fragments. Complete open reading frames were transferred for 32 of these genes, meaning they have the potential to produce functional proteins. Moreover, four transfers have evidence of life stage-specific regulation of transcription at levels similar to other nematode transcripts, strengthening the possibility that they are functional.

CONCLUSIONS

There is extensive and ongoing transfer of Wolbachia DNA to the worm genome and some transfers are transcribed in a stage-specific manner at biologically relevant levels.

The genome and transcriptome of Haemonchus contortus, a key model parasite for drug and vaccine discovery

BACKGROUND

The small ruminant parasite Haemonchus contortus is the most widely used parasitic nematode in drug discovery, vaccine development and anthelmintic resistance research. Its remarkable propensity to develop resistance threatens the viability of the sheep industry in many regions of the world and provides a cautionary example of the effect of mass drug administration to control parasitic nematodes. Its phylogenetic position makes it particularly well placed for comparison with the free-living nematode Caenorhabditis elegans and the most economically important parasites of livestock and humans.

RESULTS

Here we report the detailed analysis of a draft genome assembly and extensive transcriptomic dataset for H. contortus. This represents the first genome to be published for a strongylid nematode and the most extensive transcriptomic dataset for any parasitic nematode reported to date. We show a general pattern of conservation of genome structure and gene content between H. contortus and C. elegans, but also a dramatic expansion of important parasite gene families. We identify genes involved in parasite-specific pathways such as blood feeding, neurological function, and drug metabolism. In particular, we describe complete gene repertoires for known drug target families, providing the most comprehensive understanding yet of the action of several important anthelmintics. Also, we identify a set of genes enriched in the parasitic stages of the lifecycle and the parasite gut that provide a rich source of vaccine and drug target candidates.

CONCLUSIONS

The H. contortus genome and transcriptome provide an essential platform for postgenomic research in this and other important strongylid parasites.

TIMPs of parasitic helminths - a large-scale analysis of high-throughput sequence datasets

BACKGROUND

Tissue inhibitors of metalloproteases (TIMPs) are a multifunctional family of proteins that orchestrate extracellular matrix turnover, tissue remodelling and other cellular processes. In parasitic helminths, such as hookworms, TIMPs have been proposed to play key roles in the host-parasite interplay, including invasion of and establishment in the vertebrate animal hosts. Currently, knowledge of helminth TIMPs is limited to a small number of studies on canine hookworms, whereas no information is available on the occurrence of TIMPs in other parasitic helminths causing neglected diseases.

METHODS

In the present study, we conducted a large-scale investigation of TIMP proteins of a range of neglected human parasites including the hookworm Necator americanus, the roundworm Ascaris suum, the liver flukes Clonorchis sinensis and Opisthorchis viverrini, as well as the schistosome blood flukes. This entailed mining available transcriptomic and/or genomic sequence datasets for the presence of homologues of known TIMPs, predicting secondary structures of defined protein sequences, systematic phylogenetic analyses and assessment of differential expression of genes encoding putative TIMPs in the developmental stages of A. suum, N. americanus and Schistosoma haematobium which infect the mammalian hosts.

RESULTS

A total of 15 protein sequences with high homology to known eukaryotic TIMPs were predicted from the complement of sequence data available for parasitic helminths and subjected to in-depth bioinformatic analyses.

CONCLUSIONS

Supported by the availability of gene manipulation technologies such as RNA interference and/or transgenesis, this work provides a basis for future functional explorations of helminth TIMPs and, in particular, of their role/s in fundamental biological pathways linked to long-term establishment in the vertebrate hosts, with a view towards the development of novel approaches for the control of neglected helminthiases.

A satellite explosion in the genome of holocentric nematodes

Centromere sequences in the genome are associated with the formation of kinetochores, where spindle microtubules grow in mitosis. Centromere sequences usually have long tandem repeats (satellites). In holocentric nematodes it is not clear how kinetochores are formed during mitosis; they are distributed throughout the chromosomes. For this reason it appeared of interest to study the satellites in nematodes in order to determine if they offer any clue on how kinetochores are assembled in these species. We have studied the satellites in the genome of six nematode species. We found that the presence of satellites depends on whether the nematode chromosomes are holocentric or monocentric. It turns out that holocentric nematodes are unique because they have a large number of satellites scattered throughout their genome. Their number, length and composition are different in each species: they apparently have very little evolutionary conservation. In contrast, no scattered satellites are found in the monocentric nematode Trichinella spiralis. It appears that the absence/presence of scattered satellites in the genome distinguishes monocentric from holocentric nematodes. We conclude that the presence of satellites is related to the holocentric nature of the chromosomes of most nematodes. Satellites may stabilize a higher order structure of chromatin and facilitate the formation of kinetochores. We also present a new program, SATFIND, which is suited to find satellite sequences.

The draft genome and transcriptome of Panagrellus redivivus are shaped by the harsh demands of a free-living lifestyle

Nematodes compose an abundant and diverse invertebrate phylum with members inhabiting nearly every ecological niche. Panagrellus redivivus (the "microworm") is a free-living nematode frequently used to understand the evolution of developmental and behavioral processes given its phylogenetic distance to Caenorhabditis elegans. Here we report the de novo sequencing of the genome, transcriptome, and small RNAs of P. redivivus. Using a combination of automated gene finders and RNA-seq data, we predict 24,249 genes and 32,676 transcripts. Small RNA analysis revealed 248 microRNA (miRNA) hairpins, of which 63 had orthologs in other species. Fourteen miRNA clusters containing 42 miRNA precursors were found. The RNA interference, dauer development, and programmed cell death pathways are largely conserved. Analysis of protein family domain abundance revealed that P. redivivus has experienced a striking expansion of BTB domain-containing proteins and an unprecedented expansion of the cullin scaffold family of proteins involved in multi-subunit ubiquitin ligases, suggesting proteolytic plasticity and/or tighter regulation of protein turnover. The eukaryotic release factor protein family has also been dramatically expanded and suggests an ongoing evolutionary arms race with viruses and transposons. The P. redivivus genome provides a resource to advance our understanding of nematode evolution and biology and to further elucidate the genomic architecture leading to free-living lineages, taking advantage of the many fascinating features of this worm revealed by comparative studies.

First insights into the genetic diversity of the pinewood nematode in its native area using new polymorphic microsatellite loci

The pinewood nematode, Bursaphelenchus xylophilus, native to North America, is the causative agent of pine wilt disease and among the most important invasive forest pests in the East-Asian countries, such as Japan and China. Since 1999, it has been found in Europe in the Iberian Peninsula, where it also causes significant damage. In a previous study, 94 pairs of microsatellite primers have been identified in silico in the pinewood nematode genome. In the present study, specific PCR amplifications and polymorphism tests to validate these loci were performed and 17 microsatellite loci that were suitable for routine analysis of B. xylophilus genetic diversity were selected. The polymorphism of these markers was evaluated on nematodes from four field origins and one laboratory collection strain, all originate from the native area. The number of alleles and the expected heterozygosity varied between 2 and 11 and between 0.039 and 0.777, respectively. First insights into the population genetic structure of B. xylophilus were obtained using clustering and multivariate methods on the genotypes obtained from the field samples. The results showed that the pinewood nematode genetic diversity is spatially structured at the scale of the pine tree and probably at larger scales. The role of dispersal by the insect vector versus human activities in shaping this structure is discussed.

RNAi pathways contribute to developmental history-dependent phenotypic plasticity in C. elegans

Early environmental experiences profoundly influence adult phenotypes through complex mechanisms that are poorly understood. We previously showed that adult Caenorhabditis elegans that transiently passed through the stress-induced dauer larval stage (post-dauer adults) exhibit significant changes in gene expression profiles, chromatin states, and life history traits when compared with adults that bypassed the dauer stage (control adults). These wild-type, isogenic animals of equivalent developmental stages exhibit different signatures of molecular marks that reflect their distinct developmental trajectories. To gain insight into the mechanisms that contribute to these developmental history-dependent phenotypes, we profiled small RNAs from post-dauer and control adults by deep sequencing. RNA interference (RNAi) pathways are known to regulate genome-wide gene expression both at the chromatin and post-transcriptional level. By quantifying changes in endogenous small interfering RNA (endo-siRNA) levels in post-dauer as compared with control animals, our analyses identified a subset of genes that are likely targets of developmental history-dependent reprogramming through a complex RNAi-mediated mechanism. Mutations in specific endo-siRNA pathways affect expected gene expression and chromatin state changes for a subset of genes in post-dauer animals, as well as disrupt their increased brood size phenotype. We also find that both chromatin state and endo-siRNA distribution in dauers are unique, and suggest that remodeling in dauers provides a template for the subsequent establishment of adult post-dauer profiles. Our results indicate a role for endo-siRNA pathways as a contributing mechanism to early experience-dependent phenotypic plasticity in adults, and describe how developmental history can program adult physiology and behavior via epigenetic mechanisms.

Silencing of germline-expressed genes by DNA elimination in somatic cells

Chromatin diminution is the programmed elimination of specific DNA sequences during development. It occurs in diverse species, but the function(s) of diminution and the specificity of sequence loss remain largely unknown. Diminution in the nematode Ascaris suum occurs during early embryonic cleavages and leads to the loss of germline genome sequences and the formation of a distinct genome in somatic cells. We found that ∼43 Mb (∼13%) of genome sequence is eliminated in A. suum somatic cells, including ∼12.7 Mb of unique sequence. The eliminated sequences and location of the DNA breaks are the same in all somatic lineages from a single individual and between different individuals. At least 685 genes are eliminated. These genes are preferentially expressed in the germline and during early embryogenesis. We propose that diminution is a mechanism of germline gene regulation that specifically removes a large number of genes involved in gametogenesis and early embryogenesis.

large-scale screening for targeted knockouts in the Caenorhabditis elegans genome

The nematode Caenorhabditis elegans is a powerful model system to study contemporary biological problems. This system would be even more useful if we had mutations in all the genes of this multicellular metazoan. The combined efforts of the C. elegans Deletion Mutant Consortium and individuals within the worm community are moving us ever closer to this goal. At present, of the 20,377 protein-coding genes in this organism, 6764 genes with associated molecular lesions are either deletions or null mutations (WormBase WS220). Our three laboratories have contributed the majority of mutated genes, 6841 mutations in 6013 genes. The principal method we used to detect deletion mutations in the nematode utilizes polymerase chain reaction (PCR). More recently, we have used array comparative genome hybridization (aCGH) to detect deletions across the entire coding part of the genome and massively parallel short-read sequencing to identify nonsense, splicing, and missense defects in open reading frames. As deletion strains can be frozen and then thawed when needed, these strains will be an enduring community resource. Our combined molecular screening strategies have improved the overall throughput of our gene-knockout facilities and have broadened the types of mutations that we and others can identify. These multiple strategies should enable us to eventually identify a mutation in every gene in this multicellular organism. This knowledge will usher in a new age of metazoan genetics in which the contribution to any biological process can be assessed for all genes.

Generation of Leishmania hybrids by whole genomic DNA transformation

Genetic exchange is a powerful tool to study gene function in microorganisms. Here, we tested the feasibility of generating Leishmania hybrids by electroporating genomic DNA of donor cells into recipient Leishmania parasites. The donor DNA was marked with a drug resistance marker facilitating the selection of DNA transfer into the recipient cells. The transferred DNA was integrated exclusively at homologous locus and was as large as 45 kb. The independent generation of L. infantum hybrids with L. major sequences was possible for several chromosomal regions. Interfering with the mismatch repair machinery by inactivating the MSH2 gene enabled an increased efficiency of recombination between divergent sequences, hence favouring the selection of hybrids between species. Hybrids were shown to acquire the phenotype derived from the donor cells, as demonstrated for the transfer of drug resistance genes from L. major into L. infantum. The described method is a first step allowing the generation of in vitro hybrids for testing gene functions in a natural genomic context in the parasite Leishmania.

In vivo repair of alkylating and oxidative DNA damage in the mitochondrial and nuclear genomes of wild-type and glycosylase-deficient Caenorhabditis elegans

Base excision repair (BER) is an evolutionarily conserved DNA repair pathway that is critical for repair of many of the most common types of DNA damage generated both by endogenous metabolic pathways and exposure to exogenous stressors such as pollutants. Caenorhabditis elegans is an increasingly important model organism for the study of DNA damage-related processes including DNA repair, genotoxicity, and apoptosis, but BER is not well understood in this organism, and has not previously been measured in vivo. We report robust BER in the nuclear genome and slightly slower damage removal from the mitochondrial genome; in both cases the removal rates are comparable to those observed in mammals. However we could detect no deficiency in BER in the nth-1 strain, which carries a deletion in the only glycosylase yet described in C. elegans that repairs oxidative DNA damage. We also failed to detect increased lethality or growth inhibition in nth-1 nematodes after exposure to oxidative or alkylating damage, suggesting the existence of at least one additional as-yet undetected glycosylase.

RNA-seq analysis of the C. briggsae transcriptome

Curation of a high-quality gene set is the critical first step in genome research, enabling subsequent analyses such as ortholog assignment, cis-regulatory element finding, and synteny detection. In this project, we have reannotated the genome of Caenorhabditis briggsae, the best studied sister species of the model organism Caenorhabditis elegans. First, we applied a homology-based gene predictor genBlastG to annotate the C. briggsae genome. We then validated and further improved the C. briggsae gene annotation through RNA-seq analysis of the C. briggsae transcriptome, which resulted in the first validated C. briggsae gene set (23,159 genes), among which 7347 genes (33.9% of all genes with introns) have all of their introns confirmed. Most genes (14,812, or 68.3%) have at least one intron validated, compared with only 3.9% in the most recent WormBase release (WS228). Of all introns in the revised gene set (103,083), 61,503 (60.1%) have been confirmed. Additionally, we have identified numerous trans-splicing leaders (SL1 and SL2 variants) in C. briggsae, leading to the first genome-wide annotation of operons in C. briggsae (1105 operons). The majority of the annotated operons (564, or 51.0%) are perfectly conserved in C. elegans, with an additional 345 operons (or 31.2%) somewhat divergent. Additionally, RNA-seq analysis revealed over 10 thousand small-size assembly errors in the current C. briggsae reference genome that can be readily corrected. The revised C. briggsae genome annotation represents a solid platform for comparative genomics analysis and evolutionary studies of Caenorhabditis species.

The cell biology of schistosomes: a window on the evolution of the early metazoa

This review of schistosome cell biology has a dual purpose; its intent is to alert two separate research communities to the activities of the other. Schistosomes are by far and away the best-characterised platyhelminths, due to their medical and economic importance, but seem to be almost totally ignored by researchers on the free-living lower metazoans. Equally, in their enthusiasm for the parasitic way of life, schistosome researchers seldom pay attention to the work on free-living animals that could inform their molecular investigations. The publication of transcriptomes and/or genomes for Schistosoma mansoni and Schistosoma japonicum, the sponge Archimedon, the cnidarians Nematostella and Hydra and the planarian Schmidtea provide the raw material for comparisons. Apart from interrogation of the databases for molecular similarities, there have been differences in technical approach to these lower metazoans; widespread application of whole mount in situ hybridisation to Schmidtea contrasts with the application of targeted proteomics to schistosomes. Using schistosome cell biology as the template, the key topics of cell adhesion, development, signalling pathways, nerve and muscle, and epithelia, are reviewed, where possible interspersing comparisons with the sponge, cnidarian and planarian data. The biggest jump in the evolution of cellular capabilities appears to be in the transition from a diploblast to triploblast level of organisation associated with development of a mobile and plastic body form.

A sequence-anchored linkage map of the plant-parasitic nematode Meloidogyne hapla reveals exceptionally high genome-wide recombination

Root-knot nematodes (Meloidogyne spp.) cause major yield losses to many of the world's crops, but efforts to understand how these pests recognize and interact with their hosts have been hampered by a lack of genetic resources. Starting with progeny of a cross between inbred strains (VW8 and VW9) of Meloidogyne hapla that differed in host range and behavioral traits, we exploited the novel, facultative meiotic parthenogenic reproductive mode of this species to produce a genetic linkage map. Molecular markers were derived from SNPs identified between the sequenced and annotated VW9 genome and de novo sequence of VW8. Genotypes were assessed in 183 F2 lines. The colinearity of the genetic and physical maps supported the veracity of both. Analysis of local crossover intervals revealed that the average recombination rate is exceptionally high compared with that in other metazoans. In addition, F2 lines are largely homozygous for markers flanking crossover points, and thus resemble recombinant inbred lines. We suggest that the unusually high recombination rate may be an adaptation to generate within-population genetic diversity in this organism. This work presents the most comprehensive linkage map of a parasitic nematode to date and, together with genomic and transcript sequence resources, empowers M. hapla as a tractable model. Alongside the molecular map, these progeny lines can be used for analyses of genome organization and the inheritance of phenotypic traits that have key functions in modulating parasitism, behavior, and survival and for the eventual identification of the responsible genes.