Systematic immune cell dysregulation and molecular subtypes revealed by single-cell RNA-seq of subjects with type 1 diabetes.

BACKGROUND: Type 1 diabetes mellitus (T1DM) is a prototypic endocrine autoimmune disease resulting from an immune-mediated destruction of pancreatic insulin-secreting ß  cells. A comprehensive immune cell phenotype evaluation in T1DM has not been performed thus far at the single-cell level. METHODS: In this cross-sectional analysis, we generated a single-cell transcriptomic dataset of peripheral blood mononuclear cells (PBMCs) from 46 manifest T1DM (stage 3) cases and 31 matched controls. RESULTS: We surprisingly detected profound alterations in circulatory immune cells (1784 dysregulated genes in 13 immune cell types), far exceeding the count in the comparator systemic autoimmune disease SLE. Genes upregulated in T1DM were involved in WNT signaling, interferon signaling and migration of T/NK cells, antigen presentation by B cells, and monocyte activation. A significant fraction of these differentially expressed genes were also altered in T1DM pancreatic islets. We used the single-cell data to construct a T1DM metagene z-score (TMZ score) that distinguished cases and controls and classified patients into molecular subtypes. This score correlated with known prognostic immune markers of T1DM, as well as with drug response in clinical trials. CONCLUSIONS: Our study reveals a surprisingly strong systemic dimension at the level of immune cell network in T1DM, defines disease-relevant molecular subtypes, and has the potential to guide non-invasive test development and patient stratification.

Kidney organoid models reveal cilium-autophagy metabolic axis as a therapeutic target for PKD both in vitro and in vivo.

Human pluripotent stem cell-derived kidney organoids offer unprecedented opportunities for studying polycystic kidney disease (PKD), which still has no effective cure. Here, we developed both in vitro and in vivo organoid models of PKD that manifested tubular injury and aberrant upregulation of renin-angiotensin aldosterone system. Single-cell analysis revealed that a myriad of metabolic changes occurred during cystogenesis, including defective autophagy. Experimental activation of autophagy via ATG5 overexpression or primary cilia ablation significantly inhibited cystogenesis in PKD kidney organoids. Employing the organoid xenograft model of PKD, which spontaneously developed tubular cysts, we demonstrate that minoxidil, a potent autophagy activator and an FDA-approved drug, effectively attenuated cyst formation in vivo. This in vivo organoid model of PKD will enhance our capability to discover novel disease mechanisms and validate candidate drugs for clinical translation.

Rodent genetically modified models of glaucoma.

Glaucoma, one of the leading causes of irreversible blindness worldwide, is a complex and heterogenous disease. While environmental factors are important, it is well-recognized that the disease has a strong heritable component. With the advent of large-cohort genome wide association studies, a myriad of genetic risk loci has been linked to different forms of glaucoma. Animal models have been an indispensable tool in characterizing these loci, especially if they lie within coding regions in the genome. Not only do these models connect genotype to phenotype, advancing our understanding of glaucoma pathogenesis in the process, they also have valuable utility as a platform for the pre-clinical testing of potential therapies. In this review, we will outline genetic models used for studying the major forms of glaucoma, including primary open angle glaucoma, normal tension glaucoma, primary angle closure glaucoma, pigmentary glaucoma, pseudoexfoliation glaucoma, and early onset glaucoma, including congenital and developmental glaucoma, and how studying these models have helped shed light on human glaucoma.

APOE Genotype Modifies the Association between Midlife Adherence to the Planetary Healthy Diet and Cognitive Function in Later Life among Chinese Adults in Singapore.

BACKGROUND: It remains unclear if adherence to the planetary healthy diet (PHD), designed to improve human and environmental health, is associated with better cognitive function in aging, and if this association differs by apolipoprotein E (APOE) genotype. OBJECTIVES: We aimed to examine the association between the PHD pattern and risk of poor cognitive function, and to further assess whether the APOE ε4 allele could modify this association. METHODS: The study included 16,736 participants from the Singapore Chinese Health Study. The PHD score was calculated using data from a validated 165-item food frequency questionnaire at baseline (1993-1998), with higher scores indicating greater adherence to the PHD. Cognitive function was assessed by the Singapore-modified Mini-Mental State Examination at follow-up 3 visits (2014-2016). A subset of 9313 participants had APOE genotype data. Logistic regression models were used to estimate the odds ratios (ORs) and 95% confidence intervals (CIs), with adjustment for potential confounders. RESULTS: We identified 2397 (14.3%) cases of poor cognitive function. In the total population, OR (95% CI) of poor cognitive function for each one-SD increment in the PHD score was 0.89 (0.85, 0.93). Carriers of APOE ε4 allele had increased risk of poor cognitive function (OR: 1.36, 95% CI: 1.15, 1.61). There was a significant interaction between the PHD score and the APOE ε4 allele (P-interaction = 0.042). Each one-SD increment in the PHD score was significantly associated with lower risk of poor cognitive function (OR: 0.89; 95% CI: 0.83, 0.96) in non-carriers of APOE ε4 allele, but not in APOE ε4 allele carriers (OR: 1.04, 95% CI: 0.89, 1.23). CONCLUSIONS: Midlife adherence to the PHD was associated with reduced risk of poor cognitive function in later life. However, this was not observed in carriers of APOE ε4 allele who had higher risk of poor cognitive function.

Genome-Wide Association Analysis of Protein-Coding Variants in IgA Nephropathy.

SIGNIFICANCE STATEMENT: Genome-wide association studies have identified nearly 20 IgA nephropathy susceptibility loci. However, most nonsynonymous coding variants, particularly ones that occur rarely or at a low frequency, have not been well investigated. The authors performed a chip-based association study of IgA nephropathy in 8529 patients with the disorder and 23,224 controls. They identified a rare variant in the gene encoding vascular endothelial growth factor A (VEGFA) that was significantly associated with a two-fold increased risk of IgA nephropathy, which was further confirmed by sequencing analysis. They also identified a novel common variant in PKD1L3 that was significantly associated with lower haptoglobin protein levels. This study, which was well-powered to detect low-frequency variants with moderate to large effect sizes, helps expand our understanding of the genetic basis of IgA nephropathy susceptibility. BACKGROUND: Genome-wide association studies have identified nearly 20 susceptibility loci for IgA nephropathy. However, most nonsynonymous coding variants, particularly those occurring rarely or at a low frequency, have not been well investigated. METHODS: We performed a three-stage exome chip-based association study of coding variants in 8529 patients with IgA nephropathy and 23,224 controls, all of Han Chinese ancestry. Sequencing analysis was conducted to investigate rare coding variants that were not covered by the exome chip. We used molecular dynamic simulation to characterize the effects of mutations of VEGFA on the protein's structure and function. We also explored the relationship between the identified variants and the risk of disease progression. RESULTS: We discovered a novel rare nonsynonymous risk variant in VEGFA (odds ratio, 1.97; 95% confidence interval [95% CI], 1.61 to 2.41; P = 3.61×10 -11 ). Further sequencing of VEGFA revealed twice as many carriers of other rare variants in 2148 cases compared with 2732 controls. We also identified a common nonsynonymous risk variant in PKD1L3 (odds ratio, 1.16; 95% CI, 1.11 to 1.21; P = 1.43×10 -11 ), which was associated with lower haptoglobin protein levels. The rare VEGFA mutation could cause a conformational change and increase the binding affinity of VEGFA to its receptors. Furthermore, this variant was associated with the increased risk of kidney disease progression in IgA nephropathy (hazard ratio, 2.99; 95% CI, 1.09 to 8.21; P = 0.03). CONCLUSIONS: Our study identified two novel risk variants for IgA nephropathy in VEGFA and PKD1L3 and helps expand our understanding of the genetic basis of IgA nephropathy susceptibility.

Disruption of Aldehyde Dehydrogenase 2 protects against bacterial infection.

The ALDH2*2 (rs671) allele is one of the most common genetic mutations in humans, yet the positive evolutionary selective pressure to maintain this mutation is unknown, despite its association with adverse health outcomes. ALDH2 is responsible for the detoxification of metabolically produced aldehydes, including lipid-peroxidation end products derived from inflammation. Here, we demonstrate that host-derived aldehydes 4-hydroxynonenal (4HNE), malondialdehyde (MDA), and formaldehyde (FA), all of which are metabolized by ALDH2, are directly toxic to the bacterial pathogens Mycobacterium tuberculosis and Francisella tularensis at physiological levels. We find that Aldh2 expression in macrophages is decreased upon immune stimulation, and that bone marrow-derived macrophages from Aldh2 -/- mice contain elevated aldehydes relative to wild-type mice. Macrophages deficient for Aldh2 exhibited enhanced control of Francisella infection. Finally , mice lacking Aldh2 demonstrated increased resistance to pulmonary infection by M. tuberculosis , including in a hypersusceptible model of tuberculosis, and were also resistant to Francisella infection. We hypothesize that the absence of ALDH2 contributes to the host's ability to control infection by pathogens such as M. tuberculosis and F. tularensis , and that host-derived aldehydes act as antimicrobial factors during intracellular bacterial infections. One sentence summary: Aldehydes produced by host cells contribute to the control of bacterial infections.

Direct inference and control of genetic population structure from RNA sequencing data.

RNAseq data can be used to infer genetic variants, yet its use for estimating genetic population structure remains underexplored. Here, we construct a freely available computational tool (RGStraP) to estimate RNAseq-based genetic principal components (RG-PCs) and assess whether RG-PCs can be used to control for population structure in gene expression analyses. Using whole blood samples from understudied Nepalese populations and the Geuvadis study, we show that RG-PCs had comparable results to paired array-based genotypes, with high genotype concordance and high correlations of genetic principal components, capturing subpopulations within the dataset. In differential gene expression analysis, we found that inclusion of RG-PCs as covariates reduced test statistic inflation. Our paper demonstrates that genetic population structure can be directly inferred and controlled for using RNAseq data, thus facilitating improved retrospective and future analyses of transcriptomic data.

Third-Generation Single-Molecule Sequencing for Preimplantation Genetic Testing of Aneuploidy and Segmental Imbalances.

BACKGROUND: Current strategies for preimplantation genetic testing for aneuploidy or structural rearrangements (PGT-A/SR) rely mainly on next-generation sequencing (NGS) and microarray platforms, which are robust but require expensive instrumentation. We explored the suitability of third-generation single-molecule sequencing as a PGT-A/SR screening platform for both aneuploidy and segmental imbalance. METHODS: Single-cell and multicell replicates from aneuploid or segmentally unbalanced cell lines (n = 208) were SurePlex-amplified, randomized, and subjected to (a) Nanopore-based single-molecule sequencing (Oxford Nanopore Technologies) and (b) NGS using a leading commercial PGT-A solution (Illumina VeriSeq PGS). Archival SurePlex-amplified trophectoderm biopsy samples (n = 96) previously analyzed using the commercial kit were blinded and reanalyzed using Nanopore. RESULTS: Nanopore-based PGT-A identified the specific aberration in 95.45% (84/88) and 97.78% (88/90) of single-/multicells with an aneuploidy or segmental imbalance (10-30.5 Mb), respectively. Comparison against the commercial kit's results revealed concordances of 98.86% (87/88) and 98.89% (89/90) for the aneuploid and segmentally unbalanced (10-30.5 Mb aberration) samples, respectively. Detection sensitivity for smaller segmental imbalances (5-5.8 Mb aberration, n = 30) decreased markedly on both platforms. Nanopore-based PGT-A reanalysis of trophectoderm biopsy samples was 97.92% (94/96) concordant with the commercial kit results. CONCLUSION: Up to 24 SurePlex-amplified single-cell, multicell, or trophectoderm samples could be sequenced in a single MinION flow-cell for subsequent preimplantation genetic testing for aneuploidy or structural rearrangements (PGT-A/SR) analysis, with results obtainable in ≤3 days and at per-sample costs that are competitive with commercial offerings. Nanopore's third-generation single-molecule sequencing represents a viable alternative to current commercial NGS-based PGT-A solutions for aneuploidy and segmental imbalance (≥10 Mb) screening of single-/multicell or trophectoderm biopsy samples.

A five-safes approach to a secure and scalable genomics data repository.

Genomic researchers increasingly utilize commercial cloud service providers (CSPs) to manage data and analytics needs. CSPs allow researchers to grow Information Technology (IT) infrastructure on demand to overcome bottlenecks when combining large datasets. However, without adequate security controls, the risk of unauthorized access may be higher for data stored on the cloud. Additionally, regulators are mandating data access patterns and specific security protocols for the storage and use of genomic data. While CSP provides tools for security and regulatory compliance, building the necessary controls required for cloud solutions is not trivial. Research Assets Provisioning and Tracking Online Repository (RAPTOR) by the Genome Institute of Singapore is a cloud-native genomics data repository and analytics platform that implements a "five-safes" framework to provide security and governance controls to data contributors and users, leveraging CSP for sharing and analysis of genomic datasets without the risk of security breaches or running afoul of regulations.

Polygenic risk scores for the prediction of common cancers in East Asians: A population-based prospective cohort study.

Background: To evaluate the utility of polygenic risk scores (PRSs) in identifying high-risk individuals, different publicly available PRSs for breast (n=85), prostate (n=37), colorectal (n=22), and lung cancers (n=11) were examined in a prospective study of 21,694 Chinese adults. Methods: We constructed PRS using weights curated in the online PGS Catalog. PRS performance was evaluated by distribution, discrimination, predictive ability, and calibration. Hazard ratios (HR) and corresponding confidence intervals (CI) of the common cancers after 20 years of follow-up were estimated using Cox proportional hazard models for different levels of PRS. Results: A total of 495 breast, 308 prostate, 332 female-colorectal, 409 male-colorectal, 181 female-lung, and 381 male-lung incident cancers were identified. The area under receiver operating characteristic curve for the best-performing site-specific PRS were 0.61 (PGS000873, breast), 0.70 (PGS00662, prostate), 0.65 (PGS000055, female-colorectal), 0.60 (PGS000734, male-colorectal), 0.56 (PGS000721, female-lung), and 0.58 (PGS000070, male-lung), respectively. Compared to the middle quintile, individuals in the highest cancer-specific PRS quintile were 64% more likely to develop cancers of the breast, prostate, and colorectal. For lung cancer, the lowest cancer-specific PRS quintile was associated with 28-34% decreased risk compared to the middle quintile. In contrast, the HR observed for quintiles 4 (female-lung: 0.95 [0.61-1.47]; male-lung: 1.14 [0.82-1.57]) and 5 (female-lung: 0.95 [0.61-1.47]) were not significantly different from that for the middle quintile. Conclusions: Site-specific PRSs can stratify the risk of developing breast, prostate, and colorectal cancers in this East Asian population. Appropriate correction factors may be required to improve calibration. Funding: This work is supported by the National Research Foundation Singapore (NRF-NRFF2017-02), PRECISION Health Research, Singapore (PRECISE) and the Agency for Science, Technology and Research (A*STAR). WP Koh was supported by National Medical Research Council, Singapore (NMRC/CSA/0055/2013). CC Khor was supported by National Research Foundation Singapore (NRF-NRFI2018-01). Rajkumar Dorajoo received a grant from the Agency for Science, Technology and Research Career Development Award (A*STAR CDA - 202D8090), and from Ministry of Health Healthy Longevity Catalyst Award (HLCA20Jan-0022).The Singapore Chinese Health Study was supported by grants from the National Medical Research Council, Singapore (NMRC/CIRG/1456/2016) and the U.S. National Institutes of Health (NIH) (R01 CA144034 and UM1 CA182876).

Although humans contain the same genes, the sequence within these DNA sites can vary from person to person. These small variations, also known as genetic variants, can increase the risk of developing certain diseases. While each variant will only have a weak effect, if multiple variations are present the odds of developing the disease becomes significantly higher. To determine which variants are linked to a disease, researchers carry out genome-wide association studies which involve analyzing the genomes of individuals with and without the condition and comparing their genetic codes. This data is then used to calculate how different combinations of variants impact a person's chance of getting the disease, also known as a polygenic risk score. Currently, most genome-wide association studies only incorporate genetic data from people with European ancestry. Consequently, polygenic risk scores performed using this information may not accurately predict the risk of developing the disease for individuals with other ethnicities, such as people with Asian ancestry. Here, Ho et al. evaluated how well previously calculated polygenic risk scores for the four most common cancers (breast, colorectal, prostate and lung) worked on individuals of East Asian descent. The scores were tested on a dataset containing the genetic sequence, medical history, diet and activity levels of over 21,000 people living in Singapore in the 1990s. Ho et al. found that the polygenic risk scores for breast, prostate and colorectal cancer were able to predict disease risk. However, the score for lung cancer did not perform as well. The polygenic risk score for breast cancer was the most accurate, and was able to stratify individuals into distinct risk bands at an earlier age than other scores. These findings shed light on which existing polygenic risk scores will be effective at assessing cancer risk in individuals with East Asian ancestry. Indeed, Ho et al. have already incorporated the polygenic risk score for breast cancer into a pilot study screening individuals in a comparable population in Singapore. However, the polygenic risk scores tested still performed better on individuals with European ancestry, highlighting the need to address the lack of Asian representation in genome-wide association studies.

Identification of new risk loci shared across systemic vasculitides points towards potential target genes for drug repurposing.

OBJECTIVES: The number of susceptibility loci currently associated with vasculitis is lower than in other immune-mediated diseases due in part to small cohort sizes, a consequence of the low prevalence of vasculitides. This study aimed to identify new genetic risk loci for the main systemic vasculitides through a comprehensive analysis of their genetic overlap. METHODS: Genome-wide data from 8467 patients with any of the main forms of vasculitis and 29 795 healthy controls were meta-analysed using ASSET. Pleiotropic variants were functionally annotated and linked to their target genes. Prioritised genes were queried in DrugBank to identify potentially repositionable drugs for the treatment of vasculitis. RESULTS: Sixteen variants were independently associated with two or more vasculitides, 15 of them representing new shared risk loci. Two of these pleiotropic signals, located close to CTLA4 and CPLX1, emerged as novel genetic risk loci in vasculitis. Most of these polymorphisms appeared to affect vasculitis by regulating gene expression. In this regard, for some of these common signals, potential causal genes were prioritised based on functional annotation, including CTLA4, RNF145, IL12B, IL5, IRF1, IFNGR1, PTK2B, TRIM35, EGR2 and ETS2, each of which has key roles in inflammation. In addition, drug repositioning analysis showed that several drugs, including abatacept and ustekinumab, could be potentially repurposed in the management of the analysed vasculitides. CONCLUSIONS: We identified new shared risk loci with functional impact in vasculitis and pinpointed potential causal genes, some of which could represent promising targets for the treatment of vasculitis.

Robust SNP-based prediction of rheumatoid arthritis through machine-learning-optimized polygenic risk score.

BACKGROUND: The popular statistics-based Genome-wide association studies (GWAS) have provided deep insights into the field of complex disorder genetics. However, its clinical applicability to predict disease/trait outcomes remains unclear as statistical models are not designed to make predictions. This study employs statistics-free machine-learning (ML)-optimized polygenic risk score (PRS) to complement existing GWAS and bring the prediction of disease/trait outcomes closer to clinical application. Rheumatoid Arthritis (RA) was selected as a model disease to demonstrate the robustness of ML in disease prediction as RA is a prevalent chronic inflammatory joint disease with high mortality rates, affecting adults at the economic prime. Early identification of at-risk individuals may facilitate measures to mitigate the effects of the disease. METHODS: This study employs a robust ML feature selection algorithm to identify single nucleotide polymorphisms (SNPs) that can predict RA from a set of training data comprising RA patients and population control samples. Thereafter, selected SNPs were evaluated for their predictive performances across 3 independent, unseen test datasets. The selected SNPs were subsequently used to generate PRS which was also evaluated for its predictive capacity as a sole feature. RESULTS: Through robust ML feature selection, 9 SNPs were found to be the minimum number of features for excellent predictive performance (AUC > 0.9) in 3 independent, unseen test datasets. PRS based on these 9 SNPs was significantly associated with (P < 1 × 10-16) and predictive (AUC > 0.9) of RA in the 3 unseen datasets. A RA ML-PRS calculator of these 9 SNPs was developed ( https://xistance.shinyapps.io/prs-ra/ ) to facilitate individualized clinical applicability. The majority of the predictive SNPs are protective, reside in non-coding regions, and are either predicted to be potentially functional SNPs (pfSNPs) or in high linkage disequilibrium (r2 > 0.8) with un-interrogated pfSNPs. CONCLUSIONS: These findings highlight the promise of this ML strategy to identify useful genetic features that can robustly predict disease and amenable to translation for clinical application.

MUC5AC Genetic Variation Is Associated With Tuberculous Meningitis Cerebral Spinal Fluid Cytokine Responses and Mortality.

BACKGROUND: The purpose of this study was to assess if single nucleotide polymorphisms (SNPs) in lung mucins MUC5B and MUC5AC are associated with Mycobacterium tuberculosis outcomes. METHODS: Independent SNPs in MUC5B and MUC5AC (genotyped by Illumina HumanOmniExpress array) were assessed for associations with tumor necrosis factor (TNF) concentrations (measured by immunoassay) in cerebral spinal fluid (CSF) from tuberculous meningitis (TBM) patients. SNPs associated with CSF TNF concentrations were carried forward for analyses of pulmonary and meningeal tuberculosis susceptibility and TBM mortality. RESULTS: MUC5AC SNP rs28737416 T allele was associated with lower CSF concentrations of TNF (P = 1.8 × 10-8) and IFN-γ (P = 2.3 × 10-6). In an additive genetic model, rs28737416 T/T genotype was associated with higher susceptibility to TBM (odds ratio [OR], 1.24; 95% confidence interval [CI], 1.03-1.49; P = .02), but not pulmonary tuberculosis (OR, 1.11, 95% CI, .98-1.25; P = .10). TBM mortality was higher among participants with the rs28737416 T/T and T/C genotypes (35/119, 30.4%) versus the C/C genotype (11/89, 12.4%; log-rank P = .005) in a Vietnam discovery cohort (n = 210), an independent Vietnam validation cohort (n = 87; 9/87, 19.1% vs 1/20, 2.5%; log-rank P = .02), and an Indonesia validation cohort (n = 468, 127/287, 44.3% vs 65/181, 35.9%; log-rank P = .06). CONCLUSIONS: MUC5AC variants may contribute to immune changes that influence TBM outcomes.

Increased BMI and late-life mobility dysfunction; overlap of genetic effects in brain regions.

BACKGROUND: How obesity earlier in life impacts upon mobility dysfunctions in late life is not well understood. Pernicious effects of excess weight on the musculoskeletal system and mobility dysfunctions are well-recognized. However, increasingly more data support the link of obesity to overall motor defects that are regulated in the brain. OBJECTIVES: To assess the causal relationship between body mass index (BMI) at midlife and performance of the Timed Up-and-Go test (TUG) in late life among a population-based longitudinal cohort of Chinese adults living in Singapore. METHODS: We evaluated genetic predispositions for BMI in 8342 participants who were followed up from measurement of BMI at average 53 years, to TUG test (as a functional mobility measure) 20 years later. RESULTS: A robust 75.83% of genetically determined BMI effects on late-life TUG scores were mediated through midlife BMI (Pindirect-effect = 9.24 × 10-21). Utilizing Mendelian randomization, we demonstrated a causal effect between BMI and functional mobility in late life (ßIVW = 0.180, PIVW = 0.001). Secondary gene enrichment evaluations highlighted down-regulation of genes at BMI risk loci that were correlated with poorer functional mobility in the substantia nigra and amygdala regions as compared to all other tissues. These genes also exhibit differential expression patterns during human brain development. CONCLUSIONS: We report a causal effect of obesity on mobility dysfunction. Our findings highlight potential neuronal dysfunctions in regulating predispositions on the causal pathway from obesity to mobility dysfunction.

Cross-Ancestry Genome-Wide Association Study Defines the Extended CYP2D6 Locus as the Principal Genetic Determinant of Endoxifen Plasma Concentrations.

The therapeutic efficacy of tamoxifen is predominantly mediated by its active metabolites 4-hydroxy-tamoxifen and endoxifen, whose formation is catalyzed by the polymorphic cytochrome P450 2D6 (CYP2D6). Yet, known CYP2D6 polymorphisms only partially determine metabolite concentrations in vivo. We performed the first cross-ancestry genome-wide association study with well-characterized patients of European, Middle-Eastern, and Asian descent (n = 497) to identify genetic factors impacting active and parent metabolite formation. Genome-wide significant variants were functionally evaluated in an independent liver cohort (n = 149) and in silico. Metabolite prediction models were validated in two independent European breast cancer cohorts (n = 287, n = 189). Within a single 1-megabase (Mb) region of chromosome 22q13 encompassing the CYP2D6 gene, 589 variants were significantly associated with tamoxifen metabolite concentrations, particularly endoxifen and metabolic ratio (MR) endoxifen/N-desmethyltamoxifen (minimal P = 5.4E-35 and 2.5E-65, respectively). Previously suggested other loci were not confirmed. Functional analyses revealed 66% of associated, mostly intergenic variants to be significantly correlated with hepatic CYP2D6 activity or expression (ρ = 0.35 to -0.52), and six hotspot regions in the extended 22q13 locus impacting gene regulatory function. Machine learning models based on hotspot variants (n = 12) plus CYP2D6 activity score (AS) increased the explained variability (~ 9%) compared with AS alone, explaining up to 49% (median R2 ) and 72% of the variability in endoxifen and MR endoxifen/N-desmethyltamoxifen, respectively. Our findings suggest that the extended CYP2D6 locus at 22q13 is the principal genetic determinant of endoxifen plasma concentration. Long-distance haplotypes connecting CYP2D6 with adjacent regulatory sites and nongenetic factors may account for the unexplained portion of variability.

Impact of Genetic Variants in the Nicotine Metabolism Pathway on Nicotine Metabolite Levels in Smokers.

BACKGROUND: Nicotine metabolism is a major factor in nicotine dependence, with approximately 70% to 80% of nicotine metabolized to cotinine in Caucasians. Cotinine formation is catalyzed primarily by CYP2A6, which also converts cotinine to trans-3'-hydroxycotinine (3HC). The goal of the present study was to examine the effects of CYP2A6 deficiency on nicotine metabolism profiles in vivo and the importance of genetic variants in nicotine-metabolizing enzyme genes on urinary nicotine metabolites levels. METHODS: Urine samples from 722 smokers who participated in the Singapore Chinese Health Study were analyzed using UPLC-MS/MS to detect nicotine and eight of its urinary metabolites, and a total of 58 variants in 12 genes involved in nicotine metabolism were investigated in 475 of these subjects with informative genotyping data. RESULTS: Urine samples stratified by the ratio of 3HC/cotinine exhibited a 7-fold increase in nicotine-N'-oxide, a 6-fold increase in nicotine-Glucuronide (Gluc), and a 5-fold decrease in 3HC-Gluc when comparing the lower versus upper 3HC/cotinine ventiles. Significant (P < 0.0001) associations were observed between functional metabolizing enzyme genotypes and levels of various urinary nicotine metabolites, including CYP2A6 genotype and levels of nicotine, nicotine-Gluc, nicotine-N'-oxide and 3HC, UGT2B10 genotype and levels of cotinine, nicotine-Gluc and cotinine-Gluc, UGT2B17 genotype and levels of 3HC-Gluc, FMO3 genotype and levels of nicotine-N'-oxide, and CYP2B6 genotype and levels of nicotine-N'-oxide and 4-hydroxy-4-(3-pyridyl)-butanoic acid. CONCLUSIONS: These data suggest that several pathways are important in nicotine metabolism. IMPACT: Genotype differences in several nicotine-metabolizing enzyme pathways may potentially lead to differences in nicotine dependence and smoking behavior and cessation.

Case report: Expanding the phenotype of ARHGEF17 mutations from increased intracranial aneurysm risk to a neurodevelopmental disease.

RhoGTPase regulators play a key role in the development of the nervous system, and their dysfunction can result in brain malformation and associated disorders. Several guanine nucleotide exchange factors (GEF) have been linked to neurodevelopmental disorders. In line with this, ARHGEF17 has been recently linked as a risk gene to intracranial aneurysms. Here we report siblings of a consanguineous Pakistani family with biallelic variants in the ARHGEF17 gene associated with a neurodevelopmental disorder with intellectual disability, speech delay and motor dysfunction but not aneurysms. Cranial MRI performed in one patient revealed generalized brain atrophy with an enlarged ventricular system, thin corpus callosum and microcephaly. Whole exome sequencing followed by Sanger sequencing in two of the affected individuals revealed a homozygous missense variant (g.11:73021307, c.1624C>T (NM_014786.4), p.R542W) in the ARHGEF17 gene. This variant is in a highly conserved DCLK1 phosphorylation consensus site (I/L/V/F/M]RRXX[pS/pT][I/L/M/V/F) of the protein. Our report expands the phenotypic spectrum of ARHGEF17 variants from increased intracranial aneurysm risk to neurodevelopmental disease and thereby add ARHGEF17 to the list of GEF genes involved in neurodevelopmental disorders.

Biallelic variants in CHST3 cause Spondyloepiphyseal dysplasia with joint dislocations in three Pakistani kindreds.

BACKGROUND: Skeletal dysplasia is a heterogeneous group of disorders. Spondyloepiphyseal dysplasias comprise one subgroup. Deficiency of carbohydrate sulfotransferase 3 has been reported in a small number of patients with recessively inherited spondyloepiphyseal dysplasia with joint dislocation, short stature and scoliosis. We report here molecular and clinical findings of affected individuals in three consanguineous Pakistani families. Affected individuals in all three families had a uniform phenotype including severe short stature, multiple dislocated joints, progressive scoliosis and facial dysmorphism. METHODS: Clinical evaluation was done for three unrelated families. Radiological survey of bones was completed for patients from two of the families. Whole exome sequencing index patients from each family was performed followed by Sanger sequencing for validation of segregation of identified variants in respective families. In-silico analysis for determining pathogenicity of identified variants and conservation was done. RESULTS: Whole-exome sequencing revealed biallelic variants c.590 T > C;p.(Leu197Pro), c.603C > A;p.(Tyr201Ter) and c.661C > T;p.(Arg221Cys) in CHST3 (NM_004273.5) in the three families with eight, five and two affected individuals, respectively. Contrary to previous reports, affected individuals in none of the families exhibited a hearing loss. CONCLUSION: We describe genotypic and phenotypic findings of three unrelated families with spondyloepiphyseal dysplasia. Our study confirms phenotypic variability and adds to the genotypic spectrum of spondyloepiphyseal dysplasia.

Germline Variants Associated with Nasopharyngeal Carcinoma Predisposition Identified through Whole-Exome Sequencing.

The current understanding of genetic susceptibility factors for nasopharyngeal carcinoma (NPC) is still incomplete. To identify novel germline variants associated with NPC predisposition, we analysed whole-exome sequencing data from 119 NPC patients from Singapore with a family history of NPC and/or with early-onset NPC, together with 1337 Singaporean participants without NPC. Variants were prioritised and filtered by selecting variants with minor allele frequencies of <1% in both local control (n = 1337) and gnomAD non-cancer (EAS) (n = 9626) cohorts and a high pathogenicity prediction (CADD score > 20). Using single-variant testing, we identified 17 rare pathogenic variants in 17 genes that were associated with NPC. Consistent evidence of enrichment in NPC patients was observed for five of these variants (in JAK2, PRDM16, LRP1B, NIN, and NKX2-1) from an independent case-control comparison of 156 NPC patients and 9770 unaffected individuals. In a family with five siblings, a FANCE variant (p. P445S) was detected in two affected members, but not in three unaffected members. Gene-based burden testing recapitulated variants in NKX2-1 and FANCE as being associated with NPC risk. Using pathway analysis, endocytosis and immune-modulating pathways were found to be enriched for mutation burden. This study has identified NPC-predisposing variants and genes which could shed new insights into the genetic predisposition of NPC.