Sch 36605, structure of a novel nucleoside

"Pure" hepatoid tumors of the pancreas harboring CTNNB1 somatic mutations: a new entity among solid pseudopapillary neoplasms

Hepatoid tumors (HTs) represent a rare group of neoplasms that are histologically similar to hepatocellular carcinoma but arise outside the liver. The current World Health Organization classification recognizes the hepatoid morphology of pancreatic tumors only as a possible variant of pancreatic ductal adenocarcinoma (PDAC). Here, we describe two cases of "pure" HT of the pancreas showing common features and characterized by indolent biological behavior. These tumors were roundish nodules with pushing borders, hyaline globules, and pure hepatoid histology; they were diffusely positive for β-catenin and LEF1 on immunohistochemistry. At next-generation sequencing, both neoplasms harbored only one pathogenic somatic mutation that affected the CTNNB1 gene at exon 3 and showed a loss of heterozygosity on chromosomes 18 and 21. By integrating macroscopic and microscopic features, along with their molecular profiles, we advocate that such tumors represent a distinct entity from PDAC and should be considered a new variant of solid pseudopapillary neoplasms. The recognition of this new neoplastic category may have immediate implications not only for tumor taxonomy but also for clinical practice.

Long-read mapping to repetitive reference sequences using Winnowmap2

Approximately 5-10% of the human genome remains inaccessible due to the presence of repetitive sequences such as segmental duplications and tandem repeat arrays. We show that existing long-read mappers often yield incorrect alignments and variant calls within long, near-identical repeats, as they remain vulnerable to allelic bias. In the presence of a nonreference allele within a repeat, a read sampled from that region could be mapped to an incorrect repeat copy. To address this limitation, we developed a new long-read mapping method, Winnowmap2, by using minimal confidently alignable substrings. Winnowmap2 computes each read mapping through a collection of confident subalignments. This approach is more tolerant of structural variation and more sensitive to paralog-specific variants within repeats. Our experiments highlight that Winnowmap2 successfully addresses the issue of allelic bias, enabling more accurate downstream variant calls in repetitive sequences.

Computational analysis of cancer genome sequencing data

Distilling biologically meaningful information from cancer genome sequencing data requires comprehensive identification of somatic alterations using rigorous computational methods. As the amount and complexity of sequencing data have increased, so has the number of tools for analysing them. Here, we describe the main steps involved in the bioinformatic analysis of cancer genomes, review key algorithmic developments and highlight popular tools and emerging technologies. These tools include those that identify point mutations, copy number alterations, structural variations and mutational signatures in cancer genomes. We also discuss issues in experimental design, the strengths and limitations of sequencing modalities and methodological challenges for the future.

Locus-specific expression of transposable elements in single cells with CELLO-seq

Transposable elements (TEs) regulate diverse biological processes, from early development to cancer. Expression of young TEs is difficult to measure with next-generation, single-cell sequencing technologies because their highly repetitive nature means that short complementary DNA reads cannot be unambiguously mapped to a specific locus. Single CELl LOng-read RNA-sequencing (CELLO-seq) combines long-read single cell RNA-sequencing with computational analyses to measure TE expression at unique loci. We used CELLO-seq to assess the widespread expression of TEs in two-cell mouse blastomeres as well as in human induced pluripotent stem cells. Across both species, old and young TEs showed evidence of locus-specific expression with simulations demonstrating that only a small number of very young elements in the mouse could not be mapped back to the reference with high confidence. Exploring the relationship between the expression of individual elements and putative regulators revealed large heterogeneity, with TEs within a class showing different patterns of correlation and suggesting distinct regulatory mechanisms.

Establishment of reference standards for multifaceted mosaic variant analysis

Detection of somatic mosaicism in non-proliferative cells is a new challenge in genome research, however, the accuracy of current detection strategies remains uncertain due to the lack of a ground truth. Herein, we sought to present a set of ultra-deep sequenced WES data based on reference standards generated by cell line mixtures, providing a total of 386,613 mosaic single-nucleotide variants (SNVs) and insertion-deletion mutations (INDELs) with variant allele frequencies (VAFs) ranging from 0.5% to 56%, as well as 35,113,417 non-variant and 19,936 germline variant sites as a negative control. The whole reference standard set mimics the cumulative aspect of mosaic variant acquisition such as in the early developmental stage owing to the progressive mixing of cell lines with established genotypes, ultimately unveiling 741 possible inter-sample relationships with respect to variant sharing and asymmetry in VAFs. We expect that our reference data will be essential for optimizing the current use of mosaic variant detection strategies and for developing algorithms to enable future improvements.

Measurement(s)	genotype
Technology Type(s)	DNA sequencing
Factor Type(s)	genotyping
Sample Characteristic - Organism	Homo sapiens
Sample Characteristic - Environment	cell line

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.16970041

Combining callers improves the detection of copy number variants from whole-genome sequencing

Copy Number Variants (CNVs) are deletions, duplications or insertions larger than 50 base pairs. They account for a large percentage of the normal genome variation and play major roles in human pathology. While array-based approaches have long been used to detect them in clinical practice, whole-genome sequencing (WGS) bears the promise to allow concomitant exploration of CNVs and smaller variants. However, accurately calling CNVs from WGS remains a difficult computational task, for which a consensus is still lacking. In this paper, we explore practical calling options to reach the best compromise between sensitivity and sensibility. We show that callers based on different signal (paired-end reads, split reads, coverage depth) yield complementary results. We suggest approaches combining four selected callers (Manta, Delly, ERDS, CNVnator) and a regenotyping tool (SV2), and show that this is applicable in everyday practice in terms of computation time and further interpretation. We demonstrate the superiority of these approaches over array-based Comparative Genomic Hybridization (aCGH), specifically regarding the lack of resolution in breakpoint definition and the detection of potentially relevant CNVs. Finally, we confirm our results on the NA12878 benchmark genome, as well as one clinically validated sample. In conclusion, we suggest that WGS constitutes a timely and economically valid alternative to the combination of aCGH and whole-exome sequencing.

Interleukin-7 receptor α mutational activation can initiate precursor B-cell acute lymphoblastic leukemia

Interleukin-7 receptor α (encoded by IL7R) is essential for lymphoid development. Whether acute lymphoblastic leukemia (ALL)-related IL7R gain-of-function mutations can trigger leukemogenesis remains unclear. Here, we demonstrate that lymphoid-restricted mutant IL7R, expressed at physiological levels in conditional knock-in mice, establishes a pre-leukemic stage in which B-cell precursors display self-renewal ability, initiating leukemia resembling PAX5 P80R or Ph-like human B-ALL. Full transformation associates with transcriptional upregulation of oncogenes such as Myc or Bcl2, downregulation of tumor suppressors such as Ikzf1 or Arid2, and major IL-7R signaling upregulation (involving JAK/STAT5 and PI3K/mTOR), required for leukemia cell viability. Accordingly, maximal signaling drives full penetrance and early leukemia onset in homozygous IL7R mutant animals. Notably, we identify 2 transcriptional subgroups in mouse and human Ph-like ALL, and show that dactolisib and sphingosine-kinase inhibitors are potential treatment avenues for IL-7R-related cases. Our model, a resource to explore the pathophysiology and therapeutic vulnerabilities of B-ALL, demonstrates that IL7R can initiate this malignancy.

Hominini-specific regulation of CBLN2 increases prefrontal spinogenesis

The similarities and differences between nervous systems of various species result from developmental constraints and specific adaptations1-4. Comparative analyses of the prefrontal cortex (PFC), a cerebral cortex region involved in higher-order cognition and complex social behaviours, have identified true and potential human-specific structural and molecular specializations4-8, such as an exaggerated PFC-enriched anterior-posterior dendritic spine density gradient5. These changes are probably mediated by divergence in spatiotemporal gene regulation9-17, which is particularly prominent in the midfetal human cortex15,18-20. Here we analysed human and macaque transcriptomic data15,20 and identified a transient PFC-enriched and laminar-specific upregulation of cerebellin 2 (CBLN2), a neurexin (NRXN) and glutamate receptor-δ GRID/GluD-associated synaptic organizer21-27, during midfetal development that coincided with the initiation of synaptogenesis. Moreover, we found that species differences in level of expression and laminar distribution of CBLN2 are, at least in part, due to Hominini-specific deletions containing SOX5-binding sites within a retinoic acid-responsive CBLN2 enhancer. In situ genetic humanization of the mouse Cbln2 enhancer drives increased and ectopic laminar Cbln2 expression and promotes PFC dendritic spine formation. These findings suggest a genetic and molecular basis for the anterior-posterior cortical gradient and disproportionate increase in the Hominini PFC of dendritic spines and a developmental mechanism that may link dysfunction of the NRXN-GRID-CBLN2 complex to the pathogenesis of neuropsychiatric disorders.

An Informatics Pipeline for Profiling and Annotating RNA Modifications

While over 150 distinct types of chemical modifications are known to occur on various cellular RNAs and can be dynamically controlled, the function of most of these modifications remains poorly defined. Collectively, these RNA modifications have been recently termed the "epitranscriptome". Identification and annotation of individual RNA modifications throughout the transcriptome are key for studying the role of the epitranscriptome in the regulation of gene expression and for elucidating the functional relevance of particular RNA modifications in diverse physiological and disease processes. In this protocol, we demonstrate how to identify and annotate RNA modifications based on the informatic analysis of methylated RNA immunoprecipitation and sequencing (MeRIP-seq) data, using RNAmod, a convenient one-stop online interactive platform for the annotation, analysis, and visualization of mRNA modifications.

Retinoblastoma tumor cell proliferation is negatively associated with an immune gene expression signature and increased immune cells

This study focuses on gene expression differences between early retinal states that ultimately lead to normal development, late onset retinoblastoma, or rapid bilateral retinoblastoma tumors. The late-onset and early-onset retinoblastoma tumor cells are remarkably similar to normally proliferating retinal progenitor cells, but they fail to properly express differentiation markers associated with normal development. Further, early-onset retinoblastoma tumor cells express a robust immune gene expression signature followed by accumulation of dendritic, monocyte, macrophage, and T-lymphocyte cells in the retinoblastoma tumors. This characteristic was not shared by either normal retinae or late-onset retinoblastomas. Comparison of our data with other human and mouse retinoblastoma tumor gene expression significantly confirmed, that the immune signature is present in tumors from each species. Strikingly, we observed that the immune signature in both mouse and human tumors was most highly evident in those with the lowest proliferative capacity. We directly assessed this relationship in human retinoblastoma tumors by co-analyzing proliferation and immune cell recruitment by immunohistochemistry, uncovering a significant inverse relationship between increased immune-cell infiltration in tumors and reduced tumor cell proliferation. Directly inhibiting proliferation with a PI3K/mTOR inhibitor significantly increased the number of CD45+ immune cells in the retina. This work establishes an in vivo model for the rapid recruitment of immune cells to tumorigenic neural tissue.

A high-resolution protein architecture of the budding yeast genome

The genome-wide architecture of chromatin-associated proteins that maintains chromosome integrity and gene regulation is not well defined. Here we use chromatin immunoprecipitation, exonuclease digestion and DNA sequencing (ChIP-exo/seq)1,2 to define this architecture in Saccharomyces cerevisiae. We identify 21 meta-assemblages consisting of roughly 400 different proteins that are related to DNA replication, centromeres, subtelomeres, transposons and transcription by RNA polymerase (Pol) I, II and III. Replication proteins engulf a nucleosome, centromeres lack a nucleosome, and repressive proteins encompass three nucleosomes at subtelomeric X-elements. We find that most promoters associated with Pol II evolved to lack a regulatory region, having only a core promoter. These constitutive promoters comprise a short nucleosome-free region (NFR) adjacent to a +1 nucleosome, which together bind the transcription-initiation factor TFIID to form a preinitiation complex. Positioned insulators protect core promoters from upstream events. A small fraction of promoters evolved an architecture for inducibility, whereby sequence-specific transcription factors (ssTFs) create a nucleosome-depleted region (NDR) that is distinct from an NFR. We describe structural interactions among ssTFs, their cognate cofactors and the genome. These interactions include the nucleosomal and transcriptional regulators RPD3-L, SAGA, NuA4, Tup1, Mediator and SWI-SNF. Surprisingly, we do not detect interactions between ssTFs and TFIID, suggesting that such interactions do not stably occur. Our model for gene induction involves ssTFs, cofactors and general factors such as TBP and TFIIB, but not TFIID. By contrast, constitutive transcription involves TFIID but not ssTFs engaged with their cofactors. From this, we define a highly integrated network of gene regulation by ssTFs.

Biology and genetic diversity of phasey bean mild yellows virus, a common virus in legumes in Australia

This study examined the natural and experimental host range and aphid and graft transmission of the tentative polerovirus phasey bean mild yellows virus (PBMYV). Eleven complete coding sequences from PBMYV isolates were determined from a range of hosts and locations. We found two genetically distinct variants of PBMYV. PBMYV-1 was the originally described variant, and PBMYV-2 had a large putative recombination in open reading frame 5 such that PBMYV-1 and PBMYV-2 shared only 65-66% amino acid sequence identity in the P5 protein. The virus was transmitted by a clonal colony of cowpea aphids (Aphis craccivora) and by grafting with infected scions but was not transmitted by a clonal colony of green peach aphids (Myzus persicae). PBMYV was found in natural infections in 11 host species with a range of symptoms and severity, including seven important grain legume crops from across a wide geographic area in Australia. PBMYV was common and widespread in the tropical weed phasey bean (Macroptilium lathyroides), but it is likely that there are other major alternative hosts for the virus in temperate regions of Australia. The experimental host range of PBMYV included the Fabaceae hosts chickpea (Cicer arietinum), faba bean (Vicia faba), pea (Pisum sativum), and phasey bean, but transmissions failed to infect several other members of the families Asteraceae, Cucurbitaceae, Fabaceae and Solanaceae. PBMYV was commonly found in grain legume crops in eastern and western Australia, sometimes at greater than 90% incidence. This new knowledge about PBMYV warrants further assessments of its economic impact on important grain legume crops.

Subcutaneous and intramuscular fat transcriptomes show large differences in network organization and associations with adipose traits in pigs

Subcutaneous fat (SCF) and intramuscular fat (IMF) deposition is relevant to health in humans, as well as meat production and quality in pigs. In this study, we generated RNA sequence data for 122 SCF, 120 IMF, and 87 longissimus dorsi muscle (LDM) samples using 155 F₆ pigs from a specially designed heterogeneous population generated by intercrossing four highly selected European commercial breeds and four indigenous Chinese pig breeds. The phenotypes including waist back fat thickness and intramuscular fat content were also measured in the 155 F₆ pigs. We found that the genes in SCF and IMF differed largely in both expression levels and network connectivity, and highlighted network modules that exhibited strongest gain of connectivity in SCF and IMF, containing genes that were associated with the immune process and DNA double-strand repair, respectively. We identified 215 SCF genes related to kinase inhibitor activity, mitochondrial fission, and angiogenesis, and 90 IMF genes related to lipolysis and fat cell differentiation, displayed a tissue-specific association with back fat thickness and IMF content, respectively. We found that cis-expression QTL for trait-associated genes in the two adipose tissues tended to have tissue-dependent predictability for the two adipose traits. Alternative splicing of genes was also found to be associated with SCF or IMF deposition, but the association was much less extensive than that based on expression levels. This study provides a better understanding of SCF and IMF gene transcription and network organization and identified critical genes and network modules that displayed tissue-specific associations with subcutaneous and intramuscular fat deposition. These features are helpful for designing breeding programs to genetically improve the two adipose traits in a balanced way.

A novel ACE2 isoform is expressed in human respiratory epithelia and is upregulated in response to interferons and RNA respiratory virus infection

Angiotensin-converting enzyme 2 (ACE2) is the main entry point in airway epithelial cells for SARS-CoV-2. ACE2 binding to the SARS-CoV-2 protein spike triggers viral fusion with the cell plasma membrane, resulting in viral RNA genome delivery into the host. Despite ACE2's critical role in SARS-CoV-2 infection, full understanding of ACE2 expression, including in response to viral infection, remains unclear. ACE2 was thought to encode five transcripts and one protein of 805 amino acids. In the present study, we identify a novel short isoform of ACE2 expressed in the airway epithelium, the main site of SARS-CoV-2 infection. Short ACE2 is substantially upregulated in response to interferon stimulation and rhinovirus infection, but not SARS-CoV-2 infection. This short isoform lacks SARS-CoV-2 spike high-affinity binding sites and, altogether, our data are consistent with a model where short ACE2 is unlikely to directly contribute to host susceptibility to SARS-CoV-2 infection.

The androgen receptor is a tumor suppressor in estrogen receptor-positive breast cancer

The role of the androgen receptor (AR) in estrogen receptor (ER)-α-positive breast cancer is controversial, constraining implementation of AR-directed therapies. Using a diverse, clinically relevant panel of cell-line and patient-derived models, we demonstrate that AR activation, not suppression, exerts potent antitumor activity in multiple disease contexts, including resistance to standard-of-care ER and CDK4/6 inhibitors. Notably, AR agonists combined with standard-of-care agents enhanced therapeutic responses. Mechanistically, agonist activation of AR altered the genomic distribution of ER and essential co-activators (p300, SRC-3), resulting in repression of ER-regulated cell cycle genes and upregulation of AR target genes, including known tumor suppressors. A gene signature of AR activity positively predicted disease survival in multiple clinical ER-positive breast cancer cohorts. These findings provide unambiguous evidence that AR has a tumor suppressor role in ER-positive breast cancer and support AR agonism as the optimal AR-directed treatment strategy, revealing a rational therapeutic opportunity.

Machine learning reveals bilateral distribution of somatic L1 insertions in human neurons and glia

Retrotransposons can cause somatic genome variation in the human nervous system, which is hypothesized to have relevance to brain development and neuropsychiatric disease. However, the detection of individual somatic mobile element insertions presents a difficult signal-to-noise problem. Using a machine-learning method (RetroSom) and deep whole-genome sequencing, we analyzed L1 and Alu retrotransposition in sorted neurons and glia from human brains. We characterized two brain-specific L1 insertions in neurons and glia from a donor with schizophrenia. There was anatomical distribution of the L1 insertions in neurons and glia across both hemispheres, indicating retrotransposition occurred during early embryogenesis. Both insertions were within the introns of genes (CNNM2 and FRMD4A) inside genomic loci associated with neuropsychiatric disorders. Proof-of-principle experiments revealed these L1 insertions significantly reduced gene expression. These results demonstrate that RetroSom has broad applications for studies of brain development and may provide insight into the possible pathological effects of somatic retrotransposition.

The T-box transcription factor Eomesodermin governs haemogenic competence of yolk sac mesodermal progenitors

Extra-embryonic mesoderm (ExM)-composed of the earliest cells that traverse the primitive streak-gives rise to the endothelium as well as haematopoietic progenitors in the developing yolk sac. How a specific subset of ExM becomes committed to a haematopoietic fate remains unclear. Here we demonstrate using an embryonic stem cell model that transient expression of the T-box transcription factor Eomesodermin (Eomes) governs haemogenic competency of ExM. Eomes regulates the accessibility of enhancers that the transcription factor stem cell leukaemia (SCL) normally utilizes to specify primitive erythrocytes and is essential for the normal development of Runx1+ haemogenic endothelium. Single-cell RNA sequencing suggests that Eomes loss of function profoundly blocks the formation of blood progenitors but not specification of Flk-1+ haematoendothelial progenitors. Our findings place Eomes at the top of the transcriptional hierarchy regulating early blood formation and suggest that haemogenic competence is endowed earlier during embryonic development than was previously appreciated.

Identification and characterization of enzymes involved in the biosynthesis of pyrimidine nucleoside antibiotics

Covering: up to September 2020 Hundreds of nucleoside-based natural products have been isolated from various microorganisms, several of which have been utilized in agriculture as pesticides and herbicides, in medicine as therapeutics for cancer and infectious disease, and as molecular probes to study biological processes. Natural products consisting of structural modifications of each of the canonical nucleosides have been discovered, ranging from simple modifications such as single-step alkylations or acylations to highly elaborate modifications that dramatically alter the nucleoside scaffold and require multiple enzyme-catalyzed reactions. A vast amount of genomic information has been uncovered the past two decades, which has subsequently allowed the first opportunity to interrogate the chemically intriguing enzymatic transformations for the latter type of modifications. This review highlights (i) the discovery and potential applications of structurally complex pyrimidine nucleoside antibiotics for which genetic information is known, (ii) the established reactions that convert the canonical pyrimidine into a new nucleoside scaffold, and (iii) the important tailoring reactions that impart further structural complexity to these molecules.

Bottleneck size and selection level reproducibly impact evolution of antibiotic resistance

During antibiotic treatment, the evolution of bacterial pathogens is fundamentally affected by bottlenecks and varying selection levels imposed by the drugs. Bottlenecks-that is, reductions in bacterial population size-lead to an increased influence of random effects (genetic drift) during bacterial evolution, and varying antibiotic concentrations during treatment may favour distinct resistance variants. Both aspects influence the process of bacterial evolution during antibiotic therapy and thereby treatment outcome. Surprisingly, the joint influence of these interconnected factors on the evolution of antibiotic resistance remains largely unexplored. Here we combine evolution experiments with genomic and genetic analyses to demonstrate that bottleneck size and antibiotic-induced selection reproducibly impact the evolutionary path to resistance in pathogenic Pseudomonas aeruginosa, one of the most problematic opportunistic human pathogens. Resistance is favoured-expectedly-under high antibiotic selection and weak bottlenecks, but-unexpectedly-also under low antibiotic selection and severe bottlenecks. The latter is likely to result from a reduced probability of losing favourable variants through drift under weak selection. Moreover, the absence of high resistance under low selection and weak bottlenecks is caused by the spread of low-resistance variants with high competitive fitness under these conditions. We conclude that bottlenecks, in combination with drug-induced selection, are currently neglected key determinants of pathogen evolution and outcome of antibiotic treatment.

Mutations in circulating tumor DNA predict primary resistance to systemic therapies in advanced hepatocellular carcinoma

Little is known about the mutational landscape of advanced hepatocellular carcinoma (HCC), and predictive biomarkers of response to systemic therapies are lacking. We aimed to describe the mutational landscape of advanced HCC and to identify predictors of primary resistance to systemic therapies using circulating tumor DNA (ctDNA). We prospectively enrolled 121 patients between October 2015 and January 2019. We performed targeted ultra-deep sequencing of 25 genes and Digital Droplet PCR of TERT promoter, including sequential samples throughout treatment. Primary endpoint was progression-free survival (PFS) stratified by mutation profiles in ctDNA. Secondary endpoints were overall survival and objective response rate. The most frequent mutations in ctDNA of advanced HCC were TERT promoter (51%), TP53 (32%), CTNNB1 (17%), PTEN (8%), AXIN1, ARID2, KMT2D, and TSC2 (each 6%). TP53 and CTNNB1 mutations were mutually exclusive. Patients with mutations in the PI3K/MTOR pathway had significantly shorter PFS than those without these mutations after tyrosine kinase inhibitors (2.1 vs 3.7 months, p < 0.001), but not after immune checkpoint inhibition (CPI). WNT pathway mutations were not associated with PFS, overall survival, or objective response after CPI. Serial profiling of ctDNA in a subset correlated with treatment response. Mutation profiling of ctDNA in advanced HCC shows similar mutation frequencies for known HCC drivers compared to early stages and identifies predictive biomarkers of response to systemic therapies.

Major driver mutations are shared between sinonasal intestinal-type adenocarcinoma and the morphologically identical colorectal adenocarcinoma

PURPOSE

The purpose of our study was to compare genomic changes in sinonasal intestinal-type adenocarcinoma (sITAC) and colorectal adenocarcinoma (CRC), as they are histomorphologically indistinguishable. This can cause diagnostic difficulties as sinonasal tumours initially diagnosed as sITAC may represent metastasis from CRC, a frequent cancer. Previous studies have not uncovered the underlying mechanism behind the histomorphological resemblance.

METHODS/PATIENTS

Tissue samples from all consecutive patients with sITAC at our facility (20 patients) were compared to samples from 20 patients with CRC as well as samples from 2 patients with both CRC and sinonasal tumours. DNA sequencing was performed using Illumina TruSight Oncology 500 panel consisting of 523 cancer-associated genes. Frequent mutations were inspected manually using the Integrative Genomics Viewer.

RESULTS

Several well-known cancer-associated genes were mutated in the CRC group, but also in the sinonasal ITAC group. These genes included APC mutated in 65% of the CRC group and 37% of the sinonasal ITAC group, and TP53 mutated in 65% of CRC samples and 58% of ITAC samples. These shared mutations may explain the histomorphological similarities. Successful DNA sequencing was performed on the colorectal sample from one of the two patients with both CRC and sinonasal tumour. Comparing mutations in these samples from one patient we have shown that the sinonasal tumour in all probability was a CRC metastasis.

CONCLUSION

We have identified several genetic similarities between sITAC and CRC. This discovery brings us closer to understanding mechanisms behind the development of sITAC-and hopefully in the future targeted therapy.

Repeat expansions confer WRN dependence in microsatellite-unstable cancers

The RecQ DNA helicase WRN is a synthetic lethal target for cancer cells with microsatellite instability (MSI), a form of genetic hypermutability that arises from impaired mismatch repair1-4. Depletion of WRN induces widespread DNA double-strand breaks in MSI cells, leading to cell cycle arrest and/or apoptosis. However, the mechanism by which WRN protects MSI-associated cancers from double-strand breaks remains unclear. Here we show that TA-dinucleotide repeats are highly unstable in MSI cells and undergo large-scale expansions, distinct from previously described insertion or deletion mutations of a few nucleotides5. Expanded TA repeats form non-B DNA secondary structures that stall replication forks, activate the ATR checkpoint kinase, and require unwinding by the WRN helicase. In the absence of WRN, the expanded TA-dinucleotide repeats are susceptible to cleavage by the MUS81 nuclease, leading to massive chromosome shattering. These findings identify a distinct biomarker that underlies the synthetic lethal dependence on WRN, and support the development of therapeutic agents that target WRN for MSI-associated cancers.

Single-gene imaging links genome topology, promoter-enhancer communication and transcription control

Transcription activation by distal enhancers is essential for cell-fate specification and maintenance of cellular identities. How long-range gene regulation is physically achieved, especially within complex regulatory landscapes of non-binary enhancer-promoter configurations, remains elusive. Recent nanoscopy advances have quantitatively linked promoter kinetics and ~100- to 200-nm-sized clusters of enhancer-associated regulatory factors (RFs) at important developmental genes. Here, we further dissect mechanisms of RF clustering and transcription activation in mouse embryonic stem cells. RF recruitment into clusters involves specific molecular recognition of cognate DNA and chromatin-binding sites, suggesting underlying cis-element clustering. Strikingly, imaging of tagged genomic loci, with ≤1 kilobase and ~20-nanometer precision, in live cells, reveals distal enhancer clusters over the extended locus in frequent close proximity to target genes-within RF-clustering distances. These high-interaction-frequency enhancer-cluster 'superclusters' create nano-environments wherein clustered RFs activate target genes, providing a structural framework for relating genome organization, focal RF accumulation and transcription activation.

Antibiotic susceptibility signatures identify potential antimicrobial targets in the Acinetobacter baumannii cell envelope

A unique, protective cell envelope contributes to the broad drug resistance of the nosocomial pathogen Acinetobacter baumannii. Here we use transposon insertion sequencing to identify A. baumannii mutants displaying altered susceptibility to a panel of diverse antibiotics. By examining mutants with antibiotic susceptibility profiles that parallel mutations in characterized genes, we infer the function of multiple uncharacterized envelope proteins, some of which have roles in cell division or cell elongation. Remarkably, mutations affecting a predicted cell wall hydrolase lead to alterations in lipooligosaccharide synthesis. In addition, the analysis of altered susceptibility signatures and antibiotic-induced morphology patterns allows us to predict drug synergies; for example, certain beta-lactams appear to work cooperatively due to their preferential targeting of specific cell wall assembly machineries. Our results indicate that the pathogen may be effectively inhibited by the combined targeting of multiple pathways critical for envelope growth.

Endogenous retroviruses drive species-specific germline transcriptomes in mammals

Gene regulation in the germline ensures the production of high-quality gametes, long-term maintenance of the species, and speciation. Male germline transcriptomes undergo dynamic changes after the mitosis-to-meiosis transition and have been subject to evolutionary divergence among mammals. However, the mechanisms underlying germline regulatory divergence remain undetermined. Here, we show that endogenous retroviruses (ERVs) influence species-specific germline transcriptomes. After the mitosis-to-meiosis transition in male mice, specific ERVs function as active enhancers to drive germline genes, including a mouse-specific gene set, and bear binding motifs for critical regulators of spermatogenesis such as A-MYB. This raises the possibility that a genome-wide transposition of ERVs rewired germline gene expression in a species-specific manner. Of note, independently evolved ERVs are associated with the expression of human-specific germline genes, demonstrating the prevalence of ERV-driven mechanisms in mammals. Together, we propose that ERVs fine-tune species-specific transcriptomes in the mammalian germline.

The ubiquitin ligase Cullin-1 associates with chromatin and regulates transcription of specific c-MYC target genes

Transcription is regulated through a dynamic interplay of DNA-associated proteins, and the composition of gene-regulatory complexes is subject to continuous adjustments. Protein alterations include post-translational modifications and elimination of individual polypeptides. Spatially and temporally controlled protein removal is, therefore, essential for gene regulation and accounts for the short half-life of many transcription factors. The ubiquitin-proteasome system is responsible for site- and target-specific ubiquitination and protein degradation. Specificity of ubiquitination is conferred by ubiquitin ligases. Cullin-RING complexes, the largest family of ligases, require multi-unit assembly around one of seven cullin proteins. To investigate the direct role of cullins in ubiquitination of DNA-bound proteins and in gene regulation, we analyzed their subcellular locations and DNA-affinities. We found CUL4A and CUL7 to be largely excluded from the nucleus, whereas CUL4B was primarily nuclear. CUL1,2,3, and 5 showed mixed cytosolic and nuclear expression. When analyzing chromatin affinity of individual cullins, we discovered that CUL1 preferentially associated with active promoter sequences and co-localized with 23% of all DNA-associated protein degradation sites. CUL1 co-distributed with c-MYC and specifically repressed nuclear-encoded mitochondrial and splicing-associated genes. These studies underscore the relevance of spatial control in chromatin-associated protein ubiquitination and define a novel role for CUL1 in gene repression.

Genomic characterization of malignant progression in neoplastic pancreatic cysts

Intraductal papillary mucinous neoplasms (IPMNs) and mucinous cystic neoplasms (MCNs) are non-invasive neoplasms that are often observed in association with invasive pancreatic cancers, but their origins and evolutionary relationships are poorly understood. In this study, we analyze 148 samples from IPMNs, MCNs, and small associated invasive carcinomas from 18 patients using whole exome or targeted sequencing. Using evolutionary analyses, we establish that both IPMNs and MCNs are direct precursors to pancreatic cancer. Mutations in SMAD4 and TGFBR2 are frequently restricted to invasive carcinoma, while RNF43 alterations are largely in non-invasive lesions. Genomic analyses suggest an average window of over three years between the development of high-grade dysplasia and pancreatic cancer. Taken together, these data establish non-invasive IPMNs and MCNs as origins of invasive pancreatic cancer, identifying potential drivers of invasion, highlighting the complex clonal dynamics prior to malignant transformation, and providing opportunities for early detection and intervention.

Developmental changes in the accessible chromatin, transcriptome and Ascl1-binding correlate with the loss in Müller Glial regenerative potential

Diseases and damage to the retina lead to losses in retinal neurons and eventual visual impairment. Although the mammalian retina has no inherent regenerative capabilities, fish have robust regeneration from Müller glia (MG). Recently, we have shown that driving expression of Ascl1 in adult mouse MG stimulates neural regeneration. The regeneration observed in the mouse is limited in the variety of neurons that can be derived from MG; Ascl1-expressing MG primarily generate bipolar cells. To better understand the limits of MG-based regeneration in mouse retinas, we used ATAC- and RNA-seq to compare newborn progenitors, immature MG (P8-P12), and mature MG. Our analysis demonstrated developmental differences in gene expression and accessible chromatin between progenitors and MG, primarily in neurogenic genes. Overexpression of Ascl1 is more effective in reprogramming immature MG, than mature MG, consistent with a more progenitor-like epigenetic landscape in the former. We also used ASCL1 ChIPseq to compare the differences in ASCL1 binding in progenitors and reprogrammed MG. We find that bipolar-specific accessible regions are more frequently linked to bHLH motifs and ASCL1 binding. Overall, our analysis indicates a loss of neurogenic gene expression and motif accessibility during glial maturation that may prevent efficient reprogramming.

Chd4 choreographs self-antigen expression for central immune tolerance

Autoreactive T cells are eliminated in the thymus to prevent autoimmunity by promiscuous expression of tissue-restricted self-antigens in medullary thymic epithelial cells. This expression is dependent on the transcription factor Fezf2, as well as the transcriptional regulator Aire, but the entire picture of the transcriptional program has been obscure. Here, we found that the chromatin remodeler Chd4, also called Mi-2β, plays a key role in the self-antigen expression in medullary thymic epithelial cells. To maximize the diversity of self-antigen expression, Fezf2 and Aire utilized completely distinct transcriptional mechanisms, both of which were under the control of Chd4. Chd4 organized the promoter regions of Fezf2-dependent genes, while contributing to the Aire-mediated induction of self-antigens via super-enhancers. Mice deficient in Chd4 specifically in thymic epithelial cells exhibited autoimmune phenotypes, including T cell infiltration. Thus, Chd4 plays a critical role in integrating Fezf2- and Aire-mediated gene induction to establish central immune tolerance.

Dual RNA-seq of Orientia tsutsugamushi informs on host-pathogen interactions for this neglected intracellular human pathogen

Studying emerging or neglected pathogens is often challenging due to insufficient information and absence of genetic tools. Dual RNA-seq provides insights into host-pathogen interactions, and is particularly informative for intracellular organisms. Here we apply dual RNA-seq to Orientia tsutsugamushi (Ot), an obligate intracellular bacterium that causes the vector-borne human disease scrub typhus. Half the Ot genome is composed of repetitive DNA, and there is minimal collinearity in gene order between strains. Integrating RNA-seq, comparative genomics, proteomics, and machine learning to study the transcriptional architecture of Ot, we find evidence for wide-spread post-transcriptional antisense regulation. Comparing the host response to two clinical isolates, we identify distinct immune response networks for each strain, leading to predictions of relative virulence that are validated in a mouse infection model. Thus, dual RNA-seq can provide insight into the biology and host-pathogen interactions of a poorly characterized and genetically intractable organism such as Ot.

Laboratory culture of the California Sea Firefly Vargula tsujii (Ostracoda: Cypridinidae): Developing a model system for the evolution of marine bioluminescence

Bioluminescence, or the production of light by living organisms via chemical reaction, is widespread across Metazoa. Laboratory culture of bioluminescent organisms from diverse taxonomic groups is important for determining the biosynthetic pathways of bioluminescent substrates, which may lead to new tools for biotechnology and biomedicine. Some bioluminescent groups may be cultured, including some cnidarians, ctenophores, and brittle stars, but those use luminescent substrates (luciferins) obtained from their diets, and therefore are not informative for determination of the biosynthetic pathways of the luciferins. Other groups, including terrestrial fireflies, do synthesize their own luciferin, but culturing them is difficult and the biosynthetic pathway for firefly luciferin remains unclear. An additional independent origin of endogenous bioluminescence is found within ostracods from the family Cypridinidae, which use their luminescence for defense and, in Caribbean species, for courtship displays. Here, we report the first complete life cycle of a luminous ostracod (Vargula tsujii Kornicker & Baker, 1977, the California Sea Firefly) in the laboratory. We also describe the late-stage embryogenesis of Vargula tsujii and discuss the size classes of instar development. We find embryogenesis in V. tsujii ranges from 25–38 days, and this species appears to have five instar stages, consistent with ontogeny in other cypridinid lineages. We estimate a complete life cycle at 3–4 months. We also present the first complete mitochondrial genome for Vargula tsujii. Bringing a luminous ostracod into laboratory culture sets the stage for many potential avenues of study, including learning the biosynthetic pathway of cypridinid luciferin and genomic manipulation of an autogenic bioluminescent system.

Dissecting myogenin-mediated retinoid X receptor signaling in myogenic differentiation

Deciphering the molecular mechanisms underpinning myoblast differentiation is a critical step in developing the best strategy to promote muscle regeneration in patients suffering from muscle-related diseases. We have previously established that a rexinoid x receptor (RXR)-selective agonist, bexarotene, enhances the differentiation and fusion of myoblasts through a direct regulation of MyoD expression, coupled with an augmentation of myogenin protein. Here, we found that RXR signaling associates with the distribution of myogenin at poised enhancers and a distinct E-box motif. We also found an association of myogenin with rexinoid-responsive gene expression and identified an epigenetic signature related to histone acetyltransferase p300. Moreover, RXR signaling augments residue-specific histone acetylation at enhancers co-occupied by p300 and myogenin. Thus, genomic distribution of transcriptional regulators is an important designate for identifying novel targets as well as developing therapeutics that modulate epigenetic landscape in a selective manner to promote muscle regeneration.

Recurrent horizontal transfer identifies mitochondrial positive selection in a transmissible cancer

Autonomous replication and segregation of mitochondrial DNA (mtDNA) creates the potential for evolutionary conflict driven by emergence of haplotypes under positive selection for 'selfish' traits, such as replicative advantage. However, few cases of this phenomenon arising within natural populations have been described. Here, we survey the frequency of mtDNA horizontal transfer within the canine transmissible venereal tumour (CTVT), a contagious cancer clone that occasionally acquires mtDNA from its hosts. Remarkably, one canine mtDNA haplotype, A1d1a, has repeatedly and recently colonised CTVT cells, recurrently replacing incumbent CTVT haplotypes. An A1d1a control region polymorphism predicted to influence transcription is fixed in the products of an A1d1a recombination event and occurs somatically on other CTVT mtDNA backgrounds. We present a model whereby 'selfish' positive selection acting on a regulatory variant drives repeated fixation of A1d1a within CTVT cells.

Targeted reprogramming of H3K27me3 resets epigenetic memory in plant paternal chromatin

Epigenetic marks are reprogrammed in the gametes to reset genomic potential in the next generation. In mammals, paternal chromatin is extensively reprogrammed through the global erasure of DNA methylation and the exchange of histones with protamines1,2. Precisely how the paternal epigenome is reprogrammed in flowering plants has remained unclear since DNA is not demethylated and histones are retained in sperm3,4. Here, we describe a multi-layered mechanism by which H3K27me3 is globally lost from histone-based sperm chromatin in Arabidopsis. This mechanism involves the silencing of H3K27me3 writers, activity of H3K27me3 erasers and deposition of a sperm-specific histone, H3.10 (ref. 5), which we show is immune to lysine 27 methylation. The loss of H3K27me3 facilitates the transcription of genes essential for spermatogenesis and pre-configures sperm with a chromatin state that forecasts gene expression in the next generation. Thus, plants have evolved a specific mechanism to simultaneously differentiate male gametes and reprogram the paternal epigenome.

Consistent RNA sequencing contamination in GTEx and other data sets

A challenge of next generation sequencing is read contamination. We use Genotype-Tissue Expression (GTEx) datasets and technical metadata along with RNA-seq datasets from other studies to understand factors that contribute to contamination. Here we report, of 48 analyzed tissues in GTEx, 26 have variant co-expression clusters of four highly expressed and pancreas-enriched genes (PRSS1, PNLIP, CLPS, and/or CELA3A). Fourteen additional highly expressed genes from other tissues also indicate contamination. Sample contamination is strongly associated with a sample being sequenced on the same day as a tissue that natively expresses those genes. Discrepant SNPs across four contaminating genes validate the contamination. Low-level contamination affects ~40% of samples and leads to numerous eQTL assignments in inappropriate tissues among these 18 genes. This type of contamination occurs widely, impacting bulk and single cell (scRNA-seq) data set analysis. In conclusion, highly expressed, tissue-enriched genes basally contaminate GTEx and other datasets impacting analyses.

Chromatin architecture reorganization in murine somatic cell nuclear transfer embryos

The oocyte cytoplasm can reprogram the somatic cell nucleus into a totipotent state, but with low efficiency. The spatiotemporal chromatin organization of somatic cell nuclear transfer (SCNT) embryos remains elusive. Here, we examine higher order chromatin structures of mouse SCNT embryos using a low-input Hi-C method. We find that donor cell chromatin transforms to the metaphase state rapidly after SCNT along with the dissolution of typical 3D chromatin structure. Intriguingly, the genome undergoes a mitotic metaphase-like to meiosis metaphase II-like transition following activation. Subsequently, weak chromatin compartments and topologically associating domains (TADs) emerge following metaphase exit. TADs are further removed until the 2-cell stage before being progressively reestablished. Obvious defects including stronger TAD boundaries, aberrant super-enhancer and promoter interactions are found in SCNT embryos. These defects are partially caused by inherited H3K9me3, and can be rescued by Kdm4d overexpression. These observations provide insight into chromatin architecture reorganization during SCNT embryo development.

The organisation of chromatin in somatic cell nuclear transfer (SCNT) embryos remains poorly understood. Here, the authors examine higher order chromatin structures of mouse SCNT embryos and provide insights into chromatin architecture reorganisation during SCNT embryo development.

Merkel cell polyomavirus activates LSD1-mediated blockade of non-canonical BAF to regulate transformation and tumorigenesis

Merkel cell carcinoma (MCC), a neuroendocrine cancer of the skin, is caused by integration of Merkel cell polyomavirus (MCV) and persistent expression of Large T antigen (LT) and Small T antigen (ST). We report that ST in complex with MYCL and the EP400 complex activates expression of LSD1 (KDM1A), RCOR2, and INSM1 to repress gene expression by the lineage transcription factor ATOH1. LSD1 inhibition reduces growth of MCC in vitro and in vivo. Through a forward-genetics CRISPR-Cas9 screen, we identified an antagonistic relationship between LSD1 and the non-canonical BAF (ncBAF) chromatin remodeling complex. Changes in gene expression and chromatin accessibility caused by LSD1 inhibition could be partially rescued by BRD9 inhibition, revealing that LSD1 and ncBAF antagonistically regulate an overlapping set of genes. Our work provides mechanistic insight into the dependence of MCC on LSD1 and a tumor suppressor role for ncBAF in cancer.

Exome sequencing of familial high-grade serous ovarian carcinoma reveals heterogeneity for rare candidate susceptibility genes

High-grade serous ovarian carcinoma (HGSOC) has a significant hereditary component, approximately half of which cannot be explained by known genes. To discover genes, we analyse germline exome sequencing data from 516 BRCA1/2-negative women with HGSOC, focusing on genes enriched with rare, protein-coding loss-of-function (LoF) variants. Overall, there is a significant enrichment of rare protein-coding LoF variants in the cases (p < 0.0001, chi-squared test). Only thirty-four (6.6%) have a pathogenic variant in a known or proposed predisposition gene. Few genes have LoF mutations in more than four individuals and the majority are detected in one individual only. Forty-three highly-ranked genes are identified with three or more LoF variants that are enriched by three-fold or more compared to GnomAD. These genes represent diverse functional pathways with relatively few involved in DNA repair, suggesting that much of the remaining heritability is explained by previously under-explored genes and pathways.

Characterization of novel genetic alterations in salivary gland secretory carcinoma

Secretory carcinoma is a salivary gland tumor with a characteristic chromosomal translocation that results in an ETV6-NTRK3 fusion gene. Secretory carcinoma shows relatively frequent rates of lymph-node metastasis and tumor recurrence and has a characteristic histology. Except for the ETV6 translocation, genomic alterations in secretory carcinoma have not been reported. In the present study, we characterized the novel recurrent genetic mutations of secretory carcinoma. On the basis of histology, immunohistochemistry, and ETV6 gene break-apart fluorescence in situ hybridization assays, 22 tumors were classified as secretory carcinomas (19 ETV6 translocation-positive and 3 ETV6 translocation-negative secretory carcinomas) and their clinicopathologic characteristics were reviewed. Targeted deep sequencing analyses were performed on 20 secretory carcinomas (17 ETV6 translocation-positive and 3 ETV6 translocation-negative secretory carcinomas) to investigate their genetic alterations. The A16V (C→T) mutation in PRSS1, which encodes a cationic trypsinogen and has a mutation associated with hereditary pancreatitis and pancreatic adenocarcinoma, was observed in 40% (8/20) (7/17 of ETV6 translocation-positive and 1/3 of ETV6 translocation-negative secretory carcinomas). Pathogenic variants of MLH1, MUTYH, and STK11 were also identified. Variants of uncertain significance included mutations in KMT5A. These novel characteristic genetic alterations may advance current understandings of secretory carcinoma tumorigenesis and progression, leading to improved diagnoses and treatment strategies.

HIV-1-induced cytokines deplete homeostatic innate lymphoid cells and expand TCF7-dependent memory NK cells

Human immunodeficiency virus 1 (HIV-1) infection is associated with heightened inflammation and excess risk of cardiovascular disease, cancer and other complications. These pathologies persist despite antiretroviral therapy. In two independent cohorts, we found that innate lymphoid cells (ILCs) were depleted in the blood and gut of people with HIV-1, even with effective antiretroviral therapy. ILC depletion was associated with neutrophil infiltration of the gut lamina propria, type 1 interferon activation, increased microbial translocation and natural killer (NK) cell skewing towards an inflammatory state, with chromatin structure and phenotype typical of WNT transcription factor TCF7-dependent memory T cells. Cytokines that are elevated during acute HIV-1 infection reproduced the ILC and NK cell abnormalities ex vivo. These results show that inflammatory cytokines associated with HIV-1 infection irreversibly disrupt ILCs. This results in loss of gut epithelial integrity, microbial translocation and memory NK cells with heightened inflammatory potential, and explains the chronic inflammation in people with HIV-1.

Engineering bacterial symbionts of nematodes improves their biocontrol potential to counter the western corn rootworm

The western corn rootworm (WCR) decimates maize crops worldwide. One potential way to control this pest is treatment with entomopathogenic nematodes (EPNs) that harbor bacterial symbionts that are pathogenic to insects. However, WCR larvae sequester benzoxazinoid secondary metabolites that are produced by maize and use them to increase their resistance to the nematodes and their symbionts. Here we report that experimental evolution and selection for bacterial symbionts that are resistant to benzoxazinoids improve the ability of a nematode-symbiont pair to kill WCR larvae. We isolated five Photorhabdus symbionts from different nematodes and increased their benzoxazinoid resistance through experimental evolution. Benzoxazinoid resistance evolved through multiple mechanisms, including a mutation in the aquaporin-like channel gene aqpZ. We reintroduced benzoxazinoid-resistant Photorhabdus strains into their original EPN hosts and identified one nematode-symbiont pair that was able to kill benzoxazinoid-sequestering WCR larvae more efficiently. Our results suggest that modification of bacterial symbionts might provide a generalizable strategy to improve biocontrol of agricultural pests.

Mutations in ASH1L confer susceptibility to Tourette syndrome

Tourette syndrome (TS) is a childhood-onset neuropsychiatric disorder characterized by repetitive motor movements and vocal tics. The clinical manifestations of TS are complex and often overlap with other neuropsychiatric disorders. TS is highly heritable; however, the underlying genetic basis and molecular and neuronal mechanisms of TS remain largely unknown. We performed whole-exome sequencing of a hundred trios (probands and their parents) with detailed records of their clinical presentations and identified a risk gene, ASH1L, that was both de novo mutated and associated with TS based on a transmission disequilibrium test. As a replication, we performed follow-up targeted sequencing of ASH1L in additional 524 unrelated TS samples and replicated the association (P value = 0.001). The point mutations in ASH1L cause defects in its enzymatic activity. Therefore, we established a transgenic mouse line and performed an array of anatomical, behavioral, and functional assays to investigate ASH1L function. The Ash1l^+/- mice manifested tic-like behaviors and compulsive behaviors that could be rescued by the tic-relieving drug haloperidol. We also found that Ash1l disruption leads to hyper-activation and elevated dopamine-releasing events in the dorsal striatum, all of which could explain the neural mechanisms for the behavioral abnormalities in mice. Taken together, our results provide compelling evidence that ASH1L is a TS risk gene.

Exome sequencing in schizophrenia-affected parent-offspring trios reveals risk conferred by protein-coding de novo mutations

Protein-coding de novo mutations (DNMs) are significant risk factors in many neurodevelopmental disorders, whereas schizophrenia (SCZ) risk associated with DNMs has thus far been shown to be modest. We analyzed DNMs from 1,695 SCZ-affected trios and 1,077 published SCZ-affected trios to better understand the contribution to SCZ risk. Among 2,772 SCZ probands, exome-wide DNM burden remained modest. Gene set analyses revealed that SCZ DNMs were significantly concentrated in genes that were highly expressed in the brain, that were under strong evolutionary constraint and/or overlapped with genes identified in other neurodevelopmental disorders. No single gene surpassed exome-wide significance; however, 16 genes were recurrently hit by protein-truncating DNMs, corresponding to a 3.15-fold higher rate than the mutation model expectation (permuted 95% confidence interval: 1-10 genes; permuted P = 3 × 10-5). Overall, DNMs explain a small fraction of SCZ risk, and larger samples are needed to identify individual risk genes, as coding variation across many genes confers risk for SCZ in the population.

Mechanisms governing the pioneering and redistribution capabilities of the non-classical pioneer PU.1

Establishing gene regulatory networks during differentiation or reprogramming requires master or pioneer transcription factors (TFs) such as PU.1, a prototype master TF of hematopoietic lineage differentiation. To systematically determine molecular features that control its activity, here we analyze DNA-binding in vitro and genome-wide in vivo across different cell types with native or ectopic PU.1 expression. Although PU.1, in contrast to classical pioneer factors, is unable to access nucleosomal target sites in vitro, ectopic induction of PU.1 leads to the extensive remodeling of chromatin and redistribution of partner TFs. De novo chromatin access, stable binding, and redistribution of partner TFs both require PU.1's N-terminal acidic activation domain and its ability to recruit SWI/SNF remodeling complexes, suggesting that the latter may collect and distribute co-associated TFs in conjunction with the non-classical pioneer TF PU.1.

Control of human hemoglobin switching by LIN28B-mediated regulation of BCL11A translation

Increased production of fetal hemoglobin (HbF) can ameliorate the severity of sickle cell disease and β-thalassemia¹. BCL11A represses the genes encoding HbF and regulates human hemoglobin switching through variation in its expression during development^2–7. However, the mechanisms underlying the developmental expression of BCL11A remain mysterious. Here we show that BCL11A is regulated at the level of mRNA translation during human hematopoietic development. Despite decreased BCL11A protein synthesis earlier in development, BCL11A mRNA continues to be associated with ribosomes. Through unbiased genomic and proteomic analyses, we demonstrate that the RNA-binding protein LIN28B, which is developmentally expressed in a reciprocal pattern to BCL11A, directly interacts with ribosomes and BCL11A mRNA. Furthermore, we show that BCL11A mRNA translation is suppressed by LIN28B through direct interactions, independent of its role in regulating let-7 microRNAs, and BCL11A is the major target of LIN28B-mediated HbF induction. Our results reveal a previously unappreciated mechanism underlying human hemoglobin switching that illuminates new therapeutic opportunities.

A transcriptome-wide antitermination mechanism sustaining identity of embryonic stem cells

Eukaryotic gene expression relies on extensive crosstalk between transcription and RNA processing. Changes in this composite regulation network may provide an important means for shaping cell type-specific transcriptomes. Here we show that the RNA-associated protein Srrt/Ars2 sustains embryonic stem cell (ESC) identity by preventing premature termination of numerous transcripts at cryptic cleavage/polyadenylation sites in first introns. Srrt interacts with the nuclear cap-binding complex and facilitates recruitment of the spliceosome component U1 snRNP to cognate intronic positions. At least in some cases, U1 recruited in this manner inhibits downstream cleavage/polyadenylation events through a splicing-independent mechanism called telescripting. We further provide evidence that the naturally high expression of Srrt in ESCs offsets deleterious effects of retrotransposable sequences accumulating in its targets. Our work identifies Srrt as a molecular guardian of the pluripotent cell state.

Besides its role in splicing, U1 snRNP can suppress pre-mRNA cleavage and polyadenylation. The authors show that the nuclear cap-binding complex component Srrt/Ars2 maintains embryonic stem cell identity by promoting U1 recruitment to first introns and preventing premature termination of multiple transcripts.

Latency and interval therapy affect the evolution in metastatic colorectal cancer

While comparison of primary tumor and metastases has highlighted genomic heterogeneity in colorectal cancer (CRC), previous studies have focused on a single metastatic site or limited genomic testing. Combining data from whole exome and ultra-deep targeted sequencing, we explored possible evolutionary trajectories beyond the status of these mutations, particularly among patient-matched metastatic tumors. Our findings confirm the persistence of known clinically-relevant mutations (e.g., those of RAS family of oncogenes) in CRC primary and metastases, yet reveal that latency and interval systemic therapy affect the course of evolutionary events within metastatic lesions. Specifically, our analysis of patient-matched primary and multiple metastatic lesions, developed over time, showed a similar genetic composition for liver metastatic tumors, which were 21-months apart. This genetic makeup was different from those identified in lung metastases developed before manifestation of the second liver metastasis. These results underscore the role of latency in the evolutionary path of metastatic CRC and may have implications for future treatment options.

Clinical relevance of somatic mutations in main driver genes detected in gastric cancer patients by next-generation DNA sequencing

Somatic mutation profiling in gastric cancer (GC) enables main driver mutations to be identified and their clinical and prognostic value to be evaluated. We investigated 77 tumour samples of GC by next-generation sequencing (NGS) with the Ion AmpliSeq Hotspot Panel v2 and a custom panel covering six hereditary gastric cancer predisposition genes (BMPR1A, SMAD4, CDH1, TP53, STK11 and PTEN). Overall, 47 somatic mutations in 14 genes were detected; 22 of these mutations were novel. Mutations were detected most frequently in the CDH1 (13/47) and TP53 (12/47) genes. As expected, somatic CDH1 mutations were positively correlated with distant metastases (p = 0.019) and tumours with signet ring cells (p = 0.043). These findings confirm the association of the CDH1 mutations with diffuse GC type. TP53 mutations were found to be significantly associated with a decrease in overall survival in patients with Lauren diffuse-type tumours (p = 0.0085), T3-T4 tumours (p = 0.037), and stage III-IV tumours (p = 0.013). Our results confirm that the detection of mutations in the main driver genes may have a significant prognostic value for GC patients and provide an independent GC-related set of clinical and molecular genetic data.

Heterozygous rare genetic variants in non-syndromic early-onset obesity

BACKGROUND

Obesity is a very heterogeneous disorder at both the clinical and molecular levels and with high heritability. Several monogenic forms and genes with strong effects have been identified for non-syndromic severe obesity. Novel therapeutic interventions are in development for some genetic forms, emphasizing the importance of determining genetic contributions.

OBJECTIVE

We aimed to define the contribution of rare single-nucleotide genetic variants (RSVs) in candidate genes to non-syndromic severe early-onset obesity (EOO; body mass index (BMI) >+3 standard deviation score, <3 years).

METHODS

Using a pooled DNA-sequencing approach, we screened for RSVs in 15 obesity candidate genes in a series of 463 EOO patients and 480 controls. We also analysed exome data from 293 EOO patients from the "Viva la Familia" (VLF) study as a replication dataset.

RESULTS

Likely or known pathogenic RSVs were identified in 23 patients (5.0%), with 7 of the 15 genes (BDNF, FTO, MC3R, MC4R, NEGR1, PPARG and SIM1) harbouring RSVs only in cases (3.67%) and none in controls. All were heterozygous changes, either de novo (one in BDNF) or inherited from obese parents (seven maternal, three paternal), and no individual carried more than one variant. Results were replicated in the VLF study, where 4.10% of probands carried RSVs in the overrepresented genes. RSVs in five genes were either absent (LEP) or more common in controls than in cases (ADRB3, LEPR, PCSK1 and PCSK2) in both obese datasets.

CONCLUSIONS

Heterozygous RSVs in several candidate genes of the melanocortin pathway are found in ~5.0% patients with EOO. These results support the clinical utility of genetic testing to identify patients who might benefit from targeted therapeutic intervention.

Guanidine N-methylation by BlsL Is Dependent on Acylation of Beta-amine Arginine in the Biosynthesis of Blasticidin S

The peptidyl nucleoside blasticidin S (BS) produced by Streptomyces griseochromogenes was the first non-mercurial fungicide used to prevent rice blast and increasingly used as a selection reagent in transgenic study. Acylation by addition of a leucine residue at the beta amine group of arginine side chain of demethylblasticidin S (DBS) has been proposed as a novel self-resistance to the cytotoxic biosynthetic intermediate. But the resultant product leucyldemethylblasticidin S (LDBS) has not been isolated as a metabolite, and LDBS synthetase activity remained to be demonstrated in S. griseochromogenes. In this study, we isolated LDBS in a BS heterologous producer S. lividans WJ2 upon the deletion of blsL, which encodes a S-Adenosyl methionine-dependent methyltransferase. Purified BlsL efficiently methylated LDBS at the delta N of beta-arginine to generate the ultimate intermediate LBS, but nearly didn’t methylate DBS to final product BS. Above experiments demonstrated that LDBS is indeed an intermediate in BS biosynthetic pathway, and acylation of beta-amino group of arginine side chain is prerequisite for efficient guanidine N-methylation in addition to being a self-resistance mechanism.