Clinical Presentation of Blastomycosis is Associated With Infecting Species, Not Host Genotype.

Objective: To determine if host genetics may be a risk factor for severe blastomycosis.Design: A cohort of patients who had contracted blastomycosis underwent targeted SNP (single nucleotide polymorphism) genotyping. The genetics of these patients were compared to a set of age and gender-matched controls and between patients with severe versus mild to moderate blastomycosis.Setting: The Marshfield Clinic Health System in central and northern WisconsinParticipants: Patients with a diagnosis of blastomycosis prior to 2017 were contacted for enrollment in this study. A phone hotline was also set up to allow interested participants from outside the Marshfield Clinic Health System to request enrollment.Methods: SNP frequency was assessed for significant differences between the patient cohort and controls and between patients with severe versus mild to moderate blastomycosis. We also tested the effect of Blastomyces species identified in clinical isolates on disease symptoms and severity.Results: No significant differences were found in SNP frequency between cases and controls or between those with severe or mild to moderate blastomycosis. We did detect significant differences in symptom frequency and disease severity by Blastomyces species.Conclusions: Our study did not identify any genetic risk factors for blastomycosis. Instead, the species of Blastomyces causing the infection had a significant effect on disease severity.

A gene-based recessive diplotype exome scan discovers FGF6, a novel hepcidin-regulating iron-metabolism gene.

Standard analyses applied to genome-wide association data are well designed to detect additive effects of moderate strength. However, the power for standard genome-wide association study (GWAS) analyses to identify effects from recessive diplotypes is not typically high. We proposed and conducted a gene-based compound heterozygosity test to reveal additional genes underlying complex diseases. With this approach applied to iron overload, a strong association signal was identified between the fibroblast growth factor-encoding gene, FGF6, and hemochromatosis in the central Wisconsin population. Functional validation showed that fibroblast growth factor 6 protein (FGF-6) regulates iron homeostasis and induces transcriptional regulation of hepcidin. Moreover, specific identified FGF6 variants differentially impact iron metabolism. In addition, FGF6 downregulation correlated with iron-metabolism dysfunction in systemic sclerosis and cancer cells. Using the recessive diplotype approach revealed a novel susceptibility hemochromatosis gene and has extended our understanding of the mechanisms involved in iron metabolism.

Public Attitudes toward Consent and Data Sharing in Biobank Research: A Large Multi-site Experimental Survey in the US.

Individuals participating in biobanks and other large research projects are increasingly asked to provide broad consent for open-ended research use and widespread sharing of their biosamples and data. We assessed willingness to participate in a biobank using different consent and data sharing models, hypothesizing that willingness would be higher under more restrictive scenarios. Perceived benefits, concerns, and information needs were also assessed. In this experimental survey, individuals from 11 US healthcare systems in the Electronic Medical Records and Genomics (eMERGE) Network were randomly allocated to one of three hypothetical scenarios: tiered consent and controlled data sharing; broad consent and controlled data sharing; or broad consent and open data sharing. Of 82,328 eligible individuals, exactly 13,000 (15.8%) completed the survey. Overall, 66% (95% CI: 63%-69%) of population-weighted respondents stated they would be willing to participate in a biobank; willingness and attitudes did not differ between respondents in the three scenarios. Willingness to participate was associated with self-identified white race, higher educational attainment, lower religiosity, perceiving more research benefits, fewer concerns, and fewer information needs. Most (86%, CI: 84%-87%) participants would want to know what would happen if a researcher misused their health information; fewer (51%, CI: 47%-55%) would worry about their privacy. The concern that the use of broad consent and open data sharing could adversely affect participant recruitment is not supported by these findings. Addressing potential participants' concerns and information needs and building trust and relationships with communities may increase acceptance of broad consent and wide data sharing in biobank research.

Probing the Virtual Proteome to Identify Novel Disease Biomarkers.

BACKGROUND: Proteomic approaches allow measurement of thousands of proteins in a single specimen, which can accelerate biomarker discovery. However, applying these technologies to massive biobanks is not currently feasible because of the practical barriers and costs of implementing such assays at scale. To overcome these challenges, we used a "virtual proteomic" approach, linking genetically predicted protein levels to clinical diagnoses in >40 000 individuals. METHODS: We used genome-wide association data from the Framingham Heart Study (n=759) to construct genetic predictors for 1129 plasma protein levels. We validated the genetic predictors for 268 proteins and used them to compute predicted protein levels in 41 288 genotyped individuals in the Electronic Medical Records and Genomics (eMERGE) cohort. We tested associations for each predicted protein with 1128 clinical phenotypes. Lead associations were validated with directly measured protein levels and either low-density lipoprotein cholesterol or subclinical atherosclerosis in the MDCS (Malmö Diet and Cancer Study; n=651). RESULTS: In the virtual proteomic analysis in eMERGE, 55 proteins were associated with 89 distinct diagnoses at a false discovery rate q<0.1. Among these, 13 associations involved lipid (n=7) or atherosclerosis (n=6) phenotypes. We tested each association for validation in MDCS using directly measured protein levels. At Bonferroni-adjusted significance thresholds, levels of apolipoprotein E isoforms were associated with hyperlipidemia, and circulating C-type lectin domain family 1 member B and platelet-derived growth factor receptor-ß predicted subclinical atherosclerosis. Odds ratios for carotid atherosclerosis were 1.31 (95% CI, 1.08-1.58; P=0.006) per 1-SD increment in C-type lectin domain family 1 member B and 0.79 (0.66-0.94; P=0.008) per 1-SD increment in platelet-derived growth factor receptor-ß. CONCLUSIONS: We demonstrate a biomarker discovery paradigm to identify candidate biomarkers of cardiovascular and other diseases.

Applying family analyses to electronic health records to facilitate genetic research.

Motivation: Pedigree analysis is a longstanding and powerful approach to gain insight into the underlying genetic factors in human health, but identifying, recruiting and genotyping families can be difficult, time consuming and costly. Development of high throughput methods to identify families and foster downstream analyses are necessary. Results: This paper describes simple methods that allowed us to identify 173 368 family pedigrees with high probability using basic demographic data available in most electronic health records (EHRs). We further developed and validate a novel statistical method that uses EHR data to identify families more likely to have a major genetic component to their diseases risk. Lastly, we showed that incorporating EHR-linked family data into genetic association testing may provide added power for genetic mapping without additional recruitment or genotyping. The totality of these results suggests that EHR-linked families can enable classical genetic analyses in a high-throughput manner. Availability and implementation: Pseudocode is provided as supplementary information. Contact: HEBBRING.SCOTT@marshfieldresearch.org. Supplementary information: Supplementary data are available at Bioinformatics online.

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.

Biology-Driven Gene-Gene Interaction Analysis of Age-Related Cataract in the eMERGE Network.

Bioinformatics approaches to examine gene-gene models provide a means to discover interactions between multiple genes that underlie complex disease. Extensive computational demands and adjusting for multiple testing make uncovering genetic interactions a challenge. Here, we address these issues using our knowledge-driven filtering method, Biofilter, to identify putative single nucleotide polymorphism (SNP) interaction models for cataract susceptibility, thereby reducing the number of models for analysis. Models were evaluated in 3,377 European Americans (1,185 controls, 2,192 cases) from the Marshfield Clinic, a study site of the Electronic Medical Records and Genomics (eMERGE) Network, using logistic regression. All statistically significant models from the Marshfield Clinic were then evaluated in an independent dataset of 4,311 individuals (742 controls, 3,569 cases), using independent samples from additional study sites in the eMERGE Network: Mayo Clinic, Group Health/University of Washington, Vanderbilt University Medical Center, and Geisinger Health System. Eighty-three SNP-SNP models replicated in the independent dataset at likelihood ratio test P < 0.05. Among the most significant replicating models was rs12597188 (intron of CDH1)-rs11564445 (intron of CTNNB1). These genes are known to be involved in processes that include: cell-to-cell adhesion signaling, cell-cell junction organization, and cell-cell communication. Further Biofilter analysis of all replicating models revealed a number of common functions among the genes harboring the 83 replicating SNP-SNP models, which included signal transduction and PI3K-Akt signaling pathway. These findings demonstrate the utility of Biofilter as a biology-driven method, applicable for any genome-wide association study dataset.

Analysis of CYP1B1 in pediatric and adult glaucoma and other ocular phenotypes.

PURPOSE: The CYP1B1 gene encodes an enzyme that is a member of the cytochrome P450 superfamily. Mutations in CYP1B1 have been mainly reported in recessive pediatric ocular phenotypes, such as primary congenital glaucoma (PCG) and congenital glaucoma with anterior segment dysgenesis (CG with ASD), with some likely pathogenic variants also identified in families affected with adult-onset primary open angle glaucoma (POAG). METHODS: We examined CYP1B1 in 158 pediatric patients affected with PCG (eight), CG with ASD (22), CG with other developmental ocular disorders (11), juvenile glaucoma with or without additional ocular anomalies (26), and ASD or other developmental ocular conditions without glaucoma (91); in addition, a large cohort of adult patients with POAG (193) and POAG-negative controls (288) was examined. RESULTS: Recessive pathogenic variants in CYP1B1 were identified in two PCG pedigrees, three cases with CG and ASD, and two families with CG and other ocular defects, such as sclerocornea in one patient and microphthalmia in another individual; neither sclerocornea nor microphthalmia has been previously associated with CYP1B1. Most of the identified causative mutations are new occurrences of previously reported pathogenic alleles with two novel variants identified: a c.1325delC, p.(Pro442Glnfs*15) frameshift allele in a family with PCG and a c.157G>A, p.(Gly53Ser) variant identified in a proband with CG, Peters anomaly, and microphthalmia. Analysis of the family history in the CYP1B1-positive families revealed POAG in confirmed or presumed heterozygous relatives in one family with PCG and two families with ASD/CG; POAG was associated with the c.1064_1076del, p.(Arg355Hisfs*69) allele in two of these pedigrees. Screening of an unrelated POAG cohort identified the same c.1064_1076del heterozygous allele in one individual with sporadic POAG but not in age- and ethnicity-matched POAG-negative individuals. Overall, there was no significant enrichment for mutant alleles in CYP1B1 within the POAG cases compared to the controls. CONCLUSIONS: In summary, these data expand the mutational and phenotypic spectra of CYP1B1 to include two novel alleles and additional developmental ocular phenotypes. The contribution of CYP1B1 to POAG is less clear, but loss-of-function variants in CYP1B1, especially c.1064_1076del, p.(Arg355Hisfs*69), may be associated with an increased risk for POAG.

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.

Use of an electronic medical record to create the marshfield clinic twin/multiple birth cohort.

Population-based genetic analyses, such as the Genome-Wide Association Study (GWAS), have proven powerful for describing the genetic complexities of common disease in epidemiologic research. However, the significant challenges faced by population-based study designs have resulted in revitalization of family-based approaches, including twin studies. Twin studies are unique in their ability to ascertain both heritable and environmental contributions to human disease. Several regional and national twin registries have been constructed using a variety of methods to identify potential twins. A significant challenge in constructing these large twin registries includes the substantial resources required to recruit participants, collect phenotypic data, and update the registries as time progresses. Here we describe the use of the Marshfield Clinic electronic medical record (EMR) to identify a cohort of 19,226 patients enriched for twins or multiples. This cohort defines the Marshfield Clinic Twin/Multiple Birth Cohort (MCTC). An EMR system provides both a mechanism to identify potential twins and a source of detailed phenotypic data in near real time without the need for patient contact outside standard medical care. To demonstrate that the MCTC can be used for genetic-based epidemiologic research, concordance rates for muscular dystrophy (MD) and fragile-X syndrome-two highly heritable diseases-were assessed. Observations indicate that both MD and fragile-X syndrome are highly correlated among affected twins in the MCTC (P ≅ 3.7 × 10(-6) and 1.1 × 10(-4) , respectively). These findings suggest that EMR systems may not only be an effective resource for predicting families of twins, but can also be rapidly applied to epidemiologic research.

Cone structure in subjects with known genetic relative risk for AMD.

PURPOSE: Utilize high-resolution imaging to examine retinal anatomy in patients with known genetic relative risk (RR) for developing age-related macular degeneration (AMD). METHODS: Forty asymptomatic subjects were recruited (9 men, 31 women; age range, 51 to 69 years; mean age, 61.4 years). Comprehensive eye examination, fundus photography, and high-resolution retinal imaging using spectral domain optical coherence tomography and adaptive optics were performed on each patient. Genetic RR scores were developed using an age-independent algorithm. Adaptive optics scanning light ophthalmoscope images were acquired in the macula extending to 10 degrees temporal and superior from fixation and were used to calculate cone density in up to 35 locations for each subject. RESULTS: Relative risk was not significantly predictive of fundus grade (p = 0.98). Only patients with a high RR displayed drusen on Cirrus or Bioptigen OCT. Compared to an eye with a grade of 0, an eye with a fundus grade equal to or greater than 1 had a 12% decrease in density (p < 0.0001) and a 5% increase in spacing (p = 0.0014). No association between genetic RR and either cone density (p = 0.435) or spacing (p = 0.538) was found. Three distinct adaptive optics scanning light ophthalmoscope phenotypical variations of photoreceptor appearance were noted in patients with grade 1 to 3 fundi. These included variable reflectivity of photoreceptors, decreased waveguiding, and altered photoreceptor mosaic overlying drusen. CONCLUSIONS: Our data demonstrate the potential of multimodal assessment in the understanding of early anatomical changes associated with AMD. Adaptive optics scanning light ophthalmoscope imaging reveals a decrease in photoreceptor density and increased spacing in patients with grade 1 to 3 fundi, as well as a spectrum of photoreceptor changes, ranging from variability in reflectivity to decreased density. Future longitudinal studies are needed in genetically characterized subjects to assess the significance of these findings with respect to the development and progression of AMD.

Longitudinal dynamics of farmer and livestock nasal and faecal microbiomes and resistomes.

Globally, half a billion people are employed in animal agriculture and are directly exposed to the associated microorganisms. However, the extent to which such exposures affect resident human microbiomes is unclear. Here we conducted a longitudinal profiling of the nasal and faecal microbiomes of 66 dairy farmers and 166 dairy cows over a year-long period. We compare farmer microbiomes to those of 60 age-, sex- and ZIP code-matched people with no occupational exposures to farm animals (non-farmers). We show that farming is associated with microbiomes containing livestock-associated microbes; this is most apparent in the nasal bacterial community, with farmers harbouring a richer and more diverse nasal community than non-farmers. Similarly, in the gut microbial communities, we identify more shared microbial lineages between cows and farmers from the same farms. Additionally, we find that shared microbes are associated with antibiotic resistance genes. Overall, our study demonstrates the interconnectedness of human and animal microbiomes.

Machine Learning Prediction of Treatment Response to Inhaled Corticosteroids in Asthma.

BACKGROUND: Although inhaled corticosteroids (ICS) are the first-line therapy for patients with persistent asthma, many patients continue to have exacerbations. We developed machine learning models to predict the ICS response in patients with asthma. METHODS: The subjects included asthma patients of European ancestry (n = 1371; 448 children; 916 adults). A genome-wide association study was performed to identify the SNPs associated with ICS response. Using the SNPs identified, two machine learning models were developed to predict ICS response: (1) least absolute shrinkage and selection operator (LASSO) regression and (2) random forest. RESULTS: The LASSO regression model achieved an AUC of 0.71 (95% CI 0.67-0.76; sensitivity: 0.57; specificity: 0.75) in an independent test cohort, and the random forest model achieved an AUC of 0.74 (95% CI 0.70-0.78; sensitivity: 0.70; specificity: 0.68). The genes contributing to the prediction of ICS response included those associated with ICS responses in asthma (TPSAB1, FBXL16), asthma symptoms and severity (ABCA7, CNN2, PTRN3, and BSG/CD147), airway remodeling (ELANE, FSTL3), mucin production (GAL3ST), leukotriene synthesis (GPX4), allergic asthma (ZFPM1, SBNO2), and others. CONCLUSIONS: An accurate risk prediction of ICS response can be obtained using machine learning methods, with the potential to inform personalized treatment decisions. Further studies are needed to examine if the integration of richer phenotype data could improve risk prediction.

Estimating the efficacy of pharmacogenomics over a lifetime.

It is well known that common variants in specific genes influence drug metabolism and response, but it is currently unknown what fraction of patients are given prescriptions over a lifetime that could be contraindicated by their pharmacogenomic profiles. To determine the clinical utility of pharmacogenomics over a lifetime in a general patient population, we sequenced the genomes of 300 deceased Marshfield Clinic patients linked to lifelong medical records. Genetic variants in 33 pharmacogenes were evaluated for their lifetime impact on drug prescribing using extensive electronic health records. Results show that 93% of the 300 deceased patients carried clinically relevant variants. Nearly 80% were prescribed approximately three medications on average that may have been impacted by these variants. Longitudinal data suggested that the optimal age for pharmacogenomic testing was prior to age 50, but the optimal age is greatly influenced by the stability of the population in the healthcare system. This study emphasizes the broad clinical impact of pharmacogenomic testing over a lifetime and demonstrates the potential application of genomic medicine in a general patient population for the advancement of precision medicine.

Genetic risk score in multiple sclerosis is associated with unique gut microbiome.

Multiple sclerosis (MS) is a complex autoimmune disease in which both the roles of genetic susceptibility and environmental/microbial factors have been investigated. More than 200 genetic susceptibility variants have been identified along with the dysbiosis of gut microbiota, both independently have been shown to be associated with MS. We hypothesize that MS patients harboring genetic susceptibility variants along with gut microbiome dysbiosis are at a greater risk of exhibiting the disease. We investigated the genetic risk score for MS in conjunction with gut microbiota in the same cohort of 117 relapsing remitting MS (RRMS) and 26 healthy controls. DNA samples were genotyped using Illumina's Infinium Immuno array-24 v2 chip followed by calculating genetic risk score and the microbiota was determined by sequencing the V4 hypervariable region of the 16S rRNA gene. We identified two clusters of MS patients, Cluster A and B, both having a higher genetic risk score than the control group. However, the MS cases in cluster B not only had a higher genetic risk score but also showed a distinct gut microbiome than that of cluster A. Interestingly, cluster A which included both healthy control and MS cases had similar gut microbiome composition. This could be due to (i) the non-active state of the disease in that group of MS patients at the time of fecal sample collection and/or (ii) the restoration of the gut microbiome post disease modifying therapy to treat the MS. Our study showed that there seems to be an association between genetic risk score and gut microbiome dysbiosis in triggering the disease in a small cohort of MS patients. The MS Cluster A who have a higher genetic risk score but microbiome profile similar to that of healthy controls could be due to the remitting phase of the disease or due to the effect of disease modifying therapies.

Dietary intake in the Personalized Medicine Research Project: a resource for studies of gene-diet interaction.

BACKGROUND: To describe the dietary intake of participants in the Personalized Medicine Research Project (PMRP), and to quantify differences in nutrient intake by smoking status and APOE4-a genetic marker that has been shown to modify the association between risk factors and outcomes. METHODS: The PMRP is a population-based DNA, plasma and serum biobank of more than 20,000 adults aged 18 years and older in central Wisconsin. A questionnaire at enrollment captures demographic information as well as self-reported smoking and alcohol intake. The protocol was amended to include the collection of dietary intake and physical activity via self-reported questionnaires: the National Cancer Institute 124-item Diet History Questionnaire and the Baecke Physical Activity Questionnaire. These questionnaires were mailed out to previously enrolled participants. APOE was genotyped in all subjects. RESULTS: The response rate to the mailed questionnaires was 68.2% for subjects who could still be contacted (alive with known address). Participants ranged in age from 18 to 98 years (mean 54.7) and 61% were female. Dietary intake is variable when comparing gender, age, smoking, and APOE4. Over 50% of females are dietary supplement users; females have higher supplement intake than males, but both have increasing supplement use as age increases. Food energy, total fat, cholesterol, protein, and alcohol intake decreases as both males and females age. Female smokers had higher macronutrient intake, whereas male nonsmokers had higher macronutrient intake. Nonsmokers in both genders use more supplements. In females, nonsmokers and smokers with APOE4 had higher supplement use. In males, nonsmokers with APOE4 had higher supplement use between ages 18-39 only, and lower supplement use at ages above 39. Male smokers with APOE4 had lower supplement use. CONCLUSION: Dietary intake in PMRP subjects is relatively consistent with data from the National Health and Nutrition Examination Survey (NHANES). Findings suggest a possible correlation between the use of supplements and APOE4. The PMRP dietary data can benefit studies of gene-environment interactions and the development of common diseases.

Development of an Integrated Platform Using Multidisciplinary Real-World Data to Facilitate Biomarker Discovery for Medical Products.

Translational multidisciplinary research is important for the Center for Devices and Radiological Health's efforts for utilizing real-world data (RWD) to enhance predictive evaluation of medical device performance in patient subpopulations. As part of our efforts for developing new RWD-based evidentiary approaches, including in silico discovery of device-related risk predictors and biomarkers, this study aims to characterize the sex/race-related trends in hip replacement outcomes and identify corresponding candidate single nucleotide polymorphisms (SNPs). Adverse outcomes were assessed by deriving RWD from a retrospective analysis of hip replacement hospital discharge data from the National Inpatient Sample (NIS). Candidate SNPs were explored using pre-existing data from the Personalized Medicine Research Project (PMRP). High-Performance Integrated Virtual Environment was used for analyzing and visualizing putative associations between SNPs and adverse outcomes. Ingenuity Pathway Analysis (IPA) was used for exploring plausibility of the sex-related candidate SNPs and characterizing gene networks associated with the variants of interest. The NIS-based epidemiologic evidence showed that periprosthetic osteolysis (PO) was most prevalent among white men. The PMRP-based genetic evidence associated the PO-related male predominance with rs7121 (odds ratio = 4.89; 95% confidence interval = 1.41-17.05) and other candidate SNPs. SNP-based IPA analysis of the expected gene expression alterations and corresponding signaling pathways suggested possible role of sex-related metabolic factors in development of PO, which was substantiated by ad hoc epidemiologic analysis identifying the sex-related differences in metabolic comorbidities in men vs. women with hip replacement-related PO. Thus, our in silico study illustrates RWD-based evidentiary approaches that may facilitate cost/time-efficient discovery of biomarkers for informing use of medical products.

Use of an electronic medical record to characterize cases of intermediate statin-induced muscle toxicity.

Statin use can be accompanied by a variety of musculoskeletal complaints. The authors describe the clinical characteristics of case patients who experienced adverse statin-induced musculoskeletal symptoms within a large population-based cohort in Central Wisconsin. Case status was determined based on elevated serum creatine kinase (CK) levels and the presence of at least 1 physician note reflecting an increased index of suspicion for statin intolerance. From the medical records of nearly 2 million unique patients, the authors identified more than 20,000 potential study patients ( approximately 1%) having CK data and at least 1 exposure to a statin drug. Manual screening was completed on 2227 patients with CK levels in the upper 10th percentile. Of those screened, 267 met inclusion criteria (12.0% eligibility) and 218 agreed to participate in a retrospective study characterizing the risk determinants of statin-induced muscle toxicity. Three categoric pain variables were graded retrospectively (distribution, location, and severity of pain). The presenting complaints of the case patients were extremely heterogeneous. The number of patients with a compelling pain syndrome (diffuse, proximal muscle pain of high intensity) increased at higher serum CK levels; the number of patients with indeterminate pain variables decreased at higher serum CK levels. The lines reflecting these relationships cross at a CK level of approximately 1175 U/L, approximately half the threshold level needed to make a clinical diagnosis of "myopathy" (ie, CK >10-fold the upper limit of normal).

Enrichment sampling for a multi-site patient survey using electronic health records and census data.

Objective: We describe a stratified sampling design that combines electronic health records (EHRs) and United States Census (USC) data to construct the sampling frame and an algorithm to enrich the sample with individuals belonging to rarer strata. Materials and Methods: This design was developed for a multi-site survey that sought to examine patient concerns about and barriers to participating in research studies, especially among under-studied populations (eg, minorities, low educational attainment). We defined sampling strata by cross-tabulating several socio-demographic variables obtained from EHR and augmented with census-block-level USC data. We oversampled rarer and historically underrepresented subpopulations. Results: The sampling strategy, which included USC-supplemented EHR data, led to a far more diverse sample than would have been expected under random sampling (eg, 3-, 8-, 7-, and 12-fold increase in African Americans, Asians, Hispanics and those with less than a high school degree, respectively). We observed that our EHR data tended to misclassify minority races more often than majority races, and that non-majority races, Latino ethnicity, younger adult age, lower education, and urban/suburban living were each associated with lower response rates to the mailed surveys. Discussion: We observed substantial enrichment from rarer subpopulations. The magnitude of the enrichment depends on the accuracy of the variables that define the sampling strata and the overall response rate. Conclusion: EHR and USC data may be used to define sampling strata that in turn may be used to enrich the final study sample. This design may be of particular interest for studies of rarer and understudied populations.

Healthcare provider education to support integration of pharmacogenomics in practice: the eMERGE Network experience.

Ten organizations within the Electronic Medical Records and Genomics Network developed programs to implement pharmacogenomic sequencing and clinical decision support into clinical settings. Recognizing the importance of informed prescribers, a variety of strategies were used to incorporate provider education to support implementation. Education experiences with pharmacogenomics are described within the context of each organization's prior involvement, including the scope and scale of implementation specific to their Electronic Medical Records and Genomics projects. We describe common and distinct education strategies, provide exemplars and share challenges. Lessons learned inform future perspectives. Future pharmacogenomics clinical implementation initiatives need to include funding toward implementing provider education and evaluating outcomes.