Name: Assisted Selection of Biomarkers by Linear Discriminant Analysis Effect Size (LEfSe) in Microbiome Data (Video) | JoVE Journal
Uploaded: 2022-05-16T14:00:00.000+00:00
Duration: 4 min 57 s
Description: Learn LEfSe analysis for biomarker discovery in microbiome data using linear discriminant analysis. This method identifies significant genomic features that differentiate biological groups for research applications.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

This work was supported by a grant from Fundamental Research Funds for the Central Public Welfare Research Institutes (TKS170205) and Foundation for Development of Science and Technology, and Tianjin Research Institute for Water Transport Engineering (TIWTE), M.O.T. (KJFZJJ170201).

NOTE: The protocol was sourced and modified from the research of Segata et al.3. The method is provided at https://bitbucket.org/biobakery/biobakery/wiki/lefse.
1. Preparation of input file for analysis
<ol>
	<li>Prepare the input file (Table 1) of LEfSe, which could be easily generated by many workflows8 or previous protocols9 with the original files (sample file and corresponding species annotation file).</li>
</ol>
2. LEfSe native analysis (limited to the Linux server)
<ol>
	<li>LEfSe Installation 
	NOTE: The LEfSe pipeline is recommended to be installed with Conda10.
	<ol>
		<li>Run the following commands to exclude the possibility of dependencies conflict. Create a conda environment for LEfSe (This step is recommended but not required.). -n stands for the environment name. 
		$ conda create -n LEfSe-env</li>
		<li>To activate the LEfSe environment that was created, run: 
		$ source activate LEfSe-env</li>
		<li>To install LEfSe with channel bioBakery where -c stands for channel name, run: 
		$ conda install -c biobakery lefse</li>
	</ol>
	</li>
	<li>Format data for LEfSe
	<ol>
		<li>Run the following command to format the original file to the internal format for LEfSe. Table.txt is the input file and Table-reformat.in is the output file. -c is used to set the feature, which is used as class (default 1) and -o is used to set the normalization value (default -1.0 meaning no normalization). 
		$ format_input.py Table.txt Table-reformat.in -c 1 -o 1000000</li>
	</ol>
	</li>
	<li>Calculation of linear discriminant analysis (LDA) effect size
	<ol>
		<li>Run the following command. The purpose of this step is to perform LDA of the previous result and generate the result file for the visualization. Table-reformat.in is generated using the previous step and is used as the input file in this step. Table-reformat.res is the result file. 
		$ run_lefse.py Table-reformat.in Table-reformat.res</li>
	</ol>
	</li>
	<li>Visualization by plots
	<ol>
		<li>Plot the LEfSe results. To plot the effect size of the biomarkers in a pdf file,.Table-reformat.res is generated using the previous step and the LDA.pdf is the plot file. &#8211;format is used to set the output file format. 
		$ plot_res.py Table-reformat.res LDA.pdf --format pdf</li>
		<li>Plot the cladogram. To draw the species tree and display the biomarkers in a cladogram. cladogram.pdf is the output file. 
		$ plot_cladogram.py Table-reformat.res cladogram.pdf --format pdf</li>
		<li>Plot one feature (optional) To plot the differences of a single biomarker among different groups. -f is used to set the features of plot. If one was set, the &#8211;feature_name must be given. 
		$ plot_features.py -f one --feature_name &#34;k__Bacteria.p__Firmicutes.c__Bacilli.o__Bacillales&#34; --format pdf Table-reformat.in Table-reformat.res Bacillales.pdf</li>
		<li>Plot the differential features (optional) to draw all the features, but there is too much to be done with caution. --archive is used to choose whether to compress the results. ./ means the path of the results. 
		$ plot_features.py -f diff --archive none --format pdf Table-reformat.in Table-reformat.res ./</li>
	</ol>
	</li>
</ol>
3. LEfSe online analysis (galaxy)
<ol>
	<li>Go to the huttenhower galaxy server11: http://huttenhower.sph.harvard.edu/galaxy.</li>
	<li>Upload the files. Press the Up arrow button on the left pane and upload the file. Click on Choose local file to select the input file and select the format tabular, and then click on the Start button. 
	NOTE: Referred to the webpage (https://bitbucket.org/biobakery/biobakery/wiki/lefse), use the script (taxonomy_summary.R) to generate the input file of LEfSe, and the format (each column with a group name, each line with a different level of annotation separated by &#8220;|&#8221;) is required as shown in Table 1. A schematic overview of the uploading process is shown in Figure 1.</li>
	<li>Format the data for LEfSe. Click on the LEfSe | Format Data for LEfSe link on the left pane, and select the specific rows for class in the file, and click on the Execute button. A schematic overview of the operational process and the parameters used are shown in Figure 2.</li>
	<li>Calculate the LDA effect size. Click on the LEfSe | LDA Effect Size (LEfSe) link on the left pane, and select parameter values according to the analysis requirements. Click on Execute. A schematic overview of the operational process and the parameters used are shown in Figure 3.</li>
	<li>Plot the LEfSe results. Click on the LEfSe | Plot LEfSe Results link on the left pane, and click on the Execute button. A schematic overview of the operational process and the parameters used are shown in Figure 4.</li>
	<li>Plot the cladogram. Click on Plot Cladogram on the left pane, and click on the Execute button after selecting the parameter values. A schematic overview of the operational process and the parameters used are shown in Figure 5.</li>
	<li>Plot one feature by clicking on Plot One Feature on the left pane, and clicking on the Execute button after selecting parameter values. A schematic overview of the operational process and the parameters used are shown in Figure 6.</li>
	<li>Plot differential features by clicking on Plot Differential Features on the left pane, and clcking on the Execute button after selecting parameter values. A schematic overview of the operational process and the parameters used are shown in Figure 7. 
	NOTE: These generated figures can be visualized and downloaded against the resulting output in the right pane.</li>
</ol>

<ol><li>Bolyen, E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Reproducible%2C+interactive%2C+scalable+and+extensible+microbiome+data+science+using+QIIME+2%5BTitle%5D)+AND+%22Nature+Biotechnology%22%5BJournal%5D)">Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2.</a> Nature Biotechnology. 37 (8), 852-857 (2019).</li>
<li>Knight, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Best+practices+for+analysing+microbiomes%5BTitle%5D)+AND+%22Nature+Reviews.+Microbiology%22%5BJournal%5D)">Best practices for analysing microbiomes.</a> Nature Reviews. Microbiology. 16 (7), 410-422 (2018).</li>
<li>Segata, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Metagenomic+biomarker+discovery+and+explanation%5BTitle%5D)+AND+%22Genome+Biology%22%5BJournal%5D)">Metagenomic biomarker discovery and explanation.</a> Genome Biology. 12 (6), 60(2011).</li>
<li>McIver, M., Sayoldin, B., Shafquat, A. Biobakery / lefse [tool]. , Available from: 
 <a target="_blank" href="https://bitbucket.org/biobakery/biobakery/wiki/lefse">https://bitbucket.org/biobakery/biobakery/wiki/lefse</a> (2019).</li>
<li>Kruskal, W. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+nonparametric+test+for+the+several+sample+problem%5BTitle%5D)+AND+%22The+Annals+of+Mathematical+Statistics%22%5BJournal%5D)">A nonparametric test for the several sample problem.</a> The Annals of Mathematical Statistics. 23 (4), 525-540 (1952).</li>
<li>Wilcoxon, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Individual+comparisons+by+ranking+methods%5BTitle%5D)+AND+%22Biometrics+Bulletin%22%5BJournal%5D)">Individual comparisons by ranking methods.</a> Biometrics Bulletin. 1 (6), 80-83 (1945).</li>
<li>Fisher, R. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+use+of+multiple+measurements+in+taxonomic+problems%5BTitle%5D)+AND+%22Annals+of+Eugenics%22%5BJournal%5D)">The use of multiple measurements in taxonomic problems.</a> Annals of Eugenics. 7 (1), 179-188 (1936).</li>
<li>Liu, Y. X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+practical+guide+to+amplicon+and+metagenomic+analysis+of+microbiome+data%5BTitle%5D)+AND+%22Protein+and+Cell%22%5BJournal%5D)">A practical guide to amplicon and metagenomic analysis of microbiome data.</a> Protein and Cell. 41 (7), 1-16 (2020).</li>
<li>Shahi, S. K., Zarei, K., Guseva, N. V., Mangalam, A. K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Microbiota+analysis+using+two-step+PCR+and+next-generation+16S+rRNA+gene+sequencing%5BTitle%5D)+AND+%22Journal+of+Visualized+Experiments%3A+JoVE%22%5BJournal%5D)">Microbiota analysis using two-step PCR and next-generation 16S rRNA gene sequencing.</a> Journal of Visualized Experiments: JoVE. (152), e59980(2019).</li>
<li>Grüning, B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Bioconda%3A+sustainable+and+comprehensive+software+distribution+for+the+life+sciences%5BTitle%5D)+AND+%22Nature+Methods%22%5BJournal%5D)">Bioconda: sustainable and comprehensive software distribution for the life sciences.</a> Nature Methods. 15 (7), 475-476 (2018).</li>
<li>Blankenberg, D., Chilton, J., Coraor, N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Galaxy+external+display+applications%3A+closing+a+dataflow+interoperability+loop%5BTitle%5D)+AND+%22Nature+Methods%22%5BJournal%5D)">Galaxy external display applications: closing a dataflow interoperability loop.</a> Nature Methods. 17 (2), 123-124 (2020).</li>
<li>Langille, M. G. I., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Predictive+functional+profiling+of+microbial+communities+using+16S+rRNA+marker+gene+sequences%5BTitle%5D)+AND+%22Nature+Biotechnology%22%5BJournal%5D)">Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences.</a> Nature Biotechnology. 31 (9), 814-821 (2013).</li>
<li>Shilei, Z., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Reservoir+water+stratification+and+mixing+affects+microbial+community+structure+and+functional+community+composition+in+a+stratified+drinking+reservoir%5BTitle%5D)+AND+%22Journal+of+Environmental+Management%22%5BJournal%5D)">Reservoir water stratification and mixing affects microbial community structure and functional community composition in a stratified drinking reservoir.</a> Journal of Environmental Management. 267, 110456(2020).</li>
</ol>

The authors have nothing to disclose.

Here, the protocol for the identification and characterization of biomarkers within different groups is described. This protocol can easily be adapted for other sample types, such as OTUs of microorganisms. The statistical method by LEfSe can find the characteristic microorganisms in each group (default is LDA &#62;2), that is, the microorganisms that are more abundant in this group relative to the others12. LEfSe is available in both native and web Linux versions where users can also perform LEfS...

Biomarkers are biological characteristics that can be measured and can indicate some phenomena such as infection, disease, or environment. Among them, functional biomarkers may be specific biological functions of single species or common to some species, such as gene, protein, metabolite and pathways. Besides, taxonomic biomarkers indicate an unusual species, a group of organisms (kingdom, phylum, class, order, family, genus, species), the Amplicon Sequence Varient (ASV)1, or the Operational Taxonomic Unit (OTU)2. In order to find biomarkers more quickly and accurately, a tool for analyzing the biological data is necessary. The differences between classes can be explained by LEfSe coupled with standard tests for statistical significance and additional tests encoding biological consistency and effect relevance3. LEfSe is available as a galaxy module, a conda formula, a docker image, and included in bioBakery (VM and cloud)4. Generally, the analysis of microbial diversity often uses a non-parametric test for the uncertain distribution of a sample community. The rank sum test is a non-parametric test method, which uses the rank of samples to replace the value of samples. According to the difference of sample groups, it can be divided into two samples with the Wilcoxon rank sum test and into multiple samples with the Kruskal-Wallis test5,6. Notably, when there are significant differences among multiple groups of samples, a rank-sum test of pairwise comparison of multiple samples should be performed. LDA (which stands for Linear Discriminant Analysis) invented by Ronald Fisher in 1936, is a type of supervised learning, also known as Fisher&#8217;s Linear Discriminant7. It is a classic and popular algorithm in the current field of machine learning data mining.
Here, the LEfSe assay has been optimized by Conda and Galaxy servers. Three groups of 16S rRNA gene sequences are analyzed to demonstrate the significant differences between different groups with LDA scores of microbial communities and visualization results.

<table><tbody><tr><td>No materials used</td><td></td><td></td><td></td></tr></tbody></table>

assisted selection of biomarkers by linear discriminant analysis effect size (lefse) in microbiome data

Assisted Selection of Biomarkers by Linear Discriminant Analysis Effect Size (LEfSe) in Microbiome Data

There is growing attention toward closed biological genomes in the environment and in health. To explore and reveal the intergroup differences among different samples or environments, it is crucial to discover biomarkers with statistical differences among groups. The application of Linear discriminant analysis Effect Size (LEfSe) can help find good biomarkers. Based on the original genome data, quality control, and quantification of different sequences based on taxa or genes are carried out. First, the Kruskal-Wallis rank test was used to distinguish between specific differences among statistical and biological groups. Then, the Wilcoxon rank test was performed between the two groups obtained in the previous step to assess whether the differences were consistent. Finally, a linear discriminant analysis (LDA) was conducted to evaluate the influence of biomarkers on significantly different groups based on LDA scores. To sum up, LEfSe provided the convenience for identifying genomic biomarkers that characterize statistical differences among biological groups.

The LDA scores of microbial communities with significant differences in each group by analyzing the 16S rRNA gene sequences of three samples is shown in Figure 8. The color of the histogram represents different groups, while the length represents the LDA score, which is the influence of the species with significant differences between different groups. The histogram shows the species with significant differences whose LDA score is greater than the preset value. The default preset value is 2....

LEfSe (LDA Effect Size) is a tool for high-dimensional biomarker mining to identify genomic features (such as genes, pathways, and taxonomies) that significantly characterize two or more groups in microbiome data.

There is growing attention toward closed biological genomes in the environment and in health. To explore and reveal the intergroup differences among ...

assisted-selection-biomarkers-linear-discriminant-analysis-effect

Research

JoVE Journal

Genetics

This video describes DNA isolation protocol from stool samples of colon tumor-bearing mice. The extracted DNA is amplified using PCR, and the product is further cleaned using a magnetic beads-based approach. This method can be used to determine the effect of environmental enrichment on colon microbiota and animal mortality.

Stool Microbial DNA Isolation: Studying the Effect of Environmental Enrichment on Microbiome Diversity in Murine Tumor Model

Existing algorithms generate one solution for a biomarker detection dataset. This protocol demonstrates the existence of multiple similarly effective solutions and presents a user-friendly software to help biomedical researchers investigate their datasets for the proposed challenge. Computer scientists may also provide this feature in their biomarker detection algorithms.

Selecting Multiple Biomarker Subsets with Similarly Effective Binary Classification Performances

Described here is a simplified standard operating procedure for microbiome profiling using 16S rRNA metagenomic sequencing and analysis using freely available tools. This protocol will help researchers who are new to the microbiome field as well as those requiring updates on methods to achieve bacterial profiling at a higher resolution.

Microbiota Analysis Using Two-step PCR and Next-generation 16S rRNA Gene Sequencing

This protocol provides a simple, cost-effective DNA isolation method for the analysis of murine gut microbiota alterations during the development of autoimmune disease.

Analysis of Fecal Microbiota Dynamics in Lupus-Prone Mice Using a Simple, Cost-Effective DNA Isolation Method

This protocol shows how to perform cytoplasmic microinjection in farm animal zygotes. This technique can be used to deliver any solution into the one-cell embryo such as genome editing tools to generate knockout animals.

Laser-assisted Cytoplasmic Microinjection in Livestock Zygotes

Galaxy and DAVID have emerged as popular tools that allow investigators without bioinformatics training to analyze and interpret RNA-Seq data. We describe a protocol for C. elegans researchers to perform RNA-Seq experiments, access and process the dataset using Galaxy and obtain meaningful biological information from the gene lists using DAVID.

Transcriptomic Analysis of C. elegans RNA Sequencing Data Through the Tuxedo Suite on the Galaxy Project

This manuscript describes the experimental procedure and software analysis for a bidirectional integration site assay that can simultaneously analyze upstream and downstream vector-host junction DNA. Bidirectional PCR products can be used for any downstream sequencing platform. The resulting data are useful for a high-throughput, quantitative comparison of integrated DNA targets.

Bidirectional Retroviral Integration Site PCR Methodology and Quantitative Data Analysis Workflow

Recent advances in the ability to genetically manipulate somatic cell lines hold great potential for basic and applied research. Here, we present two approaches for CRISPR/Cas9 generated knockout production and screening in mammalian cell lines, with and without the use of selectable markers.

Selection-dependent and Independent Generation of CRISPR/Cas9-mediated Gene Knockouts in Mammalian Cells

This protocol presents an approach for whole transcriptome analysis from zebrafish embryos, larvae, or sorted cells. We include isolation of RNA, pathway analysis of RNASeq data, and qRT-PCR-based validation of gene expression changes.

Sample Preparation and Analysis of RNASeq-based Gene Expression Data from Zebrafish

Here, we describe a protocol to obtain amplicon sequence data of soil, rhizosphere, and root endosphere microbiomes. This information can be used to investigate the composition and diversity of plant-associated microbial communities, and is suitable for the use with a wide range of plant species.

Exploring the Root Microbiome: Extracting Bacterial Community Data from the Soil, Rhizosphere, and Root Endosphere

Learn LEfSe analysis for biomarker discovery in microbiome data using linear discriminant analysis. This method identifies significant genomic features that differentiate biological groups for research applications.

Assisted Selection of Biomarkers by Linear Discriminant Analysis Effect Size (LEfSe) in Microbiome Data

Explore More Videos