The genome of a subterrestrial nematode reveals adaptations to heat

Temperature Effects on Expression Levels of hsp Genes in Eggs and Second-Stage Juveniles of Meloidogyne hapla Chitwood, 1949

Meloidogyne hapla is one of the most important nematode pathogens. It is a sedentary, biotrophic parasite of plants that overwinters in the soil or in diseased roots. The development of M. hapla is temperature dependent. Numerous studies have been performed on the effect of temperature on the development of M. hapla, but only a few of them analyzed the heat shock protein (hsp) genes. The aim of the study was to perform expression profiling of eight hsp genes (Mh-hsp90, Mh-hsp1, Mh-hsp4, Mh-hsp6, Mh-hsp60, Mh-dnj19, Mh-hsp43, and Mh-hsp12.2) at two development stages of M. hapla, i.e., in eggs and second-stage juveniles (J2). The eggs and J2 were incubated under cold stress (5 °C), heat stress (35 °C, 40 °C), and non-stress (10 °C, 20 °C, and 30 °C) conditions. Expression profiling was performed by qPCR. It was demonstrated that only two genes, Mh-hsp60 and Mh-dnj19, have been upregulated by heat and cold stress at both development stages. Heat stress upregulated the expression of more hsp genes than cold stress did. The level of upregulation of most hsp genes was more marked in J2 than in eggs. The obtained results suggest that the Mh-hsp90 and Mh-hsp1 genes can be used as bioindicators of environmental impacts on nematodes of the Meloidogyne genus.

Multidrug-resistant plasmid RP4 inhibits the nitrogen removal capacity of ammonia-oxidizing archaea, ammonia-oxidizing bacteria, and comammox in activated sludge

In wastewater treatment plants (WWTPs), ammonia oxidation is primarily carried out by three types of ammonia oxidation microorganisms (AOMs): ammonia-oxidizing archaea (AOA), ammonia-oxidizing bacteria (AOB), and comammox (CMX). Antibiotic resistance genes (ARGs), which pose an important public health concern, have been identified at every stage of wastewater treatment. However, few studies have focused on the impact of ARGs on ammonia removal performance. Therefore, our study sought to investigate the effect of the representative multidrug-resistant plasmid RP4 on the functional microorganisms involved in ammonia oxidation. Using an inhibitor-based method, we first evaluated the contributions of AOA, AOB, and CMX to ammonia oxidation in activated sludge, which were determined to be 13.7%, 41.1%, and 39.1%, respectively. The inhibitory effects of C2H2, C8H14, and 3,4-dimethylpyrazole phosphate (DMPP) were then validated by qPCR. After adding donor strains to the sludge, fluorescence in situ hybridization (FISH) imaging analysis demonstrated the co-localization of RP4 plasmids and all three AOMs, thus confirming the horizontal gene transfer (HGT) of the RP4 plasmid among these microorganisms. Significant inhibitory effects of the RP4 plasmid on the ammonia nitrogen consumption of AOA, AOB, and CMX were also observed, with inhibition rates of 39.7%, 36.2%, and 49.7%, respectively. Moreover, amoA expression in AOB and CMX was variably inhibited by the RP4 plasmid, whereas AOA amoA expression was not inhibited. These results demonstrate the adverse environmental effects of the RP4 plasmid and provide indirect evidence supporting plasmid-mediated conjugation transfer from bacteria to archaea.

Evolution of a biological thermocouple by adaptation of cytochrome c oxidase in a subterrestrial metazoan

In this study we report a naturally evolved temperature-sensing electrical regulator in the cytochrome c oxidase of the Devil Worm, Halicephalobus mephisto. This extremophile metazoan was isolated 1.3 km underground in a South African goldmine, where it adapted to heat and potentially to hypoxia, making its mitochondrial sequence a likely target of adaptational change. We obtained the full mitochondrial genome sequence of this organism, and show through dN/dS analysis statistically robust evidence of positive selection in H. mephisto cytochrome c oxidase subunits. Seventeen of these positively-selected amino acid substitutions were localized in proximity to the H- and K-pathway proton channels of the complex. Surprisingly, the H. mephisto cytochrome c oxidase proton pump completely shuts down at low temperatures (20°C) leading to approximately a 4.8-fold reduction in the transmembrane proton gradient voltage (ΔΨm) compared to optimal temperature (37°C). Direct measurement of oxygen consumption found a corresponding 4.7-fold drop at 20°C compared to 37°C. Correspondingly, the lifecycle of H. mephisto takes four-fold longer at the low temperature compared to higher. This elegant evolutionary adaptation creates a finely-tuned mitochondrial temperature sensor, allowing this ectothermic organism to maximize its reproductive success in varying environmental temperatures. Our study shows that evolutionary innovation may remodel core metabolism to make it more accurately map onto environmental variation.

Virus-like transposons cross the species barrier and drive the evolution of genetic incompatibilities

Horizontal gene transfer, the movement of genetic material between species, has been reported across all major eukaryotic lineages. However, the underlying mechanisms of transfer and their impact on genome evolution are still poorly understood. While studying the evolutionary origin of a selfish element in the nematode Caenorhabditis briggsae, we discovered that Mavericks, ancient virus-like transposons related to giant viruses and virophages, are one of the long-sought vectors of horizontal gene transfer. We found that Mavericks gained a novel herpesvirus-like fusogen in nematodes, leading to the widespread exchange of cargo genes between extremely divergent species, bypassing sexual and genetic barriers spanning hundreds of millions of years. Our results show how the union between viruses and transposons causes horizontal gene transfer and ultimately genetic incompatibilities in natural populations.

Discovering common pathogenic processes between COVID-19 and HFRS by integrating RNA-seq differential expression analysis with machine learning

Zoonotic virus spillover in human hosts including outbreaks of Hantavirus and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) imposes a serious impact on the quality of life of patients. Recent studies provide a shred of evidence that patients with Hantavirus-caused hemorrhagic fever with renal syndrome (HFRS) are at risk of contracting SARS-CoV-2. Both RNA viruses shared a higher degree of clinical features similarity including dry cough, high fever, shortness of breath, and certain reported cases with multiple organ failure. However, there is currently no validated treatment option to tackle this global concern. This study is attributed to the identification of common genes and perturbed pathways by combining differential expression analysis with bioinformatics and machine learning approaches. Initially, the transcriptomic data of hantavirus-infected peripheral blood mononuclear cells (PBMCs) and SARS-CoV-2 infected PBMCs were analyzed through differential gene expression analysis for identification of common differentially expressed genes (DEGs). The functional annotation by enrichment analysis of common genes demonstrated immune and inflammatory response biological processes enriched by DEGs. The protein-protein interaction (PPI) network of DEGs was then constructed and six genes named RAD51, ALDH1A1, UBA52, CUL3, GADD45B, and CDKN1A were identified as the commonly dysregulated hub genes among HFRS and COVID-19. Later, the classification performance of these hub genes were evaluated using Random Forest (RF), Poisson Linear Discriminant Analysis (PLDA), Voom-based Nearest Shrunken Centroids (voomNSC), and Support Vector Machine (SVM) classifiers which demonstrated accuracy >70%, suggesting the biomarker potential of the hub genes. To our knowledge, this is the first study that unveiled biological processes and pathways commonly dysregulated in HFRS and COVID-19, which could be in the next future used for the design of personalized treatment to prevent the linked attacks of COVID-19 and HFRS.

Effect of Vicia sativa L. on Motility, Mortality and Expression Levels of hsp Genes in J2 Stage of Meloidogyne hapla

Assuming that the seeds of Vicia sativa L. have a stressful effect on J2 stage Meloidogyne hapla, we undertook research on the effect of these seeds on the motility and mortality of J2 and determined the expression levels of selected hsp genes in J2. The assessment of the effect of V. sativa seeds on the motility of M. hapla specimens consisted of observing the movement of J2 immersed in a seed diffusate or in a tomato root filtrate at temperatures of 10, 17, and 21°C. In J2 treated with V. sativa (cv. Ina) seed diffusates, the expression level of hsp genes was determined by qPCR. J2 exposed to V. sativa diffusates were found to lose their motility, while their mortality did not exceed 30%. J2 in the seed diffusate were characterized by an increase in the expression levels of the Mh-hsp90, Mh-hsp1, and Mh-hsp43 genes. It is suggested that the hsp90 gene may be a potential bioindicator of the environmental impact on Meloidogyne nematodes. The impaired ability to move in J2 of M. hapla is attributable to the occurrence of V. sativa seeds in their habitat. These studies may contribute to developing methods of reducing crop damage caused by M. hapla.

Nematode gene annotation by machine-learning-assisted proteotranscriptomics enables proteome-wide evolutionary analysis

Nematodes encompass more than 24,000 described species, which were discovered in almost every ecological habitat, and make up >80% of metazoan taxonomic diversity in soils. The last common ancestor of nematodes is believed to date back to ∼650-750 million years, generating a large and phylogenetically diverse group to be explored. However, for most species high-quality gene annotations are incomprehensive or missing. Combining short-read RNA sequencing with mass spectrometry-based proteomics and machine-learning quality control in an approach called proteotranscriptomics, we improve gene annotations for nine genome-sequenced nematode species and provide new gene annotations for three additional species without genome assemblies. Emphasizing the sensitivity of our methodology, we provide evidence for two hitherto undescribed genes in the model organism Caenorhabditis elegans Extensive phylogenetic systems analysis using this comprehensive proteome annotation provides new insights into evolutionary processes of this metazoan group.

Gene family expansions in Antarctic winged midge as a strategy for adaptation to cold environments

Parochlus steinenii is the only flying insect native to Antarctica. To elucidate the molecular mechanisms underlying its adaptation to cold environments, we conducted comparative genomic analyses of P. steinenii and closely related lineages. In an analysis of gene family evolution, 68 rapidly evolving gene families, involved in the innate immune system, unfolded protein response, DNA packaging, protein folding, and unsaturated fatty acid biosynthesis were detected. Some gene families were P. steinenii-specific and showed phylogenetic instability. Acyl-CoA delta desaturase and heat shock cognate protein 70 (Hsc70) were representative gene families, showing signatures of positive selection with multiple gene duplication events. Acyl-CoA delta desaturases may play pivotal roles in membrane fluidity, and expanded Hsc70 genes may function as chaperones or thermal sensors in cold environments. These findings suggest that multiple gene family expansions contributed to the adaptation of P. steinenii to cold environments.

Chromosome-level genome assembly of Bactrocera dorsalis reveals its adaptation and invasion mechanisms

Bactrocera dorsalis is an invasive polyphagous pest causing considerable ecological and economic damage worldwide. We report a high-quality chromosome-level genome assembly and combine various transcriptome data to explore the molecular mechanisms of its rapid adaptation to new environments. The expansions of the DDE transposase superfamily and key gene families related to environmental adaptation and enrichment of the expanded and unique gene families in metabolism and defence response pathways explain its environmental adaptability. The relatively high but not significantly different expression of heat-shock proteins, regardless of the environmental conditions, suggests an intrinsic mechanism underlying its adaptation to high temperatures. The mitogen-activated protein kinase pathway plays a key role in adaptation to new environments. The prevalence of duplicated genes in its genome explains the diversity in the B. dorsalis complex. These findings provide insights into the genetic basis of the invasiveness and diversity of B. dorsalis, explaining its rapid adaptation and expansion.

Jiang et al. sequence the genome of Bactrocera dorsalis, a destructive and invasive agricultural pest. Insights from this chromosome-level assembly shed light on molecular adaptations that allow for the global invasion and expansion of this pest.

Tokorhabditis n. gen. (Rhabditida, Rhabditidae), a comparative nematode model for extremophilic living

Life in extreme environments is typically studied as a physiological problem, although the existence of extremophilic animals suggests that developmental and behavioral traits might also be adaptive in such environments. Here, we describe a new species of nematode, Tokorhabditis tufae, n. gen., n. sp., which was discovered from the alkaline, hypersaline, and arsenic-rich locale of Mono Lake, California. The new species, which offers a tractable model for studying animal-specific adaptations to extremophilic life, shows a combination of unusual reproductive and developmental traits. Like the recently described sister group Auanema, the species has a trioecious mating system comprising males, females, and self-fertilizing hermaphrodites. Our description of the new genus thus reveals that the origin of this uncommon reproductive mode is even more ancient than previously assumed, and it presents a new comparator for the study of mating-system transitions. However, unlike Auanema and almost all other known rhabditid nematodes, the new species is obligately live-bearing, with embryos that grow in utero, suggesting maternal provisioning during development. Finally, our isolation of two additional, molecularly distinct strains of the new genus-specifically from non-extreme locales-establishes a comparative system for the study of extremophilic traits in this model.

Endurance of extremely prolonged nutrient prevention across kingdoms of life

Numerous observations demonstrate that microorganisms can survive very long periods of nutrient deprivation and starvation. Moreover, it is evident that prolonged periods of starvation are a feature of many habitats, and many cells in all kingdoms of life are found in prolonged starvation conditions. Bacteria exhibit a range of responses to long-term starvation. These include genetic adaptations such as the long-term stationary phase and the growth advantage in stationary phase phenotypes characterized by mutations in stress-signaling genes and elevated mutation rates. Here, we suggest using the term "endurance of prolonged nutrient prevention" (EPNP phase), to describe this phase, which was also recently described in eukaryotes. Here, we review this literature and describe the current knowledge about the adaptations to very long-term starvation conditions in bacteria and eukaryotes, its conceptual and structural conservation across all kingdoms of life, and point out possible directions that merit further research.

The GIMAP Family Proteins: An Incomplete Puzzle

Overview: Long-term survival of T lymphocytes in quiescent state is essential to maintain their cell numbers in secondary lymphoid organs and in peripheral circulation. In the BioBreeding diabetes-prone strain of rats (BB-DP), loss of functional GIMAP5 (GTPase of the immune associated nucleotide binding protein 5) results in profound peripheral T lymphopenia. This discovery heralded the identification of a new family of proteins initially called Immune-associated nucleotide binding protein (IAN) family. In this review we will use ‘GIMAP’ to refer to this family of proteins. Recent studies suggest that GIMAP proteins may interact with each other and also be involved in the movement of the cellular cargo along the cytoskeletal network. Here we will summarize the current knowledge on the characteristics and functions of GIMAP family of proteins.

Genome-Wide Transcriptional Responses of Marine Nematode Litoditis marina to Hyposaline and Hypersaline Stresses

Maintenance of osmotic homeostasis is essential for all organisms, especially for marine animals in the ocean with 3% salinity or higher. However, the underlying molecular mechanisms that how marine animals adapt to high salinity environment compared to their terrestrial relatives, remain elusive. Here, we investigated marine animal’s genome-wide transcriptional responses to salinity stresses using an emerging marine nematode model Litoditis marina. We found that the transthyretin-like family genes were significantly increased in both hyposaline and hypersaline conditions, while multiple neurotransmitter receptor and ion transporter genes were down-regulated in both conditions, suggesting the existence of conserved strategies for response to stressful salinity environments in L. marina. Unsaturated fatty acids biosynthesis related genes, neuronal related tubulins and intraflagellar transport genes were specifically up-regulated in hyposaline treated worms. By contrast, cuticle related collagen genes were enriched and up-regulated for hypersaline response. Given a wide range of salinity tolerance of the marine nematodes, this study and further genetic analysis of key gene(s) of osmoregulation in L. marina will likely provide important insights into biological evolution and environmental adaptation mechanisms in nematodes and other invertebrate animals in general.

Direct Detection of DNA and RNA on Carbon Fiber Microelectrodes Using Fast-Scan Cyclic Voltammetry

DNA and RNA have been measured with many techniques but often with relatively long analysis times. In this study, we utilize fast-scan cyclic voltammetry (FSCV) for the subsecond codetection of adenine, guanine, and cytosine, first as free nucleosides, and then within custom synthesized oligos, plasmid DNA, and RNA from the nematode Caenorhabditis elegans. Previous studies have shown the detection of adenosine and guanosine with FSCV with high spatiotemporal resolution, while we have extended the assay to include cytidine and adenine, guanine, and cytosine in RNA and single- and double-stranded DNA (ssDNA and dSDNA). We find that FSCV testing has a higher sensitivity and yields higher peak oxidative currents when detecting shorter oligonucleotides and ssDNA samples at equivalent nucleobase concentrations. This is consistent with an electrostatic repulsion from negatively charged oxide groups on the surface of the carbon fiber microelectrode (CFME), the negative holding potential, and the negatively charged phosphate backbone. Moreover, as opposed to dsDNA, ssDNA nucleobases are not hydrogen-bonded to one another and thus are free to adsorb onto the surface of the carbon electrode. We also demonstrate that the simultaneous determination of nucleobases is not masked even in biologically complex serum samples. This is the first report demonstrating that FSCV, when used with CFMEs, is able to codetect nucleobases when polymerized into DNA or RNA and could potentially pave the way for future uses in clinical, diagnostic, or research applications.

Speciation and adaptive evolution reshape antioxidant enzymatic system diversity across the phylum Nematoda

Background

Nematodes have evolved to survive in diverse ecological niches and can be a serious burden on agricultural economy, veterinary medicine, and public health. Antioxidant enzymes in parasitic nematodes play a critical role in defending against host oxidative stress. However, the features of the evolution of antioxidant enzymes in the phylum Nematoda remain elusive.

Results

Here, we systematically investigated the evolution and gene expression of antioxidant enzymes in the genomes of 59 nematodes and transcriptomes of 20 nematodes. Catalase has been independently lost in several orders, suggesting that it is unnecessary for some nematodes. Unlike in mammals, phospholipid hydroperoxide glutathione peroxidase is widely distributed in nematodes, among which it has evolved independently. We found that superoxide dismutase (SOD) has been present throughout nematode evolutionary process, and the extracellular isoform (SOD3) is diverged from the corresponding enzyme in mammals and has undergone duplication and differentiation in several nematodes. Moreover, the evolution of intracellular and extracellular SOD isoforms in filaria strongly indicates that extracellular SOD3 originated from intracellular SOD1 and underwent rapid evolution to form the diversity of extracellular SOD3. We identify a novel putative metal-independent extracellular SOD presenting independently in Steinernema and Strongyloididae lineage that featured a high expression level in Strongyloides larvae. Sequence divergence of SOD3 between parasitic nematodes and their closest free-living nematode, the specifically high expression in the parasitic female stage, and presence in excretory-secretory proteome of Strongyloides suggest that SOD3 may be related with parasitism.

Conclusions

This study advances our understanding of the complex evolution of antioxidant enzymes across Nematoda and provides targets for controlling parasitic nematode diseases.

Regional sequence expansion or collapse in heterozygous genome assemblies

High levels of heterozygosity present a unique genome assembly challenge and can adversely impact downstream analyses, yet is common in sequencing datasets obtained from non-model organisms. Here we show that by re-assembling a heterozygous dataset with variant parameters and different assembly algorithms, we are able to generate assemblies whose protein annotations are statistically enriched for specific gene ontology categories. While total assembly length was not significantly affected by assembly methodologies tested, the assemblies generated varied widely in fragmentation level and we show local assembly collapse or expansion underlying the enrichment or depletion of specific protein functional groups. We show that these statistically significant deviations in gene ontology groups can occur in seemingly high-quality assemblies, and result from difficult-to-detect local sequence expansion or contractions. Given the unpredictable interplay between assembly algorithm, parameter, and biological sequence data heterozygosity, we highlight the need for better measures of assembly quality than N50 value, including methods for assessing local expansion and collapse.

In the genomic era, genomes must be reconstructed from fragments using computational methods, or assemblers. How do we know that a new genome assembly is correct? This is important because errors in assembly can lead to downstream problems in gene predictions and these inaccurate results can contaminate databases, affecting later comparative studies. A particular challenge occurs when a diploid organism inherits two highly divergent genome copies from its parents. While it is widely appreciated that this type of data is difficult for assemblers to handle properly, here we show that the process is prone to more errors than previously appreciated. Specifically, we document examples of regional expansion and collapse, affecting downstream gene prediction accuracy, but without changing the overall genome assembly size or other metrics of accuracy. Our results suggest that assembly evaluation methods should be altered to identify whether regional expansions and collapses are present in the genome assembly.