Cell-based screen identifies porphyrins as FGFR3 activity inhibitors with therapeutic potential for achondroplasia and cancer.

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 (FGFR3ACH 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.

A Functional Polymorphism Downstream of Vitamin A Regulator Gene CYP26B1 Is Associated with Hand Osteoarthritis.

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-7.31 × 10-6) clustered at 2p13.2 downstream of the CYP26B1 gene. The strongest association signal identified was rs883313 (p-value = 6.76 × 10-7, odds ratio (OR) = 1.76), followed by rs12713768 (p-value = 1.36 × 10-6, 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, r2 = 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.

Modifiable factors affecting renal preservation in type I glycogen storage disease after liver transplantation: a single-center propensity-match cohort study.

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.73m2, 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.

Genetic profiles of 103,106 individuals in the Taiwan Biobank provide insights into the health and history of Han Chinese.

Personalized medical care focuses on prediction of disease risk and response to medications. To build the risk models, access to both large-scale genomic resources and human genetic studies is required. The Taiwan Biobank (TWB) has generated high-coverage, whole-genome sequencing data from 1492 individuals and genome-wide SNP data from 103,106 individuals of Han Chinese ancestry using custom SNP arrays. Principal components analysis of the genotyping data showed that the full range of Han Chinese genetic variation was found in the cohort. The arrays also include thousands of known functional variants, allowing for simultaneous ascertainment of Mendelian disease-causing mutations and variants that affect drug metabolism. We found that 21.2% of the population are mutation carriers of autosomal recessive diseases, 3.1% have mutations in cancer-predisposing genes, and 87.3% carry variants that affect drug response. We highlight how TWB data provide insight into both population history and disease burden, while showing how widespread genetic testing can be used to improve clinical care.

Association of an IGHV3-66 gene variant with Kawasaki disease.

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-9). 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-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.

Genetic risk score constructed from common genetic variants is associated with cardiovascular disease risk in type 2 diabetes mellitus.

BACKGROUND: Patients with type 2 diabetes mellitus (T2DM) experience a two-fold increased risk of cardiovascular diseases. Genome-wide association studies (GWAS) have identified T2DM susceptibility genetic variants. Interestingly, the genetic variants associated with cardiovascular disease risk in T2DM Han Chinese remain to be elucidated. The present study aimed to investigate the genetic variants associated with cardiovascular disease risk in T2DM. METHODS: We performed bootstrapping, GWAS and an investigation of genetic variants associated with cardiovascular disease risk in a discovery T2DM cohort and in a replication cohort. The discovery cohort included 326 cardiovascular disease patients and 1209 noncardiovascular disease patients. The replication cohort included 68 cardiovascular disease patients and 317 noncardiovascular disease patients. The main outcome measures were genetic variants for genetic risk score (GRS) in cardiovascular disease risk in T2DM. RESULTS: In total, 35 genetic variants were associated with cardiovascular disease risk. A GRS was generated by combining risk alleles from these variants weighted by their estimated effect sizes (log odds ratio [OR]). T2DM patients with weighted GRS ≥ 12.63 had an approximately 15-fold increase in cardiovascular disease risk (odds ratio = 15.67, 95% confidence interval [CI] = 10.33-24.00) compared to patients with weighted GRS < 10.39. With the addition of weighted GRS, receiver-operating characteristic curves showed that area under the curve with conventional risk factors was improved from 0.719 (95% CI = 0.689-0.750) to 0.888 (95% CI = 0.866-0.910). CONCLUSIONS: These 35 genetic variants are associated with cardiovascular disease risk in T2DM, alone and cumulatively. T2DM patients with higher levels of weighted genetic risk score have higher cardiovascular disease risks.

Genome-Wide and Candidate Gene Association Analyses Identify a 14-SNP Combination for Hypertension in Patients With Type 2 Diabetes.

BACKGROUND: High blood pressure is common and comorbid with type 2 diabetes (T2D). Almost 50% of patients with T2D have high blood pressure. Patients with both conditions of hypertension (HTN) and T2D are at risk for cardiovascular diseases and mortality. The study aim was to investigate genetic risk factors for HTN in T2D patients. METHODS: This study included 999 T2D (cohort 1) patients for the first genome scan stage and 922 T2D (cohort 2) patients for the replication stage. Here, we investigated the genetic susceptibility and cumulative weighted genetic risk score for HTN in T2D patients of Han Chinese descent in Taiwan. RESULTS: Thirty novel genetic single nucleotide polymorphisms (SNPs) were associated with HTN in T2D after adjusting for age and body mass index (P value <1 × 10-4). Eight blood pressure-related and/or HTN-related genetic SNPs were associated with HTN in T2D after adjusting for age and body mass index (P value <0.05). Linkage disequilibrium and cumulative weighted genetic risk score analyses showed that 14 of the 38 SNPs were associated with risk of HTN in a dose-dependent manner in T2D (Cochran-Armitage trend test: P value <0.0001). The 14-SNP cumulative weighted genetic risk score was also associated with increased regression tendency of systolic blood pressure in T2D (SBP = 122.05 + 0.8 × weighted genetic risk score; P value = 0.0001). CONCLUSIONS: A cumulative weighted genetic risk score composed of 14 SNPs is important for HTN, increased tendency of systolic blood pressure, and may contribute to HTN risk in T2D in Taiwan.

Towards a reference genome that captures global genetic diversity.

The current human reference genome is predominantly derived from a single individual and it does not adequately reflect human genetic diversity. Here, we analyze 338 high-quality human assemblies of genetically divergent human populations to identify missing sequences in the human reference genome with breakpoint resolution. We identify 127,727 recurrent non-reference unique insertions spanning 18,048,877 bp, some of which disrupt exons and known regulatory elements. To improve genome annotations, we linearly integrate these sequences into the chromosomal assemblies and construct a Human Diversity Reference. Leveraging this reference, an average of 402,573 previously unmapped reads can be recovered for a given genome sequenced to ~40X coverage. Transcriptomic diversity among these non-reference sequences can also be directly assessed. We successfully map tens of thousands of previously discarded RNA-Seq reads to this reference and identify transcription evidence in 4781 gene loci, underlining the importance of these non-reference sequences in functional genomics. Our extensive datasets are important advances toward a comprehensive reference representation of global human genetic diversity.

Copy number variant hotspots in Han Taiwanese population induced pluripotent stem cell lines - lessons from establishing the Taiwan human disease iPSC Consortium Bank.

BACKGROUND: The Taiwan Human Disease iPSC Service Consortium was established to accelerate Taiwan's growing stem cell research initiatives and provide a platform for researchers interested in utilizing induced pluripotent stem cell (iPSC) technology. The consortium has generated and characterized 83 iPSC lines: 11 normal and 72 disease iPSC lines covering 21 different diseases, several of which are of high incidence in Taiwan. Whether there are any reprogramming-induced recurrent copy number variant (CNV) hotspots in iPSCs is still largely unknown. METHODS: We performed genome-wide copy number variant screening of 83 Han Taiwanese iPSC lines and compared them with 1093 control subjects using an Affymetrix genome-wide human SNP array. RESULTS: In the iPSCs, we identified ten specific CNV loci and seven "polymorphic" CNV regions that are associated with the reprogramming process. Additionally, we established several differentiation protocols for our iPSC lines. We demonstrated that our iPSC-derived cardiomyocytes respond to pharmacological agents and were successfully engrafted into the mouse myocardium demonstrating their potential application in cell therapy. CONCLUSIONS: The CNV hotspots induced by cell reprogramming have successfully been identified in the current study. This finding may be used as a reference index for evaluating iPSC quality for future clinical applications. Our aim was to establish a national iPSC resource center generating iPSCs, made available to researchers, to benefit the stem cell community in Taiwan and throughout the world.

Genetic Architecture Associated With Familial Short Stature.

CONTEXT: Human height is an inheritable, polygenic trait under complex and multilocus genetic regulation. Familial short stature (FSS; also called genetic short stature) is the most common type of short stature and is insufficiently known. OBJECTIVE: To investigate the FSS genetic profile and develop a polygenic risk predisposition score for FSS risk prediction. DESIGN AND SETTING: The FSS participant group of Han Chinese ancestry was diagnosed by pediatric endocrinologists in Taiwan. PATIENTS AND INTERVENTIONS: The genetic profiles of 1163 participants with FSS were identified by using a bootstrapping subsampling and genome-wide association studies (GWAS) method. MAIN OUTCOME MEASURES: Genetic profile, polygenic risk predisposition score for risk prediction. RESULTS: Ten novel genetic single nucleotide polymorphisms (SNPs) and 9 reported GWAS human height-related SNPs were identified for FSS risk. These 10 novel SNPs served as a polygenic risk predisposition score for FSS risk prediction (area under the curve: 0.940 in the testing group). This FSS polygenic risk predisposition score was also associated with the height reduction regression tendency in the general population. CONCLUSION: A polygenic risk predisposition score composed of 10 genetic SNPs is useful for FSS risk prediction and the height reduction tendency. Thus, it might contribute to FSS risk in the Han Chinese population from Taiwan.

Both IgM and IgG Antibodies Against Polyethylene Glycol Can Alter the Biological Activity of Methoxy Polyethylene Glycol-Epoetin Beta in Mice.

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 125I-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.

Genomic interrogation of familial short stature contributes to the discovery of the pathophysiological mechanisms and pharmaceutical drug repositioning.

BACKGROUND: Genetic factors, dysregulation in the endocrine system, cytokine and paracrine factors are implicated in the pathogenesis of familial short stature (FSS). Nowadays, the treatment choice for FSS is limited, with only recombinant human growth hormone (rhGH) being available. METHODS: Herein, starting from the identification of 122 genetic loci related to FSS, we adopted a genetic-driven drug discovery bioinformatics pipeline based on functional annotation to prioritize crucial biological FSS-related genes. These genes were suggested to be potential targets for therapeutics. RESULTS: We discovered five druggable subnetworks, which contained seven FSS-related genes and 17 druggable targerts. CONCLUSIONS: This study provides a valuable drug repositioning accompanied by corresponding targetable gene clusters for FSS therapy.

Genome-wide transcriptome analysis to further understand neutrophil activation and lncRNA transcript profiles in Kawasaki disease.

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 CD177pos 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.

Interethnic analyses of blood pressure loci in populations of East Asian and European descent.

Blood pressure (BP) is a major risk factor for cardiovascular disease and more than 200 genetic loci associated with BP are known. Here, we perform a multi-stage genome-wide association study for BP (max N = 289,038) principally in East Asians and meta-analysis in East Asians and Europeans. We report 19 new genetic loci and ancestry-specific BP variants, conforming to a common ancestry-specific variant association model. At 10 unique loci, distinct non-rare ancestry-specific variants colocalize within the same linkage disequilibrium block despite the significantly discordant effects for the proxy shared variants between the ethnic groups. The genome-wide transethnic correlation of causal-variant effect-sizes is 0.898 and 0.851 for systolic and diastolic BP, respectively. Some of the ancestry-specific association signals are also influenced by a selective sweep. Our results provide new evidence for the role of common ancestry-specific variants and natural selection in ethnic differences in complex traits such as BP.

Immunoglobulin profiling identifies unique signatures in patients with Kawasaki disease during intravenous immunoglobulin treatment.

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.

Loss of GPNMB Causes Autosomal-Recessive Amyloidosis Cutis Dyschromica in Humans.

Amyloidosis cutis dyschromica (ACD) is a distinct form of primary cutaneous amyloidosis characterized by generalized hyperpigmentation mottled with small hypopigmented macules on the trunks and limbs. Affected families and sporadic case subjects have been reported predominantly in East and Southeast Asian ethnicities; however, the genetic cause has not been elucidated. We report here that the compound heterozygosity or homozygosity of GPNMB truncating alleles is the cause of autosomal-recessive ACD. Six nonsense or frameshift mutations were identified in nine individuals diagnosed with ACD. Immunofluorescence analysis of skin biopsies showed that GPNMB is expressed in all epidermal cells, with the highest staining observed in melanocytes. GPNMB staining is significantly reduced in the lesional skin of affected individuals. Hyperpigmented lesions exhibited significantly increased amounts of DNA/keratin-positive amyloid deposits in the papillary dermis and infiltrating macrophages compared with hypo- or depigmented macules. Depigmentation of the lesions was attributable to loss of melanocytes. Intracytoplasmic fibrillary aggregates were observed in keratinocytes scattered in the lesional epidermis. Thus, our analysis indicates that loss of GPNMB, which has been implicated in melanosome formation, autophagy, phagocytosis, tissue repair, and negative regulation of inflammation, underlies autosomal-recessive ACD and provides insights into the etiology of amyloidosis and pigment dyschromia.

A genome-wide association study links small-vessel ischemic stroke to autophagy.

Genome-wide association studies (GWAS) can serve as strong evidence in correlating biological pathways with human diseases. Although ischemic stroke has been found to be associated with many biological pathways, the genetic mechanism of ischemic stroke is still unclear. Here, we performed GWAS for a major subtype of stroke-small-vessel occlusion (SVO)-to identify potential genetic factors contributing to ischemic stroke. GWAS were conducted on 342 individuals with SVO stroke and 1,731 controls from a Han Chinese population residing in Taiwan. The study was replicated in an independent Han Chinese population comprising an additional 188 SVO stroke cases and 1,265 controls. Three SNPs (rs2594966, rs2594973, rs4684776) clustered at 3p25.3 in ATG7 (encoding Autophagy Related 7), with P values between 2.52 × 10-6 and 3.59 × 10-6, were identified. Imputation analysis also supported the association between ATG7 and SVO stroke. To our knowledge, this is the first GWAS to link stroke and autophagy. ATG7, which has been implicated in autophagy, could provide novel insights into the genetic basis of ischemic stroke.

Identification of functional single nucleotide polymorphisms in the branchpoint site.

BACKGROUND: The human genome contains millions of single nucleotide polymorphisms (SNPs); many of these SNPs are intronic and have unknown functional significance. SNPs occurring within intron branchpoint sites, especially at the adenine (A), would presumably affect splicing; however, this has not been systematically studied. We employed a splicing prediction tool to identify human intron branchpoint sites and screened dbSNP for identifying SNPs located in the predicted sites to generate a genome-wide branchpoint site SNP database. RESULTS: We identified 600 SNPs located within branchpoint sites; among which, 216 showed a change in A. After scoring the SNPs by counting the As in the ± 10 nucleotide region, only four SNPs were identified without additional As (rs13296170, rs12769205, rs75434223, and rs67785924). Using minigene constructs, we examined the effects of these SNPs on splicing. The three SNPs (rs13296170, rs12769205, and rs75434223) with nucleotide substitution at the A position resulted in abnormal splicing (exon skipping and/or intron inclusion). However, rs67785924, a 5-bp deletion that abolished the branchpoint A nucleotide, exhibited normal RNA splicing pattern, presumably using two of the downstream As as alternative branchpoints. The influence of additional As on splicing was further confirmed by studying rs2733532, which contains three additional As in the ± 10 nucleotide region. CONCLUSIONS: We generated a high-confidence genome-wide branchpoint site SNP database, experimentally verified the importance of A in the branchpoint, and suggested that other nearby As can protect branchpoint A substitution from abnormal splicing.

A genome-wide association study identifies a novel susceptibility locus for the immunogenicity of polyethylene glycol.

Conjugation of polyethylene glycol (PEG) to therapeutic molecules can improve bioavailability and therapeutic efficacy. However, some healthy individuals have pre-existing anti-PEG antibodies and certain patients develop anti-PEG antibody during treatment with PEGylated medicines, suggesting that genetics might play a role in PEG immunogenicity. Here we perform genome-wide association studies for anti-PEG IgM and IgG responses in Han Chinese with 177 and 140 individuals, defined as positive for anti-PEG IgM and IgG responses, respectively, and with 492 subjects without either anti-PEG IgM or IgG as controls. We validate the association results in the replication cohort, consisting of 211 and 192 subjects with anti-PEG IgM and anti-PEG IgG, respectively, and 596 controls. We identify the immunoglobulin heavy chain (IGH) locus to be associated with anti-PEG IgM response at genome-wide significance (P = 2.23 × 10-22). Our findings may provide novel genetic markers for predicting the immunogenicity of PEG and efficacy of PEGylated therapeutics.Some individuals develop antibodies against the polyethylene glycol that is commonly used in therapeutic preparations. Here the authors conduct a GWAS in Han Chinese and find the IGH locus is associated with anti-PEG IgM.

Association of human height-related genetic variants with familial short stature in Han Chinese in Taiwan.

Human height can be described as a classical and inherited trait model. Genome-wide association studies (GWAS) have revealed susceptible loci and provided insights into the polygenic nature of human height. Familial short stature (FSS) represents a suitable trait for investigating short stature genetics because disease associations with short stature have been ruled out in this case. In addition, FSS is caused only by genetically inherited factors. In this study, we explored the correlations of FSS risk with the genetic loci associated with human height in previous GWAS, alone and cumulatively. We systematically evaluated 34 known human height single nucleotide polymorphisms (SNPs) in relation to FSS in the additive model (p < 0.00005). A cumulative effect was observed: the odds ratios gradually increased with increasing genetic risk score quartiles (p < 0.001; Cochran-Armitage trend test). Six affected genes-ZBTB38, ZNF638, LCORL, CABLES1, CDK10, and TSEN15-are located in the nucleus and have been implicated in embryonic, organismal, and tissue development. In conclusion, our study suggests that 13 human height GWAS-identified SNPs are associated with FSS risk both alone and cumulatively.