OP0110 CELL-TYPE-SPECIFIC TRANSCRIPTOME ARCHITECTURE UNDERLYING THE ESTABLISHMENT AND EXACERBATION OF SYSTEMIC LUPUS ERYTHEMATOSUS

Background

Systemic lupus erythematosus (SLE) is a complex autoimmune disease with unknown etiology involving multiple immune cells and has diverse clinical phenotypes. This heterogeneous nature has hampered a better understanding of SLE pathogenesis and the development of effective therapeutic agents. While recent single-cell RNA sequencing studies of SLE identified several important cell subpopulations, they were limited by sparse expression information at single-cell level and small sample sizes.

Objectives

This study aimed to elucidate the dysregulated gene expression pattern linked to multiple clinical statuses of SLE with a fine cellular resolution and higher sensitivity. We also attempted to resolve a complex interaction between risk variants and the transcriptome dysregulation seen in SLE patients.

Methods

We conducted a large-scale bulk transcriptome study of 6,386 RNA-sequencing data including 27 purified immune cell types in peripheral blood from 136 SLE and 89 healthy donors in the Immune Cell Gene Expression Atlas from the University of Tokyo (ImmuNexUT) cohort1. At enrollment, SLE patients had diverse clinical manifestations (disease activity, organ involvement and treatment profiles) and 22 patients were re-evaluated after belimumab treatment.

Results

We first profiled two distinct cell-type-specific transcriptomic signatures: disease-state and disease-activity signatures, reflecting disease establishment and exacerbation, respectively.

After confirming the high replicability of both signatures in independent cohorts, we identified candidates of biological processes unique to each signature: e.g., upregulated E2F transcriptional activity in Th1, CD8+ memory T-lineage and NK cells, and dynamic increase of IL21 and CXCL13 in Th1 cells in an active phase of SLE. Pathway analysis highlighted the importance of immunometabolic process for SLE (e.g., oxidative phosphorylation) in cell-type-specific resolution.

Moreover, we demonstrated cell-type-specific contributions to diverse organ involvement, e.g., Th1 for mucocutaneous, monocyte-lineage cells for musculoskeletal, neutrophil-lineage cells for renal activity, respectively.

We also observed the strong associations of disease-activity signatures with treatment effect: (i) belimumab suppressed activity signatures from B-lineage cells, especially in good responders and (ii) mycophenolate mofetil substantially suppressed activity signatures from plasmablast, Th1, and central memory CD8 cells.

However, through stratified LD score regression using large-scale SLE-GWASs, we revealed that disease-activity signatures were less enriched around SLE risk variants than disease-state signatures. Consistent with this result, the directions of SLE risk alleles’ expression quantitative trait locus (eQTL) effects were significantly concordant with the directions of disease-state signatures, but not with those of activity signatures. These findings suggested that the current genetic case-control studies may not well capture clinically vital biology linked to drug target discovery for SLE. Meanwhile, we also detected some examples of activity signatures that might contribute to the disease risk by modulating risk allele’s eQTL effects.

Figure 1.

Conclusion

We identified comprehensive gene signatures reflecting the establishment and exacerbation of SLE, which provide essential foundations for future genomic, genetic, and clinical studies.

References

[1]Ota, M. et al. Dynamic landscape of immune cell-specific gene regulation in immune-mediated diseases. Cell 2021;184:3006-21.e17.

Acknowledgements

This study was supported by Chugai Pharmaceutical Co., Ltd., Tokyo, Japan; the Ministry of Education, Culture, Sports; and the Japan Agency for Medical Research and Development (AMED) (JP21tm0424221 and JP21zf0127004).

Disclosure of Interests

Masahiro Nakano: None declared, Mineto Ota Grant/research support from: Mineto Ota belongs to the Social Cooperation Program, Department of functional genomics and immunological diseases, supported by Chugai Pharmaceutical., Yusuke Takeshima Grant/research support from: Yusuke Takeshima belonged to the Social Cooperation Program, Department of functional genomics and immunological diseases, supported by Chugai Pharmaceutical., Yukiko Iwasaki: None declared, Hiroaki Hatano: None declared, Yasuo Nagafuchi Grant/research support from: Yasuo Nagafuchi belongs to the Social Cooperation Program, Department of functional genomics and immunological diseases, supported by Chugai Pharmaceutical., Kwangwoo Kim: None declared, So-Young Bang: None declared, Hye Soon Lee: None declared, Hirofumi Shoda: None declared, Xuejun Zhang: None declared, Sang-Cheol Bae: None declared, Chikashi Terao: None declared, Kazuhiko Yamamoto: None declared, Tomohisa Okamura Grant/research support from: Tomohisa Okamura belongs to the Social Cooperation Program, Department of functional genomics and immunological diseases, supported by Chugai Pharmaceutical., Kazuyoshi Ishigaki: None declared, Keishi Fujio Speakers bureau: Keishi Fujio receives speaker fees from Chugai Pharmaceutical., Consultant of: Keishi Fujio receives consulting honoraria from Chugai Pharmaceutical., Grant/research support from: Keishi Fujio receives research support from Chugai Pharmaceutical.

Development of a high-resolution typing method for SLA-3, swine MHC class I antigen 3

We developed a high-resolution and comprehensive typing method for swine leukocyte antigen 3 (SLA-3), an MHC class I gene, employing locus-specific genomic PCR followed by subsequent direct sequencing. A total of 292 individuals from nine pure, one cross-breed and six cell lines were successfully typed. A total of 21 SLA-3 alleles were identified, of which four were found to be novel alleles. However, the allelic diversity of SLA-3 was lower than that of previously reported class I genes, SLA-1 and -2. More SLA-3 alleles were observed in the Landrace and Yorkshire breeds than the other breeds. SLA-3*04:01 was identified in seven out of nine breeds and was the most widely distributed allele across all breeds. Therefore, the typing method reported in this study completes our efforts to develop high-resolution typing methods for major SLA molecules, facilitating the combined analysis of major SLA genes from field samples, which is important to understand the relationship between the adaptive immune responses against pathogens and the immunogenetic makeup of an individual.

A population-based study for difference in contribution of type 2 diabetes mellitus to cardiovascular outcomes in adult cancer survivors and general population without history of cancer

Background/Introduction

Type 2 diabetes mellitus (T2DM) is a well-established risk factor for cardiovascular disease (CVD). However, evidence is scarce on the difference in burden of T2DM as a cardiovascular risk factor in adult cancer survivors compared to the general population without history of cancer. This evidence gap needs to be addressed as increased risk for CVD in adult cancer survivors compared to the general population has been consistently reported from observational studies.

Purpose

To investigate the association of T2DM with CVD in adult cancer survivors and general population without history of cancer and compare the magnitude of associations for each population using propensity score-matched analyses

Methods

For this population-based cohort study, we used the National Health Insurance Service-National Sample Cohort (NHIS-NSC) established in the Republic of Korea (2002–2015) to identity adult cancer survivors who were aged 18 years or older and survived more than 12 months after the first-ever cancer diagnosis and general population without history of cancer matched for age and sex in a 1:1 ratio. In both cohorts, those with history of CVD prior to the index date were excluded. We used Cox proportional hazards model to compute hazard ratios (HR) and 95% confidence intervals (95% CI) for incident CVD in those with T2DM compared to those without T2DM in both adult cancer survivors and the general population after adjusting for shared confounding factors (income status, cigarette smoking, alcohol consumption, physical inactivity, body mass index, blood pressure, lipid abnormalities, family history of CVD, and aspirin use). Furthermore, we used Q statistic to compare the magnitude of association between T2DM and CVD in this matched cohort.

Results

The overall age and sex matched cohort in the NHIS-NSC comprised of 5,163 adult cancer survivors and the equivalent number of participants categorized as general population without history of cancer. The adjusted HR for CVD among adult cancer survivors with T2DM, as compared to without T2DM was 2.27 (95% CI: 1.25–4.11). Similarly, T2DM was associated with an increased risk of CVD in the general population without cancer (HR=1.88; 95% CI: 1.03–3.42). However, no strong statistical evidence was found for difference in contribution of T2DM to CVD risk between adult cancer survivors and general population without history of cancer (Pheterogeneity = 0.662).

Conclusions

In this population-based cohort study, T2DM was associated with an increased risk for CVD in both adult cancer survivors and those without history of cancer. Difference in magnitude of T2DM-CVD associations between adult cancer survivors and the general population without history of cancer was not substantial from a statistical standpoint. However, adult cancer survivor with T2DM may still need additional clinical attention for CVD prevention compared to those without history of cancer.

Funding Acknowledgement

Type of funding sources: Public grant(s) – National budget only. Main funding source(s): National Cancer Center (Republic of Korea)

Treatment potential of LPCN 1144 on liver health and metabolic regulation in a non-genomic, high fat diet induced NASH rabbit model

PURPOSE

Low free testosterone (T) level in men is independently associated with presence and severity of Non-Alcoholic Steatohepatitis (NASH). The histological and molecular effects of oral testosterone prodrug LPCN 1144 treatment on hepatic fibrosis and NASH features are unknown. A metabolic syndrome-induced NASH model in rabbits consuming high fat diet (HFD) has been previously used to assess treatment effects of injectable T on hepatic fibrosis and NASH features. Here we present results on LPCN 1144 in this HFD-induced, NASH preclinical model.

METHODS

Male rabbits were randomly assigned to five groups: regular diet (RD), HFD, HFD + 1144 vehicle (HFD + Veh), HFD + 1144 (1144), and HFD + 1144 + α-tocopherol (1144 + ALPHA). Rabbits were sacrificed after 12 weeks for liver histological, biochemical and genetic analyses. Histological scores were obtained through Giemsa (inflammation), Masson's trichrome (steatosis and ballooning), and Picrosirius Red (fibrosis) staining.

RESULTS

Compared to RD, HFD and HFD + Veh significantly worsened NASH features and hepatic fibrosis. Considering HFD and HFD + Veh arms, histological and biomarker features were not significantly different. Both 1144 and 1144 + ALPHA arms improved mean histological scores of NASH as compared to HFD arm. Importantly, percentage of fibrosis was improved in both 1144 (p < 0.05) and 1144 + ALPHA (p = 0.05) treatment arms vs. HFD. Both treatment arms also reduced HFD-induced inflammation and fibrosis mRNA markers. Furthermore, 1144 treatments significantly improved HFD-induced metabolic dysfunctions.

CONCLUSIONS

Histological and biomarker analyses demonstrate that LPCN 1144 improved HFD-induced hepatic fibrosis and NASH biochemical, biomolecular and histochemical features. These preclinical findings support a therapeutic potential of LPCN 1144 in the treatment of NASH and of hepatic fibrosis.

Spontaneous Formation of Star-Shaped Surface Patterns in a Driven Bose-Einstein Condensate

We observe experimentally the spontaneous formation of star-shaped surface patterns in driven Bose-Einstein condensates. Two-dimensional star-shaped patterns with l-fold symmetry, ranging from quadrupole (l=2) to heptagon modes (l=7), are parametrically excited by modulating the scattering length near the Feshbach resonance. An effective Mathieu equation and Floquet analysis are utilized, relating the instability conditions to the dispersion of the surface modes in a trapped superfluid. Identifying the resonant frequencies of the patterns, we precisely measure the dispersion relation of the collective excitations. The oscillation amplitude of the surface excitations increases exponentially during the modulation. We find that only the l=6 mode is unstable due to its emergent coupling with the dipole motion of the cloud. Our experimental results are in excellent agreement with the mean-field framework. Our work opens a new pathway for generating higher-lying collective excitations with applications, such as the probing of exotic properties of quantum fluids and providing a generation mechanism of quantum turbulence.

Change in the Shrinkage Forces of Composite Resins According to Controlled Deflection

The aim of this study was to investigate how the polymerization shrinkage forces of composite resins change with change in deflection. Five composites, SDR (Dentsply Caulk, Milford, DE, USA), EcuSphere-Shape (DMG, Hamburg, Germany), Tetric N-Ceram Bulk Fill (Ivoclar Vivadent, Schaan, Liechtenstein), CLEARFIL AP-X (Kuraray Noritake Dental Inc., Sakazu, Kurashiki, Okayama, Japan), and Filtek Z350 XT (3M Dental Products, St Paul, MN, USA), were tested in this experiment. The polymerization shrinkage forces of the composites were measured using a custom-made tooth-deflection-mimicking device and software (R&B Inc, Daejon, Korea). In all measurements, six modes were tested: maximum-deflection, zero-deflection, and four deflection-controlled modes. For each deflection mode, the shrinkage forces were recorded continuously every 0.1 second for 180 seconds. Polymerization shrinkage and flexural modulus were also measured. Eight specimens of each material were allocated for each test. For each material, six groups of shrinkage force values were compared using one-way ANOVA and Tukey post hoc tests at a 95% confidence level. The polymerization shrinkage force of each material in each of the six deflection modes was analyzed with 95% confidence using one-way ANOVA and Tukey post hoc tests. The relationship between the force measured in the six deflection modes and the linear polymerization shrinkage and flexural modulus was analysed with 95% confidence using Pearson correlation analysis. For each material, the following held true: the shrinkage force was highest in zero-deflection mode, the force decreased as deflection increased, and the smallest force appeared in maximum-deflection mode (p<0.05). There was a high negative correlation between allowable deflection and shrinkage force in all materials.

Is dual-energy absorptiometry accurate in the assessment of bone status of patients with chronic kidney disease?

Several patients with chronic kidney disease (CKD) have deteriorated bone status. Estimation of bone status using DXA has limitations especially in patients with CKD accompanying aortic calcifications. Quantitative CT and the trabecular bone score could be more accurate methods to estimate bone status for patients with CKD and vascular calcifications.

INTRODUCTION

It remains unclear whether dual-energy absorptiometry (DXA) is appropriate for the assessment of bone status in patients with chronic kidney disease (CKD), a disease that impacts bone health. The aims of this study were to compare DXA and central quantitative computed tomography (cQCT) and to evaluate bone status in patients with pre-dialysis CKD.

METHODS

This retrospective study included 363 healthy control subjects whose bone mineral density (BMD) was evaluated with DXA and 117 CKD patients whose BMD was evaluated using both cQCT and DXA. Diagnostic discordance was assessed between the lumbar spine (LS) and femur neck (FN) from DXA or between two modalities. The trabecular bone score (TBS) was extracted from DXA images. The volume of abdominal aortic calcification (AAC) was calculated using CT images from cQCT.

RESULTS

Using LS DXA T-score, osteoporosis was less common in the CKD group than in controls. Patients with normal LS BMD using DXA were reclassified into osteopenia or osteoporosis using cQCT in CKD patients. Among discordant subjects between FN and LS in DXA, a higher BMD of LS was more common in CKD patients than in controls. CKD patients had lower TBS than controls despite having the same diagnosis using DXA. AAC volume negatively correlated with BMD from cQCT and with TBS but not with BMD from DXA.

CONCLUSIONS

TBS and cQCT could accurately assess bone status in CKD patients since DXA may overestimate LS BMD, likely due to an increased AAC volume.

Electronic cigarette use and suicidal behaviors among adolescents

BACKGROUND

The aim of this study was to investigate the association of the use of e-cigarettes, which are regarded as alternatives to traditional cigarettes and are easy to purchase, with suicidal behaviors in adolescents.

METHODS

Data for 5405 middle and high school students aged 13-18 years who had used e-cigarettes were extracted from the 2016 Korean Youth Risk Behavior Web-Based Survey and analyzed using multiple logistic regression.

RESULTS

Among those who used e-cigarettes for 1-30 days in the past 30 days, suicidal ideation (odds ratio [OR] = 1.58, 95% CI, 1.31-1.89), suicidal plans (OR = 2.44, 95% CI, 1.94-3.08), suicidal attempts (OR = 2.44, 95% CI, 1.85-3.22) and serious attempts (OR = 3.09, 95% CI, 1.51-6.32) were higher compared to those who did not use an e-cigarette in the past 30 days.

CONCLUSIONS

Suicidal behaviors are significantly higher among current adolescent e-cigarette smokers than adolescents who have not used an e-cigarette in the past 30 days. Therefore, comprehensive intervention is needed to protect adolescents' mental health.

Markers of vitamin D metabolism and premenstrual symptoms in healthy women with regular cycles

STUDY QUESTION

Are insufficient 25-hydroxyvitamin D (25(OH)D) concentrations, and other markers of vitamin D metabolism, associated with premenstrual symptoms in healthy women with regular menstrual cycles?

SUMMARY ANSWER

25(OH)D insufficiency was associated with specific physical premenstrual symptoms, while no associations were observed with psychological symptoms or with other markers of vitamin D metabolism.

WHAT IS KNOWN ALREADY

Prior studies evaluating vitamin D and premenstrual symptoms have yielded mixed results, and it is unknown whether 25(OH)D insufficiency and other markers of vitamin D metabolism are associated with premenstrual symptoms.

STUDY DESIGN, SIZE, DURATION

We used two cohorts of women with regular menstrual cycles; 1191 women aged 18-40 years in EAGeR (cross-sectional analysis of a prospective cohort within a randomized trial) and 76 women aged 18-44 years in BioCycle (prospective cohort). In EAGeR, premenstrual symptoms over the previous year were assessed at baseline, whereas in BioCycle, symptoms were assessed prospectively at multiple points over two menstrual cycles with symptoms queried over the previous week. In both cohorts, symptomatology was assessed via questionnaire regarding presence and severity of 14 physical and psychological symptoms the week before and after menses. Both studies measured 25(OH)D in serum. We also evaluated the association of additional markers of vitamin D metabolism and calcium homeostasis, including intact parathyroid hormone (iPTH), calcium (Ca), fibroblast growth factor 23 (FGF23), and 1,25 dihydroxyvitamin D (1,25(OH)2D) with premenstrual symptoms in the BioCycle cohort.

PARTICIPANTS/MATERIALS, SETTING, METHODS

One cohort of women actively seeking pregnancy (Effects of Aspirin in Gestation and Reproduction (EAGeR)) and one cohort not seeking pregnancy (BioCycle) were evaluated. Log-binomial regression was used to estimate risk ratios (RR) and 95% CIs for associations between insufficient 25(OH)D (<30 ng/ml) and individual premenstrual symptoms, adjusting for age, BMI, race, smoking, income, physical activity, and season of blood draw.

MAIN RESULTS AND THE ROLE OF CHANCE

25(OH)D insufficiency was associated with increased risk of breast fullness/tenderness (EAGeR RR 1.27, 95% CI 1.03, 1.55; BioCycle RR 1.37, 95% CI 0.56, 3.32) and generalized aches and pains (EAGeR RR 1.33, 95% CI 1.01, 1.78; BioCycle 1.36, 95% CI 0.41, 4.45), though results were imprecise in the BioCycle study. No associations were observed between insufficient 25(OH)D and psychological symptoms in either cohort. In BioCycle, iPTH, Ca, FGF23, and 1,25(OH) 2D were not associated with any premenstrual symptoms.

LIMITATIONS, REASONS FOR CAUTION

Results from the EAGeR study were limited by the study design, which assessed both 25(OH)D at baseline and individual premenstrual symptoms over the past year at the baseline. As such, reverse causality is a potential concern. Though premenstrual symptoms were assessed prospectively in the BioCycle cohort, the power was limited due to small sample size. However, results were fairly consistent across both studies.

WIDER IMPLICATIONS OF THE FINDINGS

Serum 25(OH)D may be associated with risk and severity of specific physical premenstrual symptoms.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the Intramural Research Program of the Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland (Contract nos. HHSN267200603423, HHSN267200603424, and HHSN267200603426). JG.R. and D.L.K. have been funded by the NIH Medical Research Scholars Program, a public-private partnership jointly supported by the NIH and generous contributions to the Foundation for the NIH by the Doris Duke Charitable Foundation (Grant #2014194), the American Association for Dental Research, the Colgate Palmolive Company, Genentech, and other private donors. For a complete list, visit the foundation website at http://www.fnih.org. The authors have no conflicts of interest to disclose.

TRIAL REGISTRATION NUMBER

Clinicaltrials.gov NCT00467363.

Constraints on Cosmic Strings Using Data from the Third Advanced LIGO-Virgo Observing Run

We search for gravitational-wave signals produced by cosmic strings in the Advanced LIGO and Virgo full O3 dataset. Search results are presented for gravitational waves produced by cosmic string loop features such as cusps, kinks, and, for the first time, kink-kink collisions. A template-based search for short-duration transient signals does not yield a detection. We also use the stochastic gravitational-wave background energy density upper limits derived from the O3 data to constrain the cosmic string tension Gμ as a function of the number of kinks, or the number of cusps, for two cosmic string loop distribution models. Additionally, we develop and test a third model that interpolates between these two models. Our results improve upon the previous LIGO-Virgo constraints on Gμ by 1 to 2 orders of magnitude depending on the model that is tested. In particular, for the one-loop distribution model, we set the most competitive constraints to date: Gμ≲4×10^{-15}. In the case of cosmic strings formed at the end of inflation in the context of grand unified theories, these results challenge simple inflationary models.

Constraints on Cosmic Strings Using Data from the Third Advanced LIGO-Virgo Observing Run

We search for gravitational-wave signals produced by cosmic strings in the Advanced LIGO and Virgo full O3 dataset. Search results are presented for gravitational waves produced by cosmic string loop features such as cusps, kinks, and, for the first time, kink-kink collisions. A template-based search for short-duration transient signals does not yield a detection. We also use the stochastic gravitational-wave background energy density upper limits derived from the O3 data to constrain the cosmic string tension Gμ as a function of the number of kinks, or the number of cusps, for two cosmic string loop distribution models. Additionally, we develop and test a third model that interpolates between these two models. Our results improve upon the previous LIGO-Virgo constraints on Gμ by 1 to 2 orders of magnitude depending on the model that is tested. In particular, for the one-loop distribution model, we set the most competitive constraints to date: Gμ≲4×10^{-15}. In the case of cosmic strings formed at the end of inflation in the context of grand unified theories, these results challenge simple inflationary models.

Cannabis use while trying to conceive: a prospective cohort study evaluating associations with fecundability, live birth and pregnancy loss

STUDY QUESTION

Is cannabis use assessed via urinary metabolites and self-report during preconception associated with fecundability, live birth and pregnancy loss?

SUMMARY ANSWER

Preconception cannabis use was associated with reduced fecundability among women with a history of pregnancy loss attempting pregnancy despite an increased frequency of intercourse.

WHAT IS KNOWN ALREADY

Cannabis use continues to rise despite limited evidence of safety during critical windows of pregnancy establishment. While existing studies suggest that self-reported cannabis use is not associated with fecundability, self-report may not be reliable.

STUDY DESIGN, SIZE, DURATION

A prospective cohort study was carried out including 1228 women followed for up to six cycles while attempting pregnancy (2006 to 2012), and throughout pregnancy if they conceived.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Women aged 18-40 years with a history of pregnancy loss (n = 1228) were recruited from four clinical centers. Women self-reported preconception cannabis use at baseline and urinary tetrahydrocannabinol metabolites were measured throughout preconception and early pregnancy (up to four times during the study: at baseline, after 6 months of follow-up or at the beginning of the conception cycle, and weeks 4 and 8 of pregnancy). Time to hCG-detected pregnancy, and incidence of live birth and pregnancy loss were prospectively assessed. Fecundability odds ratios (FOR) and 95% CI were estimated using discrete time Cox proportional hazards models, and risk ratios (RRs) and 95% CI using log-binomial regression adjusting for age, race, BMI, education level, baseline urine cotinine, alcohol use and antidepressant use.

MAIN RESULTS AND THE ROLE OF CHANCE

Preconception cannabis use was 5% (62/1228), based on combined urinary metabolite measurements and self-report, and 1.3% (11/789) used cannabis during the first 8 weeks of gestation based on urinary metabolites only. Women with preconception cannabis use had reduced fecundability (FOR 0.59; 95% CI 0.38, 0.92). Preconception cannabis use was also associated with increased frequency of intercourse per cycle (9.4 ± 7 versus 7.5 ± 7 days; P = 0.02) and higher LH (percentage change 64%, 95% CI 3, 161) and higher LH:FSH ratio (percentage change 39%, 95% CI 7, 81). There were also suggestive, though imprecise, associations with anovulation (RR 1.92, 95% CI 0.88, 4.18), and live birth (42% (19/45) cannabis users versus 55% (578/1043) nonusers; RR 0.80, 95% CI 0.57, 1.12). No associations were observed between preconception cannabis use and pregnancy loss (RR 0.81, 95% CI 0.46, 1.42). Similar results were observed after additional adjustment for parity, income, employment status and stress. We were unable to estimate associations between cannabis use during early pregnancy and pregnancy loss due to limited sample size.

LIMITATIONS, REASONS FOR CAUTION

Owing to the relatively few cannabis users in our study, we had limited ability to make conclusions regarding live birth and pregnancy loss, and were unable to account for male partner use. While results were similar after excluding smokers, alcohol use and any drug use in the past year, some residual confounding may persist due to these potential co-exposures.

WIDER IMPLICATIONS OF THE FINDINGS

These findings highlight potential risks on fecundability among women attempting pregnancy with a history of pregnancy loss and the need for expanded evidence regarding the reproductive health effects of cannabis use in the current climate of increasing legalization.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the Intramural Research Program of the Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland (Contract numbers: HHSN267200603423, HHSN267200603424, HHSN267200603426, HHSN275201300023I). Jeannie G. Radoc has been funded by the National Institutes of Health Medical Research Scholars Program, a public-private partnership supported jointly by the National Institutes of Health and generous contributions to the Foundation for the National Institutes of Health from the Doris Duke Charitable Foundation (DDCF Grant # 2014194), Genentech, Elsevier, and other private donors. The authors report no conflict of interest in this work and have nothing to disclose.

TRIAL REGISTRATION NUMBER

Clinicaltrials.gov NCT00467363.

Vitamin D and Reproductive Hormones Across the Menstrual Cycle

STUDY QUESTION

How do the calciotropic hormones (25-hydroxyvitamin D, 1,25-dihydroxyvitamin D and intact parathyroid hormone (iPTH)) vary across the menstrual cycle and do cyclic patterns of reproductive hormones (estradiol, progesterone, LH, FSH) differ by vitamin D status?

SUMMARY ANSWER

Calciotropic hormones vary minimally across the menstrual cycle; however, women with 25-hydroxyvitamin D below 30 ng/ml have lower mean estradiol across the menstrual cycle.

WHAT IS KNOWN ALREADY

Prior human studies suggest that vitamin D status is associated with fecundability, but the mechanism is unknown. Exogenous estrogens and prolonged changes in endogenous estradiol (pregnancy or menopause) influence concentrations of 25-hydroxyvitamin D. In vitro, treatment with 1,25-dihydroxyvitamin D increases steroidogenesis in ovarian granulosa cells. There are little data about changes in calciotropic hormones across the menstrual cycle or cyclic patterns of reproductive hormones by categories of vitamin D status.

STUDY DESIGN, SIZE, DURATION

A prospective cohort study of 89 self-identified white women aged 18-44, across two menstrual cycles. Participants were a subset of the BioCycle Study, a community-based study conducted at the University of Buffalo, 2005-2007.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Eligible participants had self-reported regular menstrual cycles between 21 and 35 days and were not using hormonal contraception or vitamins. Early morning fasting blood samples were drawn at up to eight study visits per cycle. Visits were timed to capture information in all cycle phases. Serum samples for 89 women (N = 163 menstrual cycles) were analyzed for estradiol, progesterone, LH, FSH and 25-hydroxyvitamin D (25(OH)D). Variability in calciotropic hormones within and across menstrual cycles was assessed using intraclass correlation coefficients and non-linear mixed models. Given the relative stability of the calciotropic hormones across the menstrual cycle, non-linear mixed models were used to examine differences in the cyclic patterns of estradiol, progesterone, LH and FSH by categories of each calciotropic hormone (split at the median). These models were conducted for all ovulatory cycles (N = 142 ovulatory menstrual cycles) and were adjusted for age, BMI (measured in clinic) and self-reported physical activity.

MAIN RESULTS AND THE ROLE OF CHANCE

Median 25(OH)D concentration was 29.5 ng/ml (SD 8.4), and only 6% of women had vitamin D deficiency (<20 ng/ml). The mean concentration of 25(OH)D did not differ between the luteal and follicular phase; however, both 1,25(OH)2D and iPTH showed small fluctuations across the menstrual cycle with the highest 1,25(OH)2D (and lowest iPTH) in the luteal phase. Compared with women who had mean 25(OH)D ≥30 ng/ml, women with lower 25(OH)D had 13.8% lower mean estradiol (95% confidence interval: -22.0, -4.7) and 10.8% lower free estradiol (95% CI: -0.07, -0.004). Additionally, compared to women with iPTH ≤36 pg/ml, women with higher concentrations of iPTH had 12.7% lower mean estradiol (95% CI: -18.7, -6.3) and 7.3% lower progesterone (95% CI: -13.3, -0.9). No differences in the cyclic pattern of any of the reproductive hormones were observed comparing cycles with higher and lower 1,25(OH)2D.

LIMITATIONS, REASONS FOR CAUTION

Women included in this study had self-reported 'regular' menstrual cycles and very few were found to have 25(OH)D deficiency. This limits our ability to examine cycle characteristics, anovulation and the effects of concentrations of the calciotropic hormones found in deficient individuals. Additionally, the results may not be generalizable to women with irregular cycles, other races, or populations with a higher prevalence of vitamin D deficiency.

WIDER IMPLICATIONS OF THE FINDINGS

These findings support current clinical practice that does not time testing for vitamin D deficiency to the menstrual cycle phase. We find that women with lower vitamin D status (lower 25(OH)D or higher iPTH) have lower mean concentrations of estradiol across the menstrual cycle. Although this study cannot identify a mechanism of action, further in vitro work or clinical trials may help elucidate the biologic mechanisms linking calciotropic and reproductive hormones.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the Intramural Research Programs of the Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health (contract numbers: HHSN275200403394C, HHSN275201100002I and Task 1 HHSN27500001) and the National Institute of Environmental Health Sciences. There are no competing interests.

Intradural Disc Herniation in the Setting of Congenital Lumbar Spinal Stenosis

INTRODUCTION

Accounting for an estimated 1.10-1.76% of all lumbar herniations, lumbar intradural disc herniation (IDH) occurs primarily in males during the fourth to fifth decades of life. While not validated, congenital lumbar spinal stenosis (CLSS) is implicated as one precipitating factor for IDH.

CASE REPORT

We report 28-year-old Hispanic female with CLSS, severe obesity, and degenerative disk disease, with a history of minimally invasive surgical (MIS) decompression for a large paracentral L4-5 disc herniation at 25. After three years, the patient developed sudden burning dysesthesias in the L4-5 dermatomes bilaterally and temporary leg weakness. Lumbar magnetic resonance imaging exhibited severe L4-5 spinal stenosis, and the patient underwent repeat MIS decompression, which again provided her with adequate symptom resolution. However, 20 days postoperatively she developed cauda equina syndrome with anal dysfunction, and bilateral leg and foot weakness. Upon open surgical exploration we discovered a tense L4-5 dural protrusion. After a dorsal durotomy, a large IDH with a ventral dural tear was identified. Subsequent to adequate debulking of the IDH, the ventral tear was repaired, and an expansile duraplasty was performed. Overall, the patient's bladder and bowel function, pain, hypoesthesia, and motor strength all improved. Two weeks after surgery she presented with a lumbar pseudomeningocele that was managed conservatively.

CONCLUSION

This report not only highlights an atypical presentation of IDH and is the first case of CLSS linked with IDH, lending support to the hypothesis that CLSS can lead to IDH, but also provides a comprehensive review of IDHs.

Parsing the microRNA genetics basis regulating skeletal muscle fiber types and meat quality traits in pigs

Muscle fibers are closely related to human diseases and livestock meat quality. However, the genetics basis of microRNAs (miRNAs) in regulating muscle fibers is not completely understood. In this study, we constructed the whole genome-wide miRNA expression profiles of porcine fast-twitch muscle [biceps femoris (Bf)] and slow-twitch muscle [soleus (Sol)], and identified hundreds of miRNAs, including four skeletal muscle-highly expressed miRNAs, ssc-miR-378, ssc-let-7f, ssc-miR-26a, and ssc-miR-27b-3p. Moreover, we identified 63 differentially expressed (DE) miRNAs between biceps femoris vs. soleus, which are the key candidate miRNAs regulating the skeletal muscle fiber types. In addition, we found that the expression of DE ssc-miR-499-5p was significantly correlated to the expression of Myoglobin (r = 0.6872, P < 0.0001) and Myosin heavy chain 7 (MYH7; r = 0.5408, P = 0.0020), and pH_45 min (r = 0.3806, P = 0.0380) and glucose content (r = -0.4382, P = 0.0154); while the expression of DE ssc-miR-499-3p was significantly correlated to the expression of Myoglobin (r = 0.5340, P = 0.0024) and pH_45 min (r = 0.4857, P = 0.0065). Taken together, our data established a sound foundation for further studies on the regulatory mechanisms of miRNAs in skeletal muscle fiber conversion and meat quality traits in livestock, and could provide a genetic explanation of the role of miRNAs in human muscular diseases.

Multilevel Risk Factors for Weight Change after Breast Cancer Diagnosis Among Black Women

Background: Weight gain after breast cancer diagnosis increases the risk of mortality. African American/Black breast cancer survivors are more likely to have excess body weight than their White counterparts, which may contribute to their higher mortality rate. Emerging evidence suggests that post-diagnosis weight gain may result from multilevel determinants. However, no study has investigated the multilevel characteristics among Black breast cancer survivors. Objective: To evaluate associations between individual-level factors and neighborhood social and built environment factors with weight change after breast cancer diagnosis among Black women. Methods: We evaluated associations of interest among 785 women enrolled in the Women's Circle of Health Follow-Up Study (WCHFS), a longitudinal study of Black breast cancer survivors in New Jersey. Weight change was primarily based on measurements at baseline and follow-up visits (Median: 10.3 and 23.2 mo. since diagnosis, respectively). Participants were grouped into categories of stable weight (52.4%), ≥3% weight loss (20.0%), and ≥3% weight gain (27.6%). Using multivariate- adjusted multinomial logistic regression and multilevel multinomial logistic regression, we evaluated relative risk ratios (RRRs) for associations between multilevel factors and post-diagnosis weight change category. Results: Black breast cancer survivors who were older at diagnosis, had higher household income, post-menopausal status, and higher baseline BMI were less likely to gain weight compared to women with stable weight. Former smoking, higher tumor stage, and chemotherapy were associated with increased relative risk of weight gain (e.g. RRR-chemo: 1.45, 95% CI: 1.01, 2.08). Black women residing in neighborhoods in the highest tertile for density of walkable destinations had a decreased relative risk of post- diagnosis weight gain (e.g. RRR-T3 highest density vs. T1 lowest: 0.39, 95% CI: 0.20, 0.75), while those residing in neighborhoods with higher density of fast food restaurants had increased relative risk of weight gain (RRR-T3 highest density vs. T1 lowest: 1.94, 95% CI: 1.23, 3.05). Conclusion: Both individual and neighborhood factors may influence the risk of weight gain among Black women after breast cancer diagnosis.

Improved survival analysis by learning shared genomic information from pan-cancer data

Motivation

Recent advances in deep learning have offered solutions to many biomedical tasks. However, there remains a challenge in applying deep learning to survival analysis using human cancer transcriptome data. As the number of genes, the input variables of survival model, is larger than the amount of available cancer patient samples, deep-learning models are prone to overfitting. To address the issue, we introduce a new deep-learning architecture called VAECox. VAECox uses transfer learning and fine tuning.

Results

We pre-trained a variational autoencoder on all RNA-seq data in 20 TCGA datasets and transferred the trained weights to our survival prediction model. Then we fine-tuned the transferred weights during training the survival model on each dataset. Results show that our model outperformed other previous models such as Cox Proportional Hazard with LASSO and ridge penalty and Cox-nnet on the 7 of 10 TCGA datasets in terms of C-index. The results signify that the transferred information obtained from entire cancer transcriptome data helped our survival prediction model reduce overfitting and show robust performance in unseen cancer patient samples.

Availability and implementation

Our implementation of VAECox is available at https://github.com/dmis-lab/VAECox.

Supplementary information

Supplementary data are available at Bioinformatics online.

FlavorGraph: a large-scale food-chemical graph for generating food representations and recommending food pairings

Food pairing has not yet been fully pioneered, despite our everyday experience with food and the large amount of food data available on the web. The complementary food pairings discovered thus far created by the intuition of talented chefs, not by scientific knowledge or statistical learning. We introduce FlavorGraph which is a large-scale food graph by relations extracted from million food recipes and information of 1,561 flavor molecules from food databases. We analyze the chemical and statistical relations of FlavorGraph and apply our graph embedding method to better represent foods in dense vectors. Our graph embedding method is a modification of metapath2vec with an additional chemical property learning layer and quantitatively outperforms other baseline methods in food clustering. Food pairing suggestions made based on the food representations of FlavorGraph help achieve better results than previous works, and the suggestions can also be used to predict relations between compounds and foods. Our research offers a new perspective on not only food pairing techniques but also food science in general.

All-cause mortality, cardiovascular mortality, and incidence of cardiovascular disease according to a screening program of cardiovascular risk in South Korea among young adults: a nationwide cohort study

OBJECTIVES

We aimed to determine whether there are any differences in all-cause and cause-specific mortality with cardiovascular disease (CVD) risk between health screening attenders and non-attenders among young adults.

STUDY DESIGN

We performed a retrospective cohort study using claim data from the Korean National Health Insurance Service database.

METHODS

Individuals aged 20-39 years who had received health screening at least once between 2002 and 2005 were classified as attenders, and the others were classified as non-attenders. After propensity score matching according to attendance of health screening, 2,060,409 attenders and 2,060,409 non-attenders were included. We estimated adjusted hazard ratios (HRs) and 95% confidence interval (CI) for all-cause mortality, cause-specific mortality, and hospitalization of CVD from 2006 to 2015.

RESULTS

Survival from all-cause mortality was greater among attenders than among non-attenders (log rank P < 0.001). Similarly, death from CVD (log rank P = 0.007) and CVD events (log rank P < 0.001) were less likely among attenders. The risk for all-cause mortality in attenders was significantly lower than that in non-attenders (HR = 0.83, 95% CI = 0.81 to 0.84). The risk for CVD mortality (HR = 0.80, 95% CI = 0.73 to 0.87) and hospitalization of CVD (HR = 0.92, 95% CI = 0.91 to 0.94) were lower in attenders. In stratified analyses, the risk for all-cause and cause-specific mortalities was lower among attenders regardless of insurance type.

CONCLUSIONS

Among young adults, the risk for all-cause mortality, CVD mortality, and hospitalization of CVD were lower for those who underwent health screenings. Future studies that evaluate the cost-effectiveness of health screening with additional consideration of psychosocial aspects are needed.

Artificial neural networks to compare the contribution of basic clinical factors, ESC SCORE, and multidimensional risk factors for cardiovascular event prediction performance: an observational study

Background/Introduction

Despite the recent increase in the availability of different data sources that can be used for prediction models for cardiovascular disease (CVD), it remains unclear to what extent such data could contribute to improving performance of the models in data-driven cardiovascular research.

Purpose

To compare the contribution of different data types in basic clinical factors, the European Society of Cardiology Systematic Coronary Risk Evaluation (ESC SCORE), and multidimensional risk factors for CVD prediction performance of artificial neural networks (ANN) using the relevant input features derived from a large-scale medical claims database.

Methods

We abstracted data through the National Health Insurance Sharing Service and collected information on 258,896 middle-aged individuals free of CVD at baseline (2009–2010) who were followed up for incident CVD until 2013. Multidimensional risk factors identifiable from the database were chosen from a systematic review of published articles. Input features in ANN were classified as follows: basic clinical factors (age, sex, and body mass index), ESC SCORE (age, sex, total cholesterol, systolic blood pressure, and cigarette smoking), and multidimensional risk factors (sociodemographic, lifestyle behavior, underlying medical conditions, dental health, medication use, etc). The data were partitioned into the training and test sets with 7:3 ratio and the performance of each ANN model was evaluated with area under the curve (AUC).

Results

The ANN model with multidimensional risk factors had higher prediction performance (AUC: 0.692) compared to the models with basic clinical factors (AUC: 0.671) and ESC SCORE (AUC: 0.684). Within the multidimensional risk factors, atrial fibrillation, family history, chronic kidney disease, retinal vein occlusion, dental caries, antipsychotics, and corticosteroid use were some of the strong predictors. However, adding multidimensional risk factors only showed marginal improvement (increase in 1.17% of AUC) compared with the ESC SCORE model.

Conclusions

Adding multidimensional risk factors as input features in the ANN only showed marginal improvement in the CVD prediction performance. When assessing cardiovascular risk from the large-scale healthcare data, variables included in the ESC SCORE should primarily be considered in the model.

Funding Acknowledgement

Type of funding source: Public grant(s) – National budget only. Main funding source(s): Kyuwoong Kim received a scholarship from the BK21-plus education program provided by the National Research Foundation of the Republic of Korea. This work is a part of Kyuwoong Kim's PhD dissertation.

Impact of arterial stiffness on diastolic function and outcomes of heart failure in patients presenting with acute myocardial infarction

Background

Relation of arterial stiffness and cardiovascular disease is well-known. However, there is no data on association between diastolic function and brachial-ankle pulse wave velocity (baPWV) and influence on heart failure outcomes in patients with acute myocardial infarction (AMI).

Methods

Among patients presenting with AMI, all subjects underwent baPWV and echocardiography were recruited. Diastolic function was categorized by 2016 guideline of ASE/EACVI left ventricular diastolic function. Heart failure outcomes were defined 1) hospitalization for heart failure, 2) cardiac death and 3) sudden cardiac arrest.

Results

Totally 1,016 subjects were enrolled (FU duration; 3.5±2.0 years, mean age 65±13, predominant male 71.3% and STEMI 40.8%) from 2012 to 2015 in the our University Hospital. Elderly, female, low BMI, higher PWV, HTN, DM and stroke were associated with higher baPWV. Increased arterial stiffness (defined as baPWV ≥1700) had higher prevalence of diastolic dysfunction compared with baPWV &lt;1700 (98.3% vs 86.2%). HF outcomes of 69 events were identified (HF; n=48, cardiac death; n=16, sudden cardiac arrest; n=2, death d/t HF; n=3). In the multivariate analysis, baPWV ≥1700, HTN and low BMI (23 kg/m2) were independent predictors for HF outcomes after adjustment with age, LVEF, DM and stroke (Table 1). Furthermore, KM curve showed that increased arterial stiffness was associated with higher prevalence of diastolic dysfunction and poor outcomes of heart failure (Figure 1).

Conclusions

In patients with AMI, arterial stiffness was associated with higher prevalence of diastolic dysfunction an independent predictor for heart failure hospitalization and cardiac death.

Figure 1. Diastolic dysfunction and HF in AMI

Funding Acknowledgement

Type of funding source: None

Biodegradable-polymer versus durable-polymer drug eluting stents for coronary artery disease: systematic review and a meta-analysis

Background

Biodegradable polymer drug-eluting stents (BP-DES) has been developed to overcome the potential drawbacks of the first-generation durable polymer drug-eluting stents (DP-DES). However, it is still under debate whether BP-DES is associated with superior efficacy and safety over DP-DES.

Purpose

We sought to compare the effects of BP-DES and DP-DES in patients with coronary artery disease.

Methods

We performed systematic review and a meta-analysis of randomized controlled trials comparing BP-DES and DP-DES on clinical outcomes in patients with coronary artery disease using CE-mark approved drug-eluting stents (DES) with at least 1-year follow-up. We included 32 studies involving 39,686 patients (BP-DES: 21,439 patients, and DP-DES: 18,247 patients). Primary outcome measure was target vessel failure (TVF; equivalent to the composite of cardiac death, target-vessel myocardial infarction, or ischemia-driven target vessel revascularization). We performed subgroup analysis according to the DP-DES generations (newer-generation DP-DES: 15,179patients, and first-generation DP-DES: 3,068 patients), and the effects of newer-generation DP-DES was compared with the BP-DES according to the BP-DES strut thickness (Ultra-thin strut [&lt;80μm]: 7,572 patients, Thin-strut [80–100μm]: 5,465 patients, and Thick-strut [≥80μm]: 5,876 patients).

Results

The odds for TVF was not significantly different between the BP-DES group and the DP-DES group in the entire study population (odds ratio [OR] 0.96, 95% confidence interval [CI] [0.90–1.02], P=0.20). The odds for TVF was significantly low in the BP-DES group relative to the first-generation DP-DES group, however the odds were comparable between the BP-DES group and the newer-generation DP-DES group (BP-DES versus first-generation DP-DES: OR 0.82, 95% CI [0.73–0.92], P&lt;0.001, and BP-DES versus newer-generation DP-DES: OR 1.00, 95% CI [0.93–1.08], P=0.99). We also found no significant differences between the BP-DES and newer-generation DP-DES, in all subgroups stratified by the BP-DES strut thickness (Ultra-thin strut BP-DES versus newer-generation DP-DES: OR 0.88, 95% CI [0.76–1.02], P=0.10, Thin-strut BP-DES versus newer-generation DP-DES: OR 1.01, 95% CI [0.90–1.13], P=0.89, and Thick strut BP-DES versus newer-generation DP-DES: OR 1.11, 95% CI [0.99–1.25], P=0.08).

Conclusions

In this meta-analysis of randomized controlled trials evaluating clinical outcomes, there was no significant differences between BP-DES and DP-DES. We found beneficial effects of BP-DES relative to the first-generation DP-DES, however, there was no statistical differences between BP-DES and newer-generation DP-DES, irrespective of the BP-DES strut thickness.

Pooled odds ratios for clinical outcomes

Funding Acknowledgement

Type of funding source: None

Risk factors for major bleeding during prolonged anticoagulation therapy in cancer-associated venous thromboembolisms: from the COMMAND VTE registry

Background/Introduction

Patients with cancer-associated venous thromboembolisms (VTEs) are at a high risk for recurrent VTEs and are recommended to receive prolonged anticoagulation therapy if they are at a low risk for bleeding. However, there are no established risk factors for bleeding during prolonged anticoagulation therapy.

Purpose

We aimed to identify the risk factors for major bleeding during prolonged anticoagulation therapy in cancer-associated VTE patients.

Methods

The COMMAND VTE Registry is a multicenter retrospective registry enrolling 3027 consecutive patients with acute symptomatic VTEs among 29 Japanese centers between January 2010 and August 2014. After excluding those without active cancer (N=2332), patients with major bleeding (N=15), death (N=17), and lost to follow-up (N=10) within 10 days after the diagnosis, and those without anticoagulation therapy beyond 10 days after the diagnosis (N=61), the present study population consisted of 592 cancer-associated VTE patients with anticoagulation therapy beyond 10 days after the diagnosis. The outcome measurement was International Society of Thrombosis and Hemostasis (ISTH) major bleeding during anticoagulation therapy beyond 10 days, which occurred before the first discontinuation of the anticoagulation therapy. We constructed a multivariable Cox proportional hazard model to estimate the hazard ratio (HR) and 95% confidence intervals (CIs) of the potential risk factors for major bleeding. As a sensitivity analysis, we used Fine and Gray's method to estimate the HR and 95% CI, taking into account the competing risk of all-cause death.

Results

During a median follow-up period of 199 days, major bleeding occurred in 72 patients (31 patients within 3 months; 41 beyond 3 months). The cumulative incidence of major bleeding was 5.8% at 3-months, 13.8% at 1-year, 17.5% at 2-year, and 28.1% at 5-years. The most frequent major bleeding site was gastrointestinal (47%), followed by intracranial (17%) and genitourinary (11%). Major bleeding tended to occur from the sites of the cancer, however, the sites of the cancer and sites of major bleeding were not necessarily concordant. The multivariable Cox regression model demonstrated that terminal cancer (adjusted HR, 4.17; 95% CI, 2.22–7.85, P&lt;0.001), chronic kidney disease (adjusted HR, 1.89; 95% CI 1.06–3.37, P=0.031), and gastrointestinal cancer (adjusted HR, 1.78; 95% CI, 1.04–3.04, P=0.037) were independently associated with an increased risk of major bleeding. After taking into account the competing risk of all-cause death, the multivariable Cox regression model demonstrated almost consistent results with the main analysis.

Conclusions

Major bleeding events were common during prolonged anticoagulation therapy in real-world cancer-associated VTE patients. Terminal cancer, chronic kidney disease, and gastrointestinal cancer were the independent risk factors for major bleeding.

Funding Acknowledgement

Type of funding source: Foundation. Main funding source(s): Research Institute for Production Development, Mitsubishi Tanabe Pharma Corporation