Longitudinal Multi-omics Analyses Identify Responses of Megakaryocytes, Erythroid Cells, and Plasmablasts as Hallmarks of Severe COVID-19.

Temporal resolution of cellular features associated with a severe COVID-19 disease trajectory is needed for understanding skewed immune responses and defining predictors of outcome. Here, we performed a longitudinal multi-omics study using a two-center cohort of 14 patients. We analyzed the bulk transcriptome, bulk DNA methylome, and single-cell transcriptome (>358,000 cells, including BCR profiles) of peripheral blood samples harvested from up to 5 time points. Validation was performed in two independent cohorts of COVID-19 patients. Severe COVID-19 was characterized by an increase of proliferating, metabolically hyperactive plasmablasts. Coinciding with critical illness, we also identified an expansion of interferon-activated circulating megakaryocytes and increased erythropoiesis with features of hypoxic signaling. Megakaryocyte- and erythroid-cell-derived co-expression modules were predictive of fatal disease outcome. The study demonstrates broad cellular effects of SARS-CoV-2 infection beyond adaptive immune cells and provides an entry point toward developing biomarkers and targeted treatments of patients with COVID-19.

Blood transcriptome profiling reveals distinct gene networks induced by mRNA vaccination against COVID-19.

Messenger RNA (mRNA) vaccines represent a new class of vaccines that has been shown to be highly effective during the COVID-19 pandemic and that holds great potential for other preventative and therapeutic applications. While it is known that the transcriptional activity of various genes is altered following mRNA vaccination, identifying and studying gene networks could reveal important scientific insights that might inform future vaccine designs. In this study, we conducted an in-depth weighted gene correlation network analysis of the blood transcriptome before and 24 h after the second and third vaccination with licensed mRNA vaccines against COVID-19 in humans, following a prime vaccination with either mRNA or ChAdOx1 vaccines. Utilizing this unsupervised gene network analysis approach, we identified distinct modular networks of co-varying genes characterized by either an expressional up- or downregulation in response to vaccination. Downregulated networks were associated with cell metabolic processes and regulation of transcription factors, while upregulated networks were associated with myeloid differentiation, antigen presentation, and antiviral, interferon-driven pathways. Within this interferon-associated network, we identified highly connected hub genes such as STAT2 and RIGI and associated upstream transcription factors, potentially playing important regulatory roles in the vaccine-induced immune response. The expression profile of this network significantly correlated with S1-specific IgG levels at the follow-up visit in vaccinated individuals. Those findings could be corroborated in a second, independent cohort of mRNA vaccine recipients. Collectively, results from this modular gene network analysis enhance the understanding of mRNA vaccines from a systems immunology perspective. Influencing specific gene networks could lead to optimized vaccines that elicit augmented vaccine responses.

Metabolic and transcriptomic analyses of nectaries reveal differences in the mechanism of nectar production between monocots (Ananas comosus) and dicots (Nicotiana tabacum).

BACKGROUND: Nectar is offered by numerous flowering plants to attract pollinators. To date, the production and secretion of nectar have been analyzed mainly in eudicots, particularly rosids such as Arabidopsis. However, due to the enormous diversity of flowering plants, further research on other plant species, especially monocots, is needed. Ananas comosus (monocot) is an economically important species that is ideal for such analyses because it produces easily accessible nectar in sufficient quantities. In addition, the analyses were also carried out with Nicotiana tabacum (dicot, asterids) for comparison. RESULTS: We performed transcriptome sequencing (RNA-Seq) analyses of the nectaries of Ananas comosus and Nicotiana tabacum, to test whether the mechanisms described for nectar production and secretion in Arabidopsis are also present in these plant species. The focus of these analyses is on carbohydrate metabolism and transport (e.g., sucrose-phosphate synthases, invertases, sucrose synthases, SWEETs and further sugar transporters). In addition, the metabolites were analyzed in the nectar, nectaries and leaves of both plant species to address the question of whether concentration gradients for different metabolites exist between the nectaries and nectar The nectar of N. tabacum contains large amounts of glucose, fructose and sucrose, and the sucrose concentration in the nectar appears to be similar to the sucrose concentration in the nectaries. Nectar production and secretion in this species closely resemble corresponding processes in some other dicots, including sucrose synthesis in nectaries and sucrose secretion by SWEET9. The nectar of A. comosus also contains large amounts of glucose, fructose and sucrose and in this species the sucrose concentration in the nectar appears to be higher than the sucrose concentration in the nectaries. Furthermore, orthologs of SWEET9 generally appear to be absent in A. comosus and other monocots. Therefore, sucrose export by SWEETs from nectaries into nectar can be excluded; rather, other mechanisms, such as active sugar export or exocytosis, are more likely. CONCLUSION: The mechanisms of nectar production and secretion in N. tabacum appear to be largely similar to those in other dicots, whereas in the monocotyledonous species A. comosus, different synthesis and transport processes are involved.

18S rDNA gene metabarcoding of microeukaryotes and epi-endophytes in the holobiome of seven species of large brown algae.

Brown algae (Phaeophyceae) are habitat-forming species in coastal ecosystems and include kelp forests and seaweed beds that support a wide diversity of marine life. Host-associated microbial communities are an integral part of phaeophyte biology, and whereas the bacterial microbial partners have received considerable attention, the microbial eukaryotes associated with brown algae have hardly been studied. Here, we used broadly targeted "pan-eukaryotic" primers (metabarcoding) to investigate brown algal-associated eukaryotes (the eukaryome). Using this approach, we aimed to investigate the eukaryome of seven large brown algae that are important and common species in coastal ecosystems. We also aimed to assess whether these macroalgae harbor novel eukaryotic diversity and to ascribe putative functional roles to the host-associated eukaryome based on taxonomic affiliation and phylogenetic placement. We detected a significant diversity of microeukaryotic and algal lineages associated with the brown algal species investigated. The operational taxonomic units (OTUs) were taxonomically assigned to 10 of the eukaryotic major supergroups, including taxonomic groups known to be associated with seaweeds as epibionts, endobionts, parasites, and commensals. Additionally, we revealed previously unrecorded sequence types, including novel phaeophyte OTUs, particularly in the Fucus spp. samples, that may represent fucoid genomic variants, sequencing artifacts, or undescribed epi-/endophytes. Our results provide baseline data and technical insights that will be useful for more comprehensive seaweed eukaryome studies investigating the evidently lineage-rich and functionally diverse symbionts of brown algae.

Meta-transcriptomic comparison of two sponge holobionts feeding on coral- and macroalgal-dissolved organic matter.

BACKGROUND: Sponge holobionts (i.e., the host and its associated microbiota) play a key role in the cycling of dissolved organic matter (DOM) in marine ecosystems. On coral reefs, an ecological shift from coral-dominated to algal-dominated ecosystems is currently occurring. Given that benthic corals and macroalgae release different types of DOM, in different abundances and with different bioavailability to sponge holobionts, it is important to understand how the metabolic activity of the host and associated microbiota change in response to the exposure to both DOM sources. Here, we look at the differential gene expression of two sponge holobionts 6 hours after feeding on naturally sourced coral- and macroalgal-DOM using RNA sequencing and meta-transcriptomic analysis. RESULTS: We found a slight, but significant differential gene expression in the comparison between the coral- and macroalgal-DOM treatments in both the high microbial abundance sponge Plakortis angulospiculatus and the low microbial abundance sponge Haliclona vansoesti. In the hosts, processes that regulate immune response, signal transduction, and metabolic pathways related to cell proliferation were elicited. In the associated microbiota carbohydrate metabolism was upregulated in both treatments, but coral-DOM induced further lipid and amino acids biosynthesis, while macroalgal-DOM caused a stress response. These differences could be driven by the presence of distinct organic macronutrients in the two DOM sources and of small pathogens or bacterial virulence factors in the macroalgal-DOM. CONCLUSIONS: This work provides two new sponge meta-transcriptomes and a database of putative genes and genetic pathways that are involved in the differential processing of coral- versus macroalgal-DOM as food source to sponges with high and low abundances of associated microbes. These pathways include carbohydrate metabolism, signaling pathways, and immune responses. However, the differences in the meta-transcriptomic responses of the sponge holobionts after 6 hours of feeding on the two DOM sources were small. Longer-term responses to both DOM sources should be assessed to evaluate how the metabolism and the ecological function of sponges will be affected when reefs shift from coral towards algal dominance.

The temperature-regulated DEAD-box RNA helicase CrhR interactome: Autoregulation and photosynthesis-related transcripts.

RNA helicases play crucial functions in RNA biology. In plants, RNA helicases are encoded by large gene families, performing roles in abiotic stress responses, development, the post-transcriptional regulation of gene expression as well as house-keeping functions. Several of these RNA helicases are targeted to the organelles, mitochondria and chloroplasts. Cyanobacteria are the direct evolutionary ancestors of plant chloroplasts. The cyanobacterium Synechocystis 6803 encodes a single DEAD-box RNA helicase, CrhR, that is induced by a range of abiotic stresses, including low temperature. Though the ΔcrhR mutant exhibits a severe cold-sensitive phenotype, the physiological function(s) performed by CrhR have not been described. To identify transcripts interacting with CrhR, we performed RNA co-immunoprecipitation with extracts from a Synechocystis crhR deletion mutant expressing the FLAG-tagged native CrhR or a K57A mutated version with an anticipated enhanced RNA binding. The composition of the interactome was strikingly biased towards photosynthesis-associated and redox-controlled transcripts. A transcript highly enriched in all experiments was the crhR mRNA, suggesting an auto-regulatory molecular mechanism. The identified interactome explains the described physiological role of CrhR in response to the redox poise of the photosynthetic electron transport chain and characterizes CrhR as an enzyme with a diverse range of transcripts as molecular targets.

The Inhalable Mycobiome of Sawmill Workers: Exposure Characterization and Diversity.

Exposure to fungal spores has been associated with respiratory symptoms and allergic alveolitis among sawmill workers, but the complexity of sawmill workers' fungal exposure has been poorly studied. We characterized the fungal diversity in air samples from sawmill workers' breathing zones and identified differences in the richness, diversity, and taxonomic composition between companies, departments, wood types, and seasons. Full-shift personal inhalable dust samples (n = 86) collected from 11 industrial sawmill, sorting mill, and planer mill companies processing spruce and/or pine were subjected to DNA metabarcoding using the fungal internal transcribed spacer (ITS) region 2. The workers were exposed to a higher total number of operational taxonomic units (OTUs) in summer than in winter and when processing spruce than when processing pine. Workers in the saw department had the richest fungal exposure, followed by workers in the planing department and sorting of dry timber department. Sawmills explained 11% of the variation in the fungal community composition of the exposure, followed by season (5%) and department (3%). The fungal compositions of the exposures also differed between seasons, sawmills, wood types, and departments at the taxonomic level, ranging from the phylum to the species level. The differences in exposure diversity suggest that the potential health effects of fungal inhalation may also be different; hence, a risk assessment based on the fungal diversity differences should be performed. This study may serve as a basis for establishing a fungal profile of signature species that are specific for sawmills and that can be measured quantitatively in future risk assessments of sawmill workers.IMPORTANCE To gain more knowledge about exposure-response relationships, it is important to improve exposure characterization by comprehensively identifying the temporal and spatial fungal composition and diversity of inhalable dust at workplaces. The variation in the diverse fungal communities to which individuals are exposed in different seasons and sawmills suggests that variations in exposure-related health effects between seasons and companies can be expected. More importantly, the distinct fungal profiles between departments across companies indicate that workers in different job groups are differently exposed and that health risks can be department specific. DNA metabarcoding provides insight into a broad spectrum of airborne fungi that may serve as a basis for obtaining important knowledge about the fungi to which workers are exposed.

The RNA recognition motif protein CP33A is a global ligand of chloroplast mRNAs and is essential for plastid biogenesis and plant development.

Chloroplast RNA metabolism depends on a multitude of nuclear-encoded RNA-binding proteins (RBPs). Most known chloroplast RBPs address specific RNA targets and RNA-processing functions. However, members of the small chloroplast ribonucleoprotein family (cpRNPs) play a global role in processing and stabilizing chloroplast RNAs. Here, we show that the cpRNP CP33A localizes to a distinct sub-chloroplastic domain and is essential for chloroplast development. The loss of CP33A yields albino seedlings that exhibit aberrant leaf development and can only survive in the presence of an external carbon source. Genome-wide RNA association studies demonstrate that CP33A associates with all chloroplast mRNAs. For a given transcript, quantification of CP33A-bound versus free RNAs demonstrates that CP33A associates with the majority of most mRNAs analyzed. Our results further show that CP33A is required for the accumulation of a number of tested mRNAs, and is particularly relevant for unspliced and unprocessed precursor mRNAs. Finally, CP33A fails to associate with polysomes or to strongly co-precipitate with ribosomal RNA, suggesting that it defines a ribodomain that is separate from the chloroplast translation machinery. Collectively, these findings suggest that CP33A contributes to globally essential RNA processes in the chloroplasts of higher plants.

Stenotrophomonas lactitubi sp. nov. and Stenotrophomonas indicatrix sp. nov., isolated from surfaces with food contact.

Five Gram-stain-negative, rod-shaped, none-spore-forming isolates were obtained from biofilms on different sites of a milking machine in Germany. Another strain with similar morphological characteristics was isolated from dirty dishes. Based on phylogenetic analysis of the 16S rRNA and gyrB genes, all isolates were assigned to the genus Stenotrophomonas, but were divided into three different groups. Chemotaxonomic characterization of the isolates led to the detection of iso-C15 : 0 and anteiso-C15 : 0 as the predominant cellular fatty acids, as well as small amounts of the hydroxyl fatty acids iso-C11 : 0 3-OH, C12 : 0 3-OH and iso-C13 : 0 3-OH. One group could be assigned to the species Stenotrophomonas maltophilia, while the genome sequences of two groups displayed average nucleotide identity values of less than 94 % between each other and the genome sequences of the next related type strains Stenotrophomonas maltophilia ATCC 13637T and Stenotrophomonas rhizophila DSM 14405T. Further phylogenetic, phenotypic and chemotaxonomic analyses enabled the differentiation of these strains from these closely related species. They are therefore considered to represent two novel species, for which the names Stenotrophomonaslactitubi and Stenotrophomonasindicatrix are proposed, with strains M15T (=DSM 104152T=LMG29943T) and WS40T (=DSM28278T=LMG29942T) as type strains.

A draft genome of the neritid snail Theodoxus fluviatilis.

The neritid snail Theodoxus fluviatilis is found across habitats differing in salinity, from shallow waters along the coast of the Baltic Sea to lakes throughout Europe. Living close to the water surface makes this species vulnerable to changes in salinity in their natural habitat, and the lack of a free-swimming larval stage limits this species' dispersal. Together, these factors have resulted in a patchy distribution of quite isolated populations differing in their salinity tolerances. In preparation for investigating the mechanisms underlying the physiological differences in osmoregulation between populations that cannot be explained solely by phenotypic plasticity, we present here an annotated draft genome assembly for T. fluviatilis, generated using PacBio long reads, Illumina short reads, and transcriptomic data. While the total assembly size (1045 kb) is similar to those of related species, it remains highly fragmented (N scaffolds = 35,695; N50 = 74 kb) though moderately high in complete gene content (BUSCO single copy complete: 74.3%, duplicate: 2.6%, fragmented: 10.6%, missing: 12.5% using metazoa n = 954). Nevertheless, we were able to generate gene annotations of 21,220 protein-coding genes (BUSCO single copy complete: 65.1%, duplicate: 16.7%, fragmented: 9.1%, missing: 9.1% using metazoa n = 954). Not only will this genome facilitate comparative evolutionary studies across Gastropoda, as this is the first genome assembly for the basal snail family Neritidae, it will also greatly facilitate the study of salinity tolerance in this species. Additionally, we discuss the challenges of working with a species where high molecular weight DNA isolation is very difficult.

Oxidation of sulfur, hydrogen, and iron by metabolically versatile Hydrogenovibrio from deep sea hydrothermal vents.

Chemolithoautotrophic Hydrogenovibrio are ubiquitous and abundant at hydrothermal vents. They can oxidize sulfur, hydrogen, or iron, but none are known to use all three energy sources. This ability though would be advantageous in vents hallmarked by highly dynamic environmental conditions. We isolated three Hydrogenovibrio strains from vents along the Indian Ridge, which grow on all three electron donors. We present transcriptomic data from strains grown on iron, hydrogen, or thiosulfate with respective oxidation and autotrophic carbon dioxide (CO2) fixation rates, RubisCO activity, SEM, and EDX. Maximum estimates of one strain's oxidation potential were 10, 24, and 952 mmol for iron, hydrogen, and thiosulfate oxidation and 0.3, 1, and 84 mmol CO2 fixation, respectively, per vent per hour indicating their relevance for element cycling in-situ. Several genes were up- or downregulated depending on the inorganic electron donor provided. Although no known genes of iron-oxidation were detected, upregulated transcripts suggested iron-acquisition and so far unknown iron-oxidation-pathways.

Distinct patterns of genetic variation at low-recombining genomic regions represent haplotype structure.

Genomic regions sometimes show patterns of genetic variation distinct from the genome-wide population structure. Such deviations have often been interpreted to represent effects of selection. However, systematic investigation of whether and how non-selective factors, such as recombination rates, can affect distinct patterns has been limited. Here, we associate distinct patterns of genetic variation with reduced recombination rates in a songbird, the Eurasian blackcap (Sylvia atricapilla), using a new reference genome assembly, whole-genome resequenc- ing data and recombination maps. We find that distinct patterns of genetic variation reflect haplotype structure at genomic regions with different prevalence of reduced recombination rate across populations. At low-recombining regions shared in most populations, distinct patterns reflect conspicuous haplotypes segregating in multiple populations. At low-recombining regions found only in a few populations, distinct patterns represent variance among cryptic haplotypes within the low-recombining populations. With simulations, we confirm that these distinct patterns evolve neutrally by reduced recombination rate, on which the effects of selection can be overlaid. Our results highlight that distinct patterns of genetic variation can emerge through evolutionary reduction of local recombination rate. The recombination landscape as an evolvable trait therefore plays an important role determining the heterogeneous distribution of genetic variation along the genome.

Assessing the risk to develop a growing teratoma syndrome based on molecular and epigenetic subtyping as well as novel secreted biomarkers.

In germ cell tumors (GCT), a growing teratoma during chemotherapy with decreasing tumor markers was defined as 'growing teratoma syndrome' (GTS) by Logothetis et al. in 1982. So far, its pathogenesis and specific treatment options remain elusive. We aimed at updating the GTS definition based on molecular and epigenetic features as well as identifying circulating biomarkers. We selected 50 GTS patients for clinical characterization and subsequently 12 samples were molecularly analyzed. We further included 7 longitudinal samples of 2 GTS patients. Teratomas (TER) showing no features of GTS served as controls. GTS were stratified based on growth rates into a slow (<0.5 cm/month), medium (0.5-1.5) and rapid (>1.5) group. By analyzing DNA methylation, microRNA expression and the secretome, we identified putative epigenetic and secreted biomarkers for the GTS subgroups. We found that proteins enriched in the GTS groups compared to TER were involved in proliferation, DNA replication and the cell cycle, while proteins interacting with the immune system were depleted. Additionally, GTSrapid seem to interact more strongly with the surrounding microenvironment than GTSslow. Expression of pluripotency- and yolk-sac tumor-associated genes in GTS and formation of a yolk-sac tumor or somatic-type malignancy in the longitudinal GTS samples, pointed at an additional occult non-seminomatous component after chemotherapy. Thus, updating the Logothetis GTS definition is necessary, which we propose as follows: The GTS describes a continuously growing teratoma that might harbor occult non-seminomatous components considerably reduced during therapy but outgrowing over time again.

Green targeting predictor and ambiguous targeting predictor 2: the pitfalls of plant protein targeting prediction and of transient protein expression in heterologous systems.

The challenges of plant protein targeting prediction are the existence of dual subcellular targets and the bias of experimentally confirmed data towards few and mostly nonplant model species. To assess whether training with proteins from evolutionarily distant species has a negative impact on prediction accuracy, we developed the Green Targeting Predictor tool, which was trained with a species-specific data set for Physcomitrella patens. Its performance was compared with that of the same tool trained with a mixed data set. In addition, we updated the Ambiguous Targeting Predictor. We found that predictions deviated from in vivo observations predominantly for proteins diverging within the green lineage, as well as for dual targeted proteins. To evaluate the usefulness of heterologous expression systems, selected proteins were subjected to localization studies in P. patens, Arabidopsis thaliana and Nicotiana tabacum. Four out of six proteins that show dual targeting in the original plant system were located only in a single compartment in one or both heterologous systems. We conclude that targeting signals of divergent plant species exhibit differences, calling for custom in silico and in vivo approaches when aiming to unravel the actual distribution patterns of proteins within a plant cell.

Antiviral RNAi Response in Culex quinquefasciatus-Derived HSU Cells.

Culex spp. mosquitoes are important vectors of viruses, such as West Nile virus, Eastern equine encephalitis virus and Rift valley fever virus. However, their interactions with innate antiviral immunity, especially RNA interference (RNAi), are not well known. Most research on RNAi pathways in mosquitoes is focused on the tropical vector mosquito Aedes aegypti. Here, we investigated the production of arbovirus-specific small RNAs in Cx. quinquefasciatus-derived HSU cells. Furthermore, by silencing RNAi-related proteins, we investigated the antiviral role of these proteins for two different arboviruses: Semliki Forest virus (SFV) and Bunyamwera orthobunyavirus (BUNV). Our results showed an expansion of Ago2 and Piwi6 in Cx. quinquefasciatus compared to Ae. aegypti. While silencing Ago2a and Ago2b increased BUNV replication, only Ago2b showed antiviral activity against SFV. Our results suggest differences in the function of Cx. quinquefasciatus and Ae. aegypti RNAi proteins and highlight the virus-specific function of these proteins in Cx. quinquefasciatus.

Phylogroup-specific variation shapes the clustering of antimicrobial resistance genes and defence systems across regions of genome plasticity in Pseudomonas aeruginosa.

BACKGROUND: Pseudomonas aeruginosa is an opportunistic pathogen consisting of three phylogroups (hereafter named A, B, and C). Here, we assessed phylogroup-specific evolutionary dynamics across available and also new P. aeruginosa genomes. METHODS: In this genomic analysis, we first generated new genome assemblies for 18 strains of the major P. aeruginosa clone type (mPact) panel, comprising a phylogenetically diverse collection of clinical and environmental isolates for this species. Thereafter, we combined these new genomes with 1991 publicly available P. aeruginosa genomes for a phylogenomic and comparative analysis. We specifically explored to what extent antimicrobial resistance (AMR) genes, defence systems, and virulence genes vary in their distribution across regions of genome plasticity (RGPs) and "masked" (RGP-free) genomes, and to what extent this variation differs among the phylogroups. FINDINGS: We found that members of phylogroup B possess larger genomes, contribute a comparatively larger number of pangenome families, and show lower abundance of CRISPR-Cas systems. Furthermore, AMR and defence systems are pervasive in RGPs and integrative and conjugative/mobilizable elements (ICEs/IMEs) from phylogroups A and B, and the abundance of these cargo genes is often significantly correlated. Moreover, inter- and intra-phylogroup interactions occur at the accessory genome level, suggesting frequent recombination events. Finally, we provide here the mPact panel of diverse P. aeruginosa strains that may serve as a valuable reference for functional analyses. INTERPRETATION: Altogether, our results highlight distinct pangenome characteristics of the P. aeruginosa phylogroups, which are possibly influenced by variation in the abundance of CRISPR-Cas systems and are shaped by the differential distribution of other defence systems and AMR genes. FUNDING: German Science Foundation, Max-Planck Society, Leibniz ScienceCampus Evolutionary Medicine of the Lung, BMBF program Medical Infection Genomics, Kiel Life Science Postdoc Award.

Agua Salud Alphavirus Infection, Dissemination and Transmission in Aedes aegypti Mosquitoes.

Mosquitoes are competent vectors for many important arthropod-borne viruses (arboviruses). In addition to arboviruses, insect-specific viruses (ISV) have also been discovered in mosquitoes. ISVs are viruses that replicate in insect hosts but are unable to infect and replicate in vertebrates. They have been shown to interfere with arbovirus replication in some cases. Despite the increase in studies on ISV-arbovirus interactions, ISV interactions with their hosts and how they are maintained in nature are still not well understood. In the present study, we investigated the infection and dissemination of the Agua Salud alphavirus (ASALV) in the important mosquito vector Aedes aegypti through different infection routes (per oral infection, intrathoracic injection) and its transmission. We show here that ASALV infects the female Ae. aegypti and replicates when mosquitoes are infected intrathoracically or orally. ASALV disseminated to different tissues, including the midgut, salivary glands and ovaries. However, we observed a higher virus load in the brain than in the salivary glands and carcasses, suggesting a tropism towards brain tissues. Our results show that ASALV is transmitted horizontally during adult and larval stages, although we did not observe vertical transmission. Understanding ISV infection and dissemination dynamics in Ae. aegypti and their transmission routes could help the use of ISVs as an arbovirus control strategy in the future.

The Aedes aegypti RNA interference response against Zika virus in the context of co-infection with dengue and chikungunya viruses.

Since its detection in 2015 in Brazil, Zika virus (ZIKV) has remained in the spotlight of international public health and research as an emerging arboviral pathogen. In addition to single infection, ZIKV may occur in co-infection with dengue (DENV) and chikungunya (CHIKV) viruses, with whom ZIKV shares geographic distribution and the mosquito Aedes aegypti as a vector. The main mosquito immune response against arboviruses is RNA interference (RNAi). It is unknown whether or not the dynamics of the RNAi response differ between single arboviral infections and co-infections. In this study, we investigated the interaction of ZIKV and DENV, as well as ZIKV and CHIKV co-infections with the RNAi response in Ae. aegypti. Using small RNA sequencing, we found that the efficiency of small RNA production against ZIKV -a hallmark of antiviral RNAi-was mostly similar when comparing single and co-infections with either DENV or CHIKV. Silencing of key antiviral RNAi proteins, showed no change in effect on ZIKV replication when the cell is co-infected with ZIKV and DENV or CHIKV. Interestingly, we observed a negative effect on ZIKV replication during CHIKV co-infection in the context of Ago2-knockout cells, though his effect was absent during DENV co-infection. Overall, this study provides evidence that ZIKV single or co-infections with CHIKV or DENV are equally controlled by RNAi responses. Thus, Ae. aegypti mosquitoes and derived cells support co-infections of ZIKV with either CHIKV or DENV to a similar level than single infections, as long as the RNAi response is functional.

Bacterial c-di-GMP has a key role in establishing host-microbe symbiosis.

Most microbes evolve faster than their hosts and should therefore drive evolution of host-microbe interactions. However, relatively little is known about the characteristics that define the adaptive path of microbes to host association. Here we identified microbial traits that mediate adaptation to hosts by experimentally evolving the free-living bacterium Pseudomonas lurida with the nematode Caenorhabditis elegans as its host. After ten passages, we repeatedly observed the evolution of beneficial host-specialist bacteria, with improved persistence in the nematode being associated with increased biofilm formation. Whole-genome sequencing revealed mutations that uniformly upregulate the bacterial second messenger, cyclic diguanylate (c-di-GMP). We subsequently generated mutants with upregulated c-di-GMP in different Pseudomonas strains and species, which consistently increased host association. Comparison of pseudomonad genomes from various environments revealed that c-di-GMP underlies adaptation to a variety of hosts, from plants to humans. This study indicates that c-di-GMP is fundamental for establishing host association.

Entorhinal cortex epigenome-wide association study highlights four novel loci showing differential methylation in Alzheimer's disease.

BACKGROUND: Studies on DNA methylation (DNAm) in Alzheimer's disease (AD) have recently highlighted several genomic loci showing association with disease onset and progression. METHODS: Here, we conducted an epigenome-wide association study (EWAS) using DNAm profiles in entorhinal cortex (EC) from 149 AD patients and control brains and combined these with two previously published EC datasets by meta-analysis (total n = 337). RESULTS: We identified 12 cytosine-phosphate-guanine (CpG) sites showing epigenome-wide significant association with either case-control status or Braak's tau-staging. Four of these CpGs, located in proximity to CNFN/LIPE, TENT5A, PALD1/PRF1, and DIRAS1, represent novel findings. Integrating DNAm levels with RNA sequencing-based mRNA expression data generated in the same individuals showed significant DNAm-mRNA correlations for 6 of the 12 significant CpGs. Lastly, by calculating rates of epigenetic age acceleration using two recently proposed "epigenetic clock" estimators we found a significant association with accelerated epigenetic aging in the brains of AD patients vs. controls. CONCLUSION: In summary, our study represents the hitherto most comprehensive EWAS in AD using EC and highlights several novel differentially methylated loci with potential effects on gene expression.