Jer-Yuarn Wu

Han Chinese people comprise nearly 20% of the global population but remain under-represented in genetic studies ^1,2 , so there is an urgent need for large-scale cohorts to advance precision medicine. Here we present the Taiwan Precision Medicine Initiative (TPMI), established by Academia Sinica in collaboration with 16 major medical centres around Taiwan, which has recruited 565,390 participants who consent to provide DNA samples for genetic profiling and grant access to their electronic medical records (EMRs) for research. EMR access is both retrospective and prospective, allowing longitudinal studies. Genetic profiling is done with population-optimized arrays of single-nucleotide polymorphisms for people of Han Chinese ancestry, which enable genome-wide association ^3,4 , phenome-wide association ^5,6 and polygenic risk score ^7,8 studies to be performed to evaluate common disease risk and pharmacogenetic response. Participants also agreed to be re-contacted for future research and receive personalized genetic risk profiles with health management recommendations. The TPMI has established the TPMI Data Access Platform, a central database and analysis platform that both safeguards the security of the data and facilitates academic research. As a large cohort of individuals with non-European ancestry that merges genetic profiles with EMR data and enables longitudinal follow-up, TPMI provides a unique resource that could be used to validate genetic risk prediction models, perform clinical trials of risk-based health management and inform health policies. Ultimately, the TPMI cohort will contribute to global genetic research and serve as a model for population-based precision medicine.

The Taiwan Precision Medicine Initiative (TPMI), a project initiated by the Academia Sinica in collaboration with 16 major medical centers around Taiwan, has recruited 565,390 participants who consented to provide DNA samples for genetic profiling and grant access to their electronic medical records (EMR) for studies to develop precision medicine. Access to the EMR is both retrospective and prospective, allowing researchers to conduct prospective studies over time. Genetic profiling is done with population-optimized SNP arrays for the Han Chinese populations that enable genetic analyses such as genome-wide association, phenome-wide association, and polygenic risk score studies to evaluate common disease risk and pharmacogenetic response. Furthermore, the TPMI participants agree to be contacted for future research opportunities related to their genetic risks and receive personalized genetic risk profiles with health management recommendations. TPMI has established the TPMI Data Access Platform (TDAP), a central database and analysis platform that both safeguards the security of the data and facilitates academic research. The TPMI is the largest non-European cohort that merges genetic profiles with EMR in the world. With a cohort that can be followed over time, it can be utilized to validate genetic risk prediction models, conduct clinical trials to show the efficacy of risk-based health management, and optimize health policies based on genetic risks. In this report, we describe the TPMI study design, the population and genetic characteristics of the TPMI cohort, and the power it provides to conduct crucial studies in developing precision medicine on a population and personal level. As Han Chinese represent almost 20% of the world’s population, the results of TPMI studies will benefit >1.4 billion people around the world and serve as a model for developing population-based precision medicine.

Lumbar spinal stenosis (LSS) is a degenerative spinal condition characterized by the narrowing of the spinal canal, resulting in low back pain (LBP) and limited leg mobility. Twin and family studies have suggested that genetics contributes to disease progression. However, the genetic causes of familial LSS remain unclear. We performed whole-exome and direct sequencing on seven female patients from a Han Chinese family with LBP, among whom four developed LSS. Based on our genetic findings, we performed gene knockdown studies in human chondrocytes to study possible pathological mechanisms underlying LSS. We found a novel nonsense mutation, c.417C > G (NM_002183, p.Y139X), in IL3RA, shared by all the LBP/LSS cases. Knockdown of IL3RA led to a reduction in the total collagen content of 81.6% in female chondrocytes and 21% in male chondrocytes. The expression of MMP-1, -3, and/or -10 significantly increased, with a more pronounced effect observed in females than in males. Furthermore, EsRb expression significantly decreased following IL3RA knockdown. Moreover, the knockdown of EsRb resulted in increased MMP-1 and -10 expression in chondrocytes from females. We speculate that IL3RA deficiency could lead to a reduction in collagen content and intervertebral disk (IVD) strength, particularly in females, thereby accelerating IVD degeneration and promoting LSS occurrence. Our results illustrate, for the first time, the association between IL3RA and estrogen receptor beta, highlighting their importance and impact on MMPs and collagen in degenerative spines in women.

Overactive fibroblast growth factor receptor 3 (FGFR3) signaling drives pathogenesis in a variety of cancers and a spectrum of short-limbed bone dysplasias, including the most common form of human dwarfism, achondroplasia (ACH). Targeting FGFR3 activity holds great promise as a therapeutic approach for treatment of these diseases. Here, we established a receptor/adaptor translocation assay system that can specifically monitor FGFR3 activation, and we applied it to identify FGFR3 modulators from complex natural mixtures. An FGFR3-suppressing plant extract of Amaranthus viridis was identified from the screen, and 2 bioactive porphyrins, pheophorbide a (Pa) and pyropheophorbide a, were sequentially isolated from the extract and functionally characterized. Further analysis showed that Pa reduced excessive FGFR3 signaling by decreasing its half-life in FGFR3-overactivated multiple myeloma cells and chondrocytes. In an ex vivo culture system, Pa alleviated defective long bone growth in humanized ACH mice ( FGFR3 ^ACH mice). Overall, our study presents an approach to discovery and validation of plant extracts or drug candidates that target FGFR3 activation. The compounds identified by this approach may have applications as therapeutics for FGFR3-associated cancers and skeletal dysplasias.

While genetic analyses have revealed ~100 risk loci associated with osteoarthritis (OA), only eight have been linked to hand OA. Besides, these studies were performed in predominantly European and Caucasian ancestries. Here, we conducted a genome-wide association study in the Han Chinese population to identify genetic variations associated with the disease. We recruited a total of 1136 individuals (n = 420 hand OA-affected; n = 716 unaffected control subjects) of Han Chinese ancestry. We carried out genotyping using Axiom Asia Precisi on Medicine Research Array, and we employed the RegulomeDB database and RoadMap DNase I Hypersensitivity Sites annotations to further narrow down our potential candidate variants. Genetic variants identified were tested in the Geisinger’s hand OA cohort selected from the Geisinger MyCode community health initiative (MyCode^®). We also performed a luciferase reporter assay to confirm the potential impact of top candidate single-nucleotide polymorphisms (SNPs) on hand OA. We identified six associated SNPs (p-value = 6.76 × 10⁻⁷–7.31 × 10⁻⁶) clustered at 2p13.2 downstream of the CYP26B1 gene. The strongest association signal identified was rs883313 (p-value = 6.76 × 10⁻⁷, odds ratio (OR) = 1.76), followed by rs12713768 (p-value = 1.36 × 10⁻⁶, OR = 1.74), near or within the enhancer region closest to the CYP26B1 gene. Our findings showed that the major risk-conferring CC haplotype of SNPs rs12713768 and rs10208040 [strong linkage disequilibrium (LD); D’ = 1, r² = 0.651] drives 18.9% of enhancer expression activity. Our findings highlight that the SNP rs12713768 is associated with susceptibility to and severity of hand OA in the Han Chinese population and that the suggested retinoic acid signaling pathway may play an important role in its pathogenesis.

Background and aims: Glycogen storage disease type I (GSD-I) is an autosomal recessive disorder of carbohydrate metabolism, resulting in limited production of glucose and excessive glycogen storage in the liver and kidneys. These patients are characterized by life-threatening hypoglycemia, metabolic derangements, hepatomegaly, chronic kidney disease, and failure to thrive. Liver transplantation (LT) has been performed for poor metabolic control and delayed growth. However, renal outcome was diverse in pediatric GSD patients after LT. The aim of this study was to investigate the long-term outcome of renal function in pediatric GSD-I patients after living donor LT (LDLT), and to identify modifiable variables that potentially permits LT to confer native renal preservation. Methods: The study included eight GSD-Ia and one GSD-Ib children with a median age of 9.0 (range 4.2–15.7) years at the time of LT. Using propensity score matching, 20 children with biliary atresia (BA) receiving LT were selected as the control group by matching for age, sex, pre-operative serum creatinine (SCr) and pediatric end-stage liver disease (PELD) score. Renal function was evaluated based on the SCr, estimated glomerular filtration rate (eGFR), microalbuminuria, and morphological changes in the kidneys. Comparability in long-term renal outcome in terms of anatomic and functional parameters will help to identify pre-LT factors of GSD-I that affect renal prognosis. Results: The clinical and biochemical characteristics of the GSD and BA groups were similar, including immunosuppressive regimens and duration of follow-up (median 15 years) after LT. Overall, renal function, including eGFR and microalbuminuria was comparable in the GSD-I and BA groups (median eGFR: 111 vs. 123 ml/min/1.73m², P = 0.268; median urine microalbuminuria to creatinine ratio: 16.0 vs. 7.2 mg/g, P = 0.099, respectively) after LT. However, in the subgroups of the GSD cohort, patients starting cornstarch therapy at an older age (≥ 6-year-old) before transplantation demonstrated a worse renal outcome in terms of eGFR change over years (P < 0.001). In addition, the enlarged kidney in GSD-I returned to within normal range after LT. Conclusions: Post-LT renal function was well-preserved in most GSD-I patients. Early initiation of cornstarch therapy before preschool age, followed by LT, achieved a good renal prognosis.

In a meta-analysis of three GWAS for susceptibility to Kawasaki disease (KD) conducted in Japan, Korea, and Taiwan and follow-up studies with a total of 11,265 subjects (3428 cases and 7837 controls), a significantly associated SNV in the immunoglobulin heavy variable gene (IGHV) cluster in 14q33.32 was identified (rs4774175; OR = 1.20, P = 6.0 × 10⁻⁹). Investigation of nonsynonymous SNVs of the IGHV cluster in 9335 Japanese subjects identified the C allele of rs6423677, located in IGHV3-66, as the most significant reproducible association (OR = 1.25, P = 6.8 × 10⁻¹⁰ in 3603 cases and 5731 controls). We observed highly skewed allelic usage of IGHV3-66, wherein the rs6423677 A allele was nearly abolished in the transcripts in peripheral blood mononuclear cells of both KD patients and healthy adults. Association of the high-expression allele with KD strongly indicates some active roles of B-cells or endogenous immunoglobulins in the disease pathogenesis. Considering that significant association of SNVs in the IGHV region with disease susceptibility was previously known only for rheumatic heart disease (RHD), a complication of acute rheumatic fever (ARF), these observations suggest that common B-cell related mechanisms may mediate the symptomology of KD and ARF as well as RHD.

Pre-existing antibodies that bind polyethylene glycol are present in about 40% of healthy individuals. It is currently unknown if pre-existing anti-polyethylene glycol (PEG) antibodies can alter the bioactivity of pegylated drugs with a single long PEG chain, which represents the majority of newly developed pegylated medicines. Methoxy polyethylene glycol-epoetin beta (PEG-EPO) contains a single 30 kDa PEG chain and is used to treat patients suffering from anemia. We find that the pre-existing human anti-PEG IgM and IgG antibodies from normal donors can bind to PEG-EPO. The prevalence and concentrations of anti-PEG IgM and IgG antibodies were also higher in patients that responded poorly to PEG-EPO. Monoclonal anti-PEG IgM and IgG antibodies at concentrations found in normal donors blocked the biological activity of PEG-EPO to stimulate the production of new erythrocytes in mice and accelerated the clearance of ¹²⁵I-PEG-EPO, resulting in PEG-EPO accumulation primarily in the liver and spleen. Accelerated clearance by the anti-PEG IgG antibody was mediated by the Fc portion of the antibody. Importantly, infusing higher doses of PEG-EPO could compensate for the inhibitory effects of anti-PEG antibodies, suggesting that pre-existing anti-PEG antibodies can be “dosed through.” Our study indicates that the bioactivity and therapeutic activity of PEG-EPO may be reduced in patients with elevated levels of pre-existing anti-PEG antibodies. New pegylated medicines with a single long PEG chain may also be affected in patients with high levels of anti-PEG antibodies.

Kawasaki disease (KD) is the most common cause of acquired cardiac disease in children in developed countries. However, little is known regarding the role of transcriptomic targets of KD in the disease progression and development of complications, especially coronary artery aneurysms (CAA). The aim of our study was to identify transcripts affected by KD and their potential role in the disease. We enrolled 37 KD patients and collected blood samples along a comprehensive time-course. mRNA profiling demonstrated an abundance of CD177 transcript in acute KD, and in the intravenous immunoglobulin (IVIG)-resistant group compared to in the IVIG-sensitive group. lncRNA profiling identified XLOC_006277 as the most highly expressed molecule. XLOC_006277 expression in patients at acute stage was 3.3-fold higher relative to patients with convalescent KD. Moreover, XLOC_006277 abundance increased significantly in patients with CAA. XLOC_006277 knockdown suppressed MMP-8 and MMP-9 expression, both associated with heart lesions. Our result suggested that the increase of CD177^pos neutrophils was associated with KD. Moreover, this study provided global long non-coding RNA transcripts in the blood of patients with KD, IVIG-resistant KD, or CAA. Notably, XLOC_006277 abundance was associated with CAA, which might contribute to further understanding of CAA pathogenesis in KD.

Identifying the causes of high fever syndromes such as Kawasaki disease (KD) remains challenging. To investigate pathogen exposure signatures in suspected pathogen-mediated diseases such as KD, we performed immunoglobulin (Ig) profiling using a next-generation sequencing method. After intravenous Ig (IVIG) treatment, we observed disappearance of clonally expanded IgM clonotypes, which were dominantly observed in acute-phase patients. The complementary-determining region 3 (CDR3) sequences of dominant IgM clonotypes in acute-phase patients were commonly observed in other Ig isotypes. In acute-phase KD patients, we identified 32 unique IgM CDR3 clonotypes shared in three or more cases. Furthermore, before the IVIG treatment, the sums of dominant IgM clonotypes in IVIG-resistant KD patients were significantly higher than those of IVIG-sensitive KD patients. Collectively, we demonstrate a novel approach for identifying certain Ig clonotypes for potentially interacting with pathogens involved in KD; this approach could be applied for a wide variety of fever-causing diseases of unknown origin.