Transcriptome analysis of 20 taxonomically related benzylisoquinoline alkaloid-producing plants.

BACKGROUND: Benzylisoquinoline alkaloids (BIAs) represent a diverse class of plant specialized metabolites sharing a common biosynthetic origin beginning with tyrosine. Many BIAs have potent pharmacological activities, and plants accumulating them boast long histories of use in traditional medicine and cultural practices. The decades-long focus on a select number of plant species as model systems has allowed near or full elucidation of major BIA pathways, including those of morphine, sanguinarine and berberine. However, this focus has created a dearth of knowledge surrounding non-model species, which also are known to accumulate a wide-range of BIAs but whose biosynthesis is thus far entirely unexplored. Further, these non-model species represent a rich source of catalyst diversity valuable to plant biochemists and emerging synthetic biology efforts. RESULTS: In order to access the genetic diversity of non-model plants accumulating BIAs, we selected 20 species representing 4 families within the Ranunculales. RNA extracted from each species was processed for analysis by both 1) Roche GS-FLX Titanium and 2) Illumina GA/HiSeq platforms, generating a total of 40 deep-sequencing transcriptome libraries. De novo assembly, annotation and subsequent full-length coding sequence (CDS) predictions indicated greater success for most species using the Illumina-based platform. Assembled data for each transcriptome were deposited into an established web-based BLAST portal ( www.phytometasyn.ca) to allow public access. Homology-based mining of libraries using BIA-biosynthetic enzymes as queries yielded ~850 gene candidates potentially involved in alkaloid biosynthesis. Expression analysis of these candidates was performed using inter-library FPKM normalization methods. These expression data provide a basis for the rational selection of gene candidates, and suggest possible metabolic bottlenecks within BIA metabolism. Phylogenetic analysis was performed for each of 15 different enzyme/protein groupings, highlighting many novel genes with potential involvement in the formation of one or more alkaloid types, including morphinan, aporphine, and phthalideisoquinoline alkaloids. Transcriptome resources were used to design and execute a case study of candidate N-methyltransferases (NMTs) from Glaucium flavum, which revealed predicted and novel enzyme activities. CONCLUSIONS: This study establishes an essential resource for the isolation and discovery of 1) functional homologues and 2) entirely novel catalysts within BIA metabolism. Functional analysis of G. flavum NMTs demonstrated the utility of this resource and underscored the importance of empirical determination of proposed enzymatic function. Publically accessible, fully annotated, BLAST-accessible transcriptomes were not previously available for most species included in this report, despite the rich repertoire of bioactive alkaloids found in these plants and their importance to traditional medicine. The results presented herein provide essential sequence information and inform experimental design for the continued elucidation of BIA metabolism.

Roles of Thermophiles and Fungi in Bitumen Degradation in Mostly Cold Oil Sands Outcrops.

Oil sands are surface exposed in river valley outcrops in northeastern Alberta, where flat slabs (tablets) of weathered, bitumen-saturated sandstone can be retrieved from outcrop cliffs or from riverbeds. Although the average yearly surface temperature of this region is low (0.7°C), we found that the temperatures of the exposed surfaces of outcrop cliffs reached 55 to 60°C on sunny summer days, with daily maxima being 27 to 31°C. Analysis of the cooccurrence of taxa derived from pyrosequencing of 16S/18S rRNA genes indicated that an aerobic microbial network of fungi and hydrocarbon-, methane-, or acetate-oxidizing heterotrophic bacteria was present in all cliff tablets. Metagenomic analyses indicated an elevated presence of fungal cytochrome P450 monooxygenases in these samples. This network was distinct from the heterotrophic community found in riverbeds, which included fewer fungi. A subset of cliff tablets had a network of anaerobic and/or thermophilic taxa, including methanogens, Firmicutes, and Thermotogae, in the center. Long-term aerobic incubation of outcrop samples at 55°C gave a thermophilic microbial community. Analysis of residual bitumen with a Fourier transform ion cyclotron resonance mass spectrometer indicated that aerobic degradation proceeded at 55°C but not at 4°C. Little anaerobic degradation was observed. These results indicate that bitumen degradation on outcrop surfaces is a largely aerobic process with a minor anaerobic contribution and is catalyzed by a consortium of bacteria and fungi. Bitumen degradation is stimulated by periodic high temperatures on outcrop cliffs, which cause significant decreases in bitumen viscosity.

SnowyOwl: accurate prediction of fungal genes by using RNA-Seq and homology information to select among ab initio models.

BACKGROUND: Locating the protein-coding genes in novel genomes is essential to understanding and exploiting the genomic information but it is still difficult to accurately predict all the genes. The recent availability of detailed information about transcript structure from high-throughput sequencing of messenger RNA (RNA-Seq) delineates many expressed genes and promises increased accuracy in gene prediction. Computational gene predictors have been intensively developed for and tested in well-studied animal genomes. Hundreds of fungal genomes are now or will soon be sequenced. The differences of fungal genomes from animal genomes and the phylogenetic sparsity of well-studied fungi call for gene-prediction tools tailored to them. RESULTS: SnowyOwl is a new gene prediction pipeline that uses RNA-Seq data to train and provide hints for the generation of Hidden Markov Model (HMM)-based gene predictions and to evaluate the resulting models. The pipeline has been developed and streamlined by comparing its predictions to manually curated gene models in three fungal genomes and validated against the high-quality gene annotation of Neurospora crassa; SnowyOwl predicted N. crassa genes with 83% sensitivity and 65% specificity. SnowyOwl gains sensitivity by repeatedly running the HMM gene predictor Augustus with varied input parameters and selectivity by choosing the models with best homology to known proteins and best agreement with the RNA-Seq data. CONCLUSIONS: SnowyOwl efficiently uses RNA-Seq data to produce accurate gene models in both well-studied and novel fungal genomes. The source code for the SnowyOwl pipeline (in Python) and a web interface (in PHP) is freely available from http://sourceforge.net/projects/snowyowl/.

Genomic compartmentalization of gene families encoding core components of metazoan signaling systems.

To investigate the role of gene localization and genome organization in cell-cell signalling and regulation, we mapped the distribution pattern of gene families that comprise core components of intercellular communication networks. Our study is centered on the distinct evolutionarily conserved metazoan signalling pathways that employ proteins in the receptor tyrosine kinase, WNT, hedgehog, NOTCH, Janus kinase/STAT, transforming growth factor beta, and nuclear hormone receptor protein families. Aberrant activity of these signalling pathways is closely associated with the promotion and maintenance of human cancers. The cataloguing and mapping of genes encoding these signalling proteins and comparisons across species has led us to propose that the genome can be subdivided into six genome-wide primary linkage groups (PLGs). PLGs are composed of assemblages of gene families that are often mutually exclusive, raising the possibility of unique functional identities for each group. Examination of the localization patterns of genes with distinct functions in signal transduction demonstrates dichotomous segregation patterns. For example, gene families of cell-surface receptors localize to genomic compartments that are distinct from the locations of their cognate ligand gene families. Additionally, genes encoding negative-acting components of signalling pathways (inhibitors and antagonists) are topologically separated from their positive regulators and other signal transducer genes. We, therefore, propose the existence of conserved genomic territories that encode key proteins required for the proper activity of metazoan signaling and regulatory systems. Disruption in this pattern of topologic genomic organization may contribute to aberrant regulation in hereditary or acquired diseases such as cancer. We further propose that long-range looping genomic regulatory interactions may provide a mechanism favouring the remarkable retention of these conserved gene clusters during chordate evolution.

Metagenomics of hydrocarbon resource environments indicates aerobic taxa and genes to be unexpectedly common.

Oil in subsurface reservoirs is biodegraded by resident microbial communities. Water-mediated, anaerobic conversion of hydrocarbons to methane and CO2, catalyzed by syntrophic bacteria and methanogenic archaea, is thought to be one of the dominant processes. We compared 160 microbial community compositions in ten hydrocarbon resource environments (HREs) and sequenced twelve metagenomes to characterize their metabolic potential. Although anaerobic communities were common, cores from oil sands and coal beds had unexpectedly high proportions of aerobic hydrocarbon-degrading bacteria. Likewise, most metagenomes had high proportions of genes for enzymes involved in aerobic hydrocarbon metabolism. Hence, although HREs may have been strictly anaerobic and typically methanogenic for much of their history, this may not hold today for coal beds and for the Alberta oil sands, one of the largest remaining oil reservoirs in the world. This finding may influence strategies to recover energy or chemicals from these HREs by in situ microbial processes.

ING1 inhibits Twist1 expression to block EMT and is antagonized by the HDAC inhibitor vorinostat.

ING1 is a chromatin targeting subunit of the Sin3a histone deacetylase (HDAC) complex that alters chromatin structure to subsequently regulate gene expression. We find that ING1 knockdown increases expression of Twist1, Zeb 1&2, Snai1, Bmi1 and TSHZ1 drivers of EMT, promoting EMT and cell motility. ING1 expression had the opposite effect, promoting epithelial cell morphology and inhibiting basal and TGF-ß-induced motility in 3D organoid cultures. ING1 binds the Twist1 promoter and Twist1 was largely responsible for the ability of ING1 to reduce cell migration. Consistent with ING1 inhibiting Twist1 expression in vivo, an inverse relationship between ING1 and Twist1 levels was seen in breast cancer samples from The Cancer Genome Atlas (TCGA). The HDAC inhibitor vorinostat is approved for treatment of multiple myeloma and cutaneous T cell lymphoma and is in clinical trials for solid tumours as adjuvant therapy. One molecular target of vorinostat is INhibitor of Growth 2 (ING2), that together with ING1 serve as targeting subunits of the Sin3a HDAC complex. Treatment with sublethal (LD25-LD50) levels of vorinostat promoted breast cancer cell migration several-fold, which increased further upon ING1 knockout. These observations indicate that correct targeting of the Sin3a HDAC complex, and HDAC activity in general decreases luminal and basal breast cancer cell motility, suggesting that use of HDAC inhibitors as adjuvant therapies in breast cancers that are prone to metastasize may not be optimal and requires further investigation.

Generation of connectivity-preserving surface models of multiple sclerosis lesions.

Progression of multiple sclerosis (MS) results in brain lesions caused by white matter inflammation. MS lesions have various shapes, sizes and locations, affecting cognitive abilities of patients to different extents. To facilitate the visualization of the brain lesion distribution, we have developed a software tool to build 3D surface models of MS lesions. This tool allows users to create 3D models of lesions quickly and to visualize the lesions and brain tissues using various visual attributes and configurations. The software package is based on breadth-first search based 3D connected component analysis and a 3D flood-fill based region growing algorithm to generate 3D models from binary or non-binary segmented medical image stacks.

Building generic anatomical models using virtual model cutting and iterative registration.

BACKGROUND: Using 3D generic models to statistically analyze trends in biological structure changes is an important tool in morphometrics research. Therefore, 3D generic models built for a range of populations are in high demand. However, due to the complexity of biological structures and the limited views of them that medical images can offer, it is still an exceptionally difficult task to quickly and accurately create 3D generic models (a model is a 3D graphical representation of a biological structure) based on medical image stacks (a stack is an ordered collection of 2D images). We show that the creation of a generic model that captures spatial information exploitable in statistical analyses is facilitated by coupling our generalized segmentation method to existing automatic image registration algorithms. METHODS: The method of creating generic 3D models consists of the following processing steps: (i) scanning subjects to obtain image stacks; (ii) creating individual 3D models from the stacks; (iii) interactively extracting sub-volume by cutting each model to generate the sub-model of interest; (iv) creating image stacks that contain only the information pertaining to the sub-models; (v) iteratively registering the corresponding new 2D image stacks; (vi) averaging the newly created sub-models based on intensity to produce the generic model from all the individual sub-models. RESULTS: After several registration procedures are applied to the image stacks, we can create averaged image stacks with sharp boundaries. The averaged 3D model created from those image stacks is very close to the average representation of the population. The image registration time varies depending on the image size and the desired accuracy of the registration. Both volumetric data and surface model for the generic 3D model are created at the final step. CONCLUSIONS: Our method is very flexible and easy to use such that anyone can use image stacks to create models and retrieve a sub-region from it at their ease. Java-based implementation allows our method to be used on various visualization systems including personal computers, workstations, computers equipped with stereo displays, and even virtual reality rooms such as the CAVE Automated Virtual Environment. The technique allows biologists to build generic 3D models of their interest quickly and accurately.

Circulating cell-free DNA is predominantly composed of retrotransposable elements and non-telomeric satellite DNA.

Circulating cell-free DNAs (cfDNAs) are DNA fragments which can be isolated from mammalian blood serum or plasma. In order to gain deeper insight into their origin(s), we have characterized the composition of human and cattle cfDNA via large-scale analyses of high-throughput sequencing data. We observed significant differences between the composition of cfDNA in serum/plasma and the corresponding DNA sequence composition of the human genome. Retrotransposable elements and non-telomeric satellite DNA were particularly overrepresented in the cfDNA population, while telomeric satellite DNA was underrepresented. This was consistently observed for human plasma, bovine serum and for the supernatant of human cancer cell cultures. Our results suggest that reverse transcription of retrotransposable elements and secondary-structure formation during the replication of satellite DNA are contributing to the composition of the cfDNA molecules in the mammalian blood stream. We believe that our work is an important step towards the understanding of the biogenesis of cfDNAs and thus may also facilitate the future exploitation of their diagnostic potential.

Evaluation of host-based molecular markers for the early detection of human sepsis.

We have identified 24 molecular markers, based on circulating nucleic acids (CNA) originating from the human genome, which in combination can be used in a quantitative real-time PCR (qPCR) assay to identify the presence of human sepsis, starting two to three days before the first clinical signs develop and including patients who meet the SEPSIS-3 criteria. The accuracy was more than 87 % inside of the same patient cohort for which the markers were developed and up to 81 % in blind studies of patient cohorts which were not included in the marker development. As our markers are host-based, they can be used to capture bacterial as well as fungal sepsis, unlike the current PCR-based tests, which require species-specific primer sets for each organism causing human sepsis. Our assay directly uses an aliquot of cell-free blood as the substrate for the PCR reaction, thus allowing to obtain the diagnostic results in three to four hours after the collection of the blood samples.

Enzymes revolutionize the bioproduction of value-added compounds: From enzyme discovery to special applications.

Competitive sustainable production in industry demands new and better biocatalysts, optimized bioprocesses and cost-effective product recovery. Our review sheds light on the progress made for the individual steps towards these goals, starting with the discovery of new enzymes and their corresponding genes. The enzymes are subsequently engineered to improve their performance, combined in reaction cascades to expand the reaction scope and integrated in whole cells to provide an optimal environment for the bioconversion. Strain engineering using synthetic biology methods tunes the host for production, reaction design optimizes the reaction conditions and downstream processing ensures the efficient recovery of commercially viable products. Selected examples illustrate how modified enzymes can revolutionize future-oriented applications ranging from the bioproduction of bulk-, specialty- and fine chemicals, active pharmaceutical ingredients and carbohydrates, over the conversion of the greenhouse-gas CO2 into valuable products and biocontrol in agriculture, to recycling of synthetic polymers and recovery of precious metals.

Manipulation of fluorescent and colorimetric changes of fluorescein derivatives and applications for sensing silver ions.

Three structurally similar compounds, , bearing a fluorescence chromophore to which are appended morpholine, thiomorpholine and methylpiperazine substituents, display opposite fluorescence responses to pH changes, in contrast to that observed for fluorescein; and have extremely high binding selectivity towards Ag(+) ions and show completely different fluorescent and colorimetric changes upon addition of Ag(+), and the differences are proposed to be associated with different binding modes of and to this metal ion.

An efficient virtual dissection tool to create generic models for anatomical atlases.

We have developed an efficient virtual dissection tool to create generic 3D models for anatomical atlases without the need for artistic drawings. Our custom-developed tool can be used to extract 3D models from 2D medical image stacks, cut the models and align the sub-models. Corresponding 2D medical image portions of the sub-models can then be registered and averaged. In the end, a generic model can be obtained from the averaged 2D images of several subjects. The technique optimizes the functionalities of existing toolkits and the resulting software package will allow biologists to build their atlases more quickly and accurately.

Fast interactive integration of cross-sectional image datasets and surface data for morphometric analysis.

To investigate external facial morphology and cell proliferation patterns and their relationship with cleft lip malformation in mice, we need to compare samples of mice tissue photographs and surface reconstructions from micro-CT scans obtained from mouse embryos. Tissue samples obtained through digital photography are typically misaligned with respect to each other, which prevents further analysis. We have developed a system for fast interactive alignment of these image stacks for volume reconstruction and data visualization and analysis in 3D. The system is designed to work in multiprocessor environments and can utilize an arbitrary number of processors, cutting down significantly the turnaround time and allowing users to quickly process sets of hundreds of high resolution images using a combination of automated and interactive tools. Additional modules are used to reconstruct the shape of the original subject. Our system is interactive, fully scalable and can be applied to any photographic sliced dataset, regardless of subject and reduces significantly the processing time for stack alignment.

Enterobacteriaceae dominate the core microbiome and contribute to the resistome of arugula (Eruca sativa Mill.).

BACKGROUND: Arugula is a traditional medicinal plant and popular leafy green today. It is mainly consumed raw in the Western cuisine and known to contain various bioactive secondary metabolites. However, arugula has been also associated with high-profile outbreaks causing severe food-borne human diseases. A multiphasic approach integrating data from metagenomics, amplicon sequencing, and arugula-derived bacterial cultures was employed to understand the specificity of the indigenous microbiome and resistome of the edible plant parts. RESULTS: Our results indicate that arugula is colonized by a diverse, plant habitat-specific microbiota. The indigenous phyllosphere bacterial community was shown to be dominated by Enterobacteriaceae, which are well-equipped with various antibiotic resistances. Unexpectedly, the prevalence of specific resistance mechanisms targeting therapeutic antibiotics (fluoroquinolone, chloramphenicol, phenicol, macrolide, aminocoumarin) was only surpassed by efflux pump assignments. CONCLUSIONS: Enterobacteria, being core microbiome members of arugula, have a substantial implication in the overall resistome. Detailed insights into the natural occurrence of antibiotic resistances in arugula-associated microorganisms showed that the plant is a hotspot for distinctive defense mechanisms. The specific functioning of microorganisms in this unusual ecosystem provides a unique model to study antibiotic resistances in an ecological context.

Plasticity of a holobiont: desiccation induces fasting-like metabolism within the lichen microbiota.

The role of host-associated microbiota in enduring dehydration and drought is largely unknown. We have used lichens to study this increasingly important problem because they are the organisms that are optimally adapted to reoccurring hydration/dehydration cycles, and they host a defined and persistent bacterial community. The analysis of metatranscriptomic datasets from bacterial communities of the lung lichen (Lobaria pulmonaria (L.) Hoffm.), sampled under representative hydration stages, revealed significant structural shifts and functional specialization to host conditions. The hydrated samples showed upregulated transcription of transport systems, tRNA modification and various porins (Omp2b by Rhizobiales), whereas the desiccated samples showed different functions related to stress adaption prominently. Carbohydrate metabolism was activated under both conditions. Under dry conditions, upregulation of a specialized ketone metabolism indicated a switch to lipid-based nutrition. Several bacterial lineages were involved in a functional transition that was reminiscent of a 'fasting metaorganism'. Similar functional adaptions were assigned to taxonomically unrelated groups, indicating hydration-related specialization of the microbiota. We were able to show that host-associated bacterial communities are well adapted to dehydration by stress protection and changes of the metabolism. Moreover, our results indicate an intense interplay in holobiont functioning under drought stress.

Bluejay 1.0: genome browsing and comparison with rich customization provision and dynamic resource linking.

BACKGROUND: The Bluejay genome browser has been developed over several years to address the challenges posed by the ever increasing number of data types as well as the increasing volume of data in genome research. Beginning with a browser capable of rendering views of XML-based genomic information and providing scalable vector graphics output, we have now completed version 1.0 of the system with many additional features. Our development efforts were guided by our observation that biologists who use both gene expression profiling and comparative genomics gain functional insights above and beyond those provided by traditional per-gene analyses. RESULTS: Bluejay 1.0 is a genome viewer integrating genome annotation with: (i) gene expression information; and (ii) comparative analysis with an unlimited number of other genomes in the same view. This allows the biologist to see a gene not just in the context of its genome, but also its regulation and its evolution. Bluejay now has rich provision for personalization by users: (i) numerous display customization features; (ii) the availability of waypoints for marking multiple points of interest on a genome and subsequently utilizing them; and (iii) the ability to take user relevance feedback of annotated genes or textual items to offer personalized recommendations. Bluejay 1.0 also embeds the Seahawk browser for the Moby protocol, enabling users to seamlessly invoke hundreds of Web Services on genomic data of interest without any hard-coding. CONCLUSION: Bluejay offers a unique set of customizable genome-browsing features, with the goal of allowing biologists to quickly focus on, analyze, compare, and retrieve related information on the parts of the genomic data they are most interested in. We expect these capabilities of Bluejay to benefit the many biologists who want to answer complex questions using the information available from completely sequenced genomes.

Dietary restriction in Drosophila is dependent on mitochondrial efficiency and constrained by pre-existing extended longevity.

We tested the effects of dietary restriction (DR) on the standard w(1118) strain as well as on our previously described Ra and La strains and their reciprocal isogenic 'cybrid' lines containing heterologous nuclear-mitochondrial combinations. The w(1118) and Ra strains, but neither the La strains nor their 'cybrid' derivatives, are DR-inducible. The Ra and La animals are more robust than the w(1118) animals, and the Ra strain shows an upward shift in its DR threshold. The non-inducibility of the La strain suggests that it expresses a constituitive version of the inducible DR phenotype of its predecessor Ra strain. The difference in mitochondrial efficiency observed between the Ra and La mitochondria (Ross, 2000) has real effects on longevity in the adults: the presence of the more efficient La mitochondria enhances the longevity of an otherwise Ra animal; while the maximum longevity inherent in the La nuclear information cannot be fully expressed in the presence of the Ra mitochondria. Despite the absence of any demographic mathematical parameters robustly and uniquely associated with extended longevity, there are some common trends in that longer lived cohorts have a longer period of low and often constant early life mortality rates (comprising the "health span"), leading to a delayed onset of senescence as noted by the late life increase in age-specific mortality rates (comprising the "senescence span"). There is a genetic basis to this phenotype, but the context-dependence of the demographic data suggests that there is not likely to be only one mechanism or pathway involved in the DR response. In addition, the presence of live yeast had systemic strain specific effects such that it increased longevity in the w(1118) animals but decreased it in the Ra and derived lines. Higher density (4-10x) foods yielded a decreased longevity in all strains at the highest level, showing that malnutrition occurs at both low and high caloric levels.

In-plane three-step needle insertion technique for ultrasound-guided continuous femoral nerve block after total knee arthroplasty: a retrospective review of 488 cases.

BACKGROUND: Continuous femoral nerve block (CFNB) improves postoperative analgesia after total knee arthroplasty (TKA). The aim of this study was to investigate the clinical efficacy and complications of our in-plane three-step needle insertion technique that was devised to reduce the risk of direct femoral nerve injury during CFNB in anesthetized patients. METHODS: This retrospective study included 488 patients who had undergone TKA. Ultrasound (US)-guided CFNB was performed under general or spinal anesthesia using an in-plane, three-step needle insertion technique. The success rate and difficulties of catheter placement, clinical efficacy of analgesia, and complications were recorded. RESULTS: Femoral catheters were placed with a 100% success rate. In 488 patients, real-time US imaging revealed easy separation of the fascia iliaca and the femoral nerve following injection of local anesthetic through a Tuohy needle. Verbal numerical rating scale pain scores (0-10) were 2.0 ± 1.2, 3.5 ± 1.9, 3.2 ± 1.7, 2.9 ± 1.3, and 2.5 ± 1.1 at 1, 6, 12, 24 and 48 h postoperatively. No femoral hematoma, femoral abscess, or neurologic complications, including paresthesia or neurologic deficits, were observed during the 8-week follow-up period. CONCLUSIONS: This retrospective study suggests that an in-plane three-step needle insertion technique for CFNB may reduce the risk of femoral nerve injury in anesthetized patients.

Rhizobiales as functional and endosymbiontic members in the lichen symbiosis of Lobaria pulmonaria L.

Rhizobiales (Alphaproteobacteria) are well-known beneficial partners in plant-microbe interactions. Less is known about the occurrence and function of Rhizobiales in the lichen symbiosis, although it has previously been shown that Alphaproteobacteria are the dominating group in growing lichen thalli. We have analyzed the taxonomic structure and assigned functions to Rhizobiales within a metagenomic dataset of the lung lichen Lobaria pulmonaria L. One third (32.2%) of the overall bacteria belong to the Rhizobiales, in particular to the families Methylobacteriaceae, Bradyrhizobiaceae, and Rhizobiaceae. About 20% of our metagenomic assignments could not be placed in any of the Rhizobiales lineages, which indicates a yet undescribed bacterial diversity. SEED-based functional analysis focused on Rhizobiales and revealed functions supporting the symbiosis, including auxin and vitamin production, nitrogen fixation and stress protection. We also have used a specifically developed probe to localize Rhizobiales by confocal laser scanning microscopy after fluorescence in situ hybridization (FISH-CLSM). Bacteria preferentially colonized fungal surfaces, but there is clear evidence that members of the Rhizobiales are able to intrude at varying depths into the interhyphal gelatinous matrix of the upper lichen cortical layer and that at least occasionally some bacteria also are capable to colonize the interior of the fungal hyphae. Interestingly, the gradual development of an endosymbiotic bacterial life was found for lichen- as well as for fungal- and plant-associated bacteria. The new tools to study Rhizobiales, FISH microscopy and comparative metagenomics, suggest a similar beneficial role for lichens than for plants and will help to better understand the Rhizobiales-host interaction and their biotechnological potential.