A Platform for Generation of Chamber-Specific Cardiac Tissues and Disease Modeling.

Tissue engineering using cardiomyocytes derived from human pluripotent stem cells holds a promise to revolutionize drug discovery, but only if limitations related to cardiac chamber specification and platform versatility can be overcome. We describe here a scalable tissue-cultivation platform that is cell source agnostic and enables drug testing under electrical pacing. The plastic platform enabled on-line noninvasive recording of passive tension, active force, contractile dynamics, and Ca2+ transients, as well as endpoint assessments of action potentials and conduction velocity. By combining directed cell differentiation with electrical field conditioning, we engineered electrophysiologically distinct atrial and ventricular tissues with chamber-specific drug responses and gene expression. We report, for the first time, engineering of heteropolar cardiac tissues containing distinct atrial and ventricular ends, and we demonstrate their spatially confined responses to serotonin and ranolazine. Uniquely, electrical conditioning for up to 8 months enabled modeling of polygenic left ventricular hypertrophy starting from patient cells.

Thymic regulatory T cells arise via two distinct developmental programs.

The developmental programs that generate a broad repertoire of regulatory T cells (Treg cells) able to respond to both self antigens and non-self antigens remain unclear. Here we found that mature Treg cells were generated through two distinct developmental programs involving CD25+ Treg cell progenitors (CD25+ TregP cells) and Foxp3lo Treg cell progenitors (Foxp3lo TregP cells). CD25+ TregP cells showed higher rates of apoptosis and interacted with thymic self antigens with higher affinity than did Foxp3lo TregP cells, and had a T cell antigen receptor repertoire and transcriptome distinct from that of Foxp3lo TregP cells. The development of both CD25+ TregP cells and Foxp3lo TregP cells was controlled by distinct signaling pathways and enhancers. Transcriptomics and histocytometric data suggested that CD25+ TregP cells and Foxp3lo TregP cells arose by coopting negative-selection programs and positive-selection programs, respectively. Treg cells derived from CD25+ TregP cells, but not those derived from Foxp3lo TregP cells, prevented experimental autoimmune encephalitis. Our findings indicate that Treg cells arise through two distinct developmental programs that are both required for a comprehensive Treg cell repertoire capable of establishing immunotolerance.

SARS-CoV-2 pathogenesis in an angiotensin II-induced heart-on-a-chip disease model and extracellular vesicle screening.

Adverse cardiac outcomes in COVID-19 patients, particularly those with preexisting cardiac disease, motivate the development of human cell-based organ-on-a-chip models to recapitulate cardiac injury and dysfunction and for screening of cardioprotective therapeutics. Here, we developed a heart-on-a-chip model to study the pathogenesis of SARS-CoV-2 in healthy myocardium established from human induced pluripotent stem cell (iPSC)-derived cardiomyocytes and a cardiac dysfunction model, mimicking aspects of preexisting hypertensive disease induced by angiotensin II (Ang II). We recapitulated cytopathic features of SARS-CoV-2-induced cardiac damage, including progressively impaired contractile function and calcium handling, apoptosis, and sarcomere disarray. SARS-CoV-2 presence in Ang II-treated hearts-on-a-chip decreased contractile force with earlier onset of contractile dysfunction and profoundly enhanced inflammatory cytokines compared to SARS-CoV-2 alone. Toward the development of potential therapeutics, we evaluated the cardioprotective effects of extracellular vesicles (EVs) from human iPSC which alleviated the impairment of contractile force, decreased apoptosis, reduced the disruption of sarcomeric proteins, and enhanced beta-oxidation gene expression. Viral load was not affected by either Ang II or EV treatment. We identified MicroRNAs miR-20a-5p and miR-19a-3p as potential mediators of cardioprotective effects of these EVs.

Discovery and fine-mapping of height loci via high-density imputation of GWASs in individuals of African ancestry.

Although many loci have been associated with height in European ancestry populations, very few have been identified in African ancestry individuals. Furthermore, many of the known loci have yet to be generalized to and fine-mapped within a large-scale African ancestry sample. We performed sex-combined and sex-stratified meta-analyses in up to 52,764 individuals with height and genome-wide genotyping data from the African Ancestry Anthropometry Genetics Consortium (AAAGC). We additionally combined our African ancestry meta-analysis results with published European genome-wide association study (GWAS) data. In the African ancestry analyses, we identified three novel loci (SLC4A3, NCOA2, ECD/FAM149B1) in sex-combined results and two loci (CRB1, KLF6) in women only. In the African plus European sex-combined GWAS, we identified an additional three novel loci (RCCD1, G6PC3, CEP95) which were equally driven by AAAGC and European results. Among 39 genome-wide significant signals at known loci, conditioning index SNPs from European studies identified 20 secondary signals. Two of the 20 new secondary signals and none of the 8 novel loci had minor allele frequencies (MAF) < 5%. Of 802 known European height signals, 643 displayed directionally consistent associations with height, of which 205 were nominally significant (p < 0.05) in the African ancestry sex-combined sample. Furthermore, 148 of 241 loci contained ≤20 variants in the credible sets that jointly account for 99% of the posterior probability of driving the associations. In summary, trans-ethnic meta-analyses revealed novel signals and further improved fine-mapping of putative causal variants in loci shared between African and European ancestry populations.

Incorporating G6PD genotyping to identify patients with G6PD deficiency.

Glucose-6-phosphate-dehydrogenase (G6PD) deficiency is a common X-linked enzyme disorder associated with hemolytic anemia after exposure to fava beans or certain medications. Activity testing is the gold standard for detecting G6PD deficiency; however, this test is affected by various hematologic parameters. Clinical G6PD genotyping is now included in pharmacogenetic arrays and clinical sequencing efforts and may be reconciled with activity results. Patients (n = 1391) enrolled on an institutional pharmacogenetic testing protocol underwent clinical G6PD genotyping for 164 G6PD variants. An algorithm accounting for known interferences with the activity assay is proposed. We developed clinical decision support alerts to inform prescribers when high-risk medications were prescribed, warning of gene-drug interactions and recommending therapy alteration. Of 1391 patients with genotype results, 1334 (95.9%) patients were predicted to have normal G6PD activity, 30 (2.1%) were predicted to have variable G6PD activity and 27 (2%) were predicted to have deficient G6PD activity. Of the 417 patients with a normal genotype and an activity result, 415 (99.5%) had a concordant normal G6PD phenotype. Of the 21 patients with a deficient genotype and an activity result, 18 (85.7%) had a concordant deficient activity result. Genotyping reassigned phenotype in five patients with discordant genotype and activity results: three switched from normal to deficient, and two switched from deficient to normal. G6PD activity and genotyping are two independent testing methods that can be used in conjunction to assign a more informed G6PD phenotype than either method alone.

Associations of Mitochondrial and Nuclear Mitochondrial Variants and Genes with Seven Metabolic Traits.

Mitochondria (MT), the major site of cellular energy production, are under dual genetic control by 37 mitochondrial DNA (mtDNA) genes and numerous nuclear genes (MT-nDNA). In the CHARGEmtDNA+ Consortium, we studied genetic associations of mtDNA and MT-nDNA associations with body mass index (BMI), waist-hip-ratio (WHR), glucose, insulin, HOMA-B, HOMA-IR, and HbA1c. This 45-cohort collaboration comprised 70,775 (insulin) to 170,202 (BMI) pan-ancestry individuals. Validation and imputation of mtDNA variants was followed by single-variant and gene-based association testing. We report two significant common variants, one in MT-ATP6 associated (p ≤ 5E-04) with WHR and one in the D-loop with glucose. Five rare variants in MT-ATP6, MT-ND5, and MT-ND6 associated with BMI, WHR, or insulin. Gene-based meta-analysis identified MT-ND3 associated with BMI (p ≤ 1E-03). We considered 2,282 MT-nDNA candidate gene associations compiled from online summary results for our traits (20 unique studies with 31 dataset consortia's genome-wide associations [GWASs]). Of these, 109 genes associated (p ≤ 1E-06) with at least 1 of our 7 traits. We assessed regulatory features of variants in the 109 genes, cis- and trans-gene expression regulation, and performed enrichment and protein-protein interactions analyses. Of the identified mtDNA and MT-nDNA genes, 79 associated with adipose measures, 49 with glucose/insulin, 13 with risk for type 2 diabetes, and 18 with cardiovascular disease, indicating for pleiotropic effects with health implications. Additionally, 21 genes related to cholesterol, suggesting additional important roles for the genes identified. Our results suggest that mtDNA and MT-nDNA genes and variants reported make important contributions to glucose and insulin metabolism, adipocyte regulation, diabetes, and cardiovascular disease.

Donor-specific phenotypic variation in hiPSC cardiomyocyte-derived exosomes impacts endothelial cell function.

Exosomes are an important mechanism of cell-cell interaction in the cardiovascular system, both in maintaining homeostasis and in stress response. Interindividual differences that alter content in exosomes may play a role in cardiovascular disease pathology. To study the effect of interindividual cardiomyocyte (CM) variation, we characterized exosomal content in phenotypically diverse human induced pluripotent stem cell-derived CMs (hiPSC-CMs). Cell lines were generated from six participants in the HyperGEN cohort: three with left ventricular hypertrophy (LVH) and three with normal left ventricular mass (LVM). Sequence analysis of the intracellular and exosomal RNA populations showed distinct expression pattern differences between hiPSC-CM lines derived from individuals with LVH and those with normal LVM. Functional analysis of hiPSC-endothelial cells (hiPSC-ECs) treated with exosomes from both hiPSC-CM groups showed significant variation in response, including differences in tube formation, migration, and proliferation. Overall, treatment of hiPSC-ECs with exosomes resulted in significant expression changes associated with angiogenesis and endothelial cell vasculogenesis. However, the hiPSC-ECs treated with exosomes from the LVH-affected donors exhibited significantly increased proliferation but decreased tube formation and migration, suggesting angiogenic dysregulation.NEW & NOTEWORTHY The intracellular RNA and the miRNA content in exosomes are significantly different in hiPSC-CMs derived from LVH-affected individuals compared with those from unaffected individuals. Treatment of endothelial cells with these exosomes functionally affects cellular phenotypes in a donor-specific manner. These findings provide novel insight into underlying mechanisms of hypertrophic cell signaling between different cell types. With a growing interest in stem cells and exosomes for cardiovascular therapeutic use, this also provides information important for regenerative medicine.

Human IL12RB1 expression is allele-biased and produces a novel IL12 response regulator.

Human IL12RB1 is an autosomal gene that is essential for mycobacterial disease resistance and T cell differentiation. Using primary human tissue and PBMCs, we demonstrate that lung and T cell IL12RB1 expression is allele-biased, and the extent to which cells express one IL12RB1 allele is unaffected by activation. Furthermore following its expression the IL12RB1 pre-mRNA is processed into either IL12RB1 Isoform 1 (IL12Rß1, a positive regulator of IL12 responsiveness) or IL12RB1 Isoform 2 (a protein of heretofore unknown function). T cells choice to process pre-mRNA into Isoform 1 or Isoform 2 is controlled by intragenic competition of IL12RB1 exon 9-10 splicing with IL12RB1 exon 9b splicing, as well as an IL12RB1 exon 9b-associated polyadenylation site. Heterogeneous nuclear ribonucleoprotein H (hnRNP H) binds near the regulated polyadenylation site, but is not required for exon 9b polyadenylation. Finally, microRNA-mediated knockdown experiments demonstrated that IL12RB1 Isoform 2 promotes T cell IL12 responses. Collectively, our data support a model wherein tissue expression of human IL12RB1 is allele-biased and produces an hnRNP H-bound pre-mRNA, the processing of which generates a novel IL12 response regulator.

Concordance between glucose-6-phosphate dehydrogenase (G6PD) genotype and phenotype and rasburicase use in patients with hematologic malignancies.

Phenotypic rather than genotypic tests remain the gold standard for diagnosing glucose-6-phosphate dehydrogenase (G6PD) deficiency. However, with increasing use of genomic arrays and whole exome or genome sequencing, G6PD genetic data are increasingly available. We examined the utility of G6PD genetic data in patients with hematologic malignancies and the association of G6PD genotype and phenotype with rasburicase-induced methemoglobinemia. We analyzed G6PD activity for 990 patients. Genotype data were available from the Affymetrix DMET array (n = 379), whole exome sequencing (n = 374), and/or the Illumina exome array (n = 634) for 645 patients. Medical records of 341 patients with methemoglobin measures were assessed for the administration of rasburicase. We observed 5 non-synonymous SNPs, 4 of which were known to be associated with deficient G6PD activity (WHO Class I-III). Genotyping 367 males resulted in a positive predictive value of 81.8% (47.8-96.8%), and two males with a Class I-III allele having normal activity both received a red blood cell transfusion prior to the activity assay. However, genotyping males had only 39.1% (20.5-61.2%) sensitivity. Two of the 12 heterozygous females had deficient G6PD activity. Rasburicase-induced methemoglobinemia occurred in 6 patients, 5 of whom had at least one Class I-III allele, despite 2 of these having normal G6PD activity. We conclude that although an apparent nondeficient genotype does not necessarily imply a normal phenotype, a deficient genotype result indicates a deficient phenotype in those without transfusions, and may be a useful adjuct to phenotype to prevent adverse drug reactions.

Preemptive clinical pharmacogenetics implementation: current programs in five US medical centers.

Although the field of pharmacogenetics has existed for decades, practioners have been slow to implement pharmacogenetic testing in clinical care. Numerous publications describe the barriers to clinical implementation of pharmacogenetics. Recently, several freely available resources have been developed to help address these barriers. In this review, we discuss current programs that use preemptive genotyping to optimize the pharmacotherapy of patients. Array-based preemptive testing includes a large number of relevant pharmacogenes that impact multiple high-risk drugs. Using a preemptive approach allows genotyping results to be available prior to any prescribing decision so that genomic variation may be considered as an inherent patient characteristic in the planning of therapy. This review describes the common elements among programs that have implemented preemptive genotyping and highlights key processes for implementation, including clinical decision support.

Hepatic Abundance and Activity of Androgen- and Drug-Metabolizing Enzyme UGT2B17 Are Associated with Genotype, Age, and Sex.

The major objective of this study was to investigate the association of genetic and nongenetic factors with variability in protein abundance and in vitro activity of the androgen-metabolizing enzyme UGT2B17 in human liver microsomes (n = 455). UGT2B17 abundance was quantified by liquid chromatography-tandem mass spectrometry proteomics, and enzyme activity was determined by using testosterone and dihydrotestosterone as in vitro probe substrates. Genotyping or gene resequencing and mRNA expression were also evaluated. Multivariate analysis was used to test the association of UGT2B17 copy number variation, single nucleotide polymorphisms (SNPs), age, and sex with its mRNA expression, abundance, and activity. UGT2B17 gene copy number and SNPs (rs7436962, rs9996186, rs28374627, and rs4860305) were associated with gene expression, protein levels, and androgen glucuronidation rates in a gene dose-dependent manner. UGT2B17 protein (mean ± S.D. picomoles per milligram of microsomal protein) is sparsely expressed in children younger than 9 years (0.12 ± 0.24 years) but profoundly increases from age 9 years to adults (∼10-fold) with ∼2.6-fold greater abundance in males than in females (1.2 vs. 0.47). Association of androgen glucuronidation with UGT2B15 abundance was observed only in the low UGT2B17 expressers. These data can be used to predict variability in the metabolism of UGT2B17 substrates. Drug companies should include UGT2B17 in early phenotyping assays during drug discovery to avoid late clinical failures.

The ENTPD1 promoter polymorphism -860 A > G (rs3814159) is associated with increased gene transcription, protein expression, CD39/NTPDase1 enzymatic activity, and thromboembolism risk.

Ectonucleoside triphosphate diphosphohydrolase 1 (NTPDase1) degrades the purines ATP and ADP that are key regulators of inflammation and clotting. We hypothesized that NTPDase1 polymorphisms exist and that they regulate this pathway. We sequenced the ENTPD1 gene (encoding NTPDase1) in 216 subjects then assessed genotypes in 2 cohorts comprising 2213 humans to identify ENTPD1 polymorphisms associated with venous thromboembolism (VTE). The G allele of the intron 1 polymorphism rs3176891 was more common in VTE vs. controls (odds ratio 1.26-1.9); it did not affect RNA splicing, but it was in strong linkage disequilibrium with the G allele of the promoter polymorphism rs3814159, which increased transcriptional activity by 8-fold. Oligonucleotides containing the G allele of this promoter region bound nuclear extracts more avidly. Carriers of rs3176891 G had endothelial cells with increased NTPDase1 activity and protein expression, and had platelets with enhanced aggregation. Thus, the G allele of rs3176891 marks a haplotype associated with increased clotting and platelet aggregation attributable to a promoter variant associated with increased transcription, expression, and activity of NTPDase1. We term this gain-of-function phenotype observed with rs3814159 G "CD39 Denver."-Maloney, J. P., Branchford, B. R., Brodsky, G. L., Cosmic, M. S., Calabrese, D. W., Aquilante, C. L., Maloney, K. W., Gonzalez, J. R., Zhang, W., Moreau, K. L., Wiggins, K. L., Smith, N. L., Broeckel, U., Di Paola, J. The ENTPD1 promoter polymorphism -860 A > G (rs3814159) is associated with increased gene transcription, protein expression, CD39/NTPDase1 enzymatic activity, and thromboembolism risk.

The introduction of RNA-DNA differences underlies interindividual variation in the human IL12RB1 mRNA repertoire.

Human interleukin 12 and interleukin 23 (IL12/23) influence susceptibility or resistance to multiple diseases. However, the reasons underlying individual differences in IL12/23 sensitivity remain poorly understood. Here we report that in human peripheral blood mononuclear cells (PBMCs) and inflamed lungs, the majority of interleukin-12 receptor ß1 (IL12RB1) mRNAs contain a number of RNA-DNA differences (RDDs) that concentrate in sequences essential to IL12Rß1's binding of IL12p40, the protein subunit common to both IL-12 and IL-23. IL12RB1 RDDs comprise multiple RDD types and are detectable by next-generation sequencing and classic Sanger sequencing. As a consequence of these RDDs, the resulting IL12Rß1 proteins have an altered amino acid sequence that could not be predicted on the basis of genomic DNA sequencing alone. Importantly, the introduction of RDDs into IL12RB1 mRNAs negatively regulates IL12Rß1's binding of IL12p40 and is sensitive to activation. Collectively, these results suggest that the introduction of RDDs into an individual's IL12RB1 mRNA repertoire is a novel determinant of IL12/23 sensitivity.

Incorporating a Genetic Risk Score Into Coronary Heart Disease Risk Estimates: Effect on Low-Density Lipoprotein Cholesterol Levels (the MI-GENES Clinical Trial).

BACKGROUND: Whether knowledge of genetic risk for coronary heart disease (CHD) affects health-related outcomes is unknown. We investigated whether incorporating a genetic risk score (GRS) in CHD risk estimates lowers low-density lipoprotein cholesterol (LDL-C) levels. METHODS AND RESULTS: Participants (n=203, 45-65 years of age, at intermediate risk for CHD, and not on statins) were randomly assigned to receive their 10-year probability of CHD based either on a conventional risk score (CRS) or CRS + GRS ((+)GRS). Participants in the (+)GRS group were stratified as having high or average/low GRS. Risk was disclosed by a genetic counselor followed by shared decision making regarding statin therapy with a physician. We compared the primary end point of LDL-C levels at 6 months and assessed whether any differences were attributable to changes in dietary fat intake, physical activity levels, or statin use. Participants (mean age, 59.4±5 years; 48% men; mean 10-year CHD risk, 8.5±4.1%) were allocated to receive either CRS (n=100) or (+)GRS (n=103). At the end of the study period, the (+)GRS group had a lower LDL-C than the CRS group (96.5±32.7 versus 105.9±33.3 mg/dL; P=0.04). Participants with high GRS had lower LDL-C levels (92.3±32.9 mg/dL) than CRS participants (P=0.02) but not participants with low GRS (100.9±32.2 mg/dL; P=0.18). Statins were initiated more often in the (+)GRS group than in the CRS group (39% versus 22%, P<0.01). No significant differences in dietary fat intake and physical activity levels were noted. CONCLUSIONS: Disclosure of CHD risk estimates that incorporated genetic risk information led to lower LDL-C levels than disclosure of CHD risk based on conventional risk factors alone. CLINICAL TRIAL REGISTRATION: URL: http://www.clinicaltrials.gov. Unique identifier: NCT01936675.

Genetic and Nongenetic Factors Associated with Protein Abundance of Flavin-Containing Monooxygenase 3 in Human Liver.

Hepatic flavin-containing mono-oxygenase 3 (FMO3) metabolizes a broad array of nucleophilic heteroatom (e.g., N or S)-containing xenobiotics (e.g., amphetamine, sulindac, benzydamine, ranitidine, tamoxifen, nicotine, and ethionamide), as well as endogenous compounds (e.g., catecholamine and trimethylamine). To predict the effect of genetic and nongenetic factors on the hepatic metabolism of FMO3 substrates, we quantified FMO3 protein abundance in human liver microsomes (HLMs; n = 445) by liquid chromatography-tandem mass chromatography proteomics. Genotyping/gene resequencing, mRNA expression, and functional activity (with benzydamine as probe substrate) of FMO3 were also evaluated. FMO3 abundance increased 2.2-fold (13.0 ± 11.4 pmol/mg protein vs. 28.0 ± 11.8 pmol/mg protein) from neonates to adults. After 6 years of age, no significant difference in FMO3 abundance was found between children and adults. Female donors exhibited modestly higher mRNA fragments per kilobase per million reads values (139.9 ± 76.9 vs. 105.1 ± 73.1; P < 0.001) and protein FMO3 abundance (26.7 ± 12.0 pmol/mg protein vs. 24.1 ± 12.1 pmol/mg protein; P < 0.05) compared with males. Six single nucleotide polymorphisms (SNPs), including rs2064074, rs28363536, rs2266782 (E158K), rs909530 (N285N), rs2266780 (E308G), and rs909531, were associated with significantly decreased protein abundance. FMO3 abundance in individuals homozygous and heterozygous for haplotype 3 (H3), representing variant alleles for all these SNPs (except rs2066534), were 50.8% (P < 0.001) and 79.5% (P < 0.01), respectively, of those with the reference homozygous haplotype (H1, representing wild-type). In summary, FMO3 protein abundance is significantly associated with age, gender, and genotype. These data are important in predicting FMO3-mediated heteroatom-oxidation of xenobiotics and endogenous biomolecules in the human liver.

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.

Comprehensive evaluation of AmpliSeq transcriptome, a novel targeted whole transcriptome RNA sequencing methodology for global gene expression analysis.

BACKGROUND: Whole transcriptome sequencing (RNA-seq) represents a powerful approach for whole transcriptome gene expression analysis. However, RNA-seq carries a few limitations, e.g., the requirement of a significant amount of input RNA and complications led by non-specific mapping of short reads. The Ion AmpliSeq Transcriptome Human Gene Expression Kit (AmpliSeq) was recently introduced by Life Technologies as a whole-transcriptome, targeted gene quantification kit to overcome these limitations of RNA-seq. To assess the performance of this new methodology, we performed a comprehensive comparison of AmpliSeq with RNA-seq using two well-established next-generation sequencing platforms (Illumina HiSeq and Ion Torrent Proton). We analyzed standard reference RNA samples and RNA samples obtained from human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs). RESULTS: Using published data from two standard RNA reference samples, we observed a strong concordance of log2 fold change for all genes when comparing AmpliSeq to Illumina HiSeq (Pearson's r = 0.92) and Ion Torrent Proton (Pearson's r = 0.92). We used ROC, Matthew's correlation coefficient and RMSD to determine the overall performance characteristics. All three statistical methods demonstrate AmpliSeq as a highly accurate method for differential gene expression analysis. Additionally, for genes with high abundance, AmpliSeq outperforms the two RNA-seq methods. When analyzing four closely related hiPSC-CM lines, we show that both AmpliSeq and RNA-seq capture similar global gene expression patterns consistent with known sources of variations. CONCLUSIONS: Our study indicates that AmpliSeq excels in the limiting areas of RNA-seq for gene expression quantification analysis. Thus, AmpliSeq stands as a very sensitive and cost-effective approach for very large scale gene expression analysis and mRNA marker screening with high accuracy.

Whole exome sequencing identifies a POLRID mutation segregating in a father and two daughters with findings of Klippel-Feil and Treacher Collins syndromes.

We report on a father and his two daughters diagnosed with Klippel-Feil syndrome (KFS) but with craniofacial differences (zygomatic and mandibular hypoplasia and cleft palate) and external ear abnormalities suggestive of Treacher Collins syndrome (TCS). The diagnosis of KFS was favored, given that the neck anomalies were the predominant manifestations, and that the diagnosis predated later recognition of the association between spinal segmentation abnormalities and TCS. Genetic heterogeneity and the rarity of large families with KFS have limited the ability to identify mutations by traditional methods. Whole exome sequencing identified a nonsynonymous mutation in POLR1D (subunit of RNA polymerase I and II): exon2:c.T332C:p.L111P. Mutations in POLR1D are present in about 5% of individuals diagnosed with TCS. We propose that this mutation is causal in this family, suggesting a pathogenetic link between KFS and TCS.

PG4KDS: a model for the clinical implementation of pre-emptive pharmacogenetics.

Pharmacogenetics is frequently cited as an area for initial focus of the clinical implementation of genomics. Through the PG4KDS protocol, St. Jude Children's Research Hospital pre-emptively genotypes patients for 230 genes using the Affymetrix Drug Metabolizing Enzymes and Transporters (DMET) Plus array supplemented with a CYP2D6 copy number assay. The PG4KDS protocol provides a rational, stepwise process for implementing gene/drug pairs, organizing data, and obtaining consent from patients and families. Through August 2013, 1,559 patients have been enrolled, and four gene tests have been released into the electronic health record (EHR) for clinical implementation: TPMT, CYP2D6, SLCO1B1, and CYP2C19. These genes are coupled to 12 high-risk drugs. Of the 1,016 patients with genotype test results available, 78% of them had at least one high-risk (i.e., actionable) genotype result placed in their EHR. Each diplotype result released to the EHR is coupled with an interpretive consult that is created in a concise, standardized format. To support-gene based prescribing at the point of care, 55 interruptive clinical decision support (CDS) alerts were developed. Patients are informed of their genotyping result and its relevance to their medication use through a letter. Key elements necessary for our successful implementation have included strong institutional support, a knowledgeable clinical laboratory, a process to manage any incidental findings, a strategy to educate clinicians and patients, a process to return results, and extensive use of informatics, especially CDS. Our approach to pre-emptive clinical pharmacogenetics has proven feasible, clinically useful, and scalable.