跳到正文
cccskills
免费在 GitHub 公开

bioconductor-immunotation

MHC (major histocompatibility complex) molecules are cell surface complexes that present antigens to T cells. The repertoire of antigens presented in a given genetic background largely depends on the sequence of the encoded MHC molecules, and thus, in humans, on the highly variable HLA (human leukocyte antigen) genes of the hyperpolymorphic HLA locus. More than 28,000 different HLA alleles have been reported, with significant differences in allele frequencies between human populations worldwide

中文简介正在准备中。当前显示原文说明。

查看安装方法

包含的文件(1)

  • SKILL.md4.0 KB

SKILL.md(原文)

安装前可以查看交给 AI 代理的指令内容。

immunotation

Workflows

Standard Workflow

library(immunotation)

# 1. Retrieve valid organisms and chain lookup table for humans
organisms <- get_valid_organisms()
human_chains <- retrieve_chain_lookup_table(organism = "human")

# 2. Map HLA alleles to serotypes and NetMHCpan inputs
allele_list <- c("A*01:01:01", "A*02:01:01", "B*39:01:01")
serotypes <- get_serotypes(allele_list, mhc_type = "MHC-I")
mhcpan_input <- get_mhcpan_input(allele_list, mhc_class = "MHC-I")

# 3. Query allele frequencies in global populations
freq_data <- query_allele_frequencies(
  hla_selection = "A*02:01", 
  hla_sample_size_pattern = "bigger_than", 
  hla_sample_size = 10000, 
  standard = "g"
)

Input/Output Note: Inputs are HLA allele lists and population query parameters; outputs are serotypes, tool-compatible HLA strings, and population frequency tables.

When to Use

  • Converting HLA allele names to serotypes using get_serotypes.
  • Formatting HLA alleles for immunoinformatics tools like NetMHCpan using get_mhcpan_input.
  • Mapping HLA alleles to G groups (get_G_group) or P groups (get_P_group).
  • Querying HLA allele and haplotype frequencies from the Allele Frequency Net Database (AFND) using query_allele_frequencies and query_haplotype_frequencies.
  • Visualizing global allele frequency distributions on a world map using plot_allele_frequency.

When NOT to Use

  • For predicting peptide-MHC binding directly (use tools like NetMHCpan; immunotation only formats the inputs).
  • For non-MHC/HLA gene annotations.

Data Requirements

  • Character vectors of HLA allele names in standard WHO nomenclature (e.g., "A*01:01:01").

Key Parameters

  • organism: Species name (e.g., "human", "mouse") used in retrieve_chain_lookup_table.
  • mhc_type: MHC class type ("MHC-I" or "MHC-II") used in get_serotypes.
  • mhc_class: MHC class ("MHC-I" or "MHC-II") used in get_mhcpan_input.
  • hla_selection: HLA allele or allele group to query.
  • hla_sample_size_pattern: Pattern for filtering sample size (e.g., "bigger_than").
  • hla_sample_size: Minimum sample size for population queries.
  • standard: AFND data standard quality filter (e.g., "g" for Gold, "s" for Silver, "b" for Bronze).

Best Practices

  • Use get_valid_organisms to check if a species is supported before building a lookup table.
  • Convert alleles to G or P groups using get_G_group or get_P_group to resolve ambiguous HLA typings.
  • Filter AFND queries using standard = "g" (Gold standard) to ensure high-quality allele frequency data.

Common Pitfalls

  • Querying MHC-II serotypes with incomplete alpha/beta chain annotations, which returns NA. Ensure both chains are annotated in the MHC restriction ontology (MRO).

Alternatives

  • HLAtools: For general HLA data manipulation and analysis.

Citations

  • Robinson J, Barker DJ, Georgiou X et al. IPD-IMGT/HLA Database. Nucleic Acids Research (2020)
  • Gonzalez-Galarza FF, McCabe A, Santos EJ at al. Allele frequency net database (AFND) 2020 update. Nucleic Acids Research (2020)

References

评价

还没有评价。欢迎分享你的使用感受。

评论

无需登录即可发表。欢迎分享使用感受、提问,或说明在哪些环境下无法运行。

还没有评论。来写第一条吧。

提交后会立即公开。请勿填写邮箱、电话等个人信息。URL 最多 2 个。

0 / 1000 字

此表单受 reCAPTCHA 保护,适用 Google 的隐私权政策和服务条款。

同一仓库的技能

The package contains functions for exploratory oligonucleotide array analysis. The dependence on tkWidgets only concerns few convenience functions. 'affy' is fully functional without it.

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

This tool uses the label-tree function from HyPhy to annotate a phylogenetic tree. It allows users to select a subset of leaves using either a regular expression or a list of sequence names, and then apply a specified label to these selecte

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

This package provides a client for the Bioconductor AnnotationHub web resource. The AnnotationHub web resource provides a central location where genomic files (e.g., VCF, bed, wig) and other resources from standard locations (e.g., UCSC, En

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

apeglm provides Bayesian shrinkage estimators for effect sizes for a variety of GLM models, using approximation of the posterior for individual coefficients.

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

ASURAT is a software for single-cell data analysis. Using ASURAT, one can simultaneously perform unsupervised clustering and biological interpretation in terms of cell type, disease, biological process, and signaling pathway activity. Inputting a single-cell RNA-seq data and knowledge-based databases, such as Cell Ontology, Gene Ontology, KEGG, etc., ASURAT transforms gene expression tables into original multivariate tables, termed sign-by-sample matrices (SSMs).

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

AUCell allows to identify cells with active gene sets (e.g. signatures, gene modules...) in single-cell RNA-seq data. AUCell uses the "Area Under the Curve" (AUC) to calculate whether a critical subset of the input gene set is enriched with

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

查看 bioMate-AI 的全部技能

报告此技能的问题