宏基因组+宏转录组联合分析

humann2分析部分

具体运行脚本： 

source /PUBLIC/software/MICRO/share/MetaGenome_pipeline/Meta_Readsmapping_pipeline/HUMAnN2_pipeline_v2.0/bin/source.txt

humann2 --minpath off --gap-fill on --threads 5 --protein-database /PUBLIC/software/MICRO/share/MetaGenome_pipeline/Meta_Readsmapping_pipeline/HUMAnN2_pipeline_v2.0/bin/../database/uniref50 -i /TJPROJ6/NC_BG_SH/shouhou/202402/X101SC23126055-Z01-J001_special/normal0/normal0.fq.gz -o metacyc --output-basename normal0 --taxonomic-profile ../new/result/02.MetaPhlAn/normal0/normal0.txt

humann2_rna_dna_norm脚本运行部分： 

source /PUBLIC/software/MICRO/share/MetaGenome_pipeline/Meta_Readsmapping_pipeline/HUMAnN2_pipeline_v2.0/bin/source.txt

humann2_rna_dna_norm --input_dna DNA_T1D0_genefamilies.tsv --input_rna RNA_T1D0_genefamilies.tsv --output_basename T1D0

Readme

The HUMAnN 2.0 utility script humann2_rna_dna_norm will Laplace (or Witten-Bell) smooth paired RNA and DNA output files and then return the smoothed DNA, smoothed RNA, and relative expression. HUMAnN 2.0实用程序脚本humann2_rna_dna_norm使用Laplace的方法(或Witten-Bell)平滑成对的RNA和DNA注释信息和丰度文件，然后返回平滑的DNA、平滑的RNA和相对表达。

The relative expression of gene i in sample j is the smoothed value of RNA(i,j) divided by the smoothed value of DNA(i,j) and adjusted for differences in sequencing depth. 基因i在样品j中的相对表达量为RNA(i,j)的平滑值除以DNA(i,j)的平滑值，并根据测序深度的差异进行调整。

*-smoothed_dna.tsv DNA平滑值 *-smoothed_rna.tsv RNA平滑值 *-relative_expression.tsv 相对表达丰度