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.

An eco-epidemiological modeling approach to investigate dilution effect in two different tick-borne pathosystems.

Disease (re)emergence appears to be driven by biodiversity decline and environmental change. As a result, it is increasingly important to study host-pathogen interactions within the context of their ecology and evolution. The dilution effect is the concept that higher biodiversity decreases pathogen transmission. It has been observed especially in zoonotic vector-borne pathosystems, yet evidence against it has been found. In particular, it is still debated how the community (dis)assembly assumptions and the degree of generalism of vectors and pathogens affect the direction of the biodiversity-pathogen transmission relationship. The aim of this study was to use empirical data and mechanistic models to investigate dilution mechanisms in two rodent-tick-pathogen systems differing in their vector degree of generalism. A community was assembled to include ecological interactions that expand from purely additive to purely substitutive. Such systems are excellent candidates to analyze the link between vector ecology, community (dis)assembly dynamics, and pathogen transmission. To base our mechanistic models on empirical data, rodent live-trapping, including tick sampling, was conducted in Wales across two seasons for three consecutive years. We have developed a deterministic single-vector, multi-host compartmental model that includes ecological relationships with non-host species, uniquely integrating theoretical and observational approaches. To describe pathogen transmission across a gradient of community diversity, the model was populated with parameters describing five different scenarios differing in ecological complexity; each based around one of the pathosystems: Ixodes ricinus (generalist tick)-Borrelia burgdorferi and I. trianguliceps (small mammals specialist tick)-Babesia microti. The results suggested that community composition and interspecific dynamics affected pathogen transmission with different dilution outcomes depending on the vector degree of generalism. The model provides evidence that dilution and amplification effects are not mutually exclusive in the same community but depend on vector ecology and the epidemiological output considered (i.e., the "risk" of interest). In our scenarios, more functionally diverse communities resulted in fewer infectious rodents, supporting the dilution effect. In the pathosystem with generalist vector we identified a hump shaped relationship between diversity and infections in hosts, while for that characterized by specialist tick, this relationship was more complex and more dependent upon specific parameter values.

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.

Genotyping and Whole-Genome Resequencing of Welsh Sheep Breeds Reveal Candidate Genes and Variants for Adaptation to Local Environment and Socioeconomic Traits.

BACKGROUND: Advances in genetic tools applied to livestock breeding has prompted research into the previously neglected breeds adapted to harsh local environments. One such group is the Welsh mountain sheep breeds, which can be farmed at altitudes of 300 m above sea level but are considered to have a low productive value because of their poor wool quality and small carcass size. This is contrary to the lowland breeds which are more suited to wool and meat production qualities, but do not fare well on upland pasture. Herein, medium-density genotyping data from 317 individuals representing 15 Welsh sheep breeds were used alongside the whole-genome resequencing data of 14 breeds from the same set to scan for the signatures of selection and candidate genetic variants using haplotype- and SNP-based approaches. RESULTS: Haplotype-based selection scan performed on the genotyping data pointed to a strong selection in the regions of GBA3, PPARGC1A, APOB, and PPP1R16B genes in the upland breeds, and RNF24, PANK2, and MUC15 in the lowland breeds. SNP-based selection scan performed on the resequencing data pointed to the missense mutations under putative selection relating to a local adaptation in the upland breeds with functions such as angiogenesis (VASH1), anti-oxidation (RWDD1), cell stress (HSPA5), membrane transport (ABCA13 and SLC22A7), and insulin signaling (PTPN1 and GIGFY1). By contrast, genes containing candidate missense mutations in the lowland breeds are related to cell cycle (CDK5RAP2), cell adhesion (CDHR3), and coat color (MC1R). CONCLUSION: We found new variants in genes with potentially functional consequences to the adaptation of local sheep to their environments in Wales. Knowledge of these variations is important for improving the adaptative qualities of UK and world sheep breeds through a marker-assisted selection.

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.

Identifying and validating the presence of Guanine-Quadruplexes (G4) within the blood fluke parasite Schistosoma mansoni.

Schistosomiasis is a neglected tropical disease that currently affects over 250 million individuals worldwide. In the absence of an immunoprophylactic vaccine and the recognition that mono-chemotherapeutic control of schistosomiasis by praziquantel has limitations, new strategies for managing disease burden are urgently needed. A better understanding of schistosome biology could identify previously undocumented areas suitable for the development of novel interventions. Here, for the first time, we detail the presence of G-quadruplexes (G4) and putative quadruplex forming sequences (PQS) within the Schistosoma mansoni genome. We find that G4 are present in both intragenic and intergenic regions of the seven autosomes as well as the sex-defining allosome pair. Amongst intragenic regions, G4 are particularly enriched in 3´ UTR regions. Gene Ontology (GO) term analysis evidenced significant G4 enrichment in the wnt signalling pathway (p<0.05) and PQS oligonucleotides synthetically derived from wnt-related genes resolve into parallel and anti-parallel G4 motifs as elucidated by circular dichroism (CD) spectroscopy. Finally, utilising a single chain anti-G4 antibody called BG4, we confirm the in situ presence of G4 within both adult female and male worm nuclei. These results collectively suggest that G4-targeted compounds could be tested as novel anthelmintic agents and highlights the possibility that G4-stabilizing molecules could be progressed as candidates for the treatment of schistosomiasis.

Genetic and Methylome Variation in Turkish Brachypodium Distachyon Accessions Differentiate Two Geographically Distinct Subpopulations.

Brachypodium distachyon (Brachypodium) is a non-domesticated model grass species that can be used to test if variation in genetic sequence or methylation are linked to environmental differences. To assess this, we collected seeds from 12 sites within five climatically distinct regions of Turkey. Seeds from each region were grown under standardized growth conditions in the UK to preserve methylated sequence variation. At six weeks following germination, leaves were sampled and assessed for genomic and DNA methylation variation. In a follow-up experiment, phenomic approaches were used to describe plant growth and drought responses. Genome sequencing and population structure analysis suggested three ancestral clusters across the Mediterranean, two of which were geographically separated in Turkey into coastal and central subpopulations. Phenotypic analyses showed that the coastal subpopulation tended to exhibit relatively delayed flowering and the central, increased drought tolerance as indicated by reduced yellowing. Genome-wide methylation analyses in GpC, CHG and CHH contexts also showed variation which aligned with the separation into coastal and central subpopulations. The climate niche modelling of both subpopulations showed a significant influence from the "Precipitation in the Driest Quarter" on the central subpopulation and "Temperature of the Coldest Month" on the coastal subpopulation. Our work demonstrates genetic diversity and variation in DNA methylation in Turkish accessions of Brachypodium that may be associated with climate variables and the molecular basis of which will feature in ongoing analyses.

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.

Transcriptional and Metabolomic Analyses Indicate that Cell Wall Properties are Associated with Drought Tolerance in Brachypodium distachyon.

Brachypodium distachyon is an established model for drought tolerance. We previously identified accessions exhibiting high tolerance, susceptibility and intermediate tolerance to drought; respectively, ABR8, KOZ1 and ABR4. Transcriptomics and metabolomic approaches were used to define tolerance mechanisms. Transcriptional analyses suggested relatively few drought responsive genes in ABR8 compared to KOZ1. Linking these to gene ontology (GO) terms indicated enrichment for "regulated stress response", "plant cell wall" and "oxidative stress" associated genes. Further, tolerance correlated with pre-existing differences in cell wall-associated gene expression including glycoside hydrolases, pectin methylesterases, expansins and a pectin acetylesterase. Metabolomic assessments of the same samples also indicated few significant changes in ABR8 with drought. Instead, pre-existing differences in the cell wall-associated metabolites correlated with drought tolerance. Although other features, e.g., jasmonate signaling were suggested in our study, cell wall-focused events appeared to be predominant. Our data suggests two different modes through which the cell wall could confer drought tolerance: (i) An active response mode linked to stress induced changes in cell wall features, and (ii) an intrinsic mode where innate differences in cell wall composition and architecture are important. Both modes seem to contribute to ABR8 drought tolerance. Identification of the exact mechanisms through which the cell wall confers drought tolerance will be important in order to inform development of drought tolerant crops.

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.

Polyomic tools for an emerging livestock parasite, the rumen fluke Calicophoron daubneyi; identifying shifts in rumen functionality.

BACKGROUND: Diseases caused by parasitic flatworms of rumen tissues (paramphistomosis) are a significant threat to global food security as a cause of morbidity and mortality in ruminant livestock in subtropical and tropical climates. Calicophoron daubneyi is currently the only paramphistome species commonly infecting ruminant livestock in temperate European climates. However, recorded incidences of C. daubneyi infection in European livestock have been increasing over the last decade. Whilst clinical paramphistomosis caused by adult worms has not been confirmed in Europe, fatalities have been attributed to severe haemorrhagic enteritis of the small intestine resulting from the migration of immature paramphistomes. Large numbers of mature adults can reside in the rumen, yet to date, the impact on rumen fermentation, and consequently on productivity and economic management of infected livestock, have not been resolved. Limited publicly available nucleotide and protein sequences for C. daubneyi underpin this lack of biological and economic understanding. Here we present for the first time a de novo assembled transcriptome, with functional annotations, for adult C. daubneyi, which provides a reference database for protein and nucleotide sequence identification to facilitate fundamental biology, anthelmintic, vaccine and diagnostics discoveries. RESULTS: This dataset identifies a number of genes potentially unique to C. daubneyi and, by comparison to an existing transcriptome for the related Paramphistomum cervi, identifies novel genes which may be unique to the paramphistome group of platyhelminthes. Additionally, we present the first coverage of the excretory/secretory and soluble somatic proteome profiles for adult C. daubneyi and identify the release of extracellular vesicles from adult C. daubneyi parasites during in vitro, ex-host culture. Finally, we have performed the first analysis of rumen fluke impacting upon rumen fermentation parameters using an in vitro gas production study resulting in a significant increase in propionate production. CONCLUSIONS: The resulting data provide a discovery platform (transcriptome, proteomes, EV isolation pipeline and in vitro fermentation system) to further study C. daubneyi-host interaction. In addition, the acetate: propionate ratio has been demonstrated to decrease with rumen fluke infection suggesting that acidotic conditions in the rumen may occur.

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.

Identification of Stipules reduced, a leaf morphology gene in pea (Pisum sativum).

Pea (Pisum sativum) is one of relatively few genetically amenable plant species with compound leaves. Pea leaves have a variety of specialized organs: leaflets, tendrils, pulvini and stipules, which enable the identification of mutations that transform or affect distinct parts of the leaf. Characterization of these mutations offers insights into the development and evolution of novel leaf traits. The previously characterized morphological gene Cochleata, conferring stipule identity, was known to interact with Stipules reduced (St), which conditions stipule size in pea, but the St gene remained unknown. Here we analysed Fast Neutron irradiated pea mutants by restriction site associated DNA sequencing. We identified St as a gene encoding a C2H2 zinc finger transcription factor that is regulated by Cochleata. St regulates both cell division and cell expansion in the stipule. Our approach shows how systematic genome-wide screens can be used successfully for the analysis of traits in species for which whole genome sequences are not available.