Molecular cartography uncovers evolutionary and microenvironmental dynamics in sporadic colorectal tumors.

Colorectal cancer exhibits dynamic cellular and genetic heterogeneity during progression from precursor lesions toward malignancy. Analysis of spatial multi-omic data from 31 human colorectal specimens enabled phylogeographic mapping of tumor evolution that revealed individualized progression trajectories and accompanying microenvironmental and clonal alterations. Phylogeographic mapping ordered genetic events, classified tumors by their evolutionary dynamics, and placed clonal regions along global pseudotemporal progression trajectories encompassing the chromosomal instability (CIN+) and hypermutated (HM) pathways. Integrated single-cell and spatial transcriptomic data revealed recurring epithelial programs and infiltrating immune states along progression pseudotime. We discovered an immune exclusion signature (IEX), consisting of extracellular matrix regulators DDR1, TGFBI, PAK4, and DPEP1, that charts with CIN+ tumor progression, is associated with reduced cytotoxic cell infiltration, and shows prognostic value in independent cohorts. This spatial multi-omic atlas provides insights into colorectal tumor-microenvironment co-evolution, serving as a resource for stratification and targeted treatments.

Differential pre-malignant programs and microenvironment chart distinct paths to malignancy in human colorectal polyps.

Colorectal cancers (CRCs) arise from precursor polyps whose cellular origins, molecular heterogeneity, and immunogenic potential may reveal diagnostic and therapeutic insights when analyzed at high resolution. We present a single-cell transcriptomic and imaging atlas of the two most common human colorectal polyps, conventional adenomas and serrated polyps, and their resulting CRC counterparts. Integrative analysis of 128 datasets from 62 participants reveals adenomas arise from WNT-driven expansion of stem cells, while serrated polyps derive from differentiated cells through gastric metaplasia. Metaplasia-associated damage is coupled to a cytotoxic immune microenvironment preceding hypermutation, driven partly by antigen-presentation differences associated with tumor cell-differentiation status. Microsatellite unstable CRCs contain distinct non-metaplastic regions where tumor cells acquire stem cell properties and cytotoxic immune cells are depleted. Our multi-omic atlas provides insights into malignant progression of colorectal polyps and their microenvironment, serving as a framework for precision surveillance and prevention of CRC.

Dual indexed library design enables compatibility of in-Drop single-cell RNA-sequencing with exAMP chemistry sequencing platforms.

BACKGROUND: The increasing demand of single-cell RNA-sequencing (scRNA-seq) experiments, such as the number of experiments and cells queried per experiment, necessitates higher sequencing depth coupled to high data quality. New high-throughput sequencers, such as the Illumina NovaSeq 6000, enables this demand to be filled in a cost-effective manner. However, current scRNA-seq library designs present compatibility challenges with newer sequencing technologies, such as index-hopping, and their ability to generate high quality data has yet to be systematically evaluated. RESULTS: Here, we engineered a dual-indexed library structure, called TruDrop, on top of the inDrop scRNA-seq platform to solve these compatibility challenges, such that TruDrop libraries and standard Illumina libraries can be sequenced alongside each other on the NovaSeq. On scRNA-seq libraries, we implemented a previously-documented countermeasure to the well-described problem of index-hopping, demonstrated significant improvements in base-calling accuracy on the NovaSeq, and provided an example of multiplexing twenty-four scRNA-seq libraries simultaneously. We showed favorable comparisons in transcriptional diversity of TruDrop compared with prior inDrop libraries. CONCLUSIONS: Our approach enables cost-effective, high throughput generation of sequencing data with high quality, which should enable more routine use of scRNA-seq technologies.

Detection of drug resistant Mycobacterium tuberculosis by high-throughput sequencing of DNA isolated from acid fast bacilli smears.

BACKGROUND: Drug susceptibility testing for Mycobacterium tuberculosis (MTB) is difficult to perform in resource-limited settings where Acid Fast Bacilli (AFB) smears are commonly used for disease diagnosis and monitoring. We developed a simple method for extraction of MTB DNA from AFB smears for sequencing-based detection of mutations associated with resistance to all first and several second-line anti-tuberculosis drugs. METHODS: We isolated MTB DNA by boiling smear content in a Chelex solution, followed by column purification. We sequenced PCR-amplified segments of the rpoB, katG, embB, gyrA, gyrB, rpsL, and rrs genes, the inhA, eis, and pncA promoters and the entire pncA gene. RESULTS: We tested our assay on 1,208 clinically obtained AFB smears from Ghana (n = 379), Kenya (n = 517), Uganda (n = 262), and Zambia (n = 50). Coverage depth varied by target and slide smear grade, ranging from 300X to 12000X on average. Coverage of ≥20X was obtained for all targets in 870 (72%) slides overall. Mono-resistance (5.9%), multi-drug resistance (1.8%), and poly-resistance (2.4%) mutation profiles were detected in 10% of slides overall, and in over 32% of retreatment and follow-up cases. CONCLUSION: This rapid AFB smear DNA-based method for determining drug resistance may be useful for the diagnosis and surveillance of drug-resistant tuberculosis.

CRISPR/Cas9-targeted removal of unwanted sequences from small-RNA sequencing libraries.

In small RNA (smRNA) sequencing studies, highly abundant molecules such as adapter dimer products and tissue-specific microRNAs (miRNAs) inhibit accurate quantification of lowly expressed species. We previously developed a method to selectively deplete highly abundant miRNAs. However, this method does not deplete adapter dimer ligation products that, unless removed by gel-separation, comprise most of the library. Here, we have adapted and modified recently described methods for CRISPR/Cas9-based Depletion of Abundant Species by Hybridization ('DASH') to smRNA-seq, which we have termed miRNA and Adapter Dimer-DASH (MAD-DASH). In MAD-DASH, Cas9 is complexed with single guide RNAs (sgRNAs) targeting adapter dimer ligation products, alongside highly expressed tissue-specific smRNAs, for cleavage in vitro. This process dramatically reduces adapter dimer and targeted smRNA sequences, can be multiplexed, shows minimal off-target effects, improves the quantification of lowly expressed miRNAs from human plasma and tissue derived RNA, and obviates the need for gel-separation, greatly increasing sample throughput. Additionally, the method is fully customizable to other smRNA-seq preparation methods. Like depletion of ribosomal RNA for mRNA-seq and mitochondrial DNA for ATAC-seq, our method allows for greater proportional read-depth of non-targeted sequences.

Discovery and Validation of Circulating Biomarkers of Colorectal Adenoma by High-Depth Small RNA Sequencing.

Purpose: Colorectal cancer is the third most common cancer worldwide, causing approximately 700,000 deaths each year. The majority of colorectal cancers begin as adenomas. Definitive screening for colorectal adenomas is currently accomplished through colonoscopy but, owing largely to costs and invasiveness, is typically limited to patient groups at higher risk by virtue of age or family history. We sought to determine if blood-based small RNA markers could detect colorectal adenoma.Experimental Design: We applied high-depth small RNA sequencing to plasma from a large (n = 189) cohort of patients, balanced for age, sex, and ancestry. Our analytical methodology allowed for the detection of both microRNAs and other small RNA species. We replicated sequencing results by qPCR on plasma samples from an independent cohort (n = 140).Results: We found several small RNA species with significant associations to colorectal adenoma, including both microRNAs and non-microRNA small RNAs. These associations were robust to correction for patient covariates, including age. Among the adenoma-associated small RNAs, two, a miR-335-5p isoform and an un-annotated small RNA, were validated by qPCR in an independent cohort. A classifier trained on measures of these two RNAs in the discovery cohort yields an AUC of 0.755 (0.775 with age) for adenoma detection in the independent cohort. This classifier accurately detects adenomas in patients under 50 and is robust to sex or ancestry.Conclusions: Circulating small RNAs (including but not limited to miRNAs) discovered by sequencing and validated by qPCR identify patients with colorectal adenomas effectively. Clin Cancer Res; 24(9); 2092-9. ©2018 AACR.

Genomic regulation of invasion by STAT3 in triple negative breast cancer.

Breast cancer is a heterogeneous disease comprised of four molecular subtypes defined by whether the tumor-originating cells are luminal or basal epithelial cells. Breast cancers arising from the luminal mammary duct often express estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth receptor 2 (HER2). Tumors expressing ER and/or PR are treated with anti-hormonal therapies, while tumors overexpressing HER2 are targeted with monoclonal antibodies. Immunohistochemical detection of ER, PR, and HER2 receptors/proteins is a critical step in breast cancer diagnosis and guided treatment. Breast tumors that do not express these proteins are known as "triple negative breast cancer" (TNBC) and are typically basal-like. TNBCs are the most aggressive subtype, with the highest mortality rates and no targeted therapy, so there is a pressing need to identify important TNBC tumor regulators. The signal transducer and activator of transcription 3 (STAT3) transcription factor has been previously implicated as a constitutively active oncogene in TNBC. However, its direct regulatory gene targets and tumorigenic properties have not been well characterized. By integrating RNA-seq and ChIP-seq data from 2 TNBC tumors and 5 cell lines, we discovered novel gene signatures directly regulated by STAT3 that were enriched for processes involving inflammation, immunity, and invasion in TNBC. Functional analysis revealed that STAT3 has a key role regulating invasion and metastasis, a characteristic often associated with TNBC. Our findings suggest therapies targeting STAT3 may be important for preventing TNBC metastasis.

The mammalian LINC complex regulates genome transcriptional responses to substrate rigidity.

Mechanical integration of the nucleus with the extracellular matrix (ECM) is established by linkage between the cytoskeleton and the nucleus. This integration is hypothesized to mediate sensing of ECM rigidity, but parsing the function of nucleus-cytoskeleton linkage from other mechanisms has remained a central challenge. Here we took advantage of the fact that the LINC (linker of nucleoskeleton and cytoskeleton) complex is a known molecular linker of the nucleus to the cytoskeleton, and asked how it regulates the sensitivity of genome-wide transcription to substratum rigidity. We show that gene mechanosensitivity is preserved after LINC disruption, but reversed in direction. Combined with myosin inhibition studies, we identify genes that depend on nuclear tension for their regulation. We also show that LINC disruption does not attenuate nuclear shape sensitivity to substrate rigidity. Our results show for the first time that the LINC complex facilitates mechano-regulation of expression across the genome.

Geospatial Resolution of Human and Bacterial Diversity with City-Scale Metagenomics.

The panoply of microorganisms and other species present in our environment influence human health and disease, especially in cities, but have not been profiled with metagenomics at a city-wide scale. We sequenced DNA from surfaces across the entire New York City (NYC) subway system, the Gowanus Canal, and public parks. Nearly half of the DNA (48%) does not match any known organism; identified organisms spanned 1,688 bacterial, viral, archaeal, and eukaryotic taxa, which were enriched for harmless genera associated with skin (e.g., Acinetobacter). Predicted ancestry of human DNA left on subway surfaces can recapitulate U.S. Census demographic data, and bacterial signatures can reveal a station's history, such as marine-associated bacteria in a hurricane-flooded station. Some evidence of pathogens was found (Bacillus anthracis), but a lack of reported cases in NYC suggests that the pathogens represent a normal, urban microbiome. This baseline metagenomic map of NYC could help long-term disease surveillance, bioterrorism threat mitigation, and health management in the built environment of cities.

A 502-Base Free-Solution Electrophoretic DNA Sequencing Method Using End-Attached Wormlike Micelles.

We demonstrate that the use of wormlike nonionic micelles as drag-tags in end-labeled free-solution electrophoresis ("micelle-ELFSE") provides single-base resolution of Sanger sequencing products up to 502 bases in length, a nearly 2-fold improvement over reported ELFSE separations. "CiEj" running buffers containing 48 mM C12E5, 6 mM C10E5, and 3 M urea (32.5 °C) form wormlike micelles that provide a drag equivalent to an uncharged DNA fragment with a length (α) of 509 bases (effective Rh = 27 nm). Runtime in a 40 cm capillary (30 kV) was 35 min for elution of all products down to the 26-base primer. We also show that smaller Triton X-100 micelles give a read length of 103 bases in a 4 min run, so that a combined analysis of the Sanger products using the two buffers in separate capillaries could be completed in 14 min for the full range of lengths. A van Deemter analysis shows that resolution is limited by diffusion-based peak broadening and wall adsorption. Effects of drag-tag polydispersity are not observed, despite the inherent polydispersity of the wormlike micelles. We ascribe this to a stochastic size-sampling process that occurs as micelle size fluctuates rapidly during the runtime. A theoretical model of the process suggests that fluctuations occur with a time scale less than 10 ms, consistent with the monomer exchange process in nonionic micelles. The CiEj buffer has a low viscosity (2.7 cP) and appears to be semidilute in micelle concentration. The large drag-tag size of the CiEj buffers leads to steric segregation of the DNA and tag for short fragments and attendant mobility shifts.

Exome sequencing in amyotrophic lateral sclerosis identifies risk genes and pathways.

Amyotrophic lateral sclerosis (ALS) is a devastating neurological disease with no effective treatment. We report the results of a moderate-scale sequencing study aimed at increasing the number of genes known to contribute to predisposition for ALS. We performed whole-exome sequencing of 2869 ALS patients and 6405 controls. Several known ALS genes were found to be associated, and TBK1 (the gene encoding TANK-binding kinase 1) was identified as an ALS gene. TBK1 is known to bind to and phosphorylate a number of proteins involved in innate immunity and autophagy, including optineurin (OPTN) and p62 (SQSTM1/sequestosome), both of which have also been implicated in ALS. These observations reveal a key role of the autophagic pathway in ALS and suggest specific targets for therapeutic intervention.