Faculty Bibliography

The bone marrow microenvironment has a key role in regulating haematopoiesis, but its molecular complexity and response to stress are incompletely understood. Here we map the transcriptional landscape of mouse bone marrow vascular, perivascular and osteoblast cell populations at single-cell resolution, both at homeostasis and under conditions of stress-induced haematopoiesis. This analysis revealed previously unappreciated levels of cellular heterogeneity within the bone marrow niche and resolved cellular sources of pro-haematopoietic growth factors, chemokines and membrane-bound ligands. Our studies demonstrate a considerable transcriptional remodelling of niche elements under stress conditions, including an adipocytic skewing of perivascular cells. Among the stress-induced changes, we observed that vascular Notch delta-like ligands (encoded by Dll1 and Dll4) were downregulated. In the absence of vascular Dll4, haematopoietic stem cells prematurely induced a myeloid transcriptional program. These findings refine our understanding of the cellular architecture of the bone marrow niche, reveal a dynamic and heterogeneous molecular landscape that is highly sensitive to stress and illustrate the utility of single-cell transcriptomic data in evaluating the regulation of haematopoiesis by discrete niche populations.

Near full genome sequencing (NFGS) of HIV-1 is required to assess the genetic composition of HIV-1 strains comprehensively. Population-wide, it enables a determination of the heterogeneity of HIV-1 and the emergence of novel/recombinant strains, while for each individual it constitutes a diagnostic instrument to assist targeted therapeutic measures against viral components. There is still a lack of robust and adaptable techniques for efficient NFGS from miscellaneous HIV-1 subtypes. Using rational primer design, a broad primer set was developed for the amplification and sequencing of diverse HIV-1 group M variants from plasma. Using pure subtypes as well as diverse, unique recombinant forms (URF), variable amplicon approaches were developed for NFGS comprising all functional genes. Twenty-three different genomes composed of subtypes A (A1), B, F (F2), G, CRF01_AE, CRF02_AG, and CRF22_01A1 were successfully determined. The NFGS approach was robust irrespective of viral loads (≥306 copies/mL) and amplification method. Third-generation sequencing (TGS), single genome amplification (SGA), cloning, and bulk sequencing yielded similar outcomes concerning subtype composition and recombinant breakpoint patterns. The introduction of a simple and versatile near full genome amplification, sequencing, and cloning method enables broad application in phylogenetic studies of diverse HIV-1 subtypes and can contribute to personalized HIV therapy and diagnosis.

Humans are frequently exposed to acrylamide, a probable human carcinogen found in commonplace sources such as most heated starchy foods or tobacco smoke. Prior evidence has shown that acrylamide causes cancer in rodents, yet epidemiological studies conducted to date are limited and, thus far, have yielded inconclusive data on association of human cancers with acrylamide exposure. In this study, we experimentally identify a novel and unique mutational signature imprinted by acrylamide through the effects of its reactive metabolite glycidamide. We next show that the glycidamide mutational signature is found in a full one-third of approximately 1600 tumor genomes corresponding to 19 human tumor types from 14 organs. The highest enrichment of the glycidamide signature was observed in the cancers of the lung (88% of the interrogated tumors), liver (73%), kidney (>70%), bile duct (57%), cervix (50%), and, to a lesser extent, additional cancer types. Overall, our study reveals an unexpectedly extensive contribution of acrylamide-associated mutagenesis to human cancers.

In the version of this Article originally published the same blot was inadvertently presented as both p-Rb and Cyclin A in Fig. 2a. This blot corresponds to the p-Rb panel, as can be seen in the unprocessed version of these blots in Supplementary Fig. 9. The corrected version of the panel is shown below, together with a completely uncropped image of both blots. In addition, in the 'Viral transduction' section of the Methods, the pLKO.1 plasmids encoding short hairpin RNAs targeting human Rnd1 were incorrectly listed as clones TRCN0000018338 and TRCN0000039977. The correct clone numbers are TRCN0000047434 and TRCN0000047435.

Microglandular adenosis (MGA) is a rare breast lesion reported to be associated with invasive carcinoma in up to 20-30% of cases, and has been proposed as a non-obligate precursor to basal-like breast cancers. We identified a case of matrix-producing metaplastic carcinoma with morphologic and immunohistochemical evidence of progression from MGA to atypical MGA (AMGA), carcinoma in situ (CIS) and invasive carcinoma. We performed whole exome sequencing of each component (MGA, AMGA, CIS and cancer) to characterize the mutational landscape of these foci. There was significant copy number overlap between all foci, including a segmental amplification of the CCND1 locus (partial chromosome 11 trisomy) and MYC (8q24.12-13). Using a bioinformatics approach, we were able to identify three putative mutational clusters and recurrent, stop-gain non-synonymous mutations in both ZNF862 and TP53 that were shared across all foci. Finally, we identified a novel deleterious splice-acceptor site mutation of chr5:5186164G>T (chromosome 5p15) encoding the gene, ADAMTS16, in the invasive component.

Here, we report the draft genome of Streptococcus halitosis sp. nov. strain VT-4, a novel bacterium isolated from the dorsal part of the tongue of a patient with halitosis. The genome comprised 1,880,608 bp with a GC content of 41.0%. There were 1,782 predicted protein-coding genes, including those associated with virulence and antibiotic resistance.

BACKGROUND:In Latin America (LA), there is a high incidence rate of breast cancer (BC) in premenopausal women, and the genomic features of these BC remain unknown. Here, we aim to characterize the molecular features of BC in young LA women within the framework of the PRECAMA study, a multicenter population-based case-control study of BC in premenopausal women. METHODS:Pathological tumor tissues were collected from incident cases from four LA countries. Immunohistochemistry (IHC) was performed centrally for ER, PR, HER2, Ki67, EGFR, CK5/6, and p53 protein markers. Targeted deep sequencing was done on genomic DNA extracted from formalin-fixed, paraffin-embedded tumor tissues and their paired blood samples to screen for somatic mutations in eight genes frequently mutated in BC. A subset of samples was analyzed by exome sequencing to identify somatic mutational signatures. RESULTS:The majority of cases were positive for ER or PR (168/233; 72%), and 21% were triple-negative (TN), mainly of basal type. Most tumors were positive for Ki67 (189/233; 81%). In 126 sequenced cases, TP53 and PIK3CA were the most frequently mutated genes (32.5% and 21.4%, respectively), followed by AKT1 (9.5%). TP53 mutations were more frequent in HER2-enriched and TN IHC subtypes, whereas PIK3CA/AKT1 mutations were more frequent in ER-positive tumors, as expected. Interestingly, a higher proportion of G:C>T:A mutations was observed in TP53 in PRECAMA cases compared with TCGA and METABRIC BC series (27% vs 14%). Exome-wide mutational patterns in 10 TN cases revealed alterations in signal transduction pathways and major contributions of mutational signatures caused by altered DNA repair pathways. CONCLUSIONS:These pilot results on PRECAMA tumors give a preview of the molecular features of premenopausal BC in LA. Although the overall mutation burden was as expected from data in other populations, mutational patterns observed in TP53 and exome-wide suggested possible differences in mutagenic processes giving rise to these tumors compared with other populations. Further -omics analyses of a larger number of cases in the near future will enable the investigation of relationships between these molecular features and risk factors.

BACKGROUND:Chagas disease is still prevalent in rural areas of South America. In endemic areas of Bolivia, school children are screened for the program of Chagas disease eradication of the Ministry of Health, and positive children are treated. Here, we compared the fecal, oral and skin microbiomes of children with or without Chagas disease, and before and after benznidazol treatment of infected children. METHODS:A total of 543 Bolivian children (5-14 years old) were tested for Chagas disease, and 20 positive children were treated with Benznidazole. Fecal samples and oral and skin swabs were obtained before and after treatment, together with samples from a group of 35 uninfected controls. The 16S rRNA genes were sequenced and analyzed using QIIME to determine Alpha diversity differences and community distances, and linear discriminant analyses to determine marker taxa by infection status or treatment. RESULTS:Twenty out of 543 children screened were seropositive for Chagas disease (3.7%) and were included in the study, together with 35 control children that were seronegative for the disease. Fecal samples, oral and skin swabs were taken at the beginning of the study and after the anti-protozoa therapy with Benznidazole to the chagasic children. Infected children had higher fecal Firmicutes (Streptococcus, Roseburia, Butyrivibrio, and Blautia), and lower Bacteroides and also showed some skin -but not oral- microbiota differences. Treatment eliminated the fecal microbiota differences from control children, increasing Dialister (class Clostridia) and members of the Enterobacteriaceae, and decreasing Prevotella and Coprococcus, with minor effects on the oral and skin bacterial diversity. CONCLUSIONS:The results of this study show differences in the fecal microbiota associated with Chagas disease in children, and also evidence that treatment normalizes fecal microbiota (makes it more similar to that in controls), but is associated with oral and skin microbiota differences from control children. Since microbiota impacts in children, it is important to determine the effect of drugs on the children microbiota, since dysbiosis could lead to physiological effects which might be avoidable with microbiota restoration interventions.

Lung cancer remains the leading cause of cancer death worldwide1. With new treatment modalities, there has been a shift in focus to how we can predict who may respond to targeted treatments. Current data suggest that the human microbiome can affect lung cancer treatment through its effects on the systemic immune tone. Our group has shown that the lower airway microbiota of lung cancer patients is characterized by enrichment with oral commensals2 which triggers transcriptomic signatures (PI3K, MAPK) previously described in NSCLC 2,3. The impact of local lung dysbiosis on lung cancer progression and survival is unknown. Patients with suspicious nodules on imaging who underwent bronchoscopy were recruited. High-throughput sequencing of bacterial 16S rRNA-encoding gene amplicons was performed. Clustering was based on Dirichlet-Multinomial mixtures (DMM) modeling. RNAseq was performed on bronchial epithelial cells obtained by airway brushing. We focused our analysis on 83 NSCLC samples. Overall alpha-diversity showed that advanced stage (IIIb-VI) lower airway samples were more similar to buccal samples than local stage (I-IIIa), p<0.0001. In addition, worse 6-month and 1-year survival was associated with more similar alpha-diversity between lower airway and buccal samples (Figure 1A-D). Utilizing DMM two clusters were identified, Supraglottic-Predominant-Taxa (SPT) and Background-Predominant-Taxa (BPT). There was a significant increase in percentage of SPT in advance stage compared to local stage (p<0.008) Kaplan-Meir survival analysis shows worse survival in those with NSCLC who were clustered into the SPT group compared to BPT (p=0.0003, Figure 1E). With RNAseq, differentially expressed genes between advanced stage vs. local stage and 6-month vs. 1-year survival were not as pronounced as SPT vs. BPT (Figure 1F) suggesting that globally, transcriptomic changes between different stage and NSCLC survival were difficult to detect as compared to when airway microbiome were differentiated. In lung cancer, dysbiosis within the lower airway microenvironment, possibly by micro-aspiration, is associated with a worse 6-month and 1-year survival. This change is also associated with transcriptome changes in the local environment