[go: up one dir, main page]

WO2010010201A1 - Perfil de expresion genetica como marcador de la receptividad endometrial - Google Patents

Perfil de expresion genetica como marcador de la receptividad endometrial Download PDF

Info

Publication number
WO2010010201A1
WO2010010201A1 PCT/ES2008/000513 ES2008000513W WO2010010201A1 WO 2010010201 A1 WO2010010201 A1 WO 2010010201A1 ES 2008000513 W ES2008000513 W ES 2008000513W WO 2010010201 A1 WO2010010201 A1 WO 2010010201A1
Authority
WO
WIPO (PCT)
Prior art keywords
endometrium
profile
receptive
gene
abnormality
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.)
Ceased
Application number
PCT/ES2008/000513
Other languages
English (en)
French (fr)
Inventor
Carlos Simon Valles
José Antonio HORCAJADAS ALMANSA
Patricia Diaz Gimeno
Antonio Pellicer Martinez
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.)
Equipo IVI Investigacion SL
Original Assignee
Equipo IVI Investigacion SL
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 Equipo IVI Investigacion SL filed Critical Equipo IVI Investigacion SL
Priority to PCT/ES2008/000513 priority Critical patent/WO2010010201A1/es
Priority to EP09800092.0A priority patent/EP2333107B1/en
Priority to US13/057,135 priority patent/US10081840B2/en
Priority to PL09800092T priority patent/PL2333107T3/pl
Priority to ES09800092.0T priority patent/ES2484417T3/es
Priority to PCT/ES2009/000386 priority patent/WO2010010213A1/es
Priority to DK09800092.0T priority patent/DK2333107T3/da
Priority to CA2732849A priority patent/CA2732849C/en
Publication of WO2010010201A1 publication Critical patent/WO2010010201A1/es
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING 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/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6883—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING 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/00—Oligonucleotides characterized by their use
    • C12Q2600/112—Disease subtyping, staging or classification
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING 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/00—Oligonucleotides characterized by their use
    • C12Q2600/158—Expression markers

Definitions

  • the present invention relates to the determination of the receptivity of the human endometrium from a gene expression profile. More specifically, it consists in the elaboration of a specific expression microarray of endometrial receptivity (ERA or Endometrial Receptivity Array) that allows evaluating the receptive state of a human endometrium, as well as assessing the state of the endometrium for diagnostic and therapeutic purposes.
  • ERA Endometrial Receptivity Array
  • the endometrium is the mucosa that lines the inside of the uterine cavity. Its function is to house the embryo, allowing its implantation and favoring the development of the placenta. This process requires a receptive endometrium capable of responding to blast signals that is the stage of development in which the embryo is when it implants.
  • the human endometrium is a cyclic hormonally regulated tissue, the hormones that prepare it to reach this state of receptivity are estradiol that induces cell proliferation and progesterone that is involved in differentiation, causing a large number of changes in the profile of gene expression of the endometrium that reaches a receptive phenotype for a short period of time called the " implantation window " .
  • Noyes et al. Described, for the first time, an endometrial dating method based exclusively on histological criteria and morphological changes in the different compartments of the endometrium in response to the presence of estrogens and progesterone.
  • Noyes and his collaborators examined the histological features of endometrial biopsies taken during 8,000 spontaneous cycles in 300 women (Noyes et al., 1950; 1953). They were able to relate different histological patterns with particular moments of the menstrual cycle by correlating the histological changes with the basal body temperature. These morphological criteria are still used today and are considered as the "GoId Standard 'for the study of the endometrium, evaluation of endometrial receptivity and detection of endometrial abnormalities.
  • ASRM American Society of Reproductive Medicine
  • a DNA microarray consists of a large number of DNA molecules arranged on a solid substrate so that they form a matrix of sequences in two or three dimensions. These fragments of genetic material can be short sequences called oligonucleotides or larger, such as complementary DNA (cDNA) that is synthesized from mRNA, or PCR products (in vitro replication of DNA sequences through chain reaction of polymerase).
  • cDNA complementary DNA
  • probes These single-stranded nucleotide fragments immobilized on the support are called "probes.”
  • the nucleic acids of the samples to be analyzed are marked by various methods (enzymatic, fluorescent, etc.) and incubated on the panel of probes, allowing hybridization (recognition and binding between complementary molecules) of homologous sequences.
  • hybridization the samples of labeled genetic material bind to their complementary immobilized in the support of the chip, allowing the identification and quantification of the DNA present in the sample.
  • the appropriate bioinformatics scanner and tools allow interpreting and analyzing the data obtained (Al-Shahrour F et al., 2005).
  • microarray you can use those that exist commercially, or you can design a custom one.
  • Support material glass, plastic, membranes ...
  • microarrays also called bioarrays or biochips
  • manufacture of microarrays has been described in several patent documents, such as, for example, WO 2005/018796 A1, US 2005/0048554 A1, and US 2005/0046758 A1.
  • dendrimers WO 2005/040094 A1
  • large biomolecules US 2005/0042363 A1
  • collect information on samples such as identifying a carcinogenic or pathogenic cell in an individual (WO 2005/016230 A2).
  • the gene expression profile of the endometrium during the implantation window in stimulated cycles has also been analyzed (Mirkin and cois., 2004; Horcajadas and cois., 2005; Simón C and cois., 2005) and in response to drugs such as RU486 ( Catalano and cois., 2005; Sharkey and cois., 2005).
  • the refractory profile of the human endometrium in the presence of an intrauterine device (IUD) during the implantation window has also been studied (Horcajadas et al. 2006). All these works have been recently reviewed by the authors of this application (Horcajadas y cois., 2007).
  • a method is described to diagnose the endometriosis, based on the identification of the product of at least one of the genes of the group consisting of fibronectin, PTK7 transmenbrane receptor, type XVIII collagen, alpha 1 , subtilisin-like protein (PACE4), M-laminin chain (merosine), elastin, type IV collagen, alpha 2, inducible p27 interferon alpha gene, reticulocalbin, aldehyde dehydrogenase 6, gravine, nidogen and phospholipase C epsilon, in which A small amount of the control gene indicates the presence of endometriosis.
  • US 2003/0186300 A1 describes commercial methods and compositions for the diagnosis and treatment of diseases associated with reproduction.
  • the invention also relates to methods and compositions for determining and modulating the receptivity of the endometrium.
  • Ia of its isoforms in the uterus This enzyme is useful in the control of fertility, to monitor a premature pregnancy and to detect the receptivity of the uterus in mammals. New forms of 5/6 proprotein convertase are also described,
  • WO 01/89548 A2 refers to the pharmaceutical use of the polypeptide and the nucleic acid of fibulin-1 in the control of birth in women, and for the diagnosis and treatment of endometriosis.
  • the invention relates to a method and means for determining the specific conditions or changes in the mucosa of the uterus or in the epithelium of other organs.
  • the method allows to determine the overexpression of mRNA subunits of the type l- ⁇ ( ⁇ 7, ⁇ 6, ⁇ 6e) of human gonadotropin.
  • the measures of the expression of ⁇ 7, ⁇ 6, ⁇ 6e are used to indicate the receptivity of the uterine mucosa to the implantation of an embryo or to indicate neoplastic changes in the epithelium.
  • US 6,733,962 B2 describes a method to diagnose an abnormal development of a woman's endometrium based on the expression of cyclin E and p27 in a sample obtained after day 20 of the menstrual cycle of a woman that ideally lasts 28 days .
  • microarray available that encompasses a selection of genes that generate an expression profile that serves to diagnose and determine whether the state of a particular endometrium corresponds to the state of receptivity / non-receptivity.
  • the present invention consists in determining the functional state of human endometrial receptivity by means of the invention of a specific microarray that identifies the gene expression profile of the situation of receptivity / non-receptivity of the human endometrium.
  • the basis of the microarray is as follows: when a gene is active, many mRNA molecules are produced, which have the same sequence of complementary bases as the gene. When a gene is inactive, mRNA is not produced.
  • the researchers extract the whole mRNA from two cell populations that vary in the situation to be studied, in this case the receptive and non-receptive endometrium, mark it with a fluorescent substance and add it to the microarray. Since each mRNA is matched only to the probe of the gene that has its same sequence of complementary bases, those points that capture more mRNA - and therefore shine with more fluorescence - will indicate which genes were most active. If we compare the fluorescence pattern of the receptive endometrium against the non-receptive one, we will know which genes are differentially expressed in one situation with respect to the other, and that we therefore assume that they are related to the process.
  • the probes are designed to be joined by the mRNAs of the gene to which they belong, and are those that are attached to the array support.
  • the oligonucleotides that form the probe are automatically inserted in a layer of glass, nylon or plastic, being placed in boxes that act as a test microtube.
  • Oligonucleotide microarrays are manufactured automatically and inserted by robots using photolithography or piezoelectric printing. The result is an automated and standardized process that allows thousands of impressions per cm 2 and minute.
  • the distribution of the probes in the microarray is observed as a set of probes; those that have the same sequence are located in the matrix in a same point.
  • the probes have 60 nucleotide oligos, anyway, it is a technology in constant improvement.
  • what is marked and loose in the solution that we hybridize in the microarray are these mRNAs that will be the ones that bind to the probe attached to the support in the manner explained, by sequence homology, so that how much more labeled mRNA binds to a point, which corresponds to a specific probe, more light will be detected at that point, and more quantity will be from that gene.
  • this microarray for molecular diagnosis, it can also be used as a biotechnological tool for the study of the possible effect of drugs on the endometrium, such as the response to contraceptive drugs, both in vitro and in vivo tests.
  • Fig 1. List of the 293 genes and 738 probes of the genes included in the ERA. In “bold” the 567 probes corresponding to the 238 genes with an FDR ⁇ 0.05 and an FC> 3 are shown, which are the ones that have been selected and are specified in Figure 2.
  • Fig 2. List of the 238 genes selected with an FDRO, 05 and an FC> 3 that will be used by the bioinformatic predictor to classify.
  • Fig 3. Specific microarray (ERA) (agilent technologies). The figure shows how the ERA, oligo expression microarray, has an 8x15K format (8 arrays of 15,000 probes) per slide. Each array contains 15,744 points: 738 probes that include the selected genes (8 replicates, 5,904 points), 536 control points and 9,304 free (empty) points.
  • Fig 6. Different classifiers called SVM (Support Vector Machine), KNN (Neighboring Neighborhood), DLDA (Diagonal Linear Discriminant Analysis), PARM, SOM (Self-Organizing Map), which are mathematical algorithms that use different procedures to classify the profiles of expression as LH + 7 (receptive) and OTHER (non-receptive).
  • SVM Small Vector Machine
  • KNN Neighboring Neighborhood
  • DLDA Diagonal Linear Discriminant Analysis
  • PARM Spagonal Linear Discriminant Analysis
  • SOM Self-Organizing Map
  • the figure represents the classification error for each classifier in the axis of the ordinates, and the number of probes used grouped in order of significance, in the axis of the abscissa.
  • Fig 7. Representation of the error according to the number of probes. Two examples of the classification obtained with two of the classifiers are shown. In the case in which the classification of the algorithm coincides with the real class, it is represented in gray, and if it is not framed in a circle.
  • Fig 9. Shows the comparison of the results with the predictions obtained with the ERA against those obtained by morphological dating (Noyes).
  • the first phase of the project consists in the identification of the genes that are specifically regulated in the endometrium of the LH + 7 day and that will be part of the personalized microarray.
  • FDR 0.05 P-value depending on the size of the sample.
  • the value of FDR 0.05 is the significance that is taken into account at the statistical level and involves the risk of 5% that the differences are due to chance, and not to the biological process in question.
  • FC from the acronym of FoId change. It means: number of times the expression of a gene changes in one situation with respect to the other.
  • FO 3 the criterion is to assume that if it changes more than three times it is enough change to consider the gene important for the process.
  • the FDR has considered the possibility that differences in expression may be due to chance, and not to the biological process.
  • genes with an Fc above a threshold value of 3 have been selected so that the final number of genes with which we will work is feasible. In this way, we are giving more importance to the genes that change the most, because we assume a directly proportional relationship between greater change and greater importance for the process.
  • This strict criterion combines both the statistical and the biological requirement.
  • the final list of the ERA is composed of 293 genes represented by 738 probes plus controls.
  • Figure 1 shows in "bold" the 567 probes corresponding to the 238 genes with an FDR ⁇ 0.05 and an FC> 3, which are those that have been selected and are specified in Figure 2.
  • the ERA is an oligo expression microarray with an 8x15K format (8 arrays of 15,000 probes) per slide (Figure 3). Each array contains 15,744 points: 738 probes that include our selected genes (8 replicates, 5,904 points), 536 control points and 9,304 free (empty) points.
  • RNA thus extracted is DNAse for 1 hour at 37 0 C to remove traces of DNA and purified again using the Qiagen RNeasy kit following the manufacturer 's instructions.
  • the RNA that is obtained behind the columns of the RNeasy kit is analyzed to verify its quality in the Agilent 2100 bioanalyzer using the specific RNA chips of the Agilent brand, the Nano LabChip RNA.
  • RNA complementary DNA single strand cDNA
  • cDNA total RNA complementary DNA single strand
  • cRNA complementary RNA
  • That cRNA is purified by means of a purification kit based on affinity chromatography, and quantified.
  • that labeled cRNA is fragmented for 30 min at 60 ° C and hybridized in the microarray for 17 hours at 65 ° C. After that time, the microarray is washed to eliminate nonspecific hybridizations. Once hybridized and washed, the microarrays are centrifuged at 3,000 rpm for 3 minutes to dry the microarrays and are read by scanning on an Axon GenePix 4100A reading for the intensities of Cy3 (532nm).
  • the correction of the background effect has been made by subtracting half of the median from this, to the average of the intensity of the point.
  • Interarray normalization has been done using the quantile method.
  • results obtained in the ERA have been validated by quantitative PCR in order to give more consistency to our results and verify that the microarray analysis is reliable.
  • RNA in the form of cDNA was placed in the presence of 1 ⁇ g oligo (dT) (Clontech) until reaching a final volume of 12.5 ⁇ l with water treated with DEPC (diethyl pyrcarbonate). It was heated 2 minutes at 7O 0 C so that any possible secondary structure in the mRNAs was denatured and then kept on ice for 2 minutes.
  • dT oligo
  • DEPC diethyl pyrcarbonate
  • 6.5 ⁇ l was added for each of the 30 samples to be validated, of a MIX solution with 4 ⁇ l of buffer, 1 ⁇ l dNTP, 0.5 ⁇ l RNAse and 1 ⁇ l of retrotranscriptase (Rt-PCR clontech).
  • the 1 hr retrotranscription was continued in the thermal cycler.
  • 80 ⁇ l of water with DEPC are added and the concentration of the single stranded cDNA obtained by placing 2 ⁇ l of sample and 98 ⁇ l of DEPC water is measured by spectrophotometry.
  • the amount of cDNA that has been retrotranscribed must be between 80 to 120 ng / ⁇ l to start from similar concentrations, although it is normalized with the internal standard in our case GAPDH.
  • the range of cDNA to be amplified must be between 50-500 ng / ⁇ l. If any sample is not in those parameters, it is diluted.
  • a predictor is a mathematical tool that uses a data matrix, in our case those generated with the ERA, and learns to distinguish classes (Medina I 1 and cois, 2007), in our case two classes (LH + 7 (receptive) / OTHER (not receptive)).
  • the reasoning underlying this strategy is as follows: if we can distinguish between classes as a result of the different level of gene expression, then in theory, it is possible to find the characteristic gene expression of LH + 7 and use it to assign a class to the profile of Expression of the sample problem analyzed with the ERA custom microarray.
  • the set of samples that trains the classifier to define the classes is called a training set. That is, the gene expression profiles of these samples, measured with the ERA, are used by the program to know which probes are the most informative to predict (LH + 7 against OTHERS).
  • the biopsies used to generate the classification model are carefully chosen and dated in the most reliable way we have today.
  • This training set consists of 14 samples: LH + 7 (8), LH + 2 (4), and days 9-10 (2). They are all independent samples, of different, healthy women, in natural cycle and with proven fertility. These are Caucasian women with a body mass index between 19 and 25 kg / m 2 and between 19 and 34 years.
  • the classification is carried out by the bioinformatics program using different mathematical algorithms (Figure 6).
  • An algorithm is a well-defined, orderly and finite list of operations that allows us to find the solution to a problem. Given an initial state and an entry, through successive and well-defined steps a final state is reached, obtaining a solution.
  • the classifier manages to distinguish the classes according to a procedure called cross validation. The efficiency of the classifier is calculated and we generate classifications with error 0, that is, that all the samples of the training set are classified in the assigned real class.
  • the DLDA establishes a separation by virtue of a straight line
  • the KNN based on the nearest neighboring point is based on distances ... and thus each method is based on a mathematical separation criterion that will be adjusted more or less to The reality of the samples.
  • 16 endometrial samples considered as the test set have been analyzed with the ERA: 8 samples in: day 20 (2), day 21 (2 ), day 22 (1) and LH + 7 (3), and 8 of different days of the cycle: day 8 (2), day 14 (1), day 9 (1), day 10 (1), day 11 ( 1), day 16 (1), and day 26 (1). It is also about endometria of healthy women and of proven fertility, Caucasian, with an ICM between 19-25 kg / cm2, and an age between 19 and 34 years. Only those samples whose histological dating applying Noyes criteria have been chosen, coincide between the two pathologists and with the day of the menstrual cycle.
  • a value is assigned to what it predicts in the test set. If the error is 0 for all models, the prediction value is the same, but if not, the relative contribution is calculated to give the prediction greater or lesser weight. There is an inversely proportional relationship between the error and what should contribute to the decision. If the errors are relativized to one, and that value given by each method is subtracted from one, the direct relationship is achieved, since at greater error, the subtraction with respect to one, gives a lower value, and with that, the Methods with more error have a lower weighting coefficient, and the contribution to the decision to predict, LH + 7 or ANOTHER will be lower.
  • the value of 1 is given to LH + 7, and to another the value of -1, the weighting coefficient is multiplied according to the method error, by 1 or -1 as predicted. In the end all the methods for each array are compared, and the summation is made, if the value it gives is> 0, Io classifies it as LH + 7, and if the value is ⁇ 0, as ANOTHER.
  • the prediction establishes a value different from that classified according to Noyes in samples 7 (LH + 7), 14 (day 26), and 15 (day 16), but this is not considered an error, since these samples even correspond to one day of the cycle in the first receptive case and in the second and third non-receptive cases, they may have a different gene expression profile, as predicted by the model according to the expression profile.
  • Figure 9 shows the comparison of the results with the predictions obtained with the ERA against those obtained by morphological dating (Noyes).
  • Noyes dating has a higher error rate in the days that comprise
  • sample 7 of the array it is a patient LH + 7 according to Noyes, and
  • the ERA can be used for the positive identification of endometrial receptivity as well as for the diagnosis of its alteration associated with endometrial alterations typical of pathologies such as endometriosis, implantation failure, hydrosalpinx, etc.
  • this diagnostic tool would allow to detect the functional modifications induced by interceptive drugs or that aim to improve endometrial receptivity, altering the normality / abnormality situation in the receptive profile of a woman's endometrium.

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Health & Medical Sciences (AREA)
  • Organic Chemistry (AREA)
  • Wood Science & Technology (AREA)
  • Analytical Chemistry (AREA)
  • Zoology (AREA)
  • Genetics & Genomics (AREA)
  • Engineering & Computer Science (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • Microbiology (AREA)
  • Molecular Biology (AREA)
  • Biotechnology (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

La presente invencion se refiere a la determiinacion de la receptividad del endometrio humano a partir de un perfil de expresión genética. Más concretamente, consiste en la elaboración de un microarray de expresión especifico de receptividad endometrial (ERA o Endometrial Receptivity Array) que permita evaluar el estado receptivo de un emdometrio humano, así como valorar el estado del endometrio con fines diagnóstico y terapéuticos.

Description

PERFIL DE EXPRESIÓN GÉNICA COMO MARCADOR DE LA RECEPTIVIDAD
ENDOMETRIAL
Sector de Ia técnica
La presente invención se refiere a Ia determinación de Ia receptividad del endometrio humano a partir de un perfil de expresión génica. Más concretamente, consiste en Ia elaboración de un microarray de expresión específico de receptividad endometrial (ERA o Endometrial Receptivity Array) que permita evaluar el estado receptivo de un endometrio humano, así como valorar el estado del endometrio con fines diagnósticos y terapéuticos.
Estado de Ia técnica
El endometrio es Ia mucosa que recubre el interior de Ia cavidad uterina. Su función es Ia de alojar al embrión, permitiendo su implantación y favoreciendo el desarrollo de Ia placenta. Este proceso requiere de un endometrio receptivo capaz de responder a las señales del blastocito que es el estadio de desarrollo en el que se encuentra el embrión cuando implanta. El endometrio humano es un tejido regulado hormonalmente de forma cíclica, las hormonas que Io preparan para alcanzar dicho estado de receptividad son el estradiol que induce Ia proliferación celular y Ia progesterona que está implicada en Ia diferenciación, provocando un gran número de cambios en el perfil de expresión génica del endometrio que alcanza un fenotipo receptivo durante un corto periodo de tiempo llamado "ventana de implantación". Aunque no hay total acuerdo sobre el periodo de implantación en los humanos, los estudios clínicos sugieren que este proceso tiene lugar entre los días 20 y 24 de un ciclo ovulatorio normal (Wilcox y cois., 1999), considerándose crítico el día LH+7 (día 20-21).
La evolución de nuestro conocimiento sobre el endometrio humano contrasta con Ia ausencia de progreso en el desarrollo de nuevos métodos diagnósticos para su datación y estudio. En Ia actualidad, Ia valoración del endometrio todavía se realiza por medio de estudios histológicos basados en observaciones descritas hace más de 50 años (Noyes y cois., 1950) o con técnicas macroscópicas y poco resolutivas como estudios ecográficos igualmente poco objetivos, faltos de concreción y que producen resultados con grandes variaciones.
En 1950, Noyes y colaboradores describieron por primera vez, un método de datación endometrial basado exclusivamente en criterios histológicos y en los cambios morfológicos de los distintos compartimentos del endometrio en respuesta a Ia presencia de estrógenos y progesterona. Noyes y sus colaboradores examinaron los rasgos histológicos de biopsias endometriales tomadas durante 8.000 ciclos espontáneos en 300 mujeres (Noyes y cois., 1950; 1953). Fueron capaces de relacionar distintos patrones histológicos con momentos particulares del ciclo menstrual correlacionando los cambios histológicos con Ia temperatura corporal basal. Estos criterios morfológicos siguen siendo usados hoy en día y se consideran como el "GoId Standard' para el estudio del endometrio, evaluación de Ia receptividad endometrial y detección de anomalías endometriales.
Sin embargo, esta técnica no está exenta de inconvenientes. Se demuestra que el uso de las características histológicas falla a Ia hora de distinguir Ia fase del ciclo menstrual, así como forma de discriminar entre mujeres fértiles e infértiles, concluyéndose que no es útil para su uso clínico. La subjetividad que supone Ia observación visual hace que exista una variabilidad interobservador, intraobservador e interciclo que alteran Ia consistencia de los resultados obtenidos. Además, Ia estimulación ovárica propia de los tratamientos de reproducción asistida (TRA) produce modificaciones en el proceso de maduración endometrial comparado con ciclos naturales que difícilmente pueden ser explicados con los criterios de Noyes (Papanikolaou y cois., 2005). Es por esto, que existe en Ia literatura muchos trabajos que cuestionan las observaciones histológicas interpretadas por uno o varios patólogos, tanto en estudios clínicos retrospectivos (Balash y cois., 1992; Batista y cois., 1993; Shoupe y cois., 1989), prospectivos (Li y cois., 1989; Creus y cois., 2002; Ordi y cois., 2003), así como recientemente en estudios aleatorios (Murray y cois., 2004; Coutifaris y cois., 2004). Asimismo, el comité práctico de Ia Sociedad Americana de Medicina Reproductiva (ASRM) establece que aunque el criterio clásico del defecto de Ia fase lútea consiste en un retraso en Ia maduración endometrial de >2 días siguiendo los criterios de Noyes, este comité tiene serias dudas sobre Ia exactitud de dichos criterios histológicos y por Io tanto de Ia prevalencia del defecto de Ia fase lútea (DFL) e incluso de su relevancia clínica como causa de infertilidad (ASRM, 2000).
En este sentido, Balasch y cois., demostraron que Ia incidencia de DFL y los patrones endometriales histológicos eran similares en mujeres fértiles e ¡nfértiles. Más aún, una histología endometrial adecuada en el ciclo de Ia ovulación o en los anteriores no estaba relacionada con los datos de embarazo en mujeres infértiles concluyendo que Ia evaluación histológica del endometrio en Ia fase lútea no resulta útil para predecir o mejorar los resultados reproductivos (Balasch y cois., 1992). En otros estudios del mismo grupo se demostró Ia existencia de una clara disociación en Ia expresión temporal de una serie de marcadores relacionados con Ia ventana de implantación (integrinas alpha y beta 3) y Ia expresión de pinópodos. No encontraron, además, diferencias de expresión de estos marcadores entre mujeres fértiles e infértiles (Creus y cois., 2002). También demostraron una alta variabilidad entre ciclos y una baja reproducibilidad para estos marcadores (Ordi y cois., 2003).
Li y cois., realizaron Ia datación de 63 biopsias endometriales en dos ocasiones distintas por el mismo patólogo demostrando que sólo en un 24% de ellas hubo un acuerdo total. En un estudio separado observaron que entre distintos ciclos de Ia misma mujer sólo en un 4% de los casos hubo acuerdo total. Estos datos enfatizan Ia falta de precisión de los métodos de datación tradicionales y su falta de garantía para predecir el desarrollo en los ciclos siguientes (Li y cois., 1989).
Las diferencias entre patólogos varían dependiendo del día del ciclo menstrual en que Ia biopsia endometrial es recogida. Más de un 20% de las biopsias endometriales fueron datadas con una diferencia de al menos dos días entre patólogos en Ia fase lútea temprana, media y tardía. Con respecto a las variaciones interciclos, éstas llegan a ser del 60% en Ia fase lútea media (Murray y cois., 2004). Se ha demostrado que, durante Ia ventana de implantación, un porcentaje muy similar de mujeres tienen el endometrio fuera de fase, 49,4% de fértiles frente al 43,2% de infértiles (p = 0,33) y que, finalmente, Ia datación histológico no está relacionado con el estatus de fertilidad (Coutifaris y cois., 2004). Estas variaciones descritas sugieren que los criterios tradicionales no son precisos y que se requiere nuevas tecnologías para datar e identificar funcionalmente las muestras endometriales. En Ia era pre-genómica, sólo se podían llevar a cabo estudios "gen a gen" para seleccionar candidatas útiles para el estudio de Ia receptividad uterina o para determinar Ia situación endometrial en mujeres con endometriosis o sin ella.
Así, en Ia era genómica en Ia que nos encontramos, se buscan herramientas objetivas basadas en criterios moleculares que vengan a mejorar Ia capacidad diagnóstica de determinadas técnicas como Ia histológica, muy útil, sin embargo, para otro tipo de necesidades.
A mediados de los 90 (Schena y cois., 1995), se desarrolló una tecnología revolucionaria para Ia determinación y cuantificación de Ia expresión de los ARN mensajeros (ARNm) en una muestra, los microarrays de expresión génica. Su principal ventaja es que ofrecen Ia posibilidad de analizar simultáneamente miles de genes en un solo experimento. Un microarray de ADN consiste en un gran número de moléculas de ADN ordenadas sobre un sustrato sólido de manera que formen una matriz de secuencias en dos o tres dimensiones. Estos fragmentos de material genético pueden ser secuencias cortas llamadas oligonucleótidos o de mayor tamaño, como es el ADN complementario (ADNc) que es sintetizado a partir de ARNm, o bien productos de PCR (replicación in vitro de secuencias de ADN mediante Ia reacción en cadena de polimerasa). A estos fragmentos de nucleótidos de una sola hebra inmovilizados en el soporte, se les denomina "sondas". Los ácidos nucleicos de las muestras a analizar se marcan por diversos métodos (enzimáticos, fluorescentes, etc.) y se incuban sobre el panel de sondas, permitiendo Ia hibridación (reconocimiento y unión entre moléculas complementarias) de secuencias homologas. Durante Ia hibridación, las muestras de material genético marcadas se unen a sus complementarias inmovilizadas en el soporte del chip, permitiendo Ia identificación y cuantificación del ADN presente en Ia muestra. Después, el escáner y las herramientas bioinformáticas adecuadas permiten interpretar y analizar los datos obtenidos (Al-Shahrour F y cois., 2005).
Para utilizar un microarray se puede recurrir a los que existen comercialmente, o se puede diseñar uno personalizado.
Para diseñar un microarray se deben realizar las operaciones siguientes:
a) Elección del tipo de sonda, oligos, ADNc... b) Mareaje de sondas o muestras: enzimático, fluorescente...
c) Material de soporte: vidrio, plástico, membranas...
d) Inmovilización de sondas: activa, pasiva, covalente...
e) Fabricación: impresión, síntesis ¡n situ...
f) Detección de Ia hibridación: escáner, fluorimetría...
g) Procesamiento de datos: software.
Esta tecnología se está aplicando al análisis de Ia expresión génica, secuenciación, seguimiento de terapias, medicina preventiva, toxicología de fármacos y diagnóstico molecular. Se ha descrito Ia fabricación de microarrays, denominados también bioarrays o biochips en varios documentos de patente, como por ejemplo, WO 2005/018796 A1 , US 2005/0048554 A1 , y US 2005/0046758 A1. También su utilización se ha aplicado a dendrímeros (WO 2005/040094 A1) y grandes biomoléculas (US 2005/0042363 A1) o para recabar información sobre muestras, como por ejemplo, identificar una célula cancerígena o patógena en un individuo (WO 2005/016230 A2). También se conoce su utilización para inmovilizar ácidos nucleicos que son complementarios de una variedad de genes, aplicándose al campo de Ia química, Ia biología, Ia medicina y los diagnósticos en medicina (US 6,821 ,724 B1). En Ia actualidad se están utilizando los microarrays para hacer comparaciones basándose en datos genómicos y para investigar distintos sistemas.
Existen diversas publicaciones de literatura patente y no patente sobre este tema. En Ia actualidad, Ia tecnología del microarray ha permitido estudiar de forma global Ia expresión génica del endometrio bajo condiciones fisiológicas durante las diferentes fases del ciclo menstrual en ciclo natural (Ponnampalam y cois., 2004, Talbi y cois., 2005). En Io que respecta a Ia ventana de implantación humana los perfiles de expresión génica del endometrio en el ciclo natural han sido descritos (Borthwick y cois., 2003; Carson y cois., 2002; Riesewijk y cois., 2003;Mirkin y cois., 2005). También se ha analizado el perfil de expresión génica del endometrio durante Ia ventana de implantación en ciclos estimulados (Mirkin y cois., 2004; Horcajadas y cois., 2005; Simón C y cois., 2005) y en respuesta a fármacos como RU486 (Catalano y cois., 2005; Sharkey y cois., 2005). Asimismo se ha estudiado el perfil refractario del endometrio humano en presencia de un dispositivo intrauterino (DIU) durante Ia ventana de implantación (Horcajadas y col. 2006). Todos estos trabajos han sido revisados recientemente por los autores de Ia presente solicitud (Horcajadas y cois., 2007). La conclusión de dicho estudio es que aunque en los últimos 4 años, se han llevado a cabo diversos estudios genómicos del endometrio humano en distintas condiciones fisiológicas y patológicas, los cuales han generado gran cantidad de información sobre Ia regulación de los genes durante Ia ventana de implantación, tanto en mujeres fértiles como ¡nfértiles, sin embargo, las moléculas y mecanismos clave todavía están por descubrir.
En el campo de las patentes, existen varias que tratan de determinar Ia receptividad/no receptividad del endometrio, aunque ni los genes, ni Ia tecnología, ni los sistemas predictivos que postulan coinciden con Ia actual patente.
En el documento de patente US 2003/0077589 A1 se describe un método para diagnosticar Ia endometriosis, basado en Ia identificación del producto de al menos uno de los genes del grupo que consiste en fibronectina, receptor de transmenbrana PTK7, colágeno tipo XVIII, alfa 1 , proteína similar a Ia subtilisina (PACE4), cadena de laminina M (merosina), elastina, colágeno tipo IV, alfa 2, gen p27 interferón alfa inducible, reticulocalbina, aldehido deshidrogenasa 6, gravina, nidogen y fosfolipasa C epsilon, en Ia que una pequeña cantidad del gen control indica Ia presencia de endometriosis.
En Ia solicitud de patente US 2003/0125282 A1 se describen dos proteínas humanas MATER (se conocían las MATER de ratón) y su relación y utilización para las alteraciones relacionados con Ia fertilidad.
En el documento US 2003/0186300 A1 se describen métodos y composiciones comerciales para el diagnóstico y el tratamiento de enfermedades asociadas a Ia reproducción. La invención también se refiere a métodos y composiciones para Ia determinación y modulación de Ia receptividad del endometrio.
En Ia patente US 2005/0032111 A1 se utiliza Ia expresión de Ia cadherina-11 en el tejido endometrial como indicador de Ia capacidad para establecer o mantener un embarazo. El documento US 2005/0106134 A1 se refiere al papel de Ia enzima proproteína convertasa 5/6 durante el embarazo, y particularmente de su detección y
Ia de sus isoformas en el útero. Esta enzima es útil en el control de Ia fertilidad, para monitorizar un embarazo prematuro y para Ia detección de Ia receptividad del útero en los mamíferos. También se describen nuevas formas de proproteína convertasa 5/6,
En Ia patente US 2003/0228636 A1 se describe un método para detectar Ia receptividad del endometrio para Ia implantación embrionaria, que comprende: obtener una muestra del endometrio, poner en contacto el endometrio con un anticuerpo monoclonal de β3, y detectar β3 en el endometrio. También se mencionan los contraceptivos y kits de diagnóstico útiles para llevar a cabo los métodos de Ia invención.
En Ia solicitud de patente WO 2005/061725 A1 se describen métodos para detectar marcadores asociados con enfermedades endometriales o una determinada fase endometrial en una mujer, que comprenden medir los marcadores endometriales peptídicos o los polinucleótidos que codifican los marcadores, en Ia muestra estudio. La invención también proporciona métodos para detectar enfermedades endometriales, así como kits para llevar a cabo los métodos de Ia invención.
El documento WO 01/89548 A2 se refiere a Ia utilización farmacéutica del polipéptido y el ácido nucleico de Ia fibulina-1 en el control de Ia natalidad en las mujeres, y para el diagnóstico y tratamiento de Ia endometriosis.
En Ia patente WO 2004/058999 A2 Ia invención se refiere a un método y los medios para determinar las condiciones específicas o los cambios en Ia mucosa del útero o en el epitelio de otros órganos. El método permite determinar Ia sobre- expresión de subunidades ARNm del tipo l-β (β7,β6,β6e) de gonadotropina humana. Las medidas de Ia expresión de β7,β6,β6e se utilizan para indicar Ia receptividad de Ia mucosa uterina a Ia implantación de un embrión o para indicar cambios neoplásicos en el epitelio.
En Ia patente US 2004/0005612 A1 se identifican secuencias genéticas con niveles de expresión que son reprimidos o inducidos en el endometrio humano durante Ia ventana de implantación. Los genes caracterizados durante Ia ventana de implantación proporcionan material para ensayos con el objeto de determinar alteraciones del endometrio e infertilidad, así como métodos de control de Ia natalidad basados en el endometrio.
En Ia patente US 6,733,962 B2 se describe un método para diagnosticar un desarrollo anormal del endometrio de una mujer basándose en Ia expresión de Ia ciclina E y Ia p27 en una muestra obtenida después del día 20 del ciclo menstrual de una mujer que idealmente dura 28 días.
Resumiendo, durante más de 50 años se ha tratado de determinar un estándar histológico para utilizar en el diagnóstico clínico de receptividad endometrial basado en observaciones morfológicas. Hoy día, con Ia tecnología de los microarrays, mucho más precisa que las observaciones morfológicas, se han publicado trabajos que se refieren a distintos genes presentes a Io largo del ciclo menstrual, pero los resultados no coinciden ya que el diseño experimental, Ia recolección de muestras y Ia selección de los genes son cruciales a Ia hora de sacar conclusiones.
Por Io tanto, todavía y más que nunca, es necesario tener disponible un microarray que englobe una selección de genes que generen un perfil de expresión que sirva para diagnosticar y determinar si el estado de un endometrio particular corresponde al estado de receptividad/no receptividad.
Por ello, en Ia presente solicitud se ha determinado una lista de genes y sondas, los cuales una vez incorporados a un microarray, por medio del estudio de Ia expresión de estos genes en Ia muestra objeto del estudio, permite detectar el estado de receptividad/no receptividad de una muestra de endometrio obtenida 7 días después del pico de Ia LH, así como situaciones de sub-fertilidad de origen endometrial en función del perfil de expresión génica de todos ellos.
Descripción detallada de Ia invención
La presente invención consiste en determinar el estado funcional de receptividad endometrial humana por medio de Ia invención de un microarray específico que identifica el perfil de expresión génica de Ia situación de receptividad /no receptividad del endometrio humano.
Para ello, se siguen los pasos siguientes: 1. Identificar un conjunto de los genes que estén implicados en Ia receptividad endometrial para su inclusión en un microarray específico de receptividad endometrial (Endometrial Receptivity Array, ERA).
2. Creación del microarray específico.
3. Analizar el patrón de expresión del ERA durante Ia ventana de implantación mediante herramientas bioinformáticas, para poder establecer el perfil de receptividad endometrial y crear un clasificador.
4. Desarrollar con este clasificador un predictor que, basado en el perfil de expresión génica, permita evaluar y predecir, de forma cuantitativa y objetiva el estado receptivo endometrial in vivo.
El fundamento del microarray es el siguiente: cuando un gen está activo, se producen muchas moléculas de ARNm, que tienen Ia misma secuencia de bases complementarias que el gen. Cuando un gen está inactivo no se produce ARNm. Los investigadores extraen todo el ARNm de dos poblaciones celulares que varían en Ia situación a estudiar, en este caso el endometrio receptivo y no receptivo, Io marcan con una sustancia fluorescente y Io añaden sobre el microarray. Como cada ARNm se aparea sólo a Ia sonda del gen que tiene su misma secuencia de bases complementarias, aquellos puntos que capturen más ARNm -y que por tanto brillen con más fluorescencia- indicarán qué genes estaban más activos. Si comparamos el patrón de fluorescencia del endometrio receptivo contra el no receptivo, sabremos qué genes están expresados de forma diferencial en una situación con respecto de Ia otra, y que por tanto asumimos que están relacionados con el proceso.
Las sondas están diseñadas para que se unan a ellas los ARNm del gen al cual pertenecen, y son las que están fijadas en el soporte del array. Los oligonucleótidos que forman Ia sonda se insertan de forma automatizada en una capa de cristal, nylon o plástico, colocándose en unas casillas que actúan a modo de microtubo de ensayo. Los microarrays de oligonucleótidos son fabricados de manera automatizada e insertados por robots mediante fotolitografía o impresión piezoeléctríca. El resultado es un proceso automatizado y normalizado que permite miles de impresiones por cm2 y minuto.
Generalmente se observa Ia distribución de las sondas en el microarray como un conjunto de sondas; las que tienen Ia misma secuencia se sitúan en Ia matriz en un mismo punto. En nuestro array personalizado las sondas poseen oligos de 60 nucleótidos, de todos modos, se trata de una tecnología en constante perfeccionamiento. De este modo, Io que está marcado y suelto en Ia disolución que ponemos a hibridar en el microarray son estos ARNm marcados que serán los que se unan a Ia sonda fijada al soporte de Ia manera explicada, por homología de secuencia, de modo que cuanto más ARNm marcado se una a un punto, que corresponde a una sonda específica, más luz será detectada en ese punto, y más cantidad habrá de ese gen.
Una vez establecido el funcionamiento del microarray objeto de Ia invención, así como delimitado el patrón de expresión de receptividad para evaluar mediante métodos bioinformáticos el estado de receptividad/no receptividad de un endometrio, se podrá evaluar igualmente, utilizando el mismo método, otros procesos patológicos que deriven en infertilidad o subfertilidad de origen endometrial como fallos de implantación o endometriosis.
Además del uso de este microarray para el diagnóstico molecular, también se puede utilizar como herramienta biotecnológica para el estudio del posible efecto de fármacos en el endometrio, como por ejemplo, Ia respuesta a fármacos contraceptivos, tanto en ensayos in vitro como in vivo.
Breve descripción de las figuras
Fig 1. Lista de los 293 genes y de las 738 sondas de los genes que incluye el ERA. En "negrita" se muestran las 567 sondas correspondientes a los 238 genes con un FDR<0,05 y un FC>3, que son los que se han seleccionado y se especifican en Ia Figura 2.
Fig 2. Lista de los 238 genes seleccionados con un FDRO, 05 y un FC>3 que serán utilizados por el predictor bioinformático para clasificar.
Fig 3. Microarray específico (ERA) (agilent technologies). La figura muestra como el ERA, microarray de expresión de oligos, tiene un formato de 8x15K (8 arrays de 15.000 sondas) por portaobjetos. Cada array contiene 15.744 puntos: 738 sondas en las que se incluyen los genes seleccionados (8 réplicas, 5.904 puntos), 536 puntos control y 9.304 puntos libres (vacíos). Fig 4. Tabla en Ia que se muestran los cebadores directos e inversos diseñados de los genes a amplificar mediante PCR cuantitativa.
Fig 5. Media de Ia expresión de las sondas de cada gen en el array comparada con Ia expresión en Ia PCR cuantitativa.
Fig 6. Distintos clasificadores denominados SVM (Support Vector Machine), KNN (Nearest Neighbour), DLDA (Diagonal Linear Discriminant Análisis), PARM, SOM (Self-Organizing Map), que son algoritmos matemáticos que utilizan distintos procedimientos para clasificar los perfiles de expresión como LH+7 (receptivo) y OTRO (no receptivo). En Ia figura se representa el error de clasificación para cada clasificador en el eje de las ordenadas, y el número de sondas empleadas agrupadas por orden de significatividad, en el eje de las abscisas.
Fig 7. Representación del error en función del número de sondas. Se muestran dos ejemplos de Ia clasificación obtenida con dos de los clasificadores. En el caso en que Ia clasificación del algoritmo coincida con Ia clase real, se representa en gris, y si no se enmarca en un círculo. La figura 7a clasifica con un error de 0, con todas las agrupaciones de sondas, pero en el mostrado en 7b, se observa como con 200 sondas discrepa en 2 de 14 (error= 2/14 =0,14).
Fig 8. Tabla en Ia que se comparan todos los métodos para cada array y se clasifica como LH+7 o como OTRO. Se escogen los clasificadores con error 0 y se ve que existe una predicción consenso.
Fig 9. Muestra Ia comparación de los resultados con las predicciones obtenidas con el ERA contra los obtenidos por datación morfológica (Noyes).
Ejemplos
1. IDENTIFICACIÓN DE LOS GENES IMPLICADOS EN LA RECEPTIVIDAD ENDOMETRIAL PARA LA GENERACIÓN DEL MICROARRAY ESPECÍFICO DE RECEPTIVIDAD ENDOMETRIAL.
La primera fase del proyecto consiste en Ia identificación de los genes que se encuentran regulados de forma específica en el endometrio del día LH+7 y que van a formar parte del microarray personalizado.
En Ia mayoría de los trabajos publicados los genes mencionados se han seleccionado cuando se encuentran inducidos o reprimidos dos veces. En Ia presente invención se han seguido unos criterios de selección distintos y más estrictos:
Criterio de selección de genes.
La selección de genes ha sido realizada en base a las diferencias del perfil de expresión génica endometrial representado por LH+1 , LH+3 y LH+5 (no receptivo) contra LH+7 como estado receptivo. Los niveles de expresión se han obtenido a partir de un microarray de expresión de oligos de genoma completo. Se han escogido aquellos genes que muestran diferencias significativas de expresión en estas dos situaciones utilizando los criterios de FDR<0,05* y de FC>3**.
* FDR: de las siglas en inglés de False Discovery Rate. Parámetro que corrige el
P-valor en función del tamaño de Ia muestra. El valor de FDR 0,05 es Ia significatividad que de forma estándar se tiene en cuenta a nivel estadístico y supone el correr un riesgo del 5% de que las diferencias se deban al azar, y no al proceso biológico en cuestión.
** FC: de las siglas en ingles de FoId change. Significa: número de veces que cambia Ia expresión de un gen en una situación con respecto de Ia otra. En cuanto al FO 3, el criterio es suponer que si cambia más de tres veces es suficiente cambio para considerar el gen importante para el proceso.
Con el FDR se ha considerado Ia posibilidad de que las diferencias de expresión puedan deberse al azar, y no al proceso biológico. Además, se han seleccionado los genes con un Fc por encima de un valor umbral de 3 para que el número final de genes con los que vamos a trabajar sea factible. De este modo, estamos dando más importancia a los genes que cambian más, porque asumimos una relación directamente proporcional entre un mayor cambio y una mayor importancia para el proceso. Este criterio estricto auna tanto el requerimiento estadístico como el biológico. Además, hemos corroborado el sentido funcional de esta selección génica en el proceso biológico de Ia receptividad endometrial. Para ello, hemos realizado, mediante el uso de herramientas bioinformáticas, Ia clasificación ontológica de los genes utilizando FATIGO GEPAS (Al-Shahrour F y cois., 2005) viendo que los procesos biológicos representados de una manera superior a Ia esperada con una significatividad de 0,05, son Ia respuesta al estrés, Ia respuesta de defensa y Ia adhesión celular, procesos bastantes relevantes en preparar un endometrio a Ia posible implantación del blastocito.
Se han escogido aquellos genes con estas características y por medio de programas informáticos han resultando un total de 238 genes representados por 567 sondas.
La lista final del ERA esta compuesta por 293 genes representados por 738 sondas más controles. La Figura 1 muestra en "negrita" las 567 sondas correspondientes a los 238 genes con un FDR<0,05 y un FC>3, que son los que se han seleccionado y se especifican en Ia Figura 2.
2. CREACIÓN DEL MICROARRAY ESPECÍFICO (ERA) (Agilent Technologies)
El ERA es un microarray de expresión de oligos con un formato de 8x15K (8 arrays de 15.000 sondas) por portaobjetos (Figura 3). Cada array contiene 15.744 puntos: 738 sondas en las que se incluyen nuestros genes seleccionados (8 réplicas, 5.904 puntos), 536 puntos control y 9.304 puntos libres (vacíos).
Entrenamiento v predicción del ERA
En esta sección hemos utilizado los datos de expresión generados por el ERA para Ia clasificación de las muestras endometriales en dos grupos (LH+7 (receptivo) y OTRO (no receptivo) y para comprobar su capacidad de predicción.
Para ello, se obtuvieron 30 muestras de biopsias endometriales de mujeres fértiles (14 del conjunto de entrenamiento + 16 del conjunto de prueba). Todas las muestras independientes, de mujeres con fertilidad probada, en distintos días del ciclo menstrual. Se trata de mujeres caucásicas con un índice de masa corporal entre 19 y 25 kg/m2 y cuya edad oscila entre 19 y 34 años.
Dichas muestras se utilizaron para generar un clasificador (entrenamiento) y testar el valor de predicción del ERA.
Para ello, el ARN total fue extraído usando el protocolo del Trizol (Invitrogen) siguiendo las instrucciones del fabricante (Life Technologies, Inc., USA). Las muestras se homogeneizaron usando 1 mi de TRIzol por cada 75 mg de tejido, se incubaron a temperatura ambiente 5 min, se añadieron 200 μl volúmenes de cloroformo para Ia misma cantidad de tejido y se incubaron a temperatura ambiente 5 min. Posteriormente se centrifugaron 15 min a 12.000xg (4°C). La fase acuosa se precipitó con un volumen igual de 2-propanol (isopropanol), se incubaron en hielo 5 min y centrifugaron 30 min a 12.000xg (4°C). El precipitado se lavó con etanol 70% en agua tratada con dietilpirocarbonato (DEPC)1 para posteriormente resuspenderlo en agua DEPC (15 μl). Con este protocolo se suelen obtener 1-2 μg de ARN total por mg de tejido endometrial. El ARN así extraído se trata con ADNsa durante 1 hora a 370C para eliminar las trazas de ADN y purificarlo de nuevo usando el kit RNeasy de Qiagen siguiendo las instrucciones del fabricante. El ARN que se obtiene tras las columnas del kit RNeasy se analiza para comprobar su calidad en el bioanalizador Agilent 2100 usando los chips para ARN específicos de Ia marca Agilent, el ARN Nano LabChip.
Sólo se han utilizado para posteriores análisis aquellos ARNs que cumplían las características siguientes:
- No presentaban ADN genómico detectable,
- poseían una concentración superior a 200 μg/ml,
- el valor del radio de rARN era 28s/18S >1 ,2, y
- el valor de RIN>7,0, (ARN Integrity Number).
Tras los análisis, con las muestras seleccionadas por su calidad adecuada, a partir del ARN total se genera ADN complementarios de una sola cadena (ADNc) incubándolo entre una y dos horas a 400C con retrotranscriptasa, nucleótidos y un oligonucleótido polydT-T7 que porta, no solamente Ia secuencia poly T que híbrida con Ia cola poliA de los ARN mensajeros, sino además Ia secuencia de reconocimiento de Ia ARN polimerasa de T7.
A partir del ADNc obtenido en el paso anterior, se incuba durante 2 horas a 4O0C en presencia de ARN polimerasa de T7 y nucleótidos, uno de los cuales está marcado con Cy3, para producir ARN complementario llamado ARNc.
Se purifica ese ARNc por medio de un kit de purificación basado en una cromatografía de afinidad, y se cuantifica.
Una vez purificado, ese ARNc marcado se fragmenta durante 30 min a 60°C y se híbrida en el microarray durante 17 horas a 65°C. Una vez transcurrido ese tiempo, se lava el microarray para eliminar las hibridaciones inespecíficas. Una vez hibridados y lavados, los microarrays se centrifugan a 3.000 rpm durante 3 minutos para secar los microarrays y se procede a su lectura por medio de su escaneo en un Axon GenePix 4100A leyendo para las intensidades de Cy3 (532nm).
Como resultado, tras el pertinente procesado de datos que adjuntamos a continuación, se generó una matriz de expresión génica cuyas filas corresponden a las 567 sondas de los 238 genes seleccionados y cuyas columnas, a las distintas muestras.
Procesado de los datos del arrav
La corrección del efecto del fondo se ha realizado restando Ia mitad de Ia mediana de este, a Ia media de Ia intensidad del punto. La normalización interarray se ha hecho usando el método de los cuantiles.
Después es calculada Ia mediana de las ocho réplicas de cada sonda. Las diferentes sondas del mismo gen son analizadas individualmente y los resultados son procesados por herramientas bioinformáticas. Validación de los resultados del ERA mediante PCR
Los resultados obtenidos en el ERA han sido validados mediante PCR cuantitativa con objeto de dar mayor consistencia a nuestros resultados y comprobar que el análisis por microarrays es fiable.
Se lleva a cabo Ia retrotranscripción para obtener el ARN en forma de ADNc, para ello 1μg de ARN total se puso en presencia de 1μg oligo (dT) (Clontech) hasta llegar a un volumen final de 12,5 μl con agua tratada con DEPC (dietilpirocarbonato). Se calentó 2 minutos a 7O0C para que se desnaturalizara cualquier posible estructura secundaria en los ARNm y luego se mantuvo en hielo 2 minutos.
Después se añadió 6,5 μl por cada una de las 30 muestras a validar, de una solución MIX con 4 μl de tampón, 1 μl dNTP, 0,5 μl ARNsa y 1 μl de retrotranscriptasa (Rt-PCR clontech). Se prosiguió Ia retrotranscripción 1 hr en el termociclador. Se añaden 80 μl de agua con DEPC y se mide por espectrofotometría Ia concentración del ADNc de simple cadena obtenido poniendo 2 μl de muestra y 98 μl de agua DEPC. La cantidad de ADNc que se ha retrotranscrito debe estar entre 80 a 120 ng / μl para partir de concentraciones similares, aunque se normaliza con el patrón interno en nuestro caso GAPDH. De todos modos, para que Ia PCR cuantitativa funcione correctamente el rango de ADNc a amplificar debe estar entre 50-500 ng / μl. Si alguna muestra no está en esos parámetros, se diluye.
Diseñamos los cebadores directos e inversos para cinco genes aumentados en
LH+7 (Figura 4). Las secuencias de oligonucleótidos de los cebadores fueron diseñadas con el programa bioinformático Gene fisher (ver Figura 4 y listado de secuencias). El sistema de detección se realizó con el SYBR Green I de unión a ADN de doble cadena (Roche). Este sistema de detección establece un rango dinámico linear para detectar productos específicos de PCR. Todos los experimentos de Q-PCR se realizaron usando el SYBR Green PCR Master Mix (roche) y las condiciones universales de los parámetros de ciclos térmicos indicados por los fabricantes utilizando el Light cycler de Roche. Se realizaron 40 ciclos. Las temperaturas a las que funcionan bien los cebadores se pueden observar en Ia figura 4, La cuantificación relativa se realizo mediante el método estándar de Ia curva patrón. La expresión de GPX3; CLDN10; FXYD2 ; SPP1 ;y MT1G, corresponden en el ERA a los valores de expresión de las siguientes sondas:
Figure imgf000018_0001
Considerando que se trata de técnicas distintas, Ia PCR cuantitativa, cuya sensibilidad es mucho mayor pero que proporciona un solo valor de expresión, y los arrays en los que se tiene Ia expresión de distintas sondas para un mismo gen. Por ello, para poder comparar se ha realizado en el array Ia media de Ia expresión de las distintas sondas de un gen (Figura 5).
Debido a su distinta sensibilidad, se considera que Ia relación del valor de expresión entre ambas técnicas correspondería a un factor de corrección de 10 (aumentos de expresión de 1Ox en el array se admite que se corresponden con un máximo de 10Ox en Ia PCR cuantitativa (Figura 5).
3. ANALIZAR EL PATRÓN DE EXPRESIÓN DEL ERA DURANTE LA VENTANA DE IMPLANTACIÓN PARA PODER ESTABLECER EL PERFIL DE RECEPTIVIDAD ENDOMETRIAL. GENERACIÓN DE UN CLASIFICADOR .
Entrenamiento:
Un predictor es una herramienta matemática que utiliza una matriz de datos, en nuestro caso los generados con el ERA, y aprende a distinguir clases (Medina I1 y cois, 2007), en nuestro caso dos clases (LH+7 (receptivo)/OTRO (no receptivo)). El razonamiento que subyace a esta estrategia es el siguiente: si podemos distinguir entre las clases como consecuencia del distinto nivel de expresión génica, entonces en teoría, es posible encontrar Ia expresión génica característica de LH+7 y usarla para asignar una clase al perfil de expresión de Ia muestra problema analizada con el microarray personalizado ERA.
Al conjunto de muestras que entrena al clasificador para definir las clases, se denomina conjunto de entrenamiento ("training set"). Es decir, los perfiles de expresión génica de estas muestras, medidos con el ERA, son utilizados por el programa para saber que sondas son las más informativas para predecir (LH+7 contra OTROS). Las biopsias utilizadas para generar el modelo de clasificación, están cuidadosamente escogidas y datadas de Ia manera más fiable de Ia que disponemos en Ia actualidad. Este conjunto de entrenamiento consta de 14 muestras: LH+7 (8), LH+2 (4), y días 9- 10 (2). Son todas muestras independientes, de mujeres distintas, sanas, en ciclo natural y con fertilidad probada. Se trata de mujeres caucásicas con un índice de masa corporal entre 19 y 25 kg/m2 y entre 19 y 34 años.
La clasificación es realizada por el programa bioinformático utilizando distintos algoritmos matemáticos (Figura 6). Un algoritmo es una lista bien definida, ordenada y finita, de operaciones que permite hallar Ia solución a un problema. Dado un estado inicial y una entrada, a través de pasos sucesivos y bien definidos se llega a un estado final, obteniendo una solución. El clasificador consigue distinguir las clases según un procedimiento llamado validación cruzada. Se calcula Ia eficacia del clasificador y generamos clasificaciones con error 0, Es decir, que todas las muestras del conjunto de entrenamiento son clasificadas en Ia clase real asignada.
En Ia Figura 6, se muestran los distintos clasificadores denominados SVM
(Support Vector Machine), KNN (Nearest Neighbour), DLDA (Diagonal Linear Discriminant Análisis), PARM, SOM (Self-Organizing Map), que son algoritmos matemáticos que utilizan distintos procedimientos para clasificar los perfiles de expresión como LH+7 y OTRO. A priori no se puede saber como se distribuyen los datos en el espacio, sólo se sabe como se ubican en las dimensiones que se pueden distinguir, que son tres. De este modo, hay diferentes algoritmos a aplicar, que funcionarán mejor o peor según como los datos introducidos se distribuyan en el espacio. Se aplican los más usados en matemáticas, que son los arriba indicados, y vemos cual separa mejor en dos grupos (LH+7 y OTRO). De este modo el DLDA establece una separación en virtud a una recta, el KNN en función del punto vecino más cercano, se basa en distancias... y así cada método se basa en un criterio de separación matemático que se ajustará más a menos a Ia realidad de las muestras.
En Ia Figura 6 se representa el error de clasificación para cada clasificador en el eje de las ordenadas, y el número de sondas empleadas agrupadas por orden de significatividad, en el eje de las abcisas. Si de 14 muestras, clasifica todas como Ia clase real, el error es 0/14=0, si clasifica distinto 1 de las 14, el error será 1/14=0,07.
Otra manera de representar el error en función del número de sondas se muestra en las listas de figura 7 a y b. En dichas figuras se muestran dos ejemplos de Ia clasificación obtenida con dos de los clasificadores. Si Ia clasificación del algoritmo coincide con Ia clase real, se representa en gris, y si no se enmarca en un círculo. En el ejemplo que muestra 7 a de Ia izquierda clasifica con un error de 0, con todas las agrupaciones de sondas, pero en el mostrado en 7b, vemos como con 200 sondas discrepa en 2 de 14 (error= 2/14 =0,14). 4. DESARROLLO DE UN PREDICTOR QUE PERMITA EVALUAR Y PREDECIR, DE FORMA CUANTITATIVA Y OBJETIVA EL ESTADO RECEPTIVO ENDOMETRIAL BASÁNDOSE EN EL PERFIL DE EXPRESIÓN GÉNICA.
Testado:
Para comprobar el poder de predicción de los distintos modelos de clasificación, han sido analizadas con el ERA, 16 muestras endometriales consideradas como el conjunto de prueba ("test set"): 8 muestras en: día 20 (2), día 21 (2), día 22 (1) y LH+7 (3), y 8 de diferentes días del ciclo: día 8 (2), día 14 (1), día 9 (1), día 10 (1), día 11 (1), día 16 (1), y día 26 (1). También se trata de endometrios de mujeres sanas y de fertilidad probada, caucásicas, con un ICM entre 19-25 kg/cm2, y una edad entre 19 y 34 años. Sólo se han escogido aquellas muestras cuyo dataje histológico aplicando los criterios de Noyes, coincide entre los dos patólogos y con el día del ciclo menstrual.
Los resultados esperados para Ia predicción según los métodos de dataje actuales (histología y día de Ia última menstruación) serían que las ocho primeras muestras fueran LH+7 y las ocho restantes fueran clasificadas como OTRO. El análisis del perfil de expresión génica de las 16 muestras con los modelos de entrenamiento creados, mostró que las predicciones generadas eran distintas según el modelo.
Criterios para determinar Ia predicción
En función del error mínimo de clasificación del método con el conjunto de entrenamiento, se asigna un valor a Io que predice en el conjunto de prueba. Si el error es 0 para todos los modelos, el valor de Ia predicción es el mismo, pero si no, se calcula Ia contribución relativa para darle a Ia predicción mayor o menor peso. Existe una relación inversamente proporcional entre el error y Io que debería contribuir a Ia decisión. Si se relativizan los errores a uno, y se restan a uno ese valor que da cada método, se consigue hacer Ia relación directa, ya que a mayor error, Ia resta con respecto de uno, da un valor menor, y con ello, los métodos con más error tienen un coeficiente de ponderación menor, y Ia aportación a Ia decisión de predecir, LH+7 u OTRO será menor.
A LH+7 se Ie da el valor de 1, y a OTRO el valor de -1, se multiplica el coeficiente de ponderación según el error del método, por 1 o -1 según haya predicho. Al final se comparan todos los métodos para cada array, y se hace el sumatorio, si el valor que da es >0, Io clasifica como LH+7, y si el valor es <0, como OTRO.
En función de todo esto se establece una predicción consenso (ver figura 8).
La predicción establece un valor distinto al clasificado según Noyes en las muestras 7 (LH+7), 14 (día 26), y 15 (día 16), pero esto no es considerado un error, ya que estas muestras aunque correspondan a un día del ciclo en el primer caso receptivo y en el segundo y tercero no receptivo, pueden tener un perfil génico de expresión distinto, tal y como predice el modelo según el perfil de expresión.
En Ia figura 9 se muestra Ia comparación de los resultados con las predicciones obtenidas con el ERA contra los obtenidos por datación morfológica (Noyes).
La datación de Noyes posee mayor tasa de error en los días que comprenden
Ia fase lútea del ciclo ovárico, y precisamente las muestras 14 y 15, que son de otro día del ciclo que se clasifican distinto son las dos que se encuentran en este momento del ciclo. En cuanto a Ia muestra 7 del array, es una paciente LH+7 según Noyes, y
OTRO día del ciclo según el ERA.
El ERA se puede utilizar para Ia identificación positiva de Ia receptividad endometrial así como para el diagnóstico de su alteración asociado a alteraciones endometriales propias de patologías como Ia endometriosis, fallo de implantación, hidrosalpinx, etc. Así mismo esta herramienta diagnóstica permitiría detectar las modificaciones funcionales inducidas por fármacos interceptivos o que pretendan mejorar Ia receptividad endometrial, alterando Ia situación de normalidad/anormalidad en el perfil receptivo del endometrio de una mujer. BIBLIOGRAFÍA
- Al-Shahrour F y Dopazo J. In Azuaje F y Dopazo J (eds), Data analysis and visualization in genomics and proteomics. Wiley 2005; 99-112.
- Al-Shahrour F1 Minguez P, Vaquerizas JM, Conde L y Dopazo J. BABELOMICS: a suite of web tools for functional annotation and analysis of groups of genes ¡n high- throughput experiments. Nucleic Acids Res 2005; 33:460-464,
- Balasch J, Fabregues F, Creus M y Vanrell JA. The usefulness of endometrial biopsy for luteal phase evaluation in infertility. Hum Reprod 1992; 7:973-977.
- Batista MC, Cartledge TP, Merino MJ, Axiotis C, Platia MP, Merriam GR. Midluteal phase endometrial biopsy does not accurately predict luteal function. Fértil Steril 1993;
59:294-300,
-Borthwick J, Charnock-Jones S, Tom BD y cois, (2003) Determination of the transcript profile of human endometrium. Mol Hum Reprod 9, 19-33,
-Carson D, Lagow E, Thathiah A y cois, (2002) Changes in gene expression during the early to mid-luteal (receptive phase) transition in human endometrium detected by high-density microarray screening. Mol Hum Reprod 8, 971-979.
- Catalano RD, Yanaihara A, Evans AL, Rocha D, Prentice A, Saidi S, Print CG, Charnock-Jones DS, Sharkey AM and Smith SK (2003) The effect of RU486 on the gene expression profile in an endometrial explant model Mol Human Reprod 9,465- 473.
- Coutifaris C, Myers ER, Guzick DS, Diamond MP, Carson SA, Legro RS, McGovern PG, Schlaff WD, Carr BR, Steinkampf MP, Silva S, Vogel DL y Leppert PC. Histological dating of timed endometrial biopsy tissue is not related to fertility status. Fértil Steril 2004; 82:1264-72.
- Creus M, Ordi J, Fabregues F, Casamitjana R, Ferrer B, CoII E, Vanrell JA y Balasch J. Alphavbeta 3 integrin expression and pinopod formation in normal and out-of-phase endometria of fertile and ¡nfertile women. Hum Reprod 2002; 17:2279-2286, -Horcajadas JA, Sharkey AM, Catalano RD, Sherwin JRA, Domínguez F, Burgos LA, Castro A, Peraza MR, Pellicer A and Simón C (2006) Use of Gene-Expression Profiling to Identify Human Endometrial Refractoriness. J Clin Endocrinol Metabol.
-Horcajadas JA, Pellicer A y Simón C. (2007) Wide Genomic Analysis of Human Endometrial Receptivity. New times, new opportunities. Human Reprod Update 13, 77- 86,
-Kliman HJ, Honig S, WaIIs D, Luna M, McSweet JC, Copperman AB. Optimization of endometrial preparation results in a normal endometrial function test (EFT) and good reproductive outcome in donor ovum recipients. J Assist Reprod Genet 2006; 23:299- 303,
- Lessey BA, Castelbaum AJ, Sawin SW, Sun J. Integrins as markers of uterine receptivity in women with primary unexplained ¡nfertility. Fértil Steril 1995; 63:535-542.
- Li TC, Dockery P, Rogers AW y Cooke ID. (How precise is histologic dating of endometrium using the standard dating criteria?. Fértil Steril 1989; 51:759-763,
- Medina I, Montaner D, Tarraga J, Dopazo J. Prophet, a web-based tool for class prediction using microarray data. Bioinformatics. 2007; 23(3):390-1.
-Mirkin S, Arslan M, Churikov D, Corica A, Diaz Jl, Williams S, Bocea S y Oehninger S (2005) In search of candidate genes critically expressed in the human endometrium during the window of implantation Human Reprod 20:2104-2117.
-Mirkin S, Nikas G, Hsiu JG, Diaz J and Oehninger S (2004) Gene expression profiles y structural/functional features of the peri-implantation endometrium in natural and gonadotropin-stimulated eyeles. J Clin Endocrinol Metab 89:5742-5752.
- Montaner D, Tarraga J, Huerta-Cepas J, Burguet J, Vaquerizas JM, Conde L, Minguez P, Vera J, Mukherjee S, VaIIs J, Pujana MAG, Alloza E, Herrero J, Al-
Shahrour F y Dopazo J. Next station ¡n microarray data analysis: GEPAS. Accepted Nucleic Acids Res. 2006, - Murray WIJ, Meyer WR, Zaino RJ, Lessey BA, Novotny DB, lreland K, Zeng D y Fritz MA. A critical analysis of the accuracy, reproducibility, and clinical utility of histologic endometrial dating in fertile women. Fértil Steril 2004; 81:1333-1343,
- Noyes RW, Hertig AT, y Rock J. Dating the endometrial biopsy. Fértil Steril 1950; 1 :3- 17.
- Ordi J, Creus M, Quinto L, Casamitjana R, Cardesa A y Balasch J. Within-subject between-cycle variability of histological dating, alpha v beta 3 integrin expression, and pinopod formation in the human endometrium. J Clin Endocrinol Metab 2003; 88:2119-2125,
- Papanikolaou EG, TouARNye H, Verpoest W, Camus M, VeARNeve V, Van Steirteghem A, Devroey P; http://www.ncbi.nlm.nih. gov/sites/entrez?Db=pubmed&Cmd=ShowDetailView&TermToSearc h=15576388&ordinalpos=79&itool=EntrezSystem2.PEntrez.Pubmed.Pubmed_ResultsPanel.Pu bmed RVDocSum Early and late ovarían hyperstimulation syndrome: early pregnancy outcome and profile. Hum Reprod. 2005; 20(3): 636-641.
- Ponnampalam AP, Weston GC, Trajstman AC. Molecular classification of human endometrial cycle stages by transcriptional profiling. Mol Hum Reprod 2004; 10, 879- 893,
-Riesewijk A, Martin J, Horcajadas JA Polman J, Pellicer A, Mosselman S y Simón C (2003) Gene expression profiling of human endometrial receptivity on days LH+2 contra LH+7 by microarray technology. Mol Hum Reprod 9:253-264,
- Schena Wl, Shalon D, Davis RW y Brown PO. Quantitative monitoring of gene expression pattems with a complementary ADN microarray. Science 1995; 270:467- 470,
- Sharkey AWI, Catalano R, Evans A, Charnock-Jones DS and Smith SK (2005) Novel antiangiogenic agents for use in contraception. Contraception 71,263-271.
- Shoupe D, Mishell DR Jr, Lacarra M, Lobo RA, Horenstein J, d'Ablaing G. Correlation of endometrial maturation with four methods of estimating day of ovulation. Obstet Gynecol. Obstet Gynecol 1989; 73:88-92. - Talbi S, Hamilton AE, Vo KC, Tulac S, Overgaard M T, Dosiou C, Le Shay N, Nezhat, CN, Kempson R, Lessey BA, Nayak NR y Giudice LC. Molecular phenotyping of human endometrium distinguishes menstrual cycle phases and underlying biological processes in normo-ovulatory women. Endocrinology 2005; 147:1097-1121.
-Wilcox AJ, Baird DD, Weinberg CR. Time of implantation of the conceptus and loss of pregnancy. N Engl J Med. 1999;340: 1796-1799.
T-REX (http://www.gepas.org/)
FATIGO(http://babelomics.b¡oínfo.cipf.es/EntrvPoint?loadForm=fatigo)
PROPHET(http://gepas.b¡o¡nfo.cipf.es/cgi-bin/loadtool.cgi?tool=prophet)
Agilent earrav 4.5 ( https://earrav.chem.agilent.com/earray/)

Claims

REINVINDICACIONES
1. Un método para detectar en una muestra biológica Ia situación de normalidad/anormalidad en el perfil receptivo de un endometrio, que comprende:
a) obtener una biopsia endometrial del fondo uterino de una mujer 7 días después de su pico de LH endógena, Io que equivale al día 20-21 del ciclo menstrual.
b) realizar Ia extracción de ARN de Ia biopsia endometrial;
c) mediante Ia tecnología de microarrays determinar en dicha muestra el perfil de expresión de los genes del ERA ;
d) detectar en dicha biopsia el perfil de expresión de los genes del endometrio que se contemplan en Ia tabla de Ia figura 1 ; y
e) analizar dicho perfil de expresión de los genes mediante los programas bioinformáticos específicos que clasifican y determinan el estado del endometrio en función del perfil génico con los criterios fijados.
2. Un método para detectar en una muestra biológica Ia situación de normalidad/anormalidad en el perfil receptivo de un endometrio, según Ia reivindicación 1 , en el que Ia muestra endometrial obtenida en (a) se ha puesto en contacto con un agente específico para Ia evaluación de cada gen.
3. Un método para detectar en una muestra biológica Ia situación de normalidad/anormalidad en el perfil receptivo de un endometrio, según las reivindicaciones 1-2, en el que el agente específico es un oligo (sonda) que es complementario a una región del gen del que se quiere cuantificar su expresión.
4. Un método para detectar en una muestra biológica Ia situación de normalidad/anormalidad en el perfil receptivo de un endometrio, según las reivindicaciones 1-3, en el que cada gen tiene al menos una sonda.
5. Un método para detectar en una muestra biológica Ia situación de normalidad/anormalidad en el perfil receptivo de un endometrio, según las reivindicaciones 1-4, en el que se ha extraído el ARNm a partir de Ia muestra endometrial y se ha determinado el perfil receptivo del endometrio por comparación con los genes y las sondas de Ia tabla de Ia figura 1.
6. Un método para detectar en una muestra biológica Ia situación de normalidad/anormalidad en el perfil receptivo de un endometrio, según las reivindicaciones 1-5, en el que dicho ARNm normalmente es inducido (o reprimido) en esa fase del ciclo en un endometrio receptivo y es reprimido (o inducido) en esa misma fase del ciclo femenino en un endometrio no receptivo.
7. Un método para detectar en una muestra biológica Ia situación de normalidad/anormalidad en el perfil receptivo de un endometrio, según las reivindicaciones 1-6, en el que el perfil de expresión se ajusta al establecido por el clasificador para el ERA.
8. Un método para detectar en una muestra biológica Ia situación de anormalidad en el perfil receptivo de un endometrio, según las reivindicaciones 1-7, producida por situaciones de subfertilidad como Ia endometriosis o el hidrosálpinx.
9. Un método para detectar en una muestra biológica Ia situación de anormalidad en el perfil receptivo de un endometrio posiblemente producida por endometriosis, según Ia reivindicación 8, en el que al menos uno de los marcadores relacionados con Ia endometriosis se selecciona entre el grupo siguiente: un gen de Ia tabla de Ia figura 1 , un gen y al menos una de sus sondas de Ia tabla de Ia figura 1 , ARN que da lugar a ADNc derivado de dicho gen, fragmentos de dicho ARN, y proteínas o péptidos codificados por dicho gen o por el ARN o uno de sus fragmentos que produce el ADNc derivado de dicho gen.
10. Un método para detectar en una muestra biológica Ia situación de anormalidad en el perfil receptivo de un endometrio posiblemente producida por endometriosis, según Ia reivindicación 9, en el que el marcador de Ia endometriosis es cualquier molécula codificada por ese ARNm.
11. Un método de acuerdo con las reivindicaciones 1-7, para detectar en una muestra biológica el efecto de fármacos o dispositivos inertes o en combinación con fármacos que alteren Ia situación de normalidad/anormalidad en el perfil receptivo de un endometrio, según el perfil de expresión determinado por el ERA, cuyos resultados serían de un cambio genómico compatible con Ia no receptividad.
12. Un método de acuerdo con las reivindicación 11 , para detectar en una muestra biológica el efecto de fármacos que alteren Ia situación de normalidad/anormalidad en el perfil receptivo de un endometrio, según el perfil de expresión determinado por el ERA, tanto en una mujer sana como con diversos estados patológicos, como por ejemplo, endometriosis, fallo de implantación, etc
PCT/ES2008/000513 2008-07-22 2008-07-22 Perfil de expresion genetica como marcador de la receptividad endometrial Ceased WO2010010201A1 (es)

Priority Applications (8)

Application Number Priority Date Filing Date Title
PCT/ES2008/000513 WO2010010201A1 (es) 2008-07-22 2008-07-22 Perfil de expresion genetica como marcador de la receptividad endometrial
EP09800092.0A EP2333107B1 (en) 2008-07-22 2009-07-22 Gene expression profile as an endometrial receptivity marker
US13/057,135 US10081840B2 (en) 2008-07-22 2009-07-22 Gene expression profile as an endometrial receptivity marker
PL09800092T PL2333107T3 (pl) 2008-07-22 2009-07-22 profil ekspresji genów jako marker wrażliwości śluzówki macicy
ES09800092.0T ES2484417T3 (es) 2008-07-22 2009-07-22 Perfil de expresión génica como marcador de la receptividad endometrial
PCT/ES2009/000386 WO2010010213A1 (es) 2008-07-22 2009-07-22 Perfil de expresión génica como marcador de la receptividad endometrial
DK09800092.0T DK2333107T3 (da) 2008-07-22 2009-07-22 Genekspressionsprofil som markør for endometriel receptivitet
CA2732849A CA2732849C (en) 2008-07-22 2009-07-22 Gene expression profile as an endometrial receptivity marker

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/ES2008/000513 WO2010010201A1 (es) 2008-07-22 2008-07-22 Perfil de expresion genetica como marcador de la receptividad endometrial

Publications (1)

Publication Number Publication Date
WO2010010201A1 true WO2010010201A1 (es) 2010-01-28

Family

ID=41570044

Family Applications (2)

Application Number Title Priority Date Filing Date
PCT/ES2008/000513 Ceased WO2010010201A1 (es) 2008-07-22 2008-07-22 Perfil de expresion genetica como marcador de la receptividad endometrial
PCT/ES2009/000386 Ceased WO2010010213A1 (es) 2008-07-22 2009-07-22 Perfil de expresión génica como marcador de la receptividad endometrial

Family Applications After (1)

Application Number Title Priority Date Filing Date
PCT/ES2009/000386 Ceased WO2010010213A1 (es) 2008-07-22 2009-07-22 Perfil de expresión génica como marcador de la receptividad endometrial

Country Status (7)

Country Link
US (1) US10081840B2 (es)
EP (1) EP2333107B1 (es)
CA (1) CA2732849C (es)
DK (1) DK2333107T3 (es)
ES (1) ES2484417T3 (es)
PL (1) PL2333107T3 (es)
WO (2) WO2010010201A1 (es)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8323177B2 (en) 2009-08-22 2012-12-04 The Board Of Trustees Of The Leland Stanford Junior University Imaging and evaluating embryos, oocytes, and stem cells
EP3037101A1 (en) 2014-12-22 2016-06-29 Ferring B.V. Oxytocin receptor antagonist therapy in the luteal phase for implantation and pregnancy in women undergoing assisted reproductive technologies
US9482659B2 (en) 2010-09-27 2016-11-01 Progyny, Inc. Apparatus, method, and system for the automated imaging and evaluation of embryos, oocytes and stem cells
US9879307B2 (en) 2011-02-23 2018-01-30 The Board Of Trustees Of The Leland Stanford Junior University Methods of detecting aneuploidy in human embryos
US10241108B2 (en) 2013-02-01 2019-03-26 Ares Trading S.A. Abnormal syngamy phenotypes observed with time lapse imaging for early identification of embryos with lower development potential

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4953334B2 (ja) 2009-08-24 2012-06-13 国立大学法人金沢大学 遺伝子発現プロファイルによる消化器癌、胃癌、大腸癌、膵臓癌及び胆道癌の検出
EP2687851A1 (en) 2012-07-20 2014-01-22 Matricelab Innove Method for increasing implantation success in assisted fertilization
HRP20191240T1 (hr) 2014-06-17 2019-11-01 Igenomix S L Terapija matičnim stanicama kod patologija endometrija
US10918327B2 (en) 2017-02-02 2021-02-16 Coopersurgical, Inc. Compositions and methods for determining receptivity of an endometrium for embryonic implantation
EP3569718A1 (en) * 2018-05-16 2019-11-20 Integrated Genetic Lab Services SLU Kit and method for determining the receptivity status of an endometrium
CN110042156B (zh) * 2019-04-22 2021-12-28 苏州亿康医学检验有限公司 一种判断子宫内膜容受性的方法及其应用
US20220196670A1 (en) * 2019-05-17 2022-06-23 Tervisetehnoloogiate Arenduskeskus As Endometrial receptivity determination
CN111778326B (zh) * 2020-07-14 2021-10-22 和卓生物科技（上海）有限公司 用于子宫内膜容受性评估的基因标志物组合及其应用
CN120888649A (zh) * 2022-03-15 2025-11-04 苏州亿康医学检验有限公司 无创方式判断子宫内膜容受性的方法、模型和标志物
EP4311862A1 (en) 2022-07-29 2024-01-31 Ivi Rma Global, Sl. Methods for detection of embryo implantation failure of endometrial origen
GB202216420D0 (en) 2022-11-04 2022-12-21 Univ Warwick Biomarkers of pregnancy loss

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040005612A1 (en) * 2002-05-14 2004-01-08 Giudice Linda C. Endometrial genes in endometrial disorders

Family Cites Families (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0473612B1 (en) 1989-05-05 1996-07-31 Northern Sydney Area Health Service Fertility enhancement
US6979533B2 (en) 1992-06-12 2005-12-27 The Trustees Of The University Of Pennsylvania Determination of endometrial receptivity toward embryo implantation
CA2203718A1 (en) 1997-04-25 1998-10-25 University Of British Columbia Cadherin-11 as an indicator of viable pregnancy
US6821724B1 (en) 1998-09-17 2004-11-23 Affymetrix, Inc. Methods of genetic analysis using nucleic acid arrays
IL127896A0 (en) 1998-12-31 1999-10-28 Gail Laster A novel treatment for preeclampsia and related diseases
CA2400065A1 (en) 2000-03-08 2001-09-13 Harvey J. Kliman Methods of diagnosing and monitoring endometrial glandular development
JP2004504275A (ja) 2000-05-24 2004-02-12 シエーリング　アクチエンゲゼルシャフト フィブリン−１の医薬用途
US6780594B2 (en) 2000-09-25 2004-08-24 Schering Aktiengesellschaft Method for in vitro diagnosis of endometriosis
WO2005016230A2 (en) 2002-06-07 2005-02-24 President And Fellows Of Harvard College Evaluating protein signatures
AUPR673001A0 (en) 2001-07-31 2001-08-23 Prince Henry's Institute Of Medical Research Pregnancy-related enzyme activity
DE10139874A1 (de) 2001-08-10 2003-04-30 Schering Ag Humane Mater-Proteine
WO2005040094A1 (en) 2003-10-24 2005-05-06 Postech Foundation Novel dendrimer compound and a biochip using the same
GB0201284D0 (en) 2002-01-21 2002-03-06 Isis Innovation Screening methods for contraceptive and fertility agents
US20030186300A1 (en) 2002-03-20 2003-10-02 Ali Akoum Methods and products for modulation of reproductive processes and for diagnosis, prognostication and treatment of related conditions
JP2005539016A (ja) 2002-08-02 2005-12-22 シエーリング アクチエンゲゼルシャフト 女性用受胎調節及びホルモン置換療法のための高められた抗ゴナドトロピン活性を有するプロゲステロン受容体モジュレーター
DE10236405A1 (de) 2002-08-02 2004-02-19 Schering Ag Progesteronrezeptormodulatoren mit erhöhter antigonadotroper Aktivität für die weibliche Fertilitätskontrolle und Hormonersatztherapie
US20040121339A1 (en) 2002-12-19 2004-06-24 Jizhong Zhou Special film-coated substrate for bio-microarray fabrication and use thereof
WO2004058999A2 (de) 2002-12-21 2004-07-15 Universität Leipzig Verfahren und mittel zur bestimmung von bestimmten zuständen bzw. veränderungen im uterusepithel und im epithel anderer organe
JP4007269B2 (ja) 2003-07-29 2007-11-14 オムロン株式会社 生体分子の形状転写方法、チップ基板の製造方法及びバイオチップの製造方法
US20050042363A1 (en) 2003-08-18 2005-02-24 Kukhtin Alexander V. Method for fabrication of biochips with a macroporous polymer substrate
WO2005018796A1 (en) 2003-08-21 2005-03-03 Pamgene B.V. Microarray support for bioprobe synthesis
EP1670902A4 (en) 2003-09-08 2008-05-07 Univ Texas METHOD AND COMPOSITION FOR ENHANCING i IN VITRO /i EMBRYO DEVELOPMENT BY SUPPLEMENTING CULTURE MEDIUM WITH PROSTAGLANDIN OR A PROSTAGLANDIN ANALOG
JP2007504807A (ja) 2003-09-08 2007-03-08 ボード オブ リージェンツ オブ ザ ユニバーシティー オブ テキサス システム 培養培地にプロスタグランジン又はプロスタグランジン類似物を補足することによってインビトロ胚発生を高めるための方法及び組成
EP1711618B1 (en) 2003-12-23 2010-08-04 Mount Sinai Hospital Corporation Methods for detecting markers associated with endometrial disease or phase
EP1765083A4 (en) 2004-05-28 2008-08-20 Serono Lab USE OF IL-17 TO TREAT FERTILITY DISORDERS
JP2007278750A (ja) 2006-04-04 2007-10-25 Ono Pharmaceut Co Ltd プロスタグランジン類及びその類縁体の定量方法
WO2009143633A1 (en) 2008-05-30 2009-12-03 Institut De Recherches Cliniques De Montreal Pcsk9 inhibitors and methods of use thereof
EP2348318A1 (en) 2010-01-21 2011-07-27 Equipo Ivi Investigación, S.L. Diagnostic method for endometrial receptivity
WO2013057316A2 (en) 2011-10-21 2013-04-25 INSERM (Institut National de la Santé et de la Recherche Médicale) Methods for assessing endometrial receptivity of a patient after controlled ovarian hyperstimulation
US9423402B2 (en) 2011-11-30 2016-08-23 Singapore Health Services Pte. Ltd. GM1 ganglioside to annexin V microparticle polypeptide ratio for biological monitoring
EP2867670B1 (en) 2012-06-27 2018-05-30 Roche Diagnostics GmbH MEANS AND METHODS APPLYING sFlt-1/PlGF OR ENDOGLIN/PlGF RATIO TO RULE-OUT ONSET OF PREECLAMPSIA WITHIN A CERTAIN TIME PERIOD
CN106662589B (zh) 2014-03-21 2019-07-30 艾基诺米公司 先兆子痫的早期检测
HRP20191240T1 (hr) 2014-06-17 2019-11-01 Igenomix S L Terapija matičnim stanicama kod patologija endometrija

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040005612A1 (en) * 2002-05-14 2004-01-08 Giudice Linda C. Endometrial genes in endometrial disorders

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
BORTHWICK, J. M. ET AL.: "Determination of the transcript profile of human endometrium.", MOLECULAR HUMAN REPRODUCTION., vol. 9, no. 1, January 2003 (2003-01-01), pages 19 - 33 *
BURNEY, R. O. ET AL.: "Gene expression analysis of endometrium reveals progesterone resistance and candidate susceptibility genes in women with endometriosis.", ENDOCRINOLOGY., vol. 148, no. 8, August 2007 (2007-08-01), pages 3814 - 3826 *
CARSON, D. D. ET AL.: "Changes in gene expression during the early to mid-luteal (receptive phase) transition in human endometrium detected by high-density microarray screening.", MOLECULAR HUMAN REPRODUCTION., vol. 8, no. 9, September 2002 (2002-09-01), pages 871 - 879 *
HORCAJADAS, J. A. ET AL.: "Effect of an intrauterine device on the gene expression profile of the endometrium.", THE JOURNAL OF CLINICAL ENDOCRINOLOGY & METABOLISM., vol. 91, no. 8, August 2006 (2006-08-01), pages 3199 - 3207 *
HORCAJADAS, J. A. ET AL.: "Effect of controlled ovarian hyperstimulation in IVF on endometrial gene expression profiles.", MOLECULAR HUMAN REPRODUCTION., vol. 11, no. 3, March 2005 (2005-03-01), pages 195 - 205 *
RIESEWIJK, A. ET AL.: "Gene expression profiling of human endometrial receptivity on days LH+2 versus LH+7 by microarray technology.", MOLECULAR HUMAN REPRODUCTION., vol. 9, no. 5, May 2003 (2003-05-01), pages 253 - 264 *
TALBI, S. ET AL.: "Molecular phenotyping of human endometrium distinguishes menstrual cycle phases and underlying biological processes in normo- ovulatory women.", ENDOCRINOLOGY., vol. 147, no. 3, March 2006 (2006-03-01), pages 1097 - 1121 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8337387B2 (en) 2009-08-22 2012-12-25 The Board Of Trustees Of The Leland Stanford Junior University Imaging and evaluating embryos, oocytes, and stem cells
US8721521B2 (en) 2009-08-22 2014-05-13 The Board Of Trustees Of The Leland Stanford Junior University Imaging and evaluating embryos, oocytes, and stem cells
US8951184B2 (en) 2009-08-22 2015-02-10 The Board Of Trustees Of The Leland Stanford Junior University Imaging and evaluating embryos, oocytes, and stem cells
US8989475B2 (en) 2009-08-22 2015-03-24 The Board Of Trustees Of The Leland Stanford Junior University Imaging and evaluating embryos, oocytes, and stem cells
US9228931B2 (en) 2009-08-22 2016-01-05 The Board Of Trustees Of The Leland Stanford Junior University Imaging and evaluating embryos, oocytes, and stem cells
US8323177B2 (en) 2009-08-22 2012-12-04 The Board Of Trustees Of The Leland Stanford Junior University Imaging and evaluating embryos, oocytes, and stem cells
US9482659B2 (en) 2010-09-27 2016-11-01 Progyny, Inc. Apparatus, method, and system for the automated imaging and evaluation of embryos, oocytes and stem cells
US9879307B2 (en) 2011-02-23 2018-01-30 The Board Of Trustees Of The Leland Stanford Junior University Methods of detecting aneuploidy in human embryos
US10241108B2 (en) 2013-02-01 2019-03-26 Ares Trading S.A. Abnormal syngamy phenotypes observed with time lapse imaging for early identification of embryos with lower development potential
WO2016105190A1 (en) 2014-12-22 2016-06-30 Ferring B.V. Oxytocin receptor antagonist therapy in the luteal phase for implantation and pregnancy in women undergoing assisted reproductive technologies
EP3037101A1 (en) 2014-12-22 2016-06-29 Ferring B.V. Oxytocin receptor antagonist therapy in the luteal phase for implantation and pregnancy in women undergoing assisted reproductive technologies
EP3501533A1 (en) 2014-12-22 2019-06-26 Ferring B.V. Oxytocin receptor antagonist therapy in the luteal phase for implantation and pregnancy in women undergoing assisted reproductive technologies
US11752157B2 (en) 2014-12-22 2023-09-12 Ferring B.V. Oxytocin receptor antagonist therapy in the luteal phase for implantation and pregnancy in women undergoing assisted reproductive technologies

Also Published As

Publication number Publication date
US20120040849A1 (en) 2012-02-16
EP2333107A4 (en) 2012-03-14
CA2732849C (en) 2022-06-21
US10081840B2 (en) 2018-09-25
CA2732849A1 (en) 2010-01-28
ES2484417T3 (es) 2014-08-11
WO2010010213A1 (es) 2010-01-28
EP2333107B1 (en) 2014-04-23
DK2333107T3 (da) 2014-07-21
EP2333107A1 (en) 2011-06-15
PL2333107T3 (pl) 2014-11-28

Similar Documents

Publication Publication Date Title
WO2010010201A1 (es) Perfil de expresion genetica como marcador de la receptividad endometrial
US8202968B2 (en) Predicting lung cancer survival using gene expression
ES2650610T3 (es) Perfiles de expresión génica para predecir desenlaces en cáncer de mama
Herington et al. Gene profiling the window of implantation: Microarray analyses from human and rodent models
US11667974B2 (en) Diagnostic, prognostic and therapeutic uses of long noncoding RNAs for pathologies and toxicities inducing heart disorders
KR20170063519A (ko) 불임 및 관련 병리상태를 평가하기 위한 방법 및 시스템
JP2011500071A (ja) 遺伝子に基づくアルゴリズム的ガン予後及び患者の臨床結果
US20120094845A1 (en) Methods and devices for assessing infertility and/or egg quality
US20190241955A1 (en) Methods of determining endometrial receptivity and uses thereof
US20190276890A1 (en) Gene expression profile as an endometrial receptivity marker
WO2015079060A2 (en) Mirnas as advanced diagnostic tool in patients with cardiovascular disease, in particular acute myocardial infarction (ami)
JP2013510575A (ja) 卵丘細胞が差示発現する遺伝子およびそれらを用いる妊娠受容能をもつ卵母細胞の同定のためのアッセイ
TW201833335A (zh) 用於評估子宮內膜癌的發展風險或診斷子宮內膜癌的方法及套組
CN113755571B (zh) 用于胚胎着床成功率检测的生物标志物及应用
CN110129444B (zh) 一种检测放射性损伤的试剂及其应用
CN114941025B (zh) 用于诊断子痫前期的miRNA及其应用
WO2005030959A1 (ja) 神経芽細胞腫予後診断のためのマイクロアレイと神経芽細胞腫予後診断方法
Class et al. Patent application title: GENE EXPRESSION PROFILE AS AN ENDOMETRIAL RECEPTIVITY MARKER Inventors: Carlos Simón Valles (Valencia, ES) José Antonio Horcajadas Almansa (Valencia, ES) José Antonio Horcajadas Almansa (Valencia, ES) Patricia Diaz Gimeno (Valencia, ES) Antonio Pellicer Martínez (Valencia, ES) Assignees: EQUIPO IVI INVESTIGACION SL
JP2003235557A (ja) 精神分裂病により発現量が変化する遺伝子を規定する核酸を解析する方法
JP2010502179A (ja) 予後診断方法
JP6623548B2 (ja) 膵管内乳頭粘液性腫瘍検出用マーカーとキット
EP4311862A1 (en) Methods for detection of embryo implantation failure of endometrial origen
EP2840149B1 (en) Method of diagnosing down&#39;s syndrome
CN120418454A (zh) 用于诊断女性生殖道疾病或病症的生物标志物
WO2019180296A1 (es) Método de predicción de respuesta a tratamientos con cirugía bariátrica en diabetes mellitus tipo 2

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 08805342

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 08805342

Country of ref document: EP

Kind code of ref document: A1