To gain insight into the transcription programs activated during the formation of Drosophila larval structures, we carried out single cell RNA sequencing during two periods of Drosophila embryogenesis: stages 10-12, when most organs are first specified and initiate morphological and physiological specialization; and stages 13-16, when organs achieve their final mature architectures and begin to function. Our data confirm previous findings with regards to functional specialization of some organs - the salivary gland and trachea - and clarify the embryonic functions of another - the plasmatocytes. We also identify two early developmental trajectories in germ cells and uncover a potential role for proteolysis during germline stem cell specialization. We identify the likely cell type of origin for key components of the Drosophila matrisome and several commonly used Drosophila embryonic cell culture lines. Finally, we compare our findings with other recent related studies and with other modalities for identifying tissue-specific gene expression patterns. These data provide a useful community resource for identifying many new players in tissue-specific morphogenesis and functional specialization of developing organs.
Non-invasive methods of detecting radiation exposure show promise to improve upon current approaches to biological dosimetry in ease, speed, and accuracy. Here we developed a pipeline that employs Fourier transform infrared (FTIR) spectroscopy in the mid-infrared spectrum to identify a signature of low dose ionizing radiation exposure in mouse ear pinnae over time. Mice exposed to 0.1 to 2 Gy total body irradiation were repeatedly measured by FTIR at the stratum corneum of the ear pinnae. We found significant discriminative power for all doses and time-points out to 90 days after exposure. Classification accuracy was maximized when testing 14 days after exposure (specificity > 0.9 with a sensitivity threshold of 0.9) and dropped by roughly 30% sensitivity at 90 days. Infrared frequencies point towards biological changes in DNA conformation, lipid oxidation and accumulation and shifts in protein secondary structure. Since only hundreds of samples were used to learn the highly discriminative signature, developing human-relevant diagnostic capabilities is likely feasible and this non-invasive procedure points toward rapid, non-invasive, and reagent-free biodosimetry applications at population scales.
PURPOSE Tissue-agnostic biomarkers that capture the commonality in cancer biology, may provide a new avenue for treatment development and optimization across cancer types. Here, we aimed to evaluate and validate the clinical value of a tissue-agnostic cellular morphometrics biomarker (CMB) signature, which was discovered by artificial intelligence (AI) from H&E-stained whole-slide images (WSI) of diagnostic slides of colon cancers, in pan-gastrointestinal (pan-GI) pre-cancer lesions and cancers. METHODS We discovered CMBs from WSI using our well-established CMB-ML pipeline and established a CMB risk score (CMBRS) using multivariate regression models. Based on CMBRS, we assigned individual patients from The Cancer Genome Atlas Colon Adenocarcinoma Cohort (TCGA-COAD) (n=430) to CMB risk groups (CMBRG). We then extensively evaluated tissue-agnostic clinical value of CMB signature, CMBRS and CMBRG in multi-cohorts with different types of GI cancer (n=2,219) and risk assessment of precancerous lesions (n=1,016). We unraveled each CMB-related biological function using bulk RNA-sequencing, single-cell RNA-sequencing (scRNA-seq) and opal multiplex immunohistochemistry (IHC) techniques. RESULTS From the TCGA-COAD cohort, we developed a 13-CMB signature and constructed CMBRS/CMBRG that predict prognosis of colon cancer patients. Importantly, this 13-CMB signature proved prognostic and predictive values for TCGA patients with rectal, gastric and esophageal cancer independent of traditional clinical factors. These findings were independently validated using multiple cohorts from Drum Tower Hospital. Moreover, 13-CMB signature exhibited the power for risk stratification of colon adenoma and early esophageal neoplastic lesion patients for predicting cancer progression. In addition, we demonstrated and validated independent prognostic impacts of gene signatures and CMB signatures and a significant increase in predictive power by integration of CMB signature, gene signature and clinical factors. Correlations between CMBs and gene expression levels revealed the association of each CMB with biological functions including cell proliferation, epithelial-to-mesenchymal transition and immune microenvironment. The association of CMBs with the immune microenvironment was prospectively validated by scRNA-seq and was further confirmed by Opal multiplex IHC staining in colon cancer. CONCLUSION This study demonstrates the clinical value of tissue-agnostic AI-empowered CMB signature from WSI with defined biological functions, which can be used in clinical settings to assess risk, diagnose disease, and guide clinical interventions. Tissue-agnostic CMBs potentially provide a new avenue for a rapid, robust and cost-effective cross-cancer prediction that is essential for developing common treatment strategy for multiple cancers. ### Competing Interest Statement The authors have declared no competing interest. ### Funding Statement This work has been supported by the National Natural Science Foundation of China (ID Number: 82272952), the Natural Science Foundation of Jiangsu Province for Excellent Young Scholars (ID Number: BK20220094), China Postdoctoral Science Foundation (ID Number: 2022M721579) and funds for Clinical Trials from the Affiliated Drum Tower Hospital, Medical School of Nanjing University (ID Number: 2021-LCYJ-PY-21). ### Author Declarations I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained. Yes The details of the IRB/oversight body that provided approval or exemption for the research described are given below: The hospital validation study was approved by the Institutional Review Board (IRB) at the participating hospital and was independently carried out at Nanjing Drum Tower Hospital. I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals. Yes I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance). Yes I have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable. Yes Whole slide images and clinical data of the TCGA cohorts were downloaded from the TCGA GDC portal (). Processed data related to TCGA and Drum Tower cohorts have been provided with the manuscript. Raw data from the Nanjing Drum Tower Hospital is not currently permitted in public repositories because ethical and legal implications are still being discussed at an institutional level. * H&E : hematoxylin and eosin AI : artificial intelligence GI : gastrointestinal CMB : cellular morphometrics biomarker WSI : whole-slide images CMBRS : CMB risk score CMBRG : CMB risk group OS : overall survival scRNA-seq : single cell RNA sequencing IHC : immunohistochemistry TCGA-COAD : The Cancer Genome Atlas - Colon Adenocarcinoma TCGA-STAD : The Cancer Genome Atlas - Stomach Adenocarcinom TCGA-READ : The Cancer Genome Atlas - Rectum Adenocarcinoma TCGA-ESCA : The Cancer Genome Atlas - Esophageal Carcinoma LGIN : Low-Grade Intraepithelial Neoplasia HGIN : High-Grade Intraepithelial Neoplasia CAP : Colon Adenomatous Polyps EEL : Early Esophageal Lesion CRC : Colorectal Cancer DT : Drum Tower Hospital
Synthetic biology is creating genetically engineered organisms at an increasing rate for many potentially valuable applications, but this potential comes with the risk of misuse or accidental release. To begin to address this issue, we have developed a system called GUARDIAN that can automatically detect signatures of engineering in DNA sequencing data, and we have conducted a blinded test of this system using a curated Test and Evaluation (T&E) data set. GUARDIAN uses an ensemble approach based on the guiding principle that no single approach is likely to be able to detect engineering with perfect accuracy. Critically, ensembling enables GUARDIAN to detect sequence inserts in 13 target organisms with a high degree of specificity that requires no subject matter expert (SME) review.
SUMMARY A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here we present the culmination of the modERN (model organism Encyclopedia of Regulatory Networks) consortium that systematically assayed TF binding events in vivo in two major model organisms, Drosophila melanogaster (fly) and Caenorhabditis elegans (worm). We describe key features of these datasets, comprising 604 TFs identifying 3.6M sites in the fly and 350 TFs identifying 0.9 M sites in the worm. Applying a machine learning model to these data identifies sets of TFs with a prominent role in promoting target gene expression in specific cell types. TF binding data are available through the ENCODE Data Coordinating Center and at https://epic.gs.washington.edu/modERNresource , which provides access to processed and summary data, as well as widgets to probe cell type-specific TF-target relationships. These data are a rich resource that should fuel investigations into TF function during development.
A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here, we present the culmination of the efforts of the modENCODE (model organism Encyclopedia of DNA Elements) and modERN (model organism Encyclopedia of Regulatory Networks) consortia to systematically assay TF binding events in vivo in two major model organisms, Drosophila melanogaster (fly) and Caenorhabditis elegans (worm). These data sets comprise 605 TFs identifying 3.6 M sites in the fly and 356 TFs identifying 0.9 M sites in the worm, and represent the majority of the regulatory space in each genome. We demonstrate that TFs associate with chromatin in clusters termed "metapeaks," that larger metapeaks have characteristics of high-occupancy target (HOT) regions, and that the importance of consensus sequence motifs bound by TFs depends on metapeak size and complexity. Combining ChIP-seq data with single-cell RNA-seq data in a machine-learning model identifies TFs with a prominent role in promoting target gene expression in specific cell types, even differentiating between parent-daughter cells during embryogenesis. These data are a rich resource for the community that should fuel and guide future investigations into TF function. To facilitate data accessibility and utility, all strains expressing green fluorescent protein (GFP)-tagged TFs are available at the stock centers for each organism. The chromatin immunoprecipitation sequencing data are available through the ENCODE Data Coordinating Center, GEO, and through a direct interface that provides rapid access to processed data sets and summary analyses, as well as widgets to probe the cell-type-specific TF-target relationships.
Ionizing radiation induces complex changes in cells and tissues. The conventional approach to biological dosimetry has been to integrate physical and clinical measurements to optimize dose assessment. Molecular biodosimetry is an effective strategy to monitor radiation exposure and hematologic, cytogenetic, protein and transcript-based approaches have been developed to increase dose estimation accuracy. However, these approaches are invasive, time-consuming, have limited effectiveness over time, and importantly do not accurately inform on low dose radiation exposures. Therefore, novel, non-invasive, biomarkers are required that can overcome these limitations. We developed a pipeline that employs Fourier transform infrared (FTIR) spectroscopy in the mid-infrared spectrum to identify a signature of low dose ionizing radiation exposure in mouse ear pinnae over time. Two cohorts of C57BL/6J and one cohort of BALB/c mice were followed for ninety days after total body X-ray exposures of 10, 50, 100 or 200 cGy. The stratum corneum layer of the ear pinnae were repeatedly measured with attenuated total reflectance - focal plane array (ATR-FPA)-FTIR at 5, 14, 21, 49 and 90 days after radiation exposure. We found statistically significant discriminative power for all doses and all tested time-points out to 90 days after exposure. Classification accuracy was maximized when testing 14 days after exposure (specificity > 0.9 with a sensitivity threshold of 0.9) and dropped by roughly 30% sensitivity at 90 days. Since only hundreds of samples were required to learn highly discriminative signatures, developing human-relevant diagnostic capabilities is likely feasible and this non-invasive procedure points toward future non-invasive biodosimetry applications at population scales.
We describe a next-generation Drosophila protein interaction map-"DPIM2"-established from affinity purification-mass spectrometry of 5,805 baits, covering the largest fraction of the Drosophila proteome. The network contains 32,668 interactions among 3,644 proteins, organized into 632 clusters representing putative functional modules. Our analysis expands the pool of known protein interactions in Drosophila, provides annotation for poorly studied genes, and postulates previously undescribed protein interaction relationships. The predictive power and functional relevance of this network are probed through the lens of the Notch signaling pathway, and we find that newly identified members of complexes that include known Notch modifiers can also modulate Notch signaling. DPIM2 allows direct comparisons with a recently published human protein interaction network, defining the existence of functional interactions conserved across species. Thus, DPIM2 defines a valuable resource for predicting protein co-complex memberships and functional associations as well as generates functional hypotheses regarding specific protein interactions.
The insights into interactions between host genetics and gut microbiome (GM) in colorectal tumor susceptibility (CTS) remains lacking. We used Collaborative Cross mouse population model to identify genetic and microbial determinants of Azoxymethane-induced CTS. We identified 4417 CTS-associated single nucleotide polymorphisms (SNPs) containing 334 genes that were transcriptionally altered in human colorectal cancers (CRCs) and consistently clustered independent human CRC cohorts into two subgroups with different prognosis. We discovered a set of genera in early-life associated with CTS and defined a 16-genus signature that accurately predicted CTS, the majority of which were correlated with human CRCs. We identified 547 SNPs associated with abundances of these genera. Mediation analysis revealed GM as mediators partially exerting the effect of SNP UNC3869242 within Duox2 on CTS. Intestine cell-specific depletion of Duox2 altered GM composition and contribution of Duox2 depletion to CTS was significantly influenced by GM. Our findings provide potential novel targets for personalized CRC prevention and treatment.
A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here, we present the culmination of the efforts of the modENCODE (model organism Encyclopedia of DNA Elements) and modERN (model organism Encyclopedia of Regulatory Networks) consortia to systematically assay TF binding events in vivo in two major model organisms,Drosophila melanogaster(fly) andCaenorhabditis elegans(worm). These data sets comprise 605 TFs identifying 3.6 M sites in the fly and 356 TFs identifying 0.9 M sites in the worm, and represent the majority of the regulatory space in each genome. We demonstrate that TFs associate with chromatin in clusters termed “metapeaks,” that larger metapeaks have characteristics of high-occupancy target (HOT) regions, and that the importance of consensus sequence motifs bound by TFs depends on metapeak size and complexity. Combining ChIP-seq data with single-cell RNA-seq data in a machine-learning model identifies TFs with a prominent role in promoting target gene expression in specific cell types, even differentiating between parent–daughter cells during embryogenesis. These data are a rich resource for the community that should fuel and guide future investigations into TF function. To facilitate data accessibility and utility, all strains expressing green fluorescent protein (GFP)-tagged TFs are available at the stock centers for each organism. The chromatin immunoprecipitation sequencing data are available through the ENCODE Data Coordinating Center, GEO, and through a direct interface that provides rapid access to processed data sets and summary analyses, as well as widgets to probe the cell-type-specific TF–target relationships.
Abstract Transcription factors (TFs) play a key role in development and in cellular responses to the environment by activating or repressing the transcription of target genes in precise spatial and temporal patterns. In order to develop a catalog of target genes of Drosophila melanogaster TFs, the modERN consortium systematically knocked down the expression of TFs using RNAi in whole embryos followed by RNA-seq. We generated data for 45 TFs which have 18 different DNA-binding domains and are expressed in 15 of the 16 organ systems. The range of inactivation of the targeted TFs by RNAi ranged from log2fold change −3.52 to +0.49. The TFs also showed remarkable heterogeneity in the numbers of candidate target genes identified, with some generating thousands of candidates and others only tens. We present detailed analysis from five experiments, including those for three TFs that have been the focus of previous functional studies (ERR, sens, and zfh2) and two previously uncharacterized TFs (sens-2 and CG32006), as well as short vignettes for selected additional experiments to illustrate the utility of this resource. The RNA-seq datasets are available through the ENCODE DCC (http://encodeproject.org) and the Sequence Read Archive (SRA). TF and target gene expression patterns can be found here: https://insitu.fruitfly.org. These studies provide data that facilitate scientific inquiries into the functions of individual TFs in key developmental, metabolic, defensive, and homeostatic regulatory pathways, as well as provide a broader perspective on how individual TFs work together in local networks during embryogenesis.
Short polypeptides encoded by small open reading frames (smORFs) are ubiquitously found in eukaryotic genomes and are important regulators of physiology, development, and mitochondrial processes. Here, we focus on a subset of 298 smORFs that are evolutionarily conserved between Drosophila melanogaster and humans. Many of these smORFs are conserved broadly in the bilaterian lineage, and ∼182 are conserved in plants. We observe remarkably heterogeneous spatial and temporal expression patterns of smORF transcripts-indicating wide-spread tissue-specific and stage-specific mitochondrial architectures. In addition, an analysis of annotated functional domains reveals a predicted enrichment of smORF polypeptides localizing to mitochondria. We conduct an embryonic ribosome profiling experiment and find support for translation of 137 of these smORFs during embryogenesis. We further embark on functional characterization using CRISPR knockout/activation, RNAi knockdown, and cDNA overexpression, revealing diverse phenotypes. This study underscores the importance of identifying smORF function in disease and phenotypic diversity.
Generating reference maps of interactome networks illuminates genetic studies by providing a protein-centric approach to finding new components of existing pathways, complexes, and processes. We apply state-of-the-art methods to identify binary protein-protein interactions (PPIs) for Drosophila melanogaster . Four all-by-all yeast two-hybrid (Y2H) screens of > 10,000 Drosophila proteins result in the ‘FlyBi’ dataset of 8723 PPIs among 2939 proteins. Testing subsets of data from FlyBi and previous PPI studies using an orthogonal assay allows for normalization of data quality; subsequent integration of FlyBi and previous data results in an expanded binary Drosophila reference interaction network, DroRI, comprising 17,232 interactions among 6511 proteins. We use FlyBi data to generate an autophagy network, then validate in vivo using autophagy-related assays. The deformed wings ( dwg ) gene encodes a protein that is both a regulator and a target of autophagy. Altogether, these resources provide a foundation for building new hypotheses regarding protein networks and function.
Hematologic toxicity is a common side effect of multimodal cancer therapy. Nearly all animal studies investigating the causes of radiotherapy-induced hematologic toxicity use inbred strains with limited genetic diversity and do not reflect the diverse responses observed in humans. We used the population-based Collaborative Cross (CC) mouse resource to investigate the genetic architecture of the acute and persistent immune response after radiation exposure by measuring 22 immune parameters in 1,720 CC mice representing 35 strains. We determined relative acute and persistent radiation resistance scores at the individual strain level considering contributions from all immune parameters. Genome-wide association analysis identified quantitative trait loci associated with baseline and radiation responses. A cross-species radiation resistance score predicted recurrence-free survival in medulloblastoma patients. We present a community resource of immune parameters and genome-wide association analyses before and after radiation exposure for future investigations of the contributions of host genetics on radiosensitivity.
Supplementary Tables S1-S4 from Amplification of PVT1 Contributes to the Pathophysiology of Ovarian and Breast Cancer
Development of effective therapies against SARS-CoV-2 infections relies on mechanistic knowledge of virus host interface. Abundant physical interactions between viral and host proteins have been identified, but few have been functionally characterized. Harnessing the power of fly genetics, we develop a comprehensive Drosophila COVID-19 resource (DCR) consisting of publicly available strains for conditional tissue-specific expression of all SARS-CoV-2 encoded proteins, UAS-human cDNA transgenic lines encoding established host-viral interacting factors, and GAL4 insertion lines disrupting fly homologs of SARS-CoV-2 human interacting proteins. We demonstrate the utility of the DCR to functionally assess SARS-CoV-2 genes and candidate human binding partners. We show that NSP8 engages in strong genetic interactions with several human candidates, most prominently with the ATE1 arginyltransferase to induce actin arginylation and cytoskeletal disorganization, and that two ATE1 inhibitors can reverse NSP8 phenotypes. The DCR enables parallel global-scale functional analysis of SARS-CoV-2 components in a prime genetic model system.
Naturally produced peptides (<100 amino acids) are important regulators of physiology, development, and metabolism. Recent studies have predicted that thousands of peptides may be translated from transcripts containing small open-reading frames (smORFs). Here, we describe two peptides in Drosophila encoded by conserved smORFs, Sloth1 and Sloth2. These peptides are translated from the same bicistronic transcript and share sequence similarities, suggesting that they encode paralogs. Yet, Sloth1 and Sloth2 are not functionally redundant, and loss of either peptide causes animal lethality, reduced neuronal function, impaired mitochondrial function, and neurodegeneration. We provide evidence that Sloth1/2 are highly expressed in neurons, imported to mitochondria, and regulate mitochondrial complex III assembly. These results suggest that phenotypic analysis of smORF genes in Drosophila can provide a wealth of information on the biological functions of this poorly characterized class of genes.
Radiotherapy is the current standard of care for more than 50% of all cancer patients. Improvements in radiotherapy (RT) technology have increased tumor targeting and normal tissue sparing. Radiations at ultra-high dose rates required for FLASH-RT effects have sparked interest in potentially providing additional differential therapeutic benefits. We present a new experimental platform that is the first one to deliver petawatt laser-driven proton pulses of 2 MeV energy at 0.2 Hz repetition rate by means of a compact, tunable active plasma lens beamline to biological samples. Cell monolayers grown over a 10 mm diameter field were exposed to clinically relevant proton doses ranging from 7 to 35 Gy at ultra-high instantaneous dose rates of 107 Gy/s. Dose-dependent cell survival measurements of human normal and tumor cells exposed to LD protons showed significantly higher cell survival of normal-cells compared to tumor-cells for total doses of 7 Gy and higher, which was not observed to the same extent for X-ray reference irradiations at clinical dose rates. These findings provide preliminary evidence that compact LD proton sources enable a new and promising platform for investigating the physical, chemical and biological mechanisms underlying the FLASH effect.
SARS-CoV-2 spread rapidly, causing millions of deaths across the globe. As a result, demand for medical supplies and personal protective equipment (PPE) surged and supplies dwindled. Separate entirely, hospital-acquired infections have become commonplace and challenging to treat. To explore the potential of novel sterilization techniques, this study evaluated the disinfection efficacy of Fathhome's ozone-based, dry-sanitizing device by dose and time response. Inactivation of human pathogens was tested on non-porous (plastic) surfaces. 95.42-100% inactivation was observed across all types of vegetative microorganisms and 27.36% inactivation of bacterial endospores tested, with no residual ozone detectable after completion. These results strongly support the hypothesis that Fathhome's commercial implementation of gas-based disinfection is suitable for rapid decontamination of a wide variety of pathogens on PPE and other industrially relevant materials.
Azoxymethane (AOM) is a widely used carcinogen to study chemical-induced colorectal carcinogenesis and is an agent for studying fulminant hepatic failure. The inter-strain susceptibility to acute toxicity by AOM has been reported, but its association with host genetics or gut microbiota remains largely unexplored. Here a cohort of genetically diverse Collaborative Cross (CC) mice was used to assess the contribution of host genetics and the gut microbiome to AOM-induced acute toxicity. We observed variation in AOM-induced acute liver failure across CC strains. Quantitative trait loci (QTL) analysis revealed three chromosome regions significantly associated with AOM toxicity. Genes located within these QTL, including peroxisome proliferator-activated receptor alpha (Ppara), were enriched for enzyme activator and nucleoside-triphosphatase regulator activity. We further demonstrated that the protein level of PPARα in liver tissues from sensitive strains was remarkably lower compared to levels in resistant strains, consistent with protective role of PPAR family in liver injury. We discovered that the abundance levels of gut microbial families Anaeroplasmataceae, Ruminococcaceae, Lactobacillaceae, Akkermansiaceae and Clostridiaceae were significantly higher in the sensitive strains compared to the resistant strains. Using a random forest classifier method, we determined that the relative abundance levels of these microbial families predicted AOM toxicity with the area under the receiver-operating curve (AUC) of 0.75. Combining the three genetic loci and five microbial families increased the predictive accuracy of AOM toxicity (AUC of 0.99). Moreover, we found that Ruminococcaceae and Lactobacillaceae acted as mediators between host genetics and AOM toxicity. In conclusion, this study shows that host genetics and specific microbiome members play a critical role in AOM-induced acute toxicity, which provides a framework for analysis of the health effects from environmental toxicants.