Spirocerca lupi draft genome, vaccine and anthelmintic targets

Spirocerca lupi is a parasitic nematode affecting predominantly domestic dogs. It causes spirocercosis, a disease that is often fatal. The assembled draft genome of S. lupi consists of 13,627 predicted protein-coding genes and is approximately 150 Mb in length. Several known anthelmintic gene targets such as for β-Tubulin, glutamate, and GABA receptors as well as known vaccine gene targets such as cysteine protease inhibitor and cytokines were identified in S. lupi by comparing orthologs of C. elegans anthelmintic gene targets as well as orthologs to known vaccine candidates. New anthelmintic targets were predicted through an inclusion-exclusion strategy and new vaccine targets were predicted through an immunoinformatics approach. New anthelminthic targets include DNA-directed RNA polymerases, chitin synthase, polymerases, and other enzymes. New vaccine targets include cuticle collagens. These gene targets provide a starting platform for new drug identification and vaccine design.

Genome-wide RNA interference of the nhr gene family in barber's pole worm identified members crucial for larval viability in vitro

Nuclear hormone receptors (NHRs) are emerging target candidates against nematode infection and resistance. However, there is a lack of comprehensive information on NHR-coding genes in parasitic nematodes. In this study, we curated the nhr gene family for 60 major parasitic nematodes from humans and animals. Compared with the free-living model organism Caenorhabditis elegans, a remarkable contraction of the nhr family was revealed in parasitic species, with genetic diversification and conservation unveiled among nematode Clades I (10-13), III (16-42), IV (33-35) and V (25-64). Using an in vitro biosystem, we demonstrated that 40 nhr genes in a blood-feeding nematode Haemonchus contortus (clade V; barber's pole worm) were responsive to host serum and one nhr gene (i.e., nhr-64) was consistently stimulated by anthelmintics (i.e., ivermectin, thiabendazole and levamisole); Using a high-throughput RNA interference platform, we knocked down 43 nhr genes of H. contortus and identified at least two genes that are required for the viability (i.e., nhr-105) and development (i.e., nhr-17) of the infective larvae of this parasitic nematode in vitro. Harnessing this preliminary functional atlas of nhr genes for H. contortus will prime the biological studies of this gene family in nematode genetics, infection, and anthelmintic metabolism within host animals, as well as the promising discovery of novel intervention targets.

Chromosome-level genome assembly of the cereal cyst nematode Heterodera flipjevi

As an economically important plant parasitic nematode (PPN), Heterodera filipjevi causes great damage on wheat, and now it was widely recorded in many countries. While multiple genomes of PPNs have been published, high-quality genome assembly and annotation on H. filipjevi have yet to be performed. This study presents a chromosome-scale genome assembly and annotation for H. filipjevi, utilizing a combination of Illumina short-read, PacBio long-read, and Hi-C sequencing technologies. The genome consists of 9 pseudo-chromosomes that contain 134.19 Mb of sequence, with a scaffold N50 length of 11.88 Mb. In total, 10,036 genes were annotated, representing 75.20% of the total predicted protein-coding genes. Our study provides the first chromosome-scale genome for H. filipjevi, which is also the inaugural high-quality genome of cereal cyst nematodes (CCNs). It provides a valuable genomic resource for further biological research and pest management of cereal cyst nematodes disease.

Phased chromosome-scale genome assembly of an asexual, allopolyploid root-knot nematode reveals complex subgenomic structure

We present the chromosome-scale genome assembly of the allopolyploid root-knot nematode Meloidogyne javanica. We show that the M. javanica genome is predominantly allotetraploid, comprising two subgenomes, A and B, that most likely originated from hybridisation of two ancestral parental species. The assembly was annotated using full-length non-chimeric transcripts, comparison to reference databases, and ab initio prediction techniques, and the subgenomes were phased using ancestral k-mer spectral analysis. Subgenome B appears to show fission of chromosomal contigs, and while there is substantial synteny between subgenomes, we also identified regions lacking synteny that may have diverged in the ancestral genomes prior to or following hybridisation. This annotated and phased genome assembly forms a significant resource for understanding the origins and genetics of these globally important plant pathogens.

The impact of differential transposition activities of autonomous and nonautonomous hAT transposable elements on genome architecture and gene expression in Caenorhabditis inopinata

Transposable elements are DNA sequences capable of moving within genomes and significantly influence genomic evolution. The nematode Caenorhabditis inopinata exhibits a much higher transposable element copy number than its sister species, Caenorhabditis elegans. In this study, we identified a novel autonomous transposable element belonging to the hAT superfamily from a spontaneous transposable element-insertion mutant in C. inopinata and named this transposon Ci-hAT1. Further bioinformatic analyses uncovered 3 additional autonomous hAT elements-Ci-hAT2, Ci-hAT3, and Ci-hAT4-along with over 1,000 copies of 2 nonautonomous miniature inverted-repeat transposable elements, mCi-hAT1 and mCi-hAT4, likely derived from Ci-hAT1 and Ci-hAT4 through internal deletion. We tracked at least 3 sequential transpositions of Ci-hAT1 over several years. However, the transposition rates of the other 3 autonomous hAT elements were lower, suggesting varying activity levels. Notably, the distribution patterns of the 2 miniature inverted-repeat transposable element families differed significantly: mCi-hAT1 was primarily located in the chromosome arms, a pattern observed in the transposable elements of other Caenorhabditis species, whereas mCi-hAT4 was more evenly distributed across chromosomes. Additionally, interspecific transcriptome analysis indicated that C. inopinata genes with upstream or intronic these miniature inverted-repeat transposable element insertions tend to be more highly expressed than their orthologous genes in C. elegans. These findings highlight the significant role of de-silenced transposable elements in driving the evolution of genomes and transcriptomes, leading to species-specific genetic diversity.

WormBase 2024: status and transitioning to Alliance infrastructure

WormBase has been the major repository and knowledgebase of information about the genome and genetics of Caenorhabditis elegans and other nematodes of experimental interest for over 2 decades. We have 3 goals: to keep current with the fast-paced C. elegans research, to provide better integration with other resources, and to be sustainable. Here, we discuss the current state of WormBase as well as progress and plans for moving core WormBase infrastructure to the Alliance of Genome Resources (the Alliance). As an Alliance member, WormBase will continue to interact with the C. elegans community, develop new features as needed, and curate key information from the literature and large-scale projects.

Recombination in bdelloid rotifer genomes: asexuality, transfer and stress

Bdelloid rotifers constitute a class of microscopic animals living in freshwater habitats worldwide. Several strange features of bdelloids have drawn attention: their ability to tolerate desiccation and other stresses, a lack of reported males across the clade despite centuries of study, and unusually high numbers of horizontally acquired, non-metazoan genes. Genome sequencing is transforming our understanding of their lifestyle and its consequences, while in turn providing wider insights about recombination and genome organisation in animals. Many questions remain, not least how to reconcile apparent genomic signatures of sex with the continued absence of reported males, why bdelloids have so many horizontally acquired genes, and how their remarkable ability to survive stress interacts with recombination and other genomic processes.

Genome structure and population genomics of the canine heartworm Dirofilaria immitis

The heartworm, Dirofilaria immitis, is a filarial parasitic nematode responsible for significant morbidity and mortality in wild and domesticated canids. Resistance to macrocyclic lactone drug prevention represents a significant threat to parasite control and has prompted investigations to understand the genetic determinants of resistance. This study aimed to improve the genomic resources of D. immitis to enable a more precise understanding of how genetic variation is distributed within and between parasite populations worldwide, which will inform the likelihood and rate by which parasites, and in turn, resistant alleles, might spread. We have guided the scaffolding of a recently published genome assembly for D. immitis (ICBAS_JMDir_1.0) using the chromosomal-scale reference genomes of Brugia malayi and Onchocerca volvulus, resulting in an 89.5 Mb assembly composed of four autosomal- and one sex-linked chromosomal-scale scaffolds representing 99.7% of the genome. Publicly available and new whole-genome sequencing data from 32 D. immitis samples from Australia, Italy and the USA were assessed using principal component analysis, nucleotide diversity (Pi) and absolute genetic divergence (Dxy) to characterise the global genetic structure and measure within- and between-population diversity. These population genetic analyses revealed broad-scale genetic structure among globally diverse samples and differences in genetic diversity between populations; however, fine-scale subpopulation analysis was limited and biased by differences between sample types. Finally, we mapped single nucleotide polymorphisms previously associated with macrocyclic lactone resistance in the new genome assembly, revealing the physical linkage of high-priority variants on chromosome 3, and determined their frequency in the studied populations. This new chromosomal assembly for D. immitis now allows for a more precise investigation of selection on genome-wide genetic variation and will enhance our understanding of parasite transmission and the spread of genetic variants responsible for resistance to treatment.

Unraveling the Unusual Subgenomic Organization in the Neopolyploid Free-Living Flatworm Macrostomum lignano

Whole genome duplication (WGD) is an evolutionary event resulting in a redundancy of genetic material. Different mechanisms of WGD, allo- or autopolyploidization, lead to distinct evolutionary trajectories of newly formed polyploids. Genome studies on such species are important for understanding the early stages of genome evolution. However, assembling neopolyploid is a challenging task due to the presence of 2 homologous (or homeologous) chromosome sets and therefore the existence of the extended paralogous regions in its genome. Post-WGD evolution of polyploids includes cytogenetic diploidization leading to the formation of species, whose polyploid origin might be hidden by disomic inheritance. Earlier we uncovered the hidden polyploid origin of the free-living flatworms of the genus Macrostomum (Macrostomum lignano, M. janickei, and M. mirumnovem). Cytogenetic diploidization in these species is accompanied by intensive chromosomal rearrangements including chromosomes fusions. In this study, we unravel the M. lignano genome organization through generation and sequencing of 2 sublines of the commonly used inbred line of M. lignano (called DV1) differing only in a copy number of the largest chromosome (MLI1). Using nontrivial assembly free comparative analysis of their genomes, we deciphered DNA sequences belonging to MLI1 and validated them by sequencing the pool of microdissected MLI1. Here we presented the uncommon mechanism of genome rediplodization of M. lignano, which consists of (i) presence of 3 subgenomes, which emerged via formation of large fused chromosomes and its variants, and (ii) sustaining their heterozygosity through inter- and intrachromosomal rearrangements.

Expression of Hairpin-Enriched Mitochondrial DNA in Two Hairworm Species (Nematomorpha)

Nematomorpha (hairworms) is a phylum of parasitic ecdysozoans, best known for infecting arthropods and guiding their hosts toward water, where the parasite can complete its life cycle. Over 350 species of nematomorphs have been described, yet molecular data for the group remain scarce. The few available mitochondrial genomes of nematomorphs are enriched with long inverted repeats, which are embedded in the coding sequences of their genes-a remarkably unusual feature exclusive to this phylum. Here, we obtain and annotate the repeats in the mitochondrial genome of another nematomorph species-Parachordodes pustulosus. Using genomic and transcriptomic libraries, we investigate the impact of inverted repeats on the read coverage of the mitochondrial genome. Pronounced drops in the read coverage coincide with regions containing long inverted repeats, denoting the 'blind spots' of short-fragment sequencing libraries. Phylogenetic inference with the novel data reveals multiple disagreements between the traditional system of Nematomorpha and molecular data, rendering several genera paraphyletic, including Parachordodes.

Population differentiation and epidemic tracking of Bursaphelenchus xylophilus in China based on chromosome-level assembly and whole-genome sequencing data

BACKGROUND

Bursaphelenchus xylophilus, the pinewood nematode, kills millions of pine trees worldwide every year, and causes enormous economic and ecological losses. Despite extensive research on population variation, there is little understanding of the population-wide variation spectrum in China.

RESULTS

We sequenced an inbred B. xylophilus strain using Pacbio+Illumina+Bionano+Hi-C and generated a chromosome-level assembly (AH1) with six chromosomes of 77.1 Mb (chromosome N50: 12 Mb). The AH1 assembly shows very high continuity and completeness, and contains novel genes with potentially important functions compared with previous assemblies. Subsequently, we sequenced 181 strains from China and the USA and found ~7.8 million single nucleotide polymorphisms (SNPs). Analysis shows that the B. xylophilus population in China can be divided into geographically bounded subpopulations with severe cross-infection and potential migrations. In addition, distribution of B. xylophilus is dominated by temperature zones while geographically associated SNPs are mainly located on adaptation related GPCR gene families, suggesting the nematode has been evolving to adapt to different temperatures. A machine-learning based epidemic tracking method has been established to predict their geographical origins, which can be applied to any other species.

CONCLUSION

Our study provides the community with the first high-quality chromosome-level assembly which includes a comprehensive catalogue of genetic variations. It provides insights into population structure and effective tracking method for this invasive species, which facilitates future studies to address a variety of applied, genomic and evolutionary questions in B. xylophilus as well as related species.

Approaches for CRISPR/Cas9 Genome Editing in C. elegans

The clustered, regularly interspaced, short, palindromic repeat (CRISPR)-associated (CAS) nuclease Cas9 has been used in many organisms to generate specific mutations and transgene insertions. Here we describe our most up-to-date protocols using the S. pyogenes Cas9 in C. elegans that provides a convenient and effective approach for making heritable changes to the worm genome. We present several considerations when deciding which strategy best suits the needs of the experiment.

Genome-wide transcriptome analysis of the early developmental stages of Echinococcus granulosus protoscoleces reveals extensive alternative splicing events in the spliceosome pathway

BACKGROUND

The complex life cycle of Echinococcus granulosus involves numerous environmental conditions within different intermediate and definitive hosts. This requires adaptation at different levels of transcript regulation. Alternative splicing (AS) and the related cellular functions as one of the major fields of post-genomics has been poorly studied in tapeworms. In the present study, we investigated AS events and their potential biological effects in E. granulosus.

METHODS

Whole transcriptome sequencing data of four groups of protoscoleces were prepared for RNA-seq library construction. Fresh protoscoleces were either used as non-induced controls (NT group) or incubated for 15 min with pepsin (PEP group) and cultivated in a biphasic medium for 12 and 24 h (12 and 24 h groups). The frequency and different types of AS events were identified using rMATS software. Functional annotations and gene ontology of differential AS (DAS) genes were performed using Blast2GO software. AS events were experimentally validated by PCR on the protoscolex cDNAs using specific primers for each gene.

RESULTS

At least one AS event was found in 38.1% of the genes (3904 out of 10,245) in the protoscoleces during early strobilar development. The genes were associated primarily with cellular and metabolic processes and binding and catalytic activity. KEGG pathway analysis of DAS events revealed a number of genes belonging to different components of the spliceosome complex. These genes tended to belong to common SR proteins, U1-related factors, U2-related factors, complex A-specific factors and other splicing-related proteins.

CONCLUSIONS

The high number of AS events in the transcriptome regulatory mechanisms indicates the essential rapid molecular processes required by the parasite for adaptation in different environments.

High-efficiency CRISPR gene editing in C. elegans using Cas9 integrated into the genome

Gene editing in C. elegans using plasmid-based CRISPR reagents requires microinjection of many animals to produce a single edit. Germline silencing of plasmid-borne Cas9 is a major cause of inefficient editing. Here, we present a set of C. elegans strains that constitutively express Cas9 in the germline from an integrated transgene. These strains markedly improve the success rate for plasmid-based CRISPR edits. For simple, short homology arm GFP insertions, 50-100% of injected animals typically produce edited progeny, depending on the target locus. Template-guided editing from an extrachromosomal array is maintained over multiple generations. We have built strains with the Cas9 transgene on multiple chromosomes. Additionally, each Cas9 locus also contains a heatshock-driven Cre recombinase for selectable marker removal and a bright fluorescence marker for easy outcrossing. These integrated Cas9 strains greatly reduce the workload for producing individual genome edits.

Hookworm infection in central China: morphological and molecular diagnosis

BACKGROUND

Necator americanus is one of the major etiological agents of human ancylostomiasis. Historically, the epidemiology of ancylostomiasis in Henan Province of central China and the molecular characteristics of N. americanus have been poorly understood.

METHODS

In this study, we report a case of ancylostomiasis in Zhengzhou city of Henan Province. We also review the epidemiology of ancylostomiasis in Henan Province from 1949 to 2020. In addition, the complete mitochondrial (mt) genome of one clinical isolate is fully characterized using Illumina sequencing. All available mt genomes of hookworms in GenBank were included to reconstruct the phylogeny using both maximum likelihood (ML) and Bayesian inference (BI) methods.

RESULTS

A total of three worms were collected from the patient. These worms were identified as N. americanus based on morphological characteristics as well as confirmed by genotyping with the barcoding gene cox1. Although ancylostomiasis cases have dropped substantially in recent years, hookworm infection is still a public health problem in underdeveloped areas and remote rural areas in Henan Province. The mt genome features of the N. americanus contained 12 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and a major non-coding region. The nad1 gene showed high sequence variability among isolates, which is worth considering for future genetic studies of N. americanus. Phylogenetic analyses support the monophyly of hookworm isolates from different hosts and distinct geographical locations.

CONCLUSIONS

The mt genome of N. americanus presented here will serve as a useful data set for studying population genetics and phylogenetic relationships of hookworms. Positive measures for preventing and controlling ancylostomiasis are required by both health services and individuals in Henan Province.

C. elegans genome-wide analysis reveals DNA repair pathways that act cooperatively to preserve genome integrity upon ionizing radiation

Ionizing radiation (IR) is widely used in cancer therapy and accidental or environmental exposure is a major concern. However, little is known about the genome-wide effects IR exerts on germ cells and the relative contribution of DNA repair pathways for mending IR-induced lesions. Here, using C. elegans as a model system and using primary sequencing data from our recent high-level overview of the mutagenic consequences of 11 genotoxic agents, we investigate in detail the genome-wide mutagenic consequences of exposing wild-type and 43 DNA repair and damage response defective C. elegans strains to a Caesium (Cs-137) source, emitting γ-rays. Cs-137 radiation induced single nucleotide variants (SNVs) at a rate of ~1 base substitution per 3 Gy, affecting all nucleotides equally. In nucleotide excision repair mutants, this frequency increased 2-fold concurrently with increased dinucleotide substitutions. As observed for DNA damage induced by bulky DNA adducts, small deletions were increased in translesion polymerase mutants, while base changes decreased. Structural variants (SVs) were augmented with dose, but did not arise with significantly higher frequency in any DNA repair mutants tested. Moreover, 6% of all mutations occurred in clusters, but clustering was not significantly altered in any DNA repair mutant background. Our data is relevant for better understanding how DNA repair pathways modulate IR-induced lesions.

X-treme loss of sequence diversity linked to neo-X chromosomes in filarial nematodes

The sequence diversity of natural and laboratory populations of Brugia pahangi and Brugia malayi was assessed with Illumina resequencing followed by mapping in order to identify single nucleotide variants and insertions/deletions. In natural and laboratory Brugia populations, there is a lack of sequence diversity on chromosome X relative to the autosomes (πX/πA = 0.2), which is lower than the expected (πX/πA = 0.75). A reduction in diversity is also observed in other filarial nematodes with neo-X chromosome fusions in the genera Onchocerca and Wuchereria, but not those without neo-X chromosome fusions in the genera Loa and Dirofilaria. In the species with neo-X chromosome fusions, chromosome X is abnormally large, containing a third of the genetic material such that a sizable portion of the genome is lacking sequence diversity. Such profound differences in genetic diversity can be consequential, having been associated with drug resistance and adaptability, with the potential to affect filarial eradication.

High-quality reference genome of Fasciola gigantica: Insights into the genomic signatures of transposon-mediated evolution and specific parasitic adaption in tropical regions

Fasciola gigantica and Fasciola hepatica are causative pathogens of fascioliasis, with the widest latitudinal, longitudinal, and altitudinal distribution; however, among parasites, they have the largest sequenced genomes, hindering genomic research. In the present study, we used various sequencing and assembly technologies to generate a new high-quality Fasciola gigantica reference genome. We improved the integration of gene structure prediction, and identified two independent transposable element expansion events contributing to (1) the speciation between Fasciola and Fasciolopsis during the Cretaceous-Paleogene boundary mass extinction, and (2) the habitat switch to the liver during the Paleocene-Eocene Thermal Maximum, accompanied by gene length increment. Long interspersed element (LINE) duplication contributed to the second transposon-mediated alteration, showing an obvious trend of insertion into gene regions, regardless of strong purifying effect. Gene ontology analysis of genes with long LINE insertions identified membrane-associated and vesicle secretion process proteins, further implicating the functional alteration of the gene network. We identified 852 predicted excretory/secretory proteins and 3300 protein-protein interactions between Fasciola gigantica and its host. Among them, copper/zinc superoxide dismutase genes, with specific gene copy number variations, might play a central role in the phase I detoxification process. Analysis of 559 single-copy orthologs suggested that Fasciola gigantica and Fasciola hepatica diverged at 11.8 Ma near the Middle and Late Miocene Epoch boundary. We identified 98 rapidly evolving gene families, including actin and aquaporin, which might explain the large body size and the parasitic adaptive character resulting in these liver flukes becoming epidemic in tropical and subtropical regions.

Linked surveillance and genetic data uncovers programmatically relevant geographic scale of Guinea worm transmission in Chad

Background

Guinea worm (Dracunculus medinensis) was detected in Chad in 2010 after a supposed ten-year absence, posing a challenge to the global eradication effort. Initiation of a village-based surveillance system in 2012 revealed a substantial number of dogs infected with Guinea worm, raising questions about paratenic hosts and cross-species transmission.

Methodology/principal findings

We coupled genomic and surveillance case data from 2012-2018 to investigate the modes of transmission between dog and human hosts and the geographic connectivity of worms. Eighty-six variants across four genes in the mitochondrial genome identified 41 genetically distinct worm genotypes. Spatiotemporal modeling revealed worms with the same genotype (‘genetically identical’) were within a median range of 18.6 kilometers of each other, but largely within approximately 50 kilometers. Genetically identical worms varied in their degree of spatial clustering, suggesting there may be different factors that favor or constrain transmission. Each worm was surrounded by five to ten genetically distinct worms within a 50 kilometer radius. As expected, we observed a change in the genetic similarity distribution between pairs of worms using variants across the complete mitochondrial genome in an independent population.

Conclusions/significance

In the largest study linking genetic and surveillance data to date of Guinea worm cases in Chad, we show genetic identity and modeling can facilitate the understanding of local transmission. The co-occurrence of genetically non-identical worms in quantitatively identified transmission ranges highlights the necessity for genomic tools to link cases. The improved discrimination between pairs of worms from variants identified across the complete mitochondrial genome suggests that expanding the number of genomic markers could link cases at a finer scale. These results suggest that scaling up genomic surveillance for Guinea worm may provide additional value for programmatic decision-making critical for monitoring cases and intervention efficacy to achieve elimination.

The global eradication effort for Guinea worm disease has dramatically decreased the global burden of the disease and enabled 187 countries to be certified by the World Health Organization to be free of endemic transmission. Despite this progress, several countries continue to have endemic transmission. In Chad, a long absence of reported cases was interrupted with the identification of new Guinea worm cases, prompting a substantial scale up of surveillance and intervention efforts. Here, we study the value of increasing genomic surveillance as a tool for programmatic evaluation of surveillance and intervention efforts in Chad. Linking surveillance and genomic samples, parsimonious spatial models help reveal a consistent geographic clustering of similar genetic sequences across Chad. We also demonstrate that expanding the sequencing can offer better resolution for distinguishing Guinea worm samples. In this retrospective study, we found evidence that scaling up genomic surveillance can be an important monitoring and evaluation tool for the eradication program in Chad.

Genome-Wide Characterization and Identification of Long Non-Coding RNAs during the Molting Process of a Spider Mite, Panonychus citri

Molting is essential for arthropods to grow. As one of the important arthropod pests in agriculture, key spider mite species (Tetranychus and Panonychus) can normally molt three times from the larva to adult stage within a week. This physiological strategy results in the short lifecycle of spider mites and difficulties in their control in the field. Long non-coding RNAs (lncRNAs) regulate transcriptional editing, cellular function, and biological processes. Thus, analysis of the lncRNAs in the spider mite molting process may provide new insights into their roles in the molting mechanism. For this purpose, we used high-throughput RNA-seq to examine the expression dynamics of lncRNAs and mRNAs in the molting process of different development stages in Panonychus citri. We identified 9199 lncRNAs from 18 transcriptomes. Analysis of the lncRNAs suggested that they were shorter and had fewer exons and transcripts than mRNAs. Among these, 356 lncRNAs were differentially expressed during three molting processes: late larva to early protonymph, late protonymph to early deutonymph, and late deutonymph to early adult. A time series profile analysis of differentially expressed lncRNAs showed that 77 lncRNAs were clustered into two dynamic expression profiles (Pattern a and Pattern c), implying that lncRNAs were involved in the molting process of spider mites. Furthermore, the lncRNA-mRNA co-expression networks showed that several differentially expressed hub lncRNAs were predicted to be functionally associated with typical molting-related proteins, such as cuticle protein and chitin biosynthesis. These data reveal the potential regulatory function of lncRNAs in the molting process and provide datasets for further analysis of lncRNAs and mRNAs in spider mites.

DYF-4 regulates patched-related/DAF-6-mediated sensory compartment formation in C. elegans

Coordination of neurite extension with surrounding glia development is critical for neuronal function, but the underlying molecular mechanisms remain poorly understood. Through a genome-wide mutagenesis screen in C. elegans, we identified dyf-4 and daf-6 as two mutants sharing similar defects in dendrite extension. DAF-6 encodes a glia-specific patched-related membrane protein that plays vital roles in glial morphogenesis. We cloned dyf-4 and found that DYF-4 encodes a glia-secreted protein. Further investigations revealed that DYF-4 interacts with DAF-6 and functions in a same pathway as DAF-6 to regulate sensory compartment formation. Furthermore, we demonstrated that reported glial suppressors of daf-6 could also restore dendrite elongation and ciliogenesis in both dyf-4 and daf-6 mutants. Collectively, our data reveal that DYF-4 is a regulator for DAF-6 which promotes the proper formation of the glial channel and indirectly affects neurite extension and ciliogenesis.

In C. elegans sensory organ, the ciliated neuronal endings are wrapped in a luminal channel formed by glial cells, forming a specialized sensory compartment critical for sensory activity. Coordination of glial channel formation, dendritic tip anchoring and ciliogenesis are critical for the formation of a functional sensory compartment. DAF-6, a patched-related glial membrane protein, was reported to play an important role in glial channel morphogenesis, but the molecular function and regulatory mechanism of DAF-6 remain poorly understood. Here, we found that DYF-4, a glia-secreted protein, interacts and colocalizes with DAF-6, and functions in a same pathway as DAF-6 to regulate sensory compartment formation. We propose that DYF-4 is a novel regulator for DAF-6 to control sensory compartment formation.

Linking genome size variation to population phenotypic variation within the rotifer, Brachionus asplanchnoidis

Eukaryotic organisms usually contain much more genomic DNA than expected from their biological complexity. In explaining this pattern, selection-based hypotheses suggest that genome size evolves through selection acting on correlated life history traits, implicitly assuming the existence of phenotypic effects of (extra) genomic DNA that are independent of its information content. Here, we present conclusive evidence of such phenotypic effects within a well-mixed natural population that shows heritable variation in genome size. We found that genome size is positively correlated with body size, egg size, and embryonic development time in a population of the monogonont rotifer Brachionus asplanchnoidis. The effect on embryonic development time was mediated partly by an indirect effect (via egg size), and a direct effect, the latter indicating an increased replication cost of the larger amounts of DNA during mitosis. Our results suggest that selection-based change of genome size can operate in this population, provided it is strong enough to overcome drift or mutational change of genome size.

Protection of the C. elegans germ cell genome depends on diverse DNA repair pathways during normal proliferation

Maintaining genome integrity is particularly important in germ cells to ensure faithful transmission of genetic information across generations. Here we systematically describe germ cell mutagenesis in wild-type and 61 DNA repair mutants cultivated over multiple generations. ~44% of the DNA repair mutants analysed showed a >2-fold increased mutagenesis with a broad spectrum of mutational outcomes. Nucleotide excision repair deficiency led to higher base substitution rates, whereas polh-1(Polη) and rev-3(Polζ) translesion synthesis polymerase mutants resulted in 50-400 bp deletions. Signatures associated with defective homologous recombination fall into two classes: 1) brc-1/BRCA1 and rad-51/RAD51 paralog mutants showed increased mutations across all mutation classes, 2) mus-81/MUS81 and slx-1/SLX1 nuclease, and him-6/BLM, helq-1/HELQ or rtel-1/RTEL1 helicase mutants primarily accumulated structural variants. Repetitive and G-quadruplex sequence-containing loci were more frequently mutated in specific DNA repair backgrounds. Tandem duplications embedded in inverted repeats were observed in helq-1 helicase mutants, and a unique pattern of 'translocations' involving homeologous sequences occurred in rip-1 recombination mutants. atm-1/ATM checkpoint mutants harboured structural variants specifically enriched in subtelomeric regions. Interestingly, locally clustered mutagenesis was only observed for combined brc-1 and cep-1/p53 deficiency. Our study provides a global view of how different DNA repair pathways contribute to prevent germ cell mutagenesis.

Parallel genetics of regulatory sequences using scalable genome editing in vivo

How regulatory sequences control gene expression is fundamental for explaining phenotypes in health and disease. Regulatory elements must ultimately be understood within their genomic environment and development- or tissue-specific contexts. Because this is technically challenging, few regulatory elements have been characterized in vivo. Here, we use inducible Cas9 and multiplexed guide RNAs to create hundreds of mutations in enhancers/promoters and 3' UTRs of 16 genes in C. elegans. Our software crispr-DART analyzes indel mutations in targeted DNA sequencing. We quantify the impact of mutations on expression and fitness by targeted RNA sequencing and DNA sampling. When applying our approach to the lin-41 3' UTR, generating hundreds of mutants, we find that the two adjacent binding sites for the miRNA let-7 can regulate lin-41 expression independently of each other. Finally, we map regulatory genotypes to phenotypic traits for several genes. Our approach enables parallel analysis of regulatory sequences directly in animals.

Identifying and validating the presence of Guanine-Quadruplexes (G4) within the blood fluke parasite Schistosoma mansoni

Schistosomiasis is a neglected tropical disease that currently affects over 250 million individuals worldwide. In the absence of an immunoprophylactic vaccine and the recognition that mono-chemotherapeutic control of schistosomiasis by praziquantel has limitations, new strategies for managing disease burden are urgently needed. A better understanding of schistosome biology could identify previously undocumented areas suitable for the development of novel interventions. Here, for the first time, we detail the presence of G-quadruplexes (G4) and putative quadruplex forming sequences (PQS) within the Schistosoma mansoni genome. We find that G4 are present in both intragenic and intergenic regions of the seven autosomes as well as the sex-defining allosome pair. Amongst intragenic regions, G4 are particularly enriched in 3´ UTR regions. Gene Ontology (GO) term analysis evidenced significant G4 enrichment in the wnt signalling pathway (p<0.05) and PQS oligonucleotides synthetically derived from wnt-related genes resolve into parallel and anti-parallel G4 motifs as elucidated by circular dichroism (CD) spectroscopy. Finally, utilising a single chain anti-G4 antibody called BG4, we confirm the in situ presence of G4 within both adult female and male worm nuclei. These results collectively suggest that G4-targeted compounds could be tested as novel anthelmintic agents and highlights the possibility that G4-stabilizing molecules could be progressed as candidates for the treatment of schistosomiasis.