A network approach to exploring the functional basis of gene-gene epistatic interactions in disease susceptibility.

Motivation: Individual genetic variants explain only a small fraction of heritability in some diseases. Some variants have weak marginal effects on disease risk, but their joint effects are significantly stronger when occurring together. Most studies on such epistatic interactions have focused on methods for identifying the interactions and interpreting individual cases, but few have explored their general functional basis. This was due to the lack of a comprehensive list of epistatic interactions and uncertainties in associating variants to genes. Results: We conducted a large-scale survey of published research articles to compile the first comprehensive list of epistatic interactions in human diseases with detailed annotations. We used various methods to associate these variants to genes to ensure robustness. We found that these genes are significantly more connected in protein interaction networks, are more co-expressed and participate more often in the same pathways. We demonstrate using the list to discover novel disease pathways. Contact: kevinyip@cse.cuhk.edu.hk. Supplementary information: Supplementary data are available at Bioinformatics online.

The association of liver function and quality of life of patients with liver cancer.

BACKGROUND: Quality of life (QOL) assessments with the European Organization for Research and Treatment of Cancer (EORTC) QLQ-C30, QLQ-HCC18, C30 and HCC18 index scores have been shown to be prognostic factors for overall survival (OS) in patients with hepatocellular carcinoma (HCC), independent of disease stage and liver function. Liver function parameters (including bilirubin, albumin, international normalized ratio [INR], Child-Pugh class, ALBI grade, MELD, alkaline phosphatase [ALP]-to-platelet ratio, albumin-to-ALP ratio) have also been found to be independent prognostic factors for OS in HCC patients. There has been scanty data on whether QOL and baseline liver function per se are correlated in HCC patients. This study investigates the correlations between baseline QOL data and liver function variables in HCC patients. METHODS: From 2007 to 2011, 517 patients were enrolled. Baseline QOL was assessed at diagnosis using the EORTC QLQ-C30 and QLQ-HCC18; thereafter C30 and HCC18 index scores were derived. Clinical and laboratory data were collected. For liver function assessment, Child-Pugh class, ALBI grade, MELD, ALP-to-platelet ratio and albumin-to-ALP ratio were derived. Correlation analyses were performed between QOL and liver function data. RESULTS: Complete QOL data were available in 472 HCC patients. After adjusting for clinical variables, significant correlations were found between QOL (QLQ-C30 and QLQ-HCC18) and dichotomized liver function variables (including Child-Pugh class, ALBI grade and the presence of ascites). It was demonstrated that QOL had significant and potentially clinically important correlations with continuous liver function variables (albumin, bilirubin, ALP and albumin-to-ALP ratio), with the highest Spearman's rank correlation coefficient (rho) exceeding 0.4. HCC18 and C30 index scores were also significantly correlated with these liver function variables. HCC18 index score, which had rho up to 0.37, generally performed better than C30 index score, which had rho up to 0.33. CONCLUSIONS: In HCC patients, baseline QOL assessment (using EORTC QLQ-C30, QLQ-HCC18, C30 index-score or HCC18 index-score) is significantly correlated with liver function. Based on the findings of this study, future trials are warranted to assess whether treatment to enhance liver function could improve HCC patients' QOL.

Profiles of lipids, blood pressure and weight changes among premenopausal Chinese breast cancer patients after adjuvant chemotherapy.

BACKGROUND: Adjuvant chemotherapy improves outcome of patients with early breast cancer. However, chemotherapy may be associated with long term toxicities. In this retrospective cohort study, the objectives were to determine body weight, body mass index (BMI), blood pressure and fasting lipids levels of young premenopausal Chinese breast cancer patients after adjuvant chemotherapy. Potential factors associated with these parameters were identified. METHODS: Eligibility criteria include premenopausal Chinese patients who were diagnosed to have stage I-III breast cancer within 3-10 years, age < 45 and having received adjuvant chemotherapy at the time of breast cancer diagnosis. Information at initial breast cancer diagnosis were retrieved from patients' medical records and include age at diagnosis, tumor characteristics, anti-cancer treatments, blood pressure and body weight and height. At study entry, all patients had additional background demographics collected, as well as blood pressure, body weight and fasting serum lipid profiles measured. Incidence of chemotherapy-related amenorrhoea (CRA) and menopause were determined. Factors associated with weight gain, hypertension and dyslipidaemias were analyzed. RESULTS: Two hundred and eighty patients were studied. The median age at breast cancer diagnosis was 41 years (range: 24-45). The median time from breast cancer diagnosis to study entry was 5.0 years. The median age at study entry was 46.5 years (range: 28-54). 91.1% developed CRA; 48.9% had become menopausal and 10% were peri-menopausal. Between initial breast cancer diagnosis and the time of study entry, the median weight gain was 1.8 kg; 63.2% gained weight by >2%; 52.1% were overweight/obese; 30.7% had hypertension. Abnormal total-cholesterol and LDL-cholesterol occurred in 34.3% and 56.1% respectively. On multivariate analyses, older age was associated with reduced risk while occurrence of CRA and having received taxane-containing regimens were associated with increased risk of weight gain. Oestrogen-receptor positivity was associated with reduced risk while overweight/obese statuses were associated with increased risk of hypertension. Use of tamoxifen was associated with reduced risk of abnormal LDL-cholesterol. Weight gain, overweight/obese, older age, progression to post/peri-menopausal status at study entry, having received corticosteroid premedication before adjuvant chemotherapy and having received taxane-containing adjuvant chemotherapy were associated with increased risk of dyslipidaemias. CONCLUSION: Among young premenopausal Chinese breast cancer patients who had received adjuvant chemotherapy, the current study has revealed that although there was only a median weight gain of 1.8 kg, there was a nearly 60% increase in abnormal BMI. Further, a significant proportion of patients were detected to have hypertension and dyslipidaemias. Interventional studies with lifestyle modifications are warranted.

Functional interaction between SNPs and microsatellite in the transcriptional regulation of insulin-like growth factor 1.

A CA-repeat microsatellite in insulin-like growth factor 1 (IGF1) promoter was associated with interindividual variation of circulating IGF1 level. Previously, we reported that such association was due to variation of haplotype unit in a linkage disequilibrium block composed of microsatellite and single-nucleotide polymorphisms (SNPs), suggesting the presence of an interaction between them. In this study, reporter assays were performed to investigate the regulatory effect and interaction of genetic variants on gene expression. We used an in vitro system to compare the transcriptional activities of haplotypes (rs35767:T>C, the CA-repeat microsatellite, rs5742612:T>C, and rs2288377:T>A) in evolutionarily conserved region of IGF1 promoter. In haplotype C-T-T, a longer microsatellite had a lower transcriptional activity (17.6 ± 2.4-fold for 17 repeats and 8.3 ± 1.1-fold for 21 repeats), whereas in haplotype T-C-A, such trend could not be observed, as the microsatellite with 21 repeats had the highest transcriptional activity (17.5 ± 2.3-fold). Because the microsatellite and SNPs affected the transcriptional activity of each other, there may be an interaction between them in the regulation of IGF1 expression. For the first time, we demonstrated that a noncoding microsatellite polymorphism could act as a functional unit and interact with SNPs in the regulation of transcription in human genome.

Genetic adaptation of the hypoxia-inducible factor pathway to oxygen pressure among eurasian human populations.

Research into the mechanisms of human adaptation to the hypoxic environment of high altitude is of great interest to the fields of human physiology and clinical medicine. Recently, the gene EGLN1, from the hypoxia-inducible factor (HIF) pathway, was identified as being involved in the hypoxic adaptation of highland Andeans and Tibetans. Both highland Andeans and Tibetans have adapted to an extremely hypoxic habitat and less attention has been paid to populations living in normoxic conditions at sea level and mild-hypoxic environments of moderate altitude, thus, whether a common adaptive mechanism exists in response to quantitative variations of environmental oxygen pressure over a wide range of residing altitudes is unknown. Here, we first performed a genome-wide association study of 35 populations from the Human Genome Diversity-CEPH Panel who dwell at sea level to moderate altitude in Eurasia (N = 691; 0-2,500 m) to identify the genetic adaptation profile of normoxic and mild-hypoxic inhabitants. In addition, we systematically compared the results from the present study to six previously published genome-wide scans of highland Andeans and Tibetans to identify shared adaptive signals in response to quantitative variations of oxygen pressure. For normoxic and mild-hypoxic populations, the strongest adaptive signal came from the mu opioid receptor-encoding gene (OPRM1, 2.54 × 10(-9)), which has been implicated in the stimulation of respiration, while in the systematic survey the EGLN1-DISC1 locus was identified in all studies. A replication study performed with highland Tibetans (N = 733) and sea level Han Chinese (N = 748) confirmed the association between altitude and SNP allele frequencies in OPRM1 (in Tibetans only, P < 0.01) and in EGLN1-DISC1 (in Tibetans and Han Chinese, P < 0.01). Taken together, identification of the OPRM1 gene suggests that cardiopulmonary adaptation mechanisms are important and should be a focus in future studies of hypoxia adaptation. Furthermore, the identification of the EGLN1 gene from the HIF pathway suggests a common adaptive mechanism for Eurasian human populations residing at different altitudes with different oxygen pressures.

BOOST: A fast approach to detecting gene-gene interactions in genome-wide case-control studies.

Gene-gene interactions have long been recognized to be fundamentally important for understanding genetic causes of complex disease traits. At present, identifying gene-gene interactions from genome-wide case-control studies is computationally and methodologically challenging. In this paper, we introduce a simple but powerful method, named "BOolean Operation-based Screening and Testing" (BOOST). For the discovery of unknown gene-gene interactions that underlie complex diseases, BOOST allows examination of all pairwise interactions in genome-wide case-control studies in a remarkably fast manner. We have carried out interaction analyses on seven data sets from the Wellcome Trust Case Control Consortium (WTCCC). Each analysis took less than 60 hr to completely evaluate all pairs of roughly 360,000 SNPs on a standard 3.0 GHz desktop with 4G memory running the Windows XP system. The interaction patterns identified from the type 1 diabetes data set display significant difference from those identified from the rheumatoid arthritis data set, although both data sets share a very similar hit region in the WTCCC report. BOOST has also identified some disease-associated interactions between genes in the major histocompatibility complex region in the type 1 diabetes data set. We believe that our method can serve as a computationally and statistically useful tool in the coming era of large-scale interaction mapping in genome-wide case-control studies.

CISH and susceptibility to infectious diseases.

BACKGROUND: The interleukin-2-mediated immune response is critical for host defense against infectious pathogens. Cytokine-inducible SRC homology 2 (SH2) domain protein (CISH), a suppressor of cytokine signaling, controls interleukin-2 signaling. METHODS: Using a case-control design, we tested for an association between CISH polymorphisms and susceptibility to major infectious diseases (bacteremia, tuberculosis, and severe malaria) in blood samples from 8402 persons in Gambia, Hong Kong, Kenya, Malawi, and Vietnam. We had previously tested 20 other immune-related genes in one or more of these sample collections. RESULTS: We observed associations between variant alleles of multiple CISH polymorphisms and increased susceptibility to each infectious disease in each of the study populations. When all five single-nucleotide polymorphisms (SNPs) (at positions -639, -292, -163, +1320, and +3415 [all relative to CISH]) within the CISH-associated locus were considered together in a multiple-SNP score, we found an association between CISH genetic variants and susceptibility to bacteremia, malaria, and tuberculosis (P=3.8x10(-11) for all comparisons), with -292 accounting for most of the association signal (P=4.58x10(-7)). Peripheral-blood mononuclear cells obtained from adult subjects carrying the -292 variant, as compared with wild-type cells, showed a muted response to the stimulation of interleukin-2 production--that is, 25 to 40% less CISH expression. CONCLUSIONS: Variants of CISH are associated with susceptibility to diseases caused by diverse infectious pathogens, suggesting that negative regulators of cytokine signaling have a role in immunity against various infectious diseases. The overall risk of one of these infectious diseases was increased by at least 18% among persons carrying the variant CISH alleles.

Lack of association between polymorphisms from genome-wide association studies and tuberculosis in the Chinese population.

Genome-wide association studies (GWAS) have identified rs4331426 and rs2057178 as being associated with tuberculosis (TB) in African populations. Both are common single nucleotide polymorphisms (SNPs) in Africans, but they are much rarer in Eurasian populations. In order to corroborate these results, we carried out a case-control study in the Chinese population; these 2 SNPs were genotyped in 600 pulmonary TB patients and 618 healthy controls. The results showed that neither of the SNPs was associated with TB, even after stratification by gender, age, and smear status. Considering the limitation of poor coverage of variations in commercial available genotyping platforms in African populations, further GWAS should be conducted in other populations such as Indian and Chinese. Moreover, future genetic studies on host susceptibility to TB need to take into account all the variables, including host, environment, pathogen, and interactions.

Telomere length and cognitive function in southern Chinese community-dwelling male elders.

BACKGROUND: telomere attrition has been associated with an increased risk of different age-related diseases and is widely accepted as a marker of cellular ageing. On the other hand, it is well known that cognitive function declines with age. The telomere length may therefore act as a marker for the pathway associated with cognitive function. METHODS: we examined telomere length and cognitive functions in a community-dwelling Chinese male population aged 65 years and above living in Hong Kong. The telomere length was measured by quantitative real-time PCR in 976 men. Cognitive function was assessed by Chinese (Cantonese) version of Mini-Mental State Exam and Community Screening Interview for Dementia. RESULTS: our result showed there was a significant association between telomere length, delayed recall (P = 0.007) and category verbal fluency (P = 0.048). These associations remained significant after adjustment for age and education. Further analysis using a cut-off score for MMSE, three-item recall and word list generation tests suggested that the telomere length was positively correlated with performance in these areas (P = 0.015). CONCLUSION: the findings support the association of telomere length and cognitive function and suggested that the telomere length may serve as a biological marker for cognitive decline.

Identification of a Common Variant for Coronary Heart Disease at PDE1A Contributes to Individualized Treatment Goals and Risk Stratification of Cardiovascular Complications in Chinese Patients With Type 2 Diabetes.

OBJECTIVE: In this study we aim to unravel genetic determinants of coronary heart disease (CHD) in type 2 diabetes (T2D) and explore their applications. RESEARCH DESIGN AND METHODS: We performed a two-stage genome-wide association study for CHD in Chinese patients with T2D (3,596 case and 8,898 control subjects), followed by replications in European patients with T2D (764 case and 4,276 control subjects) and general populations (n = 51,442-547,261). Each identified variant was examined for its association with a wide range of phenotypes and its interactions with glycemic, blood pressure (BP), and lipid controls in incident cardiovascular diseases. RESULTS: We identified a novel variant (rs10171703) for CHD (odds ratio 1.21 [95% CI 1.13-1.30]; P = 2.4 × 10-8) and BP (ß ± SE 0.130 ± 0.017; P = 4.1 × 10-14) at PDE1A in Chinese T2D patients but found only a modest association with CHD in general populations. This variant modulated the effects of BP goal attainment (130/80 mmHg) on CHD (Pinteraction = 0.0155) and myocardial infarction (MI) (Pinteraction = 5.1 × 10-4). Patients with CC genotype of rs10171703 had >40% reduction in either cardiovascular events in response to BP control (2.9 × 10-8 < P < 3.6 × 10-5), those with CT genotype had no difference (0.0726 < P < 0.2614), and those with TT genotype had a threefold increase in MI risk (P = 6.7 × 10-3). CONCLUSIONS: We discovered a novel CHD- and BP-related variant at PDE1A that interacted with BP goal attainment with divergent effects on CHD risk in Chinese patients with T2D. Incorporating this information may facilitate individualized treatment strategies for precision care in diabetes, only when our findings are validated.

Nonsense-Mediated Decay Targeted RNA (ntRNA): Proposal of a ntRNA-miRNA-lncRNA Triple Regulatory Network Usable as Biomarker of Prognostic Risk in Patients with Kidney Cancer.

The most prevalent subtype of renal cell carcinoma (RCC), kidney renal clear cell carcinoma (KIRC) may be associated with a poor prognosis in a high number of cases, with a stage-specific prognostic stratification currently in use. No reliable biomarkers have been utilized so far in clinical practice despite the efforts in biomarker research in the last years. Nonsense-mediated mRNA decay (NMD) is a critical safeguard against erroneous transcripts, particularly mRNA transcripts containing premature termination codons (called nonsense-mediated decay targeted RNA, ntRNA). In this study, we first characterized 296 differentially expressed ntRNAs that were independent of the corresponding gene, 261 differentially expressed miRNAs, and 4653 differentially expressed lncRNAs. Then, we constructed a hub ntRNA-miRNA-lncRNA triple regulatory network associated with the prognosis of KIRC. Moreover, the results of immune infiltration analysis indicated that this network may influence the changes of the tumor immune microenvironment. A prognostic model derived from the genes and immune cells associated with the network was developed to distinguish between high- and low-risk patients, which was a better prognostic than other models, constructed using different biomarkers. Additionally, correlation of methylation and ntRNAs in the network suggested that some ntRNAs were regulated by methylation, which is helpful to further study the causes of abnormal expression of ntRNAs. In conclusion, this study highlighted the possible clinical implications of ntRNA functions in KIRC, proposing potential significant biomarkers that could be utilized to define the prognosis and design personalized treatment plans in kidney cancer management in the next future.

Clinical Predictors and Long-term Impact of Acute Kidney Injury on Progression of Diabetic Kidney Disease in Chinese Patients With Type 2 Diabetes.

We aim to assess the long-term impact of acute kidney injury (AKI) on progression of diabetic kidney disease (DKD) and all-cause mortality and investigate determinants of AKI in Chinese patients with type 2 diabetes (T2D). A consecutive cohort of 9,096 Chinese patients with T2D from the Hong Kong Diabetes Register was followed for 12 years (mean ± SD age 57 ± 13.2 years; 46.9% men; median duration of diabetes 5 years). AKI was defined based on the Kidney Disease: Improving Global Outcomes (KDIGO) criteria using serum creatinine. Estimated glomerular filtration rate measurements were used to identify the first episode with chronic kidney disease (CKD) and end-stage renal disease (ESRD). Polygenic risk score (PRS) composed of 27 single nucleotide polymorphisms (SNPs) known to be associated with serum uric acid (SUA) in European populations was used to examine the role of SUA in pathogenesis of AKI, CKD, and ESRD. Validation was sought in an independent cohort including 6,007 patients (age 61.2 ± 10.9 years; 59.5% men; median duration of diabetes 10 years). Patients with AKI had a higher risk for developing incident CKD (hazard ratio 14.3 [95% CI 12.69-16.11]), for developing ESRD (12.1 [10.74-13.62]), and for all-cause death (7.99 [7.31-8.74]) compared with those without AKI. Incidence rate for ESRD among patients with no episodes of AKI and one, two, and three or more episodes of AKI was 7.1, 24.4, 32.4, and 37.3 per 1,000 person-years, respectively. Baseline SUA was a strong independent predictor for AKI. A PRS composed of 27 SUA-related SNPs was associated with AKI and CKD in both discovery and replication cohorts but not ESRD. Elevated SUA may increase the risk of DKD through increasing AKI. The identification of SUA as a modifiable risk factor and PRS as a nonmodifiable risk factor may facilitate the identification of individuals at high risk to prevent AKI and its long-term impact in T2D.

A hidden two-locus disease association pattern in genome-wide association studies.

BACKGROUND: Recent association analyses in genome-wide association studies (GWAS) mainly focus on single-locus association tests (marginal tests) and two-locus interaction detections. These analysis methods have provided strong evidence of associations between genetics variances and complex diseases. However, there exists a type of association pattern, which often occurs within local regions in the genome and is unlikely to be detected by either marginal tests or interaction tests. This association pattern involves a group of correlated single-nucleotide polymorphisms (SNPs). The correlation among SNPs can lead to weak marginal effects and the interaction does not play a role in this association pattern. This phenomenon is due to the existence of unfaithfulness: the marginal effects of correlated SNPs do not express their significant joint effects faithfully due to the correlation cancelation. RESULTS: In this paper, we develop a computational method to detect this association pattern masked by unfaithfulness. We have applied our method to analyze seven data sets from the Wellcome Trust Case Control Consortium (WTCCC). The analysis for each data set takes about one week to finish the examination of all pairs of SNPs. Based on the empirical result of these real data, we show that this type of association masked by unfaithfulness widely exists in GWAS. CONCLUSIONS: These newly identified associations enrich the discoveries of GWAS, which may provide new insights both in the analysis of tagSNPs and in the experiment design of GWAS. Since these associations may be easily missed by existing analysis tools, we can only connect some of them to publicly available findings from other association studies. As independent data set is limited at this moment, we also have difficulties to replicate these findings. More biological implications need further investigation. AVAILABILITY: The software is freely available at http://bioinformatics.ust.hk/hidden_pattern_finder.zip.

Detecting two-locus associations allowing for interactions in genome-wide association studies.

MOTIVATION: Genome-wide association studies (GWASs) aim to identify genetic susceptibility to complex diseases by assaying and analyzing hundreds of thousands of single nucleotide polymorphisms (SNPs). Although traditional single-locus statistical tests have identified many genetic determinants of susceptibility, those findings cannot completely explain genetic contributions to complex diseases. Marchini and coauthors demonstrated the importance of testing two-locus associations allowing for interactions through a wide range of simulation studies. However, such a test is computationally demanding as we need to test hundreds of billions of SNP pairs in GWAS. Here, we provide a method to address this computational burden for dichotomous phenotypes. RESULTS: We have applied our method on nine datasets from GWAS, including the aged-related macular degeneration (AMD) dataset, the Parkinson's disease dataset and seven datasets from the Wellcome Trust Case Control Consortium (WTCCC). Our method has discovered many associations that were not identified before. The running time for the AMD dataset, the Parkinson's disease dataset and each of seven WTCCC datasets are 2.5, 82 and 90 h on a standard 3.0 GHz desktop with 4 G memory running Windows XP system. Our experiment results demonstrate that our method is feasible for the full-scale analyses of both single- and two-locus associations allowing for interactions in GWAS. AVAILABILITY: http://bioinformatics.ust.hk/SNPAssociation.zip CONTACT: nelsontang@cuhk.edu.hk; eeyu@ust.hk; SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Predictive rule inference for epistatic interaction detection in genome-wide association studies.

MOTIVATION: Under the current era of genome-wide association study (GWAS), finding epistatic interactions in the large volume of SNP data is a challenging and unsolved issue. Few of previous studies could handle genome-wide data due to the difficulties in searching the combinatorially explosive search space and statistically evaluating high-order epistatic interactions given the limited number of samples. In this work, we propose a novel learning approach (SNPRuler) based on the predictive rule inference to find disease-associated epistatic interactions. RESULTS: Our extensive experiments on both simulated data and real genome-wide data from Wellcome Trust Case Control Consortium (WTCCC) show that SNPRuler significantly outperforms its recent competitor. To our knowledge, SNPRuler is the first method that guarantees to find the epistatic interactions without exhaustive search. Our results indicate that finding epistatic interactions in GWAS is computationally attainable in practice. AVAILABILITY: http://bioinformatics.ust.hk/SNPRuler.zip

Novel mutations in PHKA2 gene in glycogen storage disease type IX patients from Hong Kong, China.

The diagnosis of glycogen storage disease (GSD) type IX is often complicated by the complexity of the phosphorylase kinase enzyme (PHK), and molecular analysis is the preferred way to provide definitive diagnosis. Here we reported two novel mutations found in two GSD type IX patients with different residual enzyme activities from Hong Kong, China using genetic analysis and, provided the molecular interpretation of the deficient PHK activity. These two newly described mutations would be useful for the study of future GSD patients.

Haplotype effect in the IGF1 promoter accounts for the association between microsatellite and serum IGF1 concentration.

OBJECTIVE: Systemic insulin-like growth factor 1 (IGF1) level is an important risk factor for various diseases. The inter-individual variation of serum IGF1 is determined by environmental and genetic factors, which are attributed to a microsatellite in IGF1 promoter. However, the exact nature of the underlying regulatory elements accounting for this association has not been characterized. Here, we defined the haplotype patterns, including both SNPs and the microsatellite, in the Chinese population, and investigated their regulatory effect on serum IGF1 level. This is the first study in which haplotype patterns of the microsatellite and SNPs in the IGF1 promoter are examined together. METHODS: The linkage disequilibrium (LD) patterns of IGF1 were examined using tagSNPs of the IGF1 regulatory region. The microsatellite, three tagSNPs and haplotypes were correlated with serum IGF1 concentration in 450 normal premenopausal Chinese women. RESULTS: Common alleles of the microsatellite were in strong LD with the three tagSNPs and were associated with particular haplotypes composed of SNPs. Neither the CA repeat number nor SNPs alone showed a robust association with serum IGF1 concentration. On the other hand, the haplotype T-19-A-T was significantly associated with serum IGF1 level. CONCLUSION: No association was found between SNPs and microsatellite alone. However, the haplotype showed better correlation with serum IGF1 level. The results indicate that the previously observed correlation with microsatellite was because of a haplotype effect in the IGF1 promoter. Microsatellite or tagSNPs alone are not the primary regulatory elements of IGF1 expression. The exact regulatory genetic variant needs to be defined by functional genetic studies.

Direct Measurement of B Lymphocyte Gene Expression Biomarkers in Peripheral Blood Transcriptomics Enables Early Prediction of Vaccine Seroconversion.

Peripheral blood transcriptome is a highly promising area for biomarker development. However, transcript abundances (TA) in these cell mixture samples are confounded by proportions of the component leukocyte subpopulations. This poses a challenge to clinical applications, as the cell of origin of any change in TA is not known without prior cell separation procedure. We developed a framework to develop a cell-type informative TA biomarkers which enable determination of TA of a single cell-type (B lymphocytes) directly in cell mixture samples of peripheral blood (e.g., peripheral blood mononuclear cells, PBMC) without the need for subpopulation separation. It is applicable to a panel of genes called B cell informative genes. Then a ratio of two B cell informative genes (a target gene and a stably expressed reference gene) obtained in PBMC was used as a new biomarker to represent the target gene expression in purified B lymphocytes. This approach, which eliminates the tedious procedure of cell separation and directly determines TA of a leukocyte subpopulation in peripheral blood samples, is called the Direct LS-TA method. This method is applied to gene expression datasets collected in influenza vaccination trials as early predictive biomarkers of seroconversion. By using TNFRSF17 or TXNDC5 as the target genes and TNFRSF13C or FCRLA as the reference genes, the Direct LS-TA B cell biomarkers were determined directly in the PBMC transcriptome data and were highly correlated with TA of the corresponding target genes in purified B lymphocytes. Vaccination responders had almost a 2-fold higher Direct LS-TA biomarker level of TNFRSF17 (log 2 SMD = 0.84, 95% CI = 0.47-1.21) on day 7 after vaccination. The sensitivity of these Direct LS-TA biomarkers in the prediction of seroconversion was greater than 0.7 and area-under curves (AUC) were over 0.8 in many datasets. In this paper, we report a straightforward approach to directly estimate B lymphocyte gene expression in PBMC, which could be used in a routine clinical setting. Moreover, the method enables the practice of precision medicine in the prediction of vaccination response. More importantly, seroconversion could now be predicted as early as day 7. As the acquired immunology pathway is common to vaccination against influenza and COVID-19, these biomarkers could also be useful to predict seroconversion for the new COVID-19 vaccines.

A Decade in Review after Idiopathic Scoliosis Was First Called a Complex Trait-A Tribute to the Late Dr. Yves Cotrel for His Support in Studies of Etiology of Scoliosis.

Adolescent Idiopathic Scoliosis (AIS) is a prevalent and important spine disorder in the pediatric age group. An increased family tendency was observed for a long time, but the underlying genetic mechanism was uncertain. In 1999, Dr. Yves Cotrel founded the Cotrel Foundation in the Institut de France, which supported collaboration of international researchers to work together to better understand the etiology of AIS. This new concept of AIS as a complex trait evolved in this setting among researchers who joined the annual Cotrel meetings. It is now over a decade since the first proposal of the complex trait genetic model for AIS. Here, we review in detail the vast information about the genetic and environmental factors in AIS pathogenesis gathered to date. More importantly, new insights into AIS etiology were brought to us through new research data under the perspective of a complex trait. Hopefully, future research directions may lead to better management of AIS, which has a tremendous impact on affected adolescents in terms of both physical growth and psychological development.

Development of genome-wide polygenic risk scores for lipid traits and clinical applications for dyslipidemia, subclinical atherosclerosis, and diabetes cardiovascular complications among East Asians.

BACKGROUND: The clinical utility of personal genomic information in identifying individuals at increased risks for dyslipidemia and cardiovascular diseases remains unclear. METHODS: We used data from Biobank Japan (n = 70,657-128,305) and developed novel East Asian-specific genome-wide polygenic risk scores (PRSs) for four lipid traits. We validated (n = 4271) and subsequently tested associations of these scores with 3-year lipid changes in adolescents (n = 620), carotid intima-media thickness (cIMT) in adult women (n = 781), dyslipidemia (n = 7723), and coronary heart disease (CHD) (n = 2374 cases and 6246 controls) in type 2 diabetes (T2D) patients. RESULTS: Our PRSs aggregating 84-549 genetic variants (0.251 < correlation coefficients (r) < 0.272) had comparably stronger association with lipid variations than the typical PRSs derived based on the genome-wide significant variants (0.089 < r < 0.240). Our PRSs were robustly associated with their corresponding lipid levels (7.5 × 10- 103 < P < 1.3 × 10- 75) and 3-year lipid changes (1.4 × 10- 6 < P < 0.0130) which started to emerge in childhood and adolescence. With the adjustments for principal components (PCs), sex, age, and body mass index, there was an elevation of 5.3% in TC (ß ± SE = 0.052 ± 0.002), 11.7% in TG (ß ± SE = 0.111 ± 0.006), 5.8% in HDL-C (ß ± SE = 0.057 ± 0.003), and 8.4% in LDL-C (ß ± SE = 0.081 ± 0.004) per one standard deviation increase in the corresponding PRS. However, their predictive power was attenuated in T2D patients (0.183 < r < 0.231). When we included each PRS (for TC, TG, and LDL-C) in addition to the clinical factors and PCs, the AUC for dyslipidemia was significantly increased by 0.032-0.057 in the general population (7.5 × 10- 3 < P < 0.0400) and 0.029-0.069 in T2D patients (2.1 × 10- 10 < P < 0.0428). Moreover, the quintile of TC-related PRS was moderately associated with cIMT in adult women (ß ± SE = 0.011 ± 0.005, Ptrend = 0.0182). Independent of conventional risk factors, the quintile of PRSs for TC [OR (95% CI) = 1.07 (1.03-1.11)], TG [OR (95% CI) = 1.05 (1.01-1.09)], and LDL-C [OR (95% CI) = 1.05 (1.01-1.09)] were significantly associated with increased risk of CHD in T2D patients (4.8 × 10- 4 < P < 0.0197). Further adjustment for baseline lipid drug use notably attenuated the CHD association. CONCLUSIONS: The PRSs derived and validated here highlight the potential for early genomic screening and personalized risk assessment for cardiovascular disease.