2025年大学《生物信息学》专业题库- 生物信息学在疾病易感基因鉴定中的作用

2025年大学《生物信息学》专业题库——生物信息学在疾病易感基因鉴定中的作用考试时间：______分钟总分：______分姓名：______一、简述利用全基因组关联分析（GWAS）研究疾病易感基因的基本流程，并说明在每个步骤中需要注意的关键质量控制点。二、比较全基因组关联分析（GWAS）和孟德尔随机化（MR）在疾病易感基因研究中的应用特点和主要区别。请分别说明两种方法试图解决的核心问题。三、在进行外显子组测序数据的变异检测时，通常需要进行多种类型的质量控制（QC）过滤。请列举至少五种重要的QC过滤指标，并简要说明每个指标过滤的意义。四、假设你进行了一项针对某种复杂疾病的GWAS研究，并在某个染色体区域检测到一个显著的关联信号。请描述你接下来会采取哪些生物信息学步骤来验证和注释这个关联位点，以推断其可能的功能和与疾病的关联机制。五、全外显子组测序（WES）相比全基因组测序（WGS）具有成本较低、覆盖度较高等优点，但也存在局限性。请列举WES在疾病易感基因鉴定中的主要局限性，并思考如何利用生物信息学方法来部分弥补这些局限性。六、简述生物信息学工具ANNOVAR或SnpEff在处理基因变异数据时，主要提供哪些类型的注释信息？选择其中一类信息，详细说明其对于理解变异潜在功能的重要性。七、在利用机器学习方法预测疾病易感基因或进行疾病分类时，通常需要构建特征集。请结合疾病易感基因鉴定的场景，列举至少三种可以用于构建特征的数据类型，并简要说明如何将这些数据转化为机器学习模型可识别的格式。八、设计一个研究方案，用于探索某种罕见遗传病的潜在易感基因。该疾病由小样本家系数据（假设包含10个家系，每个家系有2-3个患者和若干正常亲属）提供。请说明你会采用哪些生物信息学方法进行分析，并简要描述分析流程中的关键步骤和考虑因素。九、描述进行通路富集分析的基本思路。假设通过GWAS分析获得了一组与某种疾病显著关联的基因，请说明你会如何利用生物信息学工具对这些基因进行通路富集分析，并解释分析结果的生物学意义。十、孟德尔随机化（MR）分析要求选择的遗传变异作为工具变量（IV）必须满足三个核心假设。请列出这三个假设，并简要解释违反其中任何一个假设可能导致MR估计结果产生偏倚的方向。试卷答案一、GWAS基本流程：①样本采集与表型数据收集；②基因组DNA提取与高质量检测；③基因组测序（产生原始测序数据）；④数据预处理（质量控制QC，如去除低质量读段、去除重复序列、对参考基因组进行比对）；⑤变异检测（识别SNP、InDel等变异）；⑥变异注释（确定变异位置对应的基因、功能注释等）；⑦质量控制（过滤不符合标准的样本、变异，如Hardy-Weinberg平衡检验、连锁不平衡检验）；⑧统计关联分析（计算每个变异与疾病表型的关联程度，如计算P值）；⑨多重检验校正（如Bonferroni校正、GWASpermutationtest等，控制家族wise错误率）；⑩关联信号聚类与区间合并（识别潜在的风险区域）；⑪候选基因筛选与功能注释（基于关联信号所在的基因区间进行功能富集分析）；⑫结果验证（在独立样本集或通过实验方法验证）。质量控制关键点：主要包括样本QC（去除DNA质量差的样本）、测序数据QC（去除低质量读段、接头序列、重复序列）、比对质量QC（确保高比例读段成功比对、评估比对准确性）、变异检测QC（去除机器学习产生的假阳性变异）、变异注释QC（确保变异注释的准确性）、连锁不平衡（LD）校正前的样本和变异过滤、以及关联分析前的多重检验校正等。二、GWAS特点：直接检测全基因组范围内的变异与疾病的关联，可以发现新的、未知的易感基因，尤其适用于复杂疾病；可以绘制基因频率分布图；但成本高，分析复杂，易受人群分层影响，难以直接推断因果关系，且可能存在大量假阳性信号。MR特点：利用已知的遗传变异（作为工具变量IV）作为代理变量来推断暴露因素（如某基因变异）与结局（如疾病）之间的因果关系；可以间接检验因果假说，控制混杂因素；但要求工具变量满足无关联性、独立性、工具变量与结局的关联性三个核心假设，且这些假设往往难以完全满足，可能导致估计偏差；主要适用于检验已有GWAS发现的关联信号背后的潜在因果效应。核心区别：GWAS主要用于发现新的关联信号和潜在候选基因，关注的是“是什么”关联；MR主要用于检验已知的关联信号背后的因果效应，关注的是“为什么”关联。三、QC过滤指标：①测序读段质量分数（ReadQualityScore）：评估每个碱基的测序准确性，低质量碱基可能导致变异误检；②比对率（AlignmentRate）：衡量测序读段成功比对到参考基因组的比例，过低可能意味着样本DNA质量差或文库构建问题；③平均片段长度（AverageInsertSize）：反映文库片段化的大小，与目标区域覆盖度相关；④PCR重复读段比例（DuplicationRate）：高比例重复读段可能源于PCR扩增偏好，影响变异检测；⑤位点缺失率（MissingRate）：衡量样本在某个变异位点信息缺失的比例，过高可能意味着该位点测序困难或变异本身频率极低；⑥Hardy-Weinberg平衡检验（HWEP-value）：用于检测群体样本中基因型频率是否符合HWE预期，偏离HWE的位点可能存在样本混匀错误或选择偏倚；⑦引物二聚体/特异性片段比例（PrimerDimer/SpecificProductRatio）：评估PCR扩增效率，非特异性产物过多会影响定量。过滤意义（以HWE为例）：HWE检验用于检测群体样本中基因型频率是否符合随机婚配的预期。偏离HWE的位点通常提示存在样本混匀错误（如家系样本错误分池）、选择偏倚（如疾病选择导致某些基因型频率改变）或群体结构影响。过滤掉这些位点可以避免它们对关联分析的假阳性影响，保证样本质量和群体代表性。四、验证与注释步骤：①确认关联位点：精确确定显著关联SNP的位置及其周围一定范围内的基因组区域（如±1Mb或±500kb）；②检索数据库：利用数据库如dbSNP、Ensembl、UCSCGenomeBrowser等，确认该SNP的注释信息，包括其是否位于基因编码区（外显子）、非编码区（如内含子、启动子、5'UTR、3'UTR）、调控元件等；③预测功能影响：如果SNP位于外显子区，利用工具如SIFT、PolyPhen-2、MutationTaster等预测该SNP错义、无义或移码突变的可能性及其潜在的致病变异（missensevariant）；④查询表达数据：检索GEO等数据库中该基因在不同组织或细胞类型中的表达模式，判断该变异是否可能影响基因表达；⑤查询通路与功能：利用KEGG、GO、Reactome等数据库，查询该基因参与的生物学通路和功能注释，结合疾病相关知识，推断变异可能的功能影响及其与疾病的关联机制；⑥考虑连锁不平衡：检查该SNP与其他已知的功能性SNP或表达量变异（eQTL）是否存在强烈的连锁不平衡（LD），如果存在，可能需要关注LD块中其他更具功能意义的变异。五、WES局限性：①无法覆盖基因组非编码区域：大部分疾病相关变异可能位于外显子区之外的非编码区域（如启动子、增强子、长链非编码RNA区域），WES无法直接检测这些变异，可能导致遗漏重要易感基因；②无法检测所有类型的变异：WES主要关注编码区的变异，对于大片段缺失、重复序列、结构变异等检测能力有限；③数据稀疏性：对于基因组中重复序列较多或覆盖度不均的区域，测序深度可能不足，影响变异检测的准确性和全面性；④成本限制：虽然比WGS便宜，但对于全基因组分析，成本仍然较高。生物信息学弥补方法：①整合WGS数据：对于关键的候选基因或区域，可以补充进行全基因组测序（WGS），以捕获WES遗漏的非编码区变异或大片段结构变异；②利用计算预测：利用生物信息学工具，基于序列特征、进化保守性、已知调控元件位置等预测非编码区域的潜在功能变异；③多组学数据整合：整合表型数据、转录组数据（RNA-Seq）、蛋白质组数据等，利用机器学习等方法推断非编码变异的功能影响；④利用数据库资源：利用公共数据库中已知的非编码变异信息或调控元件注释进行补充注释和功能分析。六、ANNOVAR/SnpEff注释信息类型：①变异类型（VariantType）：如SNP（单核苷酸多态性）、InDel（插入缺失）、CNV（拷贝数变异）、SV（结构变异）等，区分不同类型变异；②基因组版本与坐标（GenomeBuild&Coordinates）：提供变异在参考基因组版本（如GRCh38）上的精确位置（染色体、起始位置、终止位置）；③参考碱基与变异碱基（Reference&AlternateAllele）：明确显示变异位点上的原始碱基和变异产生的碱基；④所在的基因/转录本信息（Gene/TranscriptID&Name）：标注变异影响的具体基因ID和名称，以及可能影响的转录本ID和名称；⑤变异位置（如外显子、内含子、5'UTR、3'UTR、编码区、非编码区）：指示变异发生在基因结构的哪个区域；⑥预测的功能影响（PredictedFunctionalEffect）：预测变异可能对基因功能产生的影响，如错义突变（missense）、无义突变（nonsense）、移码突变（frameshift）、沉默突变（silent）、剪接位点突变（splicesite）、启动子突变等；⑦致病性预测（DiseasePrediction）：一些工具会结合预测算法，评估变异作为致病基因的可能性（如SIFT、PolyPhen-2打分）；⑧通路与功能注释（Pathway&FunctionAnnotations）：提供变异所在基因涉及的生物学通路、分子功能等信息。重要性（以功能影响为例）：预测变异的功能影响对于理解变异与疾病的潜在关联至关重要。例如，预测为错义突变的SNP可能改变氨基酸序列，影响蛋白质结构和功能；预测为剪接位点突变的InDel可能影响mRNA的剪接过程，导致蛋白质表达异常或功能丧失。通过功能注释，可以将基因变异与具体的生物学过程联系起来，为后续的实验验证和机制研究提供方向。七、可用于构建特征的数据类型：①基因组数据：如基因型数据（SNP、InDel、CNV等）、表观遗传修饰数据（如甲基化水平）、基因组距离信息（如与已知位点或基因的距离）；②转录组数据：如基因或转录本的表达量（FPKM/TPM）、信使RNA（mRNA）丰度、小干扰RNA（siRNA）丰度；③蛋白质组数据：如蛋白质丰度、蛋白质修饰状态、蛋白质-蛋白质相互作用（PPI）网络信息；④临床表型数据：如年龄、性别、疾病严重程度、病程、家族史、生活方式（吸烟、饮酒等）、药物反应等；⑤生物通路与功能数据：如基因参与的通路信息、基因本体（GO）功能注释、KEGG通路富集分析结果等。数据转化格式：将上述数据转化为机器学习模型可识别的格式通常涉及以下步骤：①数据标准化/归一化：由于不同来源的数据量纲和数值范围差异很大，需要对特征数据进行标准化或归一化处理，如使用Z-score标准化、Min-Max归一化等，使不同特征具有可比性；②特征编码：对于分类特征（如性别、疾病分型），需要进行编码转换，如独热编码（One-HotEncoding）、标签编码（LabelEncoding）等；③矩阵构建：将处理后的特征数据组织成特征矩阵，每行代表一个样本，每列代表一个特征；④数据划分：将特征矩阵划分为训练集和测试集（或训练集、验证集、测试集），用于模型训练和评估；⑤保存格式：将处理好的特征矩阵保存为模型兼容的格式，如CSV、TXT、HDF5或直接加载到特定的机器学习库（如scikit-learn）中。八、研究方案设计：①数据准备：对提供的家系数据进行QC，包括检查DNA质量、确保样本无交叉污染、进行基因组捕获/测序的QC（比对率、覆盖度、重复率等）；②变异检测：进行家系样本的全外显子组或目标区域测序数据的变异检测，包括SNP和InDel检测、质量控制（去除低质量变异、过滤HWE不平衡变异、限制样本间LD强度）；③家系关系构建与合并：利用家系样本的表型信息，构建家系关系图，识别家系内部关联；合并来自不同家系的遗传数据；④连锁不平衡校正：在家系数据中，采用特定的连锁不平衡校正方法（如Merlin软件考虑家系结构）或无相关性的方法（如使用GCTA估计相关系数）；⑤关联分析：进行基于家系结构的关联分析（如使用Merlin进行FBAT分析或家系连锁不平衡调整的GWAS），或者进行非家系相关的关联分析（如使用PLINK进行家系数据的GWAS分析，考虑家系结构）；⑥候选基因筛选：根据关联分析结果，筛选出P值显著或效应量较大的变异位点；对这些位点所在的基因进行注释，重点关注位于外显子区的变异；⑦功能注释与通路分析：对候选基因进行功能注释（如GO、KEGG），进行通路富集分析，以揭示潜在的功能通路和生物学机制；⑧结果验证：考虑在更大规模的独立样本集中验证初步发现的关联信号；或通过细胞实验、动物模型等实验手段验证候选基因的功能。九、通路富集分析基本思路：①候选基因集合：从GWAS分析中获得一组与疾病显著关联的候选基因列表；②基因功能注释：利用生物信息学工具，为候选基因列表中的每个基因提供功能注释信息，如所属的通路、参与的生物学过程、分子功能等；③计算富集指标：基于基因功能注释信息，计算候选基因集合在特定功能集（如KEGG通路、GO术语）中的富集程度，常用的指标包括富集比（EnrichmentRatio）、超几何检验的P值、Fisher精确检验的P值等；④多重检验校正：由于同时测试大量通路或功能集，需要进行多重检验校正（如Bonferroni校正、FDR校正），以控制假阳性发现率