The RNA cargo of Myxococcus outer membrane vesicles.

The outer membrane vesicles (OMVs) secreted by some Gram-negative bacteria contain RNA cargo, which can be introduced into target cells, affecting their cellular processes. To test whether the antimicrobial OMVs secreted by predatory myxobacteria might contain cargo RNA with a role in prey killing, we purified OMVs and cells from four different strains of Myxococcus spp. for RNA-seq transcriptome sequencing. Myxobacterial OMVs contained distinct sets of RNA molecules. The abundance of major cellular transcripts correlated strongly with their abundance in OMVs, suggesting non-specific packaging into OMVs. However, many major cellular transcripts were absent entirely from OMVs and some transcripts were found exclusively in OMVs, suggesting OMV RNA cargo loading is not simply a consequence of sampling the cellular transcriptome. Despite considerable variation in OMV RNA cargo between biological replicates, a small number of transcripts were found consistently in replicate OMV preparations. These 'core' OMV transcripts were often found in the OMVs from multiple strains, and sometimes enriched relative to their abundance in cellular transcriptomes. In addition to providing the first transcriptomes for myxobacterial OMVs, and the first cellular transcriptomes for three strains of Myxococcus spp., we highlight five transcripts for further study. These transcripts are 'core' for at least two of the three strains of M. xanthus studied, and encode two alkyl hydroperoxidase proteins (AhpC and AhpD), two ribosome-associated inhibitors (RaiA-like) and a DO-family protease. It will be interesting to test whether the transcripts serve a biological function within OMVs, potentially being transported into prey cells for translation into toxic proteins.

Interpreting alignment-free sequence comparison: what makes a score a good score?

Alignment-free methods are alternatives to alignment-based methods when searching sequence data sets. The output from an alignment-free sequence comparison is a similarity score, the interpretation of which is not straightforward. We propose objective functions to interpret and calibrate outputs from alignment-free searches, noting that different objective functions are necessary for different biological contexts. This leads to advantages: visualising and comparing score distributions, including those from true positives, may be a relatively simple method to gain insight into the performance of different metrics. Using an empirical approach with both DNA and protein sequences, we characterise different similarity score distributions generated under different parameters. In particular, we demonstrate how sequence length can affect the scores. We show that scores of true positive sequence pairs may correlate significantly with their mean length; and even if the correlation is weak, the relative difference in length of the sequence pair may significantly reduce the effectiveness of alignment-free metrics. Importantly, we show how objective functions can be used with test data to accurately estimate the probability of true positives. This can significantly increase the utility of alignment-free approaches. Finally, we have developed a general-purpose software tool called KAST for use in high-throughput workflows on Linux clusters.

In silico characterisation of the complete Ly6 protein family in Fasciola gigantica supported through transcriptomics of the newly-excysted juveniles.

Fasciola gigantica is one of the aetiological trematodes associated with fascioliasis, which heavily impacts food-production systems and human and animal welfare on a global scale. In the absence of a vaccine, fascioliasis control and treatment is restricted to pasture management, such as clean grazing, and a limited array of chemotherapies, to which signs of resistance are beginning to appear. Research into novel control strategies is therefore urgently required and the advent of 'omics technologies presents considerable opportunity for novel drug and vaccine target discovery. Here, interrogation of the first available F. gigantica newly excysted juvenile (NEJ) transcriptome revealed several protein families of current interest to parasitic flatworm vaccine research, including orthologues of mammalian complement regulator CD59 of the Ly6 family. Ly6 proteins have previously been identified on the tegument of Schistosoma mansoni and induced protective immunity in vaccination trials. Incorporating the recently available F. gigantica genome, the current work revealed 20 novel Ly6 family members in F. gigantica and, in parallel, significantly extended the F. hepatica complement from 3 to 18 members. Phylogenetic analysis revealed several distinct clades within the family, some of which are unique to Fasciola spp. trematodes. Analysis of available proteomic databases also revealed three of the newly discovered FhLy6s were present in extracellular vesicles, which have previously been prioritised in studying the host-parasite interface. The presentation of this new transcriptomic resource, in addition to the Ly6 family proteins here identified, represents a wealth of opportunity for future vaccine research.

Evidence of Immune Modulators in the Secretome of the Equine Tapeworm Anoplocephala perfoliata.

Anoplocephala perfoliata is a neglected gastro-intestinal tapeworm, commonly infecting horses worldwide. Molecular investigation of A. perfoliata is hampered by a lack of tools to better understand the host-parasite interface. This interface is likely influenced by parasite derived immune modulators released in the secretome as free proteins or components of extracellular vesicles (EVs). Therefore, adult RNA was sequenced and de novo assembled to generate the first A. perfoliata transcriptome. In addition, excretory secretory products (ESP) from adult A. perfoliata were collected and EVs isolated using size exclusion chromatography, prior to proteomic analysis of the EVs, the EV surface and EV depleted ESP. Transcriptome analysis revealed 454 sequences homologous to known helminth immune modulators including two novel Sigma class GSTs, five α-HSP90s, and three α-enolases with isoforms of all three observed within the proteomic analysis of the secretome. Furthermore, secretome proteomics identified common helminth proteins across each sample with known EV markers, such as annexins and tetraspanins, observed in EV fractions. Importantly, 49 of the 454 putative immune modulators were identified across the secretome proteomics contained within and on the surface of EVs in addition to those identified in free ESP. This work provides the molecular tools for A. perfoliata to reveal key players in the host-parasite interaction within the horse host.

The Genetics of Prey Susceptibility to Myxobacterial Predation: A Review, Including an Investigation into Pseudomonas aeruginosa Mutations Affecting Predation by Myxococcus xanthus.

Bacterial predation is a ubiquitous and fundamental biological process, which influences the community composition of microbial ecosystems. Among the best characterised bacterial predators are the myxobacteria, which include the model organism Myxococcus xanthus. Predation by M. xanthus involves the secretion of antibiotic metabolites and hydrolytic enzymes, which results in the lysis of prey organisms and release of prey nutrients into the extracellular milieu. Due to the generalist nature of this predatory mechanism, M. xanthus has a broad prey range, being able to kill and consume Gram-negative/positive bacteria and fungi. Potential prey organisms have evolved a range of behaviours which protect themselves from attack by predators. In recent years, several investigations have studied the molecular responses of a broad variety of prey organisms to M. xanthus predation. It seems that the diverse mechanisms employed by prey belong to a much smaller number of general "predation resistance" strategies. In this mini-review, we present the current state of knowledge regarding M. xanthus predation, and how prey organisms resist predation. As previous molecular studies of prey susceptibility have focussed on individual genes/metabolites, we have also undertaken a genome-wide screen for genes of Pseudomonas aeruginosa which contribute to its ability to resist predation. P. aeruginosa is a World Health Organisation priority 1 antibiotic-resistant pathogen. It is metabolically versatile and has an array of pathogenic mechanisms, leading to its prevalence as an opportunistic pathogen. Using a library of nearly 5,500 defined transposon insertion mutants, we screened for "prey genes", which when mutated allowed increased predation by a fluorescent strain of M. xanthus. A set of candidate "prey proteins" were identified, which shared common functional roles and whose nature suggested that predation resistance by P. aeruginosa requires an effective metal/oxidative stress system, an intact motility system, and mechanisms for de-toxifying antimicrobial peptides.

A survey of non-coding RNAs in the social and predatory myxobacterium Myxococcus xanthus DK1622.

Prokaryotic ncRNAs are important regulators of gene expression, and can be involved in complex signalling networks. The myxobacteria are model organisms for studies into multicellular development and microbial predation, being particularly renowned for their large genomes and exceptionally sophisticated signalling networks. However, apart from two specific examples, little is known about their regulatory ncRNAs. Here, we integrate bioinformatic predictions and transcriptome sequence data to provide a comprehensive survey of the ncRNAs made by the exemplar myxobacterium M. xanthus DK1622. M. xanthus RNA-seq data from four experimental conditions was interrogated to identify transcripts mapping outside coding sequences and to known ncRNAs. The resulting 37 ncRNAs were clustered on the genome and most (30/37) were conserved across the myxobacteria. A majority of ncRNAs (22/37) were intergenic, while 13 were at least partially antisense to protein-coding genes. Predicted promoter and terminator sequences explained the start/stop sites of 18 ncRNAs. mRNA targets for the ncRNAs were predicted, including plausible candidates for a known regulatory ncRNA. 22 ncRNAs were differentially expressed by nutrient availability and expression of 25 predicted targets was found to correlate strongly with that of their regulatory ncRNAs. Sharing of predicted mRNA targets by multiple ncRNAs suggests that some ncRNAs might regulate each other within signalling networks. This genomic survey of M. xanthus ncRNA biology provides a starting point for further studies of myxobacterial ncRNAs, which are likely to have important functions in these industrially important and sophisticated organisms.

Genome-Wide Identification of Myxobacterial Predation Genes and Demonstration of Formaldehyde Secretion as a Potentially Predation-Resistant Trait of Pseudomonas aeruginosa.

Despite widespread use in human biology, genome-wide association studies (GWAS) of bacteria are few and have, to date, focused primarily on pathogens. Myxobacteria are predatory microbes with large patchwork genomes, with individual strains secreting unique cocktails of predatory proteins and metabolites. We investigated whether a GWAS strategy could be applied to myxobacteria to identify genes associated with predation. Deduced proteomes from 29 myxobacterial genomes (including eight Myxococcus genomes sequenced for this study), were clustered into orthologous groups, and the presence/absence of orthologues assessed in superior and inferior predators of ten prey organisms. 139 'predation genes' were identified as being associated significantly with predation, including some whose annotation suggested a testable predatory mechanism. Formaldehyde dismutase (fdm) was associated with superior predation of Pseudomonas aeruginosa, and predatory activity of a strain lacking fdm could be increased by the exogenous addition of a formaldehyde detoxifying enzyme, suggesting that production of formaldehyde by P. aeruginosa acts as an anti-predation behaviour. This study establishes the utility of bacterial GWAS to investigate microbial processes beyond pathogenesis, giving plausible and verifiable associations between gene presence/absence and predatory phenotype. We propose that the slow growth rate of myxobacteria, coupled with their predatory mechanism of constitutive secretion, has rendered them relatively resistant to genome streamlining. The resultant genome expansion made possible their observed accumulation of prey-specific predatory genes, without requiring them to be selected for by frequent or recent predation on diverse prey, potentially explaining both the large pan-genome and broad prey range of myxobacteria.

Mendel's pea crosses: varieties, traits and statistics.

A controversy arose over Mendel's pea crossing experiments after the statistician R.A. Fisher proposed how these may have been performed and criticised Mendel's interpretation of his data. Here we re-examine Mendel's experiments and investigate Fisher's statistical criticisms of bias. We describe pea varieties available in Mendel's time and show that these could readily provide all the material Mendel needed for his experiments; the characters he chose to follow were clearly described in catalogues at the time. The combination of character states available in these varieties, together with Eichling's report of crosses Mendel performed, suggest that two of his F3 progeny test experiments may have involved the same F2 population, and therefore that these data should not be treated as independent variables in statistical analysis of Mendel's data. A comprehensive re-examination of Mendel's segregation ratios does not support previous suggestions that they differ remarkably from expectation. The χ2 values for his segregation ratios sum to a value close to the expectation and there is no deficiency of extreme segregation ratios. Overall the χ values for Mendel's segregation ratios deviate slightly from the standard normal distribution; this is probably because of the variance associated with phenotypic rather than genotypic ratios and because Mendel excluded some data sets with small numbers of progeny, where he noted the ratios "deviate not insignificantly" from expectation.

Comparative genomic analysis of six Glossina genomes, vectors of African trypanosomes.

BACKGROUND: Tsetse flies (Glossina sp.) are the vectors of human and animal trypanosomiasis throughout sub-Saharan Africa. Tsetse flies are distinguished from other Diptera by unique adaptations, including lactation and the birthing of live young (obligate viviparity), a vertebrate blood-specific diet by both sexes, and obligate bacterial symbiosis. This work describes the comparative analysis of six Glossina genomes representing three sub-genera: Morsitans (G. morsitans morsitans, G. pallidipes, G. austeni), Palpalis (G. palpalis, G. fuscipes), and Fusca (G. brevipalpis) which represent different habitats, host preferences, and vectorial capacity. RESULTS: Genomic analyses validate established evolutionary relationships and sub-genera. Syntenic analysis of Glossina relative to Drosophila melanogaster shows reduced structural conservation across the sex-linked X chromosome. Sex-linked scaffolds show increased rates of female-specific gene expression and lower evolutionary rates relative to autosome associated genes. Tsetse-specific genes are enriched in protease, odorant-binding, and helicase activities. Lactation-associated genes are conserved across all Glossina species while male seminal proteins are rapidly evolving. Olfactory and gustatory genes are reduced across the genus relative to other insects. Vision-associated Rhodopsin genes show conservation of motion detection/tracking functions and variance in the Rhodopsin detecting colors in the blue wavelength ranges. CONCLUSIONS: Expanded genomic discoveries reveal the genetics underlying Glossina biology and provide a rich body of knowledge for basic science and disease control. They also provide insight into the evolutionary biology underlying novel adaptations and are relevant to applied aspects of vector control such as trap design and discovery of novel pest and disease control strategies.

Evolutionary genomics of anthroponosis in Cryptosporidium.

Human cryptosporidiosis is the leading protozoan cause of diarrhoeal mortality worldwide, and a preponderance of infections is caused by Cryptosporidium hominis and C. parvum. Both species consist of several subtypes with distinct geographical distributions and host preferences (that is, generalist zoonotic and specialist anthroponotic subtypes). The evolutionary processes that drive the adaptation to the human host and the population structures of Cryptosporidium remain unknown. In this study, we analyse 21 whole-genome sequences to elucidate the evolution of anthroponosis. We show that Cryptosporidium parvum splits into two subclades and that the specialist anthroponotic subtype IIc-a shares a subset of loci with C. hominis that is undergoing rapid convergent evolution driven by positive selection. C. parvum subtype IIc-a also has an elevated level of insertion and deletion mutations in the peri-telomeric genes, which is also a characteristic of other specialist subtypes. Genetic exchange between Cryptosporidium subtypes plays a prominent role throughout the evolution of the genus. Interestingly, recombinant regions are enriched for positively selected genes and potential virulence factors, which indicates adaptive introgression. Analysis of 467 gp60 sequences collected from locations across the world shows that the population genetic structure differs markedly between the main zoonotic subtype (isolation-by-distance) and the anthroponotic subtype (admixed population structure). We also show that introgression between the four anthroponotic Cryptosporidium subtypes and species included in this study has occurred recently, probably within the past millennium.

The importance of extracellular vesicle purification for downstream analysis: A comparison of differential centrifugation and size exclusion chromatography for helminth pathogens.

BACKGROUND: Robust protocols for the isolation of extracellular vesicles (EVs) from the rest of their excretory-secretory products are necessary for downstream studies and application development. The most widely used purification method of EVs for helminth pathogens is currently differential centrifugation (DC). In contrast, size exclusion chromatography (SEC) has been included in the purification pipeline for EVs from other pathogens, highlighting there is not an agreed research community 'gold standard' for EV isolation. In this case study, Fasciola hepatica from natural populations were cultured in order to collect EVs from culture media and evaluate a SEC or DC approach to pathogen helminth EV purification. METHODOLOGY/PRINCIPAL FINDINGS: Transmission electron and atomic force microscopy demonstrated that EVs prepared by SEC were both smaller in size and less diverse than EV resolved by DC. Protein quantification and Western blotting further demonstrated that SEC purification realised a higher EV purity to free excretory-secretory protein (ESP) yield ratio compared to DC approaches as evident by the reduction of soluble free cathepsin L proteases in SEC EV preparations. Proteomic analysis further highlighted DC contamination from ESP as shown by an increased diversity of protein identifications and unique peptide hits in DC EVs as compared to SEC EVs. In addition, SEC purified EVs contained less tegumental based proteins than DC purified EVs. CONCLUSIONS/SIGNIFICANCE: The data suggests that DC and SEC purification methods do not isolate equivalent EV population profiles and caution should be taken in the choice of EV purification utilised, with certain protocols for DC preparations including more free ES proteins and tegumental artefacts. We propose that SEC methods should be used for EV purification prior to downstream studies.

Direct Sequencing of Cryptosporidium in Stool Samples for Public Health.

The protozoan parasite Cryptosporidium is an important cause of diarrheal disease (cryptosporidiosis) in humans and animals, with significant morbidity and mortality especially in severely immunocompromised people and in young children in low-resource settings. Due to the sexual life cycle of the parasite, transmission is complex. There are no restrictions on sexual recombination between sub-populations, meaning that large-scale genetic recombination may occur within a host, potentially confounding epidemiological analysis. To clarify the relationships between infections in different hosts, it is first necessary to correctly identify species and genotypes, but these differentiations are not made by standard diagnostic tests and more sophisticated molecular methods have been developed. For instance, multilocus genotyping has been utilized to differentiate isolates within the major human pathogens, Cryptosporidium parvum and Cryptosporidium hominis. This has allowed mixed populations with multiple alleles to be identified: recombination events are considered to be the driving force of increased variation and the emergence of new subtypes. As yet, whole genome sequencing (WGS) is having limited impact on public health investigations, due in part to insufficient numbers of oocysts and purity of DNA derived from clinical samples. Moreover, because public health agencies have not prioritized parasites, validation has not been performed on user-friendly data analysis pipelines suitable for public health practitioners. Nonetheless, since the first whole genome assembly in 2004 there are now numerous genomes of human and animal-derived cryptosporidia publically available, spanning nine species. It has also been demonstrated that WGS from very low numbers of oocysts is possible, through the use of amplification procedures. These data and approaches are providing new insights into host-adapted infectivity, the presence and frequency of multiple sub-populations of Cryptosporidium spp. within single clinical samples, and transmission of infection. Analyses show that although whole genome sequences do indeed contain many alleles, they are invariably dominated by a single highly abundant allele. These insights are helping to better understand population structures within hosts, which will be important to develop novel prevention strategies in the fight against cryptosporidiosis.

Transcriptomic analysis of crustacean neuropeptide signaling during the moult cycle in the green shore crab, Carcinus maenas.

BACKGROUND: Ecdysis is an innate behaviour programme by which all arthropods moult their exoskeletons. The complex suite of interacting neuropeptides that orchestrate ecdysis is well studied in insects, but details of the crustacean ecdysis cassette are fragmented and our understanding of this process is comparatively crude, preventing a meaningful evolutionary comparison. To begin to address this issue we identified transcripts coding for neuropeptides and their putative receptors in the central nervous system (CNS) and Y-organs (YO) within the crab, Carcinus maenas, and mapped their expression profiles across accurately defined stages of the moult cycle using RNA-sequencing. We also studied gene expression within the epidermally-derived YO, the only defined role for which is the synthesis of ecdysteroid moulting hormones, to elucidate peptides and G protein-coupled receptors (GPCRs) that might have a function in ecdysis. RESULTS: Transcriptome mining of the CNS transcriptome yielded neuropeptide transcripts representing 47 neuropeptide families and 66 putative GPCRs. Neuropeptide transcripts that were differentially expressed across the moult cycle included carcikinin, crustacean hyperglycemic hormone-2, and crustacean cardioactive peptide, whilst a single putative neuropeptide receptor, proctolin R1, was differentially expressed. Carcikinin mRNA in particular exhibited dramatic increases in expression pre-moult, suggesting a role in ecdysis regulation. Crustacean hyperglycemic hormone-2 mRNA expression was elevated post- and pre-moult whilst that for crustacean cardioactive peptide, which regulates insect ecdysis and plays a role in stereotyped motor activity during crustacean ecdysis, was elevated in pre-moult. In the YO, several putative neuropeptide receptor transcripts were differentially expressed across the moult cycle, as was the mRNA for the neuropeptide, neuroparsin-1. Whilst differential gene expression of putative neuropeptide receptors was expected, the discovery and differential expression of neuropeptide transcripts was surprising. Analysis of GPCR transcript expression between YO and epidermis revealed 11 to be upregulated in the YO and thus are now candidates for peptide control of ecdysis. CONCLUSIONS: The data presented represent a comprehensive survey of the deduced C. maenas neuropeptidome and putative GPCRs. Importantly, we have described the differential expression profiles of these transcripts across accurately staged moult cycles in tissues key to the ecdysis programme. This study provides important avenues for the future exploration of functionality of receptor-ligand pairs in crustaceans.

The anti-fecundity effect of 5-azacytidine (5-AzaC) on Schistosoma mansoni is linked to dis-regulated transcription, translation and stem cell activities.

Uncontrolled host immunological reactions directed against tissue-trapped eggs precipitate a potentially lethal, pathological cascade responsible for schistosomiasis. Blocking schistosome egg production, therefore, presents a strategy for simultaneously reducing immunopathology as well as limiting disease transmission in endemic or emerging areas. We recently demonstrated that the ribonucleoside analogue 5-azacytidine (5-AzaC) inhibited Schistosoma mansoni oviposition, egg maturation and ovarian development. While these anti-fecundity effects were associated with a loss of DNA methylation, other molecular processes affected by 5-AzaC were not examined at the time. By comparing the transcriptomes of 5-AzaC-treated females to controls, we provide evidence that this ribonucleoside analogue also modulates other crucial aspects of schistosome egg-laying biology. For example, S. mansoni gene products associated with amino acid-, carbohydrate-, fatty acid-, nucleotide- and tricarboxylic acid (TCA)- homeostasis are all dysregulated in 5-AzaC treated females. To validate the metabolic pathway most significantly affected by 5-AzaC, amino acid metabolism, nascent protein synthesis was subsequently quantified in adult schistosomes. Here, 5-AzaC inhibited this process by 68% ±16.7% (SEM) in male- and 81% ±4.8% (SEM) in female-schistosomes. Furthermore, the transcriptome data indicated that adult female stem cells were also affected by 5-AzaC. For instance, 40% of transcripts associated with proliferating schistosome cells were significantly down-regulated by 5-AzaC. This finding correlated with a considerable reduction (95%) in the number of 5-ethynyl-2'-deoxyuridine (EdU) positive cells found in 5-AzaC-treated females. In addition to protein coding genes, the effect that 5-AzaC had on repetitive element expression was also assessed. Here, 46 repeats were found differentially transcribed between 5-AzaC-treated and control females with long terminal repeat (LTR) and DNA transposon classes being amongst the most significant. This study demonstrates that the anti-fecundity activity of 5-AzaC affects more than just DNA methylation in schistosome parasites. Further characterisation of these processes may reveal novel targets for schistosomiasis control.

Transcriptional changes when Myxococcus xanthus preys on Escherichia coli suggest myxobacterial predators are constitutively toxic but regulate their feeding.

Predation is a fundamental ecological process, but within most microbial ecosystems the molecular mechanisms of predation remain poorly understood. We investigated transcriptome changes associated with the predation of Escherichia coli by the myxobacterium Myxococcus xanthus using mRNA sequencing. Exposure to pre-killed prey significantly altered expression of 1319 predator genes. However, the transcriptional response to living prey was minimal, with only 12 genes being significantly up-regulated. The genes most induced by prey presence (kdpA and kdpB, members of the kdp regulon) were confirmed by reverse transcriptase quantitative PCR to be regulated by osmotic shock in M. xanthus, suggesting indirect sensing of prey. However, the prey showed extensive transcriptome changes when co-cultured with predator, with 40 % of its genes (1534) showing significant changes in expression. Bacteriolytic M. xanthus culture supernatant and secreted outer membrane vesicles (OMVs) also induced changes in expression of large numbers of prey genes (598 and 461, respectively). Five metabolic pathways were significantly enriched in prey genes up-regulated on exposure to OMVs, supernatant and/or predatory cells, including those for ribosome and lipopolysaccharide production, suggesting that the prey cell wall and protein production are primary targets of the predator's attack. Our data suggest a model of the myxobacterial predatome (genes and proteins associated with predation) in which the predator constitutively produces secretions which disable its prey whilst simultaneously generating a signal that prey is present. That signal then triggers a regulated feeding response in the predator.

Genome Sequencing and Comparative Transcriptomics Provide a Holistic View of 4-Nitrophenol Degradation and Concurrent Fatty Acid Catabolism by Rhodococcus sp. Strain BUPNP1.

Rhodococcus sp.strain BUPNP1 can utilize the priority environmental pollutant 4-nitrophenol (4-NP) as its sole source of carbon and energy. In this study, genome and transcriptome sequencing were used to gain mechanistic insights into 4-NP degradation. The draft BUPNP1 genome is 5.56 Mbp and encodes 4,963 proteins, which are significantly enriched in hypothetical proteins compared to other Rhodococcus sp. A novel 4-NP catabolic 43 gene cluster "nph" was identified that encodes all the genes required for the conversion of 4-NP into acetyl-CoA and succinate, via 4-nitrocatechol. The cluster also encodes pathways for the catabolism of other diverse aromatic compounds. Comparisons between BUPN1 growing on either 4-NP or glucose resulted in significant changes in the expression of many nph cluster genes, and, during 4-NP growth, a loss of lipid inclusions. Moreover, fatty acid degradation/synthesis genes were found within the nph cluster, suggesting fatty acids may be concurrently catabolised with 4-NP. A holistic model for the action of the nph gene cluster is proposed which incorporates genetic architecture, uptake and metabolism of aromatic compounds, enzymatic activities and transcriptional regulation. The model provides testable hypotheses for further biochemical investigations into the genes of the nph cluster, for potential exploitation in bioremediation.

The Biomphalaria glabrata DNA methylation machinery displays spatial tissue expression, is differentially active in distinct snail populations and is modulated by interactions with Schistosoma mansoni.

BACKGROUND: The debilitating human disease schistosomiasis is caused by infection with schistosome parasites that maintain a complex lifecycle alternating between definitive (human) and intermediate (snail) hosts. While much is known about how the definitive host responds to schistosome infection, there is comparably less information available describing the snail's response to infection. METHODOLOGY/PRINCIPLE FINDINGS: Here, using information recently revealed by sequencing of the Biomphalaria glabrata intermediate host genome, we provide evidence that the predicted core snail DNA methylation machinery components are associated with both intra-species reproduction processes and inter-species interactions. Firstly, methyl-CpG binding domain protein (Bgmbd2/3) and DNA methyltransferase 1 (Bgdnmt1) genes are transcriptionally enriched in gonadal compared to somatic tissues with 5-azacytidine (5-AzaC) treatment significantly inhibiting oviposition. Secondly, elevated levels of 5-methyl cytosine (5mC), DNA methyltransferase activity and 5mC binding in pigmented hybrid- compared to inbred (NMRI)- B. glabrata populations indicate a role for the snail's DNA methylation machinery in maintaining hybrid vigour or heterosis. Thirdly, locus-specific detection of 5mC by bisulfite (BS)-PCR revealed 5mC within an exonic region of a housekeeping protein-coding gene (Bg14-3-3), supporting previous in silico predictions and whole genome BS-Seq analysis of this species' genome. Finally, we provide preliminary evidence for parasite-mediated host epigenetic reprogramming in the schistosome/snail system, as demonstrated by the increase in Bgdnmt1 and Bgmbd2/3 transcript abundance following Bge (B. glabrata embryonic cell line) exposure to parasite larval transformation products (LTP). CONCLUSIONS/SIGNIFICANCE: The presence of a functional DNA methylation machinery in B. glabrata as well as the modulation of these gene products in response to schistosome products, suggests a vital role for DNA methylation during snail development/oviposition and parasite interactions. Further deciphering the role of this epigenetic process during Biomphalaria/Schistosoma co-evolutionary biology may reveal key factors associated with disease transmission and, moreover, enable the discovery of novel lifecycle intervention strategies.

Upgrading short-read animal genome assemblies to chromosome level using comparative genomics and a universal probe set.

Most recent initiatives to sequence and assemble new species' genomes de novo fail to achieve the ultimate endpoint to produce contigs, each representing one whole chromosome. Even the best-assembled genomes (using contemporary technologies) consist of subchromosomal-sized scaffolds. To circumvent this problem, we developed a novel approach that combines computational algorithms to merge scaffolds into chromosomal fragments, PCR-based scaffold verification, and physical mapping to chromosomes. Multigenome-alignment-guided probe selection led to the development of a set of universal avian BAC clones that permit rapid anchoring of multiple scaffolds to chromosomes on all avian genomes. As proof of principle, we assembled genomes of the pigeon (Columbia livia) and peregrine falcon (Falco peregrinus) to chromosome levels comparable, in continuity, to avian reference genomes. Both species are of interest for breeding, cultural, food, and/or environmental reasons. Pigeon has a typical avian karyotype (2n = 80), while falcon (2n = 50) is highly rearranged compared to the avian ancestor. By using chromosome breakpoint data, we established that avian interchromosomal breakpoints appear in the regions of low density of conserved noncoding elements (CNEs) and that the chromosomal fission sites are further limited to long CNE "deserts." This corresponds with fission being the rarest type of rearrangement in avian genome evolution. High-throughput multiple hybridization and rapid capture strategies using the current BAC set provide the basis for assembling numerous avian (and possibly other reptilian) species, while the overall strategy for scaffold assembly and mapping provides the basis for an approach that (provided metaphases can be generated) could be applied to any animal genome.

Identification and temporal expression of putative circadian clock transcripts in the amphipod crustacean Talitrus saltator.

BACKGROUND: Talitrus saltator is an amphipod crustacean that inhabits the supralittoral zone on sandy beaches in the Northeast Atlantic and Mediterranean. T. saltator exhibits endogenous locomotor activity rhythms and time-compensated sun and moon orientation, both of which necessitate at least one chronometric mechanism. Whilst their behaviour is well studied, currently there are no descriptions of the underlying molecular components of a biological clock in this animal, and very few in other crustacean species. METHODS: We harvested brain tissue from animals expressing robust circadian activity rhythms and used homology cloning and Illumina RNAseq approaches to sequence and identify the core circadian clock and clock-related genes in these samples. We assessed the temporal expression of these genes in time-course samples from rhythmic animals using RNAseq. RESULTS: We identified a comprehensive suite of circadian clock gene homologues in T. saltator including the 'core' clock genes period (Talper), cryptochrome 2 (Talcry2), timeless (Taltim), clock (Talclk), and bmal1 (Talbmal1). In addition we describe the sequence and putative structures of 23 clock-associated genes including two unusual, extended isoforms of pigment dispersing hormone (Talpdh). We examined time-course RNAseq expression data, derived from tissues harvested from behaviourally rhythmic animals, to reveal rhythmic expression of these genes with approximately circadian period in Talper and Talbmal1. Of the clock-related genes, casein kinase IIß (TalckIIß), ebony (Talebony), jetlag (Taljetlag), pigment dispensing hormone (Talpdh), protein phosphatase 1 (Talpp1), shaggy (Talshaggy), sirt1 (Talsirt1), sirt7 (Talsirt7) and supernumerary limbs (Talslimb) show temporal changes in expression. DISCUSSION: We report the sequences of principle genes that comprise the circadian clock of T. saltator and highlight the conserved structural and functional domains of their deduced cognate proteins. Our sequencing data contribute to the growing inventory of described comparative clocks. Expression profiling of the identified clock genes illuminates tantalising targets for experimental manipulation to elucidate the molecular and cellular control of clock-driven phenotypes in this crustacean.