Many cancers are caused by genetic dysregulation of canonical genes known to be associated with the cancer. However, identifying how genetic dysregulation is linked to cancer pathology under these circumstances. Furthermore, providing an effective therapeutic remains a challenging endeavor. Accordingly, new methods of diagnosis and treatment are needed to better understand how these genetic dysregulations cause a wide range of cancers.
In one aspect, the invention features a method of treating a cancer in a by identifying a sequence of a novel open reading frame (nORF) and a cancer associated therewith, wherein the sequence of the nORF is distinct from a canonical open reading frame (cORF) of a gene. The nORF is present in (i) an overlapping region of the cORF in an alternate reading frame, (ii) a 5′ untranslated region (UTR) of the cORF, (iii) a 3′ UTR of the cORF, (iv) an intronic region of the cORF, (v) an intergenic region of the cORF, or (vi) a region not associated with the cORF or the gene, and wherein the nORF has increased expression relative to the nORF in a noncancerous cell. The method further includes administering to the subject an inhibitor that reduces expression of the nORF to treat the cancer.
In another aspect, the invention features method of treating a cancer in a subject by administering to the subject an inhibitor that reduces expression of a nORF. The subject may have previously been identified with a sequence of the nORF and a cancer associated therewith, wherein the sequence of the nORF is distinct from a cORF of a gene, wherein the nORF is present in (i) an overlapping region of the cORF in an alternate reading frame, (ii) a 5′ UTR of the cORF, (iii) a 3′ UTR of the cORF, (iv) an intronic region of the cORF, (v) an intergenic region of the cORF, or (vi) a region not associated with the cORF or the gene, and wherein the nORF has increased expression relative to the nORF in a noncancerous cell.
In some embodiments of either of the foregoing aspects, the method reduces expression of the nORF, e.g., by at least 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 95%, 97%, or 99%. The nORF may exhibit an increase (e.g. by at least 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 100%, 200%, 300%, 400%, 500%, or more) in expression, e.g., as compared to the nORF in normal (e.g., noncancerous) tissue.
In some embodiments of either of the above aspects, the inhibitor is a small molecule, a polynucleotide, or a polypeptide. The polynucleotide may include a miRNA, an antisense RNA, an shRNA, or an siRNA. The polypeptide may include an antibody or antigen-binding fragment thereof (e.g., an scFv).
In some embodiments, the inhibitor is encoded by a vector, such as a viral vector. The viral vector may be selected, for example, from the group consisting of a Retroviridae family virus, an adenovirus, a parvovirus, a coronavirus, a rhabdovirus, a paramyxovirus, a picornavirus, an alphavirus, a herpes virus, and a poxvirus. The parvovirus viral vector may be, for example, an adeno-associated virus (AAV) vector.
In some embodiments, the viral vector is a Retroviridae family viral vector (e.g., a lentiviral vector, an alpharetroviral vector, or a gammaretroviral vector). The Retroviridae family viral vector may include one or more of the following: a central polypurine tract, a woodchuck hepatitis virus post-transcriptional regulatory element, a 5′-LTR, HIV signal sequence, HIV Psi signal 5′-splice site, delta-GAG element, 3′-splice site, and a 3′-self inactivating LTR.
In some embodiments, the viral vector is a pseudotyped viral vector. The pseudotyped viral vector may be selected, for example, from the group consisting of a pseudotyped adenovirus, a pseudotyped parvovirus, a pseudotyped coronavirus, a pseudotyped rhabdovirus, a pseudotyped paramyxovirus, a pseudotyped picornavirus, a pseudotyped alphavirus, a pseudotyped herpes virus, a pseudotyped poxvirus, and a pseudotyped Retroviridae family virus. The pseudotyped viral vector may be a lentiviral vector.
In some embodiments, the pseudotyped viral vector includes one or more envelope proteins from a virus selected from vesicular stomatitis virus (VSV), RD114 virus, murine leukemia virus (MLV), feline leukemia virus (FeLV), Venezuelan equine encephalitis virus (VEE), human foamy virus (HFV), walleye dermal sarcoma virus (WDSV), Semliki Forest virus (SFV), Rabies virus, avian leukosis virus (ALV), bovine immunodeficiency virus (BIV), bovine leukemia virus (BLV), Epstein-Barr virus (EBV), Caprine arthritis encephalitis virus (CAEV), Sin Nombre virus (SNV), Cherry Twisted Leaf virus (ChTLV), Simian T-cell leukemia virus (STLV), Mason-Pfizer monkey virus (MPMV), squirrel monkey retrovirus (SMRV), Rous-associated virus (RAV), Fujinami sarcoma virus (FuSV), avian carcinoma virus (MH2), avian encephalomyelitis virus (AEV), Alfa mosaic virus (AMV), avian sarcoma virus CT10, and equine infectious anemia virus (EIAV).
In some embodiments, the pseudotyped viral vector includes a VSV-G envelope protein.
In another aspect, the invention features a method of treating a cancer in a subject by identifying a sequence of a nORF and a cancer associated therewith, wherein the sequence of the nORF is distinct from a cORF of a gene. The nORF is present in (i) an overlapping region of the cORF in an alternate reading frame, (ii) a 5′ UTR of the cORF, (iii) a 3′ UTR of the cORF, (iv) an intronic region of the cORF, (v) an intergenic region of the cORF, or (vi) a region not associated with the cORF or the gene, and wherein the nORF has decreased expression relative to the nORF in a noncancerous cell. The method further includes administering to the subject an activator that increases expression of nORF to treat the cancer.
In another aspect, the invention features a method of treating a cancer in a subject by administering to the subject an activator that increases expression of a nORF. The subject may have previously been identified with a sequence of the nORF and a cancer associated therewith, wherein the sequence of the nORF is distinct from a cORF of a gene, wherein the nORF is present in (i) an overlapping region of the cORF in an alternate reading frame, (ii) a 5′ UTR of the cORF, (iii) a 3′ UTR of the cORF, (iv) an intronic region of the cORF, (v) an intergenic region of the cORF, or (vi) a region not associated with the cORF or the gene, and wherein the nORF has decreased expression relative to the nORF in a noncancerous cell.
In some embodiments of either of the foregoing aspects, the method increases expression of the nORF, e.g., by at least 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 100%, 200%, 300%, 400%, 500%, or more. The nORF may exhibit a decrease (e.g., by at least 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 95%, 97%, or 99%) in expression, e.g., as compared to the nORF in normal (e.g., noncancerous) tissue.
In some embodiments, the activator is a small molecule, a polynucleotide, or a polypeptide. The polynucleotide may include an antisense RNA. The polypeptide may include an antibody or antigen-binding fragment thereof (e.g., an scFv).
In some embodiments, the activator is encoded by a vector, such as a viral vector. The viral vector may be selected, for example, from the group consisting of a Retroviridae family virus, an adenovirus, a parvovirus, a coronavirus, a rhabdovirus, a paramyxovirus, a picornavirus, an alphavirus, a herpes virus, and a poxvirus. The parvovirus viral vector may be, for example, an AAV vector.
In some embodiments, the viral vector is a Retroviridae family viral vector (e.g., a lentiviral vector, an alpharetroviral vector, or a gammaretroviral vector). The Retroviridae family viral vector may include one or more of the following: a central polypurine tract, a woodchuck hepatitis virus post-transcriptional regulatory element, a 5′-LTR, HIV signal sequence, HIV Psi signal 5′-splice site, delta-GAG element, 3′-splice site, and a 3′-self inactivating LTR.
In some embodiments, the viral vector is a pseudotyped viral vector. The pseudotyped viral vector may be selected, for example, from the group consisting of a pseudotyped adenovirus, a pseudotyped parvovirus, a pseudotyped coronavirus, a pseudotyped rhabdovirus, a pseudotyped paramyxovirus, a pseudotyped picornavirus, a pseudotyped alphavirus, a pseudotyped herpes virus, a pseudotyped poxvirus, and a pseudotyped Retroviridae family virus. The pseudotyped viral vector may be a lentiviral vector.
In some embodiments, the pseudotyped viral vector includes one or more envelope proteins from a virus selected from VSV, RD114 virus, MLV, FeLV, VEE, HFV, WDSV, SFV, Rabies virus, ALV, BIV, BLV, EBV, CAEV, SNV, ChTLV, STLV, MPMV, SMRV, RAV, FuSV, MH2, AEV, AMV, avian sarcoma virus CT10, and EIAV.
In some embodiments, the pseudotyped viral vector includes a VSV-G envelope protein.
In another aspect, the invention features a method of treating a cancer in a subject by identifying a sequence of a nORF and a cancer associated therewith, wherein the sequence of the nORF is distinct from a cORF of a gene. The nORF is present in (i) an overlapping region of the cORF in an alternate reading frame, (ii) a 5′ UTR of the cORF, (iii) a 3′ UTR of the cORF, (iv) an intronic region of the cORF, (v) an intergenic region of the cORF, or (vi) a region not associated with the cORF or the gene, and wherein the nORF has decreased expression relative to the nORF in a noncancerous cell. The method further includes providing a protein encoded by the nORF to the subject treat the cancer.
In another aspect, the invention features a method of treating a cancer in a subject by providing a protein encoded by a nORF to the subject. The subject may have previously been identified with a sequence of the nORF and a cancer associated therewith, wherein the sequence of the nORF is distinct from a cORF of a gene, wherein the nORF is present in (i) an overlapping region of the cORF in an alternate reading frame, (ii) a 5′ UTR of the cORF, (iii) a 3′ UTR of the cORF, (iv) an intronic region of the cORF, (v) an intergenic region of the cORF, or (vi) a region not associated with the cORF or the gene, and wherein the nORF has decreased expression relative to the nORF in a noncancerous cell.
In some embodiments of either of the foregoing aspects, the method includes restoring the encoded protein product of the nORF. The method may include providing the protein product or a polynucleotide encoding the protein product. The method may include providing a vector (e.g., a viral vector) including the polynucleotide encoding the protein product.
In some embodiments, the viral vector may be selected, for example, from the group consisting of a Retroviridae family virus, an adenovirus, a parvovirus, a coronavirus, a rhabdovirus, a paramyxovirus, a picornavirus, an alphavirus, a herpes virus, and a poxvirus. The parvovirus viral vector may be, for example, an adeno-associated virus (AAV) vector.
In some embodiments, the viral vector is a Retroviridae family viral vector (e.g., a lentiviral vector, an alpharetroviral vector, or a gammaretroviral vector). The Retroviridae family viral vector may include one or more of the following: a central polypurine tract, a woodchuck hepatitis virus post-transcriptional regulatory element, a 5′-LTR, HIV signal sequence, HIV Psi signal 5′-splice site, delta-GAG element, 3′-splice site, and a 3′-self inactivating LTR.
In some embodiments, the viral vector is a pseudotyped viral vector. The pseudotyped viral vector may be selected, for example, from the group consisting of a pseudotyped adenovirus, a pseudotyped parvovirus, a pseudotyped coronavirus, a pseudotyped rhabdovirus, a pseudotyped paramyxovirus, a pseudotyped picornavirus, a pseudotyped alphavirus, a pseudotyped herpes virus, a pseudotyped poxvirus, and a pseudotyped Retroviridae family virus. The pseudotyped viral vector may be a lentiviral vector.
In some embodiments, the pseudotyped viral vector includes one or more envelope proteins from a virus selected from VSV, RD114 virus, MLV, FeLV, VEE, HFV, WDSV, SFV, Rabies virus, ALV, BIV, BLV, EBV, CAEV, SNV, ChTLV, STLV, MPMV, SMRV, RAV, FuSV, MH2, AEV, AMV, avian sarcoma virus CT10, and EIAV.
In some embodiments, the pseudotyped viral vector includes a VSV-G envelope protein.
In some embodiments of any of the above aspects, the encoded protein product of the nORF is less than about 100 amino acids.
In some embodiments, the method further includes performing a statistical analysis between the nORF and the cancer. The statistical analysis may measure a positive or negative association between the nORF and the cancer.
In some embodiments, the cancer is selected from the list consisting of breast invasive carcinoma, colon adenocarcinoma, esophageal carcinoma, head and neck squamous cell carcinoma, kidney chromophobe, kidney clear cell carcinoma, kidney papillary cell carcinoma, liver hepatocellular carcinoma, lung adenocarcinoma, lung squamous cell carcinoma, prostrate adenocarcinoma, stomach adenocarcinoma, thyroid carcinoma, and uterine corpus endometrioid carcinoma
In some embodiments, the nORF is selected from Table 1.
In some embodiments, the nORF is selected from Table 2.
In some embodiments, the nORF is selected from Table 3.
In some embodiments, the nORF is selected from Table 4.
In some embodiments, the nORF is selected from Table 5.
In some embodiments, the nORF is not HOXB-AS3.
In some embodiments, the cancer is not colorectal cancer.
In some embodiments, the nORF is not PINT87aa (LING-PINT).
In some embodiments, the cancer is not glioblastoma.
As used herein, a “novel open reading frame” or “nORF” refers to an open reading frame that is transcribed in a cell and consists of a sequence that is distinct from a canonical open reading frame (cORF) transcribed from a gene. The nORF may be present in (i) an overlapping region of the cORF in an alternate reading frame, (ii) a 5′ untranslated region (UTR) of the cORF, (iii) a 3′ UTR of the cORF, (iv) an intronic region of the cORF, (v) an intergenic region of the cORF, or (vi) a region not associated with the cORF or the gene. The nORF may be any unannotated genetic sequence that is transcribed in a cell.
As used herein, a “canonical open reading frame” or “cORF” refers to an open reading frame that is transcribed in a cell and its associated genetic elements, including the 5′ UTR, the 3′ UTR, the intronic regions, the exonic regions, and the intergenic regions flanking the gene comprising the cORF. A cORF includes either the primary open reading frame that is expressed from a gene, the most abundantly expressed open reading frame expressed from a gene, or an ORF that is annotated in a publicly available database as the primary and/or most abundantly expressed open reading frame from a gene.
Described herein are methods of diagnosing and treating a cancer associated with dysregulated novel open reading frames (nORFs). Many cancers are caused by dysregulation (e.g., upregulation or downregulation) in a gene or a genetic variant that is associated with the cancer. However, it was previously unclear how certain cancers are caused in which no dysregulation of a canonical gene or a canonical open reading frame (cORF) associated with the gene is present and no variant is known. The present invention is premised, in part, upon the discovery of dysregulation of certain novel open reading frames (nORFs) that are distinct from canonical open reading frames (cORF) of genes. In these instances, the dysregulation (e.g., upregulation or downregulation) imparts a deleterious effect on the nORF, in some instances, with or without substantially impacting a protein encoded by a cORF. In particular, the present invention features methods of treating cancer associated with a dysregulated nORF in which differential expression (e.g., increased or decreased expression) of the nORF is observed. With increased or decreased expression, the gene product encoded by the dysregulated nORF is increased or decreased as compared to the nORF, e.g., in a noncancerous cell. The methods of diagnosis and treatment are described in more detail below.
Genetic testing offers one avenue by which a patient may be diagnosed as having or is at risk of developing a particular cancer. For example, a genetic analysis can be used to determine whether a patient has a nORF associated with a cancer. The nORF may be present in any region of a gene, such as within the cORF, a 5′ untranslated region (UTR) of the cORF, a 3′ UTR of the cORF, an intronic region of the cORF, or an intergenic region of the cORF, The nORF may be present within an overlapping region of the cORF in an alternate reading frame, a 5′ UTR of the cORF, a 3′ UTR of the cORF, an intronic region of the cORF, or an intergenic region of the cORF. The nORF may be present in a region that is not associated with the cORF of the gene.
Exemplary genetic tests that can be used to determine whether a patient has such nORF include polymerase chain reaction (PCR) methods known in the art, such as DNA and RNA sequencing. nORF sequences may be identified de novo, e.g., using computational or statistical methods. Furthermore, nORF sequences may be identified from publicly available databases in genomic sequences in which the nORF was not previously identified and/or annotated as a sequence that was transcribed, and/or translated.
nORF sequences may be identified as being linked to a particular cancer by using a statistical analysis between the dysregulated nORF and the cancer. The statistical analysis may measure a positive or negative association between the dysregulated nORF and the cancer (see, e.g., Example 1).
To examine the functional importance of a nORF separately from a canonical coding sequence, datasets, such as the Genome Aggregation Database, may be used.
The invention features methods of treating a subject having a dysregulated nORF that has differential expression (e.g., increased or decreased expression). The dysregulated nORF may exhibit an increase (e.g., by at least 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 100%, 200%, 300%, 400%, 500%, or more) in expression, e.g., as compared to the nORF in normal (e.g., noncancerous) tissue. The dysregulated nORF may exhibit a decrease (e.g., by at least 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 95%, 97%, or 99%) in expression, e.g., as compared to the dysregulated nORF in normal (e.g., noncancerous) tissue. The subject may be first determined to have the dysregulated nORF and then may subsequently be treated for the cancer. The subject may have previously been determined to have the dysregulated nORF and is then treated for the cancer. The treatment varies according to the dysregulated nORF associated with the cancer. For example, the treatment may include an inhibitor that targets the dysregulated nORF to decrease (e.g., by at least 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 95%, 97%, or 99%) expression of an upregulated nORF. The treatment may include an activator that targets the dysregulated nORF to increase (e.g., by at least 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 100%, 200%, 300%, 400%, 500%, or more) expression of a downregulated nORF. Alternatively, or in addition, the treatment may include providing the nORF or a protein encoded by the nORF to restore levels of the nORF.
The methods of treatment and diagnosis described herein may include providing an inhibitor that targets the dysregulated nORF. The inhibitor may reduce (e.g., by at least 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 95%, 97%, or 99%) an amount or activity of the dysregulated nORF, such as to prevent the deleterious effect of the dysregulated nORF. The inhibitor may target the polynucleotide containing the nORF or the protein encoded by the nORF. The inhibitor may be a small molecule, a polynucleotide, or a polypeptide. Suitable small molecules may be determined or identified by using computational analysis based on the structure of the dysregulated nORF as determined by a protein folding algorithm. The small molecule may target any region of the dysregulated nORF. The small molecule may target the nORF or the protein encoded by the nORF. Suitable polypeptides for reducing an activity or amount of the dysregulated nORF include, for example, an antibody or antigen-binding fragment thereof that binds to the dysregulated nORF (e.g., a single chain antibody or antigen-binding fragment thereof). Suitable polynucleotides that can reduce an amount or activity of the dysregulated nORF include RNA. For example, an RNA for reducing an activity or amount of the dysregulated nORF may be, for example, a miRNA, an antisense RNA, an shRNA, or an siRNA. The miRNA, antisense RNA, shRNA, or siRNA may target a region of RNA (e.g., dysregulated nORF gene) to reduce expression of the dysregulated nORF. The polynucleotide may be an aptamer, e.g., an RNA aptamer that binds to and/or reduces an amount and/or activity of the dysregulated nORF or the protein encoded by the dysregulated nORF. The inhibitor may be provided directly or may be provided by a vector (e.g., a viral vector) encoding the inhibitor. The inhibitor may be formulated, e.g., in a pharmaceutical composition containing a pharmaceutically acceptable carrier. The composition can be administered by any suitable method known in the art to the skilled artisan. The composition (e.g., a vector, e.g., a viral vector) may be formulated in a virus or a virus-like particle.
Nucleic Acid Mediated Knockdown
Using the compositions and methods described herein, a patient with a cancer may be administered an interfering RNA molecule, a composition containing the same, or a vector encoding the same, so as to reduce or suppress the expression of a dysregulated nORF. Exemplary interfering RNA molecules that may be used in conjunction with the compositions and methods described herein are siRNA molecules, miRNA molecules, shRNA molecules, and antisense RNA molecules, among others. In the case of siRNA molecules, the siRNA may be single stranded or double stranded. miRNA molecules, in contrast, are single-stranded molecules that form a hairpin, thereby adopting a hydrogen-bonded structure reminiscent of a nucleic acid duplex. In either case, the interfering RNA may contain an antisense or “guide” strand that anneals (e.g., by way of complementarity) to the repeat-expanded mutant RNA target. The interfering RNA may also contain a “passenger” strand that is complementary to the guide strand and, thus, may have the same nucleic acid sequence as the RNA target.
siRNA is a class of short (e.g., 20-25 nt) double-stranded non-coding RNA that operates within the RNA interference pathway. siRNA may interfere with expression of the dysregulated nORF gene with complementary nucleotide sequences by degrading mRNA (via the Dicer and RISC pathways) after transcription, thereby preventing translation. miRNA is another short (e.g., about 22 nucleotides) non-coding RNA molecule that functions in RNA silencing and post-transcriptional regulation of gene expression. miRNAs function via base-pairing with complementary sequences within mRNA molecules, thereby leading to cleavage of the mRNA strand into two pieces and destabilization of the mRNA through shortening of its poly(A) tail. shRNA is an artificial RNA molecule with a tight hairpin turn that can be used to silence target gene expression via RNA interference. Antisense RNA are also short single stranded molecules that hybridize to a target RNA and prevent translation by occluding the translation machinery, thereby reducing expression of the target (e.g., the dysregulated nORF).
Antibody Mediated Knockdown
Using the compositions and methods described herein, a patient with a cancer may be provided an antibody or antigen-binding fragment thereof, a composition containing the same, a vector encoding the same, or a composition of cells containing a vector encoding the same, so as to suppress or reduce the activity of the dysregulated nORF. In some embodiments of the compositions and methods described herein, an antibody or antigen-biding fragment thereof may be used that binds to and reduces or eliminates the activity of the dysregulated nORF. The antibody may be monoclonal or polyclonal. In some embodiments, the antigen-binding fragment is an antibody that lacks the Fc portion, an F(ab′)2, a Fab, an Fv, or an scFv. The antigen-binding fragment may be an scFv.
One of ordinary skill in the art will appreciate that an antibody may include four polypeptides: two identical copies of a heavy chain polypeptide and two copies of a light chain polypeptide. Each of the heavy chains contains one N-terminal variable (VH) region and three C-terminal constant (CH1, CH2 and CH3) regions, and each light chain contains one N-terminal variable (VL) region and one C-terminal constant (CL) region. Thus, one of skill in the art would appreciate that as described herein, a vector that includes a transgene that encodes a polypeptide that is an antibody may be a single transgene that encodes a plurality of polypeptides. Also contemplated is a vector that includes a plurality of transgenes, each transgene encoding a separate polypeptide of the antibody. All variations are contemplated herein. The variable regions of each pair of light and heavy chains form the antigen binding site of an antibody. The transgene which encodes an antibody directed against the dysregulated nORF can include one or more transgene sequences, each of which encodes one or more of the heavy and/or light chain polypeptides of an antibody. In this respect, the transgene sequence which encodes an antibody directed against the dysregulated nORF can include a single transgene sequence that encodes the two heavy chain polypeptides and the two light chain polypeptides of an antibody. Alternatively, the transgene sequence which encodes an antibody directed against the dysregulated nORF can include a first transgene sequence that encodes both heavy chain polypeptides of an antibody, and a second transgene sequence that encodes both light chain polypeptides of an antibody. In yet another embodiment, the transgene sequence which encodes an antibody can include a first transgene sequence encoding a first heavy chain polypeptide of an antibody, a second transgene sequence encoding a second heavy chain polypeptide of an antibody, a third transgene sequence encoding a first light chain polypeptide of an antibody, and a fourth transgene sequence encoding a second light chain polypeptide of an antibody.
In some embodiments, the transgene that encodes the antibody includes a single open reading frame encoding a heavy chain and a light chain, and each chain is separated by a protease cleavage site.
In some embodiments, the transgene encodes a single open reading frame encoding both heavy chains and both light chains, and each chain is separate by protease cleavage site.
In some embodiments, full-length antibody expression can be achieved from a single transgene cassette using 2A peptides, such as foot-and-mouth disease virus (FMDV) equine rhinitis A, porcine teschovirus-1, and Thosea asigna virus 2A peptides, which are used to link two or more genes and allow the translated polypeptide to be self-cleaved into individual polypeptide chains (e.g., heavy chain and light chain, or two heavy chains and two light chains). Thus, in some embodiments, the transgene encodes a 2A peptide in between the heavy and light chains, optionally with a flexible linker flanking the 2A peptide (e.g., GSG linker). The transgene may further include one or more engineered cleavage sequences, e.g., a furin cleavage sequence to remove the 2A peptide residues attached to the heavy chain or light chain. Exemplary 2A peptides are described, e.g., in Chng et al MAbs 7: 403-412, 201f5, and Lin et al. Front. Plant Sci 9:1379, 2018, the disclosures of which are hereby incorporated by reference in their entirety.
In some embodiments, the antibody is a single-chain antibody or antigen-binding fragment thereof expressed from a single transgene.
The methods of treatment and diagnosis described herein may include providing an activator that targets the dysregulated nORF. The activator may increase (e.g., by at least 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 100%, 200%, 300%, 400%, 500%, or more) an amount or activity of the dysregulated nORF, such as to prevent the deleterious effect of the dysregulated nORF. The activator may target the polynucleotide containing the nORF or the protein encoded by the nORF. The activator may be a small molecule, a polynucleotide, or a polypeptide. Suitable small molecules may be determined or identified by using computational analysis based on the structure of the dysregulated nORF as determined by a protein folding algorithm. The small molecule may target any region of the dysregulated nORF. The small molecule may target the nORF or the protein encoded by the nORF. Suitable polypeptides for increasing an activity or amount of the dysregulated nORF include, for example, an antibody or antigen-binding fragment thereof that binds to the dysregulated nORF (e.g., a single chain antibody or antigen-binding fragment thereof). Suitable polynucleotides that can increase an amount or activity of the dysregulated nORF include RNA. For example, an RNA for increasing an activity or amount of the dysregulated nORF may be, for example, an antisense RNA. The antisense RNA may target a region of RNA (e.g., dysregulated nORF gene) upstream of the primary nORF open reading frame to reduce expression of the upstream nORFs, thereby dedicating the translation machinery to the primary nORF in order to increase expression of the primary nORF. The polynucleotide may be an aptamer, e.g., an RNA aptamer that binds to and/or increases an amount and/or activity of the dysregulated nORF or the protein encoded by the dysregulated nORF. The activator may be provided directly or may be provided by a vector (e.g., a viral vector) encoding the activator. The activator may be formulated, e.g., in a pharmaceutical composition containing a pharmaceutically acceptable carrier. The composition can be administered by any suitable method known in the art to the skilled artisan. The composition (e.g., a vector, e.g., a viral vector) may be formulated in a virus or a virus-like particle.
nORF Replacement
The present invention also features methods of treating a cancer by administering or providing a nORF or a protein encoded by the nORF. The therapy may restore the encoded protein product of the nORF, such as to replace the nORF that is no longer present due to downregulation. The therapy may include providing the protein product or a polynucleotide encoding the protein product. The method may include providing a vector (e.g., a viral vector) that encodes the protein product. Alternatively, the protein encoded by the nORF may be administered directly, e.g., as an enzyme replacement therapy. The nORF or a polynucleotide encoding the nORF (e.g., a vector, e.g., a viral vector) may be formulated, e.g., in a pharmaceutical composition containing a pharmaceutically acceptable carrier. The composition can be administered by any suitable method known in the art to the skilled artisan. The composition may be formulated in a virus or a virus-like particle.
In some embodiments, the length of the nORF is less than about 100 amino acids (e.g., from about 50 to 100, 50 to 90, 50 to 80, 60 to 90, 60 to 80, 70 to 100, 70 to 90, 70 to 80, 80 to 100, or 90 to 100 amino acids).
Viral genomes provide a rich source of vectors that can be used for the efficient delivery of exogenous genes into a mammalian cell. The gene to be delivered may include an activator or inhibitor that targets a dysregulated nORF, such as an RNA (e.g., an aptamer, a miRNA, an antisense RNA, an shRNA, or an siRNA). Alternatively, the gene to be delivered may include the nORF for replacement. Viral genomes are particularly useful vectors for gene delivery as the polynucleotides contained within such genomes are typically incorporated into the nuclear genome of a mammalian cell by generalized or specialized transduction. These processes occur as part of the natural viral replication cycle, and do not require added proteins or reagents in order to induce gene integration. Examples of viral vectors are a retrovirus (e.g., Retroviridae family viral vector), adenovirus (e.g., Ad5, Ad26, Ad34, Ad35, and Ad48), parvovirus (e.g., an adeno-associated viral (AAV) vector), coronavirus, negative strand RNA viruses such as orthomyxovirus (e.g., influenza virus), rhabdovirus (e.g., rabies and vesicular stomatitis virus), paramyxovirus (e.g. measles and Sendai), positive strand RNA viruses, such as picornavirus and alphavirus, and double stranded DNA viruses including adenovirus, herpesvirus (e.g., Herpes Simplex virus types 1 and 2, Epstein-Barr virus, cytomegalovirus), and poxvirus (e.g., vaccinia, modified vaccinia Ankara (MVA), fowlpox and canarypox). Other viruses include Norwalk virus, togavirus, flavivirus, reoviruses, papovavirus, hepadnavirus, human papilloma virus, human foamy virus, and hepatitis virus, for example. Examples of retroviruses are: avian leukosis-sarcoma, avian C-type viruses, mammalian C-type, B-type viruses, D-type viruses, oncoretroviruses, HTLV-BLV group, lentivirus, alpharetrovirus, gammaretrovirus, spumavirus (Coffin, J. M., Retroviridae: The viruses and their replication, Virology, Third Edition (Lippincott-Raven, Philadelphia, (1996))). Other examples are murine leukemia viruses, murine sarcoma viruses, mouse mammary tumor virus, bovine leukemia virus, feline leukemia virus, feline sarcoma virus, avian leukemia virus, human T-cell leukemia virus, baboon endogenous virus, Gibbon ape leukemia virus, Mason Pfizer monkey virus, simian immunodeficiency virus, simian sarcoma virus, Rous sarcoma virus and lentiviruses. Other examples of vectors are described, for example, in McVey et al., (U.S. Pat. No. 5,801,030), the teachings of which are incorporated herein by reference.
Retro viral Vectors
The delivery vector used in the methods described herein may be a retroviral vector. One type of retroviral vector that may be used in the methods and compositions described herein is a lentiviral vector. Lentiviral vectors (LVs), a subset of retroviruses, transduce a wide range of dividing and non-dividing cell types with high efficiency, conferring stable, long-term expression of the transgene encoding the polypeptide or RNA. An overview of optimization strategies for packaging and transducing LVs is provided in Delenda, The Journal of Gene Medicine 6: S125 (2004), the disclosure of which is incorporated herein by reference.
The use of lentivirus-based gene transfer techniques relies on the in vitro production of recombinant lentiviral particles carrying a highly deleted viral genome in which the agent of interest is accommodated. In particular, the recombinant lentivirus are recovered through the in trans coexpression in a permissive cell line of (1) the packaging constructs, i.e., a vector expressing the Gag-Pol precursors together with Rev (alternatively expressed in trans); (2) a vector expressing an envelope receptor, generally of an heterologous nature; and (3) the transfer vector, consisting in the viral cDNA deprived of all open reading frames, but maintaining the sequences required for replication, encapsidation, and expression, in which the sequences to be expressed are inserted.
A LV used in the methods and compositions described herein may include one or more of a 5′-Long terminal repeat (LTR), HIV signal sequence, HIV Psi signal 5-splice site (SD), delta-GAG element, Rev Responsive Element (RRE), 3′-splice site (SA), elongation factor (EF) 1-alpha promoter and 3′-self inactivating LTR (SIN-LTR). The lentiviral vector optionally includes a central polypurine tract (cPPT) and a woodchuck hepatitis virus post-transcriptional regulatory element (WPRE), as described in U.S. Pat. No. 6,136,597, the disclosure of which is incorporated herein by reference as it pertains to WPRE. The lentiviral vector may further include a pHR′ backbone, which may include for example as provided below.
The Lentigen LV described in Lu et al., Journal of Gene Medicine 6:963 (2004) may be used to express the DNA molecules and/or transduce cells. A LV used in the methods and compositions described herein may a 5′-Long terminal repeat (LTR), HIV signal sequence, HIV Psi signal 5′-splice site (SD), delta-GAG element, Rev Responsive Element (RRE), 3′-splice site (SA), elongation factor (EF) 1-alpha promoter and 3′-self inactivating L TR (SIN-LTR). It will be readily apparent to one skilled in the art that optionally one or more of these regions is substituted with another region performing a similar function.
Enhancer elements can be used to increase expression of modified DNA molecules or increase the lentiviral integration efficiency. The LV used in the methods and compositions described herein may include a nef sequence. The LV used in the methods and compositions described herein may include a cPPT sequence which enhances vector integration. The cPPT acts as a second origin of the (+)-strand DNA synthesis and introduces a partial strand overlap in the middle of its native HIV genome. The introduction of the cPPT sequence in the transfer vector backbone strongly increased the nuclear transport and the total amount of genome integrated into the DNA of target cells. The LV used in the methods and compositions described herein may include a Woodchuck Posttranscriptional Regulatory Element (WPRE). The WPRE acts at the transcriptional level, by promoting nuclear export of transcripts and/or by increasing the efficiency of polyadenylation of the nascent transcript, thus increasing the total amount of mRNA in the cells. The addition of the WPRE to LV results in a substantial improvement in the level of expression from several different promoters, both in vitro and in vivo. The LV used in the methods and compositions described herein may include both a cPPT sequence and WPRE sequence. The vector may also include an IRES sequence that permits the expression of multiple polypeptides from a single promoter.
In addition to IRES sequences, other elements which permit expression of multiple polypeptides are useful. The vector used in the methods and compositions described herein may include multiple promoters that permit expression more than one polypeptide. The vector used in the methods and compositions described herein may include a protein cleavage site that allows expression of more than one polypeptide. Examples of protein cleavage sites that allow expression of more than one polypeptide are described in Klunnp et al., Gene Ther.; 8:811 (2001), Osborn et al., Molecular Therapy 12:569 (2005), Szymczak and Vignali, Expert Opin Biol Ther. 5:627 (2005), and Szymczak et al., Nat Biotechnol. 22:589 (2004), the disclosures of which are incorporated herein by reference as they pertain to protein cleavage sites that allow expression of more than one polypeptide. It will be readily apparent to one skilled in the art that other elements that permit expression of multiple polypeptides identified in the future are useful and may be utilized in the vectors suitable for use with the compositions and methods described herein.
The vector used in the methods and compositions described herein may, be a clinical grade vector.
The viral vectors (e.g., retroviral vectors, e.g., lentiviral vectors) may include a promoter operably coupled to the transgene encoding the polypeptide or the polynucleotide encoding the RNA to control expression. The promoter may be a ubiquitous promoter. Alternatively, the promoter may be a tissue specific promoter, such as a myeloid cell-specific or hepatocyte-specific promoter. Suitable promoters that may be used with the compositions described herein include CD11 b promoter, sp146/p47 promoter, CD68 promoter, sp146/gp9 promoter, elongation factor 1 α (EF1α) promoter, EF1α short form (EFS) promoter, phosphoglycerate kinase (PGK) promoter, α-globin promoter, and β-globin promoter. Other promoters that may be used include, e.g., DC172 promoter, human serum albumin promoter, alpha1 antitrypsin promoter, thyroxine binding globulin promoter. The DC172 promoter is described in Jacob, et al. Gene Ther. 15:594-603, 2008, hereby incorporated by reference in its entirety.
The viral vectors (e.g., retroviral vectors, e.g., lentiviral vectors) may include an enhancer operably coupled to the transgene encoding the polypeptide or the polynucleotide encoding the RNA to control expression. The enhancer may include a β-globin locus control region ((3LCR).
Methods of Measuring nORF Gene Expression
Preferably, the compositions and methods of the disclosure are used to facilitate expression of a nORF at physiologically normal levels in a patient (e.g., a human patient), decrease expression of an upregulated nORF, or increase expression of a downregulated nORF. The therapeutic agents of the disclosure, for example, may reduce the dysregulated nORF expression in a human subject. For example, the therapeutic agents of the disclosure may reduce dysregulated nORF expression e.g., by about 1%, 2%, 3%, 4%, 5%, 10%, 15%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 95%, 97%, 98%, or 99%. Alternatively, the therapeutic agents of the disclosure may increase the dysregulated nORF expression in a human subject. For example, the therapeutic agents of the disclosure may increase dysregulated nORF expression, e.g., by at least 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 100%, 200%, 300%, 400%, 500%, or more.
The expression level of the nORF expressed in a patient can be ascertained, for example, by evaluating the concentration or relative abundance of mRNA transcripts derived from transcription of the nORF. Additionally, or alternatively, expression can be determined by evaluating the concentration or relative abundance of the nORF following transcription and/or translation of an inhibitor that decreases an amount of the dysregulated nORF. Protein concentrations can also be assessed using functional assays, such as MDP detection assays. Expression can be evaluated by a number of methodologies known in the art, including, but not limited to, nucleic acid sequencing, microarray analysis, proteomics, in-situ hybridization (e.g., fluorescence in-situ hybridization (FISH)), amplification-based assays, in situ hybridization, fluorescence activated cell sorting (FACS), northern analysis and/or PCR analysis of mRNAs.
Nucleic acid-based methods for determining expression (e.g., of an RNA inhibitor or an RNA encoding the nORF) detection that may be used in conjunction with the compositions and methods described herein include imaging-based techniques (e.g., Northern blotting or Southern blotting). Such techniques may be performed using cells obtained from a patient following administration of the polynucleotide encoding the agent. Northern blot analysis is a conventional technique well known in the art and is described, for example, in Molecular Cloning, a Laboratory Manual, second edition, 1989, Sambrook, Fritch, Maniatis, Cold Spring Harbor Press, 10 Skyline Drive, Plainview, NY 11803-2500. Typical protocols for evaluating the status of genes and gene products are found, for example in Ausubel et al., eds., 1995, Current Protocols In Molecular Biology, Units 2 (Northern Blotting), 4 (Southern Blotting), 15 (Immunoblotting) and 18 (PCR Analysis).
Detection techniques that may be used in conjunction with the compositions and methods described herein to evaluate nORF expression further include microarray sequencing experiments (e.g., Sanger sequencing and next-generation sequencing methods, also known as high-throughput sequencing or deep sequencing). Exemplary next generation sequencing technologies include, without limitation, Illumina sequencing, Ion Torrent sequencing, 454 sequencing, SOLiD sequencing, and nanopore sequencing platforms. Additional methods of sequencing known in the art can also be used. For instance, expression at the mRNA level may be determined using RNA-Seq (e.g., as described in Mortazavi et al., Nat. Methods 5:621-628 (2008) the disclosure of which is incorporated herein by reference in their entirety). RNA-Seq is a robust technology for monitoring expression by direct sequencing the RNA molecules in a sample. Briefly, this methodology may involve fragmentation of RNA to an average length of 200 nucleotides, conversion to cDNA by random priming, and synthesis of double-stranded cDNA (e.g., using the Just cDNA DoubleStranded cDNA Synthesis Kit from Agilent Technology). Then, the cDNA is converted into a molecular library for sequencing by addition of sequence adapters for each library (e.g., from Illumina®/Solexa), and the resulting 50-100 nucleotide reads are mapped onto the genome.
Expression levels of the nORF may be determined using microarray-based platforms (e.g., single-nucleotide polymorphism arrays), as microarray technology offers high resolution. Details of various microarray methods can be found in the literature. See, for example, U.S. Pat. No. 6,232,068 and Pollack et al., Nat. Genet. 23:41-46 (1999), the disclosures of each of which are incorporated herein by reference in their entirety. Using nucleic acid microarrays, mRNA samples are reverse transcribed and labeled to generate cDNA. The probes can then hybridize to one or more complementary nucleic acids arrayed and immobilized on a solid support. The array can be configured, for example, such that the sequence and position of each member of the array is known. Hybridization of a labeled probe with a particular array member indicates that the sample from which the probe was derived expresses that gene. Expression level may be quantified according to the amount of signal detected from hybridized probe-sample complexes. A typical microarray experiment involves the following steps: 1) preparation of fluorescently labeled target from RNA isolated from the sample, 2) hybridization of the labeled target to the microarray, 3) washing, staining, and scanning of the array, 4) analysis of the scanned image and 5) generation of gene expression profiles. One example of a microarray processor is the Affymetrix GENECHIP® system, which is commercially available and comprises arrays fabricated by direct synthesis of oligonucleotides on a glass surface. Other systems may be used as known to one skilled in the art.
Amplification-based assays also can be used to measure the expression level of the nORF or RNA in a target cell following delivery to a patient. In such assays, the nucleic acid sequences of the gene act as a template in an amplification reaction (for example, PCR, such as qPCR). In a quantitative amplification, the amount of amplification product is proportional to the amount of template in the original sample. Comparison to appropriate controls provides a measure of the expression level of the gene, corresponding to the specific probe used, according to the principles described herein. Methods of real-time qPCR using TaqMan probes are well known in the art. Detailed protocols for real-time qPCR are provided, for example, in Gibson et al., Genome Res. 6:995-1001 (1996), and in Heid et al., Genome Res. 6:986-994 (1996), the disclosures of each of which are incorporated herein by reference in their entirety. Levels of gene expression as described herein can be determined by RT-PCR technology. Probes used for PCR may be labeled with a detectable marker, such as, for example, a radioisotope, fluorescent compound, bioluminescent compound, a chemiluminescent compound, metal chelator, or enzyme.
Expression of the nORF can additionally be determined by measuring the concentration or relative abundance of a corresponding protein product (e.g., the nORF in a noncancerous cell or the dysregulated nORF). Protein levels can be assessed using standard detection techniques known in the art. Protein expression assays suitable for use with the compositions and methods described herein include proteomics approaches, immunohistochemical and/or western blot analysis, immunoprecipitation, molecular binding assays, ELISA, enzyme-linked immunofiltration assay (ELIFA), mass spectrometry, mass spectrometric immunoassay, and biochemical enzymatic activity assays. In particular, proteomics methods can be used to generate large-scale protein expression datasets in multiplex. Proteomics methods may utilize mass spectrometry to detect and quantify polypeptides (e.g., proteins) and/or peptide microarrays utilizing capture reagents (e.g., antibodies) specific to a panel of target proteins to identify and measure expression levels of proteins expressed in a sample (e.g., a single cell sample or a multi-cell population).
Exemplary peptide microarrays have a substrate-bound plurality of polypeptides, the binding of an oligonucleotide, a peptide, or a protein to each of the plurality of bound polypeptides being separately detectable. Alternatively, the peptide microarray may include a plurality of binders, including, but not limited to, monoclonal antibodies, polyclonal antibodies, phage display binders, yeast two-hybrid binders, aptamers, which can specifically detect the binding of specific oligonucleotides, peptides, or proteins. Examples of peptide arrays may be found in U.S. Pat. Nos. 6,268,210, 5,766,960, and 5,143,854, the disclosures of each of which are incorporated herein by reference in their entirety.
Mass spectrometry (MS) may be used in conjunction with the methods described herein to identify and characterize expression of the nORF in a cell from a patient (e.g., a human patient) following delivery of the transgene encoding the nORF. Any method of MS known in the art may be used to determine, detect, and/or measure a protein or peptide fragment of interest, e.g., LC-MS, ESI-MS, ESI-MS/MS, MALDI-TOF-MS, MALDI-TOF/TOF-MS, tandem MS, and the like. Mass spectrometers generally contain an ion source and optics, mass analyzer, and data processing electronics. Mass analyzers include scanning and ion-beam mass spectrometers, such as time-of-flight (TOF) and quadruple (Q), and trapping mass spectrometers, such as ion trap (IT), Orbitrap, and Fourier transform ion cyclotron resonance (FT-ICR), may be used in the methods described herein. Details of various MS methods can be found in the literature. See, for example, Yates et al., Annu. Rev. Biomed. Eng. 11:49-79, 2009, the disclosure of which is incorporated herein by reference in its entirety.
Prior to MS analysis, proteins in a sample obtained from the patient can be first digested into smaller peptides by chemical (e.g., via cyanogen bromide cleavage) or enzymatic (e.g., trypsin) digestion. Complex peptide samples also benefit from the use of front-end separation techniques, e.g., 2D-PAGE, HPLC, RPLC, and affinity chromatography. The digested, and optionally separated, sample is then ionized using an ion source to create charged molecules for further analysis. Ionization of the sample may be performed, e.g., by electrospray ionization (ESI), atmospheric pressure chemical ionization (APCI), photoionization, electron ionization, fast atom bombardment (FAB)/liquid secondary ionization (LSIMS), matrix assisted laser desorption/ionization (MALDI), field ionization, field desorption, thermospray/plasmaspray ionization, and particle beam ionization. Additional information relating to the choice of ionization method is known to those of skill in the art.
After ionization, digested peptides may then be fragmented to generate signature MS/MS spectra. Tandem MS, also known as MS/MS, may be particularly useful for analyzing complex mixtures. Tandem MS involves multiple steps of MS selection, with some form of ion fragmentation occurring in between the stages, which may be accomplished with individual mass spectrometer elements separated in space or using a single mass spectrometer with the MS steps separated in time. In spatially separated tandem MS, the elements are physically separated and distinct, with a physical connection between the elements to maintain high vacuum. In temporally separated tandem MS, separation is accomplished with ions trapped in the same place, with multiple separation steps taking place over time. Signature MS/MS spectra may then be compared against a peptide sequence database (e.g., SEQUEST). Post-translational modifications to peptides may also be determined, for example, by searching spectra against a database while allowing for specific peptide modifications.
A number of cancers are known in the art that are contemplated in conjunction with the methods described herein. The present invention contemplates treatment of a cancer in which a nORF exhibits increased or decreased expression, e.g., relative to a noncancerous cell.
The method may reduce the size (e.g., by 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 95%, 97%, or 99%) of a tumor (e.g., a breast tumor). The method may decrease or slow (e.g., by at least 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 95%, 97%, or 99%) the progression of cancer. The method may decrease (e.g., by at least 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 95%, 97%, or 99%) the risk of developing cancer. The method may decrease (e.g., by at least 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 95%, 97%, or 99%) the risk of developing cancer.
In some embodiments, the cancer is selected from the list consisting of breast invasive carcinoma, colon adenocarcinoma, esophageal carcinoma, head and neck squamous cell carcinoma, kidney chromophobe, kidney clear cell carcinoma, kidney papillary cell carcinoma, liver hepatocellular carcinoma, lung adenocarcinoma, lung squamous cell carcinoma, prostrate adenocarcinoma, stomach adenocarcinoma, thyroid carcinoma, and uterine corpus endometrioid carcinoma.
In some embodiments, the nORF is selected from Table 1.
In some embodiments, the nORF is selected from Table 2.
In some embodiments, the nORF is selected from Table 3.
In some embodiments, the nORF is selected from Table 4.
In some embodiments, the nORF is selected from Table 5.
In some embodiments, the nORF is not HOXB-AS3.
In some embodiments, the cancer is not colorectal cancer.
In some embodiments, the nORF is not PINT87aa (LINC-PINT).
In some embodiments, the cancer is not glioblastoma.
The following examples further illustrate the invention but should not be construed as in any way limiting its scope.
nORFs are Pervasively Translated and Important for Further Investigation
nORFs are typically smaller than canonical ORFs, the peptides or micro-proteins they encode are particularly attractive as putative allosteric cellular regulators, due to their size and the potential specificity of peptide interactions. Therefore, because the accepted nomenclature itself is inconsistent, we classified and catalogued all human nORFs from various sources, prioritizing those with strong evidence for translation and distinguishing between nORFs that are in frame and out of frame with overlapping canonical ORFs and released it as an open-source database (norfs.org/home).
Identifying and Characterizing Transcripts Encoding nORFs
To identify transcripts encoding nORFs (nORF transcripts), we extracted genomic coordinates of transcripts quantified in the UCSC Toil pipeline from the GENCODE v23 reference genome annotation and compared these with the genomic coordinates of nORFs acquired from the curated nORFs.org database, using a custom pipeline (
To determine if nORF transcripts are expressed in any tissue included in the study, we defined an expression threshold of 0.5 counts per million (CPM) across at least 10% of a single tissue. This allowed us to prioritize transcripts that are more likely to be accurately quantified and expressed at a biologically meaningful level. Using this threshold, we identified 926 expressed nORF transcripts for inclusion in this study.
We characterized the genomic properties of all nORF transcripts (
We considered genomic distribution and strand bias (
Following identification of nORF transcripts, we evaluated transcript mean expression across all GTEx normal tissues included in this study. We showed mean nORF transcript expression compared with canonical protein-coding transcripts and also compared against canonical antisense and lincRNA expression—as these are the two main transcript classifications within which nORF transcripts are identified (
Many nORF transcripts were poorly expressed, with mean CPM values below 0.5. We identified and prioritized nORF transcripts frequently expressed in cancer tissues or the corresponding NAT or GTEx normal tissue. Both cancer and reference normal tissues were considered when identifying frequently expressed nORF transcripts, as we aimed to capture nORF transcripts both up- and down-regulated between cancer and normal tissues. Frequently expressed nORF transcripts were defined as having CPM greater than 0.5 across at least 70% of samples in either cancer or corresponding reference tissue. A representative distribution of expression across samples in cancer tissue and corresponding NAT (
When comparing cancer with NAT, we determined 359 out of 926 nORF transcripts were frequently expressed in at least one cancer type; when comparing with GTEx normal tissue, 464 out of 926 nORF transcripts were frequently expressed in at least one cancer type. The number of frequently expressed nORF transcripts identified was consistent across cancer types (
A large proportion of nORF transcripts were frequently expressed across all cancer types—109 nORF transcripts for cancer and NAT; 115 nORF transcripts for cancer and GTEx normal tissue. On the other hand, comparatively few nORF transcripts were frequently expressed in any particular subset of cancer types—for example, just 14 nORF transcripts were only frequently expressed in thyroid carcinoma or thyroid NAT. This likely reflects consistent expression of nORF transcripts across tissues. A disproportionate number of nORF transcripts (79) are frequently expressed only in testicular germ cell tumor tissue or GTEx testis tissue, which is consistent with mean transcript expression patterns in testis tissue (
Identifying Differentially Expressed nORF Transcripts
To identify nORF transcripts dysregulated in cancer, we performed differential expression analysis for cancer compared with either NAT or GTEx normal tissue. We normalized RNA-Seq expected counts from the UCSC Toil dataset using the trimmed mean of M-values (TMM) method and performed differential expression analysis using the general linear model (GLM) framework provided by edgeR, as described in Materials and Methods. A fold change threshold of 2 and adjusted p value threshold of 0.001 were used to call differentially expressed nORF transcripts. Only frequently expressed nORF transcripts were considered. Corresponding analysis using a fold change threshold of 1.5 is provided in
This analysis revealed 152 nORF transcripts as dysregulated in at least a single cancer type when comparing cancer with NAT (
We observed a limited number of nORF transcripts with cancer-type specific dysregulation. In lung squamous cell carcinoma 13 nORF transcripts were uniquely upregulated, and 10 uniquely down-regulated, when compared against NAT. Kidney clear cell carcinoma, kidney chromophobe and testicular germ cell tumors also exhibited a large degree of cancer-type specific dysregulation (
To assess the reproducibility of differential expression results when comparing against NAT or GTEx normal tissue, we investigated differentially expressed nORF transcripts identified in eight cancer types with both types of reference normal tissue. Differential expression relative to GTEx normal tissue consistently revealed a larger number of dysregulated nORF transcripts. Most cancer types showed highly reproducible differential expression results between the two reference normal tissues (
We have shown that nORF transcripts are frequently expressed across multiple cancer types and reference normal tissues, and that many of these nORF transcripts are transcriptionally dysregulated in cancers. To determine whether any differentially expressed nORF transcripts can be used as prognostic marker, we investigated the relationship between nORF transcript expression and overall patient survival, for nORF transcripts differentially expressed between cancers and NAT. We used survival data for TOGA cohorts provided by the UCSC Toil Recompute Compendium and divided each cohort into high and low expression groups for each nORF transcript, as detailed in Materials and Methods. We identified 43 nORF transcripts where expression was significantly associated with patient overall survival in at least one of the 12 cancer types included in this survival analysis, with an adjusted p value threshold of 0.05 (
We investigated further nORF transcripts reproducibly differentially expressed both compared with NAT and GTEx normal tissue. For a subset of 33 nORF transcripts: (i) the transcript is reproducibly differentially expressed in cancer compared with NAT and GTEx normal tissue (ii) transcript expression is associated with prognosis (adjusted p<0.05) (iii) and transcripts up-regulated in cancer are associated with poor prognosis, and vice versa. Kaplan Meier survival curves are shown for the nORF transcripts most significantly associated with prognosis, in Kidney Clear Cell Carcinoma (
Through comprehensive analysis of RNA-Seq data from 22 cancer types, we have identified transcripts containing novel open reading frames and demonstrated that many nORF transcripts are frequently expressed in multiple cancers. Additionally, we have shown that many of these nORF transcripts are differentially expressed between cancer and normal tissue, and some of these nORF transcripts are uniquely differentially expressed in specific cancer types. Furthermore, we have shown that expression of some differentially expressed nORF transcripts have prognostic value—this is particularly convincing for four nORF transcripts reproducibly and uniquely identified as up-regulated in either liver hepatocellular carcinoma or lung adenocarcinoma, for which high expression was associated with poor prognosis.
TOGA and GTEx RNA-Seq and survival data was downloaded from the TCGA TARGET GTEx′ cohort of the UCSC Toil Recompute Compendium. Transcriptome alignment had been performed using STAR (GRCh38) and transcript expression quantified using RSEM, using transcripts present in the GENCODE v23 genome annotation. Transcript-level RSEM expected counts, TOGA survival data and phenotype data were obtained. The GENCODE v23 and corresponding Ensembl v81 genome annotations were downloaded, and transcript and coding sequence properties were extracted from the annotation files using a custom script. RSEM expected counts provided by the UCSC Toil Recompute Compendium were log 2(expected_count+1) transformed, and this transformation was removed to produce raw expected counts for use in this analysis. All data processing was performed using R, R Studio, the R package Tidyverse and unix command line tools. The Ensembl genome annotation was processed in R using ensembl db, and genomic coordinates were processed using GenomicRanges. Set diagrams were produced using UpSetR.
Mappings of TOGA cancer tissue samples to normal adjacent tissue (NAT) and GTEx normal tissue were extracted from the phenotype data provided by the UCSC Toil Recompute Compendium. We included solid tumor TOGA cancer tissues with at least 50 samples, with matched NAT or GTEx normal tissue with at least 10 or 50 samples respectively—a less stringent threshold for inclusion was used for NAT because these samples were less abundant. RSEM expected count data was filtered to retain only selected samples and expressed transcripts prior to normalization and differential expression analysis. A single sample containing missing expected count values was excluded from this analysis.
Prior to library size normalization and differential expression analysis, transcripts with poor expression were excluded from analysis. Applying a CPM threshold to identify expressed transcripts prior to TMM normalization and differential expression analysis has been shown to improve false discovery rate and is recommended practice for edgeR. Expected counts were transformed to CPM and transcripts are classified as expressed if they had expected count greater than 0.5 CPM in at least 10% of the samples of a single cancer or normal tissue. Expressed transcripts are retained. Best practices for setting thresholds for transcript-level expression are poorly established, and the thresholds used in this study were, whilst informed by the literature, largely arbitrary.
To characterize the expression of transcripts encoding nORFs across multiple cancer types and corresponding normal tissues, we obtained transcript-level RNA-Seq expression data from the UCSC Toil Recompute Compendium. This dataset includes 11,194 cancer and normal adjacent tissue samples (NAT) from TCGA and 8,003 normal tissue samples from GTEx. We used metadata provided by the UCSC Toil Recompute Compendium to match cancer, NAT and GTEx normal tissues and determine the number of samples available for each tissue. To ensure consistent and reliable results, we included solid tumor TCGA cancer tissues with at least 50 samples, with matched NAT or GTEx normal tissue containing at least 10 or 50 samples respectively—a less stringent threshold for inclusion was used for NAT because these samples are less abundant. This resulted in a total of 7,885 samples across 22 cancer types from TCGA, together with 677 NAT samples and 4,010 GTEx normal samples.
NAT and GTEx normal tissues provide non-redundant reference tissues. NAT samples closely resemble cancer samples both as a result of reduced variation in patient differences and sample processing. However, NAT is affected by changes in the tumor microenvironment and samples are less abundant than GTEx normal tissue samples. Seven cancer tissues included in this study are matched to both NAT and GTEx normal tissue which allowed us to determine whether differential expression results are reproducible across different reference tissues.
Genomic coordinates of nORFs with experimental evidence for translation were obtained from the nORFs.org database (norfs.org/home). Transcript genomic coordinates were obtained from the GENCODE v23 reference annotation. GffCompare was used to identify open reading frames and transcripts with completely matching intron chains. GffCompare performs stringent filtering to detect and remove redundant input transcripts, and this deduplication is described in detail in the documentation. Specifically, to achieve stringent deduplication of nORFs, GffCompare was run with nORF coordinates as the ‘reference set’ and transcript coordinates as the ‘query set’, with default parameters. The resultant ‘.refmap’ file containing information on overlaps between nORF and transcript coordinates was processed in R and annotated. nORF-transcript mappings identified by GffCompare were filtered to retain only those with a complete intron chain match, and for which the genomic coordinates of the nORF were completely contained within the transcript. nORFs present in multiple transcripts were excluded. Transcript biotypes were extracted from the GENCODE annotation file and open reading frames contained in protein-coding transcripts (transcripts with biotype: “protein_coding”, “IG_C_gene”, “IG_D_gene”, “IG_J_gene”, “IG_V_gene”, “TR_C_gene”, “TR_D_gene”, “TR_J_gene”, “TR_V_gene”) and rRNA transcripts were excluded. Novel and canonical ORF lengths were determined using ensembldb.
Normalization and differential expression were performed separately for comparison of cancer tissue with NAT and with GTEx normal tissue. RNA-Seq expected counts were normalized across samples using the trimmed mean of M-values (TMM) method to normalize for read depth and composition. As comparisons in differential expression were not made across transcripts, no normalization was introduced for effective transcript length.
To identify frequently expressed transcripts, CPM values were calculated across all expressed transcripts following TMM normalization using edgeR. Transcripts were classed as frequently expressed if they had CPM greater than 0.5 in at least 70% of the samples in the normal or cancer tissue of interest.
Transcript differential expression was performed using all expressed transcripts to provide correct significance testing and improve reliability of dispersion estimation. The R package edgeR was used to perform differential expression analysis using a general linear model framework—this package was chosen as it is (i) highly cited (ii) suitable for transcript-level analysis (iii) compatible with non-integer expected counts from RSEM (iv) and exhibits fast performance on large datasets. A simple additive model with no intercept was constructed, with normal reference tissues and cancer tissues each represented by a single coefficient. The process used for differential expression analysis is detailed in the edgeR manual. Briefly, transcript-wise dispersions were estimated under the general linear model framework using the Cox-Reid profile-adjusted likelihood method, which takes into account multiple factors by fitting the described model. A negative binomial model was fitted for each transcript, and thresholded hypotheses were tested to provide meaningful p values and reliable control of false discovery rate. A fold change threshold of 1.5 or 2 was used to identify differentially expressed transcripts, with an adjusted p value threshold of 0.001. Coefficients representing cancer tissues and their corresponding normal reference tissues were compared under this framework. The Benjamini and Hochberg method was used to adjust p values for multiple testing and control false discovery rate.
Overall survival (OS) analysis was performed using the R packages survival and survminer. nORF transcripts are included in survival analysis if they were differentially expressed in the cancer type of interest compared with NAT and were expressed at greater than 0.5 CPM in at least 70% of the samples in the cancer tissue cohort. For each cancer type and for the nORF transcript considered, the cohort was split into high and low expression groups. Groups were selected which were best segregated based on overall survival, using the Maximally Selected Rank Statistic, with at least 30% of patients assigned to each expression group to avoid forming groups with a small number of patients. Kaplan Meier curves were generated, and curves were compared using a log-rank test. The Benjamini and Hochberg method was used to adjust p values for multiple testing and control false discovery rate. A Cox proportional hazards regression model was fitted to overall survival data and hazard ratios were derived from the model coefficients. Both the Kaplan Meier and Cox proportional hazards regression models assume proportional hazards, where the hazard ratio between the high and low expression group remains constant over time.
While the invention has been described in connection with specific embodiments thereof, it will be understood that it is capable of further modifications and this application is intended to cover any variations, uses, or adaptations of the invention following, in general, the principles of the invention and including such departures from the invention that come within known or customary practice within the art to which the invention pertains and may be applied to the essential features hereinbefore set forth, and follows in the scope of the claims.
Other embodiments are within the claims.
This application claims the benefit of U.S. Provisional Application No. 63/126,309 filed on Dec. 16, 2020, which is incorporated herein by reference in its entirety.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/IB2021/061801 | 12/15/2021 | WO |
Number | Date | Country | |
---|---|---|---|
63126309 | Dec 2020 | US |