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

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. However, genetic and epigenetic differences between monozygotic twins have been described, challenging the accepted experimental model in disambiguating the effects of nature and nurture. 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 approximately 3.6 million single nucleotide polymorphisms (SNPs) or approximately 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 approximately 19,000 genes in CD4(+) T cells. Only 2 to 176 differences in the methylation of approximately 2 million CpG dinucleotides were detected between siblings of the three twin pairs, in contrast to approximately 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.

Identification and quantitation of multiple variants in RNA virus genomes.

The goal of the study was to identify and characterize RNA virus variants containing mutations spread over genomic distances >5 kb. As proof of concept, high-quality viral RNA of the Dengue 2 component of Takeda's tetravalent dengue vaccine candidate (TDV-2) was used to develop a reverse transcription-polymerase chain reaction protocol to amplify a ∼5.3 kb cDNA segment that contains the three genetic determinants of TDV-2 attenuation. Unique molecular identifiers were incorporated into each viral cDNA molecule for PacBio library preparation to improve the quantitative precision of the observed variants at the attenuation loci. Following assay optimization, PacBio long-read sequencing was validated with multiple clone-derived TDV-2 revertant variants and four complex revertant mixtures containing various compositions of TDV-2 and revertant viruses. PacBio sequencing analysis correctly identified and quantified variant composition in all tested samples, demonstrating that TDV-2 revertants could be identified and characterized and supporting the use of this method in the differentiation and quantification of complex variants of other RNA viruses. Long-read sequencing can identify complex RNA virus variants containing multiple mutations on a single-genome molecule, which is useful for in-depth genetic stability and revertant detection of live-attenuated viral vaccines, as well as research in virus evolution to reveal mechanisms of immune evasion and host cell adaption.

Genome-wide association genetics of an adaptive trait in lodgepole pine.

Pine cones that remain closed and retain seeds until fire causes the cones to open (cone serotiny) represent a key adaptive trait in a variety of pine species. In lodgepole pine, there is substantial geographical variation in serotiny across the Rocky Mountain region. This variation in serotiny has evolved as a result of geographically divergent selection, with consequences that extend to forest communities and ecosystems. An understanding of the genetic architecture of this trait is of interest owing to the wide-reaching ecological consequences of serotiny and also because of the repeated evolution of the trait across the genus. Here, we present and utilize an inexpensive and time-effective method for generating population genomic data. The method uses restriction enzymes and PCR amplification to generate a library of fragments that can be sequenced with a high level of multiplexing. We obtained data for more than 95,000 single nucleotide polymorphisms across 98 serotinous and nonserotinous lodgepole pines from three populations. We used a Bayesian generalized linear model (GLM) to test for an association between genotypic variation at these loci and serotiny. The probability of serotiny varied by genotype at 11 loci, and the association between genotype and serotiny at these loci was consistent in each of the three populations of pines. Genetic variation across these 11 loci explained 50% of the phenotypic variation in serotiny. Our results provide a first genome-wide association map of serotiny in pines and demonstrate an inexpensive and efficient method for generating population genomic data.

Three Strains of Tobacco etch virus Distinctly Alter the Transcriptome of Apical Stem Tissue in Capsicum annuum during Infection.

Tobacco etch virus (TEV; genus Potyvirus) is flexuous rod shaped with a single molecule of single-stranded RNA and causes serious yield losses in species in the Solanaceae. Three TEV strains (HAT, Mex21, and N) are genetically distinct and cause different disease symptoms in plants. Here, a transcriptomic RNA sequencing approach was taken for each TEV strain to evaluate gene expression of the apical stem segment of pepper plants during two stages of disease development. Distinct profiles of Differentially Expressed Genes (DEGs) were identified for each TEV strain. DEG numbers increased with degree of symptom severity: 24 from HAT, 1190 from Mex21, and 4010 from N. At 7 days post-inoculation (dpi), when systemic symptoms were similar, there were few DEGs for HAT- and Mex21-infected plants, whereas N-infected plants had 2516 DEGs. DEG patterns from 7 to 14 dpi corresponded to severity of disease symptoms: milder disease with smaller DEG changes for HAT and Mex21 and severe disease with larger DEG changes for N. Strikingly, in each of these comparisons, there are very few overlapping DEGs among the TEV strains, including no overlapping DEGs between all three strains at 7 or 14 dpi.

Identification of regenerative processes in neonatal spinal cord injury in the opossum (Monodelphis domestica): A transcriptomic study.

This study investigates the response to spinal cord injury in the gray short-tailed opossum (Monodelphis domestica). In opossums spinal injury early in development results in spontaneous axon growth through the injury, but this regenerative potential diminishes with maturity until it is lost entirely. The mechanisms underlying this regeneration remain unknown. RNA sequencing was used to identify differential gene expression in regenerating (SCI at postnatal Day 7, P7SCI) and nonregenerating (SCI at Day 28, P28SCI) cords +1d, +3d, and +7d after complete spinal transection, compared to age-matched controls. Genes showing significant differential expression (log2FC ≥ 1, Padj ≤ 0.05) were used for downstream analysis. Across all time-points 233 genes altered expression after P7SCI, and 472 genes altered expression after P28SCI. One hundred and forty-seven genes altered expression in both injury ages (63% of P7SCI data set). The majority of changes were gene upregulations. Gene ontology overrepresentation analysis in P7SCI gene-sets showed significant overrepresentations only in immune-associated categories, while P28SCI gene-sets showed overrepresentations in these same immune categories, along with other categories such as "cell proliferation," "cell adhesion," and "apoptosis." Cell-type-association analysis suggested that, regardless of injury age, injury-associated gene transcripts were most strongly associated with microglia and endothelial cells, with strikingly fewer astrocyte, oligodendrocyte and neuron-related genes, the notable exception being a cluster of mostly downregulated oligodendrocyte-associated genes in the P7SCI + 7d gene-set. Our findings demonstrate a more complex transcriptomic response in nonregenerating cords, suggesting a strong influence of non-neuronal cells in the outcome after injury and providing the largest survey yet of the transcriptomic changes occurring after SCI in this model.

Transcriptional changes and preservation of bone mass in hibernating black bears.

Physical inactivity leads to losses of bone mass and strength in most mammalian species. In contrast, hibernating bears show no bone loss over the prolonged periods (4-6 months) of immobility during winter, which suggests that they have adaptive mechanisms to preserve bone mass. To identify transcriptional changes that underlie molecular mechanisms preventing disuse osteoporosis, we conducted a large-scale gene expression screening in the trabecular bone and bone marrow, comparing hibernating and summer active bears through sequencing of the transcriptome. Gene set enrichment analysis showed a coordinated down-regulation of genes involved in bone resorption, osteoclast differentiation and signaling, and apoptosis during hibernation. These findings are consistent with previous histological findings and likely contribute to the preservation of bone during the immobility of hibernation. In contrast, no significant enrichment indicating directional changes in gene expression was detected in the gene sets of bone formation and osteoblast signaling in hibernating bears. Additionally, we revealed significant and coordinated transcriptional induction of gene sets involved in aerobic energy production including fatty acid beta oxidation, tricarboxylic acid cycle, oxidative phosphorylation, and mitochondrial metabolism. Mitochondrial oxidation was likely up-regulated by transcriptionally induced AMPK/PGC1α pathway, an upstream stimulator of mitochondrial function.

Mfd Affects Global Transcription and the Physiology of Stressed Bacillus subtilis Cells.

For several decades, Mfd has been studied as the bacterial transcription-coupled repair factor. However, recent observations indicate that this factor influences cell functions beyond DNA repair. Our lab recently described a role for Mfd in disulfide stress that was independent of its function in nucleotide excision repair and base excision repair. Because reports showed that Mfd influenced transcription of single genes, we investigated the global differences in transcription in wild-type and mfd mutant growth-limited cells in the presence and absence of diamide. Surprisingly, we found 1,997 genes differentially expressed in Mfd- cells in the absence of diamide. Using gene knockouts, we investigated the effect of genetic interactions between Mfd and the genes in its regulon on the response to disulfide stress. Interestingly, we found that Mfd interactions were complex and identified additive, epistatic, and suppressor effects in the response to disulfide stress. Pathway enrichment analysis of our RNASeq assay indicated that major biological functions, including translation, endospore formation, pyrimidine metabolism, and motility, were affected by the loss of Mfd. Further, our RNASeq findings correlated with phenotypic changes in growth in minimal media, motility, and sensitivity to antibiotics that target the cell envelope, transcription, and DNA replication. Our results suggest that Mfd has profound effects on the modulation of the transcriptome and on bacterial physiology, particularly in cells experiencing nutritional and oxidative stress.

Endless Forms: Within-Host Variation in the Structure of the West Nile Virus RNA Genome during Serial Passage in Bird Hosts.

RNA viruses are infamous for their high rates of mutation, which produce swarms of genetic variants within individual hosts. To date, analyses of intrahost genetic diversity have focused on the primary genome sequence. However, virus phenotypes are shaped not only by primary sequence but also by the secondary structures into which this sequence folds. Such structures enable viral replication, translation, and binding of small RNAs, yet within-host variation at the structural level has not been adequately explored. We characterized the structural diversity of the 5' untranslated region (UTR) of populations of West Nile virus (WNV) that had been subject to five serial passages in triplicate in each of three bird species. Viral genomes were sampled from host serum samples at each passage (n = 45 populations) and subjected to next-generation sequencing. For populations derived from passages 1, 3, and 5 (n = 9 populations), we predicted the impact of each mutation occurring at a frequency of ≥1% on the secondary structure of the 5' UTR. As expected, mutations in double-stranded (DS) regions of the 5' UTR stem structures caused structural changes of significantly greater magnitude than did mutations in single-stranded (SS) regions. Despite the greater impact of mutations in DS regions, mutations in DS and SS regions occurred at similar frequencies, with no evidence of enhanced selection against mutation in DS regions. In contrast, mutations in two regions that mediate genome cyclization and thereby regulate replication and translation, the 5' cyclization sequence and the UAR flanking stem (UFS), were suppressed in all three hosts.IMPORTANCE The enzymes that copy RNA genomes lack proofreading, and viruses that possess RNA genomes, such as West Nile virus, rapidly diversify into swarms of mutant lineages within a host. Intrahost variation of the primary genomic sequence of RNA viruses has been studied extensively because the extent of this variation shapes key virus phenotypes. However, RNA genomes also form complex secondary structures based on within-genome nucleotide complementarity, which are critical regulators of the cyclization of the virus genome that is necessary for efficient replication and translation. We sought to characterize variation in these secondary structures within populations of West Nile virus during serial passage in three bird species. Our study indicates that the intrahost population of West Nile virus is a diverse assortment of RNA secondary structures that should be considered in future analyses of intrahost viral diversity, but some regions that are critical for genome cyclization are conserved within hosts. Besides potential impacts on viral replication, structural diversity can influence the efficacy of small RNA antiviral therapies.

Transcriptome analysis reveals a stress response of Shewanella oneidensis deprived of background levels of ionizing radiation.

Natural ionizing background radiation has exerted a constant pressure on organisms since the first forms of life appeared on Earth, so that cells have developed molecular mechanisms to avoid or repair damages caused directly by radiation or indirectly by radiation-induced reactive oxygen species (ROS). In the present study, we investigated the transcriptional effect of depriving Shewanella oneidensis cultures of background levels of radiation by growing the cells in a mine 655 m underground, thus reducing the dose rate from 72.1 to 0.9 nGy h-1 from control to treatment, respectively. RNASeq transcriptome analysis showed the differential expression of 4.6 and 7.6% of the S. oneidensis genome during early- and late-exponential phases of growth, respectively. The greatest change observed in the treatment was the downregulation of ribosomal proteins (21% of all annotated ribosomal protein genes during early- and 14% during late-exponential) and tRNA genes (14% of all annotated tRNA genes in early-exponential), indicating a marked decrease in protein translation. Other significant changes were the upregulation of membrane transporters, implying an increase in the traffic of substrates across the cell membrane, as well as the up and downregulation of genes related to respiration, which could be interpreted as a response to insufficient oxidants in the cells. In other reports, there is evidence in multiple species that some ROS not just lead to oxidative stress, but act as signaling molecules to control cellular metabolism at the transcriptional level. Consistent with these reports, several genes involved in the metabolism of carbon and biosynthesis of amino acids were also regulated, lending support to the idea of a wide metabolic response. Our results indicate that S. oneidensis is sensitive to the withdrawal of background levels of ionizing radiation and suggest that a transcriptional response is required to maintain homeostasis and retain normal growth.

Salt and oxidative stresses uniquely regulate tomato cytokinin levels and transcriptomic response.

Cytokinins are well-known to be involved in processes responsible for plant growth and development. More recently, these hormones have begun to be associated with stress responses as well. However, it is unclear how changes in cytokinin biosynthesis, signaling, or transport relate to stress effects. This study examines in parallel how two different stresses, salt, and oxidative stress, affect changes in both cytokinin levels and whole plant transcriptome response. Solanum lycopersicum seedlings were given a short-term (6 hr) exposure to either salt (150 mM NaCl) or oxidative (20 mM H2O2) stress and then examined to determine both changes in cytokinin levels and transcriptome. LC-MS/MS was used to determine the levels of 22 different types of cytokinins in tomato plants including precursors, active, transported, and conjugated forms. When examining cytokinin levels we found that salt treatment caused an increase in both active and inactive cytokinin levels and oxidative stress caused a decrease in these levels. RNA-sequencing analyses of these same stress-treated tissues revealed 6,643 significantly differentially expressed genes (DEGs). Although many DEGs are similar between the two stresses, approximately one-third of the DEGs in each treatment were unique to that stress. Several cytokinin-related genes were among the DEGs. Examination of photosystem II efficiency revealed that cytokinins affect physiological response to stress in tomato, further validating the changes in cytokinin levels seen in planta.

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

[This corrects the article DOI: 10.1371/journal.pone.0171345.].

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.

Histone Citrullination Represses MicroRNA Expression, Resulting in Increased Oncogene mRNAs in Somatolactotrope Cells.

Peptidylarginine deiminase (PAD) enzymes convert histone arginine residues into citrulline to modulate chromatin organization and gene expression. Although PADs are expressed in anterior pituitary gland cells, their functional role and expression in pituitary adenomas are unknown. To begin to address these issues, we first examined normal human pituitaries and pituitary adenomas and found that PAD2, PAD4, and citrullinated histones are highest in prolactinomas and somatoprolactinomas. In the somatoprolactinoma-derived GH3 cell line, PADs citrullinate histone H3, which is attenuated by a pan-PAD inhibitor. RNA sequencing and chromatin immunoprecipitation (ChIP) studies show that the expression of microRNAs (miRNAs) let-7c-2, 23b, and 29c is suppressed by histone citrullination. Our studies demonstrate that these miRNAs directly target the mRNA of the oncogenes encoding HMGA, insulin-like growth factor 1 (IGF-1), and N-MYC, which are highly implicated in human prolactinoma/somatoprolactinoma pathogenesis. Our results are the first to define a direct role for PAD-catalyzed histone citrullination in miRNA expression, which may underlie the etiology of prolactinoma and somatoprolactinoma tumors through regulation of oncogene expression.

A Whole Genome Assembly of the Horn Fly, Haematobia irritans, and Prediction of Genes with Roles in Metabolism and Sex Determination.

Haematobia irritans, commonly known as the horn fly, is a globally distributed blood-feeding pest of cattle that is responsible for significant economic losses to cattle producers. Chemical insecticides are the primary means for controlling this pest but problems with insecticide resistance have become common in the horn fly. To provide a foundation for identification of genomic loci for insecticide resistance and for discovery of new control technology, we report the sequencing, assembly, and annotation of the horn fly genome. The assembled genome is 1.14 Gb, comprising 76,616 scaffolds with N50 scaffold length of 23 Kb. Using RNA-Seq data, we have predicted 34,413 gene models of which 19,185 have been assigned functional annotations. Comparative genomics analysis with the Dipteran flies Musca domestica L., Drosophila melanogaster, and Lucilia cuprina, show that the horn fly is most closely related to M. domestica, sharing 8,748 orthologous clusters followed by D. melanogaster and L. cuprina, sharing 7,582 and 7,490 orthologous clusters respectively. We also identified a gene locus for the sodium channel protein in which mutations have been previously reported that confers target site resistance to the most common class of pesticides used in fly control. Additionally, we identified 276 genomic loci encoding members of metabolic enzyme gene families such as cytochrome P450s, esterases and glutathione S-transferases, and several genes orthologous to sex determination pathway genes in other Dipteran species.

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.

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.

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.

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.

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.

Crude oil impairs immune function and increases susceptibility to pathogenic bacteria in southern flounder.

Exposure to crude oil or its individual constituents can have detrimental impacts on fish species, including impairment of the immune response. Increased observations of skin lesions in northern Gulf of Mexico fish during the 2010 Deepwater Horizon oil spill indicated the possibility of oil-induced immunocompromisation resulting in bacterial or viral infection. This study used a full factorial design of oil exposure and bacterial challenge to examine how oil exposure impairs southern flounder (Paralichthys lethostigma) immune function and increases susceptibility to the bacteria Vibrio anguillarum, a causative agent of vibriosis. Fish exposed to oil prior to bacterial challenge exhibited 94.4% mortality within 48 hours of bacterial exposure. Flounder challenged with V. anguillarum without prior oil exposure had <10% mortality. Exposure resulted in taxonomically distinct gill and intestine bacterial communities. Mortality strongly correlated with V. anguillarum levels, where it comprised a significantly higher percentage of the microbiome in Oil/Pathogen challenged fish and was nearly non-existent in the No Oil/Pathogen challenged fish bacterial community. Elevated V. anguillarum levels were a direct result of oil exposure-induced immunosuppression. Oil-exposure reduced expression of immunoglobulin M, the major systemic fish antibody, and resulted in an overall downregulation in transcriptome response, particularly in genes related to immune function, response to stimulus and hemostasis. Ultimately, sediment-borne oil exposure impairs immune function, leading to increased incidences of bacterial infections. This type of sediment-borne exposure may result in long-term marine ecosystem effects, as oil-bound sediment in the northern Gulf of Mexico will likely remain a contamination source for years to come.