Quality control and analytic best practices for testing genetic models of sex differences in large populations.

Phenotypic sex-based differences exist for many complex traits. In other cases, phenotypes may be similar, but underlying biology may vary. Thus, sex-aware genetic analyses are becoming increasingly important for understanding the mechanisms driving these differences. To this end, we provide a guide outlining the current best practices for testing various models of sex-dependent genetic effects in complex traits and disease conditions, noting that this is an evolving field. Insights from sex-aware analyses will not only teach us about the biology of complex traits but also aid in achieving the goals of precision medicine and health equity for all.

Racial and ethnic differences in epithelial ovarian cancer risk: an analysis from the Ovarian Cancer Association Consortium.

BACKGROUND: Limited estimates exist on risk factors for epithelial ovarian cancer (EOC) in Asian, Hispanic, and Native Hawaiian/Pacific Islander (NHPI) women. METHODS: Participants included 1734 Asian (785 cases, 949 controls), 266 NHPI (99 cases, 167 controls), 1149 Hispanic (505 cases, 644 controls), and 24,189 White (9,981 cases, 14,208 controls) women from 11 studies in the Ovarian Cancer Association Consortium. Logistic regression models estimated odds ratios (ORs) and 95% confidence intervals (CIs) for risk associations by race and ethnicity. RESULTS: Heterogeneity in EOC risk associations by race and ethnicity (p ≤ 0.02) was observed for oral contraceptive (OC) use, parity, tubal ligation and smoking. We observed inverse associations with EOC risk for OC use and parity across all groups; associations were strongest in NHPI and Asian women. The inverse association for tubal ligation with risk was most pronounced for NHPI participants (OR=0.25, 95% CI 0.13-0.48), versus Asian and White participants, respectively (OR=0.68, 95% CI 0.51-0.90; OR=0.78, 95% CI 0.73-0.85). CONCLUSIONS: Differences in EOC risk factor associations were observed across racial and ethnic groups, which could in part be due to varying prevalence of EOC histotypes. Inclusion of greater diversity in future studies is essential to inform prevention strategies.

Mendelian randomization in hepatology: A review of principles, opportunities, and challenges.

Mendelian randomization has become a popular tool to assess causal relationships using existing observational data. While randomized controlled trials are considered the gold standard for establishing causality between exposures and outcomes, it is not always feasible to conduct a trial. Mendelian randomization is a causal inference method that uses observational data to infer causal relationships by using genetic variation as a surrogate for the exposure of interest. Publications using the approach have increased dramatically in recent years, including in the field of hepatology. In this concise review, we describe the concepts, assumptions, and interpretation of Mendelian randomization as related to studies in hepatology. We focus on the strengths and weaknesses of the approach for a non-statistical audience, using an illustrative example to assess the causal relationship between body mass index and NAFLD.

The Interaction Between Brain-Derived Neurotrophic Factor Levels and Alcohol Consumption, Sleep Disturbance and Sex-Hormones in Alcohol Use Disorders.

AIMS: Brain-derived neurotrophic factor (BDNF) levels may be associated with alcohol use disorders (AUD) and alcohol consumption, correlate with sleep disturbance and be influenced by sex differences and sex hormones. These associations have not been examined in a single sample accounting for all these factors. METHODS: Data from 190 participants (29.4% female) with AUD were utilized. Sleep quality, craving intensity, depression, anxiety and alcohol consumption were assessed using the Pittsburgh Sleep Quality Index (PSQI), Penn Alcohol Craving Scale (PACS), Patient Health Questionnaire-9 (PHQ-9), Generalized Anxiety Disorder-7 (GAD-7) and Timeline Follow Back for 90 days(TLFB 90). Inventory of Drug Taking Situations (IDTS) assessed the tendency to drink in positive/negative emotional states. Serum BDNF (sBDNF) and plasma sex hormones (estrogen, progesterone, testosterone, FSH and SHBG) were measured. Pearson correlation analyses were used to examine the association between sBDNF and these measures in the entire sample and in men and women separately. Higher order interaction effects between these factors were evaluated for their association with sBDNF using a backward selection model. RESULTS: No significant correlations between sBDNF levels and sex hormones, PSQI, PHQ-9, PACS, IDTS scores and alcohol consumption were found (all P-values > 0.05). sBDNF levels were negatively correlated with GAD-7 scores in men (r = -0.1841; P = 0.03). When considering all quadratic and two-way interactions among PSQI, PHQ-9, GAD-7, mean and max drinks/day, number of drinking days, heavy drinking days, and sex no higher order moderating effects of sBDNF levels were found. CONCLUSION: Our study revealed no significant associations between sBDNF and alcohol measures, sleep, depression and sex hormones suggesting limited utility as a biomarker.

Polygenic prediction of bipolar disorder in a Latin American sample.

To date, bipolar disorder (BD) genetic studies and polygenic risk scores (PRSs) for BD are based primarily on populations of European descent (EUR) and lack representation from other ancestries including Latin American (LAT). Here, we describe a new LAT cohort from the Mayo Clinic Bipolar Biobank (MCBB), a multisite collaboration with recruitment sites in the United States (EUR; 1,443 cases and 777 controls) and Mexico and Chile (LAT; 211 cases and 161 controls) and use the sample to explore the performance of a BD-PRS in a LAT population. Using results from the largest genome-wide association study of BD in EUR individuals, PRSice2 and LDpred2 were used to compute BD-PRSs in the LAT and EUR samples from the MCBB. PRSs explained up to 1.4% (PRSice) and 4% (LDpred2) of the phenotypic variance on the liability scale in the LAT sample compared to 3.8% (PRSice2) and 3.4% (LDpred2) in the EUR samples. Future larger studies should further explore the differential performance of different PRS approaches across ancestries. International multisite studies, such as this one, have the potential to address diversity-related limitations of prior genomic studies and ultimately contribute to the reduction of health disparities.

Testing and estimation of X-chromosome SNP effects: Impact of model assumptions.

Interest in analyzing X chromosome single nucleotide polymorphisms (SNPs) is growing and several approaches have been proposed. Prior studies have compared power of different approaches, but bias and interpretation of coefficients have received less attention. We performed simulations to demonstrate the impact of X chromosome model assumptions on effect estimates. We investigated the coefficient biases of SNP and sex effects with commonly used models for X chromosome SNPs, including models with and without assumptions of X chromosome inactivation (XCI), and with and without SNP-sex interaction terms. Sex and SNP coefficient biases were observed when assumptions made about XCI and sex differences in SNP effect in the analysis model were inconsistent with the data-generating model. However, including a SNP-sex interaction term often eliminated these biases. To illustrate these findings, estimates under different genetic model assumptions are compared and interpreted in a real data example. Models to analyze X chromosome SNPs make assumptions beyond those made in autosomal variant analysis. Assumptions made about X chromosome SNP effects should be stated clearly when reporting and interpreting X chromosome associations. Fitting models with SNP × Sex interaction terms can avoid reliance on assumptions, eliminating coefficient bias even in the absence of sex differences in SNP effect.

A genome-wide association study of mammographic texture variation.

BACKGROUND: Breast parenchymal texture features, including grayscale variation (V), capture the patterns of texture variation on a mammogram and are associated with breast cancer risk, independent of mammographic density (MD). However, our knowledge on the genetic basis of these texture features is limited. METHODS: We conducted a genome-wide association study of V in 7040 European-ancestry women. V assessments were generated from digitized film mammograms. We used linear regression to test the single-nucleotide polymorphism (SNP)-phenotype associations adjusting for age, body mass index (BMI), MD phenotypes, and the top four genetic principal components. We further calculated genetic correlations and performed SNP-set tests of V with MD, breast cancer risk, and other breast cancer risk factors. RESULTS: We identified three genome-wide significant loci associated with V: rs138141444 (6q24.1) in ECT2L, rs79670367 (8q24.22) in LINC01591, and rs113174754 (12q22) near PGAM1P5. 6q24.1 and 8q24.22 have not previously been associated with MD phenotypes or breast cancer risk, while 12q22 is a known locus for both MD and breast cancer risk. Among known MD and breast cancer risk SNPs, we identified four variants that were associated with V at the Bonferroni-corrected thresholds accounting for the number of SNPs tested: rs335189 (5q23.2) in PRDM6, rs13256025 (8p21.2) in EBF2, rs11836164 (12p12.1) near SSPN, and rs17817449 (16q12.2) in FTO. We observed significant genetic correlations between V and mammographic dense area (rg = 0.79, P = 5.91 × 10-5), percent density (rg = 0.73, P = 1.00 × 10-4), and adult BMI (rg = - 0.36, P = 3.88 × 10-7). Additional significant relationships were observed for non-dense area (z = - 4.14, P = 3.42 × 10-5), estrogen receptor-positive breast cancer (z = 3.41, P = 6.41 × 10-4), and childhood body fatness (z = - 4.91, P = 9.05 × 10-7) from the SNP-set tests. CONCLUSIONS: These findings provide new insights into the genetic basis of mammographic texture variation and their associations with MD, breast cancer risk, and other breast cancer risk factors.

Towards defining morphologic parameters of normal parous and nulliparous breast tissues by artificial intelligence.

BACKGROUND: Breast terminal duct lobular units (TDLUs), the source of most breast cancer (BC) precursors, are shaped by age-related involution, a gradual process, and postpartum involution (PPI), a dramatic inflammatory process that restores baseline microanatomy after weaning. Dysregulated PPI is implicated in the pathogenesis of postpartum BCs. We propose that assessment of TDLUs in the postpartum period may have value in risk estimation, but characteristics of these tissues in relation to epidemiological factors are incompletely described. METHODS: Using validated Artificial Intelligence and morphometric methods, we analyzed digitized images of tissue sections of normal breast tissues stained with hematoxylin and eosin from donors ≤ 45 years from the Komen Tissue Bank (180 parous and 545 nulliparous). Metrics assessed by AI, included: TDLU count; adipose tissue fraction; mean acini count/TDLU; mean dilated acini; mean average acini area; mean "capillary" area; mean epithelial area; mean ratio of epithelial area versus intralobular stroma; mean mononuclear cell count (surrogate of immune cells); mean fat area proximate to TDLUs and TDLU area. We compared epidemiologic characteristics collected via questionnaire by parity status and race, using a Wilcoxon rank sum test or Fisher's exact test. Histologic features were compared between nulliparous and parous women (overall and by time between last birth and donation [recent birth: ≤ 5 years versus remote birth: > 5 years]) using multivariable regression models. RESULTS: Normal breast tissues of parous women contained significantly higher TDLU counts and acini counts, more frequent dilated acini, higher mononuclear cell counts in TDLUs and smaller acini area per TDLU than nulliparas (all multivariable analyses p < 0.001). Differences in TDLU counts and average acini size persisted for > 5 years postpartum, whereas increases in immune cells were most marked ≤ 5 years of a birth. Relationships were suggestively modified by several other factors, including demographic and reproductive characteristics, ethanol consumption and breastfeeding duration. CONCLUSIONS: Our study identified sustained expansion of TDLU numbers and reduced average acini area among parous versus nulliparous women and notable increases in immune responses within five years following childbirth. Further, we show that quantitative characteristics of normal breast samples vary with demographic features and BC risk factors.

Serum hormone levels and normal breast histology among premenopausal women.

PURPOSE: Breast terminal duct lobular units (TDLUs) are the main source of breast cancer (BC) precursors. Higher serum concentrations of hormones and growth factors have been linked to increased TDLU numbers and to elevated BC risk, with variable effects by menopausal status. We assessed associations of circulating factors with breast histology among premenopausal women using artificial intelligence (AI) and preliminarily tested whether parity modifies associations. METHODS: Pathology AI analysis was performed on 316 digital images of H&E-stained sections of normal breast tissues from Komen Tissue Bank donors ages ≤ 45 years to assess 11 quantitative metrics. Associations of circulating factors with AI metrics were assessed using regression analyses, with inclusion of interaction terms to assess effect modification. RESULTS: Higher prolactin levels were related to larger TDLU area (p < 0.001) and increased presence of adipose tissue proximate to TDLUs (p < 0.001), with less significant positive associations for acini counts (p = 0.012), dilated acini (p = 0.043), capillary area (p = 0.014), epithelial area (p = 0.007), and mononuclear cell counts (p = 0.017). Testosterone levels were associated with increased TDLU counts (p < 0.001), irrespective of parity, but associations differed by adipose tissue content. AI data for TDLU counts generally agreed with prior visual assessments. CONCLUSION: Among premenopausal women, serum hormone levels linked to BC risk were also associated with quantitative features of normal breast tissue. These relationships were suggestively modified by parity status and tissue composition. We conclude that the microanatomic features of normal breast tissue may represent a marker of BC risk.

Germline BRCA variants, lifestyle and ovarian cancer survival.

OBJECTIVE: Women with ovarian cancer who have a pathogenic germline variant in BRCA1 or BRCA2 (BRCA) have been shown to have better 5-year survival after diagnosis than women who are BRCA-wildtype (non-carriers). Modifiable lifestyle factors, including smoking, physical activity and body mass index (BMI) have previously been associated with ovarian cancer survival; however, it is unknown whether these associations differ by germline BRCA status. METHODS: We investigated measures of lifestyle prior to diagnosis in two cohorts of Australian women with invasive epithelial ovarian cancer, using Cox proportional hazards regression to calculate adjusted hazard ratios (HRs) and 95% confidence intervals (CIs). RESULTS: In the combined studies (n = 1923), there was little association between physical activity, BMI or alcohol intake and survival, and no difference by BRCA status. However, the association between current smoking status before diagnosis and poorer survival was stronger for BRCA variant carriers (HR 1.98; 95% CI 1.20-3.27) than non-carriers (HR 1.18; 95% CI 0.96-1.46; p-interaction 0.02). We saw a similar differential association with smoking when we pooled results from two additional cohorts from the USA and UK (n = 2120). Combining the results from all four studies gave a pooled-HR of 1.94 (95% CI 1.28-2.94) for current smoking among BRCA variant carriers compared to 1.08 (0.90-1.29) for non-carriers. CONCLUSIONS: Our results suggest that the adverse effect of smoking on survival may be stronger for women with a BRCA variant than those without. Thus, while smoking cessation may improve outcomes for all women with ovarian cancer, it might provide a greater benefit for BRCA variant carriers.

Into the future: A pilot study combining imaging with molecular profiling to predict resectability in ovarian cancer.

OBJECTIVE: We sought to determine the predictive value of combining tumor molecular subtype and computerized tomography (CT) imaging for surgical outcomes after primary cytoreductive surgery in advanced stage high-grade serous ovarian cancer (HGSOC) patients. METHODS: We identified 129 HGSOC patients who underwent pre-operative CT imaging and post-operative tumor mRNA profiling. A continuous CT-score indicative of overall disease burden was defined based on six imaging measurements of anatomic involvement. Molecular subtypes were derived from mRNA profiling of chemo-naïve tumors and classified as mesenchymal (MES) subtype (36%) or non-MES subtype (64%). Fischer exact tests and multivariate logistic regression examined residual disease and surgical complexity. RESULTS: Women with higher CT-scores were more likely to have MES subtype tumors (p = 0.014). MES subtypes and a high CT-score were independently predictive of macroscopic disease and high surgical complexity. In multivariate models adjusting for age, stage and American Society of Anesthesiologists (ASA) score, patients with a MES subtype and high CT-score had significantly elevated risk of macroscopic disease (OR = 26.7, 95% CI = [6.42, 187]) and were more likely to undergo high complexity surgery (OR = 9.53, 95% CI = [2.76, 40.6], compared to patients with non-MES tumor and low CT-score. CONCLUSION: Preoperative CT imaging combined with tumor molecular subtyping can identify a subset of women unlikely to have resectable disease and likely to require high complexity surgery. Along with other clinical factors, these may refine predictive scores for resection and assist treatment planning. Investigating methods for pre-surgical molecular subtyping is an important next step.

Population-based targeted sequencing of 54 candidate genes identifies PALB2 as a susceptibility gene for high-grade serous ovarian cancer.

PURPOSE: The known epithelial ovarian cancer (EOC) susceptibility genes account for less than 50% of the heritable risk of ovarian cancer suggesting that other susceptibility genes exist. The aim of this study was to evaluate the contribution to ovarian cancer susceptibility of rare deleterious germline variants in a set of candidate genes. METHODS: We sequenced the coding region of 54 candidate genes in 6385 invasive EOC cases and 6115 controls of broad European ancestry. Genes with an increased frequency of putative deleterious variants in cases versus controls were further examined in an independent set of 14 135 EOC cases and 28 655 controls from the Ovarian Cancer Association Consortium and the UK Biobank. For each gene, we estimated the EOC risks and evaluated associations between germline variant status and clinical characteristics. RESULTS: The ORs associated for high-grade serous ovarian cancer were 3.01 for PALB2 (95% CI 1.59 to 5.68; p=0.00068), 1.99 for POLK (95% CI 1.15 to 3.43; p=0.014) and 4.07 for SLX4 (95% CI 1.34 to 12.4; p=0.013). Deleterious mutations in FBXO10 were associated with a reduced risk of disease (OR 0.27, 95% CI 0.07 to 1.00, p=0.049). However, based on the Bayes false discovery probability, only the association for PALB2 in high-grade serous ovarian cancer is likely to represent a true positive. CONCLUSIONS: We have found strong evidence that carriers of PALB2 deleterious mutations are at increased risk of high-grade serous ovarian cancer. Whether the magnitude of risk is sufficiently high to warrant the inclusion of PALB2 in cancer gene panels for ovarian cancer risk testing is unclear; much larger sample sizes will be needed to provide sufficiently precise estimates for clinical counselling.

Molecular signatures of X chromosome inactivation and associations with clinical outcomes in epithelial ovarian cancer.

X chromosome inactivation (XCI) is a key epigenetic gene expression regulatory process, which may play a role in women's cancer. In particular tissues, some genes are known to escape XCI, yet patterns of XCI in ovarian cancer (OC) and their clinical associations are largely unknown. To examine XCI in OC, we integrated germline genotype with tumor copy number, gene expression and DNA methylation information from 99 OC patients. Approximately 10% of genes showed different XCI status (either escaping or being subject to XCI) compared with the studies of other tissues. Many of these genes are known oncogenes or tumor suppressors (e.g. DDX3X, TRAPPC2 and TCEANC). We also observed strong association between cis promoter DNA methylation and allele-specific expression imbalance (P = 2.0 × 10-10). Cluster analyses of the integrated data identified two molecular subgroups of OC patients representing those with regulated (N = 47) and dysregulated (N = 52) XCI. This XCI cluster membership was associated with expression of X inactive specific transcript (P = 0.002), a known driver of XCI, as well as age, grade, stage, tumor histology and extent of residual disease following surgical debulking. Patients with dysregulated XCI (N = 52) had shorter time to recurrence (HR = 2.34, P = 0.001) and overall survival time (HR = 1.87, P = 0.02) than those with regulated XCI, although results were attenuated after covariate adjustment. Similar findings were observed when restricted to high-grade serous tumors. We found evidence of a unique OC XCI profile, suggesting that XCI may play an important role in OC biology. Additional studies to examine somatic changes with paired tumor-normal tissue are needed.

Fully Automated Volumetric Breast Density Estimation from Digital Breast Tomosynthesis.

Background While digital breast tomosynthesis (DBT) is rapidly replacing digital mammography (DM) in breast cancer screening, the potential of DBT density measures for breast cancer risk assessment remains largely unexplored. Purpose To compare associations of breast density estimates from DBT and DM with breast cancer. Materials and Methods This retrospective case-control study used contralateral DM/DBT studies from women with unilateral breast cancer and age- and ethnicity-matched controls (September 19, 2011-January 6, 2015). Volumetric percent density (VPD%) was estimated from DBT using previously validated software. For comparison, the publicly available Laboratory for Individualized Breast Radiodensity Assessment software package, or LIBRA, was used to estimate area-based percent density (APD%) from raw and processed DM images. The commercial Quantra and Volpara software packages were applied to raw DM images to estimate VPD% with use of physics-based models. Density measures were compared by using Spearman correlation coefficients (r), and conditional logistic regression was performed to examine density associations (odds ratios [OR]) with breast cancer, adjusting for age and body mass index. Results A total of 132 women diagnosed with breast cancer (mean age ± standard deviation [SD], 60 years ± 11) and 528 controls (mean age, 60 years ± 11) were included. Moderate correlations between DBT and DM density measures (r = 0.32-0.75; all P < .001) were observed. Volumetric density estimates calculated from DBT (OR, 2.3 [95% CI: 1.6, 3.4] per SD for VPD%DBT) were more strongly associated with breast cancer than DM-derived density for both APD% (OR, 1.3 [95% CI: 0.9, 1.9] [P < .001] and 1.7 [95% CI: 1.2, 2.3] [P = .004] per SD for LIBRA raw and processed data, respectively) and VPD% (OR, 1.6 [95% CI: 1.1, 2.4] [P = .01] and 1.7 [95% CI: 1.2, 2.6] [P = .04] per SD for Volpara and Quantra, respectively). Conclusion The associations between quantitative breast density estimates and breast cancer risk are stronger for digital breast tomosynthesis compared with digital mammography. © RSNA, 2021 Online supplemental material is available for this article. See also the editorial by Yaffe in this issue.

Mammographic Variation Measures, Breast Density, and Breast Cancer Risk.

OBJECTIVE. Our previous work showed that variation measures, which represent breast architecture derived from mammograms, were significantly associated with breast cancer. For replication purposes, we examined the association of three variation measures (variation [V], which is measured in the image domain, and P1 and p1 [a normalized version of P1], which are derived from restricted regions in the Fourier domain) with breast cancer risk in an independent population. We also compared these measures to volumetric density measures (volumetric percent density [VPD] and dense volume [DV]) from a commercial product. MATERIALS AND METHODS. We examined 514 patients with breast cancer and 1377 control patients from a screening practice who were matched for age, date of examination, mammography unit, facility, and state of residence. Spearman rank-order correlation was used to evaluate the monotonic association between measures. Breast cancer associations were estimated using conditional logistic regression, after adjustment for age and body mass index. Odds ratios were calculated per SD increment in mammographic measure. RESULTS. These variation measures were strongly correlated with VPD (correlation, 0.68-0.80) but not with DV (correlation, 0.31-0.48). Similar to previous findings, all variation measures were significantly associated with breast cancer (odds ratio per SD: 1.30 [95% CI, 1.16-1.46] for V, 1.55 [95% CI, 1.35-1.77] for P1, and 1.51 [95% CI, 1.33-1.72] for p1). Associations of volumetric density measures with breast cancer were similar (odds ratio per SD: 1.54 [95% CI, 1.33-1.78] for VPD and 1.34 [95% CI, 1.20-1.50] for DV). When DV was included with each variation measure in the same model, all measures retained significance. CONCLUSION. Variation measures were significantly associated with breast cancer risk (comparable to the volumetric density measures) but were independent of the DV.

Understanding Verbosity: Funding Source and the Length of Consent Forms for Cancer Clinical Trials.

Consent forms are an important educational tool that helps cancer patients decide on whether or not to enroll on a clinical trial, but wordiness potentially detracts from their educational value. This single-institution study examined word counts of consent forms for all phase I, II, and III solid tumor clinical trials between 2004 and 2010. Consent forms were categorized by trial funding source: (1) pharmaceutical company; (2) National Clinical Trials Network (NCTN); (3) R01- or other non-government grants; and (4) mixed (funding from multiple sources). Three hundred fifteen consent forms were studied; these included 106 (34%) pharmaceutical company; 145 (46%) NCTN; 44 (14%) R01 type; and 20 (6%) mixed. The overall median word count was 5129 words per consent form (interquartile range (IQR) range, 4226 to 6695). The median word counts per consent form (IQR) were 5648 (4814, 6803), 5243 (4139, 6932), 4365 (3806, 5124), and 4319 (3862, 5944), respectively, based on the above funding sources, showing that pharmaceutical company trial consent forms had the highest median word count. Of note, phase of trial was associated with consent form length (phase III were wordier), and consent forms manifested a consistent increase in wordiness over time. These observations underscore a timely need to find ways to limit the verbosity of consent forms, particularly in those from pharmaceutical company trials.

Reveal, Don't Conceal: Transforming Data Visualization to Improve Transparency.

Reports highlighting the problems with the standard practice of using bar graphs to show continuous data have prompted many journals to adopt new visualization policies. These policies encourage authors to avoid bar graphs and use graphics that show the data distribution; however, they provide little guidance on how to effectively display data. We conducted a systematic review of studies published in top peripheral vascular disease journals to determine what types of figures are used, and to assess the prevalence of suboptimal data visualization practices. Among papers with data figures, 47.7% of papers used bar graphs to present continuous data. This primer provides a detailed overview of strategies for addressing this issue by (1) outlining strategies for selecting the correct type of figure depending on the study design, sample size, and the type of variable; (2) examining techniques for making effective dot plots, box plots, and violin plots; and (3) illustrating how to avoid sending mixed messages by aligning the figure structure with the study design and statistical analysis. We also present solutions to other common problems identified in the systematic review. Resources include a list of free tools and templates that authors can use to create more informative figures and an online simulator that illustrates why summary statistics are meaningful only when there are enough data to summarize. Last, we consider steps that investigators can take to improve figures in the scientific literature.

Cytotoxic T cell depletion with increasing epithelial abnormality in women with benign breast disease.

PURPOSE: We quantified cytotoxic T cells in nonmalignant breast tissues from women with and without subsequent breast cancer to assess evidence of whether immunosurveillance may be suppressed prior to tumor development. METHODS: We used an age-matched set of breast tissues from women with benign breast disease (BBD) who subsequently developed breast cancer (BBD with later BC), women with BBD who remained cancer free (BBD cancer-free), and normal Komen Tissue Bank (KTB) tissue donors (KTB controls). We evaluated terminal duct lobular units (lobules) for degree of epithelial abnormality and density of dual-positive CD8/CD103 T cells, as CD103+ cells are thought to be a subset of CD8+ cytotoxic T cells located primarily in the intraepithelial compartment. RESULTS: In 10 sets of age-matched women, 256 breast lobules were studied: 85 in BBD women with later BC, 85 in BBD cancer-free women, and 86 in KTB donors. The majority of all lobules were histologically normal (N = 143, 56%), with 65 (25%) nonproliferative fibrocystic change, and 48 (19%) proliferative epithelial change (with or without atypia). In BBD women with later BC, median CD8+/CD103+ cell density was 39.6, 31.7, and 10.5 cells/mm2 (p = 0.002) for normal, nonproliferative, and proliferative lobules. In BBD cancer-free women, median CD8+/CD103+ cell density values were 46.7, 14.3, and 0 cells/mm2 (p = 0.004) respectively. In KTB donors, CD8+/CD103+ cell density was not significantly different across the lobule types (medians 0, 5.8, 10.7, p = 0.43). CONCLUSION: In women with BBD, breast lobules with increasing epithelial abnormality show significant decreases in cytotoxic T cells as measured by CD8/CD103 staining, suggesting that impaired immunosurveillance may be a component of the earliest stages of breast cancer development.

Evaluation of LIBRA Software for Fully Automated Mammographic Density Assessment in Breast Cancer Risk Prediction.

Background The associations of density measures from the publicly available Laboratory for Individualized Breast Radiodensity Assessment (LIBRA) software with breast cancer have primarily focused on estimates from the contralateral breast at the time of diagnosis. Purpose To evaluate LIBRA measures on mammograms obtained before breast cancer diagnosis and compare their performance to established density measures. Materials and Methods For this retrospective case-control study, full-field digital mammograms in for-processing (raw) and for-presentation (processed) formats were obtained (March 2008 to December 2011) in women who developed breast cancer an average of 2 years later and in age-matched control patients. LIBRA measures included absolute dense area and area percent density (PD) from both image formats. For comparison, dense area and PD were assessed by using the research software (Cumulus), and volumetric PD (VPD) and absolute dense volume were estimated with a commercially available software (Volpara). Density measures were compared by using Spearman correlation coefficients (r), and conditional logistic regression (odds ratios [ORs] and 95% confidence intervals [CIs]) was performed to examine the associations of density measures with breast cancer by adjusting for age and body mass index. Results Evaluated were 437 women diagnosed with breast cancer (median age, 62 years ± 17 [standard deviation]) and 1225 matched control patients (median age, 61 years ± 16). LIBRA PD showed strong correlations with Cumulus PD (r = 0.77-0.84) and Volpara VPD (r = 0.85-0.90) (P < .001 for both). For LIBRA, the strongest breast cancer association was observed for PD from processed images (OR, 1.3; 95% CI: 1.1, 1.5), although the PD association from raw images was not significantly different (OR, 1.2; 95% CI: 1.1, 1.4; P = .25). Slightly stronger breast cancer associations were seen for Cumulus PD (OR, 1.5; 95% CI: 1.3, 1.8; processed images; P = .01) and Volpara VPD (OR, 1.4; 95% CI: 1.2, 1.7; raw images; P = .004) compared with LIBRA measures. Conclusion Automated density measures provided by the Laboratory for Individualized Breast Radiodensity Assessment from raw and processed mammograms correlated with established area and volumetric density measures and showed comparable breast cancer associations. © RSNA, 2020 Online supplemental material is available for this article.

fiddle: a tool to combat publication bias by getting research out of the file drawer and into the scientific community.

Statistically significant findings are more likely to be published than non-significant or null findings, leaving scientists and healthcare personnel to make decisions based on distorted scientific evidence. Continuously expanding ´file drawers' of unpublished data from well-designed experiments waste resources creates problems for researchers, the scientific community and the public. There is limited awareness of the negative impact that publication bias and selective reporting have on the scientific literature. Alternative publication formats have recently been introduced that make it easier to publish research that is difficult to publish in traditional peer reviewed journals. These include micropublications, data repositories, data journals, preprints, publishing platforms, and journals focusing on null or neutral results. While these alternative formats have the potential to reduce publication bias, many scientists are unaware that these formats exist and don't know how to use them. Our open source file drawer data liberation effort (fiddle) tool (RRID:SCR_017327 available at: http://s-quest.bihealth.org/fiddle/) is a match-making Shiny app designed to help biomedical researchers to identify the most appropriate publication format for their data. Users can search for a publication format that meets their needs, compare and contrast different publication formats, and find links to publishing platforms. This tool will assist scientists in getting otherwise inaccessible, hidden data out of the file drawer into the scientific community and literature. We briefly highlight essential details that should be included to ensure reporting quality, which will allow others to use and benefit from research published in these new formats.