Homozygote CRIM1 variant is associated with thiopurine-induced neutropenia in leukemic patients with both wildtype NUDT15 and TPMT.

BACKGROUND: NUDT15 and TPMT variants are strong genetic determinants of thiopurine-induced hematological toxicity that results in therapeutic failure in pediatric acute lymphoblastic leukemia (ALL). However, many patients with both wild-type (WT) NUDT15 and TPMT still suffer from thiopurine toxicity and therapeutic failure. METHODS: Whole-exome sequencing was done for discovery (N = 244) and replication (N = 76) cohorts. Age- and sex-adjusted multiple regression analyses of both WT patients were performed to identify (p < 0.01, N = 188 for discovery) and validate (p < 0.05, N = 52 for replication) candidate variants for the tolerated last-cycle 6-mercaptopurine (6-MP) dose intensity percentage (DIP). Both independent and additive effects of the candidate variants on well-known NUDT15 and TPMT were evaluated by multigene prediction models. RESULTS: Among the 12 candidate variants from the discovery phase, the rs3821169 variant of the gene encoding Cysteine-Rich Transmembrane BMP Regulator 1 (CRIM1) was successfully replicated (p < 0.05). It showed high interethnic variability with an impressively high allele frequency in East Asians (T = 0.255) compared to Africans (0.001), Americans (0.02), Europeans (0.009), and South Asians (0.05). Homozygote carriers of the CRIM1 rs3821169 variant (N = 12, 5%) showed significantly lower last-cycle 6-MP DIPs in the discovery, replication, and combined cohorts (p = 0.025, 0.013, and 0.001, respectively). The traditional two-gene model (NUDT15 and TPMT) for predicting 6-MP DIP < 25% was outperformed by the three-gene model that included CRIM1, in terms of the area under the receiver operating characteristic curve (0.734 vs. 0.665), prediction accuracy (0.759 vs. 0.756), sensitivity (0.636 vs. 0.523), positive predictive value (0.315 vs. 0.288), and negative predictive value (0.931 vs. 0.913). CONCLUSIONS: The CRIM1 rs3821169 variant is suggested to be an independent and/or additive genetic determinant of thiopurine toxicity beyond NUDT15 and TPMT in pediatric ALL.

Development of an After-School Day Treatment Program for Children with Severe Emotional/Behavioral Disturbances.

This article introduces Aizone, an after-school day treatment program for children with emotional/behavioral disturbances, and examines its effectiveness. Aizone targets children whose problems require intensive care, but are not severe enough to warrant hospitalization, and prioritizes those from low-income families. Based on the bio-psycho-social model, this program provides a variety of therapeutic interventions for both children and their parents, including individual and group therapy integrated under the concept of milieu therapy, as well as community linkage services. This study uses a prospective comparative method to examine the effectiveness of Aizone by comparing and validating the treatment outcomes of seven Aizone programs across Seoul. The Korean Child Behavior Checklist was administered by caregivers before and after intervention and the data was used in a repeated measures analysis. The effectiveness of the program in terms of reduction in emotional, behavioral symptoms, as well as the implications of the results are discussed.

Genetic Markers for Later Remission in Response to Early Improvement of Antidepressants.

Planning subsequent treatment strategies based on early responses rather than waiting for delayed antidepressant action can be helpful. We identified genetic markers for later non-remission in patients exhibiting poor early improvement using whole-exome sequencing data of depressive patients treated in a naturalistic manner. Among 1000 patients, early improvement at 2 weeks (reduction in Hamilton Depression Rating Scale [HAM-D] score ≥ 20%) and remission at 12 weeks (HAM-D score ≤ 7) were evaluated. Gene- and variant-level analyses were conducted to compare patients who did not exhibit early improvement and did not eventually achieve remission (n = 126) with those who exhibited early improvement and achieved remission (n = 385). Genes predicting final non-remission in patients who exhibited poor early improvement (COMT, PRNP, BRPF3, SLC25A40, and CGREF1 in males; PPFIBPI, LZTS3, MEPCE, MAP1A, and PFAS in females; ST3GAL5 in the total population) were determined. Among the significant genes, variants in the PRNP (rs1800014), COMT (rs6267), BRPF3 (rs200565609), and SLC25A40 genes (rs3213633) were identified. However, interpretations should be made cautiously, as complex pharmacotherapy involves various genes and pathways. Early detection of poor early improvement and final non-remission based on genetic risk would be helpful for decision-making in a clinical setting.

Contradiction in Star-Allele Nomenclature of Pharmacogenes between Common Haplotypes and Rare Variants.

The nomenclature of star alleles has been widely used in pharmacogenomics to enhance treatment outcomes, predict drug response variability, and reduce adverse reactions. However, the discovery of numerous rare functional variants through genome sequencing introduces complexities into the star-allele system. This study aimed to assess the nature and impact of the rapid discovery of numerous rare functional variants in the traditional haplotype-based star-allele system. We developed a new method to construct haplogroups, representing a common ancestry structure, by iteratively excluding rare and functional variants of the 25 representative pharmacogenes using the 2504 genomes from the 1000 Genomes Project. In total, 192 haplogroups and 288 star alleles were identified, with an average of 7.68 ± 4.2 cross-ethnic haplogroups per gene. Most of the haplogroups (70.8%, 136/192) were highly aligned with their corresponding classical star alleles (VI = 1.86 ± 0.78), exhibiting higher genetic diversity than the star alleles. Approximately 41.3% (N = 119) of the star alleles in the 2504 genomes did not belong to any of the haplogroups, and most of them (91.3%, 105/116) were determined by a single variant according to the allele-definition table provided by CPIC. These functional single variants had low allele frequency (MAF < 1%), high evolutionary conservation, and variant deleteriousness, which suggests significant negative selection. It is suggested that the traditional haplotype-based naming system for pharmacogenetic star alleles now needs to be adjusted by balancing both traditional haplotyping and newly emerging variant-sequencing approaches to reduce naming complexity.

Characteristics of Comorbid Physical Disease in Patients With Severe Mental Illness in South Korea: A Nationwide Population-Based Study (2014-2019).

OBJECTIVE: The purpose of this study was to identify the associations of chronic physical disease between patients with severe mental illness (SMI) and the general population of South Korea. METHODS: This study was conducted with National Health Insurance Corporation data from 2014 to 2019. A total of 848,058 people were diagnosed with SMI in this period, and the same number of controls were established by matching by sex and age. A descriptive analysis was conducted on the sociodemographic characteristics of patients with SMI. Conditional logistic regression analysis was performed to identify the associations between comorbid physical disease in patients with SMI and those of the general population. SAS Enterprise Guide 7.1 (SAS Inc, Cary, NC, USA) were used to perform all statistical tests. RESULTS: The analysis revealed significant differences in medical insurance, income level, and Charlson Comorbidity Index weighted by chronic physical disease, between patients with SMI and the general population. Conditional logistic regression analysis between the two groups also revealed significant differences in eight chronic physical diseases except hypertensive disease. CONCLUSION: This study confirmed the vulnerability of patients with SMI to chronic physical diseases and we were able to identify chronic physical disease that were highly related to patients with SMI.

Towards more accurate pharmacogenomic variant effect predictions.

Tweetable abstract Accurate variant interpretation has become a key bottleneck for the translation of an individual's pharmacogenome into actionable recommendations. We recommend an integrated use of multiplexed assays, structure-based predictions and biobank data to develop more accurate effect predictors.

Effectiveness of a Community-Based Intensive Case Management Model on Reducing Hospitalization for People With Severe Mental Illness in Seoul.

OBJECTIVE: To overcome the limited effectiveness of standard case management services, the Seoul Intensive Case Management program (S-ICM) for patients with serious mental illnesses was introduced in 2017. This study aimed to evaluate its effectiveness in reducing the length of hospital stay. METHODS: Monitoring data from April 2019 to March 2020 were retrieved from the Seoul Mental Health Welfare Center. A total of 759 participants with serious mental illnesses were included. The average length of admission per month was compared between the pre-ICM (previous year) and during-ICM periods. For post-ICM observation subgroup, average length of admission per month was compared between pre-ICM, during-ICM, and post-ICM periods. To determine the relative contributions of risk factors for during-ICM and post-ICM admission, multivariate logistic regression analyses were performed. RESULTS: The average admission stay for pre-ICM period was significantly longer than that for during-ICM period (1.47 vs. 0.26 days). Among the predictors for during-ICM admission, pre-ICM psychiatric admission was the most important risk factor, followed by medical aid beneficiary and suicidal behavior. In the subgroup analysis of the post-ICM observation period, the pre-ICM, during-ICM, and post-ICM average admission stays were 1.45, 0.29, and 0.57 days/month, respectively. There was a significant difference in the average length of stay between the pre-ICM and during-ICM periods and between the pre-ICM and post-ICM periods. Post-ICM admission risks included pre-ICM admission, S-ICM duration <3 months, and chronic unstable symptoms. CONCLUSION: The results suggest that the S-ICM effectively reduces psychiatric hospitalization duration, at least over a short-term period.

Genetic markers for depressive disorders with earlier age at onset.

Age at onset has been considered a potential indicator of underlying genetic risk in depression research. However, the variants associated with earlier age at onset of depressive disorder have not been elucidated. To evaluate the genetic architecture of depression onset, whole-exome sequencing of samples from 1000 patients with depressive disorder was performed. Cox proportional hazard models with false discovery rate-adjusted P-values were used to estimate the hazard ratios; carriers and non-carriers of individual coding variants were compared in terms of age at onset of depression with adjustment for sociodemographic and clinical characteristics. The clinical relevance of the candidate variants was also examined. Whole-exome sequencing revealed four variants in the CCL14, FYB, GPRASP1, and CTNND2 genes associated with an increased risk of depressive disorder with earlier age at onset. Although no individual variant was associated with any clinical characteristic except AAO, together they were associated with younger AAO, younger age at visit for treatment, and recurrent and atypical depression. Our data suggest novel candidate genes for depressive disorder with earlier age at onset. These genes could serve as markers allowing early identification of patients at risk of depression, and thus earlier intervention.

Interplay between IL6 and CRIM1 in thiopurine intolerance due to hematological toxicity in leukemic patients with wild-type NUDT15 and TPMT.

NUDT15 and TPMT variants are strong genetic determinants of thiopurine-induced hematological toxicity. Despite the impact of homozygous CRIM1 on thiopurine toxicity, several patients with wild-type NUDT15, TPMT, and CRIM1 experience thiopurine toxicity, therapeutic failure, and relapse of acute lymphoblastic leukemia (ALL). Novel pharmacogenetic interactions associated with thiopurine intolerance from hematological toxicities were investigated using whole-exome sequencing for last-cycle 6-mercaptopurine dose intensity percentages (DIP) tolerated by pediatric ALL patients (N = 320). IL6 rs13306435 carriers (N = 19) exhibited significantly lower DIP (48.0 ± 27.3%) than non-carriers (N = 209, 69.9 ± 29.0%; p = 0.0016 and 0.0028 by t test and multiple linear regression, respectively). Among 19 carriers, 7 with both heterozygous IL6 rs13306435 and CRIM1 rs3821169 showed significantly decreased DIP (24.7 ± 8.9%) than those with IL6 (N = 12, 61.6 ± 25.1%) or CRIM1 (N = 94, 68.1 ± 28.4%) variants. IL6 and CRIM1 variants showed marked inter-ethnic variability. Four-gene-interplay models revealed the best odds ratio (8.06) and potential population impact [relative risk (5.73), population attributable fraction (58%), number needed to treat (3.67), and number needed to genotype (12.50)]. Interplay between IL6 rs13306435 and CRIM1 rs3821169 was suggested as an independent and/or additive genetic determinant of thiopurine intolerance beyond NUDT15 and TPMT in pediatric ALL.

Gene-wise variant burden and genomic characterization of nearly every gene.

Aim: Current gene-level prioritization methods aim to provide information for further prioritization of 'disease-causing' mutations. Since, they are inherently biased toward disease genes, methods specific to pharmacogenetic (PGx) genes are required. Methods: We proposed a gene-wise variant burden (GVB) method that integrates in silico deleteriousness scores of the multitude of variants of a given gene at a personal-genome level. Results: GVB in its simplest form outperformed the two state-of-the-art methods with regard to predicting pharmacogenes and complex disease genes but not for rare Mendelian disease genes. GVB* adjusted by paralog counts robustly performed well in most of the pharmacogenetic subcategories. Seven molecular genetic features well characterized the unique genomic properties of PGx, complex, and Mendelian disease genes. Conclusion: Altogether, GVB is an individual-specific genescore, especially advantageous for PGx studies.

Clinical Genome Data Model (cGDM) provides Interactive Clinical Decision Support for Precision Medicine.

In light of recent developments in genomic technology and the rapid accumulation of genomic information, a major transition toward precision medicine is anticipated. However, the clinical applications of genomic information remain limited. This lag can be attributed to several complex factors, including the knowledge gap between medical experts and bioinformaticians, the distance between bioinformatics workflows and clinical practice, and the unique characteristics of genomic data, which can make interpretation difficult. Here we present a novel genomic data model that allows for more interactive support in clinical decision-making. Informational modelling was used as a basis to design a communication scheme between sophisticated bioinformatics predictions and the representative data relevant to a clinical decision. This study was conducted by a multidisciplinary working group who carried out clinico-genomic workflow analysis and attribute extraction, through Failure Mode and Effects Analysis (FMEA). Based on those results, a clinical genome data model (cGDM) was developed with 8 entities and 46 attributes. The cGDM integrates reliability-related factors that enable clinicians to access the reliability problem of each individual genetic test result as clinical evidence. The proposed cGDM provides a data-layer infrastructure supporting the intellectual interplay between medical experts and informed decision-making.

Sex differences in the genetic architecture of depression.

The prevalence and clinical characteristics of depressive disorders differ between women and men; however, the genetic contribution to sex differences in depressive disorders has not been elucidated. To evaluate sex-specific differences in the genetic architecture of depression, whole exome sequencing of samples from 1000 patients (70.7% female) with depressive disorder was conducted. Control data from healthy individuals with no psychiatric disorder (n = 72, 26.4% female) and East-Asian subpopulation 1000 Genome Project data (n = 207, 50.7% female) were included. The genetic variation between men and women was directly compared using both qualitative and quantitative research designs. Qualitative analysis identified five genetic markers potentially associated with increased risk of depressive disorder in females, including three variants (rs201432982 within PDE4A, and rs62640397 and rs79442975 within FDX1L) mapping to chromosome 19p13.2 and two novel variants (rs820182 and rs820148) within MYO15B at the chromosome 17p25.1 locus. Depressed patients homozygous for these variants showed more severe depressive symptoms and higher suicidality than those who were not homozygotes (i.e., heterozygotes and homozygotes for the non-associated allele). Quantitative analysis demonstrated that the genetic burden of protein-truncating and deleterious variants was higher in males than females, even after permutation testing. Our study provides novel genetic evidence that the higher prevalence of depressive disorders in women may be attributable to inherited variants.

Identifying genetic variants associated with ritodrine-induced pulmonary edema.

INTRODUCTION: Ritodrine is one of the most commonly used tocolytics in preterm labor, acting as a ß2-adrenergic agonist that reduces intracellular calcium levels and prevents myometrial activation. Ritodrine infusion can result in serious maternal complications, and pulmonary edema is a particular concern among these. The cause of pulmonary edema following ritodrine treatment is multifactorial; however, the contributing genetic factors remain poorly understood. This study investigates the genetic variants associated with ritodrine-induced pulmonary edema. METHODS: In this case-control study, 16 patients who developed pulmonary edema during ritodrine infusion [case], and 16 pregnant women who were treated with ritodrine and did not develop pulmonary edema [control] were included. The control pregnant women were selected after matching for plurality and gestational age at the time of tocolytic use. Maternal blood was collected during admission for tocolytic treatment, and whole exome sequencing was performed with the stored blood samples. RESULTS: Gene-wise variant burden (GVB) analysis resulted in a total of 71 candidate genes by comparing the cumulative effects of multiple coding variants for 19729 protein-coding genes between the patients with pulmonary edema and the matched controls. Subsequent data analysis selected only the statistically significant and deleterious variants compatible with ritodrine-induced pulmonary edema. Two final candidate variants in CPT2 and ADRA1A were confirmed by Sanger sequencing. CONCLUSIONS: We identified new potential variants in genes that play a role in cyclic adenosine monophosphate (cAMP)/protein kinase A (PKA) regulation, which supports their putative involvement in the predisposition to ritodrine-induced pulmonary edema in pregnant women.

Star Allele-Based Haplotyping versus Gene-Wise Variant Burden Scoring for Predicting 6-Mercaptopurine Intolerance in Pediatric Acute Lymphoblastic Leukemia Patients.

Nudix Hydrolase 15 (NUDT15) and Thiopurine S-Methyltransferase (TPMT) are strong genetic determinants of thiopurine toxicity in pediatric acute lymphoblastic leukemia (ALL) patients. Since patients with NUDT15 or TPMT deficiency suffer severe adverse drug reactions, star (*) allele-based haplotypes have been used to predict an optimal 6-mercaptopurine (6-MP) dosing. However, star allele haplotyping suffers from insufficient, inconsistent, and even conflicting designations with uncertain and/or unknown functional alleles. Gene-wise variant burden (GVB) scoring enables us to utilize next-generation sequencing (NGS) data to predict 6-MP intolerance in children with ALL. Whole exome sequencing was performed for 244 pediatric ALL patients under 6-MP treatments. We assigned star alleles with PharmGKB haplotype set translational table. GVB for NUDT15 and TPMT was computed by aggregating in silico deleteriousness scores of multiple coding variants for each gene. Poor last-cycle dose intensity percent (DIP < 25%) was considered as 6-MP intolerance, resulting therapeutic failure of ALL. DIPs showed significant differences ( p < 0.05) among NUDT15 poor (PM, n = 1), intermediate (IM, n = 48), and normal (NM, n = 195) metabolizers. TPMT exhibited no PM and only seven IMs. GVB showed significant differences among the different haplotype groups of both NUDT15 and TPMT ( p < 0.05). Kruskal-Wallis test for DIP values showed statistical significances for the seven different GVB score bins of NUDT15. GVB NUDT15 outperformed the star allele-based haplotypes in predicting patients with reduced last-cycle DIPs at all DIP threshold levels (i.e., 5%, 10%, 15%, and 25%). In NUDT15-and-TPMT combined interaction analyses, GVB NUDT15 , TPMT outperformed star alleles [area under the receiver operating curve (AUROC) = 0.677 vs. 0.645] in specificity (0.813 vs. 0.796), sensitivity (0.526 vs. 0.474), and positive (0.192 vs. 0.164) and negative (0.953 vs. 0.947) predictive values. Overall, GVB correctly classified five more patients (i.e., one into below and four into above 25% DIP groups) than did star allele haplotypes. GVB analysis demonstrated that 6-MP intolerance in pediatric ALL can be reliably predicted by aggregating NGS-based common, rare, and novel variants together without hampering the predictive power of the conventional haplotype analysis.

Discovery of donor genotype associated with long-term survival of patients with hematopoietic stem cell transplantation in refractory acute myeloid leukemia.

Allogeneic hematopoietic stem cell transplantation (allo-HSCT) has been the only treatment option for acute myeloid leukemia (AML) refractory to induction chemotherapy, with only 10-20% of patients achieving long-term survival. Certain donor genotypes may confer leukemia-clearing effects after allo-HSCT. We performed whole-exome sequencing of five pairs of the germ lines in AML patients who achieved long-term remission after allo-HSCT and in their donors, and found two significant variants: EGFR c.2982C > T and CDH11 c.945G > A. To validate the protective effects of these leukemia-clearing genotypes (LCGs), AML patients who received allo-HSCT in a complete-remission status were also analyzed. Twenty-two of 96 donors (22.9%) had LCGs in their genomes, and overall survival was significantly longer in patients who received allo-HSCT from donors with germ-line LCGs (hazard ratio=0.47, 95% confidence interval=0.24-0.94, p = .033). These findings indicate that donor germ-line LCGs have phenotypically leukemia-clearing effects and are biomarkers for predicting clinical outcomes in allogeneic transplantation in AML patients.

Molecular diagnosis of hereditary spherocytosis by multi-gene target sequencing in Korea: matching with osmotic fragility test and presence of spherocyte.

BACKGROUND: Current diagnostic tests for hereditary spherocytosis (HS) focus on the detection of hemolysis or indirectly assessing defects of membrane protein, whereas direct methods to detect protein defects are complicated and difficult to implement. In the present study, we investigated the patterns of genetic variation associated with HS among patients clinically diagnosed with HS. METHODS: Multi-gene targeted sequencing of 43 genes (17 RBC membrane protein-encoding genes, 20 RBC enzyme-encoding genes, and six additional genes for the differential diagnosis) was performed using the Illumina HiSeq platform. RESULTS: Among 59 patients with HS, 50 (84.7%) had one or more significant variants in a RBC membrane protein-encoding genes. A total of 54 significant variants including 46 novel mutations were detected in six RBC membrane protein-encoding genes, with the highest number of variants found in SPTB (n = 28), and followed by ANK1 (n = 19), SLC4A1 (n = 3), SPTA1 (n = 2), EPB41 (n = 1), and EPB42 (n = 1). Concurrent mutations of genes encoding RBC enzymes (ALDOB, GAPDH, and GSR) were detected in three patients. UGT1A1 mutations were present in 24 patients (40.7%). Positive rate of osmotic fragility test was 86.8% among patients harboring HS-related gene mutations. CONCLUSIONS: This constitutes the first large-scaled genetic study of Korean patients with HS. We demonstrated that multi-gene target sequencing is sensitive and feasible that can be used as a powerful tool for diagnosing HS. Considering the discrepancies of clinical and molecular diagnoses of HS, our findings suggest that molecular genetic analysis is required for accurate diagnosis of HS.

APEX1 Polymorphism and Mercaptopurine-Related Early Onset Neutropenia in Pediatric Acute Lymphoblastic Leukemia.

PURPOSE: Mercaptopurine (MP) is one of the main chemotherapeutics for acute lymphoblastic leukemia (ALL). To improve treatment outcomes, constant MP dose titration is essential to maintain steady drug exposure, while minimizing myelosuppression. We performed two-stage analyses to identify genetic determinants of MP-related neutropenia in Korean pediatric ALL patients. MATERIALS AND METHODS: Targeted sequencing of 40 patients who exhibited definite MP intolerance was conducted using a novel panel of 211 pharmacogenetic-related genes, and subsequent analysis was performed with 185 patients. RESULTS: Using bioinformatics tools and genetic data, four functionally interesting variants were selected (ABCC4, APEX1, CYP1A1, and CYP4F2). Including four variants, 23 variants in 12 genes potentially linked to MP adverse reactions were selected as final candidates for subsequent analysis in 185 patients. Ultimately, a variant allele in APEX1 rs2307486was found to be strongly associated with MP-induced neutropenia that occurred within 28 days of initiating MP (odds ratio, 3.44; p=0.02). Moreover, the cumulative incidence of MP-related neutropenia was significantly higher in patients with APEX1 rs2307486 variants, as GG genotypes were associated with the highest cumulative incidence (p < 0.01). NUDT15 rs116855232 variants were strongly associated with a higher cumulative incidence of neutropenia (p < 0.01), and a lower median dose of tolerated MP throughout maintenance treatment (p < 0.01). CONCLUSION: We have identified that APEX1 rs2307486 variants conferred an increased risk of MP-related early onset neutropenia. APEX1 and NUDT15 both contribute to cell protection from DNA damage or misincorporation, so alleles that impair the function of either gene may affect MP sensitivities, thereby inducing MP-related neutropenia.

Deleterious genetic variants in ciliopathy genes increase risk of ritodrine-induced cardiac and pulmonary side effects.

BACKGROUND: Ritodrine is a commonly used tocolytic to prevent preterm labour. However, it can cause unexpected serious adverse reactions, such as pulmonary oedema, pulmonary congestion, and tachycardia. It is unknown whether such adverse reactions are associated with pharmacogenomic variants in patients. METHODS: Whole-exome sequencing of 13 subjects with serious ritodrine-induced cardiac and pulmonary side-effects was performed to identify causal genes and variants. The deleterious impact of nonsynonymous substitutions for all genes was computed and compared between cases (n = 13) and controls (n = 30). The significant genes were annotated with Gene Ontology (GO), and the associated disease terms were categorised into four functional classes for functional enrichment tests. To assess the impact of distributed rare variants in cases with side effects, we carried out rare variant association tests with a minor allele frequency ≤ 1% using the burden test, the sequence Kernel association test (SKAT), and optimised SKAT. RESULTS: We identified 28 genes that showed significantly lower gene-wise deleteriousness scores in cases than in controls. Three of the identified genes-CYP1A1, CYP8B1, and SERPINA7-are pharmacokinetic genes. The significantly identified genes were categorized into four functional classes: ion binding, ATP binding, Ca2+-related, and ciliopathies-related. These four classes were significantly enriched with ciliary genes according to SYSCILIA Gold Standard genes (P < 0.01), thus representing ciliary genes. Furthermore, SKAT showed a marginal trend toward significance after Bonferroni correction with Joubert Syndrome ciliopathy genes (P = 0.05). With respect to the pharmacokinetic genes, rs1048943 (CYP1A1) and rs1804495 (SERPINA7) showed a significantly higher frequency in cases than controls, as determined by Fisher's exact test (P < 0.05 and P < 0.01, respectively). CONCLUSIONS: Ritodrine-induced cardiac and pulmonary side effects may be associated with deleterious genetic variants in ciliary and pharmacokinetic genes.

Evaluation of exome variants using the Ion Proton Platform to sequence error-prone regions.

The Ion Proton sequencer from Thermo Fisher accurately determines sequence variants from target regions with a rapid turnaround time at a low cost. However, misleading variant-calling errors can occur. We performed a systematic evaluation and manual curation of read-level alignments for the 675 ultrarare variants reported by the Ion Proton sequencer from 27 whole-exome sequencing data but that are not present in either the 1000 Genomes Project and the Exome Aggregation Consortium. We classified positive variant calls into 393 highly likely false positives, 126 likely false positives, and 156 likely true positives, which comprised 58.2%, 18.7%, and 23.1% of the variants, respectively. We identified four distinct error patterns of variant calling that may be bioinformatically corrected when using different strategies: simplicity region, SNV cluster, peripheral sequence read, and base inversion. Local de novo assembly successfully corrected 201 (38.7%) of the 519 highly likely or likely false positives. We also demonstrate that the two sequencing kits from Thermo Fisher (the Ion PI Sequencing 200 kit V3 and the Ion PI Hi-Q kit) exhibit different error profiles across different error types. A refined calling algorithm with better polymerase may improve the performance of the Ion Proton sequencing platform.