Molecular signatures of alternative reproductive strategies in a facultatively social hover wasp.

Social insect reproductives and non-reproductives represent ideal models with which to understand the expression and regulation of alternative phenotypes. Most research in this area has focused on the developmental regulation of reproductive phenotypes in obligately social taxa such as honey bees, while relatively few studies have addressed the molecular correlates of reproductive differentiation in species in which the division of reproductive labour is established only in plastic dominance hierarchies. To address this knowledge gap, we generate the first genome for any stenogastrine wasp and analyse brain transcriptomic data for non-reproductives and reproductives of the facultatively social species Liostenogaster flavolineata, a representative of one of the simplest forms of social living. By experimentally manipulating the reproductive 'queues' exhibited by social colonies of this species, we show that reproductive division of labour in this species is associated with transcriptomic signatures that are more subtle and variable than those observed in social taxa in which colony living has become obligate; that variation in gene expression among non-reproductives reflects their investment into foraging effort more than their social rank; and that genes associated with reproductive division of labour overlap to some extent with those underlying division of labour in the separate polistine origin of wasp sociality but only explain a small portion of overall variation in this trait. These results indicate that broad patterns of within-colony transcriptomic differentiation in this species are similar to those in Polistinae but offer little support for the existence of a strongly conserved 'toolkit' for sociality.

Mycoparasitism related targets of Tmk1 indicate stimulating regulatory functions of this MAP kinase in Trichoderma atroviride.

Mycoparasitism is a key feature of Trichoderma (Hypocreales, Ascomycota) biocontrol agents. Recent studies of intracellular signal transduction pathways of the potent mycoparasite Trichoderma atroviride revealed the involvement of Tmk1, a mitogen-activated protein kinase (MAPK), in triggering the mycoparasitic response. We previously showed that mutants missing Tmk1 exhibit reduced mycoparasitic activity against several plant pathogenic fungi. In this study, we identified the most robustly regulated targets that were governed by Tmk1 during mycoparasitism using transcriptome and proteome profiling. Tmk1 mainly exerts a stimulating function for T. atroviride during its mycoparasitic interaction with the fungal plant pathogen Rhizoctonia solani, as reflected by 89% of strongly differently responding genes in the ∆tmk1 mutant compared to the wild type. Specifically, 54% of these genes showed strong downregulation in the response with a deletion of the tmk1 gene, whereas in the wild type the same genes were strongly upregulated during the interaction with the fungal host. These included the gene encoding the mycoparasitism-related proteinase Prb1; genes involved in signal transduction pathways such as a candidate coding for a conserved 14-3-3 protein, and a gene coding for Tmk2, the T. atroviride cell-wall integrity MAP kinase; genes encoding a specific siderophore synthetase, and multiple FAD-dependent oxidoreductases and aminotransferases. Due to the phosphorylating activity of Tmk1, different (phospho-)proteomics approaches were applied and identified proteins associated with cellular metabolism, energy production, protein synthesis and fate, and cell organization. Members of FAD- and NAD/NADP-binding-domain proteins, vesicular trafficking of molecules between cellular organelles, fungal translational, as well as protein folding apparatus were among others found to be phosphorylated by Tmk1 during mycoparasitism. Outstanding downregulation in the response of the ∆tmk1 mutant to the fungal host compared to the wild type at both the transcriptome and the proteome levels was observed for nitrilase, indicating that its defense and detoxification functions might be greatly dependent on Tmk1 during T. atroviride mycoparasitism. An intersection network analysis between the identified transcripts and proteins revealed a strong involvement of Tmk1 in molecular functions with GTPase and oxidoreductase activity. These data suggest that during T. atroviride mycoparasitism this MAPK mainly governs processes regulating cell responses to extracellular signals and those involved in reactive oxygen stress.

Putting hornets on the genomic map.

Hornets are the largest of the social wasps, and are important regulators of insect populations in their native ranges. Hornets are also very successful as invasive species, with often devastating economic, ecological and societal effects. Understanding why these wasps are such successful invaders is critical to managing future introductions and minimising impact on native biodiversity. Critical to the management toolkit is a comprehensive genomic resource for these insects. Here we provide the annotated genomes for two hornets, Vespa crabro and Vespa velutina. We compare their genomes with those of other social Hymenoptera, including the northern giant hornet Vespa mandarinia. The three hornet genomes show evidence of selection pressure on genes associated with reproduction, which might facilitate the transition into invasive ranges. Vespa crabro has experienced positive selection on the highest number of genes, including those putatively associated with molecular binding and olfactory systems. Caste-specific brain transcriptomic analysis also revealed 133 differentially expressed genes, some of which are associated with olfactory functions. This report provides a spring-board for advancing our understanding of the evolution and ecology of hornets, and opens up opportunities for using molecular methods in the future management of both native and invasive populations of these over-looked insects.

Expression Analyses in the Rachis Hint towards Major Cell Wall Modifications in Grape Clusters Showing Berry Shrivel Symptoms.

Berry shrivel (BS) is one of the prominent and still unresolved ripening physiological disorders in grapevine. The causes of BS are unclear, and previous studies focused on the berry metabolism or histological studies, including cell viability staining in the rachis and berries of BS clusters. Herein, we studied the transcriptional modulation induced by BS in the rachis of pre-symptomatic and symptomatic clusters with a custom-made microarray qPCR in relation to a previous RNASeq study of BS berries. Gene set analysis of transcript expression in symptomatic rachis tissue determined suppression of cell wall biosynthesis, which could also be confirmed already in pre-symptomatic BS rachis by CESA8 qPCR analyses, while in BS berries, a high number of SWITCH genes were suppressed at veraison. Additionally, genes associated with the cell wall were differently affected by BS in berries. A high percentage of hydrolytic enzymes were induced in BS grapes in rachis and berries, while other groups such as, e.g., xyloglucan endotransglucosylase/hydrolase, were suppressed in BS rachis. In conclusion, we propose that modulated cell wall biosynthesis and cell wall assembly in pre-symptomatic BS rachis have potential consequences for cell wall strength and lead to a forced degradation of cell walls in symptomatic grape clusters. The similarity to sugar starvation transcriptional profiles provides a link to BS berries, which are low in sugar accumulation. However, further studies remain necessary to investigate the temporal and spatial coordination in both tissues.

Genomic distances reveal relationships of wild and cultivated beets.

Cultivated beets (Beta vulgaris ssp. vulgaris), including sugar beet, rank among the most important crops. The wild ancestor of beet crops is the sea beet Beta vulgaris ssp. maritima. Species and subspecies of wild beets are readily crossable with cultivated beets and are thus available for crop improvement. To study genomic relationships in the genus Beta, we sequence and analyse 606 beet genomes, encompassing sugar beet, sea beet, B. v. adanensis, B. macrocarpa, and B. patula. We observe two genetically distinct groups of sea beets, one from the Atlantic coast and the other from the Mediterranean area. Genomic comparisons based on k-mers identify sea beets from Greece as the closest wild relatives of sugar beet, suggesting that domestication of the ancestors of sugar beet may be traced to this area. Our work provides comprehensive insight into the phylogeny of wild and cultivated beets and establishes a framework for classification of further accessions of unknown (sub-)species assignment.

Quinoa genome assembly employing genomic variation for guided scaffolding.

KEY MESSAGE: We propose to use the natural variation between individuals of a population for genome assembly scaffolding. In today's genome projects, multiple accessions get sequenced, leading to variant catalogs. Using such information to improve genome assemblies is attractive both cost-wise as well as scientifically, because the value of an assembly increases with its contiguity. We conclude that haplotype information is a valuable resource to group and order contigs toward the generation of pseudomolecules. Quinoa (Chenopodium quinoa) has been under cultivation in Latin America for more than 7500 years. Recently, quinoa has gained increasing attention due to its stress resistance and its nutritional value. We generated a novel quinoa genome assembly for the Bolivian accession CHEN125 using PacBio long-read sequencing data (assembly size 1.32 Gbp, initial N50 size 608 kbp). Next, we re-sequenced 50 quinoa accessions from Peru and Bolivia. This set of accessions differed at 4.4 million single-nucleotide variant (SNV) positions compared to CHEN125 (1.4 million SNV positions on average per accession). We show how to exploit variation in accessions that are distantly related to establish a genome-wide ordered set of contigs for guided scaffolding of a reference assembly. The method is based on detecting shared haplotypes and their expected continuity throughout the genome (i.e., the effect of linkage disequilibrium), as an extension of what is expected in mapping populations where only a few haplotypes are present. We test the approach using Arabidopsis thaliana data from different populations. After applying the method on our CHEN125 quinoa assembly we validated the results with mate-pairs, genetic markers, and another quinoa assembly originating from a Chilean cultivar. We show consistency between these information sources and the haplotype-based relations as determined by us and obtain an improved assembly with an N50 size of 1079 kbp and ordered contig groups of up to 39.7 Mbp. We conclude that haplotype information in distantly related individuals of the same species is a valuable resource to group and order contigs according to their adjacency in the genome toward the generation of pseudomolecules.

Comparative genome characterization of the periodontal pathogen Tannerella forsythia.

BACKGROUND: Tannerella forsythia is a bacterial pathogen implicated in periodontal disease. Numerous virulence-associated T. forsythia genes have been described, however, it is necessary to expand the knowledge on T. forsythia's genome structure and genetic repertoire to further elucidate its role within pathogenesis. Tannerella sp. BU063, a putative periodontal health-associated sister taxon and closest known relative to T. forsythia is available for comparative analyses. In the past, strain confusion involving the T. forsythia reference type strain ATCC 43037 led to discrepancies between results obtained from in silico analyses and wet-lab experimentation. RESULTS: We generated a substantially improved genome assembly of T. forsythia ATCC 43037 covering 99% of the genome in three sequences. Using annotated genomes of ten Tannerella strains we established a soft core genome encompassing 2108 genes, based on orthologs present in > = 80% of the strains analysed. We used a set of known and hypothetical virulence factors for comparisons in pathogenic strains and the putative periodontal health-associated isolate Tannerella sp. BU063 to identify candidate genes promoting T. forsythia's pathogenesis. Searching for pathogenicity islands we detected 38 candidate regions in the T. forsythia genome. Only four of these regions corresponded to previously described pathogenicity islands. While the general protein O-glycosylation gene cluster of T. forsythia ATCC 43037 has been described previously, genes required for the initiation of glycan synthesis are yet to be discovered. We found six putative glycosylation loci which were only partially conserved in other bacteria. Lastly, we performed a comparative analysis of translational bias in T. forsythia and Tannerella sp. BU063 and detected highly biased genes. CONCLUSIONS: We provide resources and important information on the genomes of Tannerella strains. Comparative analyses enabled us to assess the suitability of T. forsythia virulence factors as therapeutic targets and to suggest novel putative virulence factors. Further, we report on gene loci that should be addressed in the context of elucidating T. forsythia's protein O-glycosylation pathway. In summary, our work paves the way for further molecular dissection of T. forsythia biology in general and virulence of this species in particular.

Benchmarking of long-read correction methods.

Third-generation sequencing technologies provided by Pacific Biosciences and Oxford Nanopore Technologies generate read lengths in the scale of kilobasepairs. However, these reads display high error rates, and correction steps are necessary to realize their great potential in genomics and transcriptomics. Here, we compare properties of PacBio and Nanopore data and assess correction methods by Canu, MARVEL and proovread in various combinations. We found total error rates of around 13% in the raw datasets. PacBio reads showed a high rate of insertions (around 8%) whereas Nanopore reads showed similar rates for substitutions, insertions and deletions of around 4% each. In data from both technologies the errors were uniformly distributed along reads apart from noisy 5' ends, and homopolymers appeared among the most over-represented kmers relative to a reference. Consensus correction using read overlaps reduced error rates to about 1% when using Canu or MARVEL after patching. The lowest error rate in Nanopore data (0.45%) was achieved by applying proovread on MARVEL-patched data including Illumina short-reads, and the lowest error rate in PacBio data (0.42%) was the result of Canu correction with minimap2 alignment after patching. Our study provides valuable insights and benchmarks regarding long-read data and correction methods.

The Gpr1-regulated Sur7 family protein Sfp2 is required for hyphal growth and cell wall stability in the mycoparasite Trichoderma atroviride.

Mycoparasites, e.g. fungi feeding on other fungi, are prominent within the genus Trichoderma and represent a promising alternative to chemical fungicides for plant disease control. We previously showed that the seven-transmembrane receptor Gpr1 regulates mycelial growth and asexual development and governs mycoparasitism-related processes in Trichoderma atroviride. We now describe the identification of genes being targeted by Gpr1 under mycoparasitic conditions. The identified gene set includes a candidate, sfp2, encoding a protein of the fungal-specific Sur7 superfamily, whose upregulation in T. atroviride upon interaction with a fungal prey is dependent on Gpr1. Sur7 family proteins are typical residents of membrane microdomains such as the membrane compartment of Can1 (MCC)/eisosome in yeast. We found that GFP-labeled Gpr1 and Sfp2 proteins show partly overlapping localization patterns in T. atroviride hyphae, which may point to shared functions and potential interaction during signal perception and endocytosis. Deletion of sfp2 caused heavily altered colony morphology, defects in polarized growth, cell wall integrity and endocytosis, and significantly reduced mycoparasitic activity, whereas sfp2 overexpression enhanced full overgrowth and killing of the prey. Transcriptional activation of a chitinase specific for hyphal growth and network formation and strong downregulation of chitin synthase-encoding genes were observed in Δsfp2. Taken together, these findings imply crucial functions of Sfp2 in hyphal morphogenesis of T. atroviride and its interaction with prey fungi.

Transcriptomic changes in CHO cells after adaptation to suspension growth in protein-free medium analysed by a species-specific microarray.

Chinese Hamster Ovary (CHO) cells are the preferred cell line for production of biopharmaceuticals. These cells are capable to grow without serum supplementation, but drastic changes in their phenotype occur during adaptation to protein-free growth, which typically include the change to a suspension phenotype with reduced growth rate. A possible approach to understand this transformation, with the intention to counteract the reduction in growth by targeted supplementation of protein-free media, is gene expression profiling. The increasing availability of genome-scale data for CHO now facilitates quests for a better understanding of metabolic pathways and gene networks. So far, systematic large-scale expression profiling in CHO cells by microarray was limited due to lack of publicly available array designs and limitations of alternative approaches. Based on the recent release of CHO and Chinese Hamster genome sequences, including an annotated RefSeq genome, we have constructed a publicly available microarray design for effective genome-scale expression profiling. The design employed microarray probes optimized for uniformity, sensitivity, and specificity, with probe properties computed using the latest thermodynamic models. We validated the platform in an analysis of gene expression changes in response to serum-free adaptation. The observed effects on the lipid metabolism as well as on nucleotide synthesis were used to successfully select media supplements that were able to increase growth rate.

Phylloxera (Daktulosphaira vitifoliae Fitch) alters the carbohydrate metabolism in root galls to allowing the compatible interaction with grapevine (Vitis ssp.) roots.

Gall forming phylloxera may compete for nutrients with meristematic tissues and develop heterotrophic structures that act as carbon sinks. In this work, we studied the underlying starch metabolism, sink-source translocation of soluble sugars towards and within root galls. We demonstrated that nodosities store carbohydrates by starch accumulation and monitored the expression of genes involved in the starch metabolic. Thereby we proved that the nodosity is symplastically connected to the source tissues through its development and that the starch metabolism is significantly affected to synthesize and degrade starch within the gall. Genes required for starch biosynthesis and degradation are up-regulated. Among the carbohydrate transporters the expression of a glucose-6-phosphate translocater, one sucrose transporter and two SWEET proteins were increases, whereas hexose transporters, tonoplast monosaccharide transporter and Erd6-like sugar transporters were decreased. We found general evidence for plant response to osmotic stress in the nodosity as previously suggested for gall induction processes. We conclude that nodosities are heterogenous plant organs that accumulate starch to serve as temporary storage structure that is gradually withdrawn by phylloxera. Phylloxera transcriptionally reprograms gall tissues beyond primary metabolism and included downstream secondary processes, including response to osmotic stress.

Assessing technical performance in differential gene expression experiments with external spike-in RNA control ratio mixtures.

There is a critical need for standard approaches to assess, report and compare the technical performance of genome-scale differential gene expression experiments. Here we assess technical performance with a proposed standard 'dashboard' of metrics derived from analysis of external spike-in RNA control ratio mixtures. These control ratio mixtures with defined abundance ratios enable assessment of diagnostic performance of differentially expressed transcript lists, limit of detection of ratio (LODR) estimates and expression ratio variability and measurement bias. The performance metrics suite is applicable to analysis of a typical experiment, and here we also apply these metrics to evaluate technical performance among laboratories. An interlaboratory study using identical samples shared among 12 laboratories with three different measurement processes demonstrates generally consistent diagnostic power across 11 laboratories. Ratio measurement variability and bias are also comparable among laboratories for the same measurement process. We observe different biases for measurement processes using different mRNA-enrichment protocols.

Analysis of methanotroph community structure using a pmoA-based microarray.

The analysis of methanotroph community composition is relevant to studies of methane oxidation in a number of environments where methane is a significant carbon source. The development and application of a microarray targeting the particulate methane monooxygenase gene (pmoA) have allowed a high-throughput, semiquantitative analysis of the major methanotroph groups in a number of different environments. Here we describe the use of a pmoA-based short oligo array for the analysis of methanotroph populations in sediment samples. The method is suitable for analysis of any type of environmental sample from which DNA can be extracted.

Earthworm activity in a simulated landfill cover soil shifts the community composition of active methanotrophs.

Landfills represent a major source of methane in the atmosphere. In a previous study, we demonstrated that earthworm activity in landfill cover soil can increase soil methane oxidation capacity. In this study, a simulated landfill cover soil mesocosm (1 m × 0.15 m) was used to observe the influence of earthworms (Eisenia veneta) on the active methanotroph community composition, by analyzing the expression of the pmoA gene, which is responsible for methane oxidation. mRNA-based pmoA microarray analysis revealed that earthworm activity in landfill cover soil stimulated activity of type I methanotrophs (Methylobacter, Methylomonas, Methylosarcina spp.) compared to type II methanotrophs (particularly Methylocystis spp.). These results, along with previous studies of methanotrophs in landfill cover soil, can now be used to plan in situ field studies to integrate earthworm-induced methanotrophy with other landfill management practises in order to maximize soil methane oxidation and reduce methane emissions from landfills.

Analysis of methanotroph community composition using a pmoA-based microbial diagnostic microarray.

Microbial diagnostic microarrays (MDMs) are highly parallel hybridization platforms containing multiple sets of immobilized oligonucleotide probes used for parallel detection and identification of many different microorganisms in environmental and clinical samples. Each probe is approximately specific to a given group of organisms. Here we describe the protocol used to develop and validate an MDM method for the semiquantification of a range of functional genes--in this case, particulate methane monooxygenase (pmoA)--and we give an example of its application to the study of the community structure of methanotrophs and functionally related bacteria in the environment. The development and validation of an MDM, following this protocol, takes ∼6 months. The pmoA MDM described in detail comprises 199 probes and addresses ∼50 different species-level clades. An experiment comprising 24 samples can be completed, from DNA extraction to data acquisition, within 3 d (12-13 h bench work).

Environmental distribution and abundance of the facultative methanotroph Methylocella.

Methylocella spp. are facultative methanotrophs, which are able to grow not only on methane but also on multicarbon substrates such as acetate, pyruvate or malate. Methylocella spp. were previously thought to be restricted to acidic soils such as peatlands, in which they may have a key role in methane oxidation. There is little information on the abundance and distribution of Methylocella spp. in the environment. New primers were designed, and a real-time quantitative PCR method was developed and validated targeting Methylocella mmoX (encoding the α-subunit of the soluble methane monooxygenase) that allowed the quantification of Methylocella spp. in environmental samples. We also developed and validated specific PCR assays, which target 16S rRNA genes of known Methylocella spp. These were used to investigate the distribution of Methylocella spp. in a variety of environmental samples. It was revealed that Methylocella species are widely distributed in nature and not restricted to acidic environments.

Physical disturbance to ecological niches created by soil structure alters community composition of methanotrophs.

Aggregates of different sizes and stability in soil create a composite of ecological niches differing in terms of physico-chemical and structural characteristics. The aim of this study was to identify, using DNA-SIP and mRNA-based microarray analysis, whether shifts in activity and community composition of methanotrophs occur when ecological niches created by soil structure are physically perturbed. Landfill cover soil was subject to three treatments termed: 'control' (minimal structural disruption), 'sieved' (sieved soil using 2 mm mesh) and 'ground' (grinding using mortar and pestle). 'Sieved' and 'ground' soil treatments exhibited higher methane oxidation potentials compared with the 'control' soil treatment. Analysis of the active community composition revealed an effect of physical disruption on active methanotrophs. Type I methanotrophs were the most active methanotrophs in 'sieved' and 'ground' soil treatments, whereas both Type I and Type II methanotrophs were active in the 'control' soil treatment. The result emphasize that changes to a particular ecological niche may not result in an immediate change to the active bacterial composition and change in composition will depend on the ability of the bacterial communities to respond to the perturbation.

Spatial and temporal diversity of methanotrophs in a landfill cover soil are differentially related to soil abiotic factors.

Methanotrophs present in landfill cover soil can limit methane emissions from landfill sites by oxidizing methane produced in landfill. Understanding the spatial and temporal distribution of populations of methanotrophs and the factors influencing their activity and diversity in landfill cover soil is critical to devise better landfill cover soil management strategies. pmoA-based microarray analyses of methanotroph community structure revealed a temporal shift in methanotroph populations across different seasons. Type II methanotrophs (particularly Methylocystis sp.) were found to be present across all seasons. Minor shifts in type I methanotroph populations were observed. In the case of spatial distribution, only minor differences in methanotroph community structure were observed with no recognizable patterns (both vertical and horizontal) at a 5 m scale. Correlation analysis between soil abiotic parameters (total C, N, NH4 (+) , NO3 (-) and water content) and distribution of methanotrophs revealed a lack of conclusive evidence for any distinct correlation pattern between measured abiotic parameters and methanotroph community structure, suggesting that complex interactions of several physico-chemical parameters shape methanotroph diversity and activity in landfill cover soils.

Identification of active methylotrophic bacteria inhabiting surface sediment of a marine estuary.

Methylotrophs play an essential role in the global carbon cycle due to their participation in methane oxidation and C1 metabolism. Despite this important biogeochemical role, marine and estuarine microorganisms that consume C1 compounds are poorly characterized. In this study, we investigated the diversity of active methylotrophs and methanotrophs in sediment from the Colne Estuary (Brightlingsea, UK). Aerobic surface sediment samples were examined for the presence of C1 -utilizing communities using DNA stable-isotope probing (DNA-SIP) with (13) C-labelled methane, methanol and monomethylamine. Active methylotrophic bacteria were confirmed after DNA-SIP and denaturing gradient gel electrophoresis analyses. Clone libraries of 16S rRNA gene amplicons revealed the presence of methylotrophic bacteria affiliated with Methylophaga spp. in methanol and monomethylamine incubations. The addition of marine ammonium mineral salts medium to the microcosms increased the rate of substrate metabolism in DNA-SIP incubations, although nutrient addition did not affect the active populations contributing (13) C-labelled DNA. The (13) CH4 SIP incubations indicated the predominant activity of type I methanotrophs and microarray hybridization of amplified particulate methane monooxygenase (pmoA) genes confirmed the role of type Ia methanotrophs in SIP incubations. Type II methanotrophs (i.e. Methylocystis and Methylosinus) were only detected in the original sediment and in the unlabelled DNA fractions, which indicated that type II methanotrophs were not actively involved in C1 compound assimilation in DNA-SIP incubations with estuarine surface sediment samples.