TRAF3 loss-of-function reveals the noncanonical NF-κB pathway as a therapeutic target in diffuse large B cell lymphoma.

Here, we report recurrent focal deletions of the chr14q32.31-32 locus, including TRAF3, a negative regulator of NF-κB signaling, in de novo diffuse large B cell lymphoma (DLBCL) (24/324 cases). Integrative analysis revealed an association between TRAF3 copy number loss with accumulation of NIK, the central noncanonical (NC) NF-κB kinase, and increased NC NF-κB pathway activity. Accordingly, TRAF3 genetic ablation in isogenic DLBCL model systems caused upregulation of NIK and enhanced NC NF-κB downstream signaling. Knockdown or pharmacological inhibition of NIK in TRAF3-deficient cells differentially impaired their proliferation and survival, suggesting an acquired onco-addiction to NC NF-κB. TRAF3 ablation also led to exacerbated secretion of the immunosuppressive cytokine IL-10. Coculturing of TRAF3-deficient DLBCL cells with CD8+ T cells impaired the induction of Granzyme B and interferon (IFN) γ, which were restored following neutralization of IL-10. Our findings corroborate a direct relationship between TRAF3 genetic alterations and NC NF-κB activation, and highlight NIK as a potential therapeutic target in a defined subset of DLBCL.

Molecular determinants of clinical outcomes in a real-world diffuse large B-cell lymphoma population.

Molecular heterogeneity of diffuse large B-cell lymphoma (DLBCL) underlies the variable outcomes achieved with immunochemotherapy. However, outcomes of gene expression profiling (GEP)-defined molecular subgroups in a real-world DLBCL population remain unknown. Here we examined the prevalence and outcomes of molecular subgroups in an unselected population of 1149 patients with de novo DLBCL in British Columbia, Canada. Evaluable biopsies were profiled by fluorescence in situ hybridization (FISH), immunohistochemistry, and digital GEP to assign cell-of-origin and the so-called "double-hit signature" (DHITsig)-a signature originally described as being characteristic for high-grade B-cell lymphoma with MYC and BCL2 rearrangements (HGBCL-DH-BCL2). DHITsig was expressed in 21% of 431 germinal center B-cell-like (GCB)-DLBCL and all 55 Burkitt lymphomas examined. Reflecting this latter finding, DHITsig has been renamed the "dark zone signature" (DZsig). DZsigpos-DLBCL, non-DZsigpos GCB-DLBCL and activated B-cell-like (ABC)-DLBCL were associated with a 2 year overall survival of 57%, 89%, and 71%, respectively. 62% of DZsigpos tumors were negative for HGBCL-DH-BCL2 by FISH, but were associated with outcomes similar to HGBCL-DH-BCL2. A small group of HGBCL-DH-BCL2 that lacked DZsig expression had different molecular features compared with DZsig-expressing HGBCL-DH-BCL2 and were associated with favorable outcomes comparable to DLBCL, not otherwise specified. DZsigpos and ABC-DLBCL had a shorter diagnosis-to-treatment interval (DTI) than GCB-DLBCL, with this metric being associated with outcome. In conclusion, DZsig expression extends beyond HGBCL-DH-BCL2 and captures a poor-prognosis DLBCL subgroup with short DTI, including patients unidentifiable by routine FISH testing, that should be considered for treatment intensification or novel therapies in prospective trials.

Single-cell profiling reveals the importance of CXCL13/CXCR5 axis biology in lymphocyte-rich classic Hodgkin lymphoma.

Lymphocyte-rich classic Hodgkin lymphoma (LR-CHL) is a rare subtype of Hodgkin lymphoma. Recent technical advances have allowed for the characterization of specific cross-talk mechanisms between malignant Hodgkin Reed-Sternberg (HRS) cells and different normal immune cells in the tumor microenvironment (TME) of CHL. However, the TME of LR-CHL has not yet been characterized at single-cell resolution. Here, using single-cell RNA sequencing (scRNA-seq), we examined the immune cell profile of 8 cell suspension samples of LR-CHL in comparison to 20 samples of the mixed cellularity (MC, 9 cases) and nodular sclerosis (NS, 11 cases) subtypes of CHL, as well as 5 reactive lymph node controls. We also performed multicolor immunofluorescence (MC-IF) on tissue microarrays from the same patients and an independent validation cohort of 31 pretreatment LR-CHL samples. ScRNA-seq analysis identified a unique CD4+ helper T cell subset in LR-CHL characterized by high expression of Chemokine C-X-C motif ligand 13 (CXCL13) and PD-1. PD-1+CXCL13+ T cells were significantly enriched in LR-CHL compared to other CHL subtypes, and spatial analyses revealed that in 46% of the LR-CHL cases these cells formed rosettes surrounding HRS cells. MC-IF analysis revealed CXCR5+ normal B cells in close proximity to CXCL13+ T cells at significantly higher levels in LR-CHL. Moreover, the abundance of PD-1+CXCL13+ T cells in the TME was significantly associated with shorter progression-free survival in LR-CHL (P = 0.032). Taken together, our findings strongly suggest the pathogenic importance of the CXCL13/CXCR5 axis and PD-1+CXCL13+ T cells as a treatment target in LR-CHL.

Mutational landscape of gray zone lymphoma.

The mutational landscape of gray zone lymphoma (GZL) has not yet been established, and differences from related entities are largely unknown. Here, we studied coding sequence mutations of 50 Epstein-Barr virus (EBV)-negative GZLs and 20 polymorphic EBV+ diffuse large B-cell lymphoma (DLBCL) not otherwise specified (poly-EBV-L) in comparison with classical Hodgkin lymphoma (cHL), primary mediastinal large B-cell lymphoma (PMBCL), and DLBCL. Exomes of 21 GZL and 7 poly-EBV-L cases, along with paired constitutional DNA, were analyzed as a discovery cohort, followed by targeted sequencing of 217 genes in an extension cohort of 29 GZL and 13 poly-EBV-L cases. GZL cases with thymic niche involvement (anterior mediastinal mass) exhibited a mutation profile closely resembling cHL and PMBCL, with SOCS1 (45%), B2M (45%), TNFAIP3 (35%), GNA13 (35%), LRRN3 (32%), and NFKBIA (29%) being the most recurrently mutated genes. In contrast, GZL cases without thymic niche involvement (n = 18) had a significantly distinct pattern that was enriched in mutations related to apoptosis defects (TP53 [39%], BCL2 [28%], BIRC6 [22%]) and depleted in GNA13, XPO1, or NF-κB signaling pathway mutations (TNFAIP3, NFKBIE, IKBKB, NFKBIA). They also exhibited more BCL2/BCL6 rearrangements compared with thymic GZL. Poly-EBV-L cases presented a distinct mutational profile, including STAT3 mutations and a significantly lower coding mutation load in comparison with EBV- GZL. Our study highlights characteristic mutational patterns in GZL associated with presentation in the thymic niche, suggesting a common cell of origin and disease evolution overlapping with related anterior mediastinal lymphomas.

The impact of MYC and BCL2 structural variants in tumors of DLBCL morphology and mechanisms of false-negative MYC IHC.

When the World Health Organization defined high-grade B-cell lymphoma with MYC and BCL2 and/or BCL6 rearrangements (HGBL-DH/TH) as a clinical category, rearrangements were the only structural variant (SV) incorporated. An "atypical double-hit" category has been proposed, encompassing tumors with concurrent MYC and BCL2 SVs other than cooccurring translocations (ie, copy number variations [CNVs]). Although the identification of a gene expression signature (DHITsig) shared among tumors harboring MYC and BCL2 rearrangements (HGBL-DH/TH-BCL2) has confirmed a common underlying biology, the biological implication of MYC and BCL2 CNVs requires further elucidation. We performed a comprehensive analysis of MYC and BCL2 SVs, as determined by fluorescent in situ hybridization (FISH), in a cohort of 802 de novo tumors with diffuse large B-cell lymphoma morphology. Although BCL2 CNVs were associated with increased expression, MYC CNVs were not. Furthermore, MYC and BCL2 CNVs, in the context of atypical double-hit, did not confer a similar gene expression profile as HGBL-DH/TH-BCL2. Finally, although MYC immunohistochemistry (IHC) has been proposed as a screening tool for FISH testing, 2 mechanisms were observed that uncoupled MYC rearrangement from IHC positivity: (1) low MYC messenger RNA expression; and (2) false-negative IHC staining mediated by a single-nucleotide polymorphism resulting in an asparagine-to-serine substitution at the 11th amino acid residue of MYC (MYC-N11S). Taken together, these results support the current exclusion of MYC and BCL2 CNVs from HGBL-DH/TH and highlight the ability of a molecular-based classification system to identify tumors with shared biology that FISH and IHC fail to fully capture.

Live Virtual Surgery and Virtual Reality in Surgery: Potential Applications in Hand Surgery Education.

Accelerated in part by the coronavirus disease 2019 pandemic, medical education has increasingly moved into the virtual sphere in recent years. Virtual surgical education encompasses several domains, including live virtual surgery and virtual and augmented reality. These technologies range in complexity from streaming audio and video of surgeries in real-time to fully immersive virtual simulations of surgery. This article reviews the current use of virtual surgical education and its possible applications in hand surgery. Applications of virtual technologies for preoperative planning and intraoperative guidance, as well as care in underresourced settings, are discussed. The authors describe their experience creating a virtual surgery subinternship with live virtual surgeries. There are many roles virtual technology can have in surgery, and this review explores potential value these technologies may have in hand surgery.

Utilizing Computer-Assisted Design and Virtual Surgical Planning for Correction of Unilateral Lambdoid Synostosis.

Unilateral lambdoid synostosis is the rarest form of single-suture craniosynostosis. Although various surgical approaches have been described, cranial vault remodeling remains the predominant approach. To aid in surgical planning, preoperative virtual surgical modeling using a patient's presenting computed tomography scan can be used to increase reconstructive precision and to reduce operative time. Presented is a 7-month-old male with unilateral lambdoid synostosis who underwent medically modeled cranial vault reconstruction.

Probabilistic cell-type assignment of single-cell RNA-seq for tumor microenvironment profiling.

Single-cell RNA sequencing has enabled the decomposition of complex tissues into functionally distinct cell types. Often, investigators wish to assign cells to cell types through unsupervised clustering followed by manual annotation or via 'mapping' to existing data. However, manual interpretation scales poorly to large datasets, mapping approaches require purified or pre-annotated data and both are prone to batch effects. To overcome these issues, we present CellAssign, a probabilistic model that leverages prior knowledge of cell-type marker genes to annotate single-cell RNA sequencing data into predefined or de novo cell types. CellAssign automates the process of assigning cells in a highly scalable manner across large datasets while controlling for batch and sample effects. We demonstrate the advantages of CellAssign through extensive simulations and analysis of tumor microenvironment composition in high-grade serous ovarian cancer and follicular lymphoma.

Integrative genomic analysis identifies key pathogenic mechanisms in primary mediastinal large B-cell lymphoma.

Primary mediastinal large B-cell lymphoma (PMBL) represents a clinically and pathologically distinct subtype of large B-cell lymphomas. Furthermore, molecular studies, including global gene expression profiling, have provided evidence that PMBL is more closely related to classical Hodgkin lymphoma (cHL). Although targeted sequencing studies have revealed a number of mutations involved in PMBL pathogenesis, a comprehensive description of disease-associated genetic alterations and perturbed pathways is still lacking. Here, we performed whole-exome sequencing of 95 PMBL tumors to inform on oncogenic driver genes and recurrent copy number alterations. The integration of somatic gene mutations with gene expression signatures provides further insights into genotype-phenotype interrelation in PMBL. We identified highly recurrent oncogenic mutations in the Janus kinase-signal transducer and activator of transcription and nuclear factor κB pathways, and provide additional evidence of the importance of immune evasion in PMBL (CIITA, CD58, B2M, CD274, and PDCD1LG2). Our analyses highlight the interferon response factor (IRF) pathway as a putative novel hallmark with frequent alterations in multiple pathway members (IRF2BP2, IRF4, and IRF8). In addition, our integrative analysis illustrates the importance of JAK1, RELB, and EP300 mutations driving oncogenic signaling. The identified driver genes were significantly more frequently mutated in PMBL compared with diffuse large B-cell lymphoma, whereas only a limited number of genes were significantly different between PMBL and cHL, emphasizing the close relation between these entities. Our study, performed on a large cohort of PMBL, highlights the importance of distinctive genetic alterations for disease taxonomy with relevance for diagnostic evaluation and therapeutic decision-making.

Somatic IL4R mutations in primary mediastinal large B-cell lymphoma lead to constitutive JAK-STAT signaling activation.

Primary mediastinal large B-cell lymphoma (PMBCL) is a distinct subtype of diffuse large B-cell lymphoma thought to arise from thymic medullary B cells. Gene mutations underlying the molecular pathogenesis of the disease are incompletely characterized. Here, we describe novel somatic IL4R mutations in 15 of 62 primary cases of PMBCL (24.2%) and in all PMBCL-derived cell lines tested. The majority of mutations (11/21; 52%) were hotspot single nucleotide variants in exon 8, leading to an I242N amino acid change in the transmembrane domain. Functional analyses establish this mutation as gain of function leading to constitutive activation of the JAK-STAT pathway and upregulation of downstream cytokine expression profiles and B cell-specific antigens. Moreover, expression of I242N mutant IL4R in a mouse xenotransplantation model conferred growth advantage in vivo. The pattern of concurrent mutations within the JAK-STAT signaling pathway suggests additive/synergistic effects of these gene mutations contributing to lymphomagenesis. Our data establish IL4R mutations as novel driver alterations and provide a strong preclinical rationale for therapeutic targeting of JAK-STAT signaling in PMBCL.

BPG: Seamless, automated and interactive visualization of scientific data.

BACKGROUND: We introduce BPG, a framework for generating publication-quality, highly-customizable plots in the R statistical environment. RESULTS: This open-source package includes multiple methods of displaying high-dimensional datasets and facilitates generation of complex multi-panel figures, making it suitable for complex datasets. A web-based interactive tool allows online figure customization, from which R code can be downloaded for integration with computational pipelines. CONCLUSION: BPG provides a new approach for linking interactive and scripted data visualization and is available at http://labs.oicr.on.ca/boutros-lab/software/bpg or via CRAN at https://cran.r-project.org/web/packages/BoutrosLab.plotting.general.

Comprehensive characterization of programmed death ligand structural rearrangements in B-cell non-Hodgkin lymphomas.

Programmed death ligands (PDLs) are immune-regulatory molecules that are frequently affected by chromosomal alterations in B-cell lymphomas. Although PDL copy-number variations are well characterized, a detailed and comprehensive analysis of structural rearrangements (SRs) and associated phenotypic consequences is largely lacking. Here, we used oligonucleotide capture sequencing of 67 formalin-fixed paraffin-embedded tissues derived from primary B-cell lymphomas and 1 cell line to detect and characterize, at base-pair resolution, SRs of the PDL locus (9p24.1; harboring PDL1/CD274 and PDL2/PDCD1LG2). We describe 36 novel PDL SRs, including 17 intrachromosomal events (inversions, duplications, deletions) and 19 translocations involving BZRAP-AS1, CD44, GET4, IL4R, KIAA0226L, MID1, RCC1, PTPN1 and segments of the immunoglobulin loci. Moreover, analysis of the precise chromosomal breakpoints reveals 2 distinct cluster breakpoint regions (CBRs) within either CD274 (CBR1) or PDCD1LG2 (CBR2). To determine the phenotypic consequences of these SRs, we performed immunohistochemistry for CD274 and PDCD1LG2 on primary pretreatment biopsies and found that PDL SRs are significantly associated with PDL protein expression. Finally, stable ectopic expression of wild-type PDCD1LG2 and the PDCD1LG2-IGHV7-81 fusion showed, in coculture, significantly reduced T-cell activation. Taken together, our data demonstrate the complementary utility of fluorescence in situ hybridization and capture sequencing approaches and provide a classification scheme for PDL SRs with implications for future studies using PDL immune-checkpoint inhibitors in B-cell lymphomas.

Phase 2 study of panobinostat with or without rituximab in relapsed diffuse large B-cell lymphoma.

The majority of diffuse large B-cell lymphoma (DLBCL) tumors contain mutations in histone-modifying enzymes (HMEs), indicating a potential therapeutic benefit of histone deacetylase inhibitors (HDIs), and preclinical data suggest that HDIs augment the effect of rituximab. In this randomized phase 2 study, we evaluated the response rate and toxicity of panobinostat, a pan-HDI administered 30 mg orally 3 times weekly, with or without rituximab, in 40 patients with relapsed or refractory de novo (n = 27) or transformed (n = 13) DLBCL. Candidate genes and whole exomes were sequenced in relapse tumor biopsies to search for molecular correlates, and these data were used to quantify circulating tumor DNA (ctDNA) in serial plasma samples. Eleven of 40 patients (28%) responded to panobinostat (95% confidence interval [CI] 14.6-43.9) and rituximab did not increase responses. The median duration of response was 14.5 months (95% CI 9.4 to "not reached"). At time of data censoring, 6 of 11 patients had not progressed. Of the genes tested for mutations, only those in MEF2B were significantly associated with response. We detected ctDNA in at least 1 plasma sample from 96% of tested patients. A significant increase in ctDNA at day 15 relative to baseline was strongly associated with lack of response (sensitivity 71.4%, specificity 100%). We conclude that panobinostat induces very durable responses in some patients with relapsed DLBCL, and early responses can be predicted by mutations in MEF2B or a significant change in ctDNA level at 15 days after treatment initiation. This clinical trial was registered at www.ClinicalTrials.gov (#NCT01238692).

Altered axonal excitability properties in facial palsy.

INTRODUCTION: Axonal excitability measures give insight into the biophysical properties of peripheral nerve axons. In this study we applied these techniques to the study of facial palsy. METHODS: Thirty patients with established facial palsy due to unresolved Bell's palsy or herpes zoster (>6 months duration), tumor invasion of the facial nerve, or traumatic facial nerve injury were assessed using facial nerve excitability techniques. RESULTS: Full recordings were obtained in 23 patients (15 unrecovered Bell's palsy or herpes zoster, 5 trauma, 3 tumor-related). Compared with normal controls, the facial palsy group demonstrated changes in stimulus response properties, threshold electrotonus, refractoriness, superexcitability, and I/V slope. Depolarizing threshold electrotonus distinguished between viral and non-viral etiologies on subgroup analysis. DISCUSSION: In this cross-sectional study, established facial palsy demonstrated findings similar to those seen in studies of regenerated axons. The improved understanding of underlying axonal characteristics offered by the technique may guide future treatment. Muscle Nerve 57: 268-272, 2018.

Children's experiences of congenital heart disease: a systematic review of qualitative studies.

We aimed to describe the experiences of children and adolescents with congenital heart disease (CHD). Electronic databases were searched until August 2016. Qualitative studies of children's perspectives on CHD were included. Data was extracted using thematic synthesis. From 44 studies from 12 countries involving 995 children, we identified 6 themes: disrupting normality (denying the diagnosis, oscillating between sickness and health, destabilizing the family dynamic), powerlessness in deteriorating health (preoccupation with impending mortality, vulnerability to catastrophic complications, exhaustion from medical testing), enduring medical ordeals (traumatized by invasive procedures, disappointed by treatment failure, displaced by transition, valuing empathy and continuity in care, overcoming uncertainty with information), warring with the body (losing stamina, distressing inability to participate in sport, distorted body image, testing the limits), hampering potential and goals (feeling disabled, unfair judgment and exclusion, difficulties with academic achievement, limiting attainment and maintenance of life milestones), and establishing one's own pace (demarcating disease from life, determination to survive, taking limitations in their stride, embracing the positives, finding personal enrichment, relying on social or spiritual support). CONCLUSION: Children with CHD feel vulnerable and burdened by debilitating physical symptoms, unpredictable complications, and discrimination. Clinicians may support patients by sharing recognition of these profound psychosocial consequences. What is Known: • CHD is associated with difficulties in learning and attention, school absenteeism, decreased endurance, poor body image, and peer socialization • What is lesser known is how young patients cope with the symptoms, prognostic uncertainty, and treatment burden What is New: • We found that children are challenged by lifestyle restrictions, fear of invasive procedures, impaired body image, discrimination, and uncertainty about the future. Feelings of disempowerment are intensified by the unpredictability of disease progression • Thus, strategies to improve outcomes include improved patient education on disease and lifestyle management and partnership with school teachers and counselors for unique psychosocial support.

Compendium of TCDD-mediated transcriptomic response datasets in mammalian model systems.

BACKGROUND: 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is the most potent congener of the dioxin class of environmental contaminants. Exposure to TCDD causes a wide range of toxic outcomes, ranging from chloracne to acute lethality. The severity of toxicity is highly dependent on the aryl hydrocarbon receptor (AHR). Binding of TCDD to the AHR leads to changes in transcription of numerous genes. Studies evaluating the transcriptional changes brought on by TCDD may provide valuable insight into the role of the AHR in human health and disease. We therefore compiled a collection of transcriptomic datasets that can be used to aid the scientific community in better understanding the transcriptional effects of ligand-activated AHR. RESULTS: Specifically, we have created a datasets package - TCDD.Transcriptomics - for the R statistical environment, consisting of 63 unique experiments comprising 377 samples, including various combinations of 3 species (human derived cell lines, mouse and rat), 4 tissue types (liver, kidney, white adipose tissue and hypothalamus) and a wide range of TCDD exposure times and doses. These datasets have been fully standardized using consistent preprocessing and annotation packages (available as of September 14, 2015). To demonstrate the utility of this R package, a subset of "AHR-core" genes were evaluated across the included datasets. Ahrr, Nqo1 and members of the Cyp family were significantly induced following exposure to TCDD across the studies as expected while Aldh3a1 was induced specifically in rat liver. Inmt was altered only in liver tissue and primarily by rat-AHR. CONCLUSIONS: Analysis of the "AHR-core" genes demonstrates a continued need for studies surrounding the impact of AHR-activity on the transcriptome; genes believed to be consistently regulated by ligand-activated AHR show surprisingly little overlap across species and tissues. Until now, a comprehensive assessment of the transcriptome across these studies was challenging due to differences in array platforms, processing methods and annotation versions. We believe that this package, which is freely available for download ( http://labs.oicr.on.ca/boutros-lab/tcdd-transcriptomics ) will prove to be a highly beneficial resource to the scientific community evaluating the effects of TCDD exposure as well as the variety of functions of the AHR.

SeqControl: process control for DNA sequencing.

As high-throughput sequencing continues to increase in speed and throughput, routine clinical and industrial application draws closer. These 'production' settings will require enhanced quality monitoring and quality control to optimize output and reduce costs. We developed SeqControl, a framework for predicting sequencing quality and coverage using a set of 15 metrics describing overall coverage, coverage distribution, basewise coverage and basewise quality. Using whole-genome sequences of 27 prostate cancers and 26 normal references, we derived multivariate models that predict sequencing quality and depth. SeqControl robustly predicted how much sequencing was required to reach a given coverage depth (area under the curve (AUC) = 0.993), accurately classified clinically relevant formalin-fixed, paraffin-embedded samples, and made predictions from as little as one-eighth of a sequencing lane (AUC = 0.967). These techniques can be immediately incorporated into existing sequencing pipelines to monitor data quality in real time. SeqControl is available at http://labs.oicr.on.ca/Boutros-lab/software/SeqControl/.

Range and Heterogeneity of Outcomes in Randomized Trials of Pediatric Chronic Kidney Disease.

OBJECTIVE: To determine the range and heterogeneity of outcomes reported in randomized controlled trials of interventions for children with chronic kidney disease (CKD). STUDY DESIGN: The Cochrane Kidney and Transplant Specialized Register was searched to March 2016. Randomized trials involving children across all stages of CKD were selected. All outcome domains and measurements were extracted from included trials. The frequency and characteristics of the outcome domains and measures were evaluated. RESULTS: From 205 trials included, 6158 different measurements of 100 different outcome domains were reported, with a median of 22 domains per trial (IQR 13-41). Overall, 52 domains (52%) were surrogate, 38 (38%) were clinical, and 10 (10%) were patient-reported. The 5 most commonly reported domains were blood pressure (76 [37%] trials), relapse/remission (70 [34%]), kidney function (66 [32%]), infection (61 [30%]), and height/pubertal development (51 [25%]). Mortality (14%), cardiovascular disease (4%), and quality of life (1%) were reported infrequently. The 2 most frequently reported outcomes, blood pressure and relapse/remission, had 56 and 81 different outcome measures, respectively. CONCLUSIONS: The outcomes reported in clinical trials involving children with CKD are extremely heterogeneous and are most often surrogate outcomes, rather than clinical and patient-centered outcomes such as cardiovascular disease and quality of life. Efforts to ensure consistent reporting of outcomes that are important to patients and clinicians will improve the value of trials to guide clinical decision-making. In our study, non-English articles were excluded.