https://orcid.org/0000-0003-1522-9038
Figures
Abstract
Introduction
Ritodrine is one of the most commonly used tocolytics in preterm labor, acting as a ß2-adrenergic agonist that reduces intracellular calcium levels and prevents myometrial activation. Ritodrine infusion can result in serious maternal complications, and pulmonary edema is a particular concern among these. The cause of pulmonary edema following ritodrine treatment is multifactorial; however, the contributing genetic factors remain poorly understood. This study investigates the genetic variants associated with ritodrine-induced pulmonary edema.
Methods
In this case-control study, 16 patients who developed pulmonary edema during ritodrine infusion [case], and 16 pregnant women who were treated with ritodrine and did not develop pulmonary edema [control] were included. The control pregnant women were selected after matching for plurality and gestational age at the time of tocolytic use. Maternal blood was collected during admission for tocolytic treatment, and whole exome sequencing was performed with the stored blood samples.
Results
Gene-wise variant burden (GVB) analysis resulted in a total of 71 candidate genes by comparing the cumulative effects of multiple coding variants for 19729 protein-coding genes between the patients with pulmonary edema and the matched controls. Subsequent data analysis selected only the statistically significant and deleterious variants compatible with ritodrine-induced pulmonary edema. Two final candidate variants in CPT2 and ADRA1A were confirmed by Sanger sequencing.
Citation: Lee SM, Park Y, Kim YJ, Hwang H-S, Seo H, Min B-J, et al. (2020) Identifying genetic variants associated with ritodrine-induced pulmonary edema. PLoS ONE 15(11): e0241215. https://doi.org/10.1371/journal.pone.0241215
Editor: Frank T. Spradley, University of Mississippi Medical Center, UNITED STATES
Received: July 24, 2020; Accepted: October 9, 2020; Published: November 9, 2020
Copyright: © 2020 Lee et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The authors confirm that the genetic data supporting the findings of this study are available within the article and its supplementary materials. Clinical data cannot be shared publicly because it contains potentially identifying or sensitive patient information. The Institutional Review Boards and Ethics Committee of Seoul National University Hospital (SNUH), Ewha Woman’s University Hospital (EUMC), and Konkuk University Hospital (KUH) are imposing the restrictions on the data collected in this study. Please contact the following e-mail or phone number for data access. (SNUH Email: irb@snuh.org, Tel: +82-2-2072-0694; EUMC Email: info@ewhactc.co.kr, Tel: +82-2-2650-5114; KUH Email: 20190508@kuh.ac.kr, Tel: +82-2-2030-6542).
Funding: This research was supported by a grant (14172MFDS178) from the Ministry of Food and Drug Safety and the Korean Health Technology R&D Project by the Ministry of Health and Welfare in the Republic of Korea (HI18C2386). The funders played no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist
Abbreviations: GVB, Gene-wise variant burden; cAMP, cyclic adenosine monophosphate; PKA, protein kinase A; WES, Whole exome sequencing; CPT2, carnitine palmitoyltransferase 2; ADRA1A, Adrenoceptor Alpha 1A; 1KGP, 1000 Genomes Project, Phase 3; SIFT, Sorting Intolerant From Tolerant; CADD, Combined Annotation Dependent Depletion
Introduction
Preterm birth is one of the leading causes of perinatal mortality and/or morbidity, and preterm labor constitutes more than one third of the precursors of preterm birth [1]. In women with preterm labor, tocolytics are commonly used to delay preterm birth for antenatal corticosteroids administration and transfer to referral hospitals.
Ritodrine is one of the FDA-approved tocolytics and has been reported to delay preterm birth in patients with preterm labor [2–4]. However, ritodrine can cause several side effects, and ritodrine-induced pulmonary edema is one of the rarest yet most serious side effect that can lead to death [5, 6]. Despite the severity of pulmonary edema, the pathophysiologic mechanism for its development following ritodrine administration remains unknown.
To date, there has been a paucity of information regarding genetic variants associated with ritodrine-induced pulmonary edema. Whole exome sequencing (WES) is being used to find genetic variants associated with various kinds of human diseases. Using WES, we can identify candidate genes, not only for Mendelian disorders, common diseases, and cancer, but also rare variants for complex diseases [7, 8]. Thus, if we can identify patients who are at risk of developing serious complications, such as pulmonary edema, other tocolytics may be considered, such as calcium channel blockers, oxytocin-receptor antagonists, or prostaglandin inhibitors by WES.
In the current study, we examined the genetic variants associated with ritodrine-induced pulmonary edema. Using WES data, we compared gene-wise variant burden (GVB) [9–11] for each of the 19,729 protein-coding gene between the 16 patients suffering from ritodrine-induced pulmonary edema and their 16 matched controls without symptoms after drug use under similar gestational conditions. In silico analysis of deleterious variants was performed, producing two putative East Asian population-specific variants, rs2229291 (CPT2) and rs2229126 (ADRA1A), located in genes with roles compatible with tocolytic-associated pulmonary edema. Next generation sequencing analysis results for two final candidates were confirmed by Sanger sequencing, with 100% concordance. In summary, we identified two novel potential variants that could help identify patients at high-risk of developing ritodrine-induced pulmonary edema.
Materials and methods
Study design
A case-control study was designed to evaluate the genetic variants associated with ritodrine-induced pulmonary edema. Among pregnant women who were admitted for preterm labor between 2000 and 2015 in Seoul National University Hospital, Ewha Woman’s University Hospital, or Konkuk University Hospital, the study population consisted of 16 patients who developed pulmonary edema during ritodrine infusion [case]. As a control group, 16 pregnant women who were treated with ritodrine and did not develop pulmonary edema [control] were included. The control pregnant women were selected based on plurality and gestational age at the time of tocolytics use. This study was approved by the Institutional Review Boards of Seoul National University Hospital, Ewha Woman’s University Hospital, and Konkuk University Hospital (IRB Nos. ECT 06–127–7 and KUH1040034).
Patient and control samples for WES
Maternal blood was collected during admission for preterm labor, and the blood samples were stored at -70°C until analysis. Exome sequencing was performed with stored maternal blood using the Ion AmpliSeq™ Exome panel (Thermo Fisher, USA) to screen coding sequence regions of the entire genome. This panel included the exome of 19,072 genes, and the size of the total targeted region was 57.7 Mb. The panel contained 293,903 primer pairs that were multiplexed into 12 pools to avoid primer-dimer formation and interference during polymerase chain reaction (PCR). The range of amplicons amplified by these oligo primer pairs were from 125 to 275 bp, and the rate of ‘on target’ coverage for this panel was 95.69%. PCR assays were performed directly to amplify 100 ng of genomic DNA samples extracted from peripheral blood to collect target regions using oligo primer pairs of the panel. Reaction parameters were as follows: 99°C for 2 min, followed by 10 cycles of 99°C for 15 s, 60°C for 16 min, and 10°C for 1 min. Following amplification, library construction was performed using the Ion Ampliseq Library Kit Plus (Thermo Scientific, Waltham, MA) as described in the manufacturer's instructions. The libraries were quantified using an Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA) and subsequently diluted to ~10 pM. Thereafter, 33.3 μL of the barcoded libraries were combined in sets of three barcodes. These combined libraries were sequenced using the Ion Proton platform with PI chip V3 (Thermo Scientific), following the manufacturer's instructions.
Data analysis
Reads were mapped to the human genome reference sequence (hg19/GRCh37) using Torrent Suite (v4.4). Aligned reads were sorted and indexed via SAMTools (v0.1.19) [12]. Genome Analysis Toolkit (GATK, v2.8.1) [13–15] was used for local realignment, and base recalibration (dbSNP137, and Mills and 1000 Genome Project gold-standard indels from the hg19 sites). Variant calling was performed using HaplotypeCaller in GATK. ANNOVAR (http://annovar.openbioinformatics.org/) [16] was used to annotate variants. To remove likely false positives, only variants previously identified in 1KGP (1000 Genomes Project, Phase 3) [17] were included for further analysis. Using SnpEff 4.1 (build 2015-01-07) [18], functional variants (missense, stop gained, stop lost, start lost) were selected. A gene-wise variant burden (GVB) score was calculated to predict the cumulative genetic effect for all coding variants of a gene. For each individual, the quantitative value for 19,729 protein-coding genes calculated by the geometric mean of SIFT (Sorting Intolerant From Tolerant) [19] annotated variants in a gene was obtained. Average GVB scores between patients and matched controls were compared by using Student’s t-test. To pinpoint susceptibility loci, Fisher’s exact test was conducted for variants in genes that showed statistically significant associations. We used two in silico bioinformatics tools to select for deleterious variants that fulfilled one or more criteria with SIFT ≤ 0.05 and Combined Annotation Dependent Depletion (CADD) ≥ 15 [20]. The association of biological function, drugs and pathways of genes containing the candidate variants were reviewed via KEGG Drug (https://www.genome.jp/kegg/drug/) [21], KEGG Pathway (https://www.genome.jp/kegg/ pathway.html), and Drug2Gene database [22].
Subsequent validation
For two final candidate loci, we performed independent validation using Sanger sequencing in the 32 individuals. Exon 4 of CPT2, including the rs2229291 position, and Exon 2 of ADRA1A, including the rs2229126 position, were amplified for Sanger sequencing. PCR assays were performed directly to amplify 20 ng genomic DNA samples to collect target regions using oligo-primer pairs. Reaction parameters were as follows: 95°C for 5 min, followed by 35 cycles of 95°C for 30 s, 58°C for 30 s, 72°C for 1 min, and 72°C for 10 min. The RBC HiYield Gel/PCR DNA Mini Kit (Taipai County 220, Taiwan) was used to purify DNA in PCR products. Following purification, PCR samples were directly sequenced using an ABI 3100 semi-automated sequencing analyzer (Applied Biosystems, Lincoln Center Drive Foster City, CA, USA). DNA sequences were analyzed using FinchTV version 1.4.0 (Geospiza, Inc., Seattle, WA, USA).
Statistical analysis
Variant analysis was conducted using three association tests, including Fisher’s exact test under both dominant and recessive models, and a Cochran-Armitage trend test. A comparison of continuous variables between groups was performed using the independent t-test or the Mann-Whitney U-test. A p-value of < .05 was considered statistically significant. SPSS version 21.0 software (IBM Inc., Armonk, NY) was used for the analysis.
Results
Subject population
Table 1 shows the clinical characteristics of the study population. As a result of matching, the frequency of twin pregnancy and gestational age at the time of tocolytic treatment were not different between both groups. In addition, parity and gestational age at delivery were not different between both groups. Cases with pulmonary edema had a higher maternal age and a higher risk of cesarean delivery; however, this difference did not reach statistical significance (p = .077).
Candidate variant prioritization and validation
To prioritize potential variants associated with ritodrine-induced pulmonary edema, we systematically applied several filtering steps (S1 Fig). Of the 229,813 variants identified in 32 individuals, only variants previously reported in 1KGP were included for further analysis. For 19,729 protein-coding genes, we calculated the GVB score, which represents the aggregated effects of genetic variants in a gene, and 71 genes showed statistically significant differences in the average GVB scores between the pulmonary edema group and normal controls. To identify exact variants with a potentially high functional impact on ritodrine-induced pulmonary edema, 43 variants showing an association signal for Fisher’s exact test with a p < 0.05 were selected. Among these, 16 nonsynonymous variants predicted as deleterious by meeting at least one out of two in silico prediction tools were identified (S1 Table). Candidate variants within genes previously implicated in ritodrine-induced adverse drug reactions or included in functionally interesting pathways were prioritized. Finally, the two final candidate variants located in CPT2 and ADRA1A were selected to be responsible for ritodrine-induced pulmonary edema (Table 2). Of these associations, CPT2 rs2229291, known as Perhexiline (calcium channel blocker) target gene on KEGG Drug, was specific for East Asian individuals, in whom the G allele is common (17%) with a substantial effect size (OR = 4.6 (1.3–19.5), Fisher’s exact test p-value = 0.014). ADRA1A rs2229126, which has been implicated to have a relationship with ritodrine via the Drug2Gene database, had a marginal association signal (OR = Inf (1.0-Inf), Fisher’s exact test p-value = 0.053). CPT2 rs2229291 remained statistically significant even after adjustment for the cesarean delivery effect (multiple linear regression p-value = 0.048, S2 Table). Individuals carrying a cumulative number of risk alleles were at increased risk for experiencing ritodrine-induced pulmonary edema compared to the control group with zero or one risk allele (Table 3). The distribution of nonsynonymous variants in ADRA1A and CPT2 genes for 32 individuals is summarized in S3 Table.
Subsequent validation by Sanger sequencing
The two final candidate variants in CPT2 and ADRA1A were confirmed by Sanger sequencing. Sanger sequences for each of the two loci in 32 individuals were perfectly concordant with NGS sequences (S2–S5 Figs). The results showed that individuals with ritodrine-induced pulmonary edema had a significantly higher number of homozygous (CPT2 rs2229291, GG) or heterozygous (CPT2 rs2229291, TG) variants than those without. Variants for ADRA1A (rs2229126, TA) were only identified in individuals with pulmonary edema in the heterozygous state (Fig 1), which implies that life-threatening toxicity could occur following ritodrine use in homozygous ADRA1A variant carriers.
Discussion
In the current study, we demonstrated that two variants in CPT2 and ADRA1A were selected as candidate variants for pulmonary edema. These new potential variants play a role in cAMP/PKA regulation, which supports their putative involvement in the predisposition to ritodrine-induced pulmonary edema in pregnant women.
In clinical practice, ritodrine is one of the most commonly prescribed medications to inhibit preterm labor, especially in Europe and Asia [2, 23]. Ritodrine is a β-adrenergic receptor agonist that acts by reducing intracellular ionized calcium levels and preventing the activation of myometrial contractile proteins. Ritodrine acts at the β2-adrenergic receptor, which activates adenylyl cyclase, thereby increasing cyclic adenosine monophosphate (cAMP) and activating protein kinase A (PKA) [24]. PKA inactivates the myosin light chain kinase, resulting in muscle relaxation.
Pulmonary edema is one of the most serious side effects of ritodrine. The mechanism of pulmonary edema is multifactorial. The conditions related to preterm labor, such as multiple gestation, antenatal corticosteroids, and intravenous infusion of large amounts of crystalloid, can contribute to the development of pulmonary edema. In addition, ritodrine has been shown to decrease total pulmonary resistance and cause retention of sodium and water, resulting in volume overload [25].
In the current study, we identified two final candidate variants in CPT2 and ADRA1A associated with ritodrine-induced pulmonary edema. ADRA1A is a gene encoding the β-adrenergic receptor, which is associated with a calcium channel. The receptor-channel complex contains a G-protein, adenylyl cyclase, and cAMP-dependent kinase, which are related to the mechanism of ritodrine-induced tocolysis. β-adrenergic receptors are also present, not only in the smooth muscle of the uterus but also throughout the lung tissue [26]. Recent reports suggest that β-adrenergic receptors participate in the maintenance of alveolar fluid balance in lung tissue. β2-knockout mice showed decreased ability in the removal of excess water from the airspace and had pulmonary edema following acute lung injury [27]. It is suggested that β-adrenergic receptors participate in the regulation of active Na+ transport for removal of excessive fluid out of the alveolar airspace [28]. In addition, β-adrenergic receptors can enhance endothelial barrier function via inhibition of endothelial contraction and intracellular gap formation, and can regulate surfactant secretion from alveolar cells [29–31]. Thus, decreased function in the β-adrenergic receptor can result in abnormalities in maintenance of alveolar fluid, endothelial barrier function, and surfactant secretion, all of which can contribute to the development of pulmonary edema.
CPT2 encodes carnitine palmitoyltransferase 2 (CPT2), an enzyme required for oxidation of long-chain fatty acids within mitochondria. The heart is a highly active organ that heavily relies on mitochondrial oxidative metabolism, and if disrupted, can lead to cardiac hypertrophic remodeling and altered cardiac function. The CPT2-deficient mouse model developed cardiac hypertrophy and systolic dysfunction, with a significant reduction in blood ejection fraction [32]. In addition, CPT2-deficient mice showed increased vascular permeability in the kidney, spleen, and lung, probably due to impaired junctional properties and decreased barrier function in endothelial cells [33]. Overall, the altered dysfunction in CPT2 may contribute to the pathogenesis of pulmonary edema through decreased function in the right ventricle and increased vascular permeability in lung tissue. There has been a report of decreased activity of CPT2 with rs2229291 in hyperthermic conditions, resulting in impaired fatty acid oxidation and a mitochondrial energy crisis in vascular endothelial cells [34, 35]; however, the direct functional impact of the variant on ritodrine-induced pulmonary edema should be further evaluated in future studies.
To our knowledge, this is the first study to examine the possible variants associated with ritodrine-induced pulmonary edema using WES. Recently, previous studies discussed possible genetic contributions to the pathogenesis of side effects following ritodrine use; however, these studies evaluated patients not only with pulmonary edema (n = 4 in the study of Seo), but also with nonspecific symptoms including dyspnea, tachycardia, and palpitation [10, 36]. From the database of multiple referral hospitals, we identified 16 patients with pulmonary edema, although its occurrence is one of the rarest side effects. However, the current study only evaluated patients in Korea, and further studies in other ethnicities or races are necessary to validate our results. In addition, the small number of this study population is another limitation of the current study. Future directions should focus on appropriate experimental evaluation of how candidate variants affect the pharmacokinetics of ritodrine and pulmonary edema.
The National Institutes of Health defines precision medicine as the tailoring of medical treatment to the individual characteristics of each patients [37]. Precision medicine implies the ability to classify patients into subpopulations according to their susceptibility to a particular disease/prognosis, or their potential response to a specific treatment. With these strategies, preventive or therapeutic intervention can be concentrated on specific populations who will benefit from these interventions and not experience side effects. Although precision medicine has not been universally introduced in obstetrics, the current study shows that we may identify subpopulations who will experience severe side effects following ritodrine administration by WES.
Conclusions
In conclusion, we identified new potential variants in genes that play a role in cAMP/PKA regulation, which supports their putative involvement in the predisposition to ritodrine-induced pulmonary edema in pregnant women.
Supporting information
S1 Fig. Schematic of the data analysis steps.
https://doi.org/10.1371/journal.pone.0241215.s001
(DOCX)
S2 Fig. Results of Sanger sequencing of CPT2 (rs2229291) for the individuals with ritodrine-induced pulmonary edema.
https://doi.org/10.1371/journal.pone.0241215.s002
(DOCX)
S3 Fig. Results of Sanger sequencing of CPT2 (rs2229291) for the individuals without ritodrine-induced pulmonary edema.
https://doi.org/10.1371/journal.pone.0241215.s003
(DOCX)
S4 Fig. Results of Sanger sequencing of ADRA1A (rs2229126) for the individuals with ritodrine-induced pulmonary edema.
https://doi.org/10.1371/journal.pone.0241215.s004
(DOCX)
S5 Fig. Results of Sanger sequencing of ADRA1A (rs2229126) for the individuals without ritodrine-induced pulmonary edema.
https://doi.org/10.1371/journal.pone.0241215.s005
(DOCX)
S1 Table. Frequency distribution for 16 deleterious variants significantly associated with ritodrine induced pulmonary edema.
https://doi.org/10.1371/journal.pone.0241215.s006
(DOCX)
S2 Table. Variants significantly associated with ritodrine-induced pulmonary edema.
https://doi.org/10.1371/journal.pone.0241215.s007
(DOCX)
S3 Table. Distribution of coding variants in ADRA1A and CPT2.
https://doi.org/10.1371/journal.pone.0241215.s008
(DOCX)
References
- 1. Goldenberg RL, Culhane JF, Iams JD, Romero R. Epidemiology and causes of preterm birth. Lancet. 2008;371(9606):75–84. Epub 2008/01/08. pmid:18177778.
- 2. Canadian Preterm Labor Investigators G. Treatment of preterm labor with the beta-adrenergic agonist ritodrine. The New England journal of medicine. 1992;327(5):308–12. Epub 1992/07/30. pmid:1620169.
- 3. Leveno KJ, Klein VR, Guzick DS, Young DC, Hankins GD, Williams ML. Single-centre randomised trial of ritodrine hydrochloride for preterm labour. Lancet. 1986;1(8493):1293–6. Epub 1986/06/07. pmid:2872429.
- 4. Li X, Zhang Y, Shi Z. Ritodrine in the treatment of preterm labour: a meta-analysis. The Indian journal of medical research. 2005;121(2):120–7. Epub 2005/03/10. pmid:15756046.
- 5. Shinohara S, Sunami R, Uchida Y, Hirata S, Suzuki K. Association between total dose of ritodrine hydrochloride and pulmonary oedema in twin pregnancy: a retrospective cohort study in Japan. BMJ open. 2017;7(12):e018118. Epub 2018/01/01. pmid:29289935; PubMed Central PMCID: PMC5778295.
- 6. Lamont RF. The pathophysiology of pulmonary oedema with the use of beta-agonists. BJOG: an international journal of obstetrics and gynaecology. 2000;107(4):439–44. Epub 2000/04/12. pmid:10759259.
- 7. Kaname T, Yanagi K, Naritomi K. A commentary on the promise of whole-exome sequencing in medical genetics. Journal of human genetics. 2014;59(3):117–8. Epub 2014/02/07. pmid:24500682.
- 8. Lee S, Abecasis GR, Boehnke M, Lin X. Rare-variant association analysis: study designs and statistical tests. American journal of human genetics. 2014;95(1):5–23. Epub 2014/07/06. pmid:24995866; PubMed Central PMCID: PMC4085641.
- 9. Gaedigk A, Sangkuhl K, Whirl-Carrillo M, Klein T, Leeder JS. Prediction of CYP2D6 phenotype from genotype across world populations. Genetics in medicine: official journal of the American College of Medical Genetics. 2016. Epub 2016/07/09. pmid:27388693.
- 10. Seo H, Kwon EJ, You YA, Park Y, Min BJ, Yoo K, et al. Deleterious genetic variants in ciliopathy genes increase risk of ritodrine-induced cardiac and pulmonary side effects. BMC medical genomics. 2018;11(1):4. Epub 2018/01/26. pmid:29368655; PubMed Central PMCID: PMC5784607.
- 11. Park Y, Kim H, Choi JY, Yun S, Min BJ, Seo ME, et al. Star Allele-Based Haplotyping versus Gene-Wise Variant Burden Scoring for Predicting 6-Mercaptopurine Intolerance in Pediatric Acute Lymphoblastic Leukemia Patients. Frontiers in pharmacology. 2019;10:654. Epub 2019/06/28. pmid:31244663; PubMed Central PMCID: PMC6580331.
- 12. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9. Epub 2009/06/10. pmid:19505943; PubMed Central PMCID: PMC2723002.
- 13. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome research. 2010;20(9):1297–303. Epub 2010/07/21. pmid:20644199; PubMed Central PMCID: PMC2928508.
- 14. DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nature genetics. 2011;43(5):491–8. Epub 2011/04/12. pmid:21478889; PubMed Central PMCID: PMC3083463.
- 15. Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, Del Angel G, Levy-Moonshine A, et al. From FastQ data to high confidence variant calls: the Genome Analysis Toolkit best practices pipeline. Current protocols in bioinformatics. 2013;43:11 0 1–33. Epub 2014/11/29. pmid:25431634; PubMed Central PMCID: PMC4243306.
- 16. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic acids research. 2010;38(16):e164. Epub 2010/07/06. pmid:20601685; PubMed Central PMCID: PMC2938201.
- 17. Genomes Project C, Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, et al. A global reference for human genetic variation. Nature. 2015;526(7571):68–74. Epub 2015/10/04. pmid:26432245; PubMed Central PMCID: PMC4750478.
- 18. Cingolani P, Platts A, Wang le L, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly. 2012;6(2):80–92. Epub 2012/06/26. pmid:22728672; PubMed Central PMCID: PMC3679285.
- 19. Ng PC, Henikoff S. SIFT: Predicting amino acid changes that affect protein function. Nucleic acids research. 2003;31(13):3812–4. Epub 2003/06/26. pmid:12824425; PubMed Central PMCID: PMC168916.
- 20. Kircher M, Witten DM, Jain P, O'Roak BJ, Cooper GM, Shendure J. A general framework for estimating the relative pathogenicity of human genetic variants. Nature genetics. 2014;46(3):310–5. Epub 2014/02/04. pmid:24487276; PubMed Central PMCID: PMC3992975.
- 21. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic acids research. 2000;28(1):27–30. Epub 1999/12/11. pmid:10592173; PubMed Central PMCID: PMC102409.
- 22. Roider HG, Pavlova N, Kirov I, Slavov S, Slavov T, Uzunov Z, et al. Drug2Gene: an exhaustive resource to explore effectively the drug-target relation network. BMC bioinformatics. 2014;15:68. Epub 2014/03/13. pmid:24618344; PubMed Central PMCID: PMC4234465.
- 23. Hwang HS, Na SH, Hur SE, Lee SA, Lee KA, Cho GJ, et al. Practice patterns in the management of threatened preterm labor in Korea: A multicenter retrospective study. Obstetrics & gynecology science. 2015;58(3):203–9. Epub 2015/05/30. pmid:26023669; PubMed Central PMCID: PMC4444516.
- 24. Arrowsmith S, Kendrick A, Wray S. Drugs acting on the pregnant uterus. Obstetrics, gynaecology and reproductive medicine. 2010;20(8):241–7. Epub 2010/08/01. pmid:24443652; PubMed Central PMCID: PMC3892726.
- 25. Hankins GD, Hauth JC, Cissik JH, Kuehl TJ. Effects of ritodrine hydrochloride on arteriovenous blood gas and shunt in healthy pregnant yellow baboons. American journal of obstetrics and gynecology. 1988;158(3 Pt 1):658–63. Epub 1988/03/01. pmid:3126659.
- 26. Mutlu GM, Factor P. Alveolar epithelial beta2-adrenergic receptors. American journal of respiratory cell and molecular biology. 2008;38(2):127–34. Epub 2007/08/22. pmid:17709598; PubMed Central PMCID: PMC2214676.
- 27. Mutlu GM, Dumasius V, Burhop J, McShane PJ, Meng FJ, Welch L, et al. Upregulation of alveolar epithelial active Na+ transport is dependent on beta2-adrenergic receptor signaling. Circulation research. 2004;94(8):1091–100. Epub 2004/03/16. pmid:15016730.
- 28. Mutlu GM, Koch WJ, Factor P. Alveolar epithelial beta 2-adrenergic receptors: their role in regulation of alveolar active sodium transport. American journal of respiratory and critical care medicine. 2004;170(12):1270–5. Epub 2004/09/28. pmid:15447942.
- 29. Liu F, Verin AD, Borbiev T, Garcia JG. Role of cAMP-dependent protein kinase A activity in endothelial cell cytoskeleton rearrangement. American journal of physiology Lung cellular and molecular physiology. 2001;280(6):L1309–17. Epub 2001/05/15. pmid:11350812.
- 30. Patterson CE, Lum H, Schaphorst KL, Verin AD, Garcia JG. Regulation of endothelial barrier function by the cAMP-dependent protein kinase. Endothelium: journal of endothelial cell research. 2000;7(4):287–308. Epub 2001/02/24. pmid:11201526.
- 31. Chen Q, Bates SR, Fisher AB. Secretagogues increase the expression of surfactant protein A receptors on lung type II cells. The Journal of biological chemistry. 1996;271(41):25277–83. Epub 1996/10/11. pmid:8810290.
- 32. Pereyra AS, Hasek LY, Harris KL, Berman AG, Damen FW, Goergen CJ, et al. Loss of cardiac carnitine palmitoyltransferase 2 results in rapamycin-resistant, acetylation-independent hypertrophy. The Journal of biological chemistry. 2017;292(45):18443–56. Epub 2017/09/17. pmid:28916721; PubMed Central PMCID: PMC5682957.
- 33. Xiong J, Kawagishi H, Yan Y, Liu J, Wells QS, Edmunds LR, et al. A Metabolic Basis for Endothelial-to-Mesenchymal Transition. Molecular cell. 2018;69(4):689–98 e7. Epub 2018/02/13. pmid:29429925; PubMed Central PMCID: PMC5816688.
- 34. Chen Y, Mizuguchi H, Yao D, Ide M, Kuroda Y, Shigematsu Y, et al. Thermolabile phenotype of carnitine palmitoyltransferase II variations as a predisposing factor for influenza-associated encephalopathy. FEBS Lett. 2005;579(10):2040–4. Epub 2005/04/07. pmid:15811315.
- 35. Yao D, Mizuguchi H, Yamaguchi M, Yamada H, Chida J, Shikata K, et al. Thermal instability of compound variants of carnitine palmitoyltransferase II and impaired mitochondrial fuel utilization in influenza-associated encephalopathy. Human mutation. 2008;29(5):718–27. Epub 2008/02/29. pmid:18306170.
- 36. Baek MY, Hwang HS, Park JY, Chung JE, Lee KE, Lee GY, et al. Association between CACNA1C gene polymorphisms and ritodrine-induced adverse events in preterm labor patients. European journal of clinical pharmacology. 2017;73(7):837–42. Epub 2017/04/10. pmid:28391406.
- 37. Gluck C, Ko YA, Susztak K. Precision Medicine Approaches to Diabetic Kidney Disease: Tissue as an Issue. Current diabetes reports. 2017;17(5):30. Epub 2017/04/02. pmid:28364354; PubMed Central PMCID: PMC5772762.