Cranberry-derived proanthocyanidins induce a differential transcriptomic response within Candida albicans urinary biofilms

Candida albicans is one of the most common causes of hospital-acquired urinary tract infections (UTIs). However, azoles are poorly active against biofilms, echinocandins do not achieve clinically useful urinary concentrations, and amphotericin B exhibits severe toxicities. Thus, novel strategies are needed to prevent Candida UTIs, which are often associated with urinary catheter biofilms. We previously demonstrated that cranberry-derived proanthocyanidins (PACs) prevent C. albicans biofilm formation in an in vitro urinary model. To elucidate functional pathways unique to urinary biofilm development and PAC inhibition, we investigated the transcriptome of C. albicans in artificial urine (AU), with and without PACs. C. albicans biofilm and planktonic cells were cultivated with or without PACs. Genome-wide expression analysis was performed by RNA sequencing. Differentially expressed genes were determined using DESeq2 software; pathway analysis was performed using Cytoscape. Approximately 2,341 of 6,444 total genes were significantly expressed in biofilm relative to planktonic cells. Functional pathway analysis revealed that genes involved in filamentation, adhesion, drug response and transport were up-regulated in urinary biofilms. Genes involved in carbon and nitrogen metabolism and nutrient response were down-regulated. In PAC-treated urinary biofilms compared to untreated control biofilms, 557 of 6,444 genes had significant changes in gene expression. Genes downregulated in PAC-treated biofilms were implicated in iron starvation and adhesion pathways. Although urinary biofilms share key features with biofilms formed in other environments, many genes are uniquely expressed in urinary biofilms. Cranberry-derived PACs interfere with the expression of iron acquisition and adhesion genes within urinary biofilms.

Transcriptomic responses of Biomphalaria pfeifferi to Schistosoma mansoni: Investigation of a neglected African snail that supports more S. mansoni transmission than any other snail species

Background

Biomphalaria pfeifferi is highly compatible with the widespread human-infecting blood fluke Schistosoma mansoni and transmits more cases of this parasite to people than any other snail species. For these reasons, B. pfeifferi is the world’s most important vector snail for S. mansoni, yet we know relatively little at the molecular level regarding the interactions between B. pfeifferi and S. mansoni from early-stage sporocyst transformation to the development of cercariae.

Methodology/Principal findings

We sought to capture a portrait of the response of B. pfeifferi to S. mansoni as it occurs in nature by undertaking Illumina dual RNA-Seq on uninfected control B. pfeifferi and three intramolluscan developmental stages (1- and 3-days post infection and patent, cercariae-producing infections) using field-derived west Kenyan specimens. A high-quality, well-annotated de novo B. pfeifferi transcriptome was assembled from over a half billion non-S. mansoni paired-end reads. Reads associated with potential symbionts were noted. Some infected snails yielded fewer normalized S. mansoni reads and showed different patterns of transcriptional response than others, an indication that the ability of field-derived snails to support and respond to infection is variable. Alterations in transcripts associated with reproduction were noted, including for the oviposition-related hormone ovipostatin and enzymes involved in metabolism of bioactive amines like dopamine or serotonin. Shedding snails exhibited responses consistent with the need for tissue repair. Both generalized stress and immune factors immune factors (VIgLs, PGRPs, BGBPs, complement C1q-like, chitinases) exhibited complex transcriptional responses in this compatible host-parasite system.

Significance

This study provides for the first time a large sequence data set to help in interpreting the important vector role of the neglected snail B. pfeifferi in transmission of S. mansoni, including with an emphasis on more natural, field-derived specimens. We have identified B. pfeifferi targets particularly responsive during infection that enable further dissection of the functional role of these candidate molecules.

Biomphalaria pfeifferi is the world’s most important snail vector for the widespread human-infecting blood fluke Schistosoma mansoni. Despite this, we know relatively little about the biology of this highly compatible African snail host of S. mansoni, especially for specimens from the field. Using an Illumina-based dual-seq approach, we captured a portrait of the transcriptional responses of Kenyan snails that were either uninfected with S. mansoni, or that harbored 1-day, 3-day, or cercariae-producing infections. Responses to infection were influenced both by the extent of schistosome gene expression and infection duration. We note and discuss several alterations in transcriptional activity in immune, stress and reproduction related genes in infected snails and the B. pfeifferi symbionts detected. Several host genes were highly up-regulated following infection and these might comprise excellent candidates for disruption to diminish compatibility. This study provides for the first time a large sequence dataset to help in interpreting the important vector role of B. pfeifferi in transmission of S. mansoni, including with an emphasis on more natural, field-derived specimens.

RNA-Seq Comparison of Larval and Adult Malpighian Tubules of the Yellow Fever Mosquito Aedes aegypti Reveals Life Stage-Specific Changes in Renal Function

Introduction: The life history of Aedes aegypti presents diverse challenges to its diuretic system. During the larval and pupal life stages mosquitoes are aquatic. With the emergence of the adult they become terrestrial. This shifts the organism within minutes from an aquatic environment to a terrestrial environment where dehydration has to be avoided. In addition, female mosquitoes take large blood meals, which present an entirely new set of challenges to salt and water homeostasis.

Methods: To determine differences in gene expression associated with these different life stages, we performed an RNA-seq analysis of the main diuretic tissue in A. aegypti, the Malpighian tubules. We compared transcript abundance in 4th instar larvae to that of adult females and analyzed the data with a focus on transcripts that encode proteins potentially involved in diuresis, like water and solute channels as well as ion transporters. We compared our results against the model of potassium- and sodium chloride excretion in the Malpighian tubules proposed by Hine et al. (2014), which involves at least eight ion transporters and a proton-pump.

Results: We found 3,421 of a total number of 17,478 (19.6%) unique transcripts with a P < 0.05 and at least a 2.5 fold change in expression levels between the two groups. We identified two novel transporter genes that are highly expressed in the adult Malpighian tubules, which have not previously been part of the transport model in this species and may play important roles in diuresis. We also identified candidates for hypothesized sodium and chloride channels. Detoxification genes were generally higher expressed in larvae.

Significance: This study represents the first comparison of Malpighian tubule transcriptomes between larval and adult A. aegypti mosquitoes, highlighting key differences in their renal systems that arise as they transform from an aquatic filter-feeding larval stage to a terrestrial, blood-feeding adult stage.

Zinc-Dependent Transcriptional Regulation in Paracoccus denitrificans

Zinc homeostasis is critical for bacterial survival and is mediated largely at the transcriptional level by the regulation of zinc uptake and efflux genes. Here we use RNA-seq to assess transcriptional changes as a result of zinc limitation in the denitrifying bacterium Paracoccus denitrificans. The results identify the differential expression of 147 genes, most of which were upregulated in zinc-depleted medium. Included in this set of genes are a large number of transition metal transporters, several transcription factors, and hypothetical proteins. Intriguingly, genes encoding nitric oxide reductase (norCB) and nitrite reductase (nirS) were also upregulated. A Zur consensus binding motif was identified in the promoters of the most highly upregulated genes. The zinc uptake regulator (Zur) from this organism was also characterized and shown to bind to the Zur motif in a zinc-dependent manner. This work expands our current understanding of the transcriptional response of gram-negative bacteria to zinc limitation and identifies genes involved in denitrification as part of the Zur regulon.

Dengue virus serotype 2 infection alters midgut and carcass gene expression in the Asian tiger mosquito, Aedes albopictus

Background

The Asian tiger mosquito, Aedes albopictus is currently an important vector for dengue, chikungunya and Zika virus, and its role in transmission of arthropod-borne viruses (arboviruses) may increase in the future due to its ability to colonize temperate regions. In contrast to Aedes aegypti, the dominant vector of dengue, chikungunya and Zika virus, genetic responses of Ae. albopictus upon infection with an arbovirus are not well characterized. Here we present a study of the changes in transcript expression in Ae. albopictus exposed to dengue virus serotype 2 via feeding on an artificial bloodmeal.

Methodology/Principal findings

We isolated midguts and midgut-free carcasses of Ae. albopictus fed on bloodmeals containing dengue virus as well as controls fed on virus-free control meals at day 1 and day 5 post-feeding. We confirmed infection of midguts from mosquitoes sampled on day 5 post-feeding via RT-PCR. RNAseq analysis revealed dynamic modulation of the expression of several putative immunity and dengue virus-responsive genes, some of whose expression was verified by qRT-PCR. For example, a serine protease gene was up-regulated in the midgut at 1 day post infection, which may potentially enhance mosquito susceptibility to dengue infection, while 14 leucine-rich repeat genes, previously shown to be involved in mosquito antiviral defenses, were down-regulated in the carcass at 5 days post infection. The number of significantly modulated genes decreased over time in midguts and increased in carcasses.

Conclusion/Significance

Dengue virus exposure results in the modulation of genes in a time- and site-specific manner. Previous literature on the interaction between mosquitoes and mosquito-borne pathogens suggests that most of the changes that occurred in Ae. albopictus exposed to DENV would favor virus infection. Many genes identified in this study warrant further characterization to understand their role in viral manipulation of and antiviral response of Ae. albopictus.

De novo assembly of a tadpole shrimp (Triops newberryi) transcriptome and preliminary differential gene expression analysis

Next-generation sequencing techniques, such as RNA sequencing, have provided a wealth of genomic information for nonmodel species. Transcriptomic information can be used to quantify the patterns of gene expression, which can identify how environmental differences invoke organismal stress responses and provide a gauge in predicting species adaptability. In our study, we used RNA sequencing to characterize the first transcriptome from a naupliar tadpole shrimp (Triops newberryi) to identify the genes expressed during the early life history stages and which could be important for future genomic studies. RNA was extracted from naupliar T. newberryi that were reared in a laboratory-controlled setting and in two different water types, a native and a non-native condition. A total of six replicates, three per condition, were sequenced with the Illumina Hi-Seq 2000 achieving 365 M 50-nt reads. High-quality reads were produced and de novo assembly was used to construct a T. newberryi transcriptome that was approximately 24.8 M base pairs. More than 10 000 peptides were predicted from the assembly, and genes were sorted into gene ontology categories. The use of different water conditions allowed for a preliminary differential gene expression analysis in order to compare the changes in gene expression between conditions. There were 299 differentially expressed genes between water conditions that might serve as a focal point for future genomic studies of Triops acclimation to different environments. The Triops transcriptome could serve as vital genomic information for additional studies on Branchiopod crustaceans.

Small mosquitoes, large implications: crowding and starvation affects gene expression and nutrient accumulation in Aedes aegypti

Background

Environmental factors such as temperature, nutrient availability, and larval density determine the outcome of postembryonic development in mosquitoes. Suboptimal temperatures, crowding, and starvation during the larval phase reduce adult mosquito size, nutrient stores and affect vectorial capacity.

Methods

In this study we compared adult female Aedes aegypti, Rockefeller strain, raised under standard laboratory conditions (Large) with those raised under crowded and nutritionally deprived conditions (Small). To compare the gene expression and nutritional state of the major energy storage and metabolic organ, the fat body, we performed transcriptomics using Illumina based RNA-seq and metabolomics using GC/MS on females before and 24 hours following blood feeding.

Results

Analysis of fat body gene expression between the experimental groups revealed a large number of significantly differentially expressed genes. Transcripts related to immunity, reproduction, autophagy, several metabolic pathways; including amino acid degradation and metabolism; and membrane transport were differentially expressed. Metabolite profiling identified 60 metabolites within the fat body to be significantly affected between small and large mosquitoes, with the majority of detected free amino acids at a higher level in small mosquitoes compared to large.

Conclusions

Gene expression and metabolites in the adult fat body reflect the individual post-embryonic developmental history of a mosquito larva. These changes affect nutritional storage and utilization, immunity, and reproduction. Therefore, it is apparent that changes in larval environment due to weather conditions, nutrition availability, vector control efforts, and other factors can affect adult vectorial capacity in the field.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-0863-9) contains supplementary material, which is available to authorized users.

Comparative genome analysis of non-toxigenic non-O1 versus toxigenic O1 Vibrio cholerae

Pathogenic strains of Vibrio cholerae are responsible for endemic and pandemic outbreaks of the disease cholera. The complete toxigenic mechanisms underlying virulence in Vibrio strains are poorly understood. The hypothesis of this work was that virulent versus non-virulent strains of V. cholerae harbor distinctive genomic elements that encode virulence. The purpose of this study was to elucidate genomic differences between the O1 serotypes and non-O1 V. cholerae PS15, a non-toxigenic strain, in order to identify novel genes potentially responsible for virulence. In this study, we compared the whole genome of the non-O1 PS15 strain to the whole genomes of toxigenic serotypes at the phylogenetic level, and found that the PS15 genome was distantly related to those of toxigenic V. cholerae. Thus we focused on a detailed gene comparison between PS15 and the distantly related O1 V. cholerae N16961. Based on sequence alignment we tentatively assigned chromosome numbers 1 and 2 to elements within the genome of non-O1 V. cholerae PS15. Further, we found that PS15 and O1 V. cholerae N16961 shared 98% identity and 766 genes, but of the genes present in N16961 that were missing in the non-O1 V. cholerae PS15 genome, 56 were predicted to encode not only for virulence-related genes (colonization, antimicrobial resistance, and regulation of persister cells) but also genes involved in the metabolic biosynthesis of lipids, nucleosides and sulfur compounds. Additionally, we found 113 genes unique to PS15 that were predicted to encode other properties related to virulence, disease, defense, membrane transport, and DNA metabolism. Here, we identified distinctive and novel genomic elements between O1 and non-O1 V. cholerae genomes as potential virulence factors and, thus, targets for future therapeutics. Modulation of such novel targets may eventually enhance eradication efforts of endemic and pandemic disease cholera in afflicted nations.

Analysis of the transcriptomes downstream of Eyeless and the Hedgehog, Decapentaplegic and Notch signaling pathways in Drosophila melanogaster

Tissue-specific transcription factors are thought to cooperate with signaling pathways to promote patterned tissue specification, in part by co-regulating transcription. The Drosophila melanogaster Pax6 homolog Eyeless forms a complex, incompletely understood regulatory network with the Hedgehog, Decapentaplegic and Notch signaling pathways to control eye-specific gene expression. We report a combinatorial approach, including mRNAseq and microarray analyses, to identify targets co-regulated by Eyeless and Hedgehog, Decapentaplegic or Notch. Multiple analyses suggest that the transcriptomes resulting from co-misexpression of Eyeless+signaling factors provide a more complete picture of eye development compared to previous efforts involving Eyeless alone: (1) Principal components analysis and two-way hierarchical clustering revealed that the Eyeless+signaling factor transcriptomes are closer to the eye control transcriptome than when Eyeless is misexpressed alone; (2) more genes are upregulated at least three-fold in response to Eyeless+signaling factors compared to Eyeless alone; (3) based on gene ontology analysis, the genes upregulated in response to Eyeless+signaling factors had a greater diversity of functions compared to Eyeless alone. Through a secondary screen that utilized RNA interference, we show that the predicted gene CG4721 has a role in eye development. CG4721 encodes a neprilysin family metalloprotease that is highly up-regulated in response to Eyeless+Notch, confirming the validity of our approach. Given the similarity between D. melanogaster and vertebrate eye development, the large number of novel genes identified as potential targets of Ey+signaling factors will provide novel insights to our understanding of eye development in D. melanogaster and humans.

Carrier testing for severe childhood recessive diseases by next-generation sequencing

Of 7028 disorders with suspected Mendelian inheritance, 1139 are recessive and have an established molecular basis. Although individually uncommon, Mendelian diseases collectively account for ~20% of infant mortality and ~10% of pediatric hospitalizations. Preconception screening, together with genetic counseling of carriers, has resulted in remarkable declines in the incidence of several severe recessive diseases including Tay-Sachs disease and cystic fibrosis. However, extension of preconception screening to most severe disease genes has hitherto been impractical. Here, we report a preconception carrier screen for 448 severe recessive childhood diseases. Rather than costly, complete sequencing of the human genome, 7717 regions from 437 target genes were enriched by hybrid capture or microdroplet polymerase chain reaction, sequenced by next-generation sequencing (NGS) to a depth of up to 2.7 gigabases, and assessed with stringent bioinformatic filters. At a resultant 160x average target coverage, 93% of nucleotides had at least 20x coverage, and mutation detection/genotyping had ~95% sensitivity and ~100% specificity for substitution, insertion/deletion, splicing, and gross deletion mutations and single-nucleotide polymorphisms. In 104 unrelated DNA samples, the average genomic carrier burden for severe pediatric recessive mutations was 2.8 and ranged from 0 to 7. The distribution of mutations among sequenced samples appeared random. Twenty-seven percent of mutations cited in the literature were found to be common polymorphisms or misannotated, underscoring the need for better mutation databases as part of a comprehensive carrier testing strategy. Given the magnitude of carrier burden and the lower cost of testing compared to treating these conditions, carrier screening by NGS made available to the general population may be an economical way to reduce the incidence of and ameliorate suffering associated with severe recessive childhood disorders.

C6/36 Aedes albopictus cells have a dysfunctional antiviral RNA interference response

Mosquitoes rely on RNA interference (RNAi) as their primary defense against viral infections. To this end, the combination of RNAi and invertebrate cell culture systems has become an invaluable tool in studying virus-vector interactions. Nevertheless, a recent study failed to detect an active RNAi response to West Nile virus (WNV) infection in C6/36 (Aedes albopictus) cells, a mosquito cell line frequently used to study arthropod-borne viruses (arboviruses). Therefore, we sought to determine if WNV actively evades the host's RNAi response or if C6/36 cells have a dysfunctional RNAi pathway. C6/36 and Drosophila melanogaster S2 cells were infected with WNV (Flaviviridae), Sindbis virus (SINV, Togaviridae) and La Crosse virus (LACV, Bunyaviridae) and total RNA recovered from cell lysates. Small RNA (sRNA) libraries were constructed and subjected to high-throughput sequencing. In S2 cells, virus-derived small interfering RNAs (viRNAs) from all three viruses were predominantly 21 nt in length, a hallmark of the RNAi pathway. However, in C6/36 cells, viRNAs were primarily 17 nt in length from WNV infected cells and 26–27 nt in length in SINV and LACV infected cells. Furthermore, the origin (positive or negative viral strand) and distribution (position along viral genome) of S2 cell generated viRNA populations was consistent with previously published studies, but the profile of sRNAs isolated from C6/36 cells was altered. In total, these results suggest that C6/36 cells lack a functional antiviral RNAi response. These findings are analogous to the type-I interferon deficiency described in Vero (African green monkey kidney) cells and suggest that C6/36 cells may fail to accurately model mosquito-arbovirus interactions at the molecular level.

Cell culture systems are invaluable tools for studying virus-host interactions. These systems are typically easy to maintain and manipulate; however, they can fail to accurately mimic the host environment encountered by viruses. Therefore, defining the limitations of each system is critical to properly interpreting the results. C6/36 Aedes albopictus cells are commonly used to study arthropod-borne viruses (arboviruses), such as West Nile virus (WNV). Recent evidence suggests that the RNA interference (RNAi) pathway, a critical aspect of the cellular innate antiviral immune response in invertebrates, may not actively target WNV in C6/36 cells. However, it is unknown whether this observation is limited to WNV. Therefore, we examined small RNA populations from C6/36 and Drosophila melanogastor S2 cells infected with WNV, Sindbis virus and La Crosse virus by high-throughput sequencing. We demonstrate that the RNAi pathway actively targets each of the three viruses in S2 cells, but does not in C6/36 cells. These findings suggest that C6/36 cells may fail to accurately model mosquito-arbovirus interactions.

Genome, epigenome and RNA sequences of monozygotic twins discordant for multiple sclerosis

Identical (or more correctly 'monozygotic') twins are widely used to study the contributions of genetics and environment to human disease. A study that focused on three pairs of monozygotic twins, in which one twin had multiple sclerosis and the other did not, has brought the latest techniques of genome sequencing and analysis to this field, and incidentally published the first female human genome sequences. Full sequences were determined for one pair of twins, and for these and the other two pairs the mRNA transcriptome and epigenome sequences of CD4⁺ lymphocytes were determined. The striking result is that no genetic, epigenetic or transcriptome differences were found that explained why one twin had the disease and the other did not. Digging deeper into the data, eQTL (expression quantitative trait locus) mapping revealed tantalizing differences within twin pairs that merit closer examination. And some possible causes can be ruled out. Future work might usefully concentrate on studies of other cell types and epigenetic modifications.

Studies of identical twins are widely used to dissect the contributions of genes and the environment to human diseases. In multiple sclerosis, an autoimmune demyelinating disease, identical twins often show differences. This might suggest that environmental effects are most significant in this case, but genetic and epigenetic differences between identical twins have been described. Here, however, studies of identical twins show no evidence for genetic, epigenetic or transcriptome differences that could explain disease discordance.

Monozygotic or ‘identical’ twins have been widely studied to dissect the relative contributions of genetics and environment in human diseases. In multiple sclerosis (MS), an autoimmune demyelinating disease and common cause of neurodegeneration and disability in young adults, disease discordance in monozygotic twins has been interpreted to indicate environmental importance in its pathogenesis^{1,2,3,4,5,6,7,8}. However, genetic and epigenetic differences between monozygotic twins have been described, challenging the accepted experimental model in disambiguating the effects of nature and nurture^9,10,11,12. Here we report the genome sequences of one MS-discordant monozygotic twin pair, and messenger RNA transcriptome and epigenome sequences of CD4⁺ lymphocytes from three MS-discordant, monozygotic twin pairs. No reproducible differences were detected between co-twins among ∼3.6 million single nucleotide polymorphisms (SNPs) or ∼0.2 million insertion-deletion polymorphisms. Nor were any reproducible differences observed between siblings of the three twin pairs in HLA haplotypes, confirmed MS-susceptibility SNPs, copy number variations, mRNA and genomic SNP and insertion-deletion genotypes, or the expression of ∼19,000 genes in CD4⁺ T cells. Only 2 to 176 differences in the methylation of ∼2 million CpG dinucleotides were detected between siblings of the three twin pairs, in contrast to ∼800 methylation differences between T cells of unrelated individuals and several thousand differences between tissues or between normal and cancerous tissues. In the first systematic effort to estimate sequence variation among monozygotic co-twins, we did not find evidence for genetic, epigenetic or transcriptome differences that explained disease discordance. These are the first, to our knowledge, female, twin and autoimmune disease individual genome sequences reported.

Genomic convergence analysis of schizophrenia: mRNA sequencing reveals altered synaptic vesicular transport in post-mortem cerebellum

Schizophrenia (SCZ) is a common, disabling mental illness with high heritability but complex, poorly understood genetic etiology. As the first phase of a genomic convergence analysis of SCZ, we generated 16.7 billion nucleotides of short read, shotgun sequences of cDNA from post-mortem cerebellar cortices of 14 patients and six, matched controls. A rigorous analysis pipeline was developed for analysis of digital gene expression studies. Sequences aligned to approximately 33,200 transcripts in each sample, with average coverage of 450 reads per gene. Following adjustments for confounding clinical, sample and experimental sources of variation, 215 genes differed significantly in expression between cases and controls. Golgi apparatus, vesicular transport, membrane association, Zinc binding and regulation of transcription were over-represented among differentially expressed genes. Twenty three genes with altered expression and involvement in presynaptic vesicular transport, Golgi function and GABAergic neurotransmission define a unifying molecular hypothesis for dysfunction in cerebellar cortex in SCZ.

Management of High-Throughput DNA Sequencing Projects: Alpheus

High-throughput DNA sequencing has enabled systems biology to begin to address areas in health, agricultural and basic biological research. Concomitant with the opportunities is an absolute necessity to manage significant volumes of high-dimensional and inter-related data and analysis. Alpheus is an analysis pipeline, database and visualization software for use with massively parallel DNA sequencing technologies that feature multi-gigabase throughput characterized by relatively short reads, such as Illumina-Solexa (sequencing-by-synthesis), Roche-454 (pyrosequencing) and Applied Biosystem's SOLiD (sequencing-by-ligation). Alpheus enables alignment to reference sequence(s), detection of variants and enumeration of sequence abundance, including expression levels in transcriptome sequence. Alpheus is able to detect several types of variants, including non-synonymous and synonymous single nucleotide polymorphisms (SNPs), insertions/deletions (indels), premature stop codons, and splice isoforms. Variant detection is aided by the ability to filter variant calls based on consistency, expected allele frequency, sequence quality, coverage, and variant type in order to minimize false positives while maximizing the identification of true positives. Alpheus also enables comparisons of genes with variants between cases and controls or bulk segregant pools. Sequence-based differential expression comparisons can be developed, with data export to SAS JMP Genomics for statistical analysis.