Stroke genetics informs drug discovery and risk prediction across ancestries.

Previous genome-wide association studies (GWASs) of stroke - the second leading cause of death worldwide - were conducted predominantly in populations of European ancestry1,2. Here, in cross-ancestry GWAS meta-analyses of 110,182 patients who have had a stroke (five ancestries, 33% non-European) and 1,503,898 control individuals, we identify association signals for stroke and its subtypes at 89 (61 new) independent loci: 60 in primary inverse-variance-weighted analyses and 29 in secondary meta-regression and multitrait analyses. On the basis of internal cross-ancestry validation and an independent follow-up in 89,084 additional cases of stroke (30% non-European) and 1,013,843 control individuals, 87% of the primary stroke risk loci and 60% of the secondary stroke risk loci were replicated (P < 0.05). Effect sizes were highly correlated across ancestries. Cross-ancestry fine-mapping, in silico mutagenesis analysis3, and transcriptome-wide and proteome-wide association analyses revealed putative causal genes (such as SH3PXD2A and FURIN) and variants (such as at GRK5 and NOS3). Using a three-pronged approach4, we provide genetic evidence for putative drug effects, highlighting F11, KLKB1, PROC, GP1BA, LAMC2 and VCAM1 as possible targets, with drugs already under investigation for stroke for F11 and PROC. A polygenic score integrating cross-ancestry and ancestry-specific stroke GWASs with vascular-risk factor GWASs (integrative polygenic scores) strongly predicted ischaemic stroke in populations of European, East Asian and African ancestry5. Stroke genetic risk scores were predictive of ischaemic stroke independent of clinical risk factors in 52,600 clinical-trial participants with cardiometabolic disease. Our results provide insights to inform biology, reveal potential drug targets and derive genetic risk prediction tools across ancestries.

Brugada syndrome in Japan and Europe: a genome-wide association study reveals shared genetic architecture and new risk loci.

BACKGROUND AND AIMS: Brugada syndrome (BrS) is an inherited arrhythmia with a higher disease prevalence and more lethal arrhythmic events in Asians than in Europeans. Genome-wide association studies (GWAS) have revealed its polygenic architecture mainly in European populations. The aim of this study was to identify novel BrS-associated loci and to compare allelic effects across ancestries. METHODS: A GWAS was conducted in Japanese participants, involving 940 cases and 1634 controls, followed by a cross-ancestry meta-analysis of Japanese and European GWAS (total of 3760 cases and 11 635 controls). The novel loci were characterized by fine-mapping, gene expression, and splicing quantitative trait associations in the human heart. RESULTS: The Japanese-specific GWAS identified one novel locus near ZSCAN20 (P = 1.0 × 10-8), and the cross-ancestry meta-analysis identified 17 association signals, including six novel loci. The effect directions of the 17 lead variants were consistent (94.1%; P for sign test = 2.7 × 10-4), and their allelic effects were highly correlated across ancestries (Pearson's R = .91; P = 2.9 × 10-7). The genetic risk score derived from the BrS GWAS of European ancestry was significantly associated with the risk of BrS in the Japanese population [odds ratio 2.12 (95% confidence interval 1.94-2.31); P = 1.2 × 10-61], suggesting a shared genetic architecture across ancestries. Functional characterization revealed that a lead variant in CAMK2D promotes alternative splicing, resulting in an isoform switch of calmodulin kinase II-δ, favouring a pro-inflammatory/pro-death pathway. CONCLUSIONS: This study demonstrates novel susceptibility loci implicating potentially novel pathogenesis underlying BrS. Despite differences in clinical expressivity and epidemiology, the polygenic architecture of BrS was substantially shared across ancestries.

Development and validation of genome-wide polygenic risk scores for predicting breast cancer incidence in Japanese females: a population-based case-cohort study.

PURPOSE: This study aimed to develop an ancestry-specific polygenic risk scores (PRSs) for the prediction of breast cancer events in Japanese females and validate it in a longitudinal cohort study. METHODS: Using publicly available summary statistics of female breast cancer genome-wide association study (GWAS) of Japanese and European ancestries, we, respectively, developed 31 candidate genome-wide PRSs using pruning and thresholding (P + T) and LDpred methods with varying parameters. Among the candidate PRS models, the best model was selected using a case-cohort dataset (63 breast cancer cases and 2213 sub-cohorts of Japanese females during a median follow-up of 11.9 years) according to the maximal predictive ability by Harrell's C-statistics. The best-performing PRS for each derivation GWAS was evaluated in another independent case-cohort dataset (260 breast cancer cases and 7845 sub-cohorts of Japanese females during a median follow-up of 16.9 years). RESULTS: For the best PRS model involving 46,861 single nucleotide polymorphisms (SNPs; P + T method with PT = 0.05 and R2 = 0.2) derived from Japanese-ancestry GWAS, the Harrell's C-statistic was 0.598 ± 0.018 in the evaluation dataset. The age-adjusted hazard ratio for breast cancer in females with the highest PRS quintile compared with those in the lowest PRS quintile was 2.47 (95% confidence intervals, 1.64-3.70). The PRS constructed using Japanese-ancestry GWAS demonstrated better predictive performance for breast cancer in Japanese females than that using European-ancestry GWAS (Harrell's C-statistics 0.598 versus 0.586). CONCLUSION: This study developed a breast cancer PRS for Japanese females and demonstrated the usefulness of the PRS for breast cancer risk stratification.

Body mass index and colorectal cancer risk: A Mendelian randomization study.

Traditional observational studies have reported a positive association between higher body mass index (BMI) and the risk of colorectal cancer (CRC). However, evidence from other approaches to pursue the causal relationship between BMI and CRC is sparse. A two-sample Mendelian randomization (MR) study was undertaken using 68 single nucleotide polymorphisms (SNPs) from the Japanese genome-wide association study (GWAS) and 654 SNPs from the GWAS catalogue for BMI as sets of instrumental variables. For the analysis of SNP-BMI associations, we undertook a meta-analysis with 36 303 participants in the Japanese Consortium of Genetic Epidemiology studies (J-CGE), comprising normal populations. For the analysis of SNP-CRC associations, we utilized 7636 CRC cases and 37 141 controls from five studies in Japan, and undertook a meta-analysis. Mendelian randomization analysis of inverse-variance weighted method indicated that a one-unit (kg/m2 ) increase in genetically predicted BMI was associated with an odds ratio of 1.13 (95% confidence interval, 1.06-1.20; P value <.001) for CRC using the set of 68 SNPs, and an odds ratio of 1.07 (1.03-1.11, 0.001) for CRC using the set of 654 SNPs. Sensitivity analyses robustly showed increased odds ratios for CRC for every one-unit increase in genetically predicted BMI. Our MR analyses strongly support the evidence that higher BMI influences the risk of CRC. Although Asians are generally leaner than Europeans and North Americans, avoiding higher BMI seems to be important for the prevention of CRC in Asian populations.

Frequent post-operative monitoring of colorectal cancer using individualised ctDNA validated by multiregional molecular profiling.

BACKGROUND: Circulating tumour DNA (ctDNA) is known as a tumour-specific personalised biomarker, but the mutation-selection criteria from heterogeneous tumours remain a challenge. METHODS: We conducted multiregional sequencing of 42 specimens from 14 colorectal tumours of 12 patients, including two double-cancer cases, to identify mutational heterogeneity to develop personalised ctDNA assays using 175 plasma samples. RESULTS: "Founder" mutations, defined as a mutation that is present in all regions of the tumour in a binary manner (i.e., present or absent), were identified in 12/14 tumours. In contrast, "truncal" mutations, which are the first mutation that occurs prior to the divergence of branches in the phylogenetic tree using variant allele frequency (VAF) as continuous variables, were identified in 12/14 tumours. Two tumours without founder and truncal mutations were hypermutators. Most founder and truncal mutations exhibited higher VAFs than "non-founder" and "branch" mutations, resulting in a high chance to be detected in ctDNA. In post-operative long-term observation for 10/12 patients, early relapse prediction, treatment efficacy and non-relapse corroboration were achievable from frequent ctDNA monitoring. CONCLUSIONS: A single biopsy is sufficient to develop custom dPCR probes for monitoring tumour burden in most CRC patients. However, it may not be effective for those with hypermutated tumours.

Evaluation of clinical formalin-fixed paraffin-embedded tissue quality for targeted-bisulfite sequencing.

Formalin-fixed paraffin-embedded (FFPE) tissues are promising biological resources for genetic research. Recent improvements in DNA extraction from FFPE samples allowed the use of these tissues for multiple sequencing methods. However, fundamental research addressing the application of FFPE-derived DNA for targeted-bisulfite sequencing (TB-seq) is lacking. Here, we evaluated the suitability of FFPE-derived DNA for TB-seq. We conducted TB-seq using FFPE-derived DNA and corresponding fresh frozen (FF) tissues of patients with kidney cancer and compared the quality of DNA, libraries, and TB-seq statistics between the two preservation methods. The approximately 600-bp average fragment size of the FFPE-derived DNA was significantly shorter than that of the FF-derived DNA. The sequencing libraries constructed using FFPE-derived DNA and the mapping ratio were approximately 10 times and 10% lower, respectively, than those constructed using FF-derived DNA. In the mapped data of FFPE-derived DNA, duplicated reads accounted for > 60% of the obtained sequence reads, with lower mean on-target coverage. Therefore, the standard TB-seq protocol is inadequate for obtaining high-quality data for epigenetic analysis from FFPE-derived DNA, and technical improvements are necessary for enabling the use of archived FFPE resources.

The interaction between ABCA1 polymorphism and physical activity on the HDL-cholesterol levels in a Japanese population.

Few studies have investigated the interactions between HDL-C-related SNPs identified by genome-wide association (GWA) study and physical activity (PA) on HDL-C. First, we conducted a sex-stratified GWA study in a discovery sample (2,231 men and 2,431 women) and replication sample (2,599 men and 3,109 women) to identify SNPs influencing log-transformed HDL-C in Japanese participants in the baseline survey of the Japan Multi-Institutional Collaborative Cohort Study. We also replicated previously reported HDL-C-related SNPs in a combined (discovery plus replication) sample (4,830 men and 5,540 women). We then analyzed the interactions of the HDL-C-related SNPs with PA on HDL-C. The sex-stratified GWA analyses identified 11 and 10 HDL-C-related SNPs in men and women as targets for an interaction analysis. Among these, only one interaction of ABCA1 rs1883025 with PA was statistically significant in men, after Bonferroni correction [P-interaction = 0.001 (α = 0.05/21 = 0.002)]. The per-major-allele (C allele) increase in log-transformed HDL-C was lost in men with low PA (ß = 0.008) compared with those with medium (ß = 0.032) or high PA (ß = 0.034). These findings suggest that the benefit of carrying a C allele of ABCA1 rs1883025 on enhancing HDL-C may be attenuated in inactive men.

Genome-Wide Polygenic Score and the Risk of Ischemic Stroke in a Prospective Cohort: The Hisayama Study.

Background and Purpose- Environmental and genetic factors contribute to the development of ischemic stroke (IS). We recently developed a genome-wide polygenic risk score (PRS) for IS using case-control datasets from 4 large-scale observational studies conducted in Japan. Our objective in the present study was to confirm the association between the PRS and the risk of IS with data from an independent prospective cohort recruited from the general Japanese population. Methods- A total of 3038 subjects aged ≥40 years were followed up for 10 years (2002-2012). The genome-wide PRS was calculated using genotype data from >350 000 single-nucleotide polymorphisms. The PRS levels were divided into quintiles. High and low genetic risk groups were defined as top 60% and bottom 40% of PRS, respectively. The hazard ratio (HR) for the development of IS was estimated using a Cox proportional hazards model. Results- During the follow-up period, 91 cases developed first-ever IS. The age- and sex-adjusted HR for IS increased with higher PRS levels (P for trend, 0.03). Subjects with the highest quintile level of PRS had a 2.44-fold (95% CI, 1.16-5.12) greater risk for IS than those with the lowest quintile level after adjusting for age and sex. A similar association was observed after adjusting for environmental risk factors (P for trend, 0.03). As compared with low genetic risk group, the age- and sex-adjusted HR in high genetic risk group was 1.63 (95% CI, 1.04-2.55), which was comparable to the HR of hypertension (HR, 1.41), diabetes mellitus (HR, 1.72), and smoking (HR, 1.54). The age- and sex-adjusted HR increased with the number of environmental risk factors in both high and low genetic risk groups without significant interaction. Conclusions- A high genome-wide PRS was a significant risk factor for IS independent of environmental risk factors in a general Japanese population. This finding suggests that PRS may be useful to identify individuals at a high risk of IS.

Strong association between the 12q24 locus and sweet taste preference in the Japanese population revealed by genome-wide meta-analysis.

The sweet taste preference of humans is an important adaptation to ensure the acquisition of carbohydrate nutrition; however, overconsumption of sweet foods can potentially lead to diseases such as obesity and diabetes. Although previous studies have suggested that interindividual variation of human sweet taste preference is heritable, genetic loci associated with the trait have yet to be fully elucidated. Here, we genotyped 12,312 Japanese participants using the HumanCore-12+ Custom BeadChip or the HumanCore-24 Custom BeadChip microarrays. The sweet taste preference of the participants was surveyed via an internet-based questionnaire, resulting in a five-point scale of sweet taste preference. The genome-wide meta-analysis of the Japanese participants revealed a strong association between the 12q24 locus and sweet taste preference scale (P = 2.8 × 10-70). The lead variant rs671 is monoallelic in non-East Asian populations and is located in the aldehyde dehydrogenase (ALDH2) gene, encoding an enzyme involved in alcohol metabolism. The association between the minor allele of rs671 and sweet taste preference was attenuated by adjusting for alcohol drinking. The subgroup analysis showed that the effect of rs671 on sweet taste preference was greater in males than in females. In conclusion, we found an association between the 12q24 locus and sweet taste preference in the Japanese population, and showed that the adjustment for drinking habits attenuated the association. This novel genetic association may provide new clues to elucidate mechanisms determining sweet taste preferences.

GWAS of habitual coffee consumption reveals a sex difference in the genetic effect of the 12q24 locus in the Japanese population.

BACKGROUND: Studies on genetic effects of coffee consumption are scarce for Asian populations. We conducted a genome-wide association study (GWAS) of habitual coffee consumption in Japan using a self-reporting online survey. RESULTS: Candidate genetic loci associated with habitual coffee consumption were searched within a discovery cohort (N = 6,264) and confirmed in a replication cohort (N = 5,975). Two loci achieved genome-wide significance (P < 5 × 10- 8) in a meta-analysis of the discovery and replication cohorts: an Asian population-specific 12q24 (rs79105258; P = 9.5 × 10- 15), which harbors CUX2, and 7p21 (rs10252701; P = 1.0 × 10- 14), in the upstream region of the aryl hydrocarbon receptor (AHR) gene, involved in caffeine metabolism. Subgroup analysis revealed a stronger genetic effect of the 12q24 locus in males (P for interaction = 8.2 × 10- 5). Further, rs79105258 at the 12q24 locus exerted pleiotropic effects on body mass index (P = 3.5 × 10- 4) and serum triglyceride levels (P = 8.7 × 10- 3). CONCLUSIONS: Our results consolidate the association of habitual coffee consumption with the 12q24 and 7p21 loci. The different effects of the 12q24 locus between males and females are a novel finding that improves our understanding of genetic influences on habitual coffee consumption.

Genetic Predisposition to Ischemic Stroke: A Polygenic Risk Score.

BACKGROUND AND PURPOSE: The prediction of genetic predispositions to ischemic stroke (IS) may allow the identification of individuals at elevated risk and thereby prevent IS in clinical practice. Previously developed weighted multilocus genetic risk scores showed limited predictive ability for IS. Here, we investigated the predictive ability of a newer method, polygenic risk score (polyGRS), based on the idea that a few strong signals, as well as several weaker signals, can be collectively informative to determine IS risk. METHODS: We genotyped 13 214 Japanese individuals with IS and 26 470 controls (derivation samples) and generated both multilocus genetic risk scores and polyGRS, using the same derivation data set. The predictive abilities of each scoring system were then assessed using 2 independent sets of Japanese samples (KyushuU and JPJM data sets). RESULTS: In both validation data sets, polyGRS was shown to be significantly associated with IS, but weighted multilocus genetic risk scores was not. Comparing the highest with the lowest polyGRS quintile, the odds ratios for IS were 1.75 (95% confidence interval, 1.33-2.31) and 1.99 (95% confidence interval, 1.19-3.33) in the KyushuU and JPJM samples, respectively. Using the KyushuU samples, the addition of polyGRS to a nongenetic risk model resulted in a significant improvement of the predictive ability (net reclassification improvement=0.151; P<0.001). CONCLUSIONS: The polyGRS was shown to be superior to weighted multilocus genetic risk scores as an IS prediction model. Thus, together with the nongenetic risk factors, polyGRS will provide valuable information for individual risk assessment and management of modifiable risk factors.

f-treeGC: a questionnaire-based family tree-creation software for genetic counseling and genome cohort studies.

BACKGROUND: The Tohoku Medical Megabank project aims to create a next-generation personalized healthcare system by conducting large-scale genome-cohort studies involving three generations of local residents in the areas affected by the Great East Japan Earthquake. We collected medical and genomic information for developing a biobank to be used for this healthcare system. We designed a questionnaire-based pedigree-creation software program named "f-treeGC," which enables even less experienced medical practitioners to accurately and rapidly collect family health history and create pedigree charts. RESULTS: f-treeGC may be run on Adobe AIR. Pedigree charts are created in the following manner: 1) At system startup, the client is prompted to provide required information on the presence or absence of children; f-treeGC is capable of creating a pedigree up to three generations. 2) An interviewer fills out a multiple-choice questionnaire on genealogical information. 3) The information requested includes name, age, gender, general status, infertility status, pregnancy status, fetal status, and physical features or health conditions of individuals over three generations. In addition, information regarding the client and the proband, and birth order information, including multiple gestation, custody, multiple individuals, donor or surrogate, adoption, and consanguinity may be included. 4) f-treeGC shows only marriages between first cousins via the overlay function. 5) f-treeGC automatically creates a pedigree chart, and the chart-creation process is visible for inspection on the screen in real time. 6) The genealogical data may be saved as a file in the original format. The created/modified date and time may be changed as required, and the file may be password-protected and/or saved in read-only format. To enable sorting or searching from the database, the file name automatically contains the terms typed into the entry fields, including physical features or health conditions, by default. 7) Alternatively, family histories are collected using a completed foldable interview paper sheet named "f-sheet", which is identical to the questionnaire in f-treeGC. CONCLUSIONS: We developed a questionnaire-based family tree-creation software, named f-treeGC, which is fully compliant with international recommendations for standardized human pedigree nomenclature. The present software simplifies the process of collecting family histories and pedigrees, and has a variety of uses, from genome cohort studies or primary care to genetic counseling.

Population-based biobank participants' preferences for receiving genetic test results.

There are ongoing debates on issues relating to returning individual research results (IRRs) and incidental findings (IFs) generated by genetic research in population-based biobanks. To understand how to appropriately return genetic results from biobank studies, we surveyed preferences for returning IRRs and IFs among participants of the Tohoku Medical Megabank Project (TMM). We mailed a questionnaire to individuals enrolled in the TMM cohort study (Group 1; n=1031) and a group of Tohoku region residents (Group 2; n=2314). The respondents were required to be over 20 years of age. Nearly 90% of Group 1 participants and over 80% of Group 2 participants expressed a preference for receiving their genetic test results. Furthermore, over 60% of both groups preferred to receive their genetic results 'from a genetic specialist.' A logistic regression analysis revealed that engaging in 'health-conscious behaviors' (such as regular physical activity, having a healthy diet, intentionally reducing alcohol intake and/or smoking and so on) was significant, positively associated with preferring to receive their genetic test results (odds ratio=2.397 (Group 1) and 1.897 (Group 2)). Our findings provided useful information and predictors regarding the return of IRRs and IFs in a population-based biobank.

MetaVelvet: an extension of Velvet assembler to de novo metagenome assembly from short sequence reads.

An important step in 'metagenomics' analysis is the assembly of multiple genomes from mixed sequence reads of multiple species in a microbial community. Most conventional pipelines use a single-genome assembler with carefully optimized parameters. A limitation of a single-genome assembler for de novo metagenome assembly is that sequences of highly abundant species are likely misidentified as repeats in a single genome, resulting in a number of small fragmented scaffolds. We extended a single-genome assembler for short reads, known as 'Velvet', to metagenome assembly, which we called 'MetaVelvet', for mixed short reads of multiple species. Our fundamental concept was to first decompose a de Bruijn graph constructed from mixed short reads into individual sub-graphs, and second, to build scaffolds based on each decomposed de Bruijn sub-graph as an isolate species genome. We made use of two features, the coverage (abundance) difference and graph connectivity, for the decomposition of the de Bruijn graph. For simulated datasets, MetaVelvet succeeded in generating significantly higher N50 scores than any single-genome assemblers. MetaVelvet also reconstructed relatively low-coverage genome sequences as scaffolds. On real datasets of human gut microbial read data, MetaVelvet produced longer scaffolds and increased the number of predicted genes.

Evaluating Leukocyte Telomere Length and Myeloid-Derived Suppressor Cells as Biomarkers for Prostate Cancer.

BACKGROUND: Leukocyte telomere length (LTL) and myeloid-derived suppressor cells (MDSC) are associated with aging and the development and progression of cancer. However, the exact nature of this relationship remains unclear. Our study aimed to investigate the potential of LTL and MDSC as diagnostic biomarkers for prostate cancer while also seeking to deepen our understanding of the relationship of these potential biomarkers to each other. METHODS: Our study involved patients undergoing a prostate biopsy. We analyzed the relative LTL in genomic DNA obtained from peripheral blood leukocytes as well as the percentage of MDSC and their subtypes in peripheral blood mononuclear cells (PBMC). Our evaluation focused on examining the relationship between LTL and MDSC and pathological diagnoses as well as investigating the correlation between LTL and MDSC levels. RESULTS: In our study of 102 participants, 56 were pathologically diagnosed with localized prostate cancer (cancer group), while 46 tested negative (control group). The cancer group exhibited significantly shorter LTL in comparison to the control group (p = 0.024). Additionally, the cancer group showed a tendency towards a higher percentage of monocytic MDSC (M-MDSC), although this difference did not reach statistical significance (p = 0.056). Our multivariate logistic regression analysis revealed that patients with shorter LTL and higher percentages of M-MDSC had a 2.98-fold (95% CI = 1.001-8.869, p = 0.049) and 3.03-fold (95% CI = 1.152-7.977, p = 0.025) increased risk of prostate cancer diagnosis, respectively. There was also a significant negative correlation between LTL and M-MDSC. (r = -0.347, p < 0.001). CONCLUSIONS: Our research has established a correlation between LTL and MDSC in patients undergoing biopsy for prostate cancer. Notably, we observed that individuals with localized prostate cancer tend to have shorter LTL and a higher percentage of M-MDSC prior to their diagnosis. These findings suggest that LTL and M-MDSC could potentially serve as adjunctive biomarkers for the early diagnosis of prostate cancer.

Association of genetically determined chronotype with circulating testosterone: a Mendelian randomization study.

Low testosterone levels in men have been linked to decreased physical and mental function, as well as a reduced quality of life. Previous prospective observational studies have suggested an association between testosterone and sleep traits, but the causality of this relationship remains unclear. We aimed to explore the potential causal link between genetically determined sleep traits and testosterone levels in men using Mendelian randomization (MR) analysis from the UK Biobank dataset. Our exposures were genetic variants associated with sleep traits (chronotype and sleep duration), whereas our outcomes were traits of sex steroid hormones (total testosterone, TT; bioavailable testosterone, BAT; and sex hormone-binding globulin, SHBG). We employed inverse variance weighted (IVW) and weighted median (WM) methods to assess the causal associations. The IVW method offers a robust estimate of causality, whereas the WM method provides reliable results even when some genetic variants are invalid instruments. Our main analysis involving sex steroid hormones and chronotype identified 155 chronotype-related variants. The primary findings from the analysis, which used chronotype as the exposure and sex steroid hormones as the outcomes, showed that a genetically predicted chronotype score was significantly associated with an increased levels of TT (association coefficient ß, 0.08; 95% confidence interval [CI], 0.02-0.14; P = 0.008) and BAT (ß, 0.08; 95% CI, 0.02-0.14; P = 0.007), whereas there was no significant association with SHBG (ß, 0.01; 95% CI, -0.02-0.03; P = 0.64). Meanwhile, MR analysis of sex steroid hormones and sleep duration was performed, and 69 variants associated with sleep duration were extracted. There were no significant association between sleep duration and sex steroid hormones (TT, P = 0.91; BAT, P = 0.82; and SHBG, P = 0.95). Our data support a causal association between chronotype and circulating testosterone levels in men. These findings underscore a potential causal relationship between chronotype and testosterone levels in men, suggesting that lifestyle adjustments are crucial for men's health. Recognizing factors that influence testosterone is essential. One limitation of this study is the use of one-sample MR, which can introduce potential bias due to non-independence of genetic associations for exposure and outcome. In conclusion, our findings indicate that a morning preference is correlated with circulating testosterone levels, emphasizing the potential impact of lifestyle habits on testosterone levels in men.

Dimension reduction of microbiome data linked Bifidobacterium and Prevotella to allergic rhinitis.

Dimension reduction has been used to visualise the distribution of multidimensional microbiome data, but the composite variables calculated by the dimension reduction methods have not been widely used to investigate the relationship of the human gut microbiome with lifestyle and disease. In the present study, we applied several dimension reduction methods, including principal component analysis, principal coordinate analysis (PCoA), non-metric multidimensional scaling (NMDS), and non-negative matrix factorization, to a microbiome dataset from 186 subjects with symptoms of  allergic rhinitis (AR) and 106 controls. All the dimension reduction methods supported that the distribution of microbial data points appeared to be continuous rather than discrete. Comparison of the composite variables calculated from the different dimension reduction methods showed that the characteristics of the composite variables differed depending on the distance matrices and the dimension reduction methods. The first composite variables calculated from PCoA and NMDS with the UniFrac distance were strongly associated with AR (FDR adjusted P = 2.4 × 10-4 for PCoA and P = 2.8 × 10-4 for NMDS), and also with the relative abundance of Bifidobacterium and Prevotella. The abundance of Bifidobacterium was also linked to intake of several nutrients, including carbohydrate, saturated fat, and alcohol via composite variables. Notably, the association between the composite variables and AR was much stronger than the association between the relative abundance of individual genera and AR. Our results highlight the usefulness of the dimension reduction methods for investigating the association of microbial composition with lifestyle and disease in clinical research.

Increased circulating polymorphonuclear myeloid-derived suppressor cells are associated with prognosis of metastatic castration-resistant prostate cancer.

Introduction: Myeloid-derived suppressor cell (MDSC) exhibits immunosuppressive functions and affects cancer progression, but its relationship with prostate cancer remains unclear. We elucidated the association of polymorphonuclear MDSC (PMN-MDSC) and monocytic MDSC (M-MDSC) levels of the total peripheral blood mononuclear cells (PBMCs) with prostate cancer progression and evaluated their roles as prognostic indicators. Methods: We enrolled 115 patients with non-metastatic hormone-sensitive prostate cancer (nmHSPC, n = 62), metastatic hormone-sensitive prostate cancer (mHSPC, n = 23), and metastatic castration-resistant prostate cancer (mCRPC, n = 30). Subsequently, the proportions of MDSCs in each disease progression were compared. Log-rank tests and multivariate Cox regression analyses were performed to ascertain the associations of overall survival. Results: The patients with mCRPC had significantly higher PMN-MDSC percentage than those with nmHSPC and mHSPC (P = 7.73 × 10-5 and 0.0014). Significantly elevated M-MDSC levels were observed in mCRPC patients aged <70 years (P = 0.016) and with a body mass index (BMI) <25 kg/m2 (P = 0.043). The high PMN-MDSC group had notably shorter median survival duration (159 days) than the low PMN-MDSC group (768 days, log-rank P = 0.018). In the multivariate analysis including age, BMI, and MDSC subset, PMN-MDSC was significantly associated with prognosis (hazard ratios, 3.48; 95% confidence interval: 1.05-11.56, P = 0.042). Discussion: PMN-MDSC levels are significantly associated with mCRPC prognosis. Additionally, we highlight the remarkable associations of age and BMI with M-MDSC levels in mCRPC, offering novel insights into MDSC dynamics in prostate cancer progression.

COPICAT: a software system for predicting interactions between proteins and chemical compounds.

UNLABELLED: Since tens of millions of chemical compounds have been accumulated in public chemical databases, fast comprehensive computational methods to predict interactions between chemical compounds and proteins are needed for virtual screening of lead compounds. Previously, we proposed a novel method for predicting protein-chemical interactions using two-layer Support Vector Machine classifiers that require only readily available biochemical data, i.e. amino acid sequences of proteins and structure formulas of chemical compounds. In this article, the method has been implemented as the COPICAT web service, with an easy-to-use front-end interface. Users can simply submit a protein-chemical interaction prediction job using a pre-trained classifier, or can even train their own classification model by uploading training data. COPICAT's fast and accurate computational prediction has enhanced lead compound discovery against a database of tens of millions of chemical compounds, implying that the search space for drug discovery is extended by >1000 times compared with currently well-used high-throughput screening methodologies. AVAILABILITY: The COPICAT server is available at http://copicat.dna.bio.keio.ac.jp. All functions, including the prediction function are freely available via anonymous login without registration. Registered users, however, can use the system more intensively.