Faecal metabarcoding reveals pervasive long-distance impacts of garden bird feeding.

Supplementary feeding of wildlife is widespread, being undertaken by more than half of households in many countries. However, the impact that these supplemental resources have is unclear, with impacts largely considered to be restricted to urban ecosystems. We reveal the pervasiveness of supplementary foodstuffs in the diet of a wild bird using metabarcoding of blue tit (Cyanistes caeruleus) faeces collected in early spring from a 220 km transect in Scotland with a large urbanization gradient. Supplementary foodstuffs were present in the majority of samples, with peanut (Arachis hypogaea) the single commonest (either natural or supplementary) dietary item. Consumption rates exhibited a distance decay from human habitation but remained high at several hundred metres from the nearest household and continued to our study limit of 1.4 km distant. Supplementary food consumption was associated with a near quadrupling of blue tit breeding density and a 5-day advancement of breeding phenology. We show that woodland bird species using supplementary food have increasing UK population trends, while species that do not, and/or are outcompeted by blue tits, are likely to be declining. We suggest that the impacts of supplementary feeding are larger and more spatially extensive than currently appreciated and could be disrupting population and ecosystem dynamics.

Gradients in richness and turnover of a forest passerine's diet prior to breeding: A mixed model approach applied to faecal metabarcoding data.

Rather little is known about the dietary richness and variation of generalist insectivorous species, including birds, due primarily to difficulties in prey identification. Using faecal metabarcoding, we provide the most comprehensive analysis of a passerine's diet to date, identifying the relative magnitudes of biogeographic, habitat and temporal trends in the richness and turnover in diet of Cyanistes caeruleus (blue tit) along a 39 site and 2° latitudinal transect in Scotland. Faecal samples were collected in 2014-2015 from adult birds roosting in nestboxes prior to nest building. DNA was extracted from 793 samples and we amplified COI and 16S minibarcodes. We identified 432 molecular operational taxonomic units that correspond to putative dietary items. Most dietary items were rare, with Lepidoptera being the most abundant and taxon-rich prey order. Here, we present a statistical approach for estimation of gradients and intersample variation in taxonomic richness and turnover using a generalised linear mixed model. We discuss the merits of this approach over existing tools and present methods for model-based estimation of repeatability, taxon richness and Jaccard indices. We found that dietary richness increases significantly as spring advances, but changes little with elevation, latitude or local tree composition. In comparison, dietary composition exhibits significant turnover along temporal and spatial gradients and among sites. Our study shows the promise of faecal metabarcoding for inferring the macroecology of food webs, but we also highlight the challenge posed by contamination and make recommendations of laboratory and statistical practices to minimise its impact on inference.

The genome of the emerging barley pathogen Ramularia collo-cygni.

BACKGROUND: Ramularia collo-cygni is a newly important, foliar fungal pathogen of barley that causes the disease Ramularia leaf spot. The fungus exhibits a prolonged endophytic growth stage before switching life habit to become an aggressive, necrotrophic pathogen that causes significant losses to green leaf area and hence grain yield and quality. RESULTS: The R. collo-cygni genome was sequenced using a combination of Illumina and Roche 454 technologies. The draft assembly of 30.3 Mb contained 11,617 predicted gene models. Our phylogenomic analysis confirmed the classification of this ascomycete fungus within the family Mycosphaerellaceae, order Capnodiales of the class Dothideomycetes. A predicted secretome comprising 1053 proteins included redox-related enzymes and carbohydrate-modifying enzymes and proteases. The relative paucity of plant cell wall degrading enzyme genes may be associated with the stealth pathogenesis characteristic of plant pathogens from the Mycosphaerellaceae. A large number of genes associated with secondary metabolite production, including homologs of toxin biosynthesis genes found in other Dothideomycete plant pathogens, were identified. CONCLUSIONS: The genome sequence of R. collo-cygni provides a framework for understanding the genetic basis of pathogenesis in this important emerging pathogen. The reduced complement of carbohydrate-degrading enzyme genes is likely to reflect a strategy to avoid detection by host defences during its prolonged asymptomatic growth. Of particular interest will be the analysis of R. collo-cygni gene expression during interactions with the host barley, to understand what triggers this fungus to switch from being a benign endophyte to an aggressive necrotroph.

Comparative transcriptomics of a complex of four European pine species.

BACKGROUND: Pinus sylvestris, P. mugo, P. uliginosa and P. uncinata are closely related but phenotypically and ecologically very distinct European pine species providing an excellent study system for analysis of the genetic basis of adaptive variation and speciation. For comparative genomic analysis of the species, transcriptome sequence was generated for 17 samples collected across the European distribution range using Illumina paired-end sequencing technology. RESULTS: De novo transcriptome assembly of a reference sample of P. sylvestris contained 40968 unigenes, of which fewer than 0.5% were identified as putative retrotransposon sequences. Based on gene annotation approaches, 19659 contigs were identified and assigned to unique genes covering a broad range of gene ontology categories. About 80% of the reads from each sample were successfully mapped to the reference transcriptome of P. sylvestris. Single nucleotide polymorphisms were identified in 22041-24096 of the unigenes providing a set of ~220-262 k SNPs identified for each species. Very similar levels of nucleotide polymorphism were observed across species (π=0.0044-0.0053) and highest pairwise nucleotide divergence (0.006) was found between P. mugo and P. sylvestris at a common set of unigenes. CONCLUSIONS: The study provides whole transcriptome sequence and a large set of SNPs to advance population and association genetic studies in pines. Our study demonstrates that transcriptome sequencing can be a very useful approach for development of novel genomic resources in species with large and complex genomes.

Genome sequence of Stenotrophomonas maltophilia PML168, which displays Baeyer-Villiger monooxygenase activity.

Stenotrophomonas maltophilia PML168 was isolated from Wembury Beach on the English Coast from a rock pool following growth and selection on agar plates. Here we present the permanent draft genome sequence, which has allowed prediction of function for several genes encoding enzymes relevant to industrial biotechnology, including a novel flavoprotein monooxygenase.

The first genome for the Cape Primrose Streptocarpus rexii (Gesneriaceae), a model plant for studying meristem-driven shoot diversity.

Cape Primroses (Streptocarpus, Gesneriaceae) are an ideal study system for investigating the genetics underlying species diversity in angiosperms. Streptocarpus rexii has served as a model species for plant developmental research for over five decades due to its unusual extended meristem activity present in the leaves. In this study, we sequenced and assembled the complete nuclear, chloroplast, and mitochondrial genomes of S. rexii using Oxford Nanopore Technologies long read sequencing. Two flow cells of PromethION sequencing resulted in 32 billion reads and were sufficient to generate a draft assembly including the chloroplast, mitochondrial and nuclear genomes, spanning 776 Mbp. The final nuclear genome assembly contained 5,855 contigs, spanning 766 Mbp of the 929-Mbp haploid genome with an N50 of 3.7 Mbp and an L50 of 57 contigs. Over 70% of the draft genome was identified as repeats. A genome repeat library of Gesneriaceae was generated and used for genome annotation, with a total of 45,045 genes annotated in the S. rexii genome. Ks plots of the paranomes suggested a recent whole genome duplication event, shared between S. rexii and Primulina huaijiensis. A new chloroplast and mitochondrial genome assembly method, based on contig coverage and identification, was developed, and successfully used to assemble both organellar genomes of S. rexii. This method was developed into a pipeline and proved widely applicable. The nuclear genome of S. rexii and other datasets generated and reported here will be invaluable resources for further research to aid in the identification of genes involved in morphological variation underpinning plant diversification.

PromethION Sequencing and Assembly of the Genome of Micropoecilia picta, a Fish with a Highly Degenerated Y Chromosome.

We here describe sequencing and assembly of both the autosomes and the sex chromosome in Micropoecilia picta, the closest related species to the guppy, Poecilia reticulata. Poecilia (Micropoecilia) picta is a close outgroup for studying the guppy, an important organism for studies in evolutionary ecology and in sex chromosome evolution. The guppy XY pair (LG12) has long been studied as a test case for the importance of sexually antagonistic variants in selection for suppressed recombination between Y and X chromosomes. The guppy Y chromosome is not degenerated, but appears to carry functional copies of all genes that are present on its X counterpart. The X chromosomes of M. picta (and its relative Micropoecilia parae) are homologous to the guppy XY pair, but their Y chromosomes are highly degenerated, and no genes can be identified in the fully Y-linked region. A complete genome sequence of a M. picta male may therefore contribute to understanding how the guppy Y evolved. These fish species' genomes are estimated to be about 750 Mb, with high densities of repetitive sequences, suggesting that long-read sequencing is needed. We evaluated several assembly approaches, and used our results to investigate the extent of Y chromosome degeneration in this species.

A telomere-to-telomere assembly of Oscheius tipulae and the evolution of rhabditid nematode chromosomes.

Eukaryotic chromosomes have phylogenetic persistence. In many taxa, each chromosome has a single functional centromere with essential roles in spindle attachment and segregation. Fusion and fission can generate chromosomes with no or multiple centromeres, leading to genome instability. Groups with holocentric chromosomes (where centromeric function is distributed along each chromosome) might be expected to show karyotypic instability. This is generally not the case, and in Caenorhabditis elegans, it has been proposed that the role of maintenance of a stable karyotype has been transferred to the meiotic pairing centers, which are found at one end of each chromosome. Here, we explore the phylogenetic stability of nematode chromosomes using a new telomere-to-telomere assembly of the rhabditine nematode Oscheius tipulae generated from nanopore long reads. The 60-Mb O. tipulae genome is resolved into six chromosomal molecules. We find the evidence of specific chromatin diminution at all telomeres. Comparing this chromosomal O. tipulae assembly with chromosomal assemblies of diverse rhabditid nematodes, we identify seven ancestral chromosomal elements (Nigon elements) and present a model for the evolution of nematode chromosomes through rearrangement and fusion of these elements. We identify frequent fusion events involving NigonX, the element associated with the rhabditid X chromosome, and thus sex chromosome-associated gene sets differ markedly between species. Despite the karyotypic stability, gene order within chromosomes defined by Nigon elements is not conserved. Our model for nematode chromosome evolution provides a platform for investigation of the tensions between local genome rearrangement and karyotypic evolution in generating extant genome architectures.

Hybrid assembly of an agricultural slurry virome reveals a diverse and stable community with the potential to alter the metabolism and virulence of veterinary pathogens.

BACKGROUND: Viruses are the most abundant biological entities on Earth, known to be crucial components of microbial ecosystems. However, there is little information on the viral community within agricultural waste. There are currently ~ 2.7 million dairy cattle in the UK producing 7-8% of their own bodyweight in manure daily, and 28 million tonnes annually. To avoid pollution of UK freshwaters, manure must be stored and spread in accordance with guidelines set by DEFRA. Manures are used as fertiliser, and widely spread over crop fields, yet little is known about their microbial composition. We analysed the virome of agricultural slurry over a 5-month period using short and long-read sequencing. RESULTS: Hybrid sequencing uncovered more high-quality viral genomes than long or short-reads alone; yielding 7682 vOTUs, 174 of which were complete viral genomes. The slurry virome was highly diverse and dominated by lytic bacteriophage, the majority of which represent novel genera (~ 98%). Despite constant influx and efflux of slurry, the composition and diversity of the slurry virome was extremely stable over time, with 55% of vOTUs detected in all samples over a 5-month period. Functional annotation revealed a diverse and abundant range of auxiliary metabolic genes and novel features present in the community, including the agriculturally relevant virulence factor VapE, which was widely distributed across different phage genera that were predicted to infect several hosts. Furthermore, we identified an abundance of phage-encoded diversity-generating retroelements, which were previously thought to be rare on lytic viral genomes. Additionally, we identified a group of crAssphages, including lineages that were previously thought only to be found in the human gut. CONCLUSIONS: The cattle slurry virome is complex, diverse and dominated by novel genera, many of which are not recovered using long or short-reads alone. Phages were found to encode a wide range of AMGs that are not constrained to particular groups or predicted hosts, including virulence determinants and putative ARGs. The application of agricultural slurry to land may therefore be a driver of bacterial virulence and antimicrobial resistance in the environment. Video abstract.

Digital gene expression analysis of two life cycle stages of the human-infective parasite, Trypanosoma brucei gambiense reveals differentially expressed clusters of co-regulated genes.

BACKGROUND: The evolutionarily ancient parasite, Trypanosoma brucei, is unusual in that the majority of its genes are regulated post-transcriptionally, leading to the suggestion that transcript abundance of most genes does not vary significantly between different life cycle stages despite the fact that the parasite undergoes substantial cellular remodelling and metabolic changes throughout its complex life cycle. To investigate this in the clinically relevant sub-species, Trypanosoma brucei gambiense, which is the causative agent of the fatal human disease African sleeping sickness, we have compared the transcriptome of two different life cycle stages, the potentially human-infective bloodstream forms with the non-human-infective procyclic stage using digital gene expression (DGE) analysis. RESULTS: Over eleven million unique tags were generated, producing expression data for 7360 genes, covering 81% of the genes in the genome. Compared to microarray analysis of the related T. b. brucei parasite, approximately 10 times more genes with a 2.5-fold change in expression levels were detected. The transcriptome analysis revealed the existence of several differentially expressed gene clusters within the genome, indicating that contiguous genes, presumably from the same polycistronic unit, are co-regulated either at the level of transcription or transcript stability. CONCLUSIONS: DGE analysis is extremely sensitive for detecting gene expression differences, revealing firstly that a far greater number of genes are stage-regulated than had previously been identified and secondly and more importantly, this analysis has revealed the existence of several differentially expressed clusters of genes present on what appears to be the same polycistronic units, a phenomenon which had not previously been observed in microarray studies. These differentially regulated clusters of genes are in addition to the previously identified RNA polymerase I polycistronic units of variant surface glycoproteins and procyclin expression sites, which encode the major surface proteins of the parasite. This raises a number of questions regarding the function and regulation of the gene clusters that clearly warrant further study.

Experimental evolution, genetic analysis and genome re-sequencing reveal the mutation conferring artemisinin resistance in an isogenic lineage of malaria parasites.

BACKGROUND: Classical and quantitative linkage analyses of genetic crosses have traditionally been used to map genes of interest, such as those conferring chloroquine or quinine resistance in malaria parasites. Next-generation sequencing technologies now present the possibility of determining genome-wide genetic variation at single base-pair resolution. Here, we combine in vivo experimental evolution, a rapid genetic strategy and whole genome re-sequencing to identify the precise genetic basis of artemisinin resistance in a lineage of the rodent malaria parasite, Plasmodium chabaudi. Such genetic markers will further the investigation of resistance and its control in natural infections of the human malaria, P. falciparum. RESULTS: A lineage of isogenic in vivo drug-selected mutant P. chabaudi parasites was investigated. By measuring the artemisinin responses of these clones, the appearance of an in vivo artemisinin resistance phenotype within the lineage was defined. The underlying genetic locus was mapped to a region of chromosome 2 by Linkage Group Selection in two different genetic crosses. Whole-genome deep coverage short-read re-sequencing (Illumina Solexa) defined the point mutations, insertions, deletions and copy-number variations arising in the lineage. Eight point mutations arise within the mutant lineage, only one of which appears on chromosome 2. This missense mutation arises contemporaneously with artemisinin resistance and maps to a gene encoding a de-ubiquitinating enzyme. CONCLUSIONS: This integrated approach facilitates the rapid identification of mutations conferring selectable phenotypes, without prior knowledge of biological and molecular mechanisms. For malaria, this model can identify candidate genes before resistant parasites are commonly observed in natural human malaria populations.

Field cricket genome reveals the footprint of recent, abrupt adaptation in the wild.

Evolutionary adaptation is generally thought to occur through incremental mutational steps, but large mutational leaps can occur during its early stages. These are challenging to study in nature due to the difficulty of observing new genetic variants as they arise and spread, but characterizing their genomic dynamics is important for understanding factors favoring rapid adaptation. Here, we report genomic consequences of recent, adaptive song loss in a Hawaiian population of field crickets (Teleogryllus oceanicus). A discrete genetic variant, flatwing, appeared and spread approximately 15 years ago. Flatwing erases sound-producing veins on male wings. These silent flatwing males are protected from a lethal, eavesdropping parasitoid fly. We sequenced, assembled and annotated the cricket genome, produced a linkage map, and identified a flatwing quantitative trait locus covering a large region of the X chromosome. Gene expression profiling showed that flatwing is associated with extensive genome-wide effects on embryonic gene expression. We found that flatwing male crickets express feminized chemical pheromones. This male feminizing effect, on a different sexual signaling modality, is genetically associated with the flatwing genotype. Our findings suggest that the early stages of evolutionary adaptation to extreme pressures can be accompanied by greater genomic and phenotypic disruption than previously appreciated, and highlight how abrupt adaptation might involve suites of traits that arise through pleiotropy or genomic hitchhiking.

High-throughput sequence analysis of variants of human cytomegalovirus strains Towne and AD169.

The genomes of commonly used variants of human cytomegalovirus (HCMV) strains Towne and AD169 each contain a substantial mutation in which a region (U(L)/b') at the right end of the long unique region has been replaced by an inverted duplication of a region from the left end of the genome. Using high-throughput technology, we have sequenced HCMV strain Towne (ATCC VR-977) and confirmed the presence of two variants, one exhibiting the replacement in U(L)/b' and the other intact in this region. Both variants are mutated in genes RL13, UL1, UL40, UL130, US1 and US9. We have also sequenced a novel AD169 variant (varUC) that is intact in U(L)/b' except for a small deletion that affects genes UL144, UL142, UL141 and UL140. Like other AD169 variants, varUC is mutated in genes RL5A, RL13, UL36 and UL131A. A subpopulation of varUC contains an additional deletion affecting genes IRS1, US1 and US2.

Lack of long-term acclimation in Antarctic encrusting species suggests vulnerability to warming.

Marine encrusting communities play vital roles in benthic ecosystems and have major economic implications with regards to biofouling. However, their ability to persist under projected warming scenarios remains poorly understood and is difficult to study under realistic conditions. Here, using heated settlement panel technologies, we show that after 18 months Antarctic encrusting communities do not acclimate to either +1 °C or +2 °C above ambient temperatures. There is significant up-regulation of the cellular stress response in warmed animals, their upper lethal temperatures decline with increasing ambient temperature and population genetic analyses show little evidence of differential survival of genotypes with treatment. By contrast, biofilm bacterial communities show no significant differences in community structure with temperature. Thus, metazoan and bacterial responses differ dramatically, suggesting that ecosystem responses to future climate change are likely to be far more complex than previously anticipated.

Signatures of the Evolution of Parthenogenesis and Cryptobiosis in the Genomes of Panagrolaimid Nematodes.

Most animal species reproduce sexually and fully parthenogenetic lineages are usually short lived in evolution. Still, parthenogenesis may be advantageous as it avoids the cost of sex and permits colonization by single individuals. Panagrolaimid nematodes have colonized environments ranging from arid deserts to Arctic and Antarctic biomes. Many are obligatory meiotic parthenogens, and most have cryptobiotic abilities, being able to survive repeated cycles of complete desiccation and freezing. To identify systems that may contribute to these striking abilities, we sequenced and compared the genomes and transcriptomes of parthenogenetic and outcrossing panagrolaimid species, including cryptobionts and non-cryptobionts. The parthenogens are triploids, most likely originating through hybridization. Adaptation to cryptobiosis shaped the genomes of panagrolaimid nematodes and is associated with the expansion of gene families and signatures of selection on genes involved in cryptobiosis. All panagrolaimids have acquired genes through horizontal gene transfer, some of which are likely to contribute to cryptobiosis.

Inter and Intraspecific Genomic Divergence in Drosophila montana Shows Evidence for Cold Adaptation.

The genomes of species that are ecological specialists will likely contain signatures of genomic adaptation to their niche. However, distinguishing genes related to ecological specialism from other sources of selection and more random changes is a challenge. Here, we describe the genome of Drosophila montana, which is the most extremely cold-adapted Drosophila species known. We use branch tests to identify genes showing accelerated divergence in contrasts between cold- and warm-adapted species and identify about 250 genes that show differences, possibly driven by a lower synonymous substitution rate in cold-adapted species. We also look for evidence of accelerated divergence between D. montana and D. virilis, a previously sequenced relative, but do not find strong evidence for divergent selection on coding sequence variation. Divergent genes are involved in a variety of functions, including cuticular and olfactory processes. Finally, we also resequenced three populations of D. montana from across its ecological and geographic range. Outlier loci were more likely to be found on the X chromosome and there was a greater than expected overlap between population outliers and those genes implicated in cold adaptation between Drosophila species, implying some continuity of selective process at these different evolutionary scales.

Differential gene expression is not required for facultative sex allocation: a transcriptome analysis of brain tissue in the parasitoid wasp Nasonia vitripennis.

Whole-transcriptome technologies have been widely used in behavioural genetics to identify genes associated with the performance of a behaviour and provide clues to its mechanistic basis. Here, we consider the genetic basis of sex allocation behaviour in the parasitoid wasp Nasonia vitripennis. Female Nasonia facultatively vary their offspring sex ratio in line with Hamilton's theory of local mate competition (LMC). A single female or 'foundress' laying eggs on a patch will lay just enough sons to fertilize her daughters. As the number of 'foundresses' laying eggs on a patch increases (and LMC declines), females produce increasingly male-biased sex ratios. Phenotypic studies have revealed the cues females use to estimate the level of LMC their sons will experience, but our understanding of the genetics underlying sex allocation is limited. Here, we exposed females to three foundress number conditions, i.e. three LMC conditions, and allowed them to oviposit. mRNA was extracted from only the heads of these females to target the brain tissue. The subsequent RNA-seq experiment confirmed that differential gene expression is not associated with the response to sex allocation cues and that we must instead turn to the underlying neuroscience to reveal the underpinnings of this impressive behavioural plasticity.

A simple and robust real-time qPCR method for the detection of PIK3CA mutations.

PIK3CA mutations are seemingly the most common driver mutations in breast cancer with H1047R and E545K being the most common of these, accounting together for around 60% of all PIK3CA mutations and have promising therapeutic implications. Given the low sensitivity and the high cost of current genotyping methods we sought to develop fast, simple and inexpensive assays for PIK3CA H1047R and E545K mutation screening in clinical material. The methods we describe are based on a real-time PCR including a mutation specific primer combined with a non-productive oligonucleotide which inhibits wild-type amplification and a parallel internal control reaction. We demonstrate consistent detection of PIK3CA H1047R mutant DNA in genomic DNA extracted from frozen breast cancer biopsies, FFPE material or cancer cell lines with a detection sensitivity of approximately 5% mutant allele fraction and validate these results using both Sanger sequencing and deep next generation sequencing methods. The detection sensitivity for PIK3CA E545K mutation was approximately 10%. We propose these methods as simple, fast and inexpensive diagnostic tools to determine PIK3CA mutation status.

Draft Genome Sequence of Magnetovibrio blakemorei Strain MV-1, a Marine Vibrioid Magnetotactic Bacterium.

We report here the genome sequence of Magnetovibrio blakemorei MV-1, a marine vibrioid magnetotactic bacterium with a single polar flagellum. The current assembly consists of 91 contigs with a combined size of 3,638,804 bp (54.3% G+C content). This genome allows for further investigations of the molecular biomineralization mechanisms of magnetosome formation.