Development of the ehive Digital Health App: Protocol for a Centralized Research Platform

BACKGROUND

The increasing use of smartphones, wearables, and connected devices has enabled the increasing application of digital technologies for research. Remote digital study platforms comprise a patient-interfacing digital application that enables multimodal data collection from a mobile app and connected sources. They offer an opportunity to recruit at scale, acquire data longitudinally at a high frequency, and engage study participants at any time of the day in any place. Few published descriptions of centralized digital research platforms provide a framework for their development.

OBJECTIVE

This study aims to serve as a road map for those seeking to develop a centralized digital research platform. We describe the technical and functional aspects of the ehive app, the centralized digital research platform of the Hasso Plattner Institute for Digital Health at Mount Sinai Hospital, New York, New York. We then provide information about ongoing studies hosted on ehive, including usership statistics and data infrastructure. Finally, we discuss our experience with ehive in the broader context of the current landscape of digital health research platforms.

METHODS

The ehive app is a multifaceted and patient-facing central digital research platform that permits the collection of e-consent for digital health studies. An overview of its development, its e-consent process, and the tools it uses for participant recruitment and retention are provided. Data integration with the platform and the infrastructure supporting its operations are discussed; furthermore, a description of its participant- and researcher-facing dashboard interfaces and the e-consent architecture is provided.

RESULTS

The ehive platform was launched in 2020 and has successfully hosted 8 studies, namely 6 observational studies and 2 clinical trials. Approximately 1484 participants downloaded the app across 36 states in the United States. The use of recruitment methods such as bulk messaging through the EPIC electronic health records and standard email portals enables broad recruitment. Light-touch engagement methods, used in an automated fashion through the platform, maintain high degrees of engagement and retention. The ehive platform demonstrates the successful deployment of a central digital research platform that can be modified across study designs.

CONCLUSIONS

Centralized digital research platforms such as ehive provide a novel tool that allows investigators to expand their research beyond their institution, engage in large-scale longitudinal studies, and combine multimodal data streams. The ehive platform serves as a model for groups seeking to develop similar digital health research programs.

INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID)

DERR1-10.2196/49204.

Economic Evaluation of Digital Therapeutic Care Apps for Unsupervised Treatment of Low Back Pain: Monte Carlo Simulation

BACKGROUND

Digital therapeutic care (DTC) programs are unsupervised app-based treatments that provide video exercises and educational material to patients with nonspecific low back pain during episodes of pain and functional disability. German statutory health insurance can reimburse DTC programs since 2019, but evidence on efficacy and reasonable pricing remains scarce. This paper presents a probabilistic sensitivity analysis (PSA) to evaluate the efficacy and cost-utility of a DTC app against treatment as usual (TAU) in Germany.

OBJECTIVE

The aim of this study was to perform a PSA in the form of a Monte Carlo simulation based on the deterministic base case analysis to account for model assumptions and parameter uncertainty. We also intend to explore to what extent the results in this probabilistic analysis differ from the results in the base case analysis and to what extent a shortage of outcome data concerning quality-of-life (QoL) metrics impacts the overall results.

METHODS

The PSA builds upon a state-transition Markov chain with a 4-week cycle length over a model time horizon of 3 years from a recently published deterministic cost-utility analysis. A Monte Carlo simulation with 10,000 iterations and a cohort size of 10,000 was employed to evaluate the cost-utility from a societal perspective. Quality-adjusted life years (QALYs) were derived from Veterans RAND 6-Dimension (VR-6D) and Short-Form 6-Dimension (SF-6D) single utility scores. Finally, we also simulated reducing the price for a 3-month app prescription to analyze at which price threshold DTC would result in being the dominant strategy over TAU in Germany.

RESULTS

The Monte Carlo simulation yielded on average a €135.97 (a currency exchange rate of EUR €1=US $1.069 is applicable) incremental cost and 0.004 incremental QALYs per person and year for the unsupervised DTC app strategy compared to in-person physiotherapy in Germany. The corresponding incremental cost-utility ratio (ICUR) amounts to an additional €34,315.19 per additional QALY. DTC yielded more QALYs in 54.96% of the iterations. DTC dominates TAU in 24.04% of the iterations for QALYs. Reducing the app price in the simulation from currently €239.96 to €164.61 for a 3-month prescription could yield a negative ICUR and thus make DTC the dominant strategy, even though the estimated probability of DTC being more effective than TAU is only 54.96%.

CONCLUSIONS

Decision-makers should be cautious when considering the reimbursement of DTC apps since no significant treatment effect was found, and the probability of cost-effectiveness remains below 60% even for an infinite willingness-to-pay threshold. More app-based studies involving the utilization of QoL outcome parameters are urgently needed to account for the low and limited precision of the available QoL input parameters, which are crucial to making profound recommendations concerning the cost-utility of novel apps.

Multi-ancestry genome- and phenome-wide association studies of diverticular disease in electronic health records with natural language processing enriched phenotyping algorithm

OBJECTIVE

Diverticular disease (DD) is one of the most prevalent conditions encountered by gastroenterologists, affecting ~50% of Americans before the age of 60. Our aim was to identify genetic risk variants and clinical phenotypes associated with DD, leveraging multiple electronic health record (EHR) data sources of 91,166 multi-ancestry participants with a Natural Language Processing (NLP) technique.

MATERIALS AND METHODS

We developed a NLP-enriched phenotyping algorithm that incorporated colonoscopy or abdominal imaging reports to identify patients with diverticulosis and diverticulitis from multicenter EHRs. We performed genome-wide association studies (GWAS) of DD in European, African and multi-ancestry participants, followed by phenome-wide association studies (PheWAS) of the risk variants to identify their potential comorbid/pleiotropic effects in clinical phenotypes.

RESULTS

Our developed algorithm showed a significant improvement in patient classification performance for DD analysis (algorithm PPVs ≥ 0.94), with up to a 3.5 fold increase in terms of the number of identified patients than the traditional method. Ancestry-stratified analyses of diverticulosis and diverticulitis of the identified subjects replicated the well-established associations between ARHGAP15 loci with DD, showing overall intensified GWAS signals in diverticulitis patients compared to diverticulosis patients. Our PheWAS analyses identified significant associations between the DD GWAS variants and circulatory system, genitourinary, and neoplastic EHR phenotypes.

DISCUSSION

As the first multi-ancestry GWAS-PheWAS study, we showcased that heterogenous EHR data can be mapped through an integrative analytical pipeline and reveal significant genotype-phenotype associations with clinical interpretation.

CONCLUSION

A systematic framework to process unstructured EHR data with NLP could advance a deep and scalable phenotyping for better patient identification and facilitate etiological investigation of a disease with multilayered data.

Digital Therapeutic Care Apps With Decision-Support Interventions for People With Low Back Pain in Germany: Cost-Effectiveness Analysis

BACKGROUND

Digital therapeutic care apps provide a new effective and scalable approach for people with nonspecific low back pain (LBP). Digital therapeutic care apps are also driven by personalized decision-support interventions that support the user in self-managing LBP, and may induce prolonged behavior change to reduce the frequency and intensity of pain episodes. However, these therapeutic apps are associated with high attrition rates, and the initial prescription cost is higher than that of face-to-face physiotherapy. In Germany, digital therapeutic care apps are now being reimbursed by statutory health insurance; however, price targets and cost-driving factors for the formation of the reimbursement rate remain unexplored.

OBJECTIVE

The aim of this study was to evaluate the cost-effectiveness of a digital therapeutic care app compared to treatment as usual (TAU) in Germany. We further aimed to explore under which circumstances the reimbursement rate could be modified to consider value-based pricing.

METHODS

We developed a state-transition Markov model based on a best-practice analysis of prior LBP-related decision-analytic models, and evaluated the cost utility of a digital therapeutic care app compared to TAU in Germany. Based on a 3-year time horizon, we simulated the incremental cost and quality-adjusted life years (QALYs) for people with nonacute LBP from the societal perspective. In the deterministic sensitivity and scenario analyses, we focused on diverging attrition rates and app cost to assess our model's robustness and conditions for changing the reimbursement rate. All costs are reported in Euro (€1=US $1.12).

RESULTS

Our base case results indicated that the digital therapeutic care strategy led to an incremental cost of €121.59, but also generated 0.0221 additional QALYs compared to the TAU strategy, with an estimated incremental cost-effectiveness ratio (ICER) of €5486 per QALY. The sensitivity analysis revealed that the reimbursement rate and the capability of digital therapeutic care to prevent reoccurring LBP episodes have a significant impact on the ICER. At the same time, the other parameters remained unaffected and thus supported the robustness of our model. In the scenario analysis, the different model time horizons and attrition rates strongly influenced the economic outcome. Reducing the cost of the app to €99 per 3 months or decreasing the app's attrition rate resulted in digital therapeutic care being significantly less costly with more generated QALYs, and is thus considered to be the dominant strategy over TAU.

CONCLUSIONS

The current reimbursement rate for a digital therapeutic care app in the statutory health insurance can be considered a cost-effective measure compared to TAU. The app's attrition rate and effect on the patient's prolonged behavior change essentially influence the settlement of an appropriate reimbursement rate. Future value-based pricing targets should focus on additional outcome parameters besides pain intensity and functional disability by including attrition rates and the app's long-term effect on quality of life.

Rare coding variants in 35 genes associate with circulating lipid levels-A multi-ancestry analysis of 170,000 exomes

Large-scale gene sequencing studies for complex traits have the potential to identify causal genes with therapeutic implications. We performed gene-based association testing of blood lipid levels with rare (minor allele frequency < 1%) predicted damaging coding variation by using sequence data from >170,000 individuals from multiple ancestries: 97,493 European, 30,025 South Asian, 16,507 African, 16,440 Hispanic/Latino, 10,420 East Asian, and 1,182 Samoan. We identified 35 genes associated with circulating lipid levels; some of these genes have not been previously associated with lipid levels when using rare coding variation from population-based samples. We prioritize 32 genes in array-based genome-wide association study (GWAS) loci based on aggregations of rare coding variants; three (EVI5, SH2B3, and PLIN1) had no prior association of rare coding variants with lipid levels. Most of our associated genes showed evidence of association among multiple ancestries. Finally, we observed an enrichment of gene-based associations for low-density lipoprotein cholesterol drug target genes and for genes closest to GWAS index single-nucleotide polymorphisms (SNPs). Our results demonstrate that gene-based associations can be beneficial for drug target development and provide evidence that the gene closest to the array-based GWAS index SNP is often the functional gene for blood lipid levels.

Digital Therapeutic Care and Decision Support Interventions for People With Low Back Pain: Systematic Review

Background

Low back pain (LBP) is the leading cause of worldwide years lost because of disability, with a tremendous economic burden for health care systems. Digital therapeutic care (DTC) programs provide a scalable, universally accessible, and low-cost approach to the multidisciplinary treatment of LBP. Moreover, novel decision support interventions such as personalized feedback messages, push notifications, and data-driven activity recommendations amplify DTC by guiding the user through the program while aiming to increase overall engagement and sustainable behavior change.

Objective

This systematic review aims to synthesize recent scientific literature on the impact of DTC apps for people with LBP and outline the implementation of add-on decision support interventions, including their effect on user retention and attrition rates.

Methods

We searched bibliographic databases, including MEDLINE, Cochrane Library, Web of Science, and the Physiotherapy Evidence Database, from March 1, 2016, to October 15, 2020, in accordance with the PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) guidelines and conducted this review based on related previously published systematic reviews. Besides randomized controlled trials (RCTs), we also included study designs with the evidence level of at least a retrospective comparative study. This enables the consideration of real-world user-generated data and provides information regarding the adoption and effectiveness of DTC apps in a real-life setting. For the appraisal of the risk of bias, we used the Risk of Bias 2 Tool and the Risk of Bias in Non-Randomized Studies of Interventions Tool for the RCTs and nonrandomized trials, respectively. The included studies were narratively synthesized regarding primary and secondary outcome measures, DTC components, applied decision support interventions, user retention, and attrition rates.

Results

We retrieved 1388 citations, of which 12 studies are included in this review. Of the 12 studies, 6 (50%) were RCTs and 6 (50%) were nonrandomized trials. In all included studies, lower pain levels and increased functionality compared with baseline values were observed in the DTC intervention group. A between-group comparison revealed significant improvements in pain and functionality levels in 67% (4/6) of the RCTs. The study population was mostly homogeneous, with predominantly female, young to middle-aged participants of normal to moderate weight. The methodological quality assessment revealed moderate to high risks of biases, especially in the nonrandomized trials.

Conclusions

This systematic review demonstrates the benefits of DTC for people with LBP. There is also evidence that decision support interventions benefit overall engagement with the app and increase participants’ ability to self-manage their recovery process. Finally, including retrospective evaluation studies of real-world user-generated data in future systematic reviews of digital health intervention trials can reveal new insights into the benefits, challenges, and real-life adoption of DTC programs.

Determinants of penetrance and variable expressivity in monogenic metabolic conditions across 77,184 exomes

Hundreds of thousands of genetic variants have been reported to cause severe monogenic diseases, but the probability that a variant carrier develops the disease (termed penetrance) is unknown for virtually all of them. Additionally, the clinical utility of common polygenetic variation remains uncertain. Using exome sequencing from 77,184 adult individuals (38,618 multi-ancestral individuals from a type 2 diabetes case-control study and 38,566 participants from the UK Biobank, for whom genotype array data were also available), we apply clinical standard-of-care gene variant curation for eight monogenic metabolic conditions. Rare variants causing monogenic diabetes and dyslipidemias display effect sizes significantly larger than the top 1% of the corresponding polygenic scores. Nevertheless, penetrance estimates for monogenic variant carriers average 60% or lower for most conditions. We assess epidemiologic and genetic factors contributing to risk prediction in monogenic variant carriers, demonstrating that inclusion of polygenic variation significantly improves biomarker estimation for two monogenic dyslipidemias.

A Resilience-Building App to Support the Mental Health of Health Care Workers in the COVID-19 Era: Design Process, Distribution, and Evaluation

Background

The COVID-19 pandemic has resulted in increased strain on health care systems and negative psychological effects on health care workers (HCWs). This is anticipated to result in long-term negative mental health effects on the population, with HCWs representing a particularly vulnerable group. The scope of the COVID-19 pandemic necessitates the development of a scalable mental health platform to provide services to large numbers of at-risk or affected individuals. The Mount Sinai Health System in New York City was at the epicenter of the pandemic in the United States.

Objective

The Center for Stress, Resilience, and Personal Growth (CSRPG) was created to address the current and anticipated psychological impact of the pandemic on the HCWs in the health system. The mission of the Center is to support the resilience and mental health of employees through educational offerings, outreach, and clinical care. Our aim was to build a mobile app to support the newly founded Center in its mission.

Methods

We built the app as a standalone digital platform that hosts a suite of tools that users can interact with on a daily basis. With consideration for the Center’s aims, we determined the overall vision, initiatives, and goals for the Wellness Hub app, followed by specific milestone tasks and deliverables for development. We defined the app’s primary features based on the mental health assessment and needs of HCWs. Feature definition was informed by the results of a resilience survey widely distributed to Mount Sinai HCWs and by the resources offered at CSRPG, including workshop content.

Results

We launched our app over the course of two phases, the first phase being a “soft” launch and the second being a broader launch to all of Mount Sinai. Of the 231 HCWs who downloaded the app, 173 (74.9%) completed our baseline assessment of all mental health screeners in the app. Results from the baseline assessment show that more than half of the users demonstrate a need for support in at least one psychological area. As of 3 months after the Phase 2 launch, approximately 55% of users re-entered the app after their first opening to explore additional features, with an average of 4 app openings per person.

Conclusions

To address the mental health needs of HCWs during the COVID-19 pandemic, the Wellness Hub app was built and deployed throughout the Mount Sinai Health System. To our knowledge, this is the first resilience app of its kind. The Wellness Hub app is a promising proof of concept, with room to grow, for those who wish to build a secure mobile health app to support their employees, communities, or others in managing and improving mental and physical well-being. It is a novel tool offering mental health support broadly.

Coronavirus 2019 and People Living With Human Immunodeficiency Virus: Outcomes for Hospitalized Patients in New York City

BACKGROUND

There are limited data regarding the clinical impact of coronavirus disease 2019 (COVID-19) on people living with human immunodeficiency virus (PLWH). In this study, we compared outcomes for PLWH with COVID-19 to a matched comparison group.

METHODS

We identified 88 PLWH hospitalized with laboratory-confirmed COVID-19 in our hospital system in New York City between 12 March and 23 April 2020. We collected data on baseline clinical characteristics, laboratory values, HIV status, treatment, and outcomes from this group and matched comparators (1 PLWH to up to 5 patients by age, sex, race/ethnicity, and calendar week of infection). We compared clinical characteristics and outcomes (death, mechanical ventilation, hospital discharge) for these groups, as well as cumulative incidence of death by HIV status.

RESULTS

Patients did not differ significantly by HIV status by age, sex, or race/ethnicity due to the matching algorithm. PLWH hospitalized with COVID-19 had high proportions of HIV virologic control on antiretroviral therapy. PLWH had greater proportions of smoking (P < .001) and comorbid illness than uninfected comparators. There was no difference in COVID-19 severity on admission by HIV status (P = .15). Poor outcomes for hospitalized PLWH were frequent but similar to proportions in comparators; 18% required mechanical ventilation and 21% died during follow-up (compared with 23% and 20%, respectively). There was similar cumulative incidence of death over time by HIV status (P = .94).

CONCLUSIONS

We found no differences in adverse outcomes associated with HIV infection for hospitalized COVID-19 patients compared with a demographically similar patient group.

Exome sequencing of 20,791 cases of type 2 diabetes and 24,440 controls

Protein-coding genetic variants that strongly affect disease risk can yield relevant clues to disease pathogenesis. Here we report exome-sequencing analyses of 20,791 individuals with type 2 diabetes (T2D) and 24,440 non-diabetic control participants from 5 ancestries. We identify gene-level associations of rare variants (with minor allele frequencies of less than 0.5%) in 4 genes at exome-wide significance, including a series of more than 30 SLC30A8 alleles that conveys protection against T2D, and in 12 gene sets, including those corresponding to T2D drug targets (P = 6.1 × 10-3) and candidate genes from knockout mice (P = 5.2 × 10-3). Within our study, the strongest T2D gene-level signals for rare variants explain at most 25% of the heritability of the strongest common single-variant signals, and the gene-level effect sizes of the rare variants that we observed in established T2D drug targets will require 75,000-185,000 sequenced cases to achieve exome-wide significance. We propose a method to interpret these modest rare-variant associations and to incorporate these associations into future target or gene prioritization efforts.

Serum 6-Bromotryptophan Levels Identified as a Risk Factor for CKD Progression

Background Metabolite levels reflect physiologic homeostasis and may serve as biomarkers of disease progression. Identifying metabolites associated with APOL1 risk alleles-genetic variants associated with CKD risk commonly present in persons of African descent-may reveal novel markers of CKD progression relevant to other populations.Methods We evaluated associations between the number of APOL1 risk alleles and 760 serum metabolites identified via untargeted profiling in participants of the African American Study of Kidney Disease and Hypertension (AASK) (n=588; Bonferroni significance threshold P<6.5×10^-5) and replicated findings in 678 black participants with CKD in BioMe, an electronic medical record-linked biobank. We tested the metabolite association with CKD progression in AASK, BioMe, and the Modification of Diet in Renal Disease (MDRD) Study.Results One metabolite, 6-bromotryptophan, was significant in AASK (P=4.7×10^-5) and replicated in BioMe (P=5.7×10^-3) participants, with lower levels associated with more APOL1 risk alleles. Lower levels of 6-bromotryptophan were associated with CKD progression in AASK and BioMe participants and in white participants in the MDRD Study, independent of demographics and clinical characteristics, including baseline GFR (adjusted hazard ratio per two-fold higher 6-bromotryptophan level, AASK, 0.76; 95% confidence interval [95% CI], 0.64 to 0.91; BioMe, 0.61; 95% CI, 0.43 to 0.85; MDRD, 0.52; 95% CI, 0.34 to 0.79). The interaction between the APOL1 risk alleles and 6-bromotryptophan was not significant. The identity of 6-bromotryptophan was confirmed in experiments comparing its molecular signature with that of authentic standards of other bromotryptophan isomers.Conclusions Serum 6-bromotryptophan is a consistent and novel risk factor for CKD progression.

Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.

Genome-wide study of resistant hypertension identified from electronic health records

Resistant hypertension is defined as high blood pressure that remains above treatment goals in spite of the concurrent use of three antihypertensive agents from different classes. Despite the important health consequences of resistant hypertension, few studies of resistant hypertension have been conducted. To perform a genome-wide association study for resistant hypertension, we defined and identified cases of resistant hypertension and hypertensives with treated, controlled hypertension among >47,500 adults residing in the US linked to electronic health records (EHRs) and genotyped as part of the electronic MEdical Records & GEnomics (eMERGE) Network. Electronic selection logic using billing codes, laboratory values, text queries, and medication records was used to identify resistant hypertension cases and controls at each site, and a total of 3,006 cases of resistant hypertension and 876 controlled hypertensives were identified among eMERGE Phase I and II sites. After imputation and quality control, a total of 2,530,150 SNPs were tested for an association among 2,830 multi-ethnic cases of resistant hypertension and 876 controlled hypertensives. No test of association was genome-wide significant in the full dataset or in the dataset limited to European American cases (n = 1,719) and controls (n = 708). The most significant finding was CLNK rs13144136 at p = 1.00x10^-6 (odds ratio = 0.68; 95% CI = 0.58–0.80) in the full dataset with similar results in the European American only dataset. We also examined whether SNPs known to influence blood pressure or hypertension also influenced resistant hypertension. None was significant after correction for multiple testing. These data highlight both the difficulties and the potential utility of EHR-linked genomic data to study clinically-relevant traits such as resistant hypertension.

Genome-Wide Association of CKD Progression: The Chronic Renal Insufficiency Cohort Study

The rate of decline of renal function varies significantly among individuals with CKD. To understand better the contribution of genetics to CKD progression, we performed a genome-wide association study among participants in the Chronic Renal Insufficiency Cohort Study. Our outcome of interest was CKD progression measured as change in eGFR over time among 1331 blacks and 1476 whites with CKD. We stratified all analyses by race and subsequently, diabetes status. Single-nucleotide polymorphisms (SNPs) that surpassed a significance threshold of P<1×10^-6 for association with eGFR slope were selected as candidates for follow-up and secondarily tested for association with proteinuria and time to ESRD. We identified 12 such SNPs among black patients and six such SNPs among white patients. We were able to conduct follow-up analyses of three candidate SNPs in similar (replication) cohorts and eight candidate SNPs in phenotype-related (validation) cohorts. Among blacks without diabetes, rs653747 in LINC00923 replicated in the African American Study of Kidney Disease and Hypertension cohort (discovery P=5.42×10^-7; replication P=0.039; combined P=7.42×10^-9). This SNP also associated with ESRD (hazard ratio, 2.0 (95% confidence interval, 1.5 to 2.7); P=4.90×10^-6). Similarly, rs931891 in LINC00923 associated with eGFR decline (P=1.44×10^-4) in white patients without diabetes. In summary, SNPs in LINC00923, an RNA gene expressed in the kidney, significantly associated with CKD progression in individuals with nondiabetic CKD. However, the lack of equivalent cohorts hampered replication for most discovery loci. Further replication of our findings in comparable study populations is warranted.

Genetic variation among 82 pharmacogenes: The PGRNseq data from the eMERGE network

Genetic variation can affect drug response in multiple ways, although it remains unclear how rare genetic variants affect drug response. The electronic Medical Records and Genomics (eMERGE) Network, collaborating with the Pharmacogenomics Research Network, began eMERGE‐PGx, a targeted sequencing study to assess genetic variation in 82 pharmacogenes critical for implementation of “precision medicine.” The February 2015 eMERGE‐PGx data release includes sequence‐derived data from ∼5,000 clinical subjects. We present the variant frequency spectrum categorized by variant type, ancestry, and predicted function. We found 95.12% of genes have variants with a scaled Combined Annotation‐Dependent Depletion score above 20, and 96.19% of all samples had one or more Clinical Pharmacogenetics Implementation Consortium Level A actionable variants. These data highlight the distribution and scope of genetic variation in relevant pharmacogenes, identifying challenges associated with implementing clinical sequencing for drug treatment at a broader level, underscoring the importance for multifaceted research in the execution of precision medicine.

The phenotypic legacy of admixture between modern humans and Neandertals

Many modern human genomes retain DNA inherited from interbreeding with archaic hominins, such as Neandertals, yet the influence of this admixture on human traits is largely unknown. We analyzed the contribution of common Neandertal variants to over 1000 electronic health record (EHR)-derived phenotypes in ~28,000 adults of European ancestry. We discovered and replicated associations of Neandertal alleles with neurological, psychiatric, immunological, and dermatological phenotypes. Neandertal alleles together explained a significant fraction of the variation in risk for depression and skin lesions resulting from sun exposure (actinic keratosis), and individual Neandertal alleles were significantly associated with specific human phenotypes, including hypercoagulation and tobacco use. Our results establish that archaic admixture influences disease risk in modern humans, provide hypotheses about the effects of hundreds of Neandertal haplotypes, and demonstrate the utility of EHR data in evolutionary analyses.

Association of Arrhythmia-Related Genetic Variants With Phenotypes Documented in Electronic Medical Records

IMPORTANCE

Large-scale DNA sequencing identifies incidental rare variants in established Mendelian disease genes, but the frequency of related clinical phenotypes in unselected patient populations is not well established. Phenotype data from electronic medical records (EMRs) may provide a resource to assess the clinical relevance of rare variants.

OBJECTIVE

To determine the clinical phenotypes from EMRs for individuals with variants designated as pathogenic by expert review in arrhythmia susceptibility genes.

DESIGN, SETTING, AND PARTICIPANTS

This prospective cohort study included 2022 individuals recruited for nonantiarrhythmic drug exposure phenotypes from October 5, 2012, to September 30, 2013, for the Electronic Medical Records and Genomics Network Pharmacogenomics project from 7 US academic medical centers. Variants in SCN5A and KCNH2, disease genes for long QT and Brugada syndromes, were assessed for potential pathogenicity by 3 laboratories with ion channel expertise and by comparison with the ClinVar database. Relevant phenotypes were determined from EMRs, with data available from 2002 (or earlier for some sites) through September 10, 2014.

EXPOSURES

One or more variants designated as pathogenic in SCN5A or KCNH2.

MAIN OUTCOMES AND MEASURES

Arrhythmia or electrocardiographic (ECG) phenotypes defined by International Classification of Diseases, Ninth Revision (ICD-9) codes, ECG data, and manual EMR review.

RESULTS

Among 2022 study participants (median age, 61 years [interquartile range, 56-65 years]; 1118 [55%] female; 1491 [74%] white), a total of 122 rare (minor allele frequency <0.5%) nonsynonymous and splice-site variants in 2 arrhythmia susceptibility genes were identified in 223 individuals (11% of the study cohort). Forty-two variants in 63 participants were designated potentially pathogenic by at least 1 laboratory or ClinVar, with low concordance across laboratories (Cohen κ = 0.26). An ICD-9 code for arrhythmia was found in 11 of 63 (17%) variant carriers vs 264 of 1959 (13%) of those without variants (difference, +4%; 95% CI, -5% to +13%; P = .35). In the 1270 (63%) with ECGs, corrected QT intervals were not different in variant carriers vs those without (median, 429 vs 439 milliseconds; difference, -10 milliseconds; 95% CI, -16 to +3 milliseconds; P = .17). After manual review, 22 of 63 participants (35%) with designated variants had any ECG or arrhythmia phenotype, and only 2 had corrected QT interval longer than 500 milliseconds.

CONCLUSIONS AND RELEVANCE

Among laboratories experienced in genetic testing for cardiac arrhythmia disorders, there was low concordance in designating SCN5A and KCNH2 variants as pathogenic. In an unselected population, the putatively pathogenic genetic variants were not associated with an abnormal phenotype. These findings raise questions about the implications of notifying patients of incidental genetic findings.

TGFβ-Induced Actin Cytoskeleton Rearrangement in Podocytes Is Associated with Compensatory Adaptation of Mitochondrial Energy Metabolism

BACKGROUND/AIMS

In podocytes, the overexpression of TGFβ ligands and receptors during glomerulosclerosis could be a causal factor for injury induction and perpetuation in glomerular tufts. Mitochondrial dysfunction and oxidative stress are emerging as potential therapeutic targets in glomerular injury, and TGFβ has been shown to modulate mitochondrial metabolism in different cell types. This study aims at investigating the role of TGFβ in podocyte energy metabolism and cytoskeleton dynamics.

METHODS

Mitochondrial function and cytoskeleton dynamics were analyzed in TGFβ-treated WT and Smad2/3 double KO podocytes.

RESULTS

TGFβ treatment in podocytes induced a significant Smad-dependent increase of mitochondrial oxygen consumption rate (OCR). ATP content was unchanged and increased respiration was not associated with increased mitochondrial mass. Increased cellular reactive oxygen species induced by Smad-mediated TGFβ signaling were reverted by NADPH oxidase inhibitor apocynin. TGFβ treatment did not induce mitochondrial oxidative stress, and Smad2/3-dependent TGFβ signaling and increased mitochondrial OCR were found to be associated with actin cytoskeleton dynamics. The role of motor proteins myosin II and dynamin in TGFβ-induced actin polymerization was demonstrated by specific inhibition, resulting in actin stabilization and normalization of mitochondrial OCR.

CONCLUSION

TGFβ-induced rearrangements of actin cytoskeleton are controlled by Smad2/3 signaling pathways and coupled with the activation of mitochondrial ATP synthesis as bioenergetic adaptation to ATP consumption by ATP- and GTP-dependent motor proteins, myosin II and dynamin.

Penetrance of Hemochromatosis in HFE Genotypes Resulting in p.Cys282Tyr and p.[Cys282Tyr];[His63Asp] in the eMERGE Network

Hereditary hemochromatosis (HH) is a common autosomal-recessive disorder associated with pathogenic HFE variants, most commonly those resulting in p.Cys282Tyr and p.His63Asp. Recommendations on returning incidental findings of HFE variants in individuals undergoing genome-scale sequencing should be informed by penetrance estimates of HH in unselected samples. We used the eMERGE Network, a multicenter cohort with genotype data linked to electronic medical records, to estimate the diagnostic rate and clinical penetrance of HH in 98 individuals homozygous for the variant coding for HFE p.Cys282Tyr and 397 compound heterozygotes with variants resulting in p.[His63Asp];[Cys282Tyr]. The diagnostic rate of HH in males was 24.4% for p.Cys282Tyr homozygotes and 3.5% for compound heterozygotes (p < 0.001); in females, it was 14.0% for p.Cys282Tyr homozygotes and 2.3% for compound heterozygotes (p < 0.001). Only males showed differences across genotypes in transferrin saturation levels (100% of homozygotes versus 37.5% of compound heterozygotes with transferrin saturation > 50%; p = 0.003), serum ferritin levels (77.8% versus 33.3% with serum ferritin > 300 ng/ml; p = 0.006), and diabetes (44.7% versus 28.0%; p = 0.03). No differences were found in the prevalence of heart disease, arthritis, or liver disease, except for the rate of liver biopsy (10.9% versus 1.8% [p = 0.013] in males; 9.1% versus 2% [p = 0.035] in females). Given the higher rate of HH diagnosis than in prior studies, the high penetrance of iron overload, and the frequency of at-risk genotypes, in addition to other suggested actionable adult-onset genetic conditions, opportunistic screening should be considered for p.[Cys282Tyr];[Cys282Tyr] individuals with existing genomic data.

Practical considerations in genomic decision support: The eMERGE experience

BACKGROUND

Genomic medicine has the potential to improve care by tailoring treatments to the individual. There is consensus in the literature that pharmacogenomics (PGx) may be an ideal starting point for real-world implementation, due to the presence of well-characterized drug-gene interactions. Clinical Decision Support (CDS) is an ideal avenue by which to implement PGx at the bedside. Previous literature has established theoretical models for PGx CDS implementation and discussed a number of anticipated real-world challenges. However, work detailing actual PGx CDS implementation experiences has been limited. Anticipated challenges include data storage and management, system integration, physician acceptance, and more.

METHODS

In this study, we analyzed the experiences of ten members of the Electronic Medical Records and Genomics (eMERGE) Network, and one affiliate, in their attempts to implement PGx CDS. We examined the resulting PGx CDS system characteristics and conducted a survey to understand the unanticipated implementation challenges sites encountered.

RESULTS

Ten sites have successfully implemented at least one PGx CDS rule in the clinical setting. The majority of sites elected to create an Omic Ancillary System (OAS) to manage genetic and genomic data. All sites were able to adapt their existing CDS tools for PGx knowledge. The most common and impactful delays were not PGx-specific issues. Instead, they were general IT implementation problems, with top challenges including team coordination/communication and staffing. The challenges encountered caused a median total delay in system go-live of approximately two months.

CONCLUSIONS

These results suggest that barriers to PGx CDS implementations are generally surmountable. Moreover, PGx CDS implementation may not be any more difficult than other healthcare IT projects of similar scope, as the most significant delays encountered were not unique to genomic medicine. These are encouraging results for any institution considering implementing a PGx CDS tool, and for the advancement of genomic medicine.

A conceptual model for translating omic data into clinical action

Genomic, proteomic, epigenomic, and other “omic” data have the potential to enable precision medicine, also commonly referred to as personalized medicine. The volume and complexity of omic data are rapidly overwhelming human cognitive capacity, requiring innovative approaches to translate such data into patient care. Here, we outline a conceptual model for the application of omic data in the clinical context, called “the omic funnel.” This model parallels the classic “Data, Information, Knowledge, Wisdom pyramid” and adds context for how to move between each successive layer. Its goal is to allow informaticians, researchers, and clinicians to approach the problem of translating omic data from bench to bedside, by using discrete steps with clearly defined needs. Such an approach can facilitate the development of modular and interoperable software that can bring precision medicine into widespread practice.

Genetic variation in the HLA region is associated with susceptibility to herpes zoster

Herpes zoster, commonly referred to as shingles, is caused by the varicella zoster virus (VZV). VZV initially manifests as chicken pox, most commonly in childhood, can remain asymptomatically latent in nerve tissues for many years and often re-emerges as shingles. Although reactivation may be related to immune suppression, aging and female sex, most inter-individual variability in re-emergence risk has not been explained to date. We performed a genome-wide association analyses in 22 981 participants (2280 shingles cases) from the electronic Medical Records and Genomics Network. Using Cox survival and logistic regression, we identified a genomic region in the combined and European ancestry groups that has an age of onset effect reaching genome-wide significance (P>1.0 × 10⁻⁸). This region tags the non-coding gene HCP5 (HLA Complex P5) in the major histocompatibility complex. This gene is an endogenous retrovirus and likely influences viral activity through regulatory functions. Variants in this genetic region are known to be associated with delay in development of AIDS in people infected by HIV. Our study provides further suggestion that this region may have a critical role in viral suppression and could potentially harbor a clinically actionable variant for the shingles vaccine.

Design and anticipated outcomes of the eMERGE-PGx project: a multicenter pilot for preemptive pharmacogenomics in electronic health record systems

We describe here the design and initial implementation of the eMERGE-PGx project. eMERGE-PGx, a partnership of the eMERGE and PGRN consortia, has three objectives : 1) Deploy PGRNseq, a next-generation sequencing platform assessing sequence variation in 84 proposed pharmacogenes, in nearly 9,000 patients likely to be prescribed drugs of interest in a 1–3 year timeframe across several clinical sites; 2) Integrate well-established clinically-validated pharmacogenetic genotypes into the electronic health record with associated clinical decision support and assess process and clinical outcomes of implementation; and 3) Develop a repository of pharmacogenetic variants of unknown significance linked to a repository of EHR-based clinical phenotype data for ongoing pharmacogenomics discovery. We describe site-specific project implementation and anticipated products, including genetic variant and phenotype data repositories, novel variant association studies, clinical decision support modules, clinical and process outcomes, approaches to manage incidental findings, and patient and clinician education methods.

Phenome-wide association studies demonstrating pleiotropy of genetic variants within FTO with and without adjustment for body mass index

Phenome-wide association studies (PheWAS) have demonstrated utility in validating genetic associations derived from traditional genetic studies as well as identifying novel genetic associations. Here we used an electronic health record (EHR)-based PheWAS to explore pleiotropy of genetic variants in the fat mass and obesity associated gene (FTO), some of which have been previously associated with obesity and type 2 diabetes (T2D). We used a population of 10,487 individuals of European ancestry with genome-wide genotyping from the Electronic Medical Records and Genomics (eMERGE) Network and another population of 13,711 individuals of European ancestry from the BioVU DNA biobank at Vanderbilt genotyped using Illumina HumanExome BeadChip. A meta-analysis of the two study populations replicated the well-described associations between FTO variants and obesity (odds ratio [OR] = 1.25, 95% Confidence Interval = 1.11-1.24, p = 2.10 × 10(-9)) and FTO variants and T2D (OR = 1.14, 95% CI = 1.08-1.21, p = 2.34 × 10(-6)). The meta-analysis also demonstrated that FTO variant rs8050136 was significantly associated with sleep apnea (OR = 1.14, 95% CI = 1.07-1.22, p = 3.33 × 10(-5)); however, the association was attenuated after adjustment for body mass index (BMI). Novel phenotype associations with obesity-associated FTO variants included fibrocystic breast disease (rs9941349, OR = 0.81, 95% CI = 0.74-0.91, p = 5.41 × 10(-5)) and trends toward associations with non-alcoholic liver disease and gram-positive bacterial infections. FTO variants not associated with obesity demonstrated other potential disease associations including non-inflammatory disorders of the cervix and chronic periodontitis. These results suggest that genetic variants in FTO may have pleiotropic associations, some of which are not mediated by obesity.

Smad2-dependent downregulation of miR-30 is required for TGF-β-induced apoptosis in podocytes

Transforming growth factors beta (TGF-β) are multi-functional cytokines capable of inducing apoptosis in epithelial cells, including glomerular podocytes. We and others have previously shown that podocyte-selective genetic deletion of the microRNA (miR)-processing enzyme, Dicer, caused glomerulosclerosis that was associated with podocyte apoptosis, and the miR-30 family was implicated in the process. Here, we report that apoptosis-associated genes were highly enriched among the predicted targets of miR-30 when compared with randomly selected miRs (26% vs. 4.5 ± 2.1%) or with the known TGF-β-regulated miR-192 (6%), miR-216a (5.1%), and miR-217 (0%). miR-30 family members were abundantly expressed in podocytes in normal mice but were downregulated in albumin/TGF-β transgenic mice with podocyte apoptosis and glomerulosclerosis. In vitro, TGF-β downregulated miR-30s in wildtype and Smad3-deficient, but not Smad2- or Smad2/Smad3-deficient, podocytes. The TGF-β-induced activation of caspase 3 and an increase in TUNEL-positive nuclei were significantly inhibited by the lentivirus-mediated overexpression of miR-30d, but not by a scrambled control miR, in podocytes. TGF-β stimulated the phosphorylation of pro-apoptotic p53 in podocytes with lentiviral expression of a scrambled miR, but not in podocytes expressing miR-30d. In contrast, miR-30d had no effect on the phosphorylation of pro-apoptotic p38 MAP kinase induced by TGF-β. Thus, we report that Smad2-dependent inhibition of miR-30s in podocytes is required for the activation of p53 and the induction of apoptosis by TGF-β. These results demonstrate a novel functional role for miR-30 in podocyte survival and indicate that the loss of miR-30 survival signaling is a novel and specific mechanism of TGF-β-induced podocyte apoptosis during glomerulosclerosis. We propose the therapeutic replacement of miR-30 as a novel strategy to prevent the podocyte apoptosis that is characteristic of progressive glomerular diseases.

Genome-wide association analysis of blood-pressure traits in African-ancestry individuals reveals common associated genes in African and non-African populations

High blood pressure (BP) is more prevalent and contributes to more severe manifestations of cardiovascular disease (CVD) in African Americans than in any other United States ethnic group. Several small African-ancestry (AA) BP genome-wide association studies (GWASs) have been published, but their findings have failed to replicate to date. We report on a large AA BP GWAS meta-analysis that includes 29,378 individuals from 19 discovery cohorts and subsequent replication in additional samples of AA (n = 10,386), European ancestry (EA) (n = 69,395), and East Asian ancestry (n = 19,601). Five loci (EVX1-HOXA, ULK4, RSPO3, PLEKHG1, and SOX6) reached genome-wide significance (p < 1.0 × 10(-8)) for either systolic or diastolic BP in a transethnic meta-analysis after correction for multiple testing. Three of these BP loci (EVX1-HOXA, RSPO3, and PLEKHG1) lack previous associations with BP. We also identified one independent signal in a known BP locus (SOX6) and provide evidence for fine mapping in four additional validated BP loci. We also demonstrate that validated EA BP GWAS loci, considered jointly, show significant effects in AA samples. Consequently, these findings suggest that BP loci might have universal effects across studied populations, demonstrating that multiethnic samples are an essential component in identifying, fine mapping, and understanding their trait variability.

A collaborative approach to developing an electronic health record phenotyping algorithm for drug-induced liver injury

OBJECTIVE

To describe a collaborative approach for developing an electronic health record (EHR) phenotyping algorithm for drug-induced liver injury (DILI).

METHODS

We analyzed types and causes of differences in DILI case definitions provided by two institutions-Columbia University and Mayo Clinic; harmonized two EHR phenotyping algorithms; and assessed the performance, measured by sensitivity, specificity, positive predictive value, and negative predictive value, of the resulting algorithm at three institutions except that sensitivity was measured only at Columbia University.

RESULTS

Although these sites had the same case definition, their phenotyping methods differed by selection of liver injury diagnoses, inclusion of drugs cited in DILI cases, laboratory tests assessed, laboratory thresholds for liver injury, exclusion criteria, and approaches to validating phenotypes. We reached consensus on a DILI phenotyping algorithm and implemented it at three institutions. The algorithm was adapted locally to account for differences in populations and data access. Implementations collectively yielded 117 algorithm-selected cases and 23 confirmed true positive cases.

DISCUSSION

Phenotyping for rare conditions benefits significantly from pooling data across institutions. Despite the heterogeneity of EHRs and varied algorithm implementations, we demonstrated the portability of this algorithm across three institutions. The performance of this algorithm for identifying DILI was comparable with other computerized approaches to identify adverse drug events.

CONCLUSIONS

Phenotyping algorithms developed for rare and complex conditions are likely to require adaptive implementation at multiple institutions. Better approaches are also needed to share algorithms. Early agreement on goals, data sources, and validation methods may improve the portability of the algorithms.

Transforming growth factor-β, bioenergetics, and mitochondria in renal disease

The transforming growth factor-β (TGF-β) family comprises more than 30 family members that are structurally related secreted dimeric cytokines, including TGF-β, activins, and bone morphogenetic proteins/growth and differentiation factors. TGF-β are pluripotent regulators of cell proliferation, differentiation, apoptosis, migration, and adhesion of many different cell types. TGF-β pathways are highly evolutionarily conserved and control embryogenesis, tissue repair, and tissue homeostasis in invertebrates and vertebrates. Aberrations in TGF-β activity and signaling underlie a broad spectrum of developmental disorders and major pathologies in human beings, including cancer, fibrosis, and autoimmune diseases. Recent observations have indicated an emerging role for TGF-β in the regulation of mitochondrial bioenergetics and oxidative stress responses characteristic of chronic degenerative diseases and aging. Conversely, energy and metabolic sensory pathways cross-regulate mediators of TGF-β signaling. Here, we review TGF-β and regulation of bioenergetic and mitochondrial functions, including energy and oxidant metabolism and apoptotic cell death, as well as their emerging relevance in renal biology and disease.

BAMBI regulates angiogenesis and endothelial homeostasis through modulation of alternative TGFβ signaling

BACKGROUND

BAMBI is a type I TGFβ receptor antagonist, whose in vivo function remains unclear, as BAMBI(-/-) mice lack an obvious phenotype.

METHODOLOGY/PRINCIPAL FINDINGS

Identifying BAMBI's functions requires identification of cell-specific expression of BAMBI. By immunohistology we found BAMBI expression restricted to endothelial cells and by electron microscopy BAMBI(-/-) mice showed prominent and swollen endothelial cells in myocardial and glomerular capillaries. In endothelial cells over-expression of BAMBI reduced, whereas knock-down enhanced capillary growth and migration in response to TGFβ. In vivo angiogenesis was enhanced in matrigel implants and in glomerular hypertrophy after unilateral nephrectomy in BAMBI(-/-) compared to BAMBI(+/+) mice consistent with an endothelial phenotype for BAMBI(-/-) mice. BAMBI's mechanism of action in endothelial cells was examined by canonical and alternative TGFβ signaling in HUVEC with over-expression or knock-down of BAMBI. BAMBI knockdown enhanced basal and TGFβ stimulated SMAD1/5 and ERK1/2 phosphorylation, while over-expression prevented both.

CONCLUSIONS/SIGNIFICANCE

Thus we provide a first description of a vascular phenotype for BAMBI(-/-) mice, and provide in vitro and in vivo evidence that BAMBI contributes to endothelial and vascular homeostasis. Further, we demonstrate that in endothelial cells BAMBI interferes with alternative TGFβ signaling, most likely through the ALK 1 receptor, which may explain the phenotype observed in BAMBI(-/-) mice. This newly described role for BAMBI in regulating endothelial function has potential implications for understanding and treating vascular disease and tumor neo-angiogenesis.

BAMBI is expressed in endothelial cells and is regulated by lysosomal/autolysosomal degradation

BACKGROUND

BAMBI (BMP and Activin Membrane Bound Inhibitor) is considered to influence TGFβ and Wnt signaling, and thereby fibrosis. Surprisingly data on cell type-specific expression of BAMBI are not available. We therefore examined the localization, gene regulation, and protein turnover of BAMBI in kidneys.

METHODOLOGY/PRINCIPAL FINDINGS

By immunofluorescence microscopy and by mRNA expression, BAMBI is restricted to endothelial cells of the glomerular and some peritubular capillaries and of arteries and veins in both murine and human kidneys. TGFβ upregulated mRNA of BAMBI in murine glomerular endothelial cells (mGEC). LPS did not downregulate mRNA for BAMBI in mGEC or in HUVECs. BAMBI mRNA had a half-life of only 60 minutes and was stabilized by cycloheximide, indicating post-transcriptional regulation due to AU-rich elements, which we identified in the 3' untranslated sequence of both the human and murine BAMBI gene. BAMBI protein turnover was studied in HUVECs with BAMBI overexpression using a lentiviral system. Serum starvation as an inducer of autophagy caused marked BAMBI degradation, which could be totally prevented by inhibition of lysosomal and autolysosomal degradation with bafilomycin, and partially by inhibition of autophagy with 3-methyladenine, but not by proteasomal inhibitors. Rapamycin activates autophagy by inhibiting TOR, and resulted in BAMBI protein degradation. Both serum starvation and rapamycin increased the conversion of the autophagy marker LC3 from LC3-I to LC3-II and also enhanced co-staining for BAMBI and LC3 in autolysosomal vesicles.

CONCLUSIONS/SIGNIFICANCE

1. BAMBI localizes to endothelial cells in the kidney and to HUVECs. 2. BAMBI mRNA is regulated by post-transcriptional mechanisms. 3. BAMBI protein is regulated by lysosomal and autolysosomal degradation. The endothelial localization and the quick turnover of BAMBI may indicate novel, yet to be defined functions of this modulator for TGFβ and Wnt protein actions in the renal vascular endothelium in health and disease.

HGF upregulation contributes to angiogenesis in mice with keratinocyte-specific Smad2 deletion

TGF-β signaling can promote tumor formation and development or suppress it, depending on the cellular context and tumor stage. A potential target of this dual effect of TGF-β is HGF, as TGF-β can inhibit or promote its expression, although the mechanisms underlying this are largely unknown. In the present study, we found that mice with keratinocyte-specific deletion of the TGF-β signaling mediator Smad2 (referred to herein as K5.Smad2(-/-) mice), which have increased susceptibility to squamous cell carcinomas (SCCs), exhibited angiogenesis associated with epithelial overexpression of HGF and endothelial activation of the HGF receptor c-Met. Application of a c-Met inhibitor abrogated angiogenesis, suggesting that HGF overexpression plays a major role in angiogenesis associated with epithelial Smad2 loss. On the Hgf promoter, Smad2 was mainly associated with transcriptional corepressors, whereas Smad4 was mainly associated with the transcriptional coactivator CREB-binding protein (CBP/p300). Smad2 loss caused increased binding of Smad4 and CBP/p300 to the Hgf promoter. Consistent with this, knocking down Smad2 in human keratinocytes caused increased levels of HGF, which were abrogated by concomitant knockdown of Smad3 and Smad4. Importantly, the incidence of HGF-positive human SCC was high in cases with Smad2 loss and lower when Smad4 was also lost. We therefore conclude that Smad2 loss causes HGF upregulation via loss of Smad2-mediated transcriptional repression and enhanced Smad3/4-mediated transactivation. Since Smad2 is often downregulated in human SCCs, our data suggest a therapeutic strategy of blocking HGF/c-Met activation for Smad2-deficient SCCs.

Human and murine kidneys show gender- and species-specific gene expression differences in response to injury

The incidence of End Stage Renal Disease (ESRD) is approximately 50% higher in men than women. In order to understand the molecular basis of this gender disparity, we examined sex specific gene expression patterns in control and diseased, human and murine kidney samples. Using the Affymetrix platform we performed comprehensive gene expression analysis on 42 microdissected human kidney samples (glomeruli and tubules). We identified 67 genes with gender biased expression in healthy human kidneys and 24 transcripts in diseased male and female human kidneys. Similar analysis performed in mice using male and female control and doxorubicin induced nephrotic syndrome kidneys identified significantly larger number of differentially expressed transcripts. The majority of genes showing gender biased expression either in diseased human and murine kidneys were different from those differentially expressed in healthy kidneys. Only 9 sexually dimorphic transcripts were common to healthy human and murine kidneys and five showed differential regulation in both human and murine diseased kidneys. In humans, sex biased genes showed statistical enrichment only to sex chromosomes while in mice they were enriched to sex chromosomes and various autosomes. Thus we present a comprehensive analysis of gender biased genes in the kidney. We show that sexually dimorphic genes in the kidney show species specific regulation. Our results also indicate that male and female kidneys respond differently to injury. These studies could provide the basis for the development of new treatment strategies for men and women with kidney disease.

Keratinocyte-specific Smad2 ablation results in increased epithelial-mesenchymal transition during skin cancer formation and progression

TGF-beta and its signaling mediators, Smad2, -3, and -4, are involved with tumor suppression and promotion functions. Smad4-/- mouse epidermis develops spontaneous skin squamous cell carcinomas (SCCs), and Smad3-/- mice are resistant to carcinogen-induced skin cancer; however, the role of Smad2 in skin carcinogenesis has not been explored. In the present study, we found that Smad2 and Smad4, but not Smad3, were frequently lost in human SCCs. Mice with keratinocyte-specific Smad2 deletion exhibited accelerated formation and malignant progression of chemically induced skin tumors compared with WT mice. Consistent with the loss of Smad2 in poorly differentiated human SCCs, Smad2-/- tumors were poorly differentiated and underwent epithelial-mesenchymal transition (EMT) prior to spontaneous Smad4 loss. Reduced E-cadherin and activation of its transcriptional repressor Snail were also found in Smad2-/- mouse epidermis and occurred more frequently in Smad2-negative human SCCs than in Smad2-positive SCCs. Knocking down Snail abrogated Smad2 loss-associated EMT, suggesting that Snail upregulation is a major mediator of Smad2 loss-associated EMT. Furthermore, Smad2 loss led to a significant increase in Smad4 binding to the Snail promoter, and knocking down either Smad3 or Smad4 in keratinocytes abrogated Smad2 loss-associated Snail overexpression. Our data suggest that enhanced Smad3/Smad4-mediated Snail transcription contributed to Smad2 loss-associated EMT during skin carcinogenesis.

Decorin deficiency enhances progressive nephropathy in diabetic mice

Decorin, a proteoglycan that inhibits active transforming growth factor-beta, is increased in diabetic nephropathy; however, its functional significance is unclear. In this study, we used low-dose streptozotocin to induce type 1 diabetes in wild-type (C57BL/6J Dcn(+/+)), Dcn(-/-), and Dcn(+/-) mice and studied the mice for up to 1 year of diabetes. Decorin gene dose had no effect on severity of diabetes; however, the Dcn(-/-) diabetic mice died significantly earlier than nondiabetic controls (57 versus 7.3% mortality). In contrast to wild-type diabetic mice, which failed to develop significant nephropathy, the Dcn(-/-) diabetic mice developed a significant increase in albuminuria and plasma creatinine and a concurrent decrease in circulating adiponectin levels. Interestingly, adiponectin levels at 6 months of diabetes were predictive of mortality in diabetic mice. Dcn(-/-) diabetic mice exhibited advanced glomerular lesions, including diffuse mesangial matrix accumulation and fibrin cap formation. By immunohistochemistry, Dcn(-/-) diabetic mice exhibited significant increases in glomerular transforming growth factor-beta, type I collagen, macrophage infiltration, and Nox4. We conclude that decorin is a natural protective factor against diabetic nephropathy and that the Dcn(-/-) diabetic mouse is a useful new model of progressive diabetic nephropathy.

Genome-wide molecular profiles of HCV-induced dysplasia and hepatocellular carcinoma

Although HCC is the third-leading cause of cancer-related deaths worldwide, there is only an elemental understanding of its molecular pathogenesis. In western countries, HCV infection is the main etiology underlying this cancer's accelerating incidence. To characterize the molecular events of the hepatocarcinogenic process, and to identify new biomarkers for early HCC, the gene expression profiles of 75 tissue samples were analyzed representing the stepwise carcinogenic process from preneoplastic lesions (cirrhosis and dysplasia) to HCC, including 4 neoplastic stages (very early HCC to metastatic tumors) from patients with HCV infection. We identified gene signatures that accurately reflect the pathological progression of disease at each stage. Eight genes distinguish between control and cirrhosis, 24 between cirrhosis and dysplasia, 93 between dysplasia and early HCC, and 9 between early and advanced HCC. Using quantitative real-time reverse-transcription PCR, we validated several novel molecular tissue markers for early HCC diagnosis, specifically induction of abnormal spindle-like, microcephaly-associated protein, hyaluronan-mediated motility receptor, primase 1, erythropoietin, and neuregulin 1. In addition, pathway analysis revealed dysregulation of the Notch and Toll-like receptor pathways in cirrhosis, followed by deregulation of several components of the Jak/STAT pathway in early carcinogenesis, then upregulation of genes involved in DNA replication and repair and cell cycle in late cancerous stages.

CONCLUSION

These findings provide a comprehensive molecular portrait of genomic changes in progressive HCV-related HCC.

TGF-beta isoform signaling regulates secondary transition and mesenchymal-induced endocrine development in the embryonic mouse pancreas

Transforming growth factor-beta (TGF-beta) superfamily signaling has been implicated in many developmental processes, including pancreatic development. Previous studies are conflicting with regard to an exact role for TGF-beta signaling in various aspects of pancreatic organogenesis. Here we have investigated the role of TGF-beta isoform signaling in embryonic pancreas differentiation and lineage selection. The TGF-beta isoform receptors (RI, RII and ALK1) were localized mainly to both the pancreatic epithelium and mesenchyme at early stages of development, but then with increasing age localized to the pancreatic islets and ducts. To determine the specific role of TGF-beta isoforms, we functionally inactivated TGF-beta signaling at different points in the signaling cascade. Disruption of TGF-beta signaling at the receptor level using mice overexpressing the dominant-negative TGF-beta type II receptor showed an increase in endocrine precursors and proliferating endocrine cells, with an abnormal accumulation of endocrine cells around the developing ducts of mid-late stage embryonic pancreas. This pattern suggested that TGF-beta isoform signaling may suppress the origination of secondary transition endocrine cells from the ducts. Secondly, TGF-beta isoform ligand inhibition with neutralizing antibody in pancreatic organ culture also led to an increase in the number of endocrine-positive cells. Thirdly, hybrid mix-and-match in vitro recombinations of transgenic pancreatic mesenchyme and wild-type epithelium also led to increased endocrine cell differentiation, but with different patterns depending on the directionality of the epithelial-mesenchymal signaling. Together these results suggest that TGF-beta signaling is important for restraining the growth and differentiation of pancreatic epithelial cells, particularly away from the endocrine lineage. Inhibition of TGF-beta signaling in the embryonic period may thus allow pancreatic epithelial cells to progress towards the endocrine lineage unchecked, particularly as part of the secondary transition of pancreatic endocrine cell development. TGF-beta RII in the ducts and islets may normally serve to downregulate the production of beta cells from embryonic ducts.

Exploiting noise in array CGH data to improve detection of DNA copy number change

Developing effective methods for analyzing array-CGH data to detect chromosomal aberrations is very important for the diagnosis of pathogenesis of cancer and other diseases. Current analysis methods, being largely based on smoothing and/or segmentation, are not quite capable of detecting both the aberration regions and the boundary break points very accurately. Furthermore, when evaluating the accuracy of an algorithm for analyzing array-CGH data, it is commonly assumed that noise in the data follows normal distribution. A fundamental question is whether noise in array-CGH is indeed Gaussian, and if not, can one exploit the characteristics of noise to develop novel analysis methods that are capable of detecting accurately the aberration regions as well as the boundary break points simultaneously? By analyzing bacterial artificial chromosomes (BACs) arrays with an average 1 mb resolution, 19 k oligo arrays with the average probe spacing <100 kb and 385 k oligo arrays with the average probe spacing of about 6 kb, we show that when there are aberrations, noise in all three types of arrays is highly non-Gaussian and possesses long-range spatial correlations, and that such noise leads to worse performance of existing methods for detecting aberrations in array-CGH than the Gaussian noise case. We further develop a novel method, which has optimally exploited the character of the noise, and is capable of identifying both aberration regions as well as the boundary break points very accurately. Finally, we propose a new concept, posteriori signal-to-noise ratio (p-SNR), to assign certain confidence level to an aberration region and boundaries detected.

Novel Regulators of Kidney Development from the Tips of the Ureteric Bud

Mammalian nephrogenesis depends on the interaction between the ureteric bud and the metanephric mesenchyme. As the ureteric bud undergoes branching and segmentation, the stalks differentiate into the collecting system of the mature kidney, while the tip cells interact with the adjacent cells of the metanephric mesenchyme, inducing their conversion into nephrons. This induction is mediated by secreted factors. For identifying novel mediators, the tips of the ureteric tree were isolated and microarray analyses were performed using manually refined, multistep gene ontology annotations. For identifying conserved factors, two databases were developed, one from mouse E12.5 and one from rat E13.5 ureteric buds. The overlap of mouse and rat data sets yielded 20 different transcripts that were enriched in the ureteric bud compared with metanephric mesenchyme and predicted to code for secreted proteins. Real-time reverse transcriptase-PCR and in situ hybridization confirmed these identifications. One of the genes that was highly specific to the ureteric bud tip was cytokine-like factor 1 (CLF-1). Recombinant CLF-1 in complex with its physiologic ligand, cardiotrophin-like cytokine (CLC), triggered phosphorylation of signal transducer and activator of transcription 3 in mesenchyme, a pathway characteristic of mesenchymal-to-epithelial conversion. Indeed, when applied to isolated rat metanephric mesenchyme, CLF-1/CLC (3 nM) induced mature nephron structures expressing glomerular and tubular markers. These results underline the power of this first comprehensive gene expression analysis of the ureteric bud tip to identify bioactive molecules.

Defects in TGF-beta signaling overcome senescence of mouse keratinocytes expressing v-Ha-ras

Previous studies have shown that TGFbeta1 expression is upregulated in mouse keratinocytes infected with a v-rasHa retrovirus, although the functional significance of this has not been clear. Here we show that v-rasHa retrovirus transduced primary mouse keratinocytes undergo hyperproliferation followed by a TGFbeta1 dependent G1 growth arrest and senescence. The growth arrest is accompanied by a 15-fold increase in total TGFbeta1 secreted and a fourfold increase in secreted active TGFbeta1. When cultured in the presence of a neutralizing antibody to TGFbeta1, the senescence response is suppressed. Levels of the TGFbeta1 target p15ink4b increase during senescence as does association of this kinase inhibitor with cyclinD/cdk4 complexes. However, p16ink4a, p53 and p19ARF expression also increase during senescence. Genetic analysis shows that TGFbeta1 null and dominant negative TbetaBRII expressing v-rasHa keratinocytes resist the G1 growth arrest and do not senescence. This resistance is associated with low expression of p15ink4b and p16ink4a, constitutive Rb phosphorylation and high levels of cdk4 and cdk2 kinase activity. In contrast, inactivation of TGFbetabeta1 secretion or response does not block the induction of p53 and p19ARF, but the level of p21waf1, a p53 target gene, is reduced in cyclin D/cdk4 and cyclin E/cdk2 complexes. Thus, although multiple senescence pathways are activated in response to a ras oncogene, inactivation of TGFbeta1 secretion or response is sufficient to block the senescence program. Since v-rasHa transduced TGFbeta1-/- keratinocytes form squamous cell carcinomas following skin grafting, these results suggest that in mouse keratinocytes, defects in TGFbeta1 signaling accelerate malignant progression by overcoming oncogene induced replicative senescence.

Defects in transforming growth factor-beta signaling cooperate with a Ras oncogene to cause rapid aneuploidy and malignant transformation of mouse keratinocytes

Genetic inactivation of the transforming growth factor-beta (TGF-beta) signaling pathway can accelerate tumor progression in the mouse epidermal model of multistage carcinogenesis. By using an in vitro model of keratinocyte transformation that parallels in vivo malignant conversion to squamous cell carcinoma, we show that v-ras(Ha) transduced primary TGF-beta1-/- keratinocytes and keratinocytes expressing a TGF-beta type II dominant-negative receptor transgene have significantly higher frequencies of spontaneous transformation than control genotypes. Malignant transformation in the TGF-beta1-/- keratinocytes is preceded by aneuploidy and accumulation of chromosomal aberrations. Similarly, transient inactivation of TGF-beta signaling with a type II dominant-negative receptor adenovirus causes rapid changes in ploidy. Exogenous TGF-beta1 can suppress aneuploidy, chromosome breaks, and malignant transformation of the TGF-beta1-/- keratinocytes at concentrations that do not significantly arrest cell proliferation. These results point to genomic instability as a mechanism by which defects in TGF-beta signaling could accelerate tumor progression in mouse multistage carcinogenesis.