This article has Open Peer Review reports available.
Gene expression in BMPR2 mutation carriers with and without evidence of Pulmonary Arterial Hypertension suggests pathways relevant to disease penetrance
© West et al; licensee BioMed Central Ltd. 2008
Received: 07 March 2008
Accepted: 29 September 2008
Published: 29 September 2008
While BMPR2 mutation strongly predisposes to pulmonary arterial hypertension (PAH), only 20% of mutation carriers develop clinical disease. This finding suggests that modifier genes contribute to FPAH clinical expression. Since modifiers are likely to be common alleles, this problem is not tractable by traditional genetic approaches. Furthermore, examination of gene expression is complicated by confounding effects attributable to drugs and the disease process itself.
To resolve these problems, B-cells were isolated, EBV-immortalized, and cultured from familial PAH patients with BMPR2 mutations, mutation positive but disease-free family members, and family members without mutation. This allows examination of differences in gene expression without drug or disease-related effects. These differences were assayed by Affymetrix array, with follow-up by quantitative RT-PCR and additional statistical analyses.
By gene array, we found consistent alterations in multiple pathways with known relationship to PAH, including actin organization, immune function, calcium balance, growth, and apoptosis. Selected genes were verified by quantitative RT-PCR using a larger sample set. One of these, CYP1B1, had tenfold lower expression than control groups in female but not male PAH patients. Analysis of overrepresented gene ontology groups suggests that risk of disease correlates with alterations in pathways more strongly than with any specific gene within those pathways.
Disease status in BMPR2 mutation carriers was correlated with alterations in proliferation, GTP signaling, and stress response pathway expression. The estrogen metabolizing gene CYP1B1 is a strong candidate as a modifier gene in female PAH patients.
Pulmonary arterial hypertension (PAH) is a lethal disorder characterized by pulmonary vascular constriction and remodeling leading to progressively worsening right ventricular hypertrophy, and eventually right heart failure. The familial form (FPAH) is usually caused by mutations in the type 2 receptor for the BMP pathway, BMPR2 [1–3].
The lifetime risk of developing PAH in a BMPR2 mutation carrier is less than 20%. This suggests the need for environmental or genetic modifiers for full expression of the disease. Understanding modifier genes would both clarify the molecular etiology of the disease, which is obscured by the myriad of BMPR2 functions, and allow the FPAH risk to be refined for asymptomatic individuals that carry BMPR2 mutations. The search for these modifiers, which has been ongoing for several years, is complicated by many factors[5, 6]. For example, since the modifier alleles are likely to be common variants, traditional genetic approaches are problematic, since the same allele could be associated with more than one line of inheritance within a family tree. This is further complicated by the possibility that predisposition to disease within BMPR2 mutation carriers is caused by a confluence of factors rather than a single modifier gene, even within a single family.
Attempts to find contributing genes through gene array approaches have been tried, using lung tissue from transplant patients and circulating cells[8, 9]. Interpretation of those data is made difficult by several factors. First, the majority of genes dysregulated are likely to be caused by the presence of end-stage disease. Second, drug effects cannot properly be controlled, and can be quite large. Finally, there is a large background in gene expression differences caused by the diverse background genetics of the patients. Finding direct evidence of modifier gene or even mutation effects within all of these sources of experimental noise is essentially impossible.
In order to overcome these issues, in this study we have used patient-derived lymphoblastoid cell lines to hunt for modifier genes. Lymphoblastoid cell lines, made by EBV-immortalizing B cells, resolve many of these issues. While culturing the cells introduces its own alterations in gene expression, these alterations are uniform across samples. Lymphoblastoid lines have been used successfully in gene expression studies for disease processes as diverse as schizophrenia, drug resistance, autism, and asthma[10–13]. Moreover, culturing the cells removes them from both the disease milieu and drug effect. While other types of cells may have changes in differentiation caused by the disease state, this is not true of B cells, whose lineage commitment is considered to be unidirectional and irreversible under physiologic conditions[14, 15]. Finally, in order to minimize variation caused by genetic background we derived all of our lines from within the same extended family. An additional advantage is that we are directly measuring differences in baseline gene expression, without regard for either the polymorphism that caused them or their hereditary origin.
In this study, we use Affymetrix arrays to compare gene expression in BMPR2 mutation carriers that are either asymptomatic carriers (unaffected) or who have PAH (affected). We confirmed selected differentially expressed genes using quantitative RT-PCR on lines derived from a larger number of patients as well as from family members without mutation (non-carriers).
We found several broad pathways with differential expression between affected and unaffected BMPR2 mutation carriers, including stress response, actin organization/g-protein, calcium balance, and cell-cycle related genes. Analysis of overrepresented gene ontology groups suggests that it is pathway-specific, not gene-specific changes that are associated with increased risk.
Age at Draw
Age at Diagnosis
Asymptomatic status was confirmed in each individual by echocardiography at the time of blood draw. Only unaffected carriers over the age of 60 were used for gene array experiments, since they had the least probability of later developing later disease, and would still have substantial difference in age of onset if they did harbor occult disease. Unaffected carriers are on average second degree relatives to each-other; they were on average sixth degree relation to FPAH patients. FPAH patients were on average eighth degree relations to each-other (they are 3rd cousins, 3rd cousins once removed, or 2nd cousins once removed). The close relation in the unaffected carriers is coincidence; there are unaffected carriers in many branches of the family tree, but not that were currently available for blood draw and in the correct age range.
An additional 24 blood samples were collected from 3 additional families for confirmatory studies related to CYP1B1 expression. These include a family with an exon 4–5 deletion (3 non-carriers, 3 affected and 2 unaffected mutation carriers), a family with an exon 9 deletion (1 non-carrier, 3 affected and 3 unaffected mutation carriers), and a family with an exon 9 frameshift mutation (2 non-carriers, 2 affected and 5 unaffected mutation carriers).
The study was approved by the institutional review board at Vanderbilt University Medical Center, and written, informed consent was obtained from all subjects included in the study. Unique identifiers to conceal identity were assigned to the samples before their receipt in the laboratory.
Lymphocytes were isolated from anticoagulated whole blood within 48 hrs of collection and exposed to Epstein-Barr Virus (EBV) to induce cell immortalization. Two ml blood was diluted with 2 ml PBS, layered on top of 3 ml of Lympho Separation Medium (MP Biomedicals) and centrifuged for 10 minutes at 1,000 × g at room temperature. Using a Pasteur pipet, the lymphocytes were removed from the serum/Lympho Sep Media interface, washed in 10 ml PBS and then resuspended in 3 ml lymphoblast media (RPMI 1640 media containing L-glutamine, and 20% fetal bovine serum) containing 2 μg/ml cyclosporine. The lymphocytes were then infected with 3 ml Epstein-Barr virus (EBV) and transferred to a T-25 vent capped flask. The cells were incubated at 37°C/5% CO2 and fed weekly with lymphoblast media + cyclosporine until signs of growth occurred.
RNA was isolated from lymphocytes using a Qiagen RNeasy mini kit (Valencia, CA). First and second strand complimentary DNA was synthesized using standard techniques. Biotin-labeled antisense complimentary RNA was produced by an in vitro transcription reaction. Human Genome U133 Plus 2.0 microarrays (Affymetrix, Foster City, CA) were hybridized with 20 μg cRNA. Target hybridization, washing, staining, and scanning probe arrays were done following an Affymetrix GeneChip Expression Analysis Manual. All array results have been submitted to the NCBI gene expression and hybridization array data repository (GEO, http://www.ncbi.nlm.nih.gov/geo/), as series GSE10767.
Affymetrix Cel files were loaded into dChip array analysis software. Overall signal strength from arrays was normalized to the median array, and expression levels determined using the perfect match/mismatch (PM/MM) algorithm. Differentially expressed genes were determined using a 95% probability of a minimum 1.4× change and a minimum absolute difference of 150 (arbitrary units). Using random reassignment of group identity, this produced a median 17.2% false discovery rate, a reasonable compromise between sensitivity and specificity. These requirements in themselves result in a list of genes in which there is no overlap in expression between groups. For every gene in the list, every affected mutation carrier has higher expression than every unaffected mutation carrier, or every affected mutation carrier has lower expression than every unaffected mutation carrier. Every gene thus has a p < .05 for significance by Wilcoxon rank-sum test.
Gene ontology was determined using the Classify Genes tool within dChip, with gene ontology files downloaded from the Gene Ontology Consortium http://www.geneontology.org, and were grouped for the purposes of this study by biological process (rather than molecular function or cellulcar component) classification[19, 20]. Genes which lacked biological process annotation in the database were assigned to a group through a brief literature review. Most genes fall into several gene ontology groups. Thus, the selection of group for each gene was somewhat arbitrary. However, in the case of the genes dysregulated in this study, different choice of group would primarily shift genes between groups rather than create new groups. For instance, matrix genes can be involved in vascular contractility, stress response, and growth. Actin organization and G-proteins can be involved in both contractility and lymphocyte recruitment. Thus, the groups of genes identified are for the most part tightly interrelated in function, with functions directly related to the etiology of PAH.
GGA TGC CTC TCT TGC TCT G
GTC TTC CCC TCC ATC GTG
AAC GTA CCG GCC ACT ATC AC
ACG ACC TGA TCC AAT TCT GC
CCA AGA GCA AGT GGA GAA GC
AGG CAT CTG AGG TGA ACA GG
TAT TCG ATG TCA AGC GAA CG
ATA TTT CCG GGT TGG AGA CC
GAG AGC TGG CCA AGA TGA AG
GCA CTC AAG GTA GCC AAA GG
CAT GAG TTC CAT GGT GAA CG
GCT CAA ATT TGG CCT GTA GC
AAT GTT CAC GCG CCT TAT TC
AAC CAG CGA TAT TTG GAT GC
Confidence intervals for fold changes in array analysis are determined by algorithms internal to dChip, as previously described. These functionally result in this case in genes for which difference between affected and unaffected groups is also p < .05 by Wilcoxon signed-rank test.
Other statistical analyses were performed using the JMP program (SAS, Cary, NC). Comparisons for quantitative PCRs comparing three groups (non-carrier, unaffected, affected) were performed using Kruskal-Wallis, although results are similar using ANOVA. Analysis for sex-specific differences in carriers was performed using two-way ANOVA on log transformed values, with comparisons between individual values by Tukey's HSD post-hoc. Analyses of overrepresentation of gene ontology groups was performed by fisher's exact test, and overrepresentation of genes within those groups were performed using the one sample z test. While a substantial body of literature is arising suggesting alternate and potentially more powerful methods of determining overrepresentation of gene ontology groups, the field still seems to be in flux and so we have continued to rely on Fisher's exact test for this purpose.
BMPR2 mutation carriers with disease show alteration in multiple PAH-related pathways compared to unaffected carriers
RNA derived from lymphoblastoid lines from four affected and three unaffected BMPR2 mutation carriers (marked with * on Table 1) were used to probe Affymetrix U133 Plus 2.0 gene expression arrays. These included every lymphoblastoid line from affected carriers within the family available at the time (one additional line was derived later), and every line from unaffected carriers over the age of sixty (who were thus most likely to remain unaffected).
Quantitative RT-PCR Confirms Gene Expression Data
To determine whether decreased CYP1B1 expression was associated with disease status in other FPAH, lymphoblastoid cells were derived from 24 individuals from three distinct patient families, each with too few patients for individual analysis. We found that log-transformed CYP1B1 expression values were significantly decreased (p < .05 by Kruskal-Wallis) in disease-affected mutation carriers compared with unaffected carriers or non-carriers, with a trend towards a sex-specific effect (not shown). When combined with the data from our original family, numbers became sufficient to show a clear sex-specific effect (p < .01 by two-way ANOVA), with CYP1B1 expression levels almost tenfold decreased in female, but not male, affected mutation carriers compared with unaffected carriers or non carriers (Figure 2B).
Analysis of Overrepresented Gene Ontology Groups Suggests Pathways, not Genes, as Significant Modifiers
There are two explanations for decreasing statistical significance of specific genes with increased samples; the first is that the results are spurious, the second, our hypothesis, was that they are representative of a group, rather than important in themselves.
GTP-related genes are overrepresented in genes changed in affected BMPR2 Mutation Carriers
We found that the number of genes that fell into the GTP-related ontology group was highly significantly overrepresented in each case. If the genes were randomly distributed across ontology groups, we would expect 4.9% of those with altered expression to fall into the GTP-related group (Table 3, second line). Instead, in each case, more than twice as many fell into this group (Table 3, third line), which was highly significant by one-sample z-statistics, with a p < .001 in each case.
Overlap between GTP-related genes changed in different patients
The purpose of this study was to identify genes which might predispose BMPR2 mutation carriers to FPAH by comparing carriers unaffected by disease to carriers who developed FPAH. We used cultured B-cell lines in order to remove both disease and drug effects, so that we could examine underlying differences in gene expression.
The strongest individual gene identified by our methodology was Cytochrome P450 1B1 (CYP1B1), in which low expression was associated with disease, but only in women. This gene qualifies as a plausible modifier gene for several reasons. Although we assayed B-cells, CYP1B1 is highly expressed in lung, likely in endothelial cells, and is recognized as a modifier gene for cancers [24–27]. CYP1B1 metabolizes environmental toxins, so it would suggest a plausible mechanism of a gene-environment interaction as a modifier. CYP1B1 also breaks down estrogen, reducing local concentration (systemic estrogen is regulated by other systems) and environmental estrogens are present in modern diets and medications. Since women have a 3:1 overrepresentation in PAH patients, starting at puberty, estrogen is the strongest modifier known for penetrance. Thus, lowered levels of CYP1B1 might result in increased local concentrations of estrogen, further increasing the risk of PAH; that this appears to be a sex-specific modifier gene supports this hypothesis. Follow-up to demonstrate correlation with genetic polymorphisms, thus confirming that this is cause not effect, and a finding of functional consequences in patients are necessary to fully interpret the meaning of this finding.
Our broader finding was that there were consistent differences in regulation of PAH-related pathways between BMPR2 mutation carriers affected and unaffected by disease. These changes are likely related to genetically determined differences in expression, rather than adaptive to disease state or drug effect, because the cells have been cultured to remove them from the immediate environment. Further, B-cells are not believed to unalterably change differentiation state in response to environmental effects. A previous study which examined B-cells freshly derived from patient blood in controls vs IPAH had substantially different findings consisting purely of upregulation of stress-response, suggesting that they were looking at disease and drug effects, and bolstering our belief that disease effects are not a feature of the current study.
Pathways changed included stress-response, actin organization, ras-related, g-protein, calcium balance, and proliferation-related pathways (Figure 2). All of these pathways have been previously seen as altered in PAH, but the current study is one of the first to present data suggesting that they are cause rather than effect in human patients[9, 31, 32]. A detailed analysis of the GTP-related gene ontology group suggested that the important change was in the pathway itself, which could be mediated through different subsets of the genes in the pathway (Figure 3, Tables 3, 4, Additional File 2). This implies that the search for modifier genes may involve either activities or pathways, rather than specific genes. We have seen this effect before; different strains of mice subjected to chronic hypoxia use different genes in the same pathway to achieve similar phenotypic ends.
These data are in agreement with our earlier examination of lymphoblastoid cells, in which our strongest candidate was the Ras pathway gene GRB2. Some of the other alterations found in our earlier examination of protein changes in lymphoblastoid cells are also present in this study, but were below the cutoff threshold we selected. For instance, PCMT1, which was increased 1.2× by protein array, was also increased 1.2× by Affymetrix array. Other genes did not have matching changes, likely representing issues of sensitivity and post-transcriptional regulation.
We found more than 80 genes consistently altered between affected and unaffected BMPR2 mutation carriers. It is not probable that there are this many modifier genes; rather, this implies that these genes are downstream of a much smaller number of modifier genes that are having consistent downstream effects. We have already published serotonin transporter(5HTT) and TGF-β promoter polymorphisms as modifier genes contributing to clinical expression in familial PAH, and the types of changes seen are plausible effects for either of these[35, 36]. However, we determined that these could not be upstream of the changes seen here, as two of the affected patients had the LL 5HTT promoter and two had the SL promoter. Correlation in gene expression even among altered genes between the two with the LL allele was no stronger than between the LL allele and SL allele (0.82 vs 0.88; correlation between these and unaffecteds averages 0.65 for comparison). Further, these cells do not appear to have a functional TGF-β pathway. Levels of both type 1 and type 2 receptors were below the threshold of detection on the arrays, and both canonical TGF-β target Pai1 expression levels and TGF-β luciferase-reporter response were very low (not shown). These data imply that there is a modifier gene upstream of the pathway changes seen, but that it is not a gene currently under consideration.
Our study design carried several inherent limitations. In order to obtain the strongest possible level of discrimination for modifier genes, we compared mutation carriers who developed disease early to mutation carriers still unaffected in old age (Table 1). This necessarily results in a mismatch of age between groups as part of the study design. However, we have previously shown that difference in age does not appear to create a difference in B-cell protein expression. Further, there was a difference in the degree of relation between unaffected and affected carriers; this could result in an underestimate of the false discovery rate. Another limitation lies in differences between B-cells and disease effector cells, including pulmonary smooth muscle and endothelium. For instance, in B-cells, excision/repair (Additional File 1, Figure 1) is tightly linked to proliferation, and so alterations in these genes are probably not meaningful to other cells. Conversely, B-cells do not express some genes likely of interest in disease, and so these cannot be studied in this system. For instance, the TGF-beta system does not appear to be functional in B-cells. However, the ideal tissue comparison, between lung tissue from BMPR2 mutation carriers that will never develop disease and those that will develop disease but have not yet, is inherently impossible; neither the information nor the tissue to conduct this study could be obtained. Finally, the act of EBV-transforming these cells is likely to have produced some changes. However, these changes should be consistent across samples; at worst, they could mask some effects. Concern about this issue explains our focus on GTP-related genes, rather than either stress-response or proliferation, although there is no particular reason to believe that these changes are not also valid.
The work presented here provides evidence that predisposition to disease in BMPR2 mutation carriers lies in alterations of pathways related to actin organization, stress response, and proliferation, likely achievable through alterations in any of several upstream modifiers, rather than primarily to specific individual genes within those pathways. The strongest single candidate modifier gene identified was CYP1B1, with tenfold lower expression on average in female BMPR2 mutation carriers than in unaffected carriers or non-carriers.
This study was funded by NIH P01 HL 72058.
- Cogan JD, Pauciulo MW, Batchman AP, Prince MA, Robbins IM, Hedges LK, Stanton KC, Wheeler LA, Phillips JA, Loyd JE, et al: High frequency of BMPR2 exonic deletions/duplications in familial pulmonary arterial hypertension. Am J Respir Crit Care Med. 2006, 174 (5): 590-598.View ArticlePubMedPubMed CentralGoogle Scholar
- Lane KB, Machado RD, Pauciulo MW, Thomson JR, Phillips JA, Loyd JE, Nichols WC, Trembath RC: Heterozygous germline mutations in BMPR2, encoding a TGF-beta receptor, cause familial primary pulmonary hypertension. The International PPH Consortium. Nat Genet. 2000, 26 (1): 81-84.View ArticlePubMedGoogle Scholar
- Deng Z, Morse JH, Slager SL, Cuervo N, Moore KJ, Venetos G, Kalachikov S, Cayanis E, Fischer SG, Barst RJ, et al: Familial primary pulmonary hypertension (gene PPH1) is caused by mutations in the bone morphogenetic protein receptor-II gene. Am J Hum Genet. 2000, 67 (3): 737-744.View ArticlePubMedPubMed CentralGoogle Scholar
- Newman JH, Trembath RC, Morse JA, Grunig E, Loyd JE, Adnot S, Coccolo F, Ventura C, Phillips JA, Knowles JA, et al: Genetic basis of pulmonary arterial hypertension: current understanding and future directions. J Am Coll Cardiol. 2004, 43 (12 Suppl S): 33S-39S.View ArticlePubMedGoogle Scholar
- Newman JH: Treatment of primary pulmonary hypertension – the next generation. N Engl J Med. 2002, 346 (12): 933-935.View ArticlePubMedGoogle Scholar
- Eddahibi S, Morrell N, d'Ortho MP, Naeije R, Adnot S: Pathobiology of pulmonary arterial hypertension. Eur Respir J. 2002, 20 (6): 1559-1572.View ArticlePubMedGoogle Scholar
- Carlson CS, Eberle MA, Kruglyak L, Nickerson DA: Mapping complex disease loci in whole-genome association studies. Nature. 2004, 429 (6990): 446-452.View ArticlePubMedGoogle Scholar
- Geraci MW, Moore M, Gesell T, Yeager ME, Alger L, Golpon H, Gao B, Loyd JE, Tuder RM, Voelkel NF: Gene expression patterns in the lungs of patients with primary pulmonary hypertension: a gene microarray analysis. Circ Res. 2001, 88 (6): 555-562.View ArticlePubMedGoogle Scholar
- Bull TM, Coldren CD, Moore M, Sotto-Santiago SM, Pham DV, Nana-Sinkam SP, Voelkel NF, Geraci MW: Gene microarray analysis of peripheral blood cells in pulmonary arterial hypertension. Am J Respir Crit Care Med. 2004, 170 (8): 911-919.View ArticlePubMedGoogle Scholar
- Kakiuchi C, Ishiwata M, Nanko S, Ozaki N, Iwata N, Umekage T, Tochigi M, Kohda K, Sasaki T, Imamura A, et al: Up-regulation of ADM and SEPX1 in the lymphoblastoid cells of patients in monozygotic twins discordant for schizophrenia. Am J Med Genet B Neuropsychiatr Genet. 2007Google Scholar
- Huang RS, Duan S, Bleibel WK, Kistner EO, Zhang W, Clark TA, Chen TX, Schweitzer AC, Blume JE, Cox NJ, et al: A genome-wide approach to identify genetic variants that contribute to etoposide-induced cytotoxicity. Proc Natl Acad Sci USA. 2007, 104 (23): 9758-9763.View ArticlePubMedPubMed CentralGoogle Scholar
- Baron CA, Liu SY, Hicks C, Gregg JP: Utilization of lymphoblastoid cell lines as a system for the molecular modeling of autism. J Autism Dev Disord. 2006, 36 (8): 973-982.View ArticlePubMedGoogle Scholar
- Dixon AL, Liang L, Moffatt MF, Chen W, Heath S, Wong KC, Taylor J, Burnett E, Gut I, Farrall M, et al: A genome-wide association study of global gene expression. Nat Genet. 2007, 39 (10): 1202-1207.View ArticlePubMedGoogle Scholar
- Lagna G, Ku MM, Nguyen PH, Neuman NA, Davis BN, Hata A: Control of phenotypic plasticity of smooth muscle cells by bone morphogenetic protein signaling through the myocardin-related transcription factors. J Biol Chem. 2007, 282 (51): 37244-37255.View ArticlePubMedPubMed CentralGoogle Scholar
- Monroe JG, Dorshkind K: Fate decisions regulating bone marrow and peripheral B lymphocyte development. Adv Immunol. 2007, 95: 1-50.View ArticlePubMedGoogle Scholar
- Newman JH, Wheeler L, Lane KB, Loyd E, Gaddipati R, Phillips JA, Loyd JE: Mutation in the gene for bone morphogenetic protein receptor II as a cause of primary pulmonary hypertension in a large kindred. N Engl J Med. 2001, 345 (5): 319-324.View ArticlePubMedGoogle Scholar
- Li C, Wong WH: Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. Proc Natl Acad Sci USA. 2001, 98 (1): 31-36.View ArticlePubMedGoogle Scholar
- Gusnanto A, Calza S, Pawitan Y: Identification of differentially expressed genes and false discovery rate in microarray studies. Curr Opin Lipidol. 2007, 18 (2): 187-193.View ArticlePubMedGoogle Scholar
- Zhong S, Li C, Wong WH: ChipInfo: Software for extracting gene annotation and gene ontology information for microarray analysis. Nucleic Acids Res. 2003, 31 (13): 3483-3486.View ArticlePubMedPubMed CentralGoogle Scholar
- The Gene Ontology (GO) project in 2006. Nucleic Acids Res. 2006, D322-326. 34 DatabaseGoogle Scholar
- Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.PubMedGoogle Scholar
- Nettleton D, Recknor J, Reecy JM: Identification of differentially expressed gene categories in microarray studies using nonparametric multivariate analysis. Bioinformatics. 2008, 24 (2): 192-201.View ArticlePubMedGoogle Scholar
- Tavazoie S, Hughes JD, Campbell MJ, Cho RJ, Church GM: Systematic determination of genetic network architecture. Nat Genet. 1999, 22 (3): 281-285.View ArticlePubMedGoogle Scholar
- Kim JH, Sherman ME, Curriero FC, Guengerich FP, Strickland PT, Sutter TR: Expression of cytochromes P450 1A1 and 1B1 in human lung from smokers, non-smokers, and ex-smokers. Toxicol Appl Pharmacol. 2004, 199 (3): 210-219.View ArticlePubMedGoogle Scholar
- Ryu DY, Hodgson E: Constitutive expression and induction of CYP1B1 mRNA in the mouse. J Biochem Mol Toxicol. 1999, 13 (5): 249-251.View ArticlePubMedGoogle Scholar
- Dekker RJ, van Soest S, Fontijn RD, Salamanca S, de Groot PG, VanBavel E, Pannekoek H, Horrevoets AJ: Prolonged fluid shear stress induces a distinct set of endothelial cell genes, most specifically lung Kruppel-like factor (KLF2). Blood. 2002, 100 (5): 1689-1698.View ArticlePubMedGoogle Scholar
- Roos PH, Bolt HM: Cytochrome P450 interactions in human cancers: new aspects considering CYP1B1. Expert Opin Drug Metab Toxicol. 2005, 1 (2): 187-202.View ArticlePubMedGoogle Scholar
- Tsuchiya Y, Nakajima M, Yokoi T: Cytochrome P450-mediated metabolism of estrogens and its regulation in human. Cancer Lett. 2005, 227 (2): 115-124.View ArticlePubMedGoogle Scholar
- Loyd JE: Genetics and pulmonary hypertension. Chest. 2002, 122 (6 Suppl): 284S-286S.View ArticlePubMedGoogle Scholar
- Ulrich S, Taraseviciene-Stewart L, Huber LC, Speich R, Voelkel N: Peripheral blood B lymphocytes derived from patients with idiopathic pulmonary arterial hypertension express a different RNA pattern compared with healthy controls: a cross sectional study. Respir Res. 2008, 9: 20.View ArticlePubMedPubMed CentralGoogle Scholar
- Austin ED, Loyd JE: Genetics and mediators in pulmonary arterial hypertension. Clin Chest Med. 2007, 28 (1): 43-57. vii-viiiView ArticlePubMedPubMed CentralGoogle Scholar
- Geraci MW, Hoshikawa Y, Yeager M, Golpon H, Gesell T, Tuder RM, Voelkel NF: Gene expression profiles in pulmonary hypertension. Chest. 2002, 121 (3 Suppl): 104S-105S.View ArticlePubMedGoogle Scholar
- Tada Y, Laudi S, Harral J, Carr M, Ivester C, Tanabe N, Takiguchi Y, Tatsumi K, Kuriyama T, Nichols WC, et al: Murine pulmonary response to chronic hypoxia is strain specific. Exp Lung Res. 2008, 34 (6): 313-323.View ArticlePubMedPubMed CentralGoogle Scholar
- Meyrick BO, Friedman DB, Billheimer DD, Cogan JD, Prince MA, Phillips JA, Loyd JE: Proteomics of transformed lymphocytes from a family with familial pulmonary arterial hypertension. Am J Respir Crit Care Med. 2008, 177 (1): 99-107.View ArticlePubMedGoogle Scholar
- Newman JH, Phillips JA, Loyd JE: Narrative review: the enigma of pulmonary arterial hypertension: new insights from genetic studies. Ann Intern Med. 2008, 148 (4): 278-283.View ArticlePubMedGoogle Scholar
- Willers ED, Newman JH, Loyd JE, Robbins IM, Wheeler LA, Prince MA, Stanton KC, Cogan JA, Runo JR, Byrne D, et al: Serotonin transporter polymorphisms in familial and idiopathic pulmonary arterial hypertension. Am J Respir Crit Care Med. 2006, 173 (7): 798-802.View ArticlePubMedGoogle Scholar
- Xu Y, See WA: Regulation of plasminogen activator activity in transitional carcinoma cell lines by wound site growth factors. Surg Oncol. 1994, 3 (3): 175-185.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/1/45/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.