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.

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.

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.

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.

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.

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.

Analysis of methanotrophic communities in landfill biofilters using diagnostic microarray.

Biofilters operated for the microbial oxidation of landfill methane at two sites in Northern Germany were analysed for the composition of their methanotrophic community by means of diagnostic microarray targeting the pmoA gene of methanotrophs. The gas emitted from site Francop (FR) contained the typical principal components (CH4, CO2, N2) only, while the gas at the second site Müggenburger Strasse (MU) was additionally charged with non-methane volatile organic compounds (NMVOCs). Methane oxidation activity measured at 22 degrees C varied between 7 and 103 microg CH4 (g dw)(-1) h(-1) at site FR and between 0.9 and 21 microg CH4 (g dw)(-1) h(-1) at site MU, depending on the depth considered. The calculated size of the active methanotrophic population varied between 3 x 10(9) and 5 x 10(11) cells (g dw)(-1) for biofilter FR and 4 x 10(8) to 1 x 10(10) cells (g dw)(-1) for biofilter MU. The methanotrophic community in both biofilters as well as the methanotrophs present in the landfill gas at site FR was strongly dominated by type II organisms, presumably as a result of high methane loads, low copper concentration and low nitrogen availability. Within each biofilter, community composition differed markedly with depth, reflecting either the different conditions of diffusive oxygen supply or the properties of the two layers of materials used in the filters or both. The two biofilter communities differed significantly. Type I methanotrophs were detected in biofilter FR but not in biofilter MU. The type II community in biofilter FR was dominated by Methylocystis species, whereas the biofilter at site MU hosted a high abundance of Methylosinus species while showing less overall methanotroph diversity. It is speculated that the differing composition of the type II population at site MU is driven by the presence of NMVOCs in the landfill gas fed to the biofilter, selecting for organisms capable of co-oxidative degradation of these compounds.

Diversity of the active methanotrophic community in acidic peatlands as assessed by mRNA and SIP-PLFA analyses.

The active methanotroph community was investigated for the first time in heather (Calluna)-covered moorlands and Sphagnum/Eriophorum-covered UK peatlands. Direct extraction of mRNA from these soils facilitated detection of expression of methane monooxygenase genes, which revealed that particulate methane monooxygenase and not soluble methane monooxygenase was probably responsible for CH(4) oxidation in situ, because only pmoA transcripts (encoding a subunit of particulate methane monooxygenase) were readily detectable. Differences in methanotroph community structures were observed between the Calluna-covered moorland and Sphagnum/Eriophorum-covered gully habitats. As with many other Sphagnum-covered peatlands, the Sphagnum/Eriophorum-covered gullies were dominated by Methylocystis. Methylocella and Methylocapsa-related species were also present. Methylobacter-related species were found as demonstrated by the use of a pmoA-based diagnostic microarray. In Calluna-covered moorlands, in addition to Methylocella and Methylocystis, a unique group of peat-associated type I methanotrophs (Gammaproteobacteria) and a group of uncultivated type II methanotrophs (Alphaproteobacteria) were also found. The pmoA sequences of the latter were only distantly related to Methylocapsa and also to the RA-14 group of methanotrophs, which are believed to be involved in oxidation of atmospheric concentrations of CH(4). Soil samples were also labelled with (13)CH(4), and subsequent analysis of the (13)C-labelled phospholipid fatty acids (PLFAs) showed that 16:1 omega 7, 18:1 omega 7 and 18:1 omega 9 were the major labelled PLFAs. The presence of (13)C-labelled 18:1 omega 9, which was not a major PLFA of any extant methanotrophs, indicated the presence of novel methanotrophs in this peatland.

Revealing the uncultivated majority: combining DNA stable-isotope probing, multiple displacement amplification and metagenomic analyses of uncultivated Methylocystis in acidic peatlands.

Peatlands represent an enormous carbon reservoir and have a potential impact on the global climate because of the active methanogenesis and methanotrophy in these soils. Uncultivated methanotrophs from seven European peatlands were studied using a combination of molecular methods. Screening for methanotroph diversity using a particulate methane monooxygenase-based diagnostic gene array revealed that Methylocystis-related species were dominant in six of the seven peatlands studied. The abundance and methane oxidation activity of Methylocystis spp. were further confirmed by DNA stable-isotope probing analysis of a sample taken from the Moor House peatland (England). After ultracentrifugation, (13)C-labelled DNA, containing genomic DNA of these Methylocystis spp., was separated from (12)C DNA and subjected to multiple displacement amplification (MDA) to generate sufficient DNA for the preparation of a fosmid metagenomic library. Potential bias of MDA was detected by fingerprint analysis of 16S rRNA using denaturing gradient gel electrophoresis for low-template amplification (0.01 ng template). Sufficient template (1-5 ng) was used in MDA to circumvent this bias and chimeric artefacts were minimized by using an enzymatic treatment of MDA-generated DNA with S1 nuclease and DNA polymerase I. Screening of the metagenomic library revealed one fosmid containing methanol dehydrogenase and two fosmids containing 16S rRNA genes from these Methylocystis-related species as well as one fosmid containing a 16S rRNA gene related to that of Methylocella/Methylocapsa. Sequencing of the 14 kb methanol dehydrogenase-containing fosmid allowed the assembly of a gene cluster encoding polypeptides involved in bacterial methanol utilization (mxaFJGIRSAC). This combination of DNA stable-isotope probing, MDA and metagenomics provided access to genomic information of a relatively large DNA fragment of these thus far uncultivated, predominant and active methanotrophs in peatland soil.

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.

RNA isolation from soil for bacterial community and functional analysis: evaluation of different extraction and soil conservation protocols.

The impact of three different RNA isolation methods on the community analysis of metabolically active bacteria was determined by reverse transcription (RT) and PCR amplification of 16S rRNA genes and subsequent terminal restriction fragment length polymorphism (T-RFLP) analysis. Furthermore, soil samples were stored at different conditions in order to evaluate the effect of soil conservation methods on the outcome of the population analysis. The quality of mRNA was assessed by reverse transcription and PCR amplification of eubacterial glutamine synthetase genes. Our results indicated that the community composition as well as the abundance of individual members were affected by the kind of RNA isolation method. Furthermore, the extraction method influenced the recovery of mRNA. Lyophilization, storage at -20 degrees C as well as storage in glycerol stocks at -80 degrees C proved to be equally appropriate for the storage of soils and subsequent RNA isolation.

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.

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 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.