Whole genome sequencing facilitates intragenic variant interpretation following modifier screening in C. elegans.

BACKGROUND: Intragenic modifiers (in-phase, second-site variants) are known to have dramatic effects on clinical outcomes, affecting disease attributes such as severity or age of onset. However, despite their clinical importance, the focus of many genetic screens in model systems is on the discovery of extragenic variants, with many labs still relying upon more traditional methods to identify modifiers. However, traditional methods such as PCR and Sanger sequencing can be time-intensive and do not permit a thorough understanding of the intragenic modifier effects in the context of non-isogenic genomic backgrounds. RESULTS: Here, we apply high throughput approaches to identify and understand intragenic modifiers using Caenorhabditis elegans. Specifically, we applied whole genome sequencing (WGS) to a mutagen-induced forward genetic screen to identify intragenic suppressors of a temperature-sensitive zyg-1(it25) allele in C. elegans. ZYG-1 is a polo kinase that is important for centriole function and cell divisions, and mutations that truncate its human orthologue, PLK4, have been associated with microcephaly. Combining WGS and CRISPR/Cas9, we rapidly identify intragenic modifiers, show that these variants are distributed non-randomly throughout zyg-1 and that genomic context plays an important role on phenotypic outcomes. CONCLUSIONS: Ultimately, our work shows that WGS facilitates high-throughput identification of intragenic modifiers in clinically relevant genes by reducing hands-on research time and overall costs and by allowing thorough understanding of the intragenic phenotypic effects in the context of different genetic backgrounds.

Deciphering complex genome rearrangements in C. elegans using short-read whole genome sequencing.

Genomic rearrangements cause congenital disorders, cancer, and complex diseases in human. Yet, they are still understudied in rare diseases because their detection is challenging, despite the advent of whole genome sequencing (WGS) technologies. Short-read (srWGS) and long-read WGS approaches are regularly compared, and the latter is commonly recommended in studies focusing on genomic rearrangements. However, srWGS is currently the most economical, accurate, and widely supported technology. In Caenorhabditis elegans (C. elegans), such variants, induced by various mutagenesis processes, have been used for decades to balance large genomic regions by preventing chromosomal crossover events and allowing the maintenance of lethal mutations. Interestingly, those chromosomal rearrangements have rarely been characterized on a molecular level. To evaluate the ability of srWGS to detect various types of complex genomic rearrangements, we sequenced three balancer strains using short-read Illumina technology. As we experimentally validated the breakpoints uncovered by srWGS, we showed that, by combining several types of analyses, srWGS enables the detection of a reciprocal translocation (eT1), a free duplication (sDp3), a large deletion (sC4), and chromoanagenesis events. Thus, applying srWGS to decipher real complex genomic rearrangements in model organisms may help designing efficient bioinformatics pipelines with systematic detection of complex rearrangements in human genomes.

Author Correction: Caenorhabditis elegans is a useful model for anthelmintic discovery.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The confounding effects of high genetic diversity on the determination and interpretation of differential gene expression analysis in the parasitic nematode Haemonchus contortus.

Differential expression analysis between parasitic nematode strains is commonly used to implicate candidate genes in anthelmintic resistance or other biological functions. We have tested the hypothesis that the high genetic diversity of an organism such as Haemonchus contortus could complicate such analyses. First, we investigated the extent to which sequence polymorphism affects the reliability of differential expression analysis between the genetically divergent H. contortus strains MHco3(ISE), MHco4(WRS) and MHco10(CAVR). Using triplicates of 20 adult female worms from each population isolated under parallel experimental conditions, we found that high rates of sequence polymorphism in RNAseq reads were associated with lower efficiency read mapping to gene models under default TopHat2 parameters, leading to biased estimates of inter-strain differential expression. We then showed it is possible to largely compensate for this bias by optimising the read mapping single nucleotide polymorphism (SNP) allowance and filtering out genes with particularly high single nucleotide polymorphism rates. Once the sequence polymorphism biases were removed, we then assessed the genuine transcriptional diversity between the strains, finding ≥824 differentially expressed genes across all three pairwise strain comparisons. This high level of inter-strain transcriptional diversity not only suggests substantive inter-strain phenotypic variation but also highlights the difficulty in reliably associating differential expression of specific genes with phenotypic differences. To provide a practical example, we analysed two gene families of potential relevance to ivermectin drug resistance; the ABC transporters and the ligand-gated ion channels (LGICs). Over half of genes identified as differentially expressed using default TopHat2 parameters were shown to be an artifact of sequence polymorphism differences. This work illustrates the need to account for sequence polymorphism in differential expression analysis. It also demonstrates that a large number of genuine transcriptional differences can occur between H. contortus strains and these must be considered before associating the differential expression of specific genes with phenotypic differences between strains.

Caenorhabditis elegans is a useful model for anthelmintic discovery.

Parasitic nematodes infect one quarter of the world's population and impact all humans through widespread infection of crops and livestock. Resistance to current anthelmintics has prompted the search for new drugs. Traditional screens that rely on parasitic worms are costly and labour intensive and target-based approaches have failed to yield novel anthelmintics. Here, we present our screen of 67,012 compounds to identify those that kill the non-parasitic nematode Caenorhabditis elegans. We then rescreen our hits in two parasitic nematode species and two vertebrate models (HEK293 cells and zebrafish), and identify 30 structurally distinct anthelmintic lead molecules. Genetic screens of 19 million C. elegans mutants reveal those nematicides for which the generation of resistance is and is not likely. We identify the target of one lead with nematode specificity and nanomolar potency as complex II of the electron transport chain. This work establishes C. elegans as an effective and cost-efficient model system for anthelmintic discovery.

Developmental plasticity and the evolution of parasitism in an unusual nematode, Parastrongyloides trichosuri.

BACKGROUND: Parasitism is an important life history strategy in many metazoan taxa. This is particularly true of the Phylum Nematoda, in which parasitism has evolved independently at least nine times. The apparent ease with which parasitism has evolved amongst nematodes may, in part, be due to a feature of nematode development acting as a pre-adaptation for the transition from a free-living to a parasitic life history. One candidate pre-adaptive feature for evolution in terrestrial nematodes is the dauer larva, a developmentally arrested morph formed in response to environmental signals. RESULTS: We investigated the role of dauer development in the nematode, Parastrongyloides trichosuri, which has retained a complete free-living life cycle in addition to a life cycle as a mammalian gastrointestinal parasite. We show that the developmental switch between these life histories is sensitive to the same environmental cues as dauer arrest in free-living nematodes, including sensitivity to a chemical cue produced by the free-living stages. Furthermore, we show that genetic variation for the sensitivity of the cue(s) exists in natural populations of P. trichosuri, such that we derived inbred lines that were largely insensitive to the cue and other lines that were supersensitive to the cue. CONCLUSIONS: For this parasitic clade, and perhaps more widely in the phylum, the evolution of parasitism co-opted the dauer switch of a free-living ancestor. This lends direct support to the hypothesis that the switch to developmental arrest in the dauer larva acted as a pre-adaptation for the evolution of parasitism, and suggests that the sensory transduction machinery downstream of the cue may have been similarly co-opted and modified.

Functional characterization of Caenorhabditis elegans heteromeric amino acid transporters.

Mammalian heteromeric amino acid transporters (HATs) are composed of a multi-transmembrane spanning catalytic protein covalently associated with a type II glycoprotein (e.g. 4F2hc, rBAT) through a disulfide bond. Caenorhabditis elegans has nine genes encoding close homologues of the HAT catalytic proteins. Three of these genes (designated AAT-1 to AAT-3) have a much higher degree of similarity to the mammalian homologues than the other six, including the presence of a cysteine residue at the position known to form a disulfide bridge to the glycoprotein partner in mammalian HATs. C. elegans also has two genes encoding homologues of the heteromeric amino acid transporter type II glycoprotein subunits (designated ATG-1 and ATG-2). Both ATG, and/or AAT-1, -2, -3 proteins were expressed in Xenopus oocytes and tested for amino acid transport function. This screen revealed that AAT-1 and AAT-3 facilitate amino acid transport when expressed together with ATG-2 but not with ATG-1 or the mammalian type II glycoproteins 4F2hc and rBAT. AAT-1 and AAT-3 covalently bind to both C. elegans ATG glycoproteins, but only the pairs with ATG-2 traffic to the oocyte surface. Both of these functional, surface-expressed C. elegans HATs transport most neutral amino acids and display the highest transport rate for l-Ala and l-Ser (apparent K(m) 100 microm range). Similar to their mammalian counterparts, the C. elegans HATs function as (near) obligatory amino acid exchangers. Taken together, this study demonstrates that the heteromeric structure and the amino acid exchange function of HATs have been conserved throughout the evolution of nematodes to mammals.

Aromatic amino acid transporter AAT-9 of Caenorhabditis elegans localizes to neurons and muscle cells.

The Caenorhabditis elegans genome encodes nine homologues of mammalian glycoprotein-associated amino acid transporters. Two of these C. elegans proteins (AAT-1 and AAT-3) have been shown to function as catalytic subunits (light chains) of heteromeric amino acid transporters. These proteins need to associate with a glycoprotein heavy chain subunit (ATG-2) to reach the cell surface in a manner similar to that of their mammalian homologues. AAT-1 and AAT-3 contain a cysteine residue in the second putative extracellular loop through which a disulfide bridge can form with a heavy chain. In contrast, six C. elegans members of this family (AAT-4 to AAT-9) lack such a cysteine residue. We show here that one of these transporter proteins, AAT-9, reaches the cell surface in Xenopus oocytes without an exogenous heavy chain and that it functions as an exchanger of aromatic amino acids. Two-electrode voltage clamp experiments demonstrate that AAT-9 displays a substrate-activated conductance. Immunofluorescence shows that it is expressed close to the pharyngeal bulbs within C. elegans neurons. The selective expression of an aat-9 promoter-green fluorescent protein construct in several neurons of this region and in wall muscle cells around the mouth supports and extends these localization data. Taken together, the results show that AAT-9 is expressed in excitable cells of the nematode head and pharynx in which it may provide a pathway for aromatic amino acid transport.