Spatial Transcriptomics of Nematodes Identifies Sperm Cells as a Source of Genomic Novelty and Rapid Evolution

Single-cell phylotranscriptomics of developmental and cell type evolution

Single-cell phylotranscriptomics is an emerging tool to reveal the molecular and cellular mechanisms of evolution. We summarize its utility in studying the hourglass pattern of ontogenetic evolution and for understanding the evolutionary history of cell types. The developmental hourglass model suggests that the mid-embryonic stage is the most conserved period of development across species, which is supported by morphological and molecular studies. Single-cell phylotranscriptomic analysis has revealed previously underappreciated heterogeneity in transcriptome ages among lineages and cell types throughout development, and has identified the lineages and tissues that drive the whole-organism hourglass pattern. Single-cell transcriptome age analyses also provide important insights into the origin of germ layers, the different selective forces on tissues during adaptation, and the evolutionary relationships between cell types.

Origin of functional de novo genes in humans from "hopeful monsters"

For a long time, it was believed that new genes arise only from modifications of preexisting genes, but the discovery of de novo protein-coding genes that originated from noncoding DNA regions demonstrates the existence of a "motherless" origination process for new genes. However, the features, distributions, expression profiles, and origin modes of these genes in humans seem to support the notion that their origin is not a purely "motherless" process; rather, these genes arise preferentially from genomic regions encoding preexisting precursors with gene-like features. In such a case, the gene loci are typically not brand new. In this short review, we will summarize the definition and features of human de novo genes and clarify their process of origination from ancestral non-coding genomic regions. In addition, we define the favored precursors, or "hopeful monsters," for the origin of de novo genes and present a discussion of the functional significance of these young genes in brain development and tumorigenesis in humans. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

Young duplicate genes show developmental stage- and cell type-specific expression and function in Caenorhabditis elegans

Gene duplication produces the material that fuels evolutionary innovation. The "out-of-testis" hypothesis suggests that sperm competition creates selective pressure encouraging the emergence of new genes in male germline, but the somatic expression and function of the newly evolved genes are not well understood. We systematically mapped the expression of young duplicate genes throughout development in Caenorhabditis elegans using both whole-organism and single-cell transcriptomic data. Based on the expression dynamics across developmental stages, young duplicate genes fall into three clusters that are preferentially expressed in early embryos, mid-stage embryos, and late-stage larvae. Early embryonic genes are involved in protein degradation and develop essentiality comparable to the genomic average. In mid-to-late embryos and L4-stage larvae, young genes are enriched in intestine, epidermal cells, coelomocytes, and amphid chemosensory neurons. Their molecular functions and inducible expression indicate potential roles in innate immune response and chemosensory perceptions, which may contribute to adaptation outside of the sperm.

Comparative Genomics of Sex, Chromosomes, and Sex Chromosomes in Caenorhabditis elegans and Other Nematodes

The nematode phylum has evolved a remarkable diversity of reproductive modes, including the repeated emergence of asexuality and hermaphroditism across divergent clades. The species-richness and small genome size of nematodes make them ideal systems for investigating the genome-wide causes and consequences of such major transitions. The availability of functional annotations for most Caenorhabditis elegans genes further allows the linking of patterns of gene content evolution with biological processes. Such gene-centric studies were recently complemented by investigations of chromosome evolution that made use of the first chromosome-scale genome assemblies outside the Caenorhabditis genus. This review highlights recent comparative genomic studies of reproductive mode evolution addressing the hybrid origin of asexuality and the parallel gene loss following the emergence of hermaphroditism. It further summarizes ongoing efforts to characterize ancient linkage blocks called Nigon elements, which form central units of chromosome evolution. Fusions between Nigon elements have been demonstrated to impact recombination and speciation. Finally, multiple recent fusions between autosomal and the sex-linked Nigon element reveal insights into the dynamic evolution of sex chromosomes across various timescales.

Co-option of an Astacin Metalloprotease Is Associated with an Evolutionarily Novel Feeding Morphology in a Predatory Nematode

Animals consume a wide variety of food sources to adapt to different environments. However, the genetic mechanisms underlying the acquisition of evolutionarily novel feeding morphology remain largely unknown. While the nematode Caenorhabditis elegans feeds on bacteria, the satellite species Pristionchus pacificus exhibits predatory feeding behavior toward other nematodes, which is an evolutionarily novel feeding habit. Here, we found that the astacin metalloprotease Ppa-NAS-6 is required for the predatory killing by P. pacificus. Ppa-nas-6 mutants were defective in predation-associated characteristics, specifically the tooth morphogenesis and tooth movement during predation. Comparison of expression patterns and rescue experiments of nas-6 in P. pacificus and C. elegans suggested that alteration of the spatial expression patterns of NAS-6 may be vital for acquiring predation-related traits. Reporter analysis of the Ppa-nas-6 promoter in C. elegans revealed that the alteration in expression patterns was caused by evolutionary changes in cis- and trans-regulatory elements. This study suggests that the co-option of a metalloprotease is involved in an evolutionarily novel feeding morphology.

Miniature spatial transcriptomics for studying parasite-endosymbiont relationships at the micro scale

Several important human infectious diseases are caused by microscale-sized parasitic nematodes like filarial worms. Filarial worms have their own spatial tissue organization; to uncover this tissue structure, we need methods that can spatially resolve these miniature specimens. Most filarial worms evolved a mutualistic association with endosymbiotic bacteria Wolbachia. However, the mechanisms underlying the dependency of filarial worms on the fitness of these bacteria remain unknown. As Wolbachia is essential for the development, reproduction, and survival of filarial worms, we spatially explored how Wolbachia interacts with the worm's reproductive system by performing a spatial characterization using Spatial Transcriptomics (ST) across a posterior region containing reproductive tissue and developing embryos of adult female Brugia malayi worms. We provide a proof-of-concept for miniature-ST to explore spatial gene expression patterns in small sample types, demonstrating the method's ability to uncover nuanced tissue region expression patterns, observe the spatial localization of key B. malayi - Wolbachia pathway genes, and co-localize the B. malayi spatial transcriptome in Wolbachia tissue regions, also under antibiotic treatment. We envision our approach will open up new avenues for the study of infectious diseases caused by micro-scale parasitic worms.

A Mediator subunit imparts robustness to a polyphenism decision

Polyphenism is a type of developmental plasticity that translates continuous environmental variability into discontinuous phenotypes. Such discontinuity likely requires a switch between alternative gene-regulatory networks, a principle that has been borne out by mechanisms found to promote morph-specific gene expression. However, whether robustness is required to execute a polyphenism decision has awaited testing at the molecular level. Here, we used a nematode model for polyphenism, Pristionchus pacificus, to identify the molecular regulatory factors that ensure the development of alternative forms. This species has a dimorphism in its adult feeding structures, specifically teeth, which are a morphological novelty that allows predation on other nematodes. Through a forward genetic screen, we determined that a duplicate homolog of the Mediator subunit MDT-15/MED15, P. pacificus MDT-15.1, is necessary for the polyphenism and the robustness of the resulting phenotypes. This transcriptional coregulator, which has a conserved role in metabolic responses to nutritional stress, coordinates these processes with its effects on this diet-induced polyphenism. Moreover, this MED15 homolog genetically interacts with two nuclear receptors, NHR-1 and NHR-40, to achieve dimorphism: Single and double mutants for these three factors result in morphologies that together produce a continuum of forms between the extremes of the polyphenism. In summary, we have identified a molecular regulator that confers discontinuity to a morphological polyphenism, while also identifying a role for MED15 as a plasticity effector.

An update on the role and potential mechanisms of clock genes regulating spermatogenesis: A systematic review of human and animal experimental studies

Circadian clocks can be traced in nearly all life kingdoms, with the male reproductive system no exception. However, our understanding of the circadian clock in spermatogenesis seems to fall behind other scenarios. The present review aims to summarize the current knowledge about the role and especially the potential mechanisms of clock genes in spermatogenesis regulation. Accumulating studies have revealed rhythmic oscillation in semen parameters and some physiological events of spermatogenesis. Disturbing the clock gene expression by genetic mutations or environmental changes will also notably damage spermatogenesis. On the other hand, the mechanisms of spermatogenetic regulation by clock genes remain largely unclear. Some recent studies, although not revealing the entire mechanisms, indeed attempted to shed light on this issue. Emerging clues hinted that gonadal hormones, retinoic acid signaling, homologous recombination, and the chromatoid body might be involved in the regulation of spermatogenesis by clock genes. Then we highlight the challenges and the promising directions for future studies so as to stimulate attention to this critical field which has not gained adequate concern.

Caenorhabditis elegans in microgravity: An omics perspective

The application of omics to study Caenorhabditis elegans (C. elegans) in the context of spaceflight is increasing, illuminating the wide-ranging biological impacts of spaceflight on physiology. In this review, we highlight the application of omics, including transcriptomics, genomics, proteomics, multi-omics, and integrated omics in the study of spaceflown C. elegans, and discuss the impact, use, and future direction of this branch of research. We highlight the variety of molecular alterations that occur in response to spaceflight, most notably changes in metabolic and neuromuscular gene regulation. These transcriptional features are reproducible and evident across many spaceflown species (e.g., mice and astronauts), supporting the use of C. elegans as a model organism to study spaceflight physiology with translational capital. Integrating tissue-specific, spatial, and multi-omics approaches, which quantitatively link molecular responses to phenotypic adaptations, will facilitate the identification of candidate regulatory molecules for therapeutic intervention and thus represents the next frontiers in C. elegans space omics research.

Thousands of Pristionchus pacificus orphan genes were integrated into developmental networks that respond to diverse environmental microbiota

Adaptation of organisms to environmental change may be facilitated by the creation of new genes. New genes without homologs in other lineages are known as taxonomically-restricted orphan genes and may result from divergence or de novo formation. Previously, we have extensively characterized the evolution and origin of such orphan genes in the nematode model organism Pristionchus pacificus. Here, we employ large-scale transcriptomics to establish potential functional associations and to measure the degree of transcriptional plasticity among orphan genes. Specifically, we analyzed 24 RNA-seq samples from adult P. pacificus worms raised on 24 different monoxenic bacterial cultures. Based on coexpression analysis, we identified 28 large modules that harbor 3,727 diplogastrid-specific orphan genes and that respond dynamically to different bacteria. These coexpression modules have distinct regulatory architecture and also exhibit differential expression patterns across development suggesting a link between bacterial response networks and development. Phylostratigraphy revealed a considerably high number of family- and even species-specific orphan genes in certain coexpression modules. This suggests that new genes are not attached randomly to existing cellular networks and that integration can happen very fast. Integrative analysis of protein domains, gene expression and ortholog data facilitated the assignments of biological labels for 22 coexpression modules with one of the largest, fast-evolving module being associated with spermatogenesis. In summary, this work presents the first functional annotation for thousands of P. pacificus orphan genes and reveals insights into their integration into environmentally responsive gene networks.

Evolution and implications of de novo genes in humans

Genes and translated open reading frames (ORFs) that emerged de novo from previously non-coding sequences provide species with opportunities for adaptation. When aberrantly activated, some human-specific de novo genes and ORFs have disease-promoting properties-for instance, driving tumour growth. Thousands of putative de novo coding sequences have been described in humans, but we still do not know what fraction of those ORFs has readily acquired a function. Here, we discuss the challenges and controversies surrounding the detection, mechanisms of origin, annotation, validation and characterization of de novo genes and ORFs. Through manual curation of literature and databases, we provide a thorough table with most de novo genes reported for humans to date. We re-evaluate each locus by tracing the enabling mutations and list proposed disease associations, protein characteristics and supporting evidence for translation and protein detection. This work will support future explorations of de novo genes and ORFs in humans.

The nematode Oscheius tipulae as a genetic model for programmed DNA elimination

Programmed DNA elimination (PDE) is a notable exception to the paradigm of genome integrity. In metazoa, PDE often occurs coincident with germline to somatic cell differentiation. During PDE, portions of genomic DNA are lost, resulting in reduced somatic genomes. Prior studies have described the sequences lost, as well as chromosome behavior, during metazoan PDE. However, a system for studying the mechanisms and consequences of PDE in metazoa is lacking. Here, we present a functional and genetic model for PDE in the free-living Rhabditidae nematode Oscheius tipulae, a family that also includes Caenorhabditis elegans. O. tipulae was recently suggested to eliminate DNA. Using staged embryos and DNA FISH, we showed that O. tipulae PDE occurs during embryogenesis at the 8-16 cell stages. We identified a conserved motif, named Sequence For Elimination (SFE), for all 12 break sites on the six chromosomes at the junctions of retained and eliminated DNA. SFE mutants exhibited a "fail-to-eliminate" phenotype only at the modified sites. END-seq revealed that breaks can occur at multiple positions within the SFE, with extensive end resection followed by telomere addition to both retained and eliminated ends. We identified many functional SFEs at the chromosome ends through END-seq in the wild-type embryos, genome sequencing of SFE mutants, and comparative genomics of 23 wild isolates. We suggest that these alternative SFEs provide flexibility in the sequences eliminated and a fail-safe mechanism for PDE. These studies establish O. tipulae as a new, attractive model for studying the mechanisms and consequences of PDE in a metazoan.

RNA-Sequencing of Heterorhabditis nematodes to identify factors involved in symbiosis with Photorhabdus bacteria

Background

Nematodes are a major group of soil inhabiting organisms. Heterorhabditis nematodes are insect-pathogenic nematodes and live in a close symbiotic association with Photorhabdus bacteria. Heterorhabditis-Photorhabdus pair offers a powerful and genetically tractable model to study animal-microbe symbiosis. It is possible to generate symbiont bacteria free (axenic) stages in Heterorhabditis. Here, we compared the transcriptome of symbiotic early-adult stage Heterorhabditis nematodes with axenic early-adult nematodes to determine the nematode genes and pathways involved in symbiosis with Photorhabdus bacteria.

Results

A de-novo reference transcriptome assembly of 95.7 Mb was created for H. bacteriophora by using all the reads. The assembly contained 46,599 transcripts with N50 value of 2,681 bp and the average transcript length was 2,054 bp. The differentially expressed transcripts were identified by mapping reads from symbiotic and axenic nematodes to the reference assembly. A total of 754 differentially expressed transcripts were identified in symbiotic nematodes as compared to the axenic nematodes. The ribosomal pathway was identified as the most affected among the differentially expressed transcripts. Additionally, 12,151 transcripts were unique to symbiotic nematodes. Endocytosis, cAMP signalling and focal adhesion were the top three enriched pathways in symbiotic nematodes, while a large number of transcripts coding for various responses against bacteria, such as bacterial recognition, canonical immune signalling pathways, and antimicrobial effectors could also be identified.

Conclusions

The symbiotic Heterorhabditis nematodes respond to the presence of symbiotic bacteria by expressing various transcripts involved in a multi-layered immune response which might represent non-systemic and evolved localized responses to maintain mutualistic bacteria at non-threatening levels. Subject to further functional validation of the identified transcripts, our findings suggest that Heterorhabditis nematode immune system plays a critical role in maintenance of symbiosis with Photorhabdus bacteria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08952-4.

NHR-23 and SPE-44 regulate distinct sets of genes during Caenorhabditis elegans spermatogenesis

Spermatogenesis is the process through which mature male gametes are formed and is necessary for the transmission of genetic information. While much work has established how sperm fate is promoted and maintained, less is known about how the sperm morphogenesis program is executed. We previously identified a novel role for the nuclear hormone receptor transcription factor, NHR-23, in promoting Caenorhabditis elegans spermatogenesis. The depletion of NHR-23 along with SPE-44, another transcription factor that promotes spermatogenesis, caused additive phenotypes. Through RNA-seq, we determined that NHR-23 and SPE-44 regulate distinct sets of genes. The depletion of both NHR-23 and SPE-44 produced yet another set of differentially regulated genes. NHR-23-regulated genes are enriched in phosphatases, consistent with the switch from genome quiescence to post-translational regulation in spermatids. In the parasitic nematode Ascaris suum, MFP1 and MFP2 control the polymerization of Major Sperm Protein, the molecule that drives sperm motility and serves as a signal to promote ovulation. NHR-23 and SPE-44 regulate several MFP2 paralogs, and NHR-23 depletion from the male germline caused defective localization of MSD/MFP1 and NSPH-2/MFP2. Although NHR-23 and SPE-44 do not transcriptionally regulate the casein kinase gene spe-6, a key regulator of sperm development, SPE-6 protein is lost following NHR-23+SPE-44 depletion. Together, these experiments provide the first mechanistic insight into how NHR-23 promotes spermatogenesis and an entry point to understanding the synthetic genetic interaction between nhr-23 and spe-44.

Multiple Pristionchus pacificus genomes reveal distinct evolutionary dynamics between de novo candidates and duplicated genes

The birth of new genes is a major molecular innovation driving phenotypic diversity across all domains of life. Although repurposing of existing protein-coding material by duplication is considered the main process of new gene formation, recent studies have discovered thousands of transcriptionally active sequences as a rich source of new genes. However, differential loss rates have to be assumed to reconcile the high birth rates of these incipient de novo genes with the dominance of ancient gene families in individual genomes. Here, we test this rapid turnover hypothesis in the context of the nematode model organism Pristionchus pacificus We extended the existing species-level phylogenomic framework by sequencing the genomes of six divergent P. pacificus strains. We used these data to study the evolutionary dynamics of different age classes and categories of origin at a population level. Contrasting de novo candidates with new families that arose by duplication and divergence from known genes, we find that de novo candidates are typically shorter, show less expression, and are overrepresented on the sex chromosome. Although the contribution of de novo candidates increases toward young age classes, multiple comparisons within the same age class showed significantly higher attrition in de novo candidates than in known genes. Similarly, young genes remain under weak evolutionary constraints with de novo candidates representing the fastest evolving subcategory. Altogether, this study provides empirical evidence for the rapid turnover hypothesis and highlights the importance of the evolutionary timescale when quantifying the contribution of different mechanisms toward new gene formation.

Laser Microdissection for Species-Agnostic Single-Tissue Applications

Single-cell methodologies have revolutionized the analysis of the transcriptomes of specific cell types. However, they often require species-specific genetic "toolkits," such as promoters driving tissue-specific expression of fluorescent proteins. Further, protocols that disrupt tissues to isolate individual cells remove cells from their native environment (e.g., signaling from neighbors) and may result in stress responses or other differences from native gene expression states. In the present protocol, laser microdissection (LMD) is optimized to isolate individual nematode tail tips for the study of gene expression during male tail tip morphogenesis. LMD allows the isolation of a portion of the animal without the need for cellular disruption or species-specific toolkits and is thus applicable to any species. Subsequently, single-cell RNA-seq library preparation protocols such as CEL-Seq2 can be applied to LMD-isolated single tissues and analyzed using standard pipelines, given that a well-annotated genome or transcriptome is available for the species. Such data can be used to establish how conserved or different the transcriptomes are that underlie the development of that tissue in different species. Limitations include the ability to cut out the tissue of interest and the sample size. A power analysis shows that as few as 70 tail tips per condition are required for 80% power. Tight synchronization of development is needed to obtain this number of animals at the same developmental stage. Thus, a method to synchronize animals at 1 h intervals is also described.

Small RNA pathways in the nematode Ascaris in the absence of piRNAs

Small RNA pathways play key and diverse regulatory roles in C. elegans, but our understanding of their conservation and contributions in other nematodes is limited. We analyzed small RNA pathways in the divergent parasitic nematode Ascaris. Ascaris has ten Argonautes with five worm-specific Argonautes (WAGOs) that associate with secondary 5’-triphosphate 22-24G-RNAs. These small RNAs target repetitive sequences or mature mRNAs and are similar to the C. elegans mutator, nuclear, and CSR-1 small RNA pathways. Even in the absence of a piRNA pathway, Ascaris CSR-1 may still function to “license” as well as fine-tune or repress gene expression. Ascaris ALG-4 and its associated 26G-RNAs target and likely repress specific mRNAs during testis meiosis. Ascaris WAGO small RNAs demonstrate target plasticity changing their targets between repeats and mRNAs during development. We provide a unique and comprehensive view of mRNA and small RNA expression throughout spermatogenesis. Overall, our study illustrates the conservation, divergence, dynamics, and flexibility of small RNA pathways in nematodes.

The parasitic nematode Ascaris lacks piRNAs. Here the authors compare Argonaute proteins and small RNAs from C. elegans and Ascaris, expanding our understanding of the conservation, divergence, and flexibility of Argonautes and small RNA pathways in nematodes.

Horizontally Acquired Cellulases Assist the Expansion of Dietary Range in Pristionchus Nematodes

Horizontal gene transfer (HGT) enables the acquisition of novel traits via non-Mendelian inheritance of genetic material. HGT plays a prominent role in the evolution of prokaryotes, whereas in animals, HGT is rare and its functional significance is often uncertain. Here, we investigate horizontally acquired cellulase genes in the free-living nematode model organism Pristionchus pacificus. We show that these cellulase genes 1) are likely of eukaryotic origin, 2) are expressed, 3) have protein products that are secreted and functional, and 4) result in endo-cellulase activity. Using CRISPR/Cas9, we generated an octuple cellulase mutant, which lacks all eight cellulase genes and cellulase activity altogether. Nonetheless, this cellulase-null mutant is viable and therefore allows a detailed analysis of a gene family that was horizontally acquired. We show that the octuple cellulase mutant has associated fitness costs with reduced fecundity and slower developmental speed. Furthermore, by using various Escherichia coli K-12 strains as a model for cellulosic biofilms, we demonstrate that cellulases facilitate the procurement of nutrients from bacterial biofilms. Together, our analysis of cellulases in Pristionchus provides comprehensive evidence from biochemistry, genetics, and phylogeny, which supports the integration of horizontally acquired genes into the complex life history strategy of this soil nematode.

Analysis and Visualization of Spatial Transcriptomic Data

Human and animal tissues consist of heterogeneous cell types that organize and interact in highly structured manners. Bulk and single-cell sequencing technologies remove cells from their original microenvironments, resulting in a loss of spatial information. Spatial transcriptomics is a recent technological innovation that measures transcriptomic information while preserving spatial information. Spatial transcriptomic data can be generated in several ways. RNA molecules are measured by in situ sequencing, in situ hybridization, or spatial barcoding to recover original spatial coordinates. The inclusion of spatial information expands the range of possibilities for analysis and visualization, and spurred the development of numerous novel methods. In this review, we summarize the core concepts of spatial genomics technology and provide a comprehensive review of current analysis and visualization methods for spatial transcriptomics.

Unconventional conservation reveals structure-function relationships in the synaptonemal complex

Functional requirements constrain protein evolution, commonly manifesting in a conserved amino acid sequence. Here, we extend this idea to secondary structural features by tracking their conservation in essential meiotic proteins with highly diverged sequences. The synaptonemal complex (SC) is a ~100-nm-wide ladder-like meiotic structure present in all eukaryotic clades, where it aligns parental chromosomes and regulates exchanges between them. Despite the conserved ultrastructure and functions of the SC, SC proteins are highly divergent within Caenorhabditis. However, SC proteins have highly conserved length and coiled-coil domain structure. We found the same unconventional conservation signature in Drosophila and mammals, and used it to identify a novel SC protein in Pristionchus pacificus, Ppa-SYP-1. Our work suggests that coiled-coils play wide-ranging roles in the structure and function of the SC, and more broadly, that expanding sequence analysis beyond measures of per-site similarity can enhance our understanding of protein evolution and function.

Programmed DNA elimination: silencing genes and repetitive sequences in somatic cells

In a multicellular organism, the genomes of all cells are in general the same. Programmed DNA elimination is a notable exception to this genome constancy rule. DNA elimination removes genes and repetitive elements in the germline genome to form a reduced somatic genome in various organisms. The process of DNA elimination within an organism is highly accurate and reproducible; it typically occurs during early embryogenesis, coincident with germline-soma differentiation. DNA elimination provides a mechanism to silence selected genes and repeats in somatic cells. Recent studies in nematodes suggest that DNA elimination removes all chromosome ends, resolves sex chromosome fusions, and may also promote the birth of novel genes. Programmed DNA elimination processes are diverse among species, suggesting DNA elimination likely has evolved multiple times in different taxa. The growing list of organisms that undergo DNA elimination indicates that DNA elimination may be more widespread than previously appreciated. These various organisms will serve as complementary and comparative models to study the function, mechanism, and evolution of programmed DNA elimination in metazoans.

New Genes Interacted with Recent Whole Genome Duplicates in the Fast Stem Growth of Bamboos

As drivers of evolutionary innovations, new genes allow organisms to explore new niches. However, clear examples of this process remain scarce. Bamboos, the unique grass lineage diversifying into the forest, have evolved with a key innovation of fast growth of woody stem, reaching up to 1 m/day. Here, we identify 1,622 bamboo-specific orphan genes that appeared in recent 46 million years, and 19 of them evolved from noncoding ancestral sequences with entire de novo origination process reconstructed. The new genes evolved gradually in exon−intron structure, protein length, expression specificity, and evolutionary constraint. These new genes, whether or not from de novo origination, are dominantly expressed in the rapidly developing shoots, and make transcriptomes of shoots the youngest among various bamboo tissues, rather than reproductive tissue in other plants. Additionally, the particularity of bamboo shoots has also been shaped by recent whole-genome duplicates (WGDs), which evolved divergent expression patterns from ancestral states. New genes and WGDs have been evolutionarily recruited into coexpression networks to underline fast-growing trait of bamboo shoot. Our study highlights the importance of interactions between new genes and genome duplicates in generating morphological innovation.

Single worm transcriptomics identifies a developmental core network of oscillating genes with deep conservation across nematodes

High-resolution spatial and temporal maps of gene expression have facilitated a comprehensive understanding of animal development and evolution. In nematodes, the small body size represented a major challenge for such studies, but recent advancements have helped overcome this limitation. Here, we have implemented single worm transcriptomics (SWT) in the nematode model organism Pristionchus pacificus to provide a high-resolution map of the developmental transcriptome. We selected 38 time points from hatching of the J2 larvae to young adults to perform transcriptome analysis over 60 h of postembryonic development. A mean sequencing depth of 4.5 million read pairs allowed the detection of more than 23,135 (80%) of all genes. Nearly 3000 (10%) genes showed oscillatory expression with discrete expression levels, phases, and amplitudes. Gene age analysis revealed an overrepresentation of ancient gene classes among oscillating genes, and around one-third of them have 1:1 orthologs in C. elegans. One important gene family overrepresented among oscillating genes is collagens. Several of these collagen genes are regulated by the developmental switch gene eud-1, indicating a potential function in the regulation of mouth-form plasticity, a key developmental process in this facultative predatory nematode. Together, our analysis provides (1) an updated protocol for SWT in nematodes that is applicable to many microscopic species, (2) a 1- to 2-h high-resolution catalog of P. pacificus gene expression throughout postembryonic development, and (3) a comparative analysis of oscillatory gene expression between the two model organisms P. pacificus and C. elegans and associated evolutionary dynamics.

MFP1/MSD-1 and MFP2/NSPH-2 co-localize with MSP during C. elegans spermatogenesis

Until recently, the only verified component of Fibrous Bodies (FBs) within Caenorhabditis elegans spermatocytes was the Major Sperm Protein (MSP), a nematode-specific cytoskeletal element. Earlier studies in the pig parasite Ascaris suum had identified accessory proteins that facilitate MSP polymerization and depolymerization within the pseudopod of crawling spermatozoa. In this study, we show that C. elegans homologs of the two Ascaris accessory proteins MFP1 and MFP2 co-localize with MSP in both the pseudopods of C. elegans sperm and the FBs of C. elegans spermatocytes.

Spatial transcriptomics of the nematode Caenorhabditis elegans using RNA tomography

RNA tomography or tomo-seq combines mRNA sequencing and cryo-sectioning to spatially resolve gene expression. We have adapted this method for the nematode Caenorhabditis elegans to generate anteroposterior gene expression maps at near-cellular resolution. Here, we provide a detailed overview of the method and present two approaches: one that includes RNA isolation for maximum sensitivity and one that is suitable for partial automatization and is therefore less time-consuming. For complete details on the use and execution of this protocol, please refer to Ebbing et al. (2018).

Genomics of the Parasitic Nematode Ascaris and Its Relatives

Nematodes of the genus Ascaris are important parasites of humans and swine, and the phylogenetically related genera (Parascaris, Toxocara, and Baylisascaris) infect mammals of veterinary interest. Over the last decade, considerable genomic resources have been established for Ascaris, including complete germline and somatic genomes, comprehensive mRNA and small RNA transcriptomes, as well as genome-wide histone and chromatin data. These datasets provide a major resource for studies on the basic biology of these parasites and the host-parasite relationship. Ascaris and its relatives undergo programmed DNA elimination, a highly regulated process where chromosomes are fragmented and portions of the genome are lost in embryonic cells destined to adopt a somatic fate, whereas the genome remains intact in germ cells. Unlike many model organisms, Ascaris transcription drives early development beginning prior to pronuclear fusion. Studies on Ascaris demonstrated a complex small RNA network even in the absence of a piRNA pathway. Comparative genomics of these ascarids has provided perspectives on nematode sex chromosome evolution, programmed DNA elimination, and host-parasite coevolution. The genomic resources enable comparison of proteins across diverse species, revealing many new potential drug targets that could be used to control these parasitic nematodes.

The community-curated Pristionchus pacificus genome facilitates automated gene annotation improvement in related nematodes

Background

The nematode Pristionchus pacificus is an established model organism for comparative studies with Caenorhabditis elegans. Over the past years, it developed into an independent animal model organism for elucidating the genetic basis of phenotypic plasticity. Community-based curations were employed recently to improve the quality of gene annotations of P. pacificus and to more easily facilitate reverse genetic studies using candidate genes from C. elegans.

Results

Here, I demonstrate that the reannotation of phylogenomic data from nine related nematode species using the community-curated P. pacificus gene set as homology data substantially improves the quality of gene annotations. Benchmarking of universal single copy orthologs (BUSCO) estimates a median completeness of 84% which corresponds to a 9% increase over previous annotations. Nevertheless, the ability to infer gene models based on homology already drops beyond the genus level reflecting the rapid evolution of nematode lineages. This also indicates that the highly curated C. elegans genome is not optimally suited for annotating non-Caenorhabditis genomes based on homology. Furthermore, comparative genomic analysis of apparently missing BUSCO genes indicates a failure of ortholog detection by the BUSCO pipeline due to the insufficient sample size and phylogenetic breadth of the underlying OrthoDB data set. As a consequence, the quality of multiple divergent nematode genomes might be underestimated.

Conclusions

This study highlights the need for optimizing gene annotation protocols and it demonstrates the benefit of a high quality genome for phylogenomic data of related species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07529-x.