Replication of the results of genome-wide and candidate gene association studies on telomere length in a Korean population

rs2736098, a synonymous polymorphism, is associated with carcinogenesis and cell count in multiple tissue types by regulating TERT expression

rs2736098 is a synonymous polymorphism in TERT (telomerase reverse transcriptase), an enzyme involved in tumor onset of multiple tissues, and should play no roles in carcinogenesis. However, a search in cancer somatic mutation database indicated that the mutation frequency at rs2736098 is much higher than the average one for TERT. Moreover, there are significant H3K4me1 and H3K27Ac signals, two universal histone modifications for active enhancers, surrounding rs2736098. Therefore, we hypothesized that rs2736098 might be within an enhancer region, regulate TERT expression and influence cancer risk. Through luciferase assay, it was verified that the enhancer activity of rs2736098C allele is significantly higher than that of T in multiple tissues. Transfection of plasmids containing TERT coding region with two different alleles indicated that rs2736098C allele can induce a significantly higher TERT expression than T. By chromatin immunoprecipitation, it was observed that the fragment spanning rs2736098 can interact with USF1 (upstream transcription factor 1). The two alleles of rs2736098 present evidently different binding affinity with nuclear proteins. Database and literature search indicated that rs2736098 is significantly associated with carcinogenesis in multiple tissues and count of multiple cell types. All these facts indicated that rs2736098 is also an oncogenic polymorphism and plays important role in cell proliferation.

Comprehensive characterization of coding and non-coding single nucleotide polymorphisms of the Myoneurin (MYNN) gene using molecular dynamics simulation and docking approaches

Genome-wide association studies (GWAS) identified a coding single nucleotide polymorphism, MYNN rs10936599, at chromosome 3q. MYNN gene encodes myoneurin protein, which has been associated with several cancer pathogenesis and disease development processes. However, there needed to be a more detailed characterization of this polymorphism's (and other coding and non-coding polymorphisms) structural, functional, and molecular impact. The current study addressed this gap and analyzed different properties of rs10936599 and non-coding SNPs of MYNN via a thorough computational method. The variant, rs10936599, was predicted functionally deleterious by nine functionality prediction approaches, like SIFT, PolyPhen-2, and REVEL, etc. Following that, structural modifications were estimated through the HOPE server and Mutation3D. Moreover, the mutation was found in a conserved and active residue, according to ConSurf and CPORT. Further, the secondary structures were predicted, followed by tertiary structures, and there was a significant deviation between the native and variant models. Similarly, molecular simulation also showed considerable differences in the dynamic pattern of the wildtype and mutant structures. Molecular docking revealed that the variant binds with better docking scores with ligand NOTCH2. In addition to that, non-coding SNPs located at the MYNN locus were retrieved from the ENSEMBL database. These were found to disrupt the transcription factor binding regulatory regions; nonetheless, only two affect miRNA target sites. Again, eight non-coding variants were detected in the testes with normalized expression, whereas HaploReg v4.1 unveiled annotations for non-coding variants. In summary, in silico comprehensive characterization of coding and non-coding single nucleotide polymorphisms of MYNN gene will assist researchers to work on MYNN gene and establish their association with certain types of cancers.

GWAS-identified telomere length associated genetic variants predict risk of recurrence of HPV-positive oropharyngeal cancer after definitive radiotherapy

BACKGROUND

Lymphocyte telomere length (LTL)-related genetic variants may modulate LTL and affect recurrence of squamous cell carcinoma of the oropharynx (SCCOP).

METHODS

A total of 1013 patients with incident SCCOP were recruited and genotyped for 16 genome-wide association study (GWAS)-identified TL-related polymorphisms. Of these patients, 489 had tumour HPV16 status determination. Univariate and multivariate analyses were performed to evaluate associations.

FINDINGS

Of the 16 TL-related polymorphisms, four were significantly associated with LTL: rs1920116, rs3027234, rs6772228, and rs11125529, and the patients with putatively favourable genotypes had approximately 1.5-3 times the likelihood of shorter LTL compared with patients with the corresponding risk genotypes. Moreover, patients with one to four favourable genotypes of the four combined polymorphisms had approximately 3-11 times the likelihood of shorter LTL compared with patients with no favourable genotype. The four LTL-related polymorphisms were significantly associated with approximately 40% reduced risk (for favourable genotypes) or doubled risk (for risk genotypes) of recurrence, and similar but more pronounced associations were observed in patients with tumour HPV16-positive SCCOP. Similarly, patients with one to four risk genotypes had significantly approximately 2.5-4 times increased recurrence risk compared with patients with no risk genotype, and similar but more pronounced associations were observed in patients with tumour HPV16-positive SCCOP.

INTERPRETATION

Four LTL-related polymorphisms individually or jointly modify LTL and risk of recurrence of SCCOP, particularly HPV-positive SCCOP. These LTL-related polymorphisms could have potential to further stratify patients with HPV-positive SCCOP for individualized treatment and better survival.

FUNDING

Not applicable.

Differential Methylation of Telomere-Related Genes Is Associated with Kidney Disease in Individuals with Type 1 Diabetes

Diabetic kidney disease (DKD) represents a major global health problem. Accelerated ageing is a key feature of DKD and, therefore, characteristics of accelerated ageing may provide useful biomarkers or therapeutic targets. Harnessing multi-omics, features affecting telomere biology and any associated methylome dysregulation in DKD were explored. Genotype data for nuclear genome polymorphisms in telomere-related genes were extracted from genome-wide case-control association data (n = 823 DKD/903 controls; n = 247 end-stage kidney disease (ESKD)/1479 controls). Telomere length was established using quantitative polymerase chain reaction. Quantitative methylation values for 1091 CpG sites in telomere-related genes were extracted from epigenome-wide case-control association data (n = 150 DKD/100 controls). Telomere length was significantly shorter in older age groups (p = 7.6 × 10^-6). Telomere length was also significantly reduced (p = 6.6 × 10^-5) in DKD versus control individuals, with significance remaining after covariate adjustment (p = 0.028). DKD and ESKD were nominally associated with telomere-related genetic variation, with Mendelian randomisation highlighting no significant association between genetically predicted telomere length and kidney disease. A total of 496 CpG sites in 212 genes reached epigenome-wide significance (p ≤ 10^-8) for DKD association, and 412 CpG sites in 193 genes for ESKD. Functional prediction revealed differentially methylated genes were enriched for Wnt signalling involvement. Harnessing previously published RNA-sequencing datasets, potential targets where epigenetic dysregulation may result in altered gene expression were revealed, useful as potential diagnostic and therapeutic targets for intervention.

Harnessing Genomic Analysis to Explore the Role of Telomeres in the Pathogenesis and Progression of Diabetic Kidney Disease

The prevalence of diabetes is increasing globally, and this trend is predicted to continue for future decades. Research is needed to uncover new ways to manage diabetes and its co-morbidities. A significant secondary complication of diabetes is kidney disease, which can ultimately result in the need for renal replacement therapy, via dialysis or transplantation. Diabetic kidney disease presents a substantial burden to patients, their families and global healthcare services. This review highlights studies that have harnessed genomic, epigenomic and functional prediction tools to uncover novel genes and pathways associated with DKD that are useful for the identification of therapeutic targets or novel biomarkers for risk stratification. Telomere length regulation is a specific pathway gaining attention recently because of its association with DKD. Researchers are employing both observational and genetics-based studies to identify telomere-related genes associated with kidney function decline in diabetes. Studies have also uncovered novel functions for telomere-related genes beyond the immediate regulation of telomere length, such as transcriptional regulation and inflammation. This review summarises studies that have revealed the potential to harness therapeutics that modulate telomere length, or the associated epigenetic modifications, for the treatment of DKD, to potentially slow renal function decline and reduce the global burden of this disease.

Genetic susceptibility to telomere shortening through the rs2293607 polymorphism is associated with a greater risk of alcohol use disorder

Telomere shortening is usually considered a biomarker of ageing. Harmful alcohol use promotes accelerated biological ageing and alcohol use disorders (AUDs) are associated with short telomere length (TL). This study was conducted to examine the relationship of TL to AUD and determine whether single nucleotide polymorphisms (SNPs) in TERC and TERT modulate this association. For this purpose, we genotyped TERC SNPs rs2293607, rs12696304, and rs16847897 and TERT SNPs rs2735940, rs2736100, and rs2736098 in 308 male patients with AUD and 255 sex-matched healthy controls and measured TL in a subset of 99 patients and 99 controls paired by age and smoking status. Our results showed that the mean TL was shorter in patients with AUD than in controls. The area under the ROC curve was 0.70 (P < 0.001). The GG genotype of TERC rs2293607 was more common among patients with AUD than among controls (9.8% vs. 5.1%; P = 0.038). No difference was found for the other SNPs. Carriers of the GG genotype of rs2293607 had shorter telomeres than did allele A carriers. In conclusion, patients with AUD had shorter telomeres. Genetic susceptibility to telomere shortening through the rs2293607 SNP is associated with a greater risk of AUD.

Myoneurin regulates BMP signaling by competing with Ppm1a for Smad binding

A delicate balance of BMP activity is critical for tissue formation and organogenesis. However, the mechanical molecular details in ensuring the proper duration and intensity of BMP signaling have yet to be fully elucidated. Here, we identified a zebrafish mutant with a disrupted gene encoding for the BTB/POZ and zinc finger protein myoneurin (Mynn). mynn^−/− mutants exhibited severe loss of pharyngeal cartilage elements, owing to poor proliferation, blocked differentiation, and low viability of cranial neural crest cells. Depletion of mynn in both zebrafish embryos and mammalian cells led to a reduction of the BMP signal activity. Mechanistically, Mynn interacts with Smad proteins in the nucleus, thereby disrupting the association between Smad protein and the phosphatase Ppm1a. Ultimately, this interaction prevents Smad dephosphorylation. More broadly, our findings may provide a new strategy to balance BMP signal activity via competitive binding of Mynn and Ppm1a to Smad proteins during pharyngeal cartilage formation.

•

mynn gene is essential for pharyngeal cartilage development

•

mynn is required for the proliferation, differentiation, and survival of the CNCCs

•

Mynn has an evolutionarily conserved function in supporting BMP signal

•

Mynn maintains BMP signal activity by competing with Ppm1a for Smad binding

The association of telomere maintenance and TERT expression with susceptibility to human papillomavirus infection in cervical epithelium

The role of telomerase reverse transcriptase (TERT) induction and telomere maintenance in carcinogenesis including cervical cancer (CC) pathogenesis has been well established. However, it remains unclear whether they affect infection of high-risk human papillomavirus (hrHPV), an initiating event for CC development. Similarly, genetic variants at the TERT locus are shown to be associated with susceptibility to CC, but it is unclear whether these SNPs modify the risk for cervical HPV infection. Here we show that in CC-derived HeLa cells, TERT overexpression inhibits, while its depletion upregulates expression of Syndecan-1 (SDC-1), a key component for HPV entry receptors. The TCGA cohort of CC analyses reveals an inverse correlation between TERT and SDC-1 expression (R = -0.23, P = 0.001). We further recruited 1330 females (520 non-HPV and 810 hrHPV-infected) without CC or high-grade cervical intraepithelial neoplasia to analyze telomeres in cervical epithelial cells and SNPs at rs2736098, rs2736100 and rs2736108, previously identified TERT SNPs for CC risk. Non-infected females exhibited age-related telomere shortening in cervical epithelial cells and their telomeres were significantly longer than those in hrHPV-infected group (1.31 ± 0.62 vs 1.19 ± 0.48, P < 0.001). There were no differences in rs2736098 and rs2736100 genotypes, but non-infected individuals had significantly a higher C-allele frequency (associated with higher TERT expression) while lower T-allele levels at rs2736108 compared with those in the hrHPV group (P = 0.020). Collectively, appropriate telomere maintenance and TERT expression in normal cervical cells may prevent CC by modulating hrHPV infection predisposition, although they are required for CC development and progression.

Large-scale cross-cancer fine-mapping of the 5p15.33 region reveals multiple independent signals

Genome-wide association studies (GWASs) have identified thousands of cancer risk loci revealing many risk regions shared across multiple cancers. Characterizing the cross-cancer shared genetic basis can increase our understanding of global mechanisms of cancer development. In this study, we collected GWAS summary statistics based on up to 375,468 cancer cases and 530,521 controls for fourteen types of cancer, including breast (overall, estrogen receptor [ER]-positive, and ER-negative), colorectal, endometrial, esophageal, glioma, head/neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancer, to characterize the shared genetic basis of cancer risk. We identified thirteen pairs of cancers with statistically significant local genetic correlations across eight distinct genomic regions. Specifically, the 5p15.33 region, harboring the TERT and CLPTM1L genes, showed statistically significant local genetic correlations for multiple cancer pairs. We conducted a cross-cancer fine-mapping of the 5p15.33 region based on eight cancers that showed genome-wide significant associations in this region (ER-negative breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, and prostate cancer). We used an iterative analysis pipeline implementing a subset-based meta-analysis approach based on cancer-specific conditional analyses and identified ten independent cross-cancer associations within this region. For each signal, we conducted cross-cancer fine-mapping to prioritize the most plausible causal variants. Our findings provide a more in-depth understanding of the shared inherited basis across human cancers and expand our knowledge of the 5p15.33 region in carcinogenesis.

A genetic variant in telomerase reverse transcriptase (TERT) modifies cancer risk in Lynch syndrome patients harbouring pathogenic MSH2 variants

Individuals with Lynch syndrome (LS), have an increased risk of developing cancer. Common genetic variants of telomerase reverse transcriptase (TERT) have been associated with a wide range of cancers, including colorectal cancer (CRC) in LS. We combined genotype data from 1881 LS patients, carrying pathogenic variants in MLH1, MSH2 or MSH6, for rs2075786 (G>A, intronic variant), 1207 LS patients for rs2736108 (C>T, upstream variant) and 1201 LS patients for rs7705526 (C>A, intronic variant). The risk of cancer was estimated by heterozygous/homozygous odds ratio (OR) with mixed-effects logistic regression to adjust for gene/gender/country of sample origin considering family identity. The AA genotype of SNP rs2075786 is associated with 85% higher odds at developing cancer compared to GG genotype in MSH2 pathogenic variant carriers (p = 0.0160). Kaplan-Meier analysis also shows an association for rs2075786; the AA allele for MSH2 variant carriers confers risk for earlier diagnosis of LS cancer (log-rank p = 0.0011). We report a polymorphism in TERT to be a possible modifier of disease risk in MSH2 pathogenic variant carriers. The rs2075786 SNP in TERT is associated with a differential risk of developing cancer for MSH2 pathogenic variant carriers. Use of this information has the potential to personalise screening protocols for LS patients.

Clinicopathological significance of the single nucleotide polymorphism, rs2853669 within the TERT promoter in papillary thyroid carcinoma

Papillary thyroid carcinoma (PTC) is the most common thyroid malignancy. Point mutations in the telomerase reverse transcriptase (TERT) promoter, C228T and C250T and oncogene BRAF^V600E have been investigated as risk factors for PTC. However, little research has been done on the single nucleotide polymorphism rs2853669 in the TERT promoter in PTC. This study aimed to clarify the clinicopathological significance of rs2853669 in Japanese patients with PTC. The genetic frequencies of rs2853669, C228T, C250T and BRAF^V600E were investigated in 58 patients with PTC and compared with the clinicopathological parameters of PTC. rs2853669, C228T, C250T and BRAF^V600E were found in 58.6%, 17.2%, 5.2% and 37.0% of the PTC patients, respectively. PTC with rs2853669 and C228T were associated only with tumor sizes larger than 2.0 cm (P < 0.05). Furthermore, the coexistence of rs2853669 and C228T was strongly associated with tumor size (P < 0.01), with an odds ratio of 6.4 (P < 0.05). We showed that rs2853669, as well as C228T, may be a risk factor for the aggressiveness of PTC, and the coexistence of these mutations might represent greater risk.

The association between polymorphisms in tankyrase gene and telomere length in omethoate-exposed workers

Peripheral blood leukocyte telomere length in omethoate-exposed workers is related to environmental exposure and single nucleotide polymorphisms (SNPs) in genes including p21, GSTM1, miR-145, etc. However, the roles of SNPs in tankyrase (TNKS) gene in telomere length are still unknown. The aim of this study was to explore the association between SNPs in TNKS gene and telomere length in omethoate-exposed workers. Telomere length in peripheral blood leukocyte DNA from 180 omethoate-exposed workers and 115 healthy controls was measured using Real-time quantitative polymerase chain reaction (PCR). Genotyping of the selected functional and susceptible SNPs was performed by the flight mass spectrometry based on PCR and single-base extension. The analysis of covariance was performed to find effects of SNPs on telomere length. Generalized linear models were used to analyze the environment, gene, and interaction on telomere length. The results showed that telomere length in the CG + CC genotypes in rs1055328 in TNKS gene was significantly longer than that in the wild homozygous GG genotype both in exposure group (P = 0.017) and in control group (P = 0.038) after adjusting the covariates. The variables kept in the generalized linear models included omethoate-exposure (β = 0.580, P = 0.001) and rs1055328 (CG + CC) in TNKS gene (β = 0.339, P = 0.002). The study suggests that the prolongation of telomere length is associated with omethoate-exposure and the CG + CC genotypes in rs1055328 in TNKS gene.

Genes associated with Type 2 Diabetes and vascular complications

Type 2 Diabetes (T2D) is a chronic disease associated with a number of micro- and macrovascular complications that increase the morbidity and mortality of patients. The risk of diabetic complications has a strong genetic component. To this end, we sought to evaluate the association of 40 single nucleotide polymorphisms (SNPs) in 21 candidate genes with T2D and its vascular complications in 503 T2D patients and 580 healthy controls. The genes were chosen because previously reported to be associated with T2D complications and/or with the aging process. We replicated the association of T2D risk with IGF2BP rs4402960 and detected novel associations with TERT rs2735940 and rs2736098. The addition of these SNPs to a model including traditional risk factors slightly improved risk prediction. After stratification of patients according to the presence/absence of vascular complications, we found significant associations of variants in the CAT, FTO, and UCP1 genes with diabetic retinopathy and nephropathy. Additionally, a variant in the ADIPOQ gene was found associated with macrovascular complications. Notably, these genes are involved in some way in mitochondrial biology and reactive oxygen species regulation. Hence, our findings strongly suggest a potential link between mitochondrial oxidative homeostasis and individual predisposition to diabetic vascular complications.

Construction of JRG (Japanese reference genome) with single-molecule real-time sequencing

In recent genome analyses, population-specific reference panels have indicated important. However, reference panels based on short-read sequencing data do not sufficiently cover long insertions. Therefore, the nature of long insertions has not been well documented. Here, we assembled a Japanese genome using single-molecule real-time sequencing data and characterized insertions found in the assembled genome. We identified 3691 insertions ranging from 100 bps to ~10,000 bps in the assembled genome relative to the international reference sequence (GRCh38). To validate and characterize these insertions, we mapped short-reads from 1070 Japanese individuals and 728 individuals from eight other populations to insertions integrated into GRCh38. With this result, we constructed JRGv1 (Japanese Reference Genome version 1) by integrating the 903 verified insertions, totaling 1,086,173 bases, shared by at least two Japanese individuals into GRCh38. We also constructed decoyJRGv1 by concatenating 3559 verified insertions, totaling 2,536,870 bases, shared by at least two Japanese individuals or by six other assemblies. This assembly improved the alignment ratio by 0.4% on average. These results demonstrate the importance of refining the reference assembly and creating a population-specific reference genome. JRGv1 and decoyJRGv1 are available at the JRG website.

Researchers in Japan have assembled a Japanese reference genome, which includes sequences missing from the international reference genome, as well as others specific to East Asian populations. A team led by Masao Nagasaki and Masayuki Yamamoto sequenced a Japanese individual using a method, which produces longer sequences than previous technologies. Using this approach, they identified thousands of sequences spanning 2.5 million bases, which were absent in the international reference genome. Many of these were sequences able to move within the genome. They showed that the majority of these sequences are also present in early humans and chimpanzees, demonstrating that their absence from the current reference is due to deletions or limitations of earlier sequencing methodologies. In addition to providing a population-specific reference, these findings demonstrate the importance of continually improving the international reference genome.

Candidate susceptibility variants in angioimmunoblastic T-cell lymphoma

Angioimmunoblastic T-cell lymphoma (AITL) is a subtype of peripheral T-cell lymphoma with a poor prognosis: the 5-year survival rate is approximately 30%. Somatic driver mutations have been found in TET2, IDH2, DNMT3A, RHOA, FYN, PLCG1, and CD28, whereas germline susceptibility to AITL has to our knowledge not been studied. The homogenous Finnish population is well suited for studies on genetic predisposition. Here, we performed an exome-wide rare variant analysis in 23 AITL patients. No germline mutations were found in the driver genes, implying that they are not frequently involved in genetic AITL predisposition. Potentially pathogenic variants present in at least two patients and showing significant (p < 0.01) enrichment in our sample set were found in ten genes: POLK, PRKCB, ZNF676, PRRC2B, PCDHGB6, GNL3L, TTC36, OTOG, OSGEPL1, and RASSF9. The most significantly enriched variants, causing p.Lys469Ter in a splice variant of POLK and p.Pro588His in PRKCB, are intriguing candidates as Polk deficient mice display a spontaneous mutator phenotype, whereas PRKCB was recently shown to be somatically mutated in 33% of another peripheral T-cell lymphoma, adult T-cell lymphoma. If validated, our findings would provide new insight into the pathogenesis of AITL, as well as tools for early detection in susceptible individuals.

The contribution of parent-to-offspring transmission of telomeres to the heritability of telomere length in humans

Leukocyte telomere length (LTL) is a heritable trait with two potential sources of heritability (h²): inherited variation in non-telomeric regions (e.g., SNPs that influence telomere maintenance) and variability in the lengths of telomeres in gametes that produce offspring zygotes (i.e., "direct" inheritance). Prior studies of LTL h² have not attempted to disentangle these two sources. Here, we use a novel approach for detecting the direct inheritance of telomeres by studying the association between identity-by-descent (IBD) sharing at chromosome ends and phenotypic similarity in LTL. We measured genome-wide SNPs and LTL for a sample of 5069 Bangladeshi adults with substantial relatedness. For each of the 6318 relative pairs identified, we used SNPs near the telomeres to estimate the number of chromosome ends shared IBD, a proxy for the number of telomeres shared IBD (T_shared). We then estimated the association between T_shared and the squared pairwise difference in LTL ((ΔLTL)²) within various classes of relatives (siblings, avuncular, cousins, and distant), adjusting for overall genetic relatedness (ϕ). The association between T_shared and (ΔLTL)² was inverse among all relative pair types. In a meta-analysis including all relative pairs (ϕ > 0.05), the association between T_shared and (ΔLTL)² (P = 0.01) was stronger than the association between ϕ and (ΔLTL)² (P = 0.43). Our results provide strong evidence that telomere length (TL) in parental germ cells impacts TL in offspring cells and contributes to LTL h² despite telomere "reprogramming" during embryonic development. Applying our method to larger studies will enable robust estimation of LTL h² attributable to direct transmission of telomeres.

Association of TERT polymorphisms with chronic hepatitis B in a Chinese Han population

In this study, we investigated the association between the polymorphisms of telomerase reverse transcriptase (TERT) gene and the risk of chronic hepatitis B (CHB) in a Chinese Han population. Four single nucleotide polymorphisms (SNPs) in TERT (rs10069690, rs2242652, rs2853677 and rs2853676) were genotyped from 224 CHB patients and 300 healthy controls using the Sequenom Mass-ARRAY platform. We used genetic model, haplotype analyses, chi-square test, logistic regression analysis to evaluate the association between SNPs and CHB risk. The relative risk was estimated by odd ratios (ORs) and 95% confidence intervals (CIs). We found that rs10069690 was significantly associated with an increased CHB risk in the dominant model (adjusted OR = 1.70, 95% CI: 1.06–2.71, P = 0.031) and additive model (adjusted OR = 1.62, 95% CI: 1.09–2.41, P = 0.018). The haplotype “TA” (rs10069690 and rs2242652) was found to be associated with an increased risk of CHB (adjusted OR = 1.58, 95% CI: 1.05–2.38, P = 0.027). Our results suggested potential genetic contributes for TERT in CHB development in a Chinese Han population. Future functional and association studies with larger sample sizes are required to confirm these findings.

Transcriptome-wide association study revealed two novel genes associated with nonobstructive azoospermia in a Chinese population

OBJECTIVE

To investigate the associations between genetically cis-regulated gene expression levels and nonobstructive azoospermia (NOA) susceptibility.

DESIGN

Transcriptome-wide association study (TWAS).

SETTING

Medical university.

INTERVENTIONS

None.

MAIN OUTCOME MEASURE(S)

The cis-hg² values for each gene were estimated with GCTA software. The effect sizes of cis-single-nucleotide polymorphisms (SNPs) on gene expression were measured using GEMMA software. Gene expression levels were entered into our existing NOA GWAS cohort using GEMMA software. The TWAS P-values were calculated using logistic regression models.

RESULT(S)

Expression levels of 1,296 cis-heritable genes were entered into our existing NOA GWAS data. The TWAS results identified two novel genes as statistically significantly associated with NOA susceptibility: PILRA and ZNF676. In addition, 6p21.32, previously reported in NOA GWAS, was further validated to be a susceptible region to NOA risk.

CONCLUSION(S)

Analysis with TWAS provides fruitful targets for follow-up functional studies.

Genome-wide association study of telomere length among South Asians identifies a second RTEL1 association signal

Background

Leucocyte telomere length (TL) is a potential biomarker of ageing and risk for age-related disease. Leucocyte TL is heritable and shows substantial differences by race/ethnicity. Recent genome-wide association studies (GWAS) report ~10 loci harbouring SNPs associated with leucocyte TL, but these studies focus primarily on populations of European ancestry.

Objective

This study aims to enhance our understanding of genetic determinants of TL across populations.

Methods

We performed a GWAS of TL using data on 5075 Bangladeshi adults. We measured TL using one of two technologies (qPCR or a Luminex-based method) and used standardised variables as TL phenotypes.

Results

Our results replicate previously reported associations in the TERC and TERT regions (P=2.2×10⁻⁸ and P=6.4×10⁻⁶, respectively). We observed a novel association signal in the RTEL1 gene (intronic SNP rs2297439; P=2.82×10⁻⁷) that is independent of previously reported TL-associated SNPs in this region. The minor allele for rs2297439 is common in South Asian populations (≥0.25) but at lower frequencies in other populations (eg, 0.07 in Northern Europeans). Among the eight other previously reported association signals, all were directionally consistent with our study, but only rs8105767 (ZNF208) was nominally significant (P=0.003). SNP-based heritability estimates were as high as 44% when analysing close relatives but much lower when analysing distant relatives only.

Conclusions

In this first GWAS of TL in a South Asian population, we replicate some, but not all, of the loci reported in prior GWAS of individuals of European ancestry, and we identify a novel second association signal at the RTEL1 locus.

Confirmation of five novel susceptibility loci for systemic lupus erythematosus (SLE) and integrated network analysis of 82 SLE susceptibility loci

We recently identified ten novel SLE susceptibility loci in Asians and uncovered several additional suggestive loci requiring further validation. This study aimed to replicate five of these suggestive loci in a Han Chinese cohort from Hong Kong, followed by meta-analysis (11,656 cases and 23,968 controls) on previously reported Asian and European populations, and to perform bioinformatic analyses on all 82 reported SLE loci to identify shared regulatory signatures. We performed a battery of analyses for these five loci, as well as joint analyses on all 82 SLE loci. All five loci passed genome-wide significance: MYNN (rs10936599, Pmeta = 1.92 × 10-13, OR = 1.14), ATG16L2 (rs11235604, Pmeta = 8.87 × 10 -12, OR = 0.78), CCL22 (rs223881, Pmeta = 5.87 × 10-16, OR = 0.87), ANKS1A (rs2762340, Pmeta = 4.93 × 10-15, OR = 0.87) and RNASEH2C (rs1308020, Pmeta = 2.96 × 10-19, OR = 0.84) and co-located with annotated gene regulatory elements. The novel loci share genetic signatures with other reported SLE loci, including effects on gene expression, transcription factor binding, and epigenetic characteristics. Most (56%) of the correlated (r2 > 0.8) SNPs from the 82 SLE loci were implicated in differential expression (9.81 × 10-198 < P < 5 × 10-3) of cis-genes. Transcription factor binding sites for p53, MEF2A and E2F1 were significantly (P < 0.05) over-represented in SLE loci, consistent with apoptosis playing a critical role in SLE. Enrichment analysis revealed common pathways, gene ontology, protein domains, and cell type-specific expression. In summary, we provide evidence of five novel SLE susceptibility loci. Integrated bioinformatics using all 82 loci revealed that SLE susceptibility loci share many gene regulatory features, suggestive of conserved mechanisms of SLE etiopathogenesis.

Genetic variants affecting telomere length are associated with the prognosis of esophageal squamous cell carcinoma in a Chinese population

Telomeres are essential for maintaining chromosomal stability and are crucial in tumor progression. Previous studies have explored the associations between telomere length and cancer prognosis, but the findings are inconclusive. Genome-wide association studies have identified several genetic variants associated with telomere length in Caucasians. However, the roles of telomere length and related genetic variants on esophageal squamous cell carcinoma (ESCC) prognosis are largely unknown. Therefore, we conducted a case-cohort study with 431 ESCC patients to assess the associations between relative telomere length (RTL), eight known telomere length related variants and the overall survival of ESCC in Chinese population. We found that as compared with the reference group, patients in the fifth (the longest) quintile had a significantly better prognosis [(adjusted hazard ratio (HR) = 0.58, 95% confidence interval (CI) = 0.34-0.98, P = 0.041]. Furthermore, A allele of rs2736108 was significantly associated with both the increased RTL (P = 0.048) and the better prognosis of ESCC (adjusted HR = 0.55, 95%CI = 0.38-0.79, P = 1.31 × 10^-3 ). Mediation analysis indicated that the effect of rs2736108 on ESCC prognosis was partly explained by RTL (1.99%). Stepwise Cox proportional hazard analysis suggested that rs2736108 played an important protective role in ESCC prognosis (HR = 0.57, 95%CI = 0.40-0.81, P = 1.97 × 10^-3 ). Our findings provide evidence that prolonged telomere length is a protective factor for ESCC patients' survival and the known telomere length related genetic variant rs2736108 can contribute to the prognosis of ESCC as well in Chinese population. © 2016 Wiley Periodicals, Inc.

TERT Polymorphism rs2853669 Influences on Lung Cancer Risk in the Korean Population

Short telomeres are known as one of the risk factors for human cancers. The present study was conducted to evaluate the association between 6 polymorphisms, which were related with short telomere length in the Korean population, and lung cancer risk using 1,100 cases and 1,096 controls. Among the 6 polymorphisms, TERT rs2853669 was significantly associated with increased lung cancer risk under a recessive model (odds ratio [OR]=1.38, 95% confidence interval [CI]=1.05-1.81, P=0.02). The effect of rs2853669 on lung cancer risk was significant in younger individuals (OR=1.73, 95% CI=1.18-2.54, P=0.005) and adenocarcinoma (OR=1.50, 95% CI=1.07-2.07, P=0.02). Our results suggest that a common functional promoter polymorphism, TERT rs2853669, may influence both telomere length and lung cancer risk in the Korean population.