CN112469836A - 子宫内膜的miRNA容受性分析 - Google Patents

子宫内膜的miRNA容受性分析 Download PDF

Info

Publication number
CN112469836A
CN112469836A CN202080004080.1A CN202080004080A CN112469836A CN 112469836 A CN112469836 A CN 112469836A CN 202080004080 A CN202080004080 A CN 202080004080A CN 112469836 A CN112469836 A CN 112469836A
Authority
CN
China
Prior art keywords
endometrial
female
receptivity
mirnas
mirna expression
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202080004080.1A
Other languages
English (en)
Other versions
CN112469836B (zh
Inventor
康诗婷
陈伟铭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Inti Biotechnology Co Ltd
Original Assignee
Quark Biosciences Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Quark Biosciences Inc filed Critical Quark Biosciences Inc
Publication of CN112469836A publication Critical patent/CN112469836A/zh
Application granted granted Critical
Publication of CN112469836B publication Critical patent/CN112469836B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6883Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6844Nucleic acid amplification reactions
    • C12Q1/686Polymerase chain reaction [PCR]
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/112Disease subtyping, staging or classification
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/158Expression markers
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/178Oligonucleotides characterized by their use miRNA, siRNA or ncRNA
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Organic Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Engineering & Computer Science (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Genetics & Genomics (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Engineering & Computer Science (AREA)
  • Microbiology (AREA)
  • Molecular Biology (AREA)
  • Immunology (AREA)
  • Biotechnology (AREA)
  • Biophysics (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Pathology (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Biology (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Bioinformatics & Computational Biology (AREA)

Abstract

本公开涉及一种使用来自女性的样本,例如子宫内膜检体来测定子宫内膜状态的方法,所述方法包括:(a)对来自女性的子宫内膜样本执行分析以测定子宫内膜样本的微小核糖核酸(miRNA)表达谱,其中miRNA表达谱包括多个miRNA,例如分别具有SEQ ID NO:1到SEQ ID NO:167的序列的167个miRNA的表达水平;和(b)使用例如计算机算法分析miRNA表达谱以获得容受性预测评分。本公开的另一方面涉及适合执行所述方法的试剂盒,以及使用试剂盒于诊断和治疗目的。

Description

子宫内膜的miRNA容受性分析
相关申请的交叉引用
本申请要求2019年7月2日提交的美国临时专利申请第62/869,574号的权益,其全部内容以引用的方式并入本文中。
技术领域
本公开是有关于一种测定女性的子宫内膜容受性的方法,所述方法使用(a)包括多个miRNA,例如167个miRNA的表达水平的微小核糖核酸(miRNA)表达谱;和(b)计算机算法,所述计算机算法基于miRNA表达谱对女性的子宫内膜状态进行分类。本公开的另一方面是有关于适合执行所述方法的试剂盒,以及使用试剂盒于诊断和治疗目的。在一些实施例中,所述方法和/或试剂盒用于分类女性对体外受精(in vitro fertilization;IVF)疗程的反应性。
背景技术
辅助生殖技术(包含IVF)为具解决无法成功生殖的潜力方法。IVF能否成功的主要因素之一为子宫内膜的容受期状态。子宫内膜接受胚胎着床的时期仅开放相对短的时间,所述时期称为着床窗口期(window of implantation;WOI)。着床窗口期通常出现在大约月经周期的第19天到第21天。长期以来,着床窗口期大部分依赖日历方法计算可能的受孕时间,但此方法往往不可靠,因此需要直接透过检查子宫内膜本身以监测子宫内膜当下的状态,以更可靠的方式预测胚胎植入的机会点。
人类子宫内膜为通过蛋白质和miRNA两者周期性地调节的组织。人类基因组包括超过2500个miRNA,其中一些已被证实在生殖周期中具有调控作用。举例来说,最近的文献指出某些miRNA可以调控与着床窗口期发展过程有关的基因。
传统上,组织学和影像方法可用于评估子宫内膜的状态。然而,上述方法极为耗时,且通常无法清楚地区分子宫内膜的容受期状态与非容受期状态。市场上也陆续研发出利用特定基因表达水平检测的方法,早期研究多集中于特定的生物标记基因。艾捷隆(Igenomix)开发出“子宫内膜容受性分析”(Endometrial Receptivity Analysis;ERA)检测,其利用微阵列芯片分析与子宫内膜容受性有关的特定238个基因表达量。然而,基于微阵列芯片技术平台的ERA检测具有某些缺点。举例来说,众所周知以微阵列芯片分析基因表达量需要较多的组织样本。另外,相较于实时定量聚合酶连锁反应(quantitativepolymerase chain reaction;qPCR)技术,微阵列芯片技术平台通常具有较低的特异性。近期,ERA检测也推出以次世代定序(next-generation sequencing;NGS)技术平台进行子宫内膜容受性分析,然而该技术平台也同样需要较大量的检体组织及以及完整性佳的RNA样品质量进行分析。
综上所述,目前市面上仍需要一种更可靠的子宫内膜容受性检测方法,其可适用于较少量的检体组织或是低质量、低样品量样品,以判定子宫内膜处于容受期状态或非容受期状态。
发明内容
本公开涉及使用来自女性的样本,例如子宫内膜检体来测定子宫内膜容受性的方法,所述方法包括:(a)对来自女性的子宫内膜样本进行分析,以测定子宫内膜样本的miRNA表达谱,其中miRNA表达谱包括多个miRNA,例如分别具有SEQ ID NO:1到SEQ ID NO:167的序列的167个miRNA的表达水平;以及(b)分析miRNA表达谱以获得容受性预测评分,其中容受性预测评分用以判断女性的子宫内膜容受性状态。本公开的另一方面是有关于适合于执行所述方法的试剂盒,以及所述试剂盒用于测定女性的子宫内膜状态。
本公开的某些实施例概括于以下段落中。本列举仅为示范性的且不详列本公开所提供的所有实施例。
实施例1.一种测定子宫内膜状态的方法,其包括:(a)对来自女性的子宫内膜样本进行分析,以测定子宫内膜样本的miRNA表达谱,其中miRNA表达谱包括多个miRNA的表达水平;和(b)分析miRNA表达谱以获得容受性预测评分,其中容受性预测评分对女性的子宫内膜状态进行分类,且其中子宫内膜状态包括容受期前、容受期或容受期后,且其中多个miRNA包括至少50个、75个、100个、125个、150个或200个miRNA,且优先选择分别具有SEQ IDNO:1到SEQ ID NO:167的序列的至少167个miRNA。
实施例2.根据实施例1所述的测定子宫内膜状态的方法,其中子宫内膜样本取自女性的子宫腔。
实施例3.根据实施例1或实施例2所述的测定子宫内膜状态的方法,其中子宫内膜样本包括子宫内膜检体、子宫内膜灌洗液或其组合。
实施例4.根据实施例1至实施例3中的任一项所述的测定子宫内膜状态的方法,其中子宫内膜样本(i)在女性内源性促黄体激素(luteinizing hormone;LH)遽增之后七天或(ii)在女性孕酮(progesterone)投药之后五天取样。
实施例5.根据实施例1至实施例4中的任一项所述的测定子宫内膜状态的方法,其中miRNA表达谱通过qPCR、定序、微阵列芯片或RNA-DNA杂交捕获技术来测定。
实施例6.根据实施例5所述的测定子宫内膜状态的方法,其中miRNA表达谱通过对由子宫内膜样本中的miRNA合成的cDNA进行qPCR来测定。
实施例7.根据实施例6所述的测定子宫内膜状态的方法,其中cDNA合成使用具有由以下通式表示的核苷酸序列的通用逆转录引物执行:5'-R-(dT)nVN-3',其中R包括SEQID NO:168,(dT)n为n个连续胸腺嘧啶残基,其中n为19,V为腺嘌呤残基、鸟嘌呤残基或胞嘧啶残基,且N为腺嘌呤残基、鸟嘌呤残基、胞嘧啶残基或胸腺嘧啶残基。
实施例8.根据实施例1至实施例7中的任一项所述的测定子宫内膜状态的方法,其中容受性预测评分为通过计算机算法产生且使用算式MIRA score=f(X∈eq(C))=Xβ+ε所计算的数值,β为系数向量,且ε为误差。
实施例9.根据实施例8所述的测定子宫内膜状态的方法,其中计算机算法通过执行以下一或多个步骤来建立:数据正规化(data normalization)、数据缩放(datascaling)、数据转换(data transformation)、预测建模(prediction modeling)以及交叉验证(cross-validation)。
实施例10.根据实施例8或实施例9所述的测定子宫内膜状态的方法,其中容受性预测评分大于1指容受期前,容受性预测评分小于-1指容受期后,且容受性预测评分介于-1到1指容受期。
实施例11.根据实施例1至实施例10中的任一项所述的测定子宫内膜状态的方法,其中若子宫内膜状态被测定为处于容受期前或容受期后,则所述方法进一步包括:重复步骤(a)和(b)至少一次或直到子宫内膜状态被测定为处于容受期。
实施例12.根据实施例1至实施例11中的任一项所述的测定子宫内膜状态的方法,其中女性遭受或曾遭受胚胎植入失败。
实施例13.根据实施例1至实施例12中的任一项所述的测定子宫内膜状态的方法,其中女性进行IVF疗程。
实施例14.根据实施例13所述的测定子宫内膜状态的方法,其中容受性预测评分进一步分类女性对IVF疗程的反应性。
实施例15.一种检测女性胚胎植入的子宫内膜容受性的方法,其包括:(a)对来自女性的子宫内膜样本执行分析以测定子宫内膜样本的miRNA表达谱,其中miRNA表达谱包括多个miRNA的表达水平;和(b)分析所述miRNA表达谱以获得容受性预测评分,其中容受性预测评分用以判定女性子宫内膜是否具有适合胚胎植入的容受性状态,且其中多个miRNA包括至少50个、75个、100个、125个、150个或200个miRNA,且优先选择分别具有SEQ ID NO:1到SEQ ID NO:167的序列的至少167个miRNA。
实施例16.根据实施例15所述的检测女性胚胎植入的子宫内膜容受性的方法,其中子宫内膜样本获自女性的子宫腔。
实施例17.根据实施例15或实施例16所述的检测女性胚胎植入的子宫内膜容受性的方法,其中子宫内膜样本包括子宫内膜检体、子宫内膜灌洗液或其组合。
实施例18.根据实施例15至实施例17中的任一项所述的检测女性胚胎植入的子宫内膜容受性的方法,其中子宫内膜样本(i)在女性内源性促黄体激素遽增之后七天或(ii)在女性孕酮投药之后五天获得。
实施例19.根据实施例15至实施例18中的任一项所述的检测女性胚胎植入的子宫内膜容受性的方法,其中miRNA表达谱通过qPCR、定序、微阵列芯片或RNA-DNA杂交捕获技术来测定。
实施例20.根据实施例19所述的检测女性胚胎植入的子宫内膜容受性的方法,其中miRNA表达谱通过对由子宫内膜样本中的miRNA合成的cDNA进行qPCR来测定。
实施例21.根据实施例20所述的检测女性胚胎植入的子宫内膜容受性的方法,其中cDNA合成使用具有由以下通式表示的核苷酸序列的通用逆转录引物执行:5'-R-(dT)nVN-3',其中R包括SEQ ID NO:168,(dT)n为n个连续胸腺嘧啶残基,n为19,V为腺嘌呤残基、鸟嘌呤残基或胞嘧啶残基,且N为腺嘌呤残基、鸟嘌呤残基、胞嘧啶残基或胸腺嘧啶残基。
实施例22.根据实施例15至实施例21中的任一项所述的检测女性胚胎植入的子宫内膜容受性的方法,其中容受性预测评分为通过计算机算法产生且使用算式MIRA score=f(X∈eq(C))=Xβ+ε所计算的数值,β为系数向量,且ε为误差。
实施例23.根据实施例22所述的检测女性胚胎植入的子宫内膜容受性的方法,其中计算机算法通过执行以下一或多个步骤中来建立:数据正规化、数据缩放、数据转换、预测建模以及交叉验证。
实施例24.根据实施例22或实施例23所述的检测女性胚胎植入的子宫内膜容受性的方法,其中容受性预测评分介于-1到1指女性子宫内膜具有适合胚胎植入的容受性状态。
实施例25.根据实施例15至实施例24中的任一项所述的检测女性胚胎植入的子宫内膜容受性的方法,其中女性遭受或曾遭受胚胎植入失败。
实施例26.一种试剂盒,其包括:(a)一或多个针对多个miRNA的miRNA表达谱分析芯片,和(b)关于(i)任选地使用一或多个miRNA表达谱分析芯片测定来自女性的子宫内膜样本的miRNA表达谱和(ii)基于miRNA表达谱使用计算机算法获得容受性预测评分的使用说明,其中多个miRNA包括至少50个、75个、100个、125个、150个或200个miRNA,且优先选择分别具有SEQ ID NO:1到SEQ ID NO:167的序列的至少167个miRNA。
实施例27.根据实施例26所述的试剂盒,其中一或多个miRNA表达谱分析芯片包括用于检测多个miRNA的表达水平的引物。
实施例28.根据实施例27所述的试剂盒,其中miRNA表达谱分析芯片适合于进行实时定量PCR(qPCR)、定序、微阵列芯片或RNA-DNA杂交捕获分析,优先选择qPCR,以检测多个miRNA的表达水平。
实施例29.一种根据实施例27或实施例28所述的试剂盒的用途,其用于测定女性的子宫内膜状态。
实施例30.根据实施例29所述的试剂盒的用途,其中女性遭受或曾遭受胚胎植入失败和/或进行体外受精(IVF)疗程。
附图说明
图1描绘在自然周期或激素替代治疗周期中女性的子宫内膜状态。LH+5:在女性内源性促黄体激素(LH)遽增之后五天;LH+7:在女性内源性LH遽增之后七天;且LH+9:在女性内源性LH遽增之后九天。P+3:在女性孕酮投药之后三天;P+5:在女性孕酮投药之后五天;且P+7:在女性孕酮投药之后七天。
图2描绘根据本公开使用特定167个miRNA的MIRA PanelChip的子宫内膜容受性检测的工作流程。
图3描绘计算机算法(MIRA模型)如何建构且MIRA模型如何产生检测结果的过程。
图4A示出子宫内膜容受性的示范性分析,其将子宫内膜状态分类为以下三种状态中的一种:容受期前状态、容受期状态或容受期后状态。
图4B示出分类为三种容受期状态的示范性女性胚胎植入结果。
图5示出使用183个子宫内膜样本所得的具167个miRNA表达水平的miRNA表达谱的10折(10-fold)交叉验证和妊娠率。SEN:灵敏性=真阳性/(真阳性+假阴性);SPE:特异性=真阴性/(真阴性+假阳性);PPV:精确性或阳性预测值=真阳性/(真阳性+假阳性);且NPV:阴性预测值=真阴性/(真阴性+假阴性)。P+6:在女性孕酮投药之后六天胚胎植入,所述女性的子宫内膜先前被测定为处于容受期前状态;P+5:在女性孕酮投药之后五天胚胎植入,所述女性的子宫内膜先前被测定为处于容受期状态;且P+4.5:在女性孕酮投药之后4.5天(即108个小时)胚胎植入,所述女性的子宫内膜先前被测定为处于容受期后状态。
图6示出MIRA评分系统,其根据容受性预测评分的值而将子宫内膜样本分类为以下三种状态中的一种:容受期前状态、容受期状态或容受期后状态。
具体实施方式
本文所阐述的公开内容和实施例将被理解为仅为示范性的且并不限制本发明的范围。虽然本文采用特定术语,除非另外指出,否则所述术语仅在通用意义和描述性意义上使用且不用于限制目的。
定义
除非上下文另外清楚地指示,否则如本文所用,单数形式“一(a/an)”和“所述(the)”意图同样包含复数形式。
术语“cDNA”是指通过使用逆转录酶对RNA进行逆转录而产生的互补DNA。在一些实施例中,RNA含有自子宫内膜组织样本提取的miRNA。参看实例1。
术语“包括”、“具有”以及“包含”是开放的连系动词。这些动词中的一或多个的任何形式或时态,例如“包括(comprises/comprising)”、“具有(has/having)”以及“包含(includes/including)”也是开放的。举例来说,“包括”、“具有”或“包含”一或多个步骤的任何方法不限于仅具有那些一或多个步骤,且也可涵盖其它未列出的步骤。类似地,“包括”、“具有”或“包含”一个或多个特征的任何组合物或试剂盒不限于仅具有那些一个或多个特征,且可涵盖其它未列出的特征。除非另外要求,否则关于本文的某些实施例提供的任何和所有实例或示范性语言(例如,“如”)的使用仅意图更好地阐明本公开,且不对本公开的范围构成限制。
术语“表达”是指生物样本,例如女性的子宫内膜组织样本中的RNA分子的转录和/或积累。在此上下文中,术语“miRNA表达”是指生物样本中的一或多个miRNA的数量,且可通过使用所属领域中已知的合适方法来检测miRNA表达。参看,例如,实例1。
术语“微小核糖核酸”(“microRNA”或“miRNA”)是指从内源基因衍生的一类长度为大约18个到25个核苷酸的非编码RNA。miRNA通过与其目标mRNA的3'非转译区(UTR)进行碱基配对来作为基因表达的转录后调控因子,以用于mRNA降解或转译抑制。
术语“核酸”、“核苷酸”以及“多核苷酸”可互换地使用且是指呈单链或双链形式的DNA或RNA的聚合物。除非另外指出,否则这些术语涵盖含有天然核苷酸的已知类似物的多核苷酸,所述多核苷酸具有与参考核酸相似的结合特性且以与天然存在的核苷酸相似的方式进行代谢。
术语“引物”是指寡核苷酸,当在诱导引物延伸产物的合成的条件下,例如在核苷酸和聚合诱导剂(如DNA或核糖核酸聚合酶)的存在下且在合适温度、pH、金属离子浓度以及盐浓度下放置所述寡核苷酸时,所述寡核苷酸用以引发互补核酸链的合成。
术语“探针”是指包括多核苷酸的结构,其含有与存在于目标核酸分析物(例如,核酸扩增产物)中的核酸序列互补的核酸序列。探针的多核苷酸区可由DNA和/或RNA和/或合成核苷酸类似物构成。探针的长度通常与其用于专一性检测目标核酸的所有或部分目标序列兼容。
术语“实时定量PCR”(“qPCR”)是指使用聚合酶链反应以同时扩增且定量目标DNA和/或RNA的实验方法。定量使用多种化学物质(包含例如
Figure BDA0002909595370000051
Green的荧光染料或塔克曼(Taqman)探针的荧光报导子寡核苷酸探针)执行,且实时定量通过测量在一或多个扩增周期之后反应中的扩增DNA和/或RNA来执行。
术语“标靶”(targeting)是指选择与所关注核酸序列杂交的合适核苷酸序列。在一些实施例中,所关注核酸序列包含具有SEQ ID NO:1到SEQ ID NO:167中的任一个的序列的miRNA。参看实例1。
用于测定子宫内膜状态的方法的概述
子宫内膜容受性是指女性的子宫内膜准备用于胚胎植入的状态。此发生在称为着床窗口期(WOI)的时间段内的所有月经周期中。如图1所示,在自然周期内,排卵在LH遽增之后发生,且WOI为LH遽增之后大约七天(LH+7)。在激素替代治疗周期中,WOI在孕酮投药之后大约五天(P+5)。这些估计给出了关于子宫内膜容受性的可能信息。然而,子宫内膜状态的最终答案只可通过检查子宫内膜本身来提供。
为此,子宫内膜样本可在激素替代治疗周期中孕酮投药之后五天(P+5)或在自然周期中内源性LH遽增之后七天(LH+7)从女性的子宫腔采集。随后样本以分子诊断工具分析子宫内膜容受性状态。在根据本公开的测定子宫内膜状态的方法中,分子诊断工具分析子宫内膜样本的miRNA表达谱。
如图2中所示,本公开提供测定子宫内膜状态的方法,其包括:(a)对子宫内膜样本执行分析以测定子宫内膜样本的miRNA表达谱,其中miRNA表达谱包括多个miRNA,例如分别具有SEQ ID NO:1到SEQ ID NO:167的序列的167个miRNA的表达水平;和(b)用计算机算法分析miRNA表达谱以获得容受性预测评分,其中容受性预测评分将子宫内膜状态分类为容受期前状态、容受期状态或容受期后状态。
容受期前状态系指子宫内膜尚未准备好接受胚胎且此时胚胎植入可能过早。容受期状态(WOI)系指子宫内膜处于胚胎植入的最佳时间。容受期后状态系指子宫内膜已经过了胚胎植入的最佳阶段。
分析miRNA表达谱以测定子宫内膜容受性
本公开测定子宫内膜样本的miRNA表达谱。在一些实施例中,miRNA表达谱包括多个miRNA,例如至少10个、25个、50个、75个、100个、125个、150个或200个miRNA的表达水平,其全部可能涉及子宫内膜容受性的调控作用。在优先选择的实施例中,本公开提供167个miRNA的选择,其表达水平涉及子宫内膜容受性的调控。参看实例1。通过首先从HumanDisease Ontology数据库中识别与生殖疾病有关的基因,且随后使用miRTARBase、TargetScan以及miRDB选择潜在的调控因子miRNA来选出这167个miRNA。
为测定子宫内膜状态,根据本公开的方法包括执行分析以测定子宫内膜样本的miRNA表达谱,其中miRNA表达谱包括表1中所示的167个miRNA的表达水平。
表1. 167个miRNA的名称和序列。
Figure BDA0002909595370000071
Figure BDA0002909595370000081
Figure BDA0002909595370000091
Figure BDA0002909595370000101
可用所属领域中已知的定量方法来分析miRNA的表达水平。在一些实施例中,为促进分析,一或多个以167个miRNA为标靶的miRNA表达谱分析芯片可作为使用。举例来说,在实例1中,两个miRNA表达谱分析芯片被设计开发以用于分析167个miRNA的表达水平。在一些实施例中,一或多个芯片另外以可作为miRNA表达分析的内源性对照的某些RNA序列为标靶,例如18s rRNA。参看实例1。
本公开提供测定子宫内膜样本的miRNA表达谱的方法。所述方法通常包括(i)获得或已获得来自女性子宫腔的子宫内膜样本,(ii)进行分析以测定子宫内膜样本的miRNA表达谱,其中miRNA表达谱包括多个miRNA,例如分别具有SEQ ID NO:1到SEQ ID NO:167的序列的167个miRNA的表达水平。
在一些实施例中,子宫内膜样本可通过侵入性方法,例如通过获取来自子宫内膜的小量检体来获得。参看实例1。在一些实施例中,子宫内膜样本可通过较不具侵入性的方法,例如通过收集存在于子宫灌洗液中的脱落细胞来获得。不希望因此受任何理论束缚,相信相较于基于微阵列芯片的mRNA表达谱分析方法,所主张的基于qPCR的miRNA表达谱方法提供更高的特异性和灵敏性,使得在根据本公开的方法中可能仅需要显著较少量的子宫内膜样本。参看王(Wang)等人,“Large scale real-time PCR validation on geneexpression measurements from two commercial long-oligonucleotidemicroarrays”,BMC Genomics,2006,7:59-75。
在一些实施例中,子宫内膜样本在女性的内源性LH遽增之后七天(LH+7)获得。在一些实施例中,子宫内膜样本在女性孕酮投药五天之后(P+5)获得。
使用所属领域中已知的方法可提取且富集子宫内膜样本中的miRNA。举例来说,可遵循制造商的说明书使用miRNeasy Micro Kit(凯杰(QIAGEN))从子宫内膜组织提取miRNA。参看实例1。富集的miRNA可在-80℃下储存。可使用所属领域中已知的方法来分析miRNA的数量和质量。举例来说,可使用商购的安捷伦生物分析仪(Agilent bioanalyzer)来分析miRNA。
可通过所属领域中已知的方法,包含qPCR、定序、微阵列芯片或RNA-DNA杂交捕获技术来定量各miRNA的表达水平。在一些实施例中,根据本公开的方法使用qPCR反应,其通常具有比北方印迹杂交和/或微阵列基因芯片分析更高的灵敏性和特异性。为此,cDNA可由逆转录反应中所提取和富集的miRNA合成,且可执行qPCR反应以定量miRNA的表达水平。因此,在一些实施例中,通过qPCR,任选地使用本文中所公开的一或多个miRNA表达谱分析芯片来测定miRNA表达谱。参看实例1。
目前,qPCR分析可分成两种类型。第一种类型使用茎-环(stem-loop)逆转录引物进行cDNA合成,且使用miRNA特异性探针或通用探针来定量miRNA。第二种方法使用线性通用逆转录引物进行cDNA合成,且使用miRNA特异性正向引物、具逆转录引物特异性的反向引物和双股DNA嵌入染料来定量miRNA。
在一些实施例中,使用如美国专利第10,590,478号中所公开的通用逆转录引物来进行cDNA合成,其以引用的方式并入本文中。在一些实施例中,使用具有由以下通式表示的核苷酸序列的通用逆转录引物进行cDNA合成:5'-R-(dT)nVN-3',其中R包括CAACTCAGGTCGTAGGCAATTCGT的序列(SEQ ID NO:168),(dT)n为n个连续胸腺嘧啶残基,其中n为19,V为腺嘌呤残基、鸟嘌呤残基或胞嘧啶残基,且N为腺嘌呤残基、鸟嘌呤残基、胞嘧啶残基或胸腺嘧啶残基。
为降低成本且易于使用,在一些实施例中,可使用根据本公开的以所有167个miRNA为标靶的一或多个miRNA表达谱分析芯片来进行qPCR反应。参看实例1。在一些实施例中,每一个miRNA表达谱分析芯片预载有合适的引物和/或探针,所述引物和/或探针能够同时分析至少20个、30个、40个、50个、60个、70个、80个、90个、95个、96个、97个、98个、99个、100个、110个、120个、130个、140个、150个、160个、170个、180个、190个或200个miRNA的表达。在一些实施例中,miRNA表达谱分析芯片含有如美国专利第9,724,692号、专利第10,415,084号、申请第16/191,451号以及申请第16/233,121号中所公开的多工试片(multiplex slide plate),其以引用的方式并入本文中。
可使用所属领域中已知的方法执行qPCR反应。在一些实施例中,可使用如美国专利第9,168,533号和申请第16/559,642号中所公开的热循环仪装置来进行qPCR反应,其以引用的方式并入本文中。同样参看实例1。
miRNA分析算法和其用于测定子宫内膜容受性的用途
根据本公开的方法,miRNA表达谱可用于运用计算机miRNA分析算法产生容受性预测评分。容受性预测评分将子宫内膜状态分类为以下三种状态中的一种:容受期前状态、容受期状态或容受期后状态。
计算机miRNA分析算法为数学预测分类器,其使用miRNA表达数据且根据不同容受性状态学习区分类别。
为建构算法,将关于miRNA表达水平的原始数据分成训练组和验证组。训练组用于训练预测分类器且验证组用于评估并改善预测分类器的性能。如图3中所示,执行以下一或多个步骤来建构且验证算法:数据正规化、数据缩放、数据转换、预测建模以及交叉验证。
为使在统计特性中分布相同,数据可通过百分位正规化(QuantileNormalization)来正规化,如博尔斯塔(Bolstad)等人,于“A comparison ofnormalization methods for high density oligonucleotide array data based onvariance and bias”,(Bioinformatics,2003,19(2):185-193中所描述。此外,为确保目标函数恰当地起作用,可标准化数据的数值范围以使数据具有零均值(zero-mean)和单位变异数(unit-variance)。
出于数据简化(data reduction)和特征撷取(feature extraction)的原因,主成分分析(principal component analysis;PCA)可用于压缩来自大量原始变量的信息且通过线性地组合原始变量来产生一小组新特征。
PCA转换的数据可用于进一步建构具有弹性网正则化(elastic netregularization)的广义线性模型,其为线性地结合lasso和ridge方法的L1和L2惩罚的正则化回归法,如邹(Zou)等人,“Regularization and variable selection via theelastic net”J.R.Statist.Soc.B,2005,67,part 2,301-320中所描述。关于glmnet的额外信息为已知且可在glmnet.stanford.edu获得。
k折交叉验证方法(例如,10折交叉验证)可用于评估计算机miRNA分析算法在最终完成之前的预测值。参看图5。在k折交叉验证中,将原始样本随机地分割成k等分的子样本。在k个子样本中,其中一个子样本保留作为测试模型的验证数据,且将其余k-1个子样本用作训练数据。随后重复进行交叉验证过程k次(折),其中k个子样本中的每一个都刚好使用一次作为验证数据。来自等分数的k个结果随后可被平均化(或以其它方式结合)以产生单一估算值。
妊娠率可用于评估计算机miRNA分析算法的预测值。参看实例2。
在验证及优化之后,产出计算机miRNA分析算法。运行算法产生容受性预测评分,所述评分将女性的子宫内膜状态分类成如下三种状态中的一个:若评分大于1,则女性的子宫内膜处于容受期前状态;若评分小于-1,则女性的子宫内膜处于容受期后状态;若评分介于-1到1,则女性的子宫内膜处于容受期状态。参见图6。
根据本公开的方法的应用
本公开提供一种使用样本,例如子宫内膜检体以测定子宫内膜状态的方法,所述方法包括:(a)对女性的子宫内膜样本进行分析以测定子宫内膜样本的miRNA表达谱,其中miRNA表达谱包括多个miRNA,例如分别具有SEQ ID NO:1到SEQ ID NO:167的序列的167个miRNA的表达水平;和(b)使用例如计算机算法来分析miRNA表达谱以获得容受性预测评分。
本公开的方法可用于各种诊断和治疗目的,包含(但不限于)IVF疗程。举例来说,在一些实施例中,基于子宫内膜的结果,可进一步包含在女性体内植入胚胎或向遭受或曾遭受胚胎植入失败的女性投予一或多种疗程的方法。在一些实施例中,本公开提供检测胚胎植入的子宫内膜容受性的方法,其包括:(a)对女性的子宫内膜样本进行分析以测定子宫内膜样本的miRNA表达谱,其中miRNA表达谱包括多个miRNA例如分别具有SEQ ID NO:1到SEQ ID NO:167的序列的167个miRNA的表达水平,(b)分析miRNA表达谱以获得容受性预测评分,其中容受性预测评分判定女性是否具有子宫内膜容受性,以及(c)将胚胎转移至被测定具有子宫内膜容受性的女性的子宫内膜。
在一些实施例中,测定子宫内膜状态的方法可用于判定女性胚胎植入的时机。在一些实施例中,若子宫内膜状态处于容受期状态,则认为女性适合胚胎植入。若子宫内膜状态处于容受期前或容受期后状态,则认为女性不适合胚胎植入。在一些实施例中,当子宫内膜状态被判定为处于容受期前状态或容受期后状态时,本公开提供基于子宫内膜状态的信息的胚胎植入的方法。举例来说,若子宫内膜状态被判定为处于容受期前状态,则在下一周期期间,可在孕酮投药之后的5.5天到7.5天之间,例如5.5天、6天、6.5天、7天或7.5天进行胚胎植入。或者,若子宫内膜状态被判定为处于容受期后状态,则在下一周期期间,可在孕酮投药之后的2.5天到4.5天之间,例如2.5天、3天、3.5天、4天或4.5天进行胚胎植入。
在取样时子宫内膜显示非容受期状态的情况下,所获得信息为具指导性的,使得所述方法可通过在另一时间点获取子宫内膜样本来重复进行,以根据首次判定结果进行修改。借助于实例,若子宫内膜状态处于容受期前状态,则获取子宫内膜样本的下一个时间点可在内源性LH遽增之后超过七天或在孕酮投药之后超过五天。举例来说,获取子宫内膜样本的下一个时点可在内源性LH遽增之后7.5天到10.5天之间,例如7.5天、8天、8.5天、9天、9.5天、10天、或10.5天,或在孕酮投药之后5.5天到7.5天之间,例如5.5天、6天、6.5天、7天或7.5天。或者,若子宫内膜状态处于容受期后状态,则获取子宫内膜样本的下一个时间点可在内源性LH遽增之后少于七天或在孕酮投药之后少于五天。举例来说,获取子宫内膜样本的下一个时点可在内源性LH遽增之后3.5天到6.5天之间,例如3.5天、4天、4.5天、5天、5.5天、6天或6.5天,或在孕酮投药之后2.5天到4.5天之间,例如2.5天、3天、3.5天、4天或4.5天。通过遵循这些程序,可得出容受期状态,且可提高IVF疗程的成功率。对于这些用途中的任何一种,女性遭受或曾遭受胚胎植入失败。在一些实施例中,女性进行IVF疗程。
在一些实施例中,若子宫内膜状态被判定为处于容受期前状态或容受期后状态,则可将测定子宫内膜状态的方法重复至少一次或直到子宫内膜状态被判定为处于容受期状态。
在一些实施例中,根据本公开的测定子宫内膜状态的方法可用于判定女性的WOI。在一些实施例中,根据本公开的方法可用于分类女性对IVF治疗的反应性。对于这些用途中的任何一种,在一些实施例中,女性遭受或曾遭受胚胎植入失败。在一些实施例中,女性进行IVF疗程。
在一些实施例中,根据本公开的测定子宫内膜状态的方法可作为探讨妊娠药物对女性子宫内膜影响的有用工具。在这些实施例中,女性遭受或曾遭受胚胎植入失败。在一些实施例中,使女性进行IVF疗程。
试剂盒
本公开的另一方面是有关于用于实施测定子宫内膜状态方法的试剂盒。在一些实施例中,试剂盒包括适合于检测多个miRNA,例如分别具有SEQ ID NO:1到SEQ ID NO:167的序列的167个miRNA的表达水平的引物和/或探针。参看实例1。在一些实施例中,引物和/或探针适合于进行qPCR反应以检测167个miRNA的表达水平。在一些实施例中,试剂盒包括一或多个以167个miRNA为标靶的miRNA表达谱分析芯片。在一些实施例中,一或多个芯片另外以可作为miRNA表达分析的内源性对照的RNA序列为标靶,例如18s rRNA。
试剂盒可另外含有关于(i)任选地使用一或多个miRNA表达谱分析芯片测定来自女性子宫内膜样本的miRNA表达谱,和/或(ii)基于miRNA表达谱使用计算机算法获得容受性预测评分的使用说明。在一些实施例中,试剂盒含有关于如何解释且运用容受性预测评分的指示说明。
在一些实施例中,试剂盒有助于诊断和治疗目的,包含但不限于IVF疗程。
实例
实例1:产出miRNA表达谱的材料和方法。
子宫内膜检体。在激素替代治疗周期中孕酮给药之后五天(P+5)或在自然周期中内源性促黄体激素遽增之后七天(LH+7),使用Pipelle子宫内膜吸引刮匙(库柏外科公司(Cooper Surgical,Inc.))从女性子宫腔采集子宫内膜检体。紧接着将子宫内膜组织储存于RNAlater中。
RNA提取和miRNA富集。遵循制造商的说明书使用miRNeasy Micro Kit(凯杰(QIAGEN))从子宫内膜组织分离总RNA。简要来说,将五毫克子宫内膜组织用电动机和研杵破碎且在液氮中均质化。将700微升QIAzol Lysis Reagent添加到均质化组织中,且将所得样本在室温下培育五分钟以促进核蛋白复合物的分解。每700微升QIAzol Lysis Reagent添加140微升三氯甲烷到试管中,且手动剧烈摇晃试管15秒并在室温下培育2分钟到3分钟。将样本在4℃下以12,000克离心15分钟。在离心之后,将上部水相转移到新试管,将一体积的70%乙醇添加到试管中,且充分地涡动试管。将样本转移到RNeasy MinElute旋转柱且在室温下以8,000克离心15秒。将流过物移液到2毫升试管内,将0.65体积的100%乙醇添加到流过物中,且充分地涡动所得样本。随后将样本转移到RNeasy MinElute旋转柱内且在室温下以8,000克离心15秒。丢弃流过物,将700微升缓冲液RWT添加到RNeasy MinElute旋转柱中,且将柱以8000克离心15秒以洗涤柱。丢弃流过物,将500微升缓冲液RPE添加到RNeasyMinElute旋转柱内,且将柱以8,000克离心15秒以洗涤柱。丢弃流过物,将500微升的80%乙醇添加到RNeasy MinElute旋转柱内,且将柱以8,000克离心2分钟以干燥旋转柱膜片。将RNeasy MinElute旋转柱放置于新的2毫升收集试管内且以8,000克离心5分钟。将RNeasyMinElute旋转柱放置于1.5毫升收集试管内,将14微升到20微升的无核酸酶水添加到旋转柱膜片上,且将柱以8,000克离心1分钟以洗脱富集miRNA的级分。在-80℃下储存富集miRNA的级分。
cDNA合成。在20微升逆转录反应中,来自子宫内膜组织的≥2纳克富集miRNA的级分用于合成cDNA。遵循制造商的说明书使用QuarkBio microRNA Universal RT Kit(奎克生技光电股份有限公司(Quark Biosciences Taiwan,Inc.))进行逆转录。简要来说,使用poly-A聚合酶将poly-A尾部添加到miRNA,之后进行cDNA合成。随后使用以下程序执行cDNA合成:42℃持续60分钟和95℃持续5分钟,且然后4℃直到程序完成。在-20℃下储存合成的cDNA。
使用NextAmp分析系统和MIRA PanelChip组进行miRNA表达谱分析。MIRAPanelChip组含有总共167个miRNA分析。167个miRNA的序列列示于表1中。另外,RNU6B、RNU43以及18s rRNA均作为内源性对照。三个外源性插入对照用于监测miRNA提取、cDNA合成以及qPCR效能(奎克生技光电股份有限公司)。使用MIRA PanelChip组分析cDNA。将cDNA(等同于0.1纳克富集miRNA的级分)添加到含有30微升2×SYBR Master Mix的混合物中(奎克生技光电股份有限公司),且将无核酸酶的水添加到混合物中以获得60微升的最终体积。将混合物手动充分地混合且短暂地快速离心以收集底部的液体。使用Pipetman将60微升混合物沿芯片边缘分配且随后通过用玻璃载片进行刮擦动作将混合物施加在MIRAPanelChip的整个表面上。随后将每个芯片浸没到含有Channeling Solution(奎克生技光电股份有限公司)的托盘内,其中使反应孔面向托盘的底部。然后将每个托盘放入到QStation,其为MIRA PanelChip应用的热循环仪(参看图2中的PanelStation)且包含内建的样本管理数据库和分析平台,使得MIRA PanelChip检测和数据分析可便利且快速地进行。其后根据以下程序执行MIRA PanelChip分析:95℃持续36秒和60℃持续72秒,持续40个周期。
实例2:计算机miRNA分析算法和其用途。
如图3中所示,计算机miRNA分析算法(MIRA)通过执行以下一或多个步骤来建构:数据正规化、数据缩放、数据转换、预测建模以及交叉验证。
数据正规化。为使在统计特性中分布相同,通过百分位正规化来正规化数据。参看图3中的算式(A);也参看博尔斯塔(Bolstad)等人,“A comparison of normalizationmethods for high density oligonucleotide array data based on variance andbias”,Bioinformatics,2003,19(2):185-193。
数据缩放。为确保目标函数恰当地作用,可标准化数据的数值范围以使数据具有零均值和单位变异数。参看图3中的算式(B)。
数据转换。出于数据简化和特征撷取的原因,PCA压缩来自大量原始变量的信息且通过线性地组合原始变量来产生一小组新特征。请参照图3中的算式(C)。
建模。PCA转换的数据用于进一步建构具有弹性网正则化的广义线性模型,其为线性地结合lasso和ridge方法的L1和L2惩罚的正则化回归法。参看图3中的算式(D);也参看邹等人(Zou),“Regularization and variable selection via the elastic net”,J.R.Statist.Soc.B,2005,67,part 2,301-320。
在完成MIRA模型之前进行交叉验证以评估计算机miRNA分析算法的预测值。如图4A所示,使用含有具有表1中所示的SEQ ID NO:1到SEQ ID NO:167的序列的167个miRNA的表达水平的miRNA表达谱,MIRA模型能够成功地将临床样本分类成以下三种状态组中的一种:容受期前状态、容受期状态、容受期后状态。此外,如图4B所示,初步验证显示分类为容受期状态的女性(测试组)的妊娠率为100%。
将来自183位女性的数据分成10个子集以达到模型评估的10折交叉验证。图5示出使用183个子宫内膜样本所得的具167个miRNA表达水平的miRNA表达谱之10折交叉验证和妊娠率。在这些测试中,于第一周期中,测定每位女性的子宫内膜状态。若女性的子宫内膜被判定为处于容受期前状态,则在下一周期孕酮投药之后六天进行胚胎植入(P+6组;35位女性)。若女性的子宫内膜被判定为处于容受期状态,则在下一周期孕酮投药之后五天进行胚胎植入(P+5组;142位女性)。若女性的子宫内膜被判定为处于容受期后状态,则在下一周期孕酮投药之后4.5天进行胚胎植入(P+4.5组;6位女性)。另外,图5示出灵敏性、特异性、PPV、NPV和10折交叉验证结果的总体一致率。
在三个组当中,检测到137件妊娠事件,其中22件事件来自P+6组,113件事件来自P+5组以及2件事件来自P+4.5组。参看图5,关于计算机miRNA分析算法的预测评估,在所有137件妊娠事件当中,2分之1来自P+4.5组、113分之107来自P+5组以及22分之17来自P+6组,显示出可通过算法判定正确的胚胎植入时间调整且产生91.24%的妊娠率(125/137)。参看图5。
MIRA模型。考虑描述于此实例中的所有参数(参看图3,eq(A-D)且随后根据交叉验证微调其参数),产出将所有样本分类成三种不同子宫内膜状态的预测模型。运行MIRA产生容受性预测评分(MIRA评分),其使用以下算式计算:MIRA score=f(X∈eq(C))=Xβ+ε,其中β为系数向量,且ε为误差,均通过交叉验证由glmnet产出(图3)。这个模型可适用于子宫内膜的任何qPCR图谱分析以预测子宫内膜状态。
如图6所示,运行计算机miRNA分析算法产出容受性预测评分,其将女性的子宫内膜状态分类成如下三种状态中的一种:若评分大于1,则女性的子宫内膜处于容受期前状态;若评分小于-1,则女性的子宫内膜处于容受期后状态;若评分介于-1到1,则女性的子宫内膜处于容受期状态(WOI)。
尽管本公开已参照具体实施方式进行特定呈现和描述,所属领域的技术人员应理解在不脱离本公开的精神和范围的情况下,可对各种形式和细节做出改变。
参考文献
安布罗斯(Ambros),“微小核糖核酸:具有大潜能的微小调节子(microRNAs:tinyregulators with great potential)”,《细胞(Cell)》,2001,107(7):823-826。
巴特尔(Bartel),“微小核糖核酸:基因组学、生物合成、机制以及功能(MicroRNAs:genomics,biogenesis,mechanis,and function)”,《细胞(Cell)》,2004,116(2):81-97。
李(Lee)等人,“秀丽隐杆线虫异时基因lin-4编码与lin-14具有反义互补性的小RNA(The C.elegans heterochronic gene lin-4encodes small RNAs with antisensecomplementarity to lin-14)”,《细胞(Cell)》,1993,75(5):843-854。
帕斯奎内利(Pasquinelli)等人,“let-7异时调节RNA的序列保守性和瞬时表达(Conservation of the sequence and temporal expression of let-7heterochronicregulatory RNA)”,《自然(Nature)》,2000,408(6808):86-89。
怀特曼(Wightman)等人,“lin-4对异时基因lin-14的转录后调节介导秀丽隐杆线虫的暂态模式形成(Posttranscriptional regulation of the heterochronic genelin-14by lin-4mediates temporal pattern formation in C.elegans)”,《细胞(Cell)》,1993,75(5):855-862。
纳沃特(Navot)等人,“胚胎移植的窗口期和人类体外受孕的功效(The window ofembryo transfer and the efficiency of human conception in vitro)”,《生育与不孕(Fertil Steril)》,1991,55(1):114-118。
哈珀(Harper),“植入窗口期(The implantation window)”,《贝利勒临床妇产科(Baillieres Clin Obstet Gynaecol)》,1992,6(2):351-371。
赫尔(Hull)等人,“组织和循环微小核糖核酸影响子宫内膜疾病的生殖功能(Tissue and circulating microRNA influence reproductive function inendometrial disease)”,《生殖生物医学在线(Reprod Biomed Online)》,2013,27(5):515-529。
加利亚诺(Galliano)等人,“微小核糖核酸和植入(MicroRNA andimplantation)”,《生育与不孕(Fertil Steril)》,2014,101(6):1531-1544。
刘(Liu)等人,“微小核糖核酸和胚胎植入(MicroRNA and EmbryoImplantation)”,《美国生殖免疫学杂志(Am J Reprod Immunol)》,2016,75(3):263-271。
史(Shi)等人,“植入窗口期间的子宫内膜微小核糖核酸标签在反复植入失败的患者中发生改变(Endometrial MicroRNA Signature during the Window of ImplantationChanged in Patients with Implantation Failure)”,《新英格兰医学杂志(Chin Med J(Engl))》,2017,130(5):566-573。
克瑞斯威克(Kresowik)等人,“微小核糖核酸-31在植入窗口期间在人类子宫内膜和血清中显著地升高:最优容受性的潜在生物标记物(MicroRNA-31is significantlyelevated in both human endometrium and serum during the window ofimplantation:a potential biomarker for optimum receptivity)”,《生殖生物学(BiolReprod)》,2014,91(1):17。
李(Li)等人,“微小核糖核酸-30a-3p调节上皮间叶细胞跃迁以通过靶向Snai2dagger来影响胚胎植入(MicroRNA-30a-3p regulates epithelial-mesenchymaltransition to affect embryo implantation by targeting Snai2dagger)”,《生殖生物学(Biol Reprod)》,2019,100(5):1171-1179。
梁(Liang)等人,“微小核糖核酸在胚胎植入中的作用(Role of microRNAs inembryo implantation)”,《生殖生物学和内分泌学(Reprod Biol Endocrinol)》,2017,15(1):90。
保罗(Paul)等人,“微小核糖核酸在人类胚胎植入中的作用:综述(The role ofmicroRNAs in human embryo implantation:a review)”,《辅助生殖技术和遗传学杂志(JAssist Reprod Genet)》,2019,36(2):179-187。
维莱纳(Vilella)等人,“通过人类子宫内膜分泌的Hsa-miR-30d由植入前胚胎吸收且可能修改其转录物组(Hsa-miR-30d,secreted by the human endometrium,is takenup by the pre-implantation embryo and might modify its transcriptome)”,《发育(Development)》,2015,142(18):3210-3221。
郑(Zheng)等人,“微小核糖核酸-200c通过靶向FUT4和α1,3-岩藻糖基化损害子宫容受性形成(MicroRNA-200c impairs uterine receptivity formation by targetingFUT4 andα1,3-fucosylation)”,《细胞死亡与分化(Cell Death Differ)》,2017,24(12):2161-2172。
博尔斯塔(Bolstad)等人,“基于房差和偏差的高密度寡核苷酸阵列归一化方法的比较(A comparison of normalization methods for high density oligonucleotidearray data based on variance and bias)”,《生物信息(Bioinformatics)》,2003,19(2):185-193。
邹(Zou)等人,“通过弹性网进行正则化和变量选择(Regularization andvariable selection via the elastic net)”,《皇家统计学会杂志系列B(J.R.Statist.Soc.B)》,2005,67,第2部分,301-320。
施里姆儿(Schriml)等人,“人类疾病本体2018更新:分类、内容以及工作流程扩展(Human Disease Ontology 2018update:classification,content and workflowexpansion)”,《核酸研究(Nucleic Acids Res)》,2019,47:D955-D62。
黄(Huang)等人,“miRTarBase 2020:实验上地验证的微小核糖核酸-靶向交互数据库的更新(miRTarBase 2020:updates to the experimentally validated microRNA-target interaction database)”,《核酸研究(Nucleic Acids Res)》,2020,48:D148-D54。
阿加瓦尔(Agarwal)等人,“预测哺乳动物mRNA中的有效微小核糖核酸靶向部位(Predicting effective microRNA target sites in mammalian mRNAs)”,2015,《Elife》,4。
陈(Chen)等人,“miRDB:用于预测功能性微小核糖核酸目标的在线数据库(miRDB:an online database for prediction of functional microRNA targets)”,《核酸研究(Nucleic Acids Res)》,48:D127-D31。
王(Wang)等人,“关于来自两种商业长寡核苷酸微阵列的基因表达测量的大规模即时PCR验证(Large scale real-time PCR validation on gene expressionmeasurements from two commercial long-oligonucleotide microarrays)”,《BMC基因组学(BMC Genomics)》,2006,7:59-75。

Claims (30)

1.一种测定子宫内膜状态的方法,包括:
(a)对来自女性的子宫内膜样本进行分析,以测定所述子宫内膜样本的miRNA表达谱,其中所述miRNA表达谱包括多个miRNA的表达水平;以及
(b)分析所述miRNA表达谱以获得容受性预测评分,其中所述容受性预测评分对所述女性的子宫内膜状态进行分类,所述子宫内膜状态包括容受期前状态、容受期状态或容受期后状态,且所述多个miRNA包括至少50个、75个、100个、125个、150个或200个miRNA,且优先选择分别具有SEQ ID NO:1到SEQ ID NO:167的序列的至少167个miRNA。
2.根据权利要求1所述的测定子宫内膜状态的方法,其中所述子宫内膜样本取自所述女性的子宫腔。
3.根据权利要求1或权利要求2所述的测定子宫内膜状态的方法,其中所述子宫内膜样本包括子宫内膜检体、子宫内膜灌洗液或其组合。
4.根据权利要求1-3中任一项所述的测定子宫内膜状态的方法,其中所述子宫内膜样本(i)在所述女性的内源性促黄体激素(luteinizing hormone;LH)遽增之后七天或(ii)在所述女性的孕酮投药之后五天获取。
5.根据权利要求1-4中任一项所述的测定子宫内膜状态的方法,其中所述miRNA表达谱通过qPCR、定序、微阵列芯片或RNA-DNA杂交捕获技术来测定。
6.根据权利要求5所述的测定子宫内膜状态的方法,其中所述miRNA表达谱通过对由所述子宫内膜样本中的miRNA所合成的cDNA进行qPCR来测定。
7.根据权利要求6所述的测定子宫内膜状态的方法,其中cDNA合成使用具有由以下通式表示的核苷酸序列的通用逆转录引物进行:5'-R-(dT)nVN-3',其中R包括SEQ ID NO:168,(dT)n为n个连续胸腺嘧啶残基,其中n为19,V为腺嘌呤残基、鸟嘌呤残基或胞嘧啶残基,且N为腺嘌呤残基、鸟嘌呤残基、胞嘧啶残基或胸腺嘧啶残基。
8.根据权利要求1-7中任一项所述的测定子宫内膜状态的方法,其中所述容受性预测评分为通过计算机算法产出且使用算式MIRA score=f(X∈eq(C))=Xβ+ε所计算的数值,其中β为系数向量,且ε为误差。
9.根据权利要求8所述的测定子宫内膜状态的方法,其中计算机算法通过进行以下一或多个步骤来建立:数据正规化、数据缩放、数据转换、预测建模以及交叉验证。
10.根据权利要求8或权利要求9所述的测定子宫内膜状态的方法,其中所述容受性预测评分大于1指容受期前状态,所述容受性预测评分小于-1指容受期后状态,所述容受性预测评分介于-1到1指容受期状态。
11.根据权利要求1-10中任一项所述的测定子宫内膜状态的方法,其中若所述子宫内膜状态被判定为处于容受期前状态或容受期后状态,则所述测定子宫内膜状态的方法另外包括:重复步骤(a)和(b)至少一次或直到所述子宫内膜状态被判定为处于容受期状态。
12.根据权利要求1-11中任一项所述的测定子宫内膜状态的方法,其中所述女性遭受或曾遭受胚胎植入失败。
13.根据权利要求1-12中任一项所述的测定子宫内膜状态的方法,其中使所述女性进行体外受精疗程。
14.根据权利要求13所述的测定子宫内膜状态的方法,其中所述容受性预测评分进一步分类所述女性对体外授精疗程的反应性。
15.一种检测女性胚胎植入的子宫内膜容受性的方法,包括:
(a)对来自所述女性的子宫内膜样本进行分析,以测定所述子宫内膜样本的miRNA表达谱,其中所述miRNA表达谱包括多个miRNA的表达水平;以及
(b)分析所述miRNA表达谱以获得容受性预测评分,其中所述容受性预测评分判定所述女性是否具有胚胎植入的子宫内膜容受性,且其中所述多个miRNA包括至少50个、75个、100个、125个、150个或200个miRNA,且优先选择分别具有SEQ ID NO:1到SEQ ID NO:167的序列的至少167个miRNA。
16.根据权利要求15所述的检测女性胚胎植入的子宫内膜容受性的方法,其中所述子宫内膜样本取自所述女性的子宫腔。
17.根据权利要求15或权利要求16所述的检测女性胚胎植入的子宫内膜容受性的方法,其中所述子宫内膜样本包括子宫内膜检体、子宫内膜灌洗液或其组合。
18.根据权利要求15-17中任一项所述的检测女性胚胎植入的子宫内膜容受性的方法,其中所述子宫内膜样本(i)在所述女性内源性促黄体激素(LH)遽增之后七天或(ii)在所述女性孕酮投药之后五天获取。
19.根据权利要求15-18中任一项所述的检测女性胚胎植入的子宫内膜容受性的方法,其中所述miRNA表达谱通过qPCR、定序、微阵列芯片或RNA-DNA杂交捕获技术来测定。
20.根据权利要求19所述的检测女性胚胎植入的子宫内膜容受性的方法,其中所述miRNA表达谱通过对由所述子宫内膜样本中的miRNA所合成的cDNA进行qPCR来测定。
21.根据权利要求20所述的检测女性胚胎植入的子宫内膜容受性的方法,其中cDNA合成使用具有由以下通式表示的核苷酸序列的通用逆转录引物进行:5'-R-(dT)nVN-3',其中R包括SEQ ID NO:168,(dT)n为n个连续胸腺嘧啶残基,n为19,V为腺嘌呤残基、鸟嘌呤残基或胞嘧啶残基且N为腺嘌呤残基、鸟嘌呤残基、胞嘧啶残基或胸腺嘧啶残基。
22.根据权利要求15-21中任一项所述的检测女性胚胎植入的子宫内膜容受性的方法,其中所述容受性预测评分为通过计算机算法产出且使用算式MIRA score=f(X∈eq(C))=Xβ+ε所计算的数值,其中β为系数向量,且ε为误差。
23.根据权利要求22所述的检测女性胚胎植入的子宫内膜容受性的方法,其中计算机算法通过进行以下一或多个步骤来建立:数据正规化、数据缩放、数据转换、预测建模以及交叉验证。
24.根据权利要求22或权利要求23所述的检测女性胚胎植入的子宫内膜容受性的方法,其中所述容受性预测评分介于-1到1指所述女性具有胚胎植入的子宫内膜容受性。
25.根据权利要求15-24中任一项所述的检测女性胚胎植入的子宫内膜容受性的方法,其中所述女性遭受或曾遭受胚胎植入失败。
26.一种试剂盒,包括:
(a)一或多个以多个miRNA为标靶的miRNA表达谱分析芯片;以及
(b)关于(i)任选地使用一或多个所述miRNA表达谱分析芯片测定女性子宫内膜样本的miRNA表达谱及(ii)基于所述miRNA表达谱使用计算机算法获取容受性预测评分的使用说明,其中所述多个miRNA包括至少50个、75个、100个、125个、150个或200个,且优先选择分别具有SEQ ID NO:1到SEQ ID NO:167的序列的至少167个miRNA。
27.根据权利要求26所述的试剂盒,其中一或多个所述miRNA表达谱分析芯片包括用于检测所述多个miRNA的表达水平的引物。
28.根据权利要求27所述的试剂盒,其中所述miRNA表达谱分析芯片适合于进行qPCR、定序、微阵列芯片或RNA-DNA杂交捕获分析,优先选择qPCR,以检测所述多个miRNA的表达水平。
29.一种试剂盒的用途,所述试剂盒为如权利要求27或权利要求28所述的试剂盒,用于测定女性的子宫内膜状态。
30.根据权利要求29所述的试剂盒的用途,其中所述女性遭受或曾遭受胚胎植入失败和/或进行体外受精疗程。
CN202080004080.1A 2019-07-02 2020-07-01 测定子宫内膜状态的方法、检测女性胚胎植入的子宫内膜容受性的方法、试剂盒及用途 Active CN112469836B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201962869574P 2019-07-02 2019-07-02
US62/869,574 2019-07-02
PCT/CN2020/099781 WO2021000893A1 (en) 2019-07-02 2020-07-01 Mirna receptivity analysis of the endometrium

Publications (2)

Publication Number Publication Date
CN112469836A true CN112469836A (zh) 2021-03-09
CN112469836B CN112469836B (zh) 2024-03-22

Family

ID=74066350

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202080004080.1A Active CN112469836B (zh) 2019-07-02 2020-07-01 测定子宫内膜状态的方法、检测女性胚胎植入的子宫内膜容受性的方法、试剂盒及用途

Country Status (6)

Country Link
US (1) US20210002698A1 (zh)
EP (1) EP3999657A4 (zh)
JP (2) JP2022539037A (zh)
CN (1) CN112469836B (zh)
TW (1) TWI753474B (zh)
WO (1) WO2021000893A1 (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20210007275A (ko) 2019-07-10 2021-01-20 삼성전자주식회사 오버레이 보정 방법, 및 그 보정 방법을 기초로 한 포토리소그라피 방법, 반도체 소자 제조방법 및 스캐너 시스템
CN112662758B (zh) * 2021-02-07 2021-08-06 成都西囡妇科医院有限公司 一种与子宫内膜容受性辅助诊断相关的miRNA标志物及其应用
EP4183888A1 (en) * 2021-11-22 2023-05-24 Asociación Centro de Investigación Cooperativa en Biociencias - CIC bioGUNE Mirna signature for identification of the receptive endometrium
WO2023202612A1 (en) * 2022-04-19 2023-10-26 The Chinese University Of Hong Kong The identification of endometrial immune cell densities and clustering analysis in the mid-luteal phase as predictor for pregnancy outcomes

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014062442A1 (en) * 2012-10-15 2014-04-24 University Of Iowa Research Foundation Use of microrna for assessing fertility in a female patient
US20170022544A1 (en) * 2015-06-02 2017-01-26 Quark Biosciences, Inc. NUCLEOTIDE SEQUENCE, UNIVERSAL REVERSE PRIMER, UNIVERSAL RT PRIMER, METHOD FOR DESIGNING PRIMER AND miRNA DETECTION METHOD
WO2018178171A1 (en) * 2017-03-29 2018-10-04 INSERM (Institut National de la Santé et de la Recherche Médicale) Methods for assessing pregnancy outcome

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2009281969A1 (en) * 2008-08-12 2010-02-18 The Ohio State University Research Foundation Micro-RNA-based compositions and methods for the diagnosis, prognosis and treatment of multiple myeloma

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014062442A1 (en) * 2012-10-15 2014-04-24 University Of Iowa Research Foundation Use of microrna for assessing fertility in a female patient
US20170022544A1 (en) * 2015-06-02 2017-01-26 Quark Biosciences, Inc. NUCLEOTIDE SEQUENCE, UNIVERSAL REVERSE PRIMER, UNIVERSAL RT PRIMER, METHOD FOR DESIGNING PRIMER AND miRNA DETECTION METHOD
WO2018178171A1 (en) * 2017-03-29 2018-10-04 INSERM (Institut National de la Santé et de la Recherche Médicale) Methods for assessing pregnancy outcome

Also Published As

Publication number Publication date
WO2021000893A1 (en) 2021-01-07
JP2022539037A (ja) 2022-09-07
US20210002698A1 (en) 2021-01-07
CN112469836B (zh) 2024-03-22
JP2023139115A (ja) 2023-10-03
EP3999657A1 (en) 2022-05-25
TW202108770A (zh) 2021-03-01
TWI753474B (zh) 2022-01-21
EP3999657A4 (en) 2023-07-05

Similar Documents

Publication Publication Date Title
CN112469836B (zh) 测定子宫内膜状态的方法、检测女性胚胎植入的子宫内膜容受性的方法、试剂盒及用途
US11111541B2 (en) Diagnostic MiRNA markers for Parkinson's disease
GB2597379A (en) Methods and systems for determining a pregnancy-related state of a subject
CN107660234A (zh) 使用二代测序预测器官移植排斥的方法
US20180327856A1 (en) Diagnostic mirnas for differential diagnosis of incidental pancreatic cystic lesions
CN109504784B (zh) 用于预测人辅助生殖技术中早期胚胎质量的miRNA分子标志及其应用
WO2010069129A1 (zh) 非小细胞肺癌标记物及其检测方法、试剂盒和生物芯片
JP2015527870A (ja) 推定出生児が状態を発症する危険性を評価するための方法およびデバイス
CN111778326B (zh) 用于子宫内膜容受性评估的基因标志物组合及其应用
JPWO2020168118A5 (zh)
Chen et al. A novel platform for discovery of differentially expressed microRNAs in patients with repeated implantation failure
KR102089559B1 (ko) 난소의 fsh에 대한 반응성의 진단 또는 예측용 바이오마커 및 이의 용도
GB2614979A (en) Methods and systems for determining a pregnancy-related state of a subject
KR20200002241A (ko) 비만 진단을 위한 마이크로RNA-26b 또는 마이크로RNA-4449 바이오마커 및 이의 용도
CN113755570B (zh) 用于预测不明原因复发性流产的生物标志物及应用
CN113755571B (zh) 用于胚胎着床成功率检测的生物标志物及应用
CN113201592A (zh) 一种鉴定肿瘤类器官的方法
CN111944894A (zh) 用于唇腭裂、神经管畸形和先天性心脏病胎儿产前无创诊断的分子标志物及其应用
WO2019168971A1 (en) Methods for assessing risk of increased time-to-first-conception
EP4311862A1 (en) Methods for detection of embryo implantation failure of endometrial origen
Van Hoeck et al. Gene expression profiling by high throughput sequencing to determine signatures for the bovine receptive uterus at early gestation
KR102480430B1 (ko) 주사증 진단을 위한 마이크로RNA-31-5p 및 이의 용도
EP4183888A1 (en) Mirna signature for identification of the receptive endometrium
CN118006758A (zh) 一种用于2型糖尿病性骨质疏松症诊断的生物标志物及其应用
KR20220052474A (ko) 주사증 진단을 위한 마이크로RNA-21-3p 및 이의 용도

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20220809

Address after: 10 / F, 7 / F, No. 81, Linshui Road, 10 qianjiali, East District, Hsinchu City, Taiwan, China, China

Applicant after: Inti Biotechnology Co., Ltd.

Address before: Taiwan China jhubei City, Hsinchu Science Park Hsinchu County medical Road 2 No. 6-2 4 floor 1

Applicant before: QUARK BIOSCIENCES TAIWAN, Inc.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant