Transcriptomic Responses to Koi Herpesvirus in Isolated Blood Leukocytes from Infected Common Carp.

Koi herpesvirus (KHV, CyHV-3) causes severe economic losses in carp farms. Its eradication is challenging due to the establishment of latency in blood leukocytes and other tissues. To understand the molecular mechanisms leading to KHV infection in leukocytes, common carp were bath-exposed to KHV at 17 °C. After confirming the presence of viral transcripts in blood leukocytes at ten days post infection, RNA-Seq was performed on peripheral blood leukocytes on the Illumina NovaSeq. KHV infection triggered a robust immune response mediated by pattern recognition receptors, mainly toll-like receptors (tlr2, tlr5, tlr7, and tlr13), urokinase plasminogen activator surface receptor-like, galectin proteins, and lipid mediators such as leukotriene B4 receptor 1. Enriched pathways showed increased mitochondria oxidative phosphorylation and the activation of signalling pathways such as mitogen-activated protein kinases (MAPKs) and vascular endothelial growth factor (VEGF). KHV-infected leukocytes showed low production of reactive oxygen species (ROS) and glutathione metabolism, high iron export and phagocytosis activity, and low autophagy. Macrophage polarization was deduced from the up-regulation of genes such as arginase non-hepatic 1-like, macrophage mannose receptor-1, crem, il-10, and il-13 receptors, while markers for cytotoxic T cells were observed to be down-regulated. Further work is required to characterise these leukocyte subsets and the molecular events leading to KHV latency in blood leukocytes.

Resistance to white spot syndrome virus in the European shore crab is associated with suppressed virion trafficking and heightened immune responses.

Introduction: All decapod crustaceans are considered potentially susceptible to White Spot Syndrome Virus (WSSV) infection, but the degree of White Spot Disease (WSD) susceptibility varies widely between species. The European shore crab Carcinus maenas can be infected with the virus for long periods of time without signs of disease. Given the high mortality rate of susceptible species, the differential susceptibility of these resistant hosts offers an opportunity to investigate mechanisms of disease resistance. Methods: Here, the temporal transcriptional responses (mRNA and miRNA) of C. maenas following WSSV injection were analysed and compared to a previously published dataset for the highly WSSV susceptible Penaeus vannamei to identify key genes, processes and pathways contributing to increased WSD resistance. Results: We show that, in contrast to P. vannamei, the transcriptional response during the first 2 days following WSSV injection in C. maenas is limited. During the later time points (7 days onwards), two groups of crabs were identified, a recalcitrant group where no replication of the virus occurred, and a group where significant viral replication occurred, with the transcriptional profiles of the latter group resembling those of WSSV-susceptible species. We identify key differences in the molecular responses of these groups to WSSV injection. Discussion: We propose that increased WSD resistance in C. maenas may result from impaired WSSV endocytosis due to the inhibition of internal vesicle budding by dynamin-1, and a delay in movement to the nucleus caused by the downregulation of cytoskeletal transcripts required for WSSV cytoskeleton docking, during early stages of the infection. This response allows resistant hosts greater time to fine-tune immune responses associated with miRNA expression, apoptosis and the melanisation cascade to defend against, and clear, invading WSSV. These findings suggest that the initial stages of infection are key to resistance to WSSV in the crab and highlight possible pathways that could be targeted in farmed crustacean to enhance resistance to WSD.

The complete mitochondrial genome of the pink sea fan, Eunicella verrucosa (Pallas, 1766).

The pink sea fan, Eunicella verrucosa (Pallas, 1766), inhabits rocky substrates across the northeast Atlantic and the western Mediterranean. Across much of its range it has been detrimentally affected by fishing. DNA from 17 E. verrucosa specimens was amplified by phi29-induced rolling circle amplification. Following purification by sodium acetate-ethanol precipitation, the circular genomic DNA was sequenced on an Illumina MiSeq v2. Specimens originated from sites along the west coast of Ireland, southwest Wales, southwest/southern England, northwest France, southern Portugal, and the Mediterranean coast of northeast Spain. All samples had identical mitochondrial genome sequences of 19,267 bp and included 14 protein-coding genes (including the mutS gene), two ribosomal RNA subunits (12S and 16S) and one methionine tRNA gene. Two genes (nad2 and nad5) overlapped by 13 bp; all other genes were separated by non-coding intergenic regions. All protein-coding genes had the same start codon (ATG) and a TAA or TAG stop codon, except for cox1 that terminated with the incomplete stop codon T--. The mitochondrial genome of E. verrucosa (MW588805) showed 99.72% similarity with that of a related sea fan species, Eunicella cavolini, with six SNPs and a 49 bp deletion between nad5 and nad4 in E. verrucosa distinguishing the two.

Transcriptomic analysis of the trade-off between endurance and burst-performance in the frog Xenopus allofraseri.

BACKGROUND: Variation in locomotor capacity among animals often reflects adaptations to different environments. Despite evidence that physical performance is heritable, the molecular basis of locomotor performance and performance trade-offs remains poorly understood. In this study we identify the genes, signaling pathways, and regulatory processes possibly responsible for the trade-off between burst performance and endurance observed in Xenopus allofraseri, using a transcriptomic approach. RESULTS: We obtained a total of about 121 million paired-end reads from Illumina RNA sequencing and analyzed 218,541 transcripts obtained from a de novo assembly. We identified 109 transcripts with a significant differential expression between endurant and burst performant individuals (FDR ≤ 0.05 and logFC ≥2), and blast searches resulted in 103 protein-coding genes. We found major differences between endurant and burst-performant individuals in the expression of genes involved in the polymerization and ATPase activity of actin filaments, cellular trafficking, proteoglycans and extracellular proteins secreted, lipid metabolism, mitochondrial activity and regulators of signaling cascades. Remarkably, we revealed transcript isoforms of key genes with functions in metabolism, apoptosis, nuclear export and as a transcriptional corepressor, expressed in either burst-performant or endurant individuals. Lastly, we find two up-regulated transcripts in burst-performant individuals that correspond to the expression of myosin-binding protein C fast-type (mybpc2). This suggests the presence of mybpc2 homoeologs and may have been favored by selection to permit fast and powerful locomotion. CONCLUSION: These results suggest that the differential expression of genes belonging to the pathways of calcium signaling, endoplasmic reticulum stress responses and striated muscle contraction, in addition to the use of alternative splicing and effectors of cellular activity underlie locomotor performance trade-offs. Ultimately, our transcriptomic analysis offers new perspectives for future analyses of the role of single nucleotide variants, homoeology and alternative splicing in the evolution of locomotor performance trade-offs.

Feral populations of Brassica oleracea along Atlantic coasts in western Europe.

There has been growing emphasis on the role that crop wild relatives might play in supporting highly selected agriculturally valuable species in the face of climate change. In species that were domesticated many thousands of years ago, distinguishing wild populations from escaped feral forms can be challenging, but reintroducing variation from either source could supplement current cultivated forms. For economically important cabbages (Brassicaceae: Brassica oleracea), "wild" populations occur throughout Europe but little is known about their genetic variation or potential as resources for breeding more resilient crop varieties. The main aim of this study was to characterize the population structure of geographically isolated wild cabbage populations along the coasts of the UK and Spain, including the Atlantic range edges. Double-digest restriction-site-associated DNA sequencing was used to sample individual cabbage genomes, assess the similarity of plants from 20 populations, and explore environment-genotype associations across varying climatic conditions. Interestingly, there were no indications of isolation by distance; several geographically close populations were genetically more distinct from each other than to distant populations. Furthermore, several distant populations shared genetic ancestry, which could indicate that they were established by escapees of similar source cultivars. However, there were signals of local adaptation to different environments, including a possible relationship between genetic diversity and soil pH. Overall, these results highlight wild cabbages in the Atlantic region as an important genetic resource worthy of further research into their relationship with existing crop varieties.

Common Host Responses in Murine Aerosol Models of Infection Caused by Highly Virulent Gram-Negative Bacteria from the Genera Burkholderia, Francisella and Yersinia.

Highly virulent bacterial pathogens cause acute infections which are exceptionally difficult to treat with conventional antibiotic therapies alone. Understanding the chain of events that are triggered during an infection of a host has the potential to lead to new therapeutic strategies. For the first time, the transcriptomic responses within the lungs of Balb/C mice have been compared during an acute infection with the intracellular pathogens Burkholderia pseudomallei, Francisella tularensis and Yersinia pestis. Temporal changes were determined using RNAseq and a bioinformatics pipeline; expression of protein was also studied from the same sample. Collectively it was found that early transcriptomic responses within the infected host were associated with the (a) slowing down of critical cellular functions, (b) production of circulatory system components, (c) lung tissue integrity, and (d) intracellular regulatory processes. One common molecule was identified, Errfi1 (ErbB receptor feedback inhibitor 1); upregulated in response to all three pathogens and a potential novel marker of acute infection. Based upon the pro-inflammatory responses observed, we sought to synchronise each infection and report that 24 h p.i. of B. pseudomallei infection closely aligned with 48 h p.i. of infection with F. tularensis and Y. pestis. Post-transcriptional modulation of RANTES expression occurred across all pathogens, suggesting that these infections directly or indirectly modulate cell trafficking through chemokine expression/detection. Collectively, this unbiased NGS approach has provided an in-depth characterisation of the host transcriptome following infection with these highly virulent pathogens ultimately aiding in the development of host-directed therapies as adjuncts or alternatives to antibiotic treatment.

Palaeogenomics of the Hydrocarbon Producing Microalga Botryococcus braunii.

Botryococcus braunii is a colonial microalga that appears early in the fossil record and is a sensitive proxy of environmental and hydroclimatic conditions. Palaeozoic Botryococcus fossils which contribute up to 90% of oil shales and approximately 1% of crude oil, co-localise with diagnostic geolipids from the degradation of source-signature hydrocarbons. However more recent Holocene sediments demonstrate no such association. Consequently, Botryococcus are identified in younger sediments by morphology alone, where potential misclassifications could lead to inaccurate paleoenvironmental reconstructions. Here we show that a combination of flow cytometry and ancient DNA (aDNA) sequencing can unambiguously identify Botryococcus microfossils in Holocene sediments with hitherto unparalleled accuracy and rapidity. The application of aDNA sequencing to microfossils offers a far-reaching opportunity for understanding environmental change in the recent geological record. When allied with other high-resolution palaeoenvironmental information such as aDNA sequencing of humans and megafauna, aDNA from microfossils may allow a deeper and more precise understanding of past environments, ecologies and migrations.

Molecular epidemiology of Pseudomonas aeruginosa in an unsegregated bronchiectasis cohort sharing hospital facilities with a cystic fibrosis cohort.

While Pseudomonas aeruginosa (PA) cross-infection is well documented among patients with cystic fibrosis (CF), the equivalent risk among patients with non-CF bronchiectasis (NCFB) is unclear, particularly those managed alongside patients with CF. We performed analysis of PA within a single centre that manages an unsegregated NCFB cohort alongside a segregated CF cohort. We found no evidence of cross-infection between the two cohorts or within the segregated CF cohort. However, within the unsegregated NCFB cohort, evidence of cross-infection was found between three (of 46) patients. While we do not presently advocate any change in the management of our NCFB cohort, longitudinal surveillance is clearly warranted.

Decay of the glycolytic pathway and adaptation to intranuclear parasitism within Enterocytozoonidae microsporidia.

Glycolysis and oxidative phosphorylation are the fundamental pathways of ATP generation in eukaryotes. Yet in microsporidia, endoparasitic fungi living at the limits of cellular streamlining, oxidative phosphorylation has been lost: energy is obtained directly from the host or, during the dispersive spore stage, via glycolysis. It was therefore surprising when the first sequenced genome from the Enterocytozoonidae - a major family of human and animal-infecting microsporidians - appeared to have lost genes for glycolysis. Here, we sequence and analyse genomes from additional members of this family, shedding new light on their unusual biology. Our survey includes the genome of Enterocytozoon hepatopenaei, a major aquacultural parasite currently causing substantial economic losses in shrimp farming, and Enterospora canceri, a pathogen that lives exclusively inside epithelial cell nuclei of its crab host. Our analysis of gene content across the clade suggests that Ent. canceri's adaptation to intranuclear life is underpinned by the expansion of transporter families. We demonstrate that this entire lineage of pathogens has lost glycolysis and, uniquely amongst eukaryotes, lacks any obvious intrinsic means of generating energy. Our study provides an important resource for the investigation of host-pathogen interactions and reductive evolution in one of the most medically and economically important microsporidian lineages.

Metagenomic Analysis of Silage.

Metagenomics is defined as the direct analysis of deoxyribonucleic acid (DNA) purified from environmental samples and enables taxonomic identification of the microbial communities present within them. Two main metagenomic approaches exist; sequencing the 16S rRNA gene coding region, which exhibits sufficient variation between taxa for identification, and shotgun sequencing, in which genomes of the organisms that are present in the sample are analyzed and ascribed to "operational taxonomic units"; species, genera or families depending on the extent of sequencing coverage. In this study, shotgun sequencing was used to analyze the microbial community present in cattle silage and, coupled with a range of bioinformatics tools to quality check and filter the DNA sequence reads, perform taxonomic classification of the microbial populations present within the sampled silage, and achieve functional annotation of the sequences. These methods were employed to identify potentially harmful bacteria that existed within the silage, an indication of silage spoilage. If spoiled silage is not remediated, then upon ingestion it could be potentially fatal to the livestock.

Pitfalls of haplotype phasing from amplicon-based long-read sequencing.

The long-read sequencers from Pacific Bioscience (PacBio) and Oxford Nanopore Technologies (ONT) offer the opportunity to phase mutations multiple kilobases apart directly from sequencing reads. In this study, we used long-range PCR with ONT and PacBio sequencing to phase two variants 9 kb apart in the RET gene. We also re-analysed data from a recent paper which had apparently successfully used ONT to phase clinically important haplotypes at the CYP2D6 and HLA loci. From these analyses, we demonstrate PCR-chimera formation during PCR amplification and reference alignment bias are pitfalls that need to be considered when attempting to phase variants using amplicon-based long-read sequencing technologies. These methodological pitfalls need to be avoided if the opportunities provided by long-read sequencers are to be fully exploited.

Transcriptomic profiling of host-parasite interactions in the microsporidian Trachipleistophora hominis.

BACKGROUND: Trachipleistophora hominis was isolated from an HIV/AIDS patient and is a member of a highly successful group of obligate intracellular parasites. METHODS: Here we have investigated the evolution of the parasite and the interplay between host and parasite gene expression using transcriptomics of T. hominis-infected rabbit kidney cells. RESULTS: T. hominis has about 30% more genes than small-genome microsporidians. Highly expressed genes include those involved in growth, replication, defence against oxidative stress, and a large fraction of uncharacterised genes. Chaperones are also highly expressed and may buffer the deleterious effects of the large number of non-synonymous mutations observed in essential T. hominis genes. Host expression suggests a general cellular shutdown upon infection, but ATP, amino sugar and nucleotide sugar production appear enhanced, potentially providing the parasite with substrates it cannot make itself. Expression divergence of duplicated genes, including transporters used to acquire host metabolites, demonstrates ongoing functional diversification during microsporidian evolution. We identified overlapping transcription at more than 100 loci in the sparse T. hominis genome, demonstrating that this feature is not caused by genome compaction. The detection of additional transposons of insect origin strongly suggests that the natural host for T. hominis is an insect. CONCLUSIONS: Our results reveal that the evolution of contemporary microsporidian genomes is highly dynamic and innovative. Moreover, highly expressed T. hominis genes of unknown function include a cohort that are shared among all microsporidians, indicating that some strongly conserved features of the biology of these enormously successful parasites remain uncharacterised.

Draft Genome Sequence of Stenotrophomonas maltophilia SeITE02, a Gammaproteobacterium Isolated from Selenite-Contaminated Mining Soil.

Stenotrophomonas maltophilia strain SeITE02 was isolated from the rhizosphere of the selenium-hyperaccumulating legume Astragalus bisculcatus. In this report, we provide the 4.56-Mb draft genome sequence of S. maltophilia SeITE02, a gammaproteobacterium that can withstand high concentrations of selenite and reduce these to elemental selenium.

The stochastic silencing phenotype of Arabidopsis morc6 mutants reveals a role in efficient RNA-directed DNA methylation.

The RNA-directed DNA methylation (RdDM) pathway is of central importance to the initiation and maintenance of transcriptional gene silencing in plants. DNA methylation is directed to target sequences by a mechanism that involves production of small RNAs by RNA polymerase IV and long non-coding RNAs by RNA polymerase V. DNA methylation then leads to recruitment of histone-modifying enzymes, followed by establishment of a silenced chromatin state. Recently MORC6, a member of the microrchidia (MORC) family of adenosine triphosphatases (ATPases), has been shown to be involved in transcriptional gene silencing. However, reports differ regarding whether MORC6 is involved in RdDM itself or acts downstream of DNA methylation to enable formation of higher-order chromatin structure. Here we demonstrate that MORC6 is required for efficient RdDM at some target loci, and, using a GFP reporter system, we found that morc6 mutants show a stochastic silencing phenotype. By using cell sorting to separate silenced and unsilenced cells, we show that release of silencing at this locus is associated with a loss of DNA methylation. Thus our data support a view that MORC6 influences RdDM and that it is not acting downstream of DNA methylation. For some loci, efficient initiation or maintenance of DNA methylation may depend on the ability to form higher-order chromatin structure.