SYSTEM FOR DETERMINING A HEALTH STATUS OF A TISSUE OF INTEREST
The invention generally relates to methods for assessing the health of a tissue by characterizing circulating nucleic acids in a biological sample. According to certain embodiments, methods for assessing the health of a tissue include the steps of detecting a sample level of RNA in a biological sample, comparing the sample level of RNA to a reference level of RNA specific to the tissue, determining whether a difference exists between the sample level and the reference level, and characterizing the tissue as abnormal if a difference is detected.
This application claims benefit and is a Continuation of application Ser. No. 15/969,603 filed May 2, 2018, which is a continuation of application Ser. No. 14/861,650 filed Sep. 22, 2015, now U.S. Pat. No. 10,240,200, issued Mar. 26, 2019, which is a Divisional of application Ser. No. 13/752,131 filed Jan. 28, 2013, now abandon, which claims benefit of U.S. Provisional Patent Application No. 61/591,642, filed on Jan. 27, 2012, now expired, which applications are incorporated herein by reference in their entirety.
TECHNICAL FIELDThe present invention relates the field of nucleic acid analysis from a biological sample containing genetic material. Specifically, methods of the invention relate to quantitating tissue-specific nucleic acids in a biological sample.
BACKGROUNDIt is often challenging to gauge the health of organs within an individual's body. Physicians are often forced to use expensive imaging techniques or perform invasive biopsies for cancer screening to identify diagnostic biomarkers and monitor tumor initiation and progression. The invasive nature of biopsies makes them unsuitable for widespread screening of patients. In addition, many diagnostic biomarkers are only identified in cancer cell lines or from biopsy specimens obtained from patients with late-stage disease and metastasis.
The presence of circulating nucleic acids (DNA and RNA) detectable in the plasma and serum of cancer patients has been investigated for its potential use to serve as markers for diagnostic purposes, with the obvious benefit being a non-invasive diagnostic tool. It has been shown that markers within the plasma are identical to the ones found in the carcinogenic tissue of the patient. Circulating RNA is particularly of interest for use in early detection cancer screenings due to RNA markers close association with malignancy.
In addition to cancer detection, the discovery of fetal specific cell-free RNA present in maternal plasma has opened up new horizons on prenatal molecular diagnostics (see e.g., Poon et al., Clinical Chemistry, 46(11): 1832-1834 (2000)). Specifically, analysis of plasma RNA holds promise for noninvasive gene expression profiling of the fetus. However, only a handful of pregnancy specific cell-free RNA transcripts have been characterized to date. A comprehensive profiling of such RNA has not been performed.
A problem with analyzing cell-free RNA in non-maternal and maternal blood is the lack of suitable data to estimate the biological causes of the cell-free RNA present. For example, there lacks a reliable method for determining tissue origins of the cell-free RNA present in blood.
SUMMARYThe present invention provides methods for profiling the origin of the cell-free RNA to assess the health of an organ or tissue. Deviations in normal cell-free transcriptomes are caused when organ/tissue-specific transcripts are released in to the blood in large amounts as those organs/tissue begin to fail or are attacked by the immune system or pathogens. As a result inflammation process can occur as part of body's complex biological response to these harmful stimuli. The invention, according to certain aspects, utilizes tissue-specific RNA transcripts of healthy individuals to deduce the relative optimal contributions of different tissues in the normal cell-free transcriptome, with each tissue-specific RNA transcript of the sample being indicative of the apotopic rate of that tissue. The normal cell-free transcriptome serves as a baseline or reference level to assess tissue health of other individuals. The invention includes a comparative measurement of the cell-free transcriptome of a sample to the normal cell free transcriptome to assess the sample levels of tissue-specific transcripts circulating in plasma and to assess the health of tissues contributing to the cell-free transcriptome.
In addition to normal reference levels, methods of the invention also utilize reference levels for cell-free transcriptomes specific to other patient populations. Using methods of the invention one can determine the relative contribution of tissue-specific transcripts to the cell-free transcriptome of maternal subjects, fetus subjects, and/or subjects having a condition or disease.
By analyzing the health of tissue based on tissue-specific transcripts, methods of the invention advantageously allow one to assess the health of a tissue without relying on disease-related protein biomarkers. In certain aspects, methods of the invention assess the health of a tissue by comparing a sample level of RNA in a biological sample to a reference level of RNA specific to a tissue, determining whether a difference exists between the sample level and the reference level, and characterizing the tissue as abnormal if a difference is detected. For example, if a patient's RNA expression levels for a specific tissue differs from the RNA expression levels for the specific tissue in the normal cell-free transcriptome, this indicates that patient's tissue is not functioning properly.
In certain aspects, methods of the invention involve assessing health of a tissue by characterizing the tissue as abnormal if a specified level of RNA is present in the blood. The method may further include detecting a level of RNA in a blood sample, comparing the sample level of RNA to a reference level of RNA specific to a tissue, determining whether a difference exists between the sample level and the reference level, and characterizing the tissue as abnormal if the sample level and the reference level are the same.
The present invention also provides methods for comprehensively profiling fetal specific cell-free RNA in maternal plasma and deconvoluting the cell-free transcriptome of fetal origin with relative proportion to different fetal tissue types. Methods of the invention involve the use of next-generation sequencing technology and/or microarrays to characterize the cell-free RNA transcripts that are present in maternal plasma at different stages of pregnancy. Quantification of these transcripts allows one to deduce changes of these genes across different trimesters, and hence provides a way of quantification of temporal changes in transcripts.
Methods of the invention allow diagnosis and identification of the potential for complications during or after pregnancy. Methods also allow the identification of pregnancy-associated transcripts which, in turn, elucidates maternal and fetal developmental programs. Methods of the invention are useful for preterm diagnosis as well as elucidation of transcript profiles associated with fetal developmental pathways generally. Thus, methods of the invention are useful to characterize fetal development and are not limited to characterization only of disease states or complications associated with pregnancy. Exemplary embodiments of the methods are described in the detailed description, claims, and figures provided below.
The invention is best understood from the following detailed description when read in conjunction with the accompanying drawings. It is emphasized that, according to common practice, the various features of the drawings are not to-scale. On the contrary, the dimensions of the various features are arbitrarily expanded or reduced for clarity. Included in the drawings are the following figures.
Methods and materials described herein apply a combination of next-generation sequencing and microarray techniques for detecting, quantitating and characterizing RNA sequences present in a biological sample. In certain embodiments, the biological sample contains a mixture of genetic material from different genomic sources, i.e. pregnant female and a fetus.
Unlike other methods of digital analysis in which the nucleic acid in the sample is isolated to a nominal single target molecule in a small reaction volume, methods of the present invention are conducted without diluting or distributing the genetic material in the sample. Methods of the invention allow for simultaneous screening of multiple transcriptomes, and provide informative sequence information for each transcript at the single-nucleotide level, thus providing the capability for non-invasive, high throughput screening for a broad spectrum of diseases or conditions in a subject from a limited amount of biological sample.
In one particular embodiment, methods of the invention involve analysis of mixed fetal and maternal RNA in the maternal blood to identify differentially expressed transcripts throughout different stages of pregnancy that may be indicative of a preterm or pathological pregnancy. Differential detection of transcripts is achieved, in part, by isolating and amplifying plasma RNA from the maternal blood throughout the different stages of pregnancy, and quantitating and characterizing the isolated transcripts via microarray and RNA-Seq.
Methods and materials specific for analyzing a biological sample containing RNA (including non-maternal, maternal, maternal-fetus mixed) as described herein, are merely one example of how methods of the invention can be applied and are not intended to limit the invention. Methods of the invention are also useful to screen for the differential expression of target genes related to cancer diagnosis, progression and/or prognosis using cell-free RNA in blood, stool, sputum, urine, transvaginal fluid, breast nipple aspirate, cerebrospinal fluid, etc.
In certain embodiments, methods of the invention generally include the following steps: obtaining a biological sample containing genetic material from different genomic sources, isolating total RNA from the biological sample containing biological sample containing a mizture of genetic material from different genomic sources, preparing amplified cDNA from total RNA, sequencing amplified cDNA, and digital counting and analysis, and profiling the amplified cDNA.
Methods of the invention also involve assessing the health of a tissue contributing to the cell-free transcriptome. In certain embodiments, the invention involves assessing the cell-free transcriptome of a biological sample to determine tissue-specific contributions of individual tissues to the cell-free transcriptome. According to certain aspects, the invention assesses the health of a tissue by detecting a sample level of RNA in a biological sample, comparing the sample level of RNA to a reference level of RNA specific to the tissue, and characterizing the tissue as abnormal if a difference is detected. This method is applicable to characterize the health of a tissue in non-maternal subjects, pregnant subjects, and live fetuses.
In certain aspects, methods of the invention employ a deconvolution of a reference cell-free RNA transcriptome to determine a reference level for a tissue. Preferably, the reference cell-free RNA transcriptome is a normal, healthy transcriptome, and the reference level of a tissue is a relative level of RNA specific to the tissue present in the blood of healthy, normal individuals. Methods of the invention assume that apoptotic cells from different tissue types release their RNA into plasma of a subject. Each of these tissues expresses a specific number of genes unique to the tissue type, and the cell-free RNA transcriptome of a subject is a summation of the different tissue types. Each tissue may express one or more numbers of genes. In certain embodiments, the reference level is a level associated with one of the genes expressed by a certain tissue. In other embodiments, the reference level is a level associated with a plurality of genes expressed by a certain tissue. It should be noted that a reference level or threshold amount for a tissue-specific transcript present in circulating RNA may be zero or a positive number.
For healthy, normal subjects, the relative contributions of circulating RNA from different tissue types are relatively stable, and each tissue-specific RNA transcript of the cell-free RNA transcriptome for normal subjects can serve as a reference level for that tissue. Applying methods of the invention, a tissue is characterized as unhealthy or abnormal if a sample includes a level of RNA that differs from a reference level of RNA specific to the tissue. The tissue of the sample may be characterized as unhealthy if the actual level of RNA is statistically different from the reference level. Statistical significance can be determined by any method known in the art. These measurements can be used to screen for organ health, as diagnostic tool, and as a tool to measure response to pharmaceuticals or in clinical trials to monitor health.
If a difference is detected between the sample level of RNA and the reference level of RNA, such difference suggests that the associated tissue is not functioning properly. The change in circulating RNA may be the precursor to organ failure or indicate that the tissue is being attacked by the immune system or pathogens. If a tissue is identified as abnormal, the next step(s), according to certain embodiments, may include more extensive testing of the tissue (e.g. invasive biopsy of the tissue), prescribing course of treatment specific to the tissue, and/or routine monitoring of the tissue.
Methods of the invention can be used to infer organ health non-invasively. This non-invasive testing can be used to screen for appendicitis, incipient diabetes and pathological conditions induced by diabetes such as nephropathy, neuropathy, retinopathy etc. In addition, the invention can be used to determine the presence of graft versus host disease in organ transplants, particularly in bone marrow transplant recipients whose new immune system is attacking the skin, GI tract or liver. The invention can also be used to monitor the health of solid organ transplant recipients such as heart, lung and kidney. The methods of the invention can assess likelihood of prematurity, preeclampsia and anomalies in pregnancy and fetal development. In addition, methods of the invention could be used to identify and monitor neurological disorders (e.g. multiple sclerosis and Alzheimer's disease) that involve cell specific death (e.g. of neurons or due to demyelination) or that involve the generation of plaques or protein aggregation.
A cell-free transcriptome for purposes of determining a reference level for tissue-specific transcripts can be the cell-free transcriptome of one or more normal subjects, maternal subjects, subjects having a certain conditions and diseases, or fetus subjects. In the case of certain conditions, the reference level of a tissue is a level of RNA specific to the tissue present in blood of one or more subjects having a certain disease or condition. In such aspect, the method includes detecting a level of RNA in a blood, comparing the sample level of RNA to a reference level of RNA specific to a tissue, determining whether a difference exists between the sample level and the reference level, and characterizing the as abnormal if the sample level and the reference level are the same.
A deconvolution of a cell-free transcriptome is used to determine the relative contribution of each tissue type towards the cell-free RNA transcriptome. The following steps are employed to determine the relative RNA contributions of certain tissues in a sample. First, a panel of tissue-specific transcripts is identified. Second, total RNA in plasma from a sample is determined using methods known in the art. Third, the total RNA is assessed against the panel of tissue-specific transcripts, and the total RNA is considered a summation these different tissue-specific transcripts. Quadratic programming can be used as a constrained optimization method to deduce the relative optimal contributions of different organs/tissues towards the cell-free transcriptome of the sample.
One or more databases of genetic information can be used to identify a panel of tissue-specific transcripts. Accordingly, aspects of the invention provide systems and methods for the use and development of a database. Particularly, methods of the invention utilize databases containing existing data generated across tissue types to identify the tissue-specific genes. Databases utilized for identification of tissue-specific genes include the Human 133A/GNF1H Gene Atlas and RNA-Seq Atlas, although any other database or literature can be used. In order to identify tissue-specific transcripts from one or more databases, certain embodiments employ a template-matching algorithm to the databases. Template matching algorithms used to filter data are known in the art, see e.g., Pavlidis P, Noble W S (2001) Analysis of strain and regional variation in gene expression in mouse brain. Genome Biol 2:research0042.1-0042.15.
In certain embodiments, quadratic programming is used as a constrained optimization method to deduce relative optimal contributions of different organs/tissues towards the cell-free transcriptome in a sample. Quadratic programming is known in the art and described in detail in Goldfarb and A. Idnani (1982). Dual and Primal-Dual Methods for Solving Strictly Convex Quadratic Programs. In J. P. Hennart (ed.), Numerical Analysis, Springer-Verlag, Berlin, pages 226-239, and D. Goldfarb and A. Idnani (1983). A numerically stable dual method for solving strictly convex quadratic programs. Mathematical Programming, 27,1-33.
Method of the invention for determining the relative contributions for each tissue can be used to determine the reference level for the tissue. That is, a certain population of subjects (e.g., maternal, normal, and cancerous) can be subject to the deconvolution process outlined in
According to certain embodiments, an unknown sample of blood can be subject to process outlined in
The biological sample can be blood, saliva, sputum, urine, semen, transvaginal fluid, cerebrospinal fluid, sweat, breast milk, breast fluid (e.g., breast nipple aspirate), stool, a cell or a tissue biopsy. In certain embodiments, the samples of the same biological sample are obtained at multiple different time points in order to analyze differential transcript levels in the biological sample over time. For example, maternal plasma may be analyzed in each trimester. In some embodiments, the biological sample is drawn blood and circulating nucleic acids, such as cell-free RNA. The cell-free RNA may be from different genomic sources is found in the blood or plasma, rather than in cells.
In a particular embodiment, the drawn blood is maternal blood. In order to obtain a sufficient amount of nucleic acids for testing, it is preferred that approximately 10-50 mL of blood be drawn. However, less blood may be drawn for a genetic screen in which less statistical significance is required, or in which the RNA sample is enriched for fetal RNA.
Methods of the invention involve isolating total RNA from a biological sample. Total RNA can be isolated from the biological sample using any methods known in the art. In certain embodiments, total RNA is extracted from plasma. Plasma RNA extraction is described in Enders et al., “The Concentration of Circulating Corticotropin-releasing Hormone mRNA in Maternal Plasma Is Increased in Preeclampsia,” Clinical Chemistry 49: 727-731, 2003. As described there, plasma harvested after centrifugation steps is mixed Trizol LS reagent (Invitrogen) and chloroform. The mixture is centrifuged, and the aqueous layer transferred to new tubes. Ethanol is added to the aqueous layer. The mixture is then applied to an RNeasy mini column (Qiagen) and processed according to the manufacturer's recommendations.
In the embodiments where the biological sample is maternal blood, the maternal blood may optionally be processed to enrich the fetal RNA concentration in the total RNA. For example, after extraction, the RNA can be separated by gel electrophoresis and the gel fraction containing circulatory RNA with a size of corresponding to fetal RNA (e.g., <300 bp) is carefully excised. The RNA is extracted from this gel slice and eluted using methods known in the art.
Alternatively, fetal specific RNA may be concentrated by known methods, including centrifugation and various enzyme inhibitors. The RNA is bound to a selective membrane (e.g., silica) to separate it from contaminants. The RNA is preferably enriched for fragments circulating in the plasma, which are less than less 300 bp. This size selection is done on an RNA size separation medium, such as an electrophoretic gel or chromatography material.
Flow cytometry techniques can also be used to enrich for fetal cells in maternal blood (Herzenberg et al., PNAS 76: 1453-1455 (1979); Bianchi et al., PNAS 87: 3279-3283 (1990); Bruch et al., Prenatal Diagnosis 11: 787-798 (1991)). U.S. Pat. No. 5,432,054 also describes a technique for separation of fetal nucleated red blood cells, using a tube having a wide top and a narrow, capillary bottom made of polyethylene. Centrifugation using a variable speed program results in a stacking of red blood cells in the capillary based on the density of the molecules. The density fraction containing low-density red blood cells, including fetal red blood cells, is recovered and then differentially hemolyzed to preferentially destroy maternal red blood cells. A density gradient in a hypertonic medium is used to separate red blood cells, now enriched in the fetal red blood cells from lymphocytes and ruptured maternal cells. The use of a hypertonic solution shrinks the red blood cells, which increases their density, and facilitates purification from the more dense lymphocytes. After the fetal cells have been isolated, fetal RNA can be purified using standard techniques in the art.
Further, an agent that stabilizes cell membranes may be added to the maternal blood to reduce maternal cell lysis including but not limited to aldehydes, urea formaldehyde, phenol formaldehyde, DMAE (dimethylaminoethanol), cholesterol, cholesterol derivatives, high concentrations of magnesium, vitamin E, and vitamin E derivatives, calcium, calcium gluconate, taurine, niacin, hydroxylamine derivatives, bimoclomol, sucrose, astaxanthin, glucose, amitriptyline, isomer A hopane tetral phenylacetate, isomer B hopane tetral phenylacetate, citicoline, inositol, vitamin B, vitamin B complex, cholesterol hemisuccinate, sorbitol, calcium, coenzyme Q, ubiquinone, vitamin K, vitamin K complex, menaquinone, zonegran, zinc, ginkgo biloba extract, diphenylhydantoin, perftoran, polyvinylpyrrolidone, phosphatidylserine, tegretol, PABA, disodium cromglycate, nedocromil sodium, phenyloin, zinc citrate, mexitil, dilantin, sodium hyaluronate, or polaxamer 188.
An example of a protocol for using this agent is as follows: The blood is stored at 4° C. until processing. The tubes are spun at 1000 rpm for ten minutes in a centrifuge with braking power set at zero. The tubes are spun a second time at 1000 rpm for ten minutes. The supernatant (the plasma) of each sample is transferred to a new tube and spun at 3000 rpm for ten minutes with the brake set at zero. The supernatant is transferred to a new tube and stored at −80° C. Approximately two milliliters of the “buffy coat,” which contains maternal cells, is placed into a separate tube and stored at −80° C.
Methods of the invention also involve preparing amplified cDNA from total RNA. cDNA is prepared and indiscriminately amplified without diluting the isolated RNA sample or distributing the mixture of genetic material in the isolated RNA into discrete reaction samples. Preferably, amplification is initiated at the 3′ end as well as randomly throughout the whole transcriptome in the sample to allow for amplification of both mRNA and non-polyadenylated transcripts. The double-stranded cDNA amplification products are thus optimized for the generation of sequencing libraries for Next Generation Sequencing platforms. Suitable kits for amplifying cDNA in accordance with the methods of the invention include, for example, the Ovation® RNA-Seq System.
Methods of the invention also involve sequencing the amplified cDNA. While any known sequencing method can be used to sequence the amplified cDNA mixture, single molecule sequencing methods are preferred. Preferably, the amplified cDNA is sequenced by whole transcriptome shotgun sequencing (also referred to herein as (“RNA-Seq”). Whole transcriptome shotgun sequencing (RNA-Seq) can be accomplished using a variety of next-generation sequencing platforms such as the Illumina Genome Analyzer platform, ABI Solid Sequencing platform, or Life Science's 454 Sequencing platform.
Methods of the invention further involve subjecting the cDNA to digital counting and analysis. The number of amplified sequences for each transcript in the amplified sample can be quantitated via sequence reads (one read per amplified strand). Unlike previous methods of digital analysis, sequencing allows for the detection and quantitation at the single nucleotide level for each transcript present in a biological sample containing a genetic material from different genomic sources and therefore multiple transcriptomes.
After digital counting, the ratios of the various amplified transcripts can compared to determine relative amounts of differential transcript in the biological sample. Where multiple biological samples are obtained at different time-points, the differential transcript levels can be characterized over the course of time.
Differential transcript levels within the biological sample can also be analyzed using via microarray techniques. The amplified cDNA can be used to probe a microarray containing gene transcripts associated with one or conditions or diseases, such as any prenatal condition, or any type of cancer, inflammatory, or autoimmune disease.
It will be understood that methods and any flow diagrams disclosed herein can be implemented by computer program instructions. These program instructions may be provided to a computer processor, such that the instructions, which execute on the processor, create means for implementing the actions specified in the flowchart blocks or described in methods for assessing tissue disclosed herein. The computer program instructions may be executed by a processor to cause a series of operational steps to be performed by the processor to produce a computer implemented process. The computer program instructions may also cause at least some of the operational steps to be performed in parallel. Moreover, some of the steps may also be performed across more than one processor, such as might arise in a multi-processor computer system. In addition, one or more processes may also be performed concurrently with other processes or even in a different sequence than illustrated without departing from the scope or spirit of the invention.
The computer program instructions can be, stored on any suitable computer-readable medium including, but not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, digital versatile disks (DVD) or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can be accessed by a computing device.
EXAMPLES Example 1: Profiling Maternal Plasma Cell-Free RNA by RNA Sequencing-A Comprehensive Approach Overview:The plasma RNA profiles of 5 pregnant women were collected during the first trimester, second trimester, post-partum, as well as those of 2 non-pregnant female donors and 2 male donors using both microarray and RNA-Seq.
Among these pregnancies, there were 2 pregnancies with clinical complications such as premature birth and one pregnancy with bi-lobed placenta. Comparison of these pregnancies against normal cases reveals genes that exhibit significantly different gene expression pattern across different temporal stages of pregnancy. Application of such technique to samples associated with complicated pregnancies may help identify transcripts that can be used as molecular markers that are predictive of these pathologies.
Study Design and Methods: SubjectsSamples were collected from 5 pregnant women were during the first trimester, second trimester, third trimester, and post-partum. As a control, blood plasma samples were also collected from 2 non-pregnant female donors and 2 male donors.
Blood Collection and ProcessingBlood samples were collected in EDTA tube and centrifuged at 1600 g for 10 min at 4° C. Supernatant were placed in 1 ml aliquots in a 1.5 ml microcentrifuge tube which were then centrifuged at 16000 g for 10 min at 4° C. to remove residual cells. Supernatants were then stored in 1.5 ml microcentrifuge tubes at −80° C. until use.
RNA Extraction and AmplificationThe cell-free maternal plasma RNAs was extracted by Trizol LS reagent. The extracted and purified total RNA was converted to cDNA and amplified using the RNA-Seq Ovation Kit (NuGen). (The above steps were the same for both Microarray and RNA-Seq sample preparation).
The cDNA was fragmented using DNase I and labeled with Biotin, following by hybridization to Affymetrix GeneChip ST 1.0 microarrays. The Illumina sequencing platform and standard Illumina library preparation protocols were used for sequencing.
Data Analysis: Correlation Between Microarray and RNA-SeqThe RMA algorithm was applied to process the raw microarray data for background correction and normalization. RPKM values of the sequenced transcripts were obtained using the CASAVA 1.7 pipeline for RNA-seq. The RPKM in the RNA-Seq and the probe intensities in the microarray were converted to log 2 scale. For the RNA-Seq data, to avoid taking the log of 0, the gene expressions with RPKM of 0 were set to 0.01 prior to taking logs. Correlation coefficients between these two platforms ranges were then calculated.
Differential Expression of RNA Transcripts Levels Ssing RNA-seqDifferential gene expression analysis was performed using edgeR, a set of library functions which are specifically written to analyze digital gene expression data. Gene Ontology was then performed using DAVID to identify for significantly enriched GO terms.
Principle Component Analysis & Identification of Significant Time Varying GenesPrinciple component analysis was carried out using a custom script in R. To identify time varying genes, the time course library of functions in R were used to implement empirical Bayes methods for assessing differential expression in experiments involving time course which in our case are the different trimesters and post-partum for each individual patients.
Results and DiscussionRNA-Seq Reveals That Pregnancy-Associated Transcripts are Detected at Significantly Different Levels Between Pregnant and Non Pregnant Subjects.
A comparison of the transcripts level derived using RNA-Seq and Gene Ontology Analysis between pregnant and non-pregnant subjects revealed that transcripts exhibiting differential transcript levels are significantly associated with female pregnancy, suggesting that RNA-Seq are enabling observation of real differences between these two class of transcriptome due to pregnancy. The top rank significantly expressed gene is PLAC4 which has also been known as a target in previous studies for developing RNA based test for trisomy 21. A listing of the top detected female pregnancy associated differentially expressed transcripts is shown in
Principle Component Analysis (PCA) on Plasma Cell Free RNA Transcripts Levels in Maternal Plasma Distinguishes Between Pre-Mature and Normal Pregnancy
Using the plasma cell free transcript level profiles as inputs for Principle Component Analysis, the profile from each patient at different time points clustered into different pathological clusters suggesting that cell free plasma RNA transcript profile in maternal plasma may be used to distinguish between pre-term and non-preterm pregnancy.
Plasma Cell free RNA levels were quantified using both microarray and RNA-Seq. Transcripts expression levels profile from microarray and RNA-Seq from each patient are correlated with a Pearson correlation of approximately 0.7. Plots of the two main principal components for cell free RNA transcript levels is shown in
Identification of Cell Free RNA Transcripts in Maternal Plasma Exhibiting Significantly Different Time Varying Trends Between Pre-Term and Normal Pregnancy Across All Three Trimesters and Post Partum
A heatmap of the top 100 cell free transcript levels exhibiting different temporal levels in preterm and normal pregnancy using microarrays is shown in
Common Cell Free RNA Transcripts Identified by Microarray and RNA-Seq Which Exhibit Significantly Different Time Varying Trends Between Pre-Term and Normal Pregnancy Across All Three Trimesters and Post-Partum
A ranking of the top 20 transcripts differentially expressed between pre-term and normal pregnancy is shown in
The protein encoded by PVALB gene is a high affinity calcium ion-binding protein that is structurally and functionally similar to calmodulin and troponin C. The encoded protein is thought to be involved in muscle relaxation. As shown in
Results from quantification and characterization of maternal plasma cell-free RNA using RNA-Seq strongly suggest that pregnancy associated transcripts can be detected.
Furthermore, both RNA-Seq and microarray methods can detect considerable gene transcripts whose level showed differential time trends that has a high probability of being associated with premature births.
The methods described herein can be modified to investigate pregnancies of different pathological situations and can also be modified to investigate temporal changes at more frequent time points.
Example 2: Quantification of Tissue-Specific Cell-Free RNA Exhibiting Temporal Variation During Pregnancy Overview:Cell-free fetal DNA found in maternal plasma has been exploited extensively for non-invasive diagnostics. In contrast, cell-free fetal RNA which has been shown to be similarly detected in maternal circulation has yet been applied widely as a form of diagnostics. Both fetal cell-free RNA and DNA face similar challenges in distinguishing the fetal from maternal component because in both cases the maternal component dominates. To detect cell-free RNA of fetal origin, focus can be placed on genes that are highly expressed only during fetal development, which are subsequently inferred to be of fetal in origin and easily distinguished from background maternal RNA. Such a perspective is collaborated by studies that has established that cell-free fetal RNA derived from genes that are highly expressed in the placenta are detectable in maternal plasma during pregnancy.
A significant characteristic that set RNA apart from DNA can be attributed to RNA transcripts dynamic nature which is well reflected during fetal development. Life begins as a series of well-orchestrated events that starts with fertilization to form a single-cell zygote and ends with a multi-cellular organism with diverse tissue types. During pregnancy, majority of fetal tissues undergoes extensive remodeling and contain functionally diverse cell types. This underlying diversity can be generated as a result of differential gene expression from the same nuclear repertoire; where the quantity of RNA transcripts dictate that different cell types make different amount of proteins, despite their genomes being identical. The human genome comprises approximately 30,000 genes. Only a small set of genes are being transcribed to RNA within a particular differentiated cell type. These tissue specific RNA transcripts have been identified through many studies and databases involving developing fetuses of classical animal models. Combining known literature available with high throughput data generated from samples via sequencing, the entire collection of RNA transcripts contained within maternal plasma can be characterized.
Fetal organ formation during pregnancy depends on successive programs of gene expression. Temporal regulation of RNA quantity is necessary to generate this progression of cell differentiation events that accompany fetal organ genesis. To unravel similar temporal dynamics for cell free RNA, the expression profile of maternal plasma cell free RNA, especially the selected fetal tissue specific panel of genes, as a function across all three trimesters during pregnancy and post-partum were analyzed. Leveraging high throughput qPCR and sequencing technologies capability for simultaneous quantification of cell free fetal tissue specific RNA transcripts, a system level view of the spectrum of RNA transcripts with fetal origins in maternal plasma was obtained. In addition, maternal plasma was analyzed to deconvolute the heterogeneous cell free transcriptome of fetal origin a relative proportion of the different fetal tissue types. This approach incorporated physical constraints regarding the fetal contributions in maternal plasma, specifically the fraction of contribution of each fetal tissues were required to be non-negative and sum to one during all three trimesters of the pregnancy. These constraints on the data set enabled the results to be interpreted as relative proportions from different fetal organs. That is, a panel of previously selected fetal tissue-specific RNA transcripts exhibiting temporal variation can be used as a foundation for applying quadratic programing in order to determine the relative tissue-specific RNA contribution in one or more samples.
When considered individually, quantification of each of these fetal tissue specific transcripts within the maternal plasma can be used as a measure for the apoptotic rate of that particular fetal tissue during pregnancy. Normal fetal organ development is tightly regulated by cell division and apoptotic cell death. Developing tissues compete to survive and proliferate, and organ size is the result of a balance between cell proliferation and death. Due to the close association between aberrant cell death and developmental diseases, therapeutic modulation of apoptosis has become an area of intense research, but with this comes the demand for monitoring the apoptosis rate of specific. Quantification of fetal cell-free RNA transcripts provide such prognostic value, especially in premature births where the incidence of apoptosis in various organs of these preterm infants has been have been shown to contribute to neurodevelopmental deficits and cerebral palsy of preterm infants.
Sample Collection and Study Design
Selection of Fetal Tissue Specific Transcript Panel
To detect the presence of these fetal tissue-specific transcripts, a list of known fetal tissue specific genes was prepared from known literature and databases. The specificity for fetal tissues was validated by cross referencing between two main databases: TISGeD (Xiao, S.-J., Zhang, C. & Ji, Z.-L. TiSGeD: a Database for Tissue-Specific Genes. Bioinformatics (Oxford, England) 26, 1273-1275 (2010)) and BioGPS (Wu, C. et al. BioGPS: an extensible and customizable portal for querying and organizing gene annotation resources. Genome biology 10, R130 (2009); Su, A. I. et al. A gene atlas of the mouse and human protein-encoding transcriptomes. Proceedings of the National Academy of Sciences of the United States of America 101, 6062-7 (2004)). Most of these selected transcripts are associated with known fetal developmental processes. This list of genes was overlapped with RNA sequencing and microarray data to generate the panel of selected fetal tissue-specic transcripts shown in
Subjects
Samples of maternal blood were collected from normal pregnant women during the first trimester, second trimester, third trimester, and post-partum. For positive controls, fetal tissue specific RNA from the various fetal tissue types were bought from Agilent. Negative controls for the experiments were performed with the entire process with water, as well as with samples that did not undergoes the reverse transcription process.
Blood Collection and Processing
At each time-point, 7 to 15 mL of peripheral blood was drawn from each subject. Blood was centrifuged at 1600 g for 10 mins and transferred to microcentrifuge tubes for further centrifugation at 16000 g for 10 mins to remove residual cells. The above steps were carried out within 24 hours of the blood draw. Resulting plasma is stored at −80 Celsius for subsequent RNA extractions.
RNA Extraction
Cell free RNA extractions were carried using Trizol followed by Qiagen's RNeasy Mini Kit. To ensure that there are no contaminating DNA, DNase digestion is performed after RNA elution using RNase free DNase from Qiagen. Resulting cell free RNA from the pregnant subjects was then processed using standard microarrays and Illumina RNA-seq protocols. These steps generate the sequencing library that we used to generate RNA-seq data as well as the microarray expression data. The remaining cell free RNA are then used for parallel qPCR.
Parallel qPCR of Selected Transcripts
Accurate quantification of these fetal tissue specific transcripts was carried out using the Fluidigm BioMark system (See e.g. Spurgeon, S. L., Jones, R. C. & Ramakrishnan, R. High throughput gene expression measurement with real time PCR in a microfluidic dynamic array. PloS one 3, e1662 (2008)). This system allows for simultaneous query of a panel of fetal tissue specific transcripts. Two parallel forms of inquiry were conducted using different starting source of material. One was using the cDNA library from the Illumina sequencing protocol and the other uses the eluted RNA directly. Both sources of material were amplified with evagreen primers targeting the genes of interest. Both sources, RNA and cDNA, were preamplified. cDNA is preamplifed using evagreen PCR supermix and primers. RNA source is preamplified using the CellsDirect One-Step qRT-PCR kit from Invitrogen. Modifications were made to the default One-Step qRT-PCR protocol to accomodate a longer incubation time for reverse transcription. 19 cycles of preamplfication were conducted for both sources and the collected PCR products were cleaned up using Exonuclease I Treatment. To increase the dynamic range and the ability to quantify the efficiency of the later qPCR steps, serial dilutions were performed on the PCR products from 5 fold, 10 fold and 10 fold dilutions. Each of the collected maternal plasma from individual pregnant women across the time points went through the same procedures and was loaded onto 48×48 Dynamic Arrary Chips from Fluidigm to perform the qPCR. For positive control, fetal tissue specific RNA from the various fetal tissue types were bought from Agilent. Each of these RNA from fetal tissues went through the same preamplification and clean-up steps. A pool sample with equal proportions of different fetal tissues was created as well for later analysis to deconvolute the relative contribution of each tissue type in the pooled samples. All collected data from the Fluidigm BioMark system were pre-processed using Fluidigm Real Time PCR Analysis software to obtain the respective Ct values for each of the transcript across all samples. Negative controls of the experiments were performed with the entire process with water, as well as with samples that did not undergoes the reverse transcription process.
Data Analysis:
Fetal tissue specific RNA transcripts clear from the maternal peripheral bloodstream within a short period after birth. That is, the post-partum cell-free RNA transcriptome of maternal blood lacks fetal tissue specific RNA transcripts. As a result, it is expected that the quantity of these fetal tissue-specific transcripts to be higher before than after birth. The data of interest were the relative quantitative changes of the tissue specific transcripts across all three trimesters of pregnancy as compared to this baseline level after the baby is born. As described the methods, the fetal tissue-specific transcripts were quantified in parallel both using the actual cell-free RNA as well as the cDNA library of the same cell-free RNA. An example of the raw data obtained is shown in
Results and Discussion:
As shown in
Biological Significance of Quantified Fetal Tissue-Specific RNA: Most of the transcripts in the panel were involved in fetal organ development and many are also found within the amniotic fluid. Once such example is ZNF238. This transcript is specific to fetal brain tissue and is known to be vital for cerebral cortex expansion during embryogenesis when neuronal layers are formed. Loss of ZNF238 in the central nervous system leads to severe disruption of neurogenesis, resulting in a striking postnatal small-brain phenotype. Using methods of the invention, one can determine whether ZNF238 is presenting in healthy, normal levels according to the stage of development.
Known defects due to the loss of ZNF238 include a striking postnatal small-brain phenotype: microcephaly, agenesis of the corpus callosum and cerebellar hypoplasia. Microcephaly can sometimes be diagnosed before birth by prenatal ultrasound. In many cases, however, it might not be evident by ultrasound until the third trimester. Typically, diagnosis is not made until birth or later in infancy upon finding that the baby's head circumference is much smaller than normal. Microcephaly is a life-long condition and currently untreatable. A child born with microcephaly will require frequent examinations and diagnostic testing by a doctor to monitor the development of the head as he or she grows. Early detection of ZNf238 differential expression using methods of the invention provides for prenatal diagnosis and may hold prognostic value for drug treatments and dosing during course of treatment.
Beyond ZNF238, many of the characterized transcripts may hold diagnostic value in developmental diseases involving apoptosis, i.e., diseases caused by removal of unnecessary neurons during neural development. Seeing that apoptosis of neurons is essential during development, one could extrapolate that similar apoptosis might be activated in neurodegenerative diseases such as Alzheimer's disease, Huntington's disease, and amyotrophic lateral sclerosis. In such a scenario, the methodology described herein will allow for close monitoring for disease progression and possibly an ideal dosage according to the progression.
Deducing relative contributions of different fetal tissue types: Differential rate of apoptosis of specific tissues may directly correlate with certain developmental diseases. That is, certain developmental diseases may increase the levels of a particular specific RNA transcripts being observed in the maternal transcriptome. Knowledge of the relative contribution from various tissue types will allow for observations of these types of changes during the progression of these diseases. The quantified panel of fetal tissue specific transcripts during pregnancy can be considered as a summation of the contributions from the various fetal tissues.
Expressing
Yi=Σjπixij+ε
where Y is the observed transcript quantity in maternal plasma for gene i, X is the known transcript quantity for gene i in known fetal tissue j and ε the normally distributed error. Additional physical constraints includes:
-
- 1. Summation of all fraction contributing to the observed quantification is 1, given by the condition: Σπi=1
- 2. All the contribution from each tissue type has to greater than or equal zero. There is no physical meaning to having a negative contribution. This is given by πi≥0, since π is defined as the fractional contribution of each tissue types.
Consequently to obtain the optimal fractional contribution of each tissue type, the least-square error is minimized. The above equations are then solved using quadratic programming in R to obtain the optimal relative contributions of the tissue types towards the maternal cell free RNA transcripts. In the workflow, the quantity of RNA transcripts are given relative to the housekeeping genes in terms of Ct values obtained from qPCR. Therefore, the Ct value can be considered as a proxy of the measured transcript quantity. An increase in Ct value of one is similar to a two-fold change in transcript quantity, i.e. 2 raised to the power of 1. The process beings with normalizing all of the data in CT relative to the housekeeping gene, and is followed by quadratic programming.
As a proof of concept for the above scheme, different fetal tissue types (Brain, Placenta, Liver, Thymus, Lung) were mixed in equal proportions to generate a pool sample. Each fetal tissue types (Brain, Placenta, Liver, Thymus, Lung) along with the pooled sample were quantified using the same Fluidigm Biomark System to obtain the Ct values from qPCR for each fetal tissue specific transcript across all tissues and the pooled sample. These values were used to perform the same deconvolution. The resulting fetal fraction of each of the fetal tissue organs (Brain, Placenta, Liver, Thymus, Lung) was 0.109, 0.206, 0.236, 0.202 & 0.245 respectively.
Conclusion:
In summary, the panel of fetal specific cell free transcripts provides valuable biological information across different fetal tissues at once. Most particularly, the method can deduce the different relative proportions of fetal tissue-specific transcripts to total RNA, and, when considered individually, each transcript can be indicative of the apoptotic rate of the fetal tissue. Such measurements have numerous potential applications for developmental and fetal medicine. Most human fetal development studies have relied mainly on postnatal tissue specimens or aborted fetuses. Methods described herein provide quick and rapid assay of the rate of fetal tissue/organ growth or death on live fetuses with minimal risk to the pregnant mother and fetus. Similar methods may be employed to monitor major adult organ tissue systems that exhibit specific cell free RNA transcripts in the plasma.
Example 3: Deconvolution of Adult Cell-Free Transcriptome OverviewThe plasma RNA profiles of 4 healthy, normal adults were analyzed. Based on the gene expression profile of different tissue types, the methods described quantify the relative contributions of each tissue type towards the cell-free RNA component in a donor's plasma. For quantification, apoptotic cells from different tissue types are assumed to release their RNA into the plasma. Each of these tissues expressed a specific number of genes unique to the tissue type, and the observed cell-free RNA transcriptome is a summation of these different tissue types.
Study Design and Methods:To determine the contribution of tissue-specific transcripts to the cell-free adult transriptome, a list of known tissue-specific genes was prepared from known literature and databases. Two database sources were utilized: Human U133A/GNF1H Gene Atlas and RNA-Seq Atlas. Using the raw data from these two database, tissue-specific genes were identified by the following method. A template-matching process was applied to data obtained from the two databases for the purpose of identifying tissue-specific gene. The list of tissue specific genes identified by the method is given in
The adult cell-free transcriptome can be considered as a summation of the tissue-specific transcripts obtained from the two databases. To quantitatively deduce the relative proportions of the different tissues in an adult cell-free transcriptome, quadratic programming is performed as a constrained optimization method to deduce the relative optimal contributions of different organs/tissues towards the cell free-transcriptome. The specificity and accuracy of this process is dependent on the table of genes provided in Figure X and the extent by which that they are detectable in RNA-seq and microarray.
Subjects: Plasma Samples Were Collected From 4 Healthy, Normal Adults.
Initial Results:Deconvolution of our adult cell-free RNA transcriptome from microarray using the above methods revealed the relative contributions of the different tissue and organs are tabulated in
An interesting tissue type contributing to circulating RNA is the hypothalamus. The hypothalamus is bounded by specialized brain regions that lack an effective blood-brain barrier; the capillary endothelium at these sites is fenestrated to allow free passage of even large proteins and other molecules which in our case we believed that RNA transcripts from apoptotic cells in that region could be released into the plasma cell free RNA component.
The same methods were performed on the subjects using RNA-seq. The results described herein are limited due to the amount of tissue-specific RNA-Seq data available. However, it is understood that tissue-specific data is expanding with the increasing rate of sequencing of various tissue rates, and future analysis will be able to leverage those datasets. For RNA-seq data (as compared to microarray), whole blood nor the bone marrow samples are not available. The cell free transcriptome can only be decomposed to the available 11 different tissue types of RNA-seq data. Of which, only relative contributions from the hypothalamus and spleen were observed, as shown in
A list of 94 tissue-specific genes (as shown in
References and citations to other documents, such as patents, patent applications, patent publications, journals, books, papers, web contents, have been made throughout this disclosure. All such documents are hereby incorporated herein by reference in their entirety for all purposes.
EQUIVALENTSThe invention may be embodied in other specific forms without departing from the spirit or essential characteristics thereof. The foregoing embodiments are therefore to be considered in all respects illustrative rather than limiting on the invention described herein. Scope of the invention is thus indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein.
Claims
1. A method for determining that a pregnant subject or a fetus of said pregnant subject is at an elevated risk of a pregnancy-related health condition, comprising:
- (a) obtaining a biological sample from said pregnant subject;
- (b) obtaining sequence information of nucleic acid molecules derived from said biological sample to identify a set of pregnancy-related markers in said biological sample; and
- (c) using a programmed computer to deconvolve said set of pregnancy-related markers identified in (b) to determine tissue-specific contributions of one or more tissue types among said set of pregnancy-related markers, thereby determining that said pregnant subject or said fetus of said pregnant subject is at said elevated risk of said pregnancy-related health condition.
2. The method of claim 1, wherein said biological sample is selected from the group consisting of a blood sample, serum sample, plasma sample, saliva sample, stool sample, sputum sample, urine sample, semen sample, transvaginal fluid sample, breast milk sample, breast fluid sample, breast nipple aspirate sample, cerebrospinal fluid sample, sweat sample, a cell biopsy sample, and a tissue biopsy sample.
3. The method of claim 2, wherein said biological sample is a blood sample or a plasma sample.
4. The method of claim 1, wherein said set of pregnancy-related markers comprises a set of ribonucleic acid (RNA) transcripts.
5. The method of claim 4, wherein obtaining said sequence information of said nucleic acid molecules comprises reverse transcribing RNA molecules derived from said biological sample to produce complementary deoxyribonucleic acid (cDNA) molecules, and sequencing said cDNA molecules or derivatives thereof to identify said set of RNA transcripts.
6. The method of claim 5, further comprising amplifying said cDNA molecules to produce amplified products, and sequencing said amplified products or derivatives thereof to identify said set of RNA transcripts.
7. The method of claim 1, wherein said pregnancy-related health condition comprises pre-term birth.
8. The method of claim 7, wherein said pre-term birth is indicated by pre-eclampsia.
9. The method of claim 1, wherein said pregnancy-related health condition comprises pre-eclampsia.
10. The method of claim 1, wherein said pregnancy-related health condition comprises one or more anomalies in pregnancy
11. The method of claim 1, wherein said pregnancy-related health condition comprises one or more anomalies in fetal development.
12. The method of claim 1, wherein said set of pregnancy-related markers comprises one or more markers corresponding to tissue-specific differentially expressed genes.
13. The method of claim 12, wherein said one or more tissue-specific differentially expressed genes comprise one or more maternal-specific genes.
14. The method of claim 12, wherein said one or more tissue-specific differentially expressed genes comprise one or more placental-specific genes.
15. The method of claim 12, wherein said one or more tissue-specific differentially expressed genes comprise one or more fetal-specific genes.
16. The method of claim 15, wherein said one or more fetal-specific genes comprise one or more genes corresponding to a fetal tissue type selected from the group consisting of: brain, liver, thymus, and lung.
17. The method of claim 1, further comprising identifying a clinical intervention for said pregnant subject or said fetus of said pregnant subject based at least in part on said elevated risk of pregnancy-related health condition determined in (c).
18. The method of claim 17, wherein said clinical intervention comprises a drug treatment.
19. The method of claim 1, wherein said pregnant subject is in a first trimester of pregnancy.
20. The method of claim 1, wherein said pregnant subject is in a second trimester of pregnancy.
21. The method of claim 1, wherein said pregnant subject is in a third trimester of pregnancy.
22. The method of claim 1, wherein deconvolving said set of pregnancy-related markers comprises comparing processing said set of pregnancy-related markers with a reference.
23. The method of claim 22, wherein said reference corresponds to a set of pregnancy-related markers from one or more non-pregnant subjects or one or more pregnant subjects.
24. The method of claim 22, wherein deconvolving said set of pregnancy-related markers comprises identifying a difference between said set of pregnancy-related markers and said reference, and using said difference to determine said elevated risk of said pregnancy-related health condition.
25. The method of claim 24, further comprising determining a level of fold change in quantitative polymerase chain reaction (qPCR) measurements based at least in part on data corresponding to said set of pregnancy-related markers and said reference to identify said difference.
26. The method of claim 24, further comprising performing principle component analysis on data corresponding to said set of pregnancy-related markers and said reference to identify said difference.
27. The method of claim 1, wherein deconvolving said set of pregnancy-related markers comprises determining a relative contribution of a plurality of distinct fetal tissue types, and determining that said pregnant subject or said fetus of said pregnant subject is at said elevated risk of said pregnancy-related health condition based at least in part on said relative contribution of said plurality of distinct fetal tissue types.
28. The method of claim 27, wherein said deconvolving comprises performing a constrained optimization.
29. The method of claim 1, further comprising using a programmed computer to deconvolve sets of pregnancy-related markers identified in biological samples obtained from said pregnant subject at two or more different time points to determine sets of tissue-specific contributions of one or more tissue types among said sets of pregnancy-related markers, and comparing said sets of tissue-specific contributions to each other to determine that said pregnant subject or said fetus of said pregnant subject is at said elevated risk of said pregnancy-related health condition.
30. A method for monitoring a pregnant subject or a fetus of said pregnant subject, comprising:
- (a) obtaining a biological sample from said pregnant subject;
- (b) obtaining sequence information of nucleic acid molecules derived from said biological sample to identify a set of pregnancy-related markers in said biological sample; and
- (c) using a programmed computer to deconvolve said set of pregnancy-related markers identified in (b) to determine tissue-specific contributions of one or more tissue types among said set of pregnancy-related markers, thereby monitoring a health of said pregnant subject or a health or developmental state of said fetus of said pregnant subject.
Type: Application
Filed: Aug 5, 2020
Publication Date: Nov 19, 2020
Inventors: Lian Chye Winston Koh (Stanford, CA), Stephen R. Quake (Stanford, CA), Hei-Mun Christina Fan (Fremont, CA), Wenying Pan (Stanford, CA)
Application Number: 16/986,042
