Automation of high-throughput mRNA-seq library preparation: a robust, hands-free and time efficient methodology.

Over the last decade, whole transcriptome profiling, also known as RNA-sequencing (RNA-seq), has quickly gained traction as a reliable method for unbiased assessment of gene expression. Integration of RNA-seq expression data into other omics datasets (e.g., proteomics, metabolomics, or epigenetics) solidifies our understanding of cell-specific regulatory patterns, yielding pathways to investigate the key rules of gene regulation. A limitation to efficient, at-scale utilization of RNA-seq is the time-demanding library preparation workflows, which is a 2-day or longer endeavor per cohort/sample size. To tackle this bottleneck, we designed an automated workflow that increases throughput capacity, while minimizing human error to enhance reproducibility. To this end, we converted the manual protocol of the NEBNext Directional Ultra II RNA Library Prep Kit for Illumina on the Beckman Coulter liquid handler, Biomek i7 Hybrid workstation. A total of 84 RNA samples were isolated from two human cell lines and subjected to comparative manual and automated library preparation methods. Qualitative and quantitative results indicated a high degree of similarity between libraries generated manually or through automation. Yet, there was a significant reduction in both hands-on and assay time from a 2-day manual to a 9-hour automated workflow. Using linear regression analysis, we found the Pearson correlation coefficient between libraries generated manually or by automation to be almost identical to a sample being sequenced twice (R²= 0.985 vs 0.983). This demonstrates that high-throughput automated workflows can be of great benefit to genomic laboratories by enhancing efficiency of library preparation, reducing hands-on time and increasing throughput potential.

PBRM1 loss in kidney cancer unbalances the proximal tubule master transcription factor hub to repress proximal tubule differentiation.

PBRM1, a subunit of the PBAF coactivator complex that transcription factors use to activate target genes, is genetically inactivated in almost all clear cell renal cell cancers (RCCs). Using unbiased proteomic analyses, we find that PAX8, a master transcription factor driver of proximal tubule epithelial fates, recruits PBRM1/PBAF. Reverse analyses of the PAX8 interactome confirm recruitment specifically of PBRM1/PBAF and not functionally similar BAF. More conspicuous in the PAX8 hub in RCC cells, however, are corepressors, which functionally oppose coactivators. Accordingly, key PAX8 target genes are repressed in RCC versus normal kidneys, with the loss of histone lysine-27 acetylation, but intact lysine-4 trimethylation, activation marks. Re-introduction of PBRM1, or depletion of opposing corepressors using siRNA or drugs, redress coregulator imbalance and release RCC cells to terminal epithelial fates. These mechanisms thus explain RCC resemblance to the proximal tubule lineage but with suppression of the late-epithelial program that normally terminates lineage-precursor proliferation.

Differentiation therapy and the mechanisms that terminate cancer cell proliferation without harming normal cells.

Chemotherapeutic drugs have a common intent to activate apoptosis in tumor cells. However, master regulators of apoptosis (e.g., p53, p16/CDKN2A) are frequently genetically inactivated in cancers, resulting in multidrug resistance. An alternative, p53-independent method for terminating malignant proliferation is to engage terminal-differentiation. Normally, the exponential proliferation of lineage-committed progenitors, coordinated by the master transcription factor (TF) MYC, is self-limited by forward-differentiation to terminal lineage-fates. In cancers, however, this exponential proliferation is disengaged from terminal-differentiation. The mechanisms underlying this decoupling are mostly unknown. We performed a systematic review of published literature (January 2007-June 2018) to identify gene pathways linked to differentiation-failure in three treatment-recalcitrant cancers: hepatocellular carcinoma (HCC), ovarian cancer (OVC), and pancreatic ductal adenocarcinoma (PDAC). We analyzed key gene alterations in various apoptosis, proliferation and differentiation pathways to determine whether it is possible to predict treatment outcomes and suggest novel therapies. Poorly differentiated tumors were linked to poorer survival across histologies. Our analyses suggested loss-of-function events to master TF drivers of lineage-fates and their cofactors as being linked to differentiation-failure: genomic data in TCGA and ICGC databases demonstrated frequent haploinsufficiency of lineage master TFs (e.g., GATA4/6) in poorly differentiated tumors; the coactivators that these TFs use to activate genes (e.g. ARID1A, PBRM1) were also frequently inactivated by genetic mutation and/or deletion. By contrast, corepressor components (e.g., DNMT1, EED, UHRF1, and BAZ1A/B), that oppose coactivators to repress or turn off genes, were frequently amplified instead, and the level of amplification was highest in poorly differentiated lesions. This selection by neoplastic evolution towards unbalanced activity of transcriptional corepressors suggests these enzymes as candidate targets for inhibition aiming to re-engage forward-differentiation. This notion is supported by both pre-clinical and clinical trial literature.

Leukemogenic nucleophosmin mutation disrupts the transcription factor hub that regulates granulomonocytic fates.

Nucleophosmin (NPM1) is among the most frequently mutated genes in acute myeloid leukemia (AML). It is not known, however, how the resulting oncoprotein mutant NPM1 is leukemogenic. To reveal the cellular machinery in which NPM1 participates in myeloid cells, we analyzed the endogenous NPM1 protein interactome by mass spectrometry and discovered abundant amounts of the master transcription factor driver of monocyte lineage differentiation PU.1 (also known as SPI1). Mutant NPM1, which aberrantly accumulates in cytoplasm, dislocated PU.1 into cytoplasm with it. CEBPA and RUNX1, the master transcription factors that collaborate with PU.1 to activate granulomonocytic lineage fates, remained nuclear; but without PU.1, their coregulator interactions were toggled from coactivators to corepressors, repressing instead of activating more than 500 granulocyte and monocyte terminal differentiation genes. An inhibitor of nuclear export, selinexor, by locking mutant NPM1/PU.1 in the nucleus, activated terminal monocytic fates. Direct depletion of the corepressor DNA methyltransferase 1 (DNMT1) from the CEBPA/RUNX1 protein interactome using the clinical drug decitabine activated terminal granulocytic fates. Together, these noncytotoxic treatments extended survival by more than 160 days versus vehicle in a patient-derived xenotransplant model of NPM1/FLT3-mutated AML. In sum, mutant NPM1 represses monocyte and granulocyte terminal differentiation by disrupting PU.1/CEBPA/RUNX1 collaboration, a transforming action that can be reversed by pharmacodynamically directed dosing of clinical small molecules.

GATA4 loss of function in liver cancer impedes precursor to hepatocyte transition.

The most frequent chromosomal structural loss in hepatocellular carcinoma (HCC) is of the short arm of chromosome 8 (8p). Genes on the remaining homologous chromosome, however, are not recurrently mutated, and the identity of key 8p tumor-suppressor genes (TSG) is unknown. In this work, analysis of minimal commonly deleted 8p segments to identify candidate TSG implicated GATA4, a master transcription factor driver of hepatocyte epithelial lineage fate. In a murine model, liver-conditional deletion of 1 Gata4 allele to model the haploinsufficiency seen in HCC produced enlarged livers with a gene expression profile of persistent precursor proliferation and failed hepatocyte epithelial differentiation. HCC mimicked this gene expression profile, even in cases that were morphologically classified as well differentiated. HCC with intact chromosome 8p also featured GATA4 loss of function via GATA4 germline mutations that abrogated GATA4 interactions with a coactivator, MED12, or by inactivating mutations directly in GATA4 coactivators, including ARID1A. GATA4 reintroduction into GATA4-haploinsufficient HCC cells or ARID1A reintroduction into ARID1A-mutant/GATA4-intact HCC cells activated hundreds of hepatocyte genes and quenched the proliferative precursor program. Thus, disruption of GATA4-mediated transactivation in HCC suppresses hepatocyte epithelial differentiation to sustain replicative precursor phenotype.