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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Oviposition but Not Sex Allocation Is Associated with Transcriptomic Changes in Females of the Parasitoid Wasp Nasonia vitripennis.

Linking the evolution of the phenotype to the underlying genotype is a key aim of evolutionary genetics and is crucial to our understanding of how natural selection shapes a trait. Here, we consider the genetic basis of sex allocation behavior in the parasitoid wasp Nasonia vitripennis using a transcriptomics approach. Females allocate offspring sex in line with the local mate competition (LMC) theory. Female-biased sex ratios are produced when one or a few females lay eggs on a patch. As the number of females contributing offspring to a patch increases, less female-biased sex ratios are favored. We contrasted the transcriptomic responses of females as they oviposit under conditions known to influence sex allocation: foundress number (a social cue) and the state of the host (parasitized or not). We found that when females encounter other females on a patch or assess host quality with their ovipositors, the resulting changes in sex allocation is not associated with significant changes in whole-body gene expression. We also found that the gene expression changes produced by females as they facultatively allocate sex in response to a host cue and a social cue are very closely correlated. We expanded the list of candidate genes associated with oviposition behavior in Nasonia, some of which may be involved in fundamental processes underlying the ability to facultatively allocate sex, including sperm storage and utilization.

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.

Quality control of next-generation sequencing data without a reference.

Next-generation sequencing (NGS) technologies have dramatically expanded the breadth of genomics. Genome-scale data, once restricted to a small number of biomedical model organisms, can now be generated for virtually any species at remarkable speed and low cost. Yet non-model organisms often lack a suitable reference to map sequence reads against, making alignment-based quality control (QC) of NGS data more challenging than cases where a well-assembled genome is already available. Here we show that by generating a rapid, non-optimized draft assembly of raw reads, it is possible to obtain reliable and informative QC metrics, thus removing the need for a high quality reference. We use benchmark datasets generated from control samples across a range of genome sizes to illustrate that QC inferences made using draft assemblies are broadly equivalent to those made using a well-established reference, and describe QC tools routinely used in our production facility to assess the quality of NGS data from non-model organisms.

The transcriptomic basis of oviposition behaviour in the parasitoid wasp Nasonia vitripennis.

Linking behavioural phenotypes to their underlying genotypes is crucial for uncovering the mechanisms that underpin behaviour and for understanding the origins and maintenance of genetic variation in behaviour. Recently, interest has begun to focus on the transcriptome as a route for identifying genes and gene pathways associated with behaviour. For many behavioural traits studied at the phenotypic level, we have little or no idea of where to start searching for "candidate" genes: the transcriptome provides such a starting point. Here we consider transcriptomic changes associated with oviposition in the parasitoid wasp Nasonia vitripennis. Oviposition is a key behaviour for parasitoids, as females are faced with a variety of decisions that will impact offspring fitness. These include choosing between hosts of differing quality, as well as making decisions regarding clutch size and offspring sex ratio. We compared the whole-body transcriptomes of resting or ovipositing female Nasonia using a "DeepSAGE" gene expression approach on the Illumina sequencing platform. We identified 332 tags that were significantly differentially expressed between the two treatments, with 77% of the changes associated with greater expression in resting females. Oviposition therefore appears to focus gene expression away from a number of physiological processes, with gene ontologies suggesting that aspects of metabolism may be down-regulated during egg-laying. Nine of the most abundant differentially expressed tags showed greater expression in ovipositing females though, including the genes purity-of-essence (associated with behavioural phenotypes in Drosophila) and glucose dehydrogenase (GLD). The GLD protein has been implicated in sperm storage and release in Drosophila and so provides a possible candidate for the control of sex allocation by female Nasonia during oviposition. Oviposition in Nasonia therefore clearly modifies the transcriptome, providing a starting point for the genetic dissection of oviposition.

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.