Comprehensive Molecular Characterization of Muscle-Invasive Bladder Cancer.

We report a comprehensive analysis of 412 muscle-invasive bladder cancers characterized by multiple TCGA analytical platforms. Fifty-eight genes were significantly mutated, and the overall mutational load was associated with APOBEC-signature mutagenesis. Clustering by mutation signature identified a high-mutation subset with 75% 5-year survival. mRNA expression clustering refined prior clustering analyses and identified a poor-survival "neuronal" subtype in which the majority of tumors lacked small cell or neuroendocrine histology. Clustering by mRNA, long non-coding RNA (lncRNA), and miRNA expression converged to identify subsets with differential epithelial-mesenchymal transition status, carcinoma in situ scores, histologic features, and survival. Our analyses identified 5 expression subtypes that may stratify response to different treatments.

The Immune Landscape of Cancer.

We performed an extensive immunogenomic analysis of more than 10,000 tumors comprising 33 diverse cancer types by utilizing data compiled by TCGA. Across cancer types, we identified six immune subtypes-wound healing, IFN-γ dominant, inflammatory, lymphocyte depleted, immunologically quiet, and TGF-ß dominant-characterized by differences in macrophage or lymphocyte signatures, Th1:Th2 cell ratio, extent of intratumoral heterogeneity, aneuploidy, extent of neoantigen load, overall cell proliferation, expression of immunomodulatory genes, and prognosis. Specific driver mutations correlated with lower (CTNNB1, NRAS, or IDH1) or higher (BRAF, TP53, or CASP8) leukocyte levels across all cancers. Multiple control modalities of the intracellular and extracellular networks (transcription, microRNAs, copy number, and epigenetic processes) were involved in tumor-immune cell interactions, both across and within immune subtypes. Our immunogenomics pipeline to characterize these heterogeneous tumors and the resulting data are intended to serve as a resource for future targeted studies to further advance the field.

Expression patterns of mismatch repair proteins in cervical cancer uncover independent prognostic value of MSH-2.

OBJECTIVE: Although early-detected cervical cancer is associated with good survival, the prognosis for late-stage disease is poor and treatment options are sparse. Mismatch repair deficiency (MMR-D) has surfaced as a predictor of prognosis and response to immune checkpoint inhibitor(s) in several cancer types, but its value in cervical cancer remains unclear. This study aimed to define the prevalence of MMR-D in cervical cancer and assess the prognostic value of MMR protein expression. METHODS: Expression of the MMR proteins MLH-1, PMS-2, MSH-2, and MSH-6 was investigated by immunohistochemical staining in a prospectively collected cervical cancer cohort (n=508) with corresponding clinicopathological and follow-up data. Sections were scored as either loss or intact expression to define MMR-D, and by a staining index, based on staining intensity and area, evaluating the prognostic potential. RNA and whole exome sequencing data were available for 72 and 75 of the patients and were used for gene set enrichment and mutational analyses, respectively. RESULTS: Five (1%) tumors were MMR-deficient, three of which were of neuroendocrine histology. MMR status did not predict survival (HR 1.93, p=0.17). MSH-2 low (n=48) was associated with poor survival (HR 1.94, p=0.02), also when adjusting for tumor stage, tumor type, and patient age (HR 2.06, p=0.013). MSH-2 low tumors had higher tumor mutational burden (p=0.003) and higher frequency of (frameshift) mutations in the double-strand break repair gene RAD50 (p<0.01). CONCLUSION: MMR-D is rare in cervical cancer, yet low MSH-2 expression is an independent predictor of poor survival.

Partnership status and retention in care among cisgender heterosexual newly diagnosed people with HIV: a cohort study.

This study examined the role of partnership status (married, unmarried-partnered, and unpartnered) on retention in care among newly diagnosed (2013-2017), cisgender heterosexual people with HIV in Birmingham, Alabama (n = 152). This study evaluated all scheduled HIV primary care provider visits for two years following diagnosis date. A kept-visit measure was calculated such that, if an individual attended ≥1 visit in each of the four 6-month intervals, they were considered to have high visit constancy. A missed-visit measure was categorized as ≥1 no-show vs. 0 no-show for first and second year after diagnosis. Odds ratios (OR) and 95% confidence intervals (CI) were calculated using logistic regression models. Models were adjusted for confounding sociodemographic and clinical characteristics. The study population was 76% Black, 57% male, median age of 37 years. Overall, 65% had high visit consistency and 34.5% had ≥1 no-show in both years. Compared to unpartnered, married individuals had higher visit constancy [AOR (95% CI): 2.88 (1.02, 8.16)]; no differences were observed among unmarried-partnered individuals. No differences in having ≥1 no-show among partnership status groups were observed for either year. These findings suggest potential success of interventions involving a social confidant in optimizing retention in care among newly diagnosed, heterosexual PWH.

Treatment-related Correlates of Growth in Children With Sickle Cell Disease in the DISPLACE Cohort.

Reduced growth and delayed maturation have been described in children with sickle cell disease (SCD). This study investigated growth and hemolysis in children with SCD in the DISPLACE (Dissemination and Implementation of Stroke Prevention Looking at the Care Environment) cohort. The database includes 5287 children, of which, 3305 had at least 2 growth measurements over a 5-year period. Body mass index was converted to z-scores (zBMI), and 19.8%, 66.1%, 14.2% of children were classified as underweight, normal, and overweight/obese, respectively. Multivariable analysis of growth was conducted and included variables: age, sex, blood pressure, hemoglobin, reticulocyte count, treatment with chronic red cell transfusion therapy (CRCT), or hydroxyurea therapy. Baseline hemoglobin levels were associated with the lower odds of being underweight (odds ratio [OR]=0.93, 95% confidence interval [CI]: 0.86-0.99), and higher odds of being overweight/obese (OR: 1.26, 95% CI: 1.17-1.36) compared with normal zBMI. CRCT was associated with being overweight/obese at baseline (OR: 1.85, 95% CI: 1.31-2.60). Overall, results showed that children who were underweight improved regardless of therapy over the 2-year time period. However, children on CRCT are at higher risk for being overweight and should be monitored closely.

A 10-gene prognostic signature points to LIMCH1 and HLA-DQB1 as important players in aggressive cervical cancer disease.

BACKGROUND: Advanced cervical cancer carries a particularly poor prognosis, and few treatment options exist. Identification of effective molecular markers is vital to improve the individualisation of treatment. We investigated transcriptional data from cervical carcinomas related to patient survival and recurrence to identify potential molecular drivers for aggressive disease. METHODS: Primary tumour RNA-sequencing profiles from 20 patients with recurrence and 53 patients with cured disease were compared. Protein levels and prognostic impact for selected markers were identified by immunohistochemistry in a population-based patient cohort. RESULTS: Comparison of tumours relative to recurrence status revealed 121 differentially expressed genes. From this gene set, a 10-gene signature with high prognostic significance (p = 0.001) was identified and validated in an independent patient cohort (p = 0.004). Protein levels of two signature genes, HLA-DQB1 (n = 389) and LIMCH1 (LIM and calponin homology domain 1) (n = 410), were independent predictors of survival (hazard ratio 2.50, p = 0.007 for HLA-DQB1 and 3.19, p = 0.007 for LIMCH1) when adjusting for established prognostic markers. HLA-DQB1 protein expression associated with programmed death ligand 1 positivity (p < 0.001). In gene set enrichment analyses, HLA-DQB1high tumours associated with immune activation and response to interferon-γ (IFN-γ). CONCLUSIONS: This study revealed a 10-gene signature with high prognostic power in cervical cancer. HLA-DQB1 and LIMCH1 are potential biomarkers guiding cervical cancer treatment.

Landscape of genomic alterations in cervical carcinomas.

Cervical cancer is responsible for 10-15% of cancer-related deaths in women worldwide. The aetiological role of infection with high-risk human papilloma viruses (HPVs) in cervical carcinomas is well established. Previous studies have also implicated somatic mutations in PIK3CA, PTEN, TP53, STK11 and KRAS as well as several copy-number alterations in the pathogenesis of cervical carcinomas. Here we report whole-exome sequencing analysis of 115 cervical carcinoma-normal paired samples, transcriptome sequencing of 79 cases and whole-genome sequencing of 14 tumour-normal pairs. Previously unknown somatic mutations in 79 primary squamous cell carcinomas include recurrent E322K substitutions in the MAPK1 gene (8%), inactivating mutations in the HLA-B gene (9%), and mutations in EP300 (16%), FBXW7 (15%), NFE2L2 (4%), TP53 (5%) and ERBB2 (6%). We also observe somatic ELF3 (13%) and CBFB (8%) mutations in 24 adenocarcinomas. Squamous cell carcinomas have higher frequencies of somatic nucleotide substitutions occurring at cytosines preceded by thymines (Tp*C sites) than adenocarcinomas. Gene expression levels at HPV integration sites were statistically significantly higher in tumours with HPV integration compared with expression of the same genes in tumours without viral integration at the same site. These data demonstrate several recurrent genomic alterations in cervical carcinomas that suggest new strategies to combat this disease.

GATK PathSeq: a customizable computational tool for the discovery and identification of microbial sequences in libraries from eukaryotic hosts.

Summary: We present an updated version of our computational pipeline, PathSeq, for the discovery and identification of microbial sequences in genomic and transcriptomic libraries from eukaryotic hosts. This pipeline is available in the Genome Analysis Toolkit (GATK) as a suite of configurable tools that can report the microbial composition of DNA or RNA short-read sequencing samples and identify unknown sequences for downstream assembly of novel organisms. GATK PathSeq enables sample analysis in minutes at low cost. In addition, these tools are built with the GATK engine and Apache Spark framework, providing robust, rapid parallelization of read quality filtering, host subtraction and microbial alignment in workstation, cluster and cloud environments. Availability and implementation: These tools are available as a part of the GATK at https://github.com/broadinstitute/gatk. Supplementary information: Supplementary data are available at Bioinformatics online.

The tumor virus landscape of AIDS-related lymphomas.

Immunodeficiency dramatically increases susceptibility to cancer as a result of reduced immune surveillance and enhanced opportunities for virus-mediated oncogenesis. Although AIDS-related lymphomas (ARLs) are frequently associated with known oncogenic viruses, many cases contain no known transforming virus. To discover novel transforming viruses, we profiled a set of ARL samples using whole transcriptome sequencing. We determined that Epstein-Barr virus (EBV) was the only virus detected in the tumor samples of this cohort, suggesting that if unidentified pathogens exist in this disease, they are present in <10% of cases or undetectable by our methods. To evaluate the role of EBV in ARL pathogenesis, we analyzed viral gene expression and found highly heterogeneous patterns of viral transcription across samples. We also found significant heterogeneity of viral antigen expression across a large cohort, with many patient samples presenting with restricted type I viral latency, indicating that EBV latency proteins are under increased immunosurveillance in the post-combined antiretroviral therapies era. Furthermore, EBV infection of lymphoma cells in HIV-positive individuals was associated with a distinct host gene expression program. These findings provide insight into the joint host-virus regulatory network of primary ARL tumor samples and expand our understanding of virus-associated oncogenesis. Our findings may also have therapeutic implications, as treatment may be personalized to target specific viral and virus-associated host processes that are only present in a subset of patients.

Clinicopathologic and molecular markers in cervical carcinoma: a prospective cohort study.

BACKGROUND: Cervical cancer is a major health problem worldwide. Identification of effective clinicopathologic and molecular markers is vital to improve treatment stratification. OBJECTIVES: The purpose of this study was to validate a set of well-defined clinicopathologic features in a large population-based, prospectively collected cervical cancer cohort to support their use in the clinic. Further, we explored p53 and human epidermal growth factor receptor 2 as potential prognostic markers in cervical cancer. STUDY DESIGN: Tissue was collected from 401 patients with cervical cancer. Clinical data that included follow-up evaluations were collected from patient journals. Histopathologic data were evaluated and revised by an expert pathologist. The prognostic impact of selected clinicopathologic variables was analyzed in the whole cohort. Tissue microarrays were prepared from 292 carcinomas, and p53 and human epidermal growth factor receptor 2 protein levels were evaluated by immunohistochemistry. Fresh frozen samples from overlapping cervical carcinomas previously were subjected to human papilloma virus typing (n=94), whole exome (n=100) and RNA (n=79) sequencing; the results were available for our analyses. RESULTS: Among the clinicopathologic variables, vascular space invasion, histologic type, and tumor size were verified as strong independent prognostic markers. High p53 protein levels were associated significantly with markers for aggressive phenotype and survival, also in multivariate survival analysis, but did not reflect TP53 mutational status. High human epidermal growth factor receptor 2 protein levels were identified in 21% of all tumors. ERBB2 amplification was associated with poor outcome (P=.003); human epidermal growth factor receptor 2 protein level was not. CONCLUSIONS: Our findings support that the Féderation Internationale de Gynécologie et d'Obstétrique s guidelines should include vascular space invasion and tumor size 2-4 cm and that careful selection of histologic type is essential for stratification of patient risk groups. High p53 levels independently predict poor survival yet do not reflect mutational status in cervical cancer. Amplified ERBB2 significantly links to poor survival, while HercepTest does not. With optimal stratification, human epidermal growth factor receptor 2-based therapy may improve cervical cancer treatment.

Characterization of HPV and host genome interactions in primary head and neck cancers.

Previous studies have established that a subset of head and neck tumors contains human papillomavirus (HPV) sequences and that HPV-driven head and neck cancers display distinct biological and clinical features. HPV is known to drive cancer by the actions of the E6 and E7 oncoproteins, but the molecular architecture of HPV infection and its interaction with the host genome in head and neck cancers have not been comprehensively described. We profiled a cohort of 279 head and neck cancers with next generation RNA and DNA sequencing and show that 35 (12.5%) tumors displayed evidence of high-risk HPV types 16, 33, or 35. Twenty-five cases had integration of the viral genome into one or more locations in the human genome with statistical enrichment for genic regions. Integrations had a marked impact on the human genome and were associated with alterations in DNA copy number, mRNA transcript abundance and splicing, and both inter- and intrachromosomal rearrangements. Many of these events involved genes with documented roles in cancer. Cancers with integrated vs. nonintegrated HPV displayed different patterns of DNA methylation and both human and viral gene expressions. Together, these data provide insight into the mechanisms by which HPV interacts with the human genome beyond expression of viral oncoproteins and suggest that specific integration events are an integral component of viral oncogenesis.

Sequence-based discovery of Bradyrhizobium enterica in cord colitis syndrome.

BACKGROUND: Immunosuppression is associated with a variety of idiopathic clinical syndromes that may have infectious causes. It has been hypothesized that the cord colitis syndrome, a complication of umbilical-cord hematopoietic stem-cell transplantation, is infectious in origin. METHODS: We performed shotgun DNA sequencing on four archived, paraffin-embedded endoscopic colon-biopsy specimens obtained from two patients with cord colitis. Computational subtraction of human and known microbial sequences and assembly of residual sequences into a bacterial draft genome were performed. We used polymerase-chain-reaction (PCR) assays and fluorescence in situ hybridization to determine whether the corresponding bacterium was present in additional patients and controls. RESULTS: DNA sequencing of the biopsy specimens revealed more than 2.5 million sequencing reads that did not match known organisms. These sequences were computationally assembled into a 7.65-Mb draft genome showing a high degree of homology with genomes of bacteria in the bradyrhizobium genus. The corresponding newly discovered bacterium was provisionally named Bradyrhizobium enterica. PCR identified B. enterica nucleotide sequences in biopsy specimens from all three additional patients with cord colitis whose samples were tested, whereas B. enterica sequences were absent in samples obtained from healthy controls and patients with colon cancer or graft-versus-host disease. CONCLUSIONS: We assembled a novel bacterial draft genome from the direct sequencing of tissue specimens from patients with cord colitis. Association of these sequences with cord colitis suggests that B. enterica may be an opportunistic human pathogen. (Funded by the National Cancer Institute and others.)

Racial disparities in individual breast cancer outcomes by hormone-receptor subtype, area-level socio-economic status and healthcare resources.

The aim of the study is to determine the influence of area-level socio-economic status and healthcare access in addition to tumor hormone-receptor subtype on individual breast cancer stage, treatment, and mortality among Non-Hispanic (NH)-Black, NH-White, and Hispanic US adults. Analysis was based on 456,217 breast cancer patients in the SEER database from 2000 to 2010. Multilevel and multivariable-adjusted logistic and Cox proportional hazards regression analysis was conducted to account for clustering by SEER registry of diagnosis. NH-Black women had greater area-level access to healthcare resources compared with women of other races. For instance, the average numbers of oncology hospitals per million population in counties with NH-Black, NH-White, and Hispanic women were 8.1, 7.7, and 5.0 respectively; average numbers of medical doctors per million in counties with NH-Black, NH-White, and Hispanic women were 100.7, 854.0, and 866.3 respectively; and average number of Ob/Gyn in counties with NH-Black, NH-White, and Hispanic women was 155.6, 127.4, and 127.3, respectively (all p values <0.001). Regardless, NH-Black women (HR 1.39, 95 % CI 1.36-1.43) and Hispanic women (HR 1.05, 95 % CI 1.03-1.08) had significantly higher breast cancer mortality compared with NH-White women even after adjusting for hormone-receptor subtype, area-level socio-economic status, and area-level healthcare access. In addition, lower county-level socio-economic status and healthcare access measures were significantly and independently associated with stage at presentation, surgery, and radiation treatment as well as mortality after adjusting for age, race/ethnicity, and HR subtype. Although breast cancer HR subtype is a strong, important, and consistent predictor of breast cancer outcomes, we still observed significant and independent influences of area-level SES and HCA on breast cancer outcomes that deserve further study and may be critical to eliminating breast cancer outcome disparities.

Genomic aberrations in cervical adenocarcinomas in Hong Kong Chinese women.

Although the rates of cervical squamous cell carcinoma have been declining, the rates of cervical adenocarcinoma are increasing in some countries. Outcomes for advanced cervical adenocarcinoma remain poor. Precision mapping of genetic alterations in cervical adenocarcinoma may enable better selection of therapies and deliver improved outcomes when combined with new sequencing diagnostics. We present whole-exome sequencing results from 15 cervical adenocarcinomas and paired normal samples from Hong Kong Chinese women. These data revealed a heterogeneous mutation spectrum and identified several frequently altered genes including FAT1, ARID1A, ERBB2 and PIK3CA. Exome sequencing identified human papillomavirus (HPV) sequences in 13 tumors in which the HPV genome might have integrated into and hence disrupted the functions of certain exons, raising the possibility that HPV integration can alter pathways other than p53 and pRb. Together, these provisionary data suggest the potential for individualized therapies for cervical adenocarcinoma based on genomic information.

Genomic analysis identifies association of Fusobacterium with colorectal carcinoma.

The tumor microenvironment of colorectal carcinoma is a complex community of genomically altered cancer cells, nonneoplastic cells, and a diverse collection of microorganisms. Each of these components may contribute to carcinogenesis; however, the role of the microbiota is the least well understood. We have characterized the composition of the microbiota in colorectal carcinoma using whole genome sequences from nine tumor/normal pairs. Fusobacterium sequences were enriched in carcinomas, confirmed by quantitative PCR and 16S rDNA sequence analysis of 95 carcinoma/normal DNA pairs, while the Bacteroidetes and Firmicutes phyla were depleted in tumors. Fusobacteria were also visualized within colorectal tumors using FISH. These findings reveal alterations in the colorectal cancer microbiota; however, the precise role of Fusobacteria in colorectal carcinoma pathogenesis requires further investigation.

Discovery and prioritization of somatic mutations in diffuse large B-cell lymphoma (DLBCL) by whole-exome sequencing.

To gain insight into the genomic basis of diffuse large B-cell lymphoma (DLBCL), we performed massively parallel whole-exome sequencing of 55 primary tumor samples from patients with DLBCL and matched normal tissue. We identified recurrent mutations in genes that are well known to be functionally relevant in DLBCL, including MYD88, CARD11, EZH2, and CREBBP. We also identified somatic mutations in genes for which a functional role in DLBCL has not been previously suspected. These genes include MEF2B, MLL2, BTG1, GNA13, ACTB, P2RY8, PCLO, and TNFRSF14. Further, we show that BCL2 mutations commonly occur in patients with BCL2/IgH rearrangements as a result of somatic hypermutation normally occurring at the IgH locus. The BCL2 point mutations are primarily synonymous, and likely caused by activation-induced cytidine deaminase-mediated somatic hypermutation, as shown by comprehensive analysis of enrichment of mutations in WRCY target motifs. Those nonsynonymous mutations that are observed tend to be found outside of the functionally important BH domains of the protein, suggesting that strong negative selection against BCL2 loss-of-function mutations is at play. Last, by using an algorithm designed to identify likely functionally relevant but infrequent mutations, we identify KRAS, BRAF, and NOTCH1 as likely drivers of DLBCL pathogenesis in some patients. Our data provide an unbiased view of the landscape of mutations in DLBCL, and this in turn may point toward new therapeutic strategies for the disease.

Deciphering the divergent transcriptomic landscapes of cervical cancer cells grown in 3D and 2D cell culture systems.

Cervical cancer remains a significant health challenge for women worldwide, with a disproportionate impact on developing regions like sub-Saharan Africa. Taking advantage of recent advancements in developing suitable preclinical models to study cell proliferation, differentiation, and gene expression, we used RNA sequencing to compare the transcriptomic profiles of SiHa cervical cancer cells grown in 3D versus 2D culture systems. Pathway analysis of 3D cultures revealed upregulation of immune activation, angiogenesis, and tissue remodeling pathways. The high expression of cytokines, chemokines, matrix metalloproteinases, and immediate early genes, suggests that 3D cultures replicate the tumor microenvironment better than 2D monolayer cultures. HPV gene expression analysis further demonstrated higher expression levels of HPV16 E1, E2, E6, and E7 genes in 3D versus 2D cultures. Further, by using a set of linear models, we identified 79 significantly differentially expressed genes in 3D culture compared to 2D culture conditions, independent of HPV16 viral gene effects. We subsequently validated five of these genes at the protein level in both the SiHa cell line and a newly developed, patient-derived cervical cancer cell line. In addition, correlation analysis identified 26 human genes positively correlated with viral genes across 2D and 3D culture conditions. The top five 3D versus 2D differentially expressed and HPV-correlated genes were validated via qRT-PCR in our patient derived cell line. Altogether, these findings suggest that 3D cultures provide superior model systems to explore mechanisms of immune evasion, cancer progression and antiviral therapeutics.

From Gut to Hormones: Unraveling the Role of Gut Microbiota in (Phyto)Estrogen Modulation in Health and Disease.

The human gut microbiota regulates estrogen metabolism through the "estrobolome," the collection of bacterial genes that encode enzymes like ß-glucuronidases and ß-glucosidases. These enzymes deconjugate and reactivate estrogen, influencing circulating levels. The estrobolome mediates the enterohepatic circulation and bioavailability of estrogen. Alterations in gut microbiota composition and estrobolome function have been associated with estrogen-related diseases like breast cancer, enometrial cancer, and polycystic ovarian syndrome (PCOS). This is likely due to dysregulated estrogen signaling partly contributed by the microbial impacts on estrogen metabolism. Dietary phytoestrogens also undergo bacterial metabolism into active metabolites like equol, which binds estrogen receptors and exhibits higher estrogenic potency than its precursor daidzein. However, the ability to produce equol varies across populations, depending on the presence of specific gut microbes. Characterizing the estrobolome and equol-producing genes across populations can provide microbiome-based biomarkers. Further research is needed to investigate specific components of the estrobolome, phytoestrogen-microbiota interactions, and mechanisms linking dysbiosis to estrogen-related pathology. However, current evidence suggests that the gut microbiota is an integral regulator of estrogen status with clinical relevance to women's health and hormonal disorders.