Matthew Hegarty

De novo assembly of red clover transcriptome based on RNA-Seq data provides insight into drought response, gene discovery and marker identification

BackgroundRed 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.ResultsPlants 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.ConclusionsThis 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.

Nonadditive changes to cytosine methylation as a consequence of hybridization and genome duplication in Senecio (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 F1Senecio (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 F1 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.

Restriction and Recruitment—Gene Duplication and the Origin and Evolution of Snake Venom Toxins

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.

Using DNA Metabarcoding to Identify the Floral Composition of Honey: A New Tool for Investigating Honey Bee Foraging Preferences

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.

A metagenomic snapshot of taxonomic and functional diversity in an alpine glacier cryoconite ecosystem

Cryoconite is a microbe‐mineral aggregate which darkens the ice surface of glaciers. Microbial process and marker gene PCR-dependent measurements reveal active and diverse cryoconite microbial communities on polar glaciers. Here, we provide the first report of a cryoconite metagenome and culture-independent study of alpine cryoconite microbial diversity. We assembled 1.2 Gbp of metagenomic DNA sequenced using an Illumina HiScanSQ from cryoconite holes across the ablation zone of Rotmoosferner in the Austrian Alps. The metagenome revealed a bacterially-dominated community, with Proteobacteria (62% of bacterialassigned contigs) and Bacteroidetes (14%) considerably more abundant than Cyanobacteria (2.5%). Streptophyte DNA dominated the eukaryotic metagenome. Functional genes linked to N, Fe, S and P cycling illustrated an acquisitive trend and a nitrogen cycle based upon efficient ammonia recycling. A comparison of 32 metagenome datasets revealed a similarity in functional profiles between the cryoconite and metagenomes characterized from other cold microbe‐mineral aggregates. Overall, the metagenomic snapshot reveals the cryoconite ecosystem of this alpine glacier as dependent on scavenging carbon and nutrients from allochthonous sources, in particular mosses transported by wind from ice-marginal habitats, consistent with net heterotrophy indicated by productivity measurements. A transition from singular snapshots of cryoconite metagenomes to comparative analyses is advocated.

Genotyping by RAD sequencing enables mapping of fatty acid composition traits in perennial ryegrass (Lolium perenne (L.))

Perennial ryegrass (Lolium perenne L.) is the most important forage crop in temperate livestock agriculture. Its nutritional quality has significant impact on the quality of meat and milk for human consumption. Evidence suggests that higher energy content in forage can assist in reducing greenhouse gas emissions from ruminants. Increasing the fatty acid content (especially α-linolenic acid, an omega-3 fatty acid) may thus contribute to better forage, but little is known about the genetic basis of variation for this trait. To this end, quantitative trait loci (QTLs) were identified associated with major fatty acid content in perennial ryegrass using a population derived from a cross between the heterozygous and outbreeding high-sugar grass variety AberMagic and an older variety, Aurora. A genetic map with 434 restriction-associated DNA (RAD) and SSR markers was generated. Significant QTLs for the content of palmitic (C16:0) on linkage groups (LGs) 2 and 7; stearic (C18:0) on LGs 3, 4 and 7; linoleic (C18:2n-6) on LGs 2 and 5; and α-linolenic acids (C18:3n-3) on LG 1 were identified. Two candidate genes (a lipase and a beta-ketoacyl CoA synthase), both associated with C16:0, and separately with C18:2n-6 and C18:0 contents, were identified. The physical positions of these genes in rice and their genetic positions in perennial ryegrass were consistent with established syntenic relationships between these two species. Validation of these associations is required, but the utility of RAD markers for rapid generation of genetic maps and QTL analysis has been demonstrated for fatty acid composition in a global forage crop.

Red clover (Trifolium pratense L.) draft genome provides a platform for trait improvement.

Red clover (Trifolium pratense L.) is a globally significant forage legume in pastoral livestock farming systems. It is an attractive component of grassland farming, because of its high yield and protein content, nutritional value and ability to fix atmospheric nitrogen. Enhancing its role further in sustainable agriculture requires genetic improvement of persistency, disease resistance, and tolerance to grazing. To help address these challenges, we have assembled a chromosome-scale reference genome for red clover. We observed large blocks of conserved synteny with Medicago truncatula and estimated that the two species diverged ~23 million years ago. Among the 40,868 annotated genes, we identified gene clusters involved in biochemical pathways of importance for forage quality and livestock nutrition. Genotyping by sequencing of a synthetic population of 86 genotypes show that the number of markers required for genomics-based breeding approaches is tractable, making red clover a suitable candidate for association studies and genomic selection.

The human salivary microbiome exhibits temporal stability in bacterial diversity

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.

Using DNA metabarcoding to investigate honey bee foraging reveals limited flower use despite high floral availability.

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.

Metagenomic Sequencing of the Chronic Obstructive Pulmonary Disease Upper Bronchial Tract Microbiome Reveals Functional Changes Associated with Disease Severity

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.