Search

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 natural genetic diversity of agricultural species is an essential genetic resource for breeding programs aiming to improve their ecosystem and production services. A large natural ecotype diversity is usually available for most grassland species. This could be used to recombine natural climatic adaptations and agronomic value to create improved populations of grassland species adapted to future regional climates. However describing natural genetic resources can be long and costly. Molecular markers may provide useful information to help this task. This opportunity was investigated for Lolium perenne L., using a set of 385 accessions from the natural diversity of this species collected right across Europe and provided by genebanks of several countries. For each of these populations, genotyping provided the allele frequencies of 189,781 SNP markers. GWAS were implemented for over 30 agronomic and/or putatively adaptive traits recorded in three climatically contrasted locations (France, Belgium, Germany). Significant associations were detected for hundreds of markers despite a strong confounding effect of the genetic background; most of them pertained to phenology traits. It is likely that genetic variability in these traits has had an important contribution to environmental adaptation and ecotype differentiation. Genomic prediction models calibrated using natural diversity were found to be highly effective to describe natural populations for almost all traits as well as commercial synthetic populations for some important traits such as disease resistance, spring growth or phenological traits. These results will certainly be valuable information to help the use of natural genetic resources of other species., (Copyright © 2020 Keep et al.)

The UK and Ireland have many native pony breeds with historical and cultural importance as well as being a source of uncharacterized genetic diversity. However, there is a lack of comprehensive research investigating their genetic diversity and phylogenetic interrelationships. Many studies contain a limited number of pony breeds or small sample sizes for these breeds. This may result in erroneous grouping of pony breeds that otherwise have intricate interrelationships with each other and are not evaluated correctly when placed as a token subset of a larger dataset. This is the first study that specifically investigates the genetic diversity within and between British and Irish native pony breeds using large sample numbers from locations of their native origin. This study used a panel of microsatellite markers and sequence analysis of the mitochondrial control region to analyze the genetic diversity within and between 11 pony breeds from Britain and Ireland. A large dataset was collected (a total of 485 animals were used for mtDNA analysis and 450 for microsatellite analysis), and previously published data were used to place the British and Irish ponies in a global context. The native ponies of Britain and Ireland were found to have had a complex history, and the interrelationships between the breeds were revealed. Overall, high levels of genetic diversity were maintained in native breeds, although some reduction was evident in small or isolated populations (Shetland, Carneddau, and Section C). Unusual mitochondrial diversity distribution patterns were apparent for the Carneddau and Dartmoor, although among breeds and global haplogroups there was a high degree of haplotype sharing evident, well-represented within British and Irish ponies. Ancestral maternal diversity was maintained by most populations, particularly the Fells and Welsh ponies, which exhibited rare and ancient lineages. The maternal and paternal histories of the breeds are distinct, with male-biased crossings between native breeds, and other shared influences, likely Arabs and Thoroughbreds, are apparent. The data generated herein provide valuable information to guide and implement the conservation of increasingly rare native genetic resources., Competing Interests: There are no conflicts of interest between authors, and permission from each author has been granted. No other institutions hold copyright over this work., (© 2020 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

Understanding the origin and maintenance of phenotypic variation, particularly across a continuous spatial distribution, represents a key challenge in evolutionary biology. For this, animal venoms represent ideal study systems: they are complex, variable, yet easily quantifiable molecular phenotypes with a clear function. Rattlesnakes display tremendous variation in their venom composition, mostly through strongly dichotomous venom strategies, which may even coexist within a single species. Here, through dense, widespread population-level sampling of the Mojave rattlesnake, Crotalus scutulatus, we show that genomic structural variation at multiple loci underlies extreme geographical variation in venom composition, which is maintained despite extensive gene flow. Unexpectedly, neither diet composition nor neutral population structure explain venom variation. Instead, venom divergence is strongly correlated with environmental conditions. Individual toxin genes correlate with distinct environmental factors, suggesting that different selective pressures can act on individual loci independently of their co-expression patterns or genomic proximity. Our results challenge common assumptions about diet composition as the key selective driver of snake venom evolution and emphasize how the interplay between genomic architecture and local-scale spatial heterogeneity in selective pressures may facilitate the retention of adaptive functional polymorphisms across a continuous space.

Bee populations are currently undergoing severe global declines driven by the interactive effects of a number of factors. Ongoing urbanisation has the potential to exacerbate bee declines, unless steps are taken to ensure appropriate floral resources are available. Sown wildflower strips are one way in which floral resources can be provided to urban bees. However, the use of these strips by pollinators in urban environments remains little studied. Here, we employ pollen metabarcoding of the rbcL gene to compare the foraging patterns of different bee species observed using urban sown wildflower strips in July 2016, with a goal of identifying which plant species are most important for bees. We also demonstrate the use of a non-destructive method of pollen collection. Bees were found to forage on a wide variety of plant genera and families, including a diverse range of plants from outside the wildflower plots, suggesting that foragers visiting sown wildflower strips also utilize other urban habitats. Particular plants within the wildflower strips dominated metabarcoding data, particularly Papaver rhoeas and Phacelia tanacetifolia . Overall, we demonstrate that pollinators observed in sown wildflower strips use certain sown foodplants as part of a larger urban matrix., Competing Interests: Natasha de Vere is an Academic Editor for PeerJ.

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.

Height is an important characteristic in the equine industry although little is known about its genetic control in native British breeds of ponies. This study aimed to map QTL data with the withers height in four pony breeds native to the British Isles, including two different sections within Welsh Cobs. In this study, a genome-wide analysis approach using the Illumina EquineSNP50 Infinium BeadChip was applied to 105 ponies and cobs. Analysis identified 222 highly significant height-associated SNPs (P ≤ 10 -5 ), among which three SNPs on ECA9 have also been previously reported elsewhere. The highest number of significant SNPs associated to height in the native British horses were located on ECA1, ECA8, and ECA16.

Three variants of the multidrug-resistant plasmid pLUH01 were assembled by deep sequencing from nasopharyngeal swabs. All have a 21-bp deletion in the RS14515 hypothetical gene. Variants 1 through 3 have 2, 6, and 3 nucleotide substitutions, respectively, compared to the pLUH01 reference genome. We named the new plasmid variants pLUH01/Lancaster/2015/1 to pLUH01/Lancaster/2015/3.

Nasopharyngeal swabs were taken from volunteers attending a general medical practice and a general hospital in Lancaster, UK, and at Lancaster University, in the winter of 2014-2015. 51 swabs were selected based on high RNA yield and allocated to deep sequencing pools as follows: patients with chronic obstructive pulmonary disease; asthmatics; adults with no respiratory symptoms; adults with feverish respiratory symptoms; adults with respiratory symptoms and presence of antibodies against influenza C; paediatric patients with respiratory symptoms (2 pools); adults with influenza C infection (2 pools), giving a total of 9 pools. Illumina sequencing was performed, with data yields per pool in the range of 345.6 megabases to 14 gigabases after removal of reads aligning to the human genome. The data were deposited in the Sequence Read Archive at NCBI, and constitute a resource for study of the viral, bacterial and fungal metagenome of the human nasopharynx in healthy and diseased states and comparison with other metagenomic studies on the human respiratory tract.

The genome sequence of human papillomavirus type 20 (HPV-20; family Papillomaviridae , genus Betapapillomavirus , species Betapapillomavirus 1 , type 20) was assembled by deep sequencing from nasopharyngeal swabs. The assembled genome is 0.37% divergent over its full length from the single complete genome of HPV-20 in GenBank (U31778). We named the strain HPV-20/Lancaster/2015., (Copyright © 2017 Atkinson et al.)

The genome of human papillomavirus type 23 (HPV-23; family Papillomaviridae , genus Betapapillomavirus , species Betapapillomavirus 2 , type 23) was assembled by deep sequencing from nasopharyngeal swabs. The assembled genome is 2.7% divergent over its full length from the single complete genome of HPV-23 in GenBank (accession no. U31781). We named the strain HPV-23/Lancaster/2015., (Copyright © 2017 Atkinson et al.)

Influenza C is not included in the annual seasonal influenza vaccine, and has historically been regarded as a minor respiratory pathogen. However, recent work has highlighted its potential role as a cause of pneumonia in infants. We performed nasopharyngeal or nasal swabbing and/or serum sampling (n = 148) in Lancaster, UK, over the winter of 2014-2015. Using enzyme-linked immunosorbent assay (ELISA), we obtain seropositivity of 77%. By contrast, only 2 individuals, both asymptomatic adults, were influenza C-positive by polymerase chain reaction (PCR). Deep sequencing of nasopharyngeal samples produced partial sequences for 4 genome segments in one of these patients. Bayesian phylogenetic analysis demonstrated that the influenza C genome from this individual is evolutionarily distant to those sampled in recent years and represents a novel genome constellation, indicating that it may be a product of a decades-old reassortment event. Although we find no evidence that influenza C was a significant respiratory pathogen during the winter of 2014-2015 in Lancaster, we confirm previous observations of seropositivity in the majority of the population. (170 words).

The genome of human rhinovirus A22 (HRV-A22) was assembled by deep sequencing RNA samples from nasopharyngeal swabs. The assembled genome is 8.7% divergent from the HRV-A22 reference strain over its full length, and it is only the second full-length genome sequence for HRV-A22. The new strain is designated strain HRV-A22/Lancaster/2015., (Copyright © 2017 Atkinson et al.)

Understanding which flowers honey bees (Apis mellifera) use for forage can help us to provide suitable plants for healthy honey bee colonies. Accordingly, honey DNA metabarcoding provides a valuable tool for investigating pollen and nectar collection. We investigated early season (April and May) floral choice by honey bees provided with a very high diversity of flowering plants within the National Botanic Garden of Wales. There was a close correspondence between the phenology of flowering and the detection of plants within the honey. Within the study area there were 437 genera of plants in flower during April and May, but only 11% of these were used. Thirty-nine plant taxa were recorded from three hives but only ten at greater than 1%. All three colonies used the same core set of native or near-native plants, typically found in hedgerows and woodlands. The major plants were supplemented with a range of horticultural species, with more variation in plant choice between the honey bee colonies. We conclude that during the spring, honey bees need access to native hedgerows and woodlands to provide major plants for foraging. Gardens provide supplementary flowers that may increase the nutritional diversity of the honey bee diet.

Understanding the relationship between ingested plant material and the attached microbiome is essential for developing methodologies to improve ruminant nutrient use efficiency. We have previously shown that perennial ryegrass (PRG) rumen bacterial colonization events follow a primary (up to 4 h) and secondary (after 4 h) pattern based on the differences in diversity of the attached bacteria. In this study, we investigated temporal niche specialization of primary and secondary populations of attached rumen microbiota using metagenomic shotgun sequencing as well as monitoring changes in the plant chemistry using mid-infrared spectroscopy (FT-IR). Metagenomic Rapid Annotation using Subsystem Technology (MG-RAST) taxonomical analysis of shotgun metagenomic sequences showed that the genera Butyrivibrio, Clostridium, Eubacterium, Prevotella , and Selenomonas dominated the attached microbiome irrespective of time. MG-RAST also showed that Acidaminococcus, Bacillus, Butyrivibrio , and Prevotella rDNA increased in read abundance during secondary colonization, whilst Blautia decreased in read abundance. MG-RAST Clusters of Orthologous Groups (COG) functional analysis also showed that the primary function of the attached microbiome was categorized broadly within "metabolism;" predominantly amino acid, carbohydrate, and lipid metabolism and transport. Most sequence read abundances (51.6, 43.8, and 50.0% of COG families pertaining to amino acid, carbohydrate and lipid metabolism, respectively) within these categories were higher in abundance during secondary colonization. Kyoto encyclopedia of genes and genomes (KEGG) pathways analysis confirmed that the PRG-attached microbiota present at 1 and 4 h of rumen incubation possess a similar functional capacity, with only a few pathways being uniquely found in only one incubation time point only. FT-IR data for the plant residues also showed that the main changes in plant chemistry between primary and secondary colonization was due to increased carbohydrate, amino acid, and lipid metabolism. This study confirmed primary and secondary colonization events and supported the hypothesis that functional changes occurred as a consequence of taxonomical changes. Sequences within the carbohydrate metabolism COG families contained only 3.2% of cellulose activities, on average across both incubation times (1 and 4 h), suggesting that degradation of the plant cell walls may be a key rate-limiting factor in ensuring the bioavailability of intra-plant nutrients in a timely manner to the microbes and ultimately the animal. This suggests that a future focus for improving ruminant nutrient use efficiency should be altering the recalcitrant plant cell wall components and/or improving the cellulolytic capacity of the rumen microbiota.

Background: Sainfoin (Onobrychis viciifolia) is a highly nutritious tannin-containing forage legume. In the diet of ruminants sainfoin can have anti-parasitic effects and reduce methane emissions under in vitro conditions. Many of these benefits have been attributed to condensed tannins or proanthocyanidins in sainfoin. A combination of increased use of industrially produced nitrogen fertilizer, issues with establishment and productivity in the first year and more reliable alternatives, such as red clover led to a decline in the use of sainfoin since the middle of the last century. In recent years there has been a resurgence of interest in sainfoin due to its potential beneficial nutraceutical and environmental attributes. However, genomic resources are scarce, thus hampering progress in genetic analysis and improvement. To address this we have used next generation RNA sequencing technology to obtain the first transcriptome of sainfoin. We used the library to identify gene-based simple sequence repeats (SSRs) and potential single nucleotide polymorphisms (SNPs)., Results: One genotype from each of five sainfoin accessions was sequenced. Paired-end (PE) sequences were generated from cDNA libraries of RNA extracted from 7 day old seedlings. A combined assembly of 92,772 transcripts was produced de novo using the Trinity programme. About 18,000 transcripts were annotated with at least one GO (gene ontology) term. A total of 63 transcripts were annotated as involved in the tannin biosynthesis pathway. We identified 3786 potential SSRs. SNPs were identified by mapping the reads of the individual assemblies against the combined assembly. After stringent filtering a total of 77,000 putative SNPs were identified. A phylogenetic analysis of single copy number genes showed that sainfoin was most closely related to red clover and Medicago truncatula, while Lotus japonicus, bean and soybean are more distant relatives., Conclusions: This work describes the first transcriptome assembly in sainfoin. The 92 K transcripts provide a rich source of SNP and SSR polymorphisms for future use in genetic studies of this crop. Annotation of genes involved in the condensed tannin biosynthesis pathway has provided the basis for further studies of the genetic control of this important trait in sainfoin.

Chronic Obstructive Pulmonary Disease (COPD) is a major source of mortality and morbidity worldwide. The microbiome associated with this disease may be an important component of the disease, though studies to date have been based on sequencing of the 16S rRNA gene, and have revealed unequivocal results. Here, we employed metagenomic sequencing of the upper bronchial tract (UBT) microbiome to allow for greater elucidation of its taxonomic composition, and revealing functional changes associated with the disease. The bacterial metagenomes within sputum samples from eight COPD patients and ten 'healthy' smokers (Controls) were sequenced, and suggested significant changes in the abundance of bacterial species, particularly within the Streptococcus genus. The functional capacity of the COPD UBT microbiome indicated an increased capacity for bacterial growth, which could be an important feature in bacterial-associated acute exacerbations. Regression analyses correlated COPD severity (FEV1% of predicted) with differences in the abundance of Streptococcus pneumoniae and functional classifications related to a reduced capacity for bacterial sialic acid metabolism. This study suggests that the COPD UBT microbiome could be used in patient risk stratification and in identifying novel monitoring and treatment methods, but study of a longitudinal cohort will be required to unequivocally relate these features of the microbiome with COPD severity.

The temporal variability of the human microbiome may be an important factor in determining its relationship with health and disease. In this study, the saliva of 40 participants was collected every 2 months over a one-year period to determine the temporal variability of the human salivary microbiome. Salivary pH and 16S rRNA gene copy number were measured for all participants, with the microbiome of 10 participants assessed through 16S rRNA amplicon sequencing. In February 2013, 16S rRNA gene copy number was significantly (P < 0.001) higher, with individual changes between time points significant (P = 0.003). Salivary pH levels were significantly (P < 0.001) higher in December 2012 than in October 2012 and February 2013, with significant (P < 0.001) individual variations seen throughout. Bacterial α-diversity showed significant differences between participants (P < 0.001), but not sampling periods (P = 0.801), and a significant positive correlation with salivary pH (R(2) = 7.8%; P = 0.019). At the phylum level, significant differences were evident between participants in the Actinobacteria (P < 0.001), Bacteroidetes (P < 0.001), Firmicutes (P = 0.008), Fusobacteria (P < 0.001), Proteobacteria (P < 0.001), Synergistetes (P < 0.001) and Spirochaetes (P = 0.003) phyla. This study charted the temporal variability of the salivary microbiome, suggesting that bacterial diversity is stable, but that 16S rRNA gene copy number may be subject to seasonal flux., (© FEMS 2015. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Background: Whole genome sequencing (WGS) of Cryptosporidium spp. has previously relied on propagation of the parasite in animals to generate enough oocysts from which to extract DNA of sufficient quantity and purity for analysis. We have developed and validated a method for preparation of genomic Cryptosporidium DNA suitable for WGS directly from human stool samples and used it to generate 10 high-quality whole Cryptosporidium genome assemblies. Our method uses a combination of salt flotation, immunomagnetic separation (IMS), and surface sterilisation of oocysts prior to DNA extraction, with subsequent use of the transposome-based Nextera XT kit to generate libraries for sequencing on Illumina platforms. IMS was found to be superior to caesium chloride density centrifugation for purification of oocysts from small volume stool samples and for reducing levels of contaminant DNA., Results: The IMS-based method was used initially to sequence whole genomes of Cryptosporidium hominis gp60 subtype IbA10G2 and Cryptosporidium parvum gp60 subtype IIaA19G1R2 from small amounts of stool left over from diagnostic testing of clinical cases of cryptosporidiosis. The C. parvum isolate was sequenced to a mean depth of 51.8X with reads covering 100 % of the bases of the C. parvum Iowa II reference genome (Bioproject PRJNA 15586), while the C. hominis isolate was sequenced to a mean depth of 34.7X with reads covering 98 % of the bases of the C. hominis TU502 v1 reference genome (Bioproject PRJNA 15585). The method was then applied to a further 17 stools, successfully generating another eight new whole genome sequences, of which two were C. hominis (gp60 subtypes IbA10G2 and IaA14R3) and six C. parvum (gp60 subtypes IIaA15G2R1 from three samples, and one each of IIaA17G1R1, IIaA18G2R1, and IIdA22G1), demonstrating the utility of this method to sequence Cryptosporidium genomes directly from clinical samples. This development is especially important as it reduces the requirement to propagate Cryptosporidium oocysts in animal models prior to genome sequencing., Conclusion: This represents the first report of high-quality whole genome sequencing of Cryptosporidium isolates prepared directly from human stool samples.

Identifying the floral composition of honey provides a method for investigating the plants that honey bees visit. We compared melissopalynology, where pollen grains retrieved from honey are identified morphologically, with a DNA metabarcoding approach using the rbcL DNA barcode marker and 454-pyrosequencing. We compared nine honeys supplied by beekeepers in the UK. DNA metabarcoding and melissopalynology were able to detect the most abundant floral components of honey. There was 92% correspondence for the plant taxa that had an abundance of over 20%. However, the level of similarity when all taxa were compared was lower, ranging from 22-45%, and there was little correspondence between the relative abundance of taxa found using the two techniques. DNA metabarcoding provided much greater repeatability, with a 64% taxa match compared to 28% with melissopalynology. DNA metabarcoding has the advantage over melissopalynology in that it does not require a high level of taxonomic expertise, a greater sample size can be screened and it provides greater resolution for some plant families. However, it does not provide a quantitative approach and pollen present in low levels are less likely to be detected. We investigated the plants that were frequently used by honey bees by examining the results obtained from both techniques. Plants with a broad taxonomic range were detected, covering 46 families and 25 orders, but a relatively small number of plants were consistently seen across multiple honey samples. Frequently found herbaceous species were Rubus fruticosus, Filipendula ulmaria, Taraxacum officinale, Trifolium spp., Brassica spp. and the non-native, invasive, Impatiens glandulifera. Tree pollen was frequently seen belonging to Castanea sativa, Crataegus monogyna and species of Malus, Salix and Quercus. We conclude that although honey bees are considered to be supergeneralists in their foraging choices, there are certain key species or plant groups that are particularly important in the honey bees environment. The reasons for this require further investigation in order to better understand honey bee nutritional requirements. DNA metabarcoding can be easily and widely used to investigate floral visitation in honey bees and can be adapted for use with other insects. It provides a starting point for investigating how we can better provide for the insects that we rely upon for pollination.

Most species exist as subdivided ex situ daughter population(s) derived from a single original group of individuals. Such subdivision occurs for many reasons both natural and manmade. Traditional British and Irish pony breeds were introduced to North America (U.S.A. and Canada) within the last 150 years, and subsequently equivalent breed societies were established. We have analyzed selected U.K. and North American equivalent pony populations as a case study for understanding the relationship between putative source and derived subpopulations. Diversity was measured using mitochondrial DNA and a panel of microsatellite markers. Genetic signatures differed between the North American subpopulations according to historical management processes. Founder effect and stochastic drift was apparent, particularly pronounced in some breeds, with evidence of admixture of imported mares of different North American breeds. This demonstrates the importance of analysis of subpopulations to facilitate understanding the genetic effects of past management practices and to lead to informed future conservation strategies.

Obligate bacterial symbionts are widespread in many invertebrates, where they are often confined to specialized host cells and are transmitted directly from mother to progeny. Increasing numbers of these bacteria are being characterized but questions remain about their population structure and evolution. Here we take a comparative genomics approach to investigate two prominent bacterial symbionts (BFo1 and BFo2) isolated from geographically separated populations of western flower thrips, Frankliniella occidentalis. Our multifaceted approach to classifying these symbionts includes concatenated multilocus sequence analysis (MLSA) phylogenies, ribosomal multilocus sequence typing (rMLST), construction of whole-genome phylogenies, and in-depth genomic comparisons. We showed that the BFo1 genome clusters more closely to species in the genus Erwinia, and is a putative close relative to Erwinia aphidicola. BFo1 is also likely to have shared a common ancestor with Erwinia pyrifoliae/Erwinia amylovora and the nonpathogenic Erwinia tasmaniensis and genetic traits similar to Erwinia billingiae. The BFo1 genome contained virulence factors found in the genus Erwinia but represented a divergent lineage. In contrast, we showed that BFo2 belongs within the Enterobacteriales but does not group closely with any currently known bacterial species. Concatenated MLSA phylogenies indicate that it may have shared a common ancestor to the Erwinia and Pantoea genera, and based on the clustering of rMLST genes, it was most closely related to Pantoea ananatis but represented a divergent lineage. We reconstructed a core genome of a putative common ancestor of Erwinia and Pantoea and compared this with the genomes of BFo bacteria. BFo2 possessed none of the virulence determinants that were omnipresent in the Erwinia and Pantoea genera. Taken together, these data are consistent with BFo2 representing a highly novel species that maybe related to known Pantoea., (© The Author(s) 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background: Understanding the evolution of the vertebrate pancreas is key to understanding its functions. The chondrichthyes (cartilaginous fish such as sharks and rays) have often been suggested to possess the most ancient example of a distinct pancreas with both hormonal (endocrine) and digestive (exocrine) roles. The lack of genetic, genomic and transcriptomic data for cartilaginous fish has hindered a more thorough understanding of the molecular-level functions of the chondrichthyan pancreas, particularly with respect to their "unusual" energy metabolism (where ketone bodies and amino acids are the main oxidative fuel source) and their paradoxical ability to both maintain stable blood glucose levels and tolerate extensive periods of hypoglycemia. In order to shed light on some of these processes, we carried out the first large-scale comparative transcriptomic survey of multiple cartilaginous fish tissues: the pancreas, brain and liver of the lesser spotted catshark, Scyliorhinus canicula., Results: We generated a mutli-tissue assembly comprising 86,006 contigs, of which 44,794 were assigned to a particular tissue or combination of tissues based on mapping of sequencing reads. We have characterised transcripts encoding genes involved in insulin regulation, glucose sensing, transcriptional regulation, signaling and digestion, as well as many peptide hormone precursors and their receptors for the first time. Comparisons to mammalian pancreas transcriptomes reveals that mechanisms of glucose sensing and insulin regulation used to establish and maintain a stable internal environment are conserved across jawed vertebrates and likely pre-date the vertebrate radiation. Conservation of pancreatic hormones and genes encoding digestive proteins support the single, early evolution of a distinct pancreatic gland with endocrine and exocrine functions in jawed vertebrates. In addition, we demonstrate that chondrichthyes lack pancreatic polypeptide (PP) and that reports of PP in the literature are likely due cross-reaction with PYY and/or NPY in the pancreas. A three hormone islet organ is therefore the ancestral jawed vertebrate condition, later elaborated upon only in the tetrapod lineage., Conclusions: The cartilaginous fish are a great untapped resource for the reconstruction of patterns and processes of vertebrate evolution and new approaches such as those described in this paper will greatly facilitate their incorporation into the rank of "model organism".

Snake venom has been hypothesized to have originated and diversified through a process that involves duplication of genes encoding body proteins with subsequent recruitment of the copy to the venom gland, where natural selection acts to develop or increase toxicity. However, gene duplication is known to be a rare event in vertebrate genomes, and the recruitment of duplicated genes to a novel expression domain (neofunctionalization) is an even rarer process that requires the evolution of novel combinations of transcription factor binding sites in upstream regulatory regions. Therefore, although this hypothesis concerning the evolution of snake venom is very unlikely and should be regarded with caution, it is nonetheless often assumed to be established fact, hindering research into the true origins of snake venom toxins. To critically evaluate this hypothesis, we have generated transcriptomic data for body tissues and salivary and venom glands from five species of venomous and nonvenomous reptiles. Our comparative transcriptomic analysis of these data reveals that snake venom does not evolve through the hypothesized process of duplication and recruitment of genes encoding body proteins. Indeed, our results show that many proposed venom toxins are in fact expressed in a wide variety of body tissues, including the salivary gland of nonvenomous reptiles and that these genes have therefore been restricted to the venom gland following duplication, not recruited. Thus, snake venom evolves through the duplication and subfunctionalization of genes encoding existing salivary proteins. These results highlight the danger of the elegant and intuitive "just-so story" in evolutionary biology., (© The Author(s) 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background: Red clover (Trifolium pratense L.) is a versatile forage crop legume, which can tolerate a variety of soils and is suitable for silage production for winter feed and for grazing. It is one of the most important forage legumes in temperate livestock agriculture. Its beneficial attributes include ability to fix nitrogen, improve soil and provide protein rich animal feed. It is however, a short-lived perennial providing good biomass yield for two or three years. Improved persistency is thus a major breeding target. Better water-stress tolerance is one of the key factors influencing persistency, but little is known about how red clover tolerates water stress., Results: Plants from a full sib mapping family were used in a drought experiment, in which the growth rate and relative water content (RWC) identified two pools of ten plants contrasting in their tolerance to drought. Key metabolites were measured and RNA-Seq analysis was carried out on four bulked samples: the two pools sampled before and after drought. Massively parallel sequencing was used to analyse the bulked RNA samples. A de novo transcriptome reconstruction based on the RNA-Seq data was made, resulting in 45181 contigs, representing 'transcript tags'. These transcript tags were annotated with gene ontology (GO) terms. One of the most striking results from the expression analysis was that the drought sensitive plants were characterised by having approximately twice the number of differentially expressed transcript tags than the tolerant plants after drought. This difference was evident in most of the major GO terms. Before onset of drought the sensitive plants overexpressed a number of genes annotated as senescence-related. Furthermore, the concentration of three metabolites, particularly pinitol, but also proline and malate increased in leaves after drought stress., Conclusions: This de novo assembly of a red clover transcriptome from leaf material of droughted and non-droughted plants provides a rich source for gene identification, single nucleotide polymorphisms (SNP) and short sequence repeats (SSR). Comparison of gene expression levels between pools and treatments identified candidate genes for further analysis of the genetic basis of drought tolerance in red clover.

The conservation of unique populations of animals is critical in order to preserve valuable genetic diversity and, where populations are free-living, maintain their irreplaceable influence upon habitat ecology. An accurate assessment of genetic diversity and structure within and between populations is crucial in order to design and implement conservation strategies in natural and domesticated species. Moreover, where it is possible to identify relic populations that are related to a structured breed an ideal opportunity presents itself to model processes that reveal historical factors that have shaped genetic diversity. The origins of native UK mountain and moorland ponies are uncertain, but they may have directly descended from prehistoric populations and potentially harbour specific adaptations to the uplands of Britain and Ireland. To date, there have been no studies of population structure and genetic diversity present within a free-living group of ponies in the Carneddau mountain range of North Wales. Herein, we describe the use of microsatellites and SNPs together with analysis of the mitochondrial control region to quantify the extent and magnitude of genetic diversity present in the feral Carneddau pony and relate this to several recognised British and Irish pony breeds. Our results establish that the feral Carneddau ponies represent a unique and distinctive population that merits recognition as a defined population and conservation priority. We discuss the implications for conservation of this population as a unique pool of genetic diversity adapted to the British uplands and potentially of particular value in maintaining the biodiversity of these habitats.

Human activity and climate change are increasingly driving species, which were once separate together, leading to the potential for gene flow. Hybridization between diverged species brings together two genomes which have evolved to meet different adaptive requirements. The unique combination of these traits in a hybrid may be beneficial or maladaptive, but either way it results in increased phenotypic variation. A percentage of hybrid individuals may, therefore, find themselves able to exploit environmental niches which their progenitors cannot, leading to invasive hybrid swarms becoming established in new habitats. Previous research into hybrids, most famously that of Loren Rieseberg and co-workers (Rieseberg et al. 1999, 2003) in sunflowers, demonstrated that hybridization can give rise to transgressive segregation of adaptive traits, wherein the combination of favourable alleles from both parents in hybrids can enable them to outperform either. However, the question still remains as to how much of the competitive ability of hybrids is a direct result of admixture and how much is the result of selection after the fact. In this issue of Molecular Ecology, (Czypionka et al. 2012) describe their study of transcriptional changes resulting from hybridization in a fish hybrid termed invasive sculpins (Cottus). Using gene expression microarray assays, they compare gene expression in both wild and lab-reared invasive hybrids to the progenitor species and experimentally produced F(2) hybrids. They demonstrate that whilst hybridization alone does result in higher variance in gene expression (some of which is transgressive), many of the transgressive changes distinguishing the invasives appear to have come about subsequent to the initial natural hybridization event. They speculate that initial success of the hybrids in their new habitat is facilitated by hybridization, but that optimization of the invasive phenotype and removal of maladaptive traits rapidly reduces the variation in gene expression seen in early hybrids., (© 2012 Blackwell Publishing Ltd.)

Background: Pollen-pistil interactions are an essential prelude to fertilization in angiosperms and determine compatibility/incompatibility. Pollen-pistil interactions have been studied at a molecular and cellular level in relatively few families. Self-incompatibility (SI) is the best understood pollen-pistil interaction at a molecular level where three different molecular mechanisms have been identified in just five families. Here we review studies of pollen-pistil interactions and SI in the Asteraceae, an important family that has been relatively understudied in these areas of reproductive biology., Scope: We begin by describing the historical literature which first identified sporophytic SI (SSI) in species of Asteraceae, the SI system later identified and characterized at a molecular level in the Brassicaceae. Early structural and cytological studies in these two families suggested that pollen-pistil interactions and SSI were similar, if not the same. Recent cellular and molecular studies in Senecio squalidus (Oxford ragwort) have challenged this belief by revealing that despite sharing the same genetic system of SSI, the Brassicaceae and Asteraceae molecular mechanisms are different. Key cellular differences have also been highlighted in pollen-stigma interactions, which may arise as a consequence of the Asteraceae possessing a 'semi-dry' stigma, rather than the 'dry' stigma typical of the Brassicaceae. The review concludes with a summary of recent transcriptomic analyses aimed at identifying proteins regulating pollen-pistil interactions and SI in S. squalidus, and by implication the Asteraceae. The Senecio pistil transcriptome contains many novel pistil-specific genes, but also pistil-specific genes previously shown to play a role in pollen-pistil interactions in other species., Conclusions: Studies in S. squalidus have shown that stigma structure and the molecular mechanism of SSI in the Asteraceae and Brassicaceae are different. The availability of a pool of pistil-specific genes for S. squalidus offers an opportunity to elucidate the molecular mechanisms of pollen-pistil interactions and SI in the Asteraceae.

The merger of two or more divergent genomes within an allopolyploid nucleus can facilitate speciation and adaptive evolution in flowering plants. Widespread changes to gene expression have been shown to result from interspecific hybridisation and polyploidy in a number of plant species, and attention has now shifted to determining the epigenetic processes that drive these changes. We present here an analysis of cytosine methylation patterns in triploid F(1) Senecio (ragwort) hybrids and their allohexaploid derivatives. We observe that, in common with similar studies in Arabidopsis, Spartina and Triticum, a small but significant proportion of loci display nonadditive methylation in the hybrids, largely resulting from interspecific hybridisation. Despite this, genome duplication results in a secondary effect on methylation, with reversion to additivity at some loci and novel methylation status at others. We also observe differences in methylation state between different allopolyploid generations, predominantly in cases of additive methylation with regard to which parental methylation state is dominant. These changes to methylation state in both F(1) triploids and their allohexaploid derivatives largely mirror the overall patterns of nonadditive gene expression observed in our previous microarray analyses and may play a causative role in generating those expression changes. These similar global changes to DNA methylation resulting from hybridisation and genome duplication may serve as a source of epigenetic variation in natural populations, facilitating adaptive evolution. Our observations that methylation state can also vary between different generations of polyploid hybrids suggests that newly formed allopolyploid species may display a high degree of epigenetic diversity upon which natural selection can act., (© 2010 Blackwell Publishing Ltd.)

Hybridization is an important cause of abrupt speciation. Hybrid speciation without a change in ploidy (homoploid hybrid speciation) is well-established in plants but has also been reported in animals and fungi. A notable example of recent homoploid hybrid speciation is Senecio squalidus (Oxford ragwort), which originated in the UK in the 18th Century following introduction of hybrid material from a hybrid zone between S. chrysanthemifolius and S. aethnensis on Mount Etna, Sicily. To investigate genetic divergence between these taxa, we used complementary DNA microarrays to compare patterns of floral gene expression. These analyses revealed major differences in gene expression between the parent species and wild and resynthesized S. squalidus. Comparisons of gene expression between S. aethnensis, S. chrysanthemifolius and natural S. squalidus identified genes potentially involved in local environmental adaptation. The analysis also revealed non-additive patterns of gene expression in the hybrid relative to its progenitors. These expression changes were more dramatic and widespread in resynthesized hybrids than in natural S. squalidus, suggesting that a unique expression pattern may have been fixed during the allopatric divergence of British S. squalidus. We speculate that hybridization-induced gene-expression change may provide an immediate source of novel phenotypic variation upon which selection can act to facilitate homoploid hybrid speciation in plants.

Interspecific hybridization is an important mechanism of speciation in higher plants. In flowering plants, hybrid speciation is usually associated with polyploidy (allopolyploidy), but hybrid speciation without genome duplication (homoploid hybrid speciation) is also possible, although it is more difficult to detect. The combination of divergent genomes within a hybrid can result in profound changes to both genome and transcriptome. Recent transcriptomic studies of wild and resynthesized homoploid and allopolyploid hybrids have revealed widespread changes to gene expression in hybrids relative to expression levels in their parents. Many of these changes to gene expression are 'non-additive', i.e. not simply the sum of the combined expression levels of parental genes. Some gene expression changes are far outside the range of gene expression in either parent, and can therefore be viewed as 'transgressive'. Such profound changes to gene expression may enable new hybrids to survive in novel habitats not accessible to their parent species. Here, we give a brief overview of hybrid speciation in plants, with an emphasis on genomic change, before focusing discussion on findings from recent transcriptomic studies. We then discuss our current work on gene expression change associated with hybrid speciation in the genus Senecio (ragworts and groundsels) focusing on the findings from a reanalysis of gene expression data obtained from recent microarray studies of wild and resynthesized allopolyploid Senecio cambrensis. These data, showing extensive non-additive and transgressive gene expression changes in Senecio hybrids, are discussed in the light of findings from other model systems, and in the context of the potential importance of gene expression change to hybrid speciation in plants.

Polyploidy, or the presence of two or more diploid parental genome sets within an organism, is found to an amazing degree in higher plants. In addition, many plant species traditionally considered to be diploid have recently been demonstrated to have undergone rounds of genome duplication in the past and are now referred to as paleopolyploids. Polyploidy and interspecific hybridisation (with which it is often associated) have long been thought to be important mechanisms of rapid species formation. The widespread occurrence of polyploids, which are frequently found in habitats different from that of their diploid progenitors, would seem to indicate that polyploidy is associated with evolutionary success in terms of the ability to colonise new environmental niches. A flurry of recent genomic studies has provided fresh insights into the potential basis of the phenotypic novelty of polyploid species. Here we review current knowledge of genetic, epigenetic, and transcriptional changes associated with polyploidy in plants and assess how these changes might contribute to the evolutionary success of polyploid plants. We conclude by stressing the need for field-based experiments to determine whether genetic changes associated with polyploidy are indeed adaptive.

Allopolyploidy, which involves genome doubling of an interspecific hybrid is an important mechanism of abrupt speciation in flowering plants [1-6]. Recent studies show that allopolyploid formation is accompanied by extensive changes to patterns of parental gene expression ("transcriptome shock") [7-15] and that this is likely the consequence of interspecific hybridization rather than polyploidization [16]. To investigate the relative impacts of hybridization and polyploidization on transcription, we compared floral gene expression in allohexaploid Senecio cambrensis with that in its parent species, S. vulgaris (tetraploid) and S. squalidus (diploid), and their triploid F1 hybrid, S. x baxteri [17]. Major changes to parental gene expression were associated principally with S. x baxteri, suggesting that the polyploidization event responsible for the formation of S. cambrensis had a widespread calming effect on altered gene expression arising from hybridization [17]. To test this hypothesis, we analyzed floral gene expression in resynthesized lines of S. cambrensis and show that, for many genes, the "transcriptome shock" observed in S. x baxteri is calmed ("ameliorated") after genome doubling in the first generation of synthetic S. cambrensis and this altered expression pattern is maintained in subsequent generations. These findings indicate that hybridization and polyploidization have immediate yet distinct effects on large-scale patterns of gene expression.

Interspecific hybridization is an important process through which abrupt speciation can occur. In recent years, genetic changes associated with hybrid speciation have been identified through a variety of techniques, including AFLP/SSR mapping, GISH/FISH and cDNA-AFLP differential display. However, progress in using microarray technology to analyse whole genome/transcriptome changes associated with hybrid speciation has been limited due to the lack of extensive sequence data for many hybrid species and the difficulties in extrapolating results from commercially available microarrays for model species onto nonmodel hybrid taxa. Increasingly therefore researchers studying nonmodel systems are turning to the development of 'anonymous' cDNA microarrays, where the time and cost of producing microarrays is reduced by printing unsequenced cDNA clones, and sequencing only those clones that display interesting expression patterns. Here we describe the creation, testing and preliminary use of anonymous cDNA microarrays to study changes in floral transcriptome associated with allopolyploid speciation in the genus Senecio. We report a comparison of gene expression between the allohexaploid hybrid, Senecio cambrensis, its parental taxa Senecio squalidus (diploid) and Senecio vulgaris (tetraploid), and the intermediate triploid (sterile) hybrid Senecioxbaxteri. Anonymous microarray analysis revealed dramatic differences in floral gene expression between these four taxa and demonstrates the power of this technique for studies of the genetic impact of hybridization in nonmodel flowering plants.

Abrupt speciation through interspecific hybridisation is an important mechanism in angiosperm evolution. Flowering plants therefore offer excellent opportunities for studying genetic processes associated with hybrid speciation. Novel molecular approaches are now available to examine these processes at the level of both genome organization and gene expression - transcriptomics. Here, we present an overview of the molecular technologies currently used to study hybrid speciation and how they are providing new insights into this mode of speciation in flowering plants. We begin with an introduction to hybrid speciation in plants, followed by a review of techniques, such as isozymes and other markers, which have been used to study hybrid species in the past. We then review advances in molecular techniques that have the potential to be applied to studies of hybrid species, followed by an overview of the main genomic and transcriptomic changes suspected, or known, to occur in newly formed hybrids, together with commentary on the application of advanced molecular tools to studying these changes.