Search

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.

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.

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.

Microbacterium sp. BDGP8 is a species of facultative anaerobic gram-positive bacterium of the family Microbacteriaceae. The complete genome consists of a single circular chromosome of 3,293,567 bp with a G + C content of 69.84% and two plasmids of 49,365 bp and 32,884 bp.

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.

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.

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.

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.

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.

The gut microbiome composition is influenced by many factors including environmental exposures. Here, we investigated the effect of thirdhand cigarette smoke (THS) and exposure age on gut microbiome diversity. C57BL/6 mice were exposed to THS at human exposure relevant levels for three weeks during three different life stages: postnatal (0-3 weeks of age), pubescent (4-7 weeks of age), and adult (9-12 weeks of age), respectively. Cecal microbiome profiles were assessed at 17 weeks of age by 16S rRNA gene sequencing. We found that age at THS exposure strongly influenced the cecal microbiome composition. Although postnatal THS exposure significantly influenced the microbial composition, pubescent and adulthood exposures only had minor effects. The microbiome of postnatally THS-exposed mice significantly increased several degradation pathways that regulate glycolysis and pyruvate decarboxylation, and significantly decreased coenzyme A biosynthesis and pyrimidine deoxyribonucleoside salvage. Our results indicate that mouse postnatal development is particularly susceptible to persistent THS exposure effects on the gut microbiome.

As a means to understand human neuropsychiatric disorders from human brain samples, we compared the transcription patterns and histological features of postmortem brain to fresh human neocortex isolated immediately following surgical removal. Compared to a number of neuropsychiatric disease-associated postmortem transcriptomes, the fresh human brain transcriptome had an entirely unique transcriptional pattern. To understand this difference, we measured genome-wide transcription as a function of time after fresh tissue removal to mimic the postmortem interval. Within a few hours, a selective reduction in the number of neuronal activity-dependent transcripts occurred with relative preservation of housekeeping genes commonly used as a reference for RNA normalization. Gene clustering indicated a rapid reduction in neuronal gene expression with a reciprocal time-dependent increase in astroglial and microglial gene expression that continued to increase for at least 24 h after tissue resection. Predicted transcriptional changes were confirmed histologically on the same tissue demonstrating that while neurons were degenerating, glial cells underwent an outgrowth of their processes. The rapid loss of neuronal genes and reciprocal expression of glial genes highlights highly dynamic transcriptional and cellular changes that occur during the postmortem interval. Understanding these time-dependent changes in gene expression in post mortem brain samples is critical for the interpretation of research studies on human brain disorders.

The gut microbiome produces vitamins, nutrients, and neurotransmitters, and helps to modulate the host immune system-and also plays a major role in the metabolism of many exogenous compounds, including drugs and chemical toxicants. However, the extent to which specific microbial species or communities modulate hazard upon exposure to chemicals remains largely opaque. Focusing on the effects of collateral dietary exposure to the widely used herbicide atrazine, we applied integrated omics and phenotypic screening to assess the role of the gut microbiome in modulating host resilience in Drosophila melanogaster. Transcriptional and metabolic responses to these compounds are sex-specific and depend strongly on the presence of the commensal microbiome. Sequencing the genomes of all abundant microbes in the fly gut revealed an enzymatic pathway responsible for atrazine detoxification unique to Acetobacter tropicalis. We find that Acetobacter tropicalis alone, in gnotobiotic animals, is sufficient to rescue increased atrazine toxicity to wild-type, conventionally reared levels. This work points toward the derivation of biotic strategies to improve host resilience to environmental chemical exposures, and illustrates the power of integrative omics to identify pathways responsible for adverse health outcomes.

BackgroundRecent evidence has linked the gut microbiome to host behavior via the gut-brain axis [1-3]; however, the underlying mechanisms remain unexplored. Here, we determined the links between host genetics, the gut microbiome and memory using the genetically defined Collaborative Cross (CC) mouse cohort, complemented with microbiome and metabolomic analyses in conventional and germ-free (GF) mice.ResultsA genome-wide association analysis (GWAS) identified 715 of 76,080 single-nucleotide polymorphisms (SNPs) that were significantly associated with short-term memory using the passive avoidance model. The identified SNPs were enriched in genes known to be involved in learning and memory functions. By 16S rRNA gene sequencing of the gut microbial community in the same CC cohort, we identified specific microorganisms that were significantly correlated with longer latencies in our retention test, including a positive correlation with Lactobacillus. Inoculation of GF mice with individual species of Lactobacillus (L. reuteri F275, L. plantarum BDGP2 or L. brevis BDGP6) resulted in significantly improved memory compared to uninoculated or E. coli DH10B inoculated controls. Untargeted metabolomics analysis revealed significantly higher levels of several metabolites, including lactate, in the stools of Lactobacillus-colonized mice, when compared to GF control mice. Moreover, we demonstrate that dietary lactate treatment alone boosted memory in conventional mice. Mechanistically, we show that both inoculation with Lactobacillus or lactate treatment significantly increased the levels of the neurotransmitter, gamma-aminobutyric acid (GABA), in the hippocampus of the mice.ConclusionTogether, this study provides new evidence for a link between Lactobacillus and memory and our results open possible new avenues for treating memory impairment disorders using specific gut microbial inoculants and/or metabolites. Video Abstract.

Lactobacillus brevis Oregon-R-modENCODE strain BDGP6 was isolated from the gut of Drosophila melanogaster for functional host-microbial interaction studies. The bacterial chromosome is a single circular DNA molecule of 2,785,111 bp with a G+C content of 46%.

Citrobacter freundii is a species of facultative anaerobic Gram-negative bacteria of the family Enterobacteriaceae The complete genome is composed of a single chromosomal circle of 4,957,773 bp with a G+C content of 52%.

ObjectiveThe Collaborative Cross (CC) is a mouse population model with diverse and reproducible genetic backgrounds used to identify novel disease models and genes that contribute to human disease. Since spontaneous tumour susceptibility in CC mice remains unexplored, we assessed tumour incidence and spectrum.DesignWe monitored 293 mice from 18 CC strains for tumour development. Genetic association analysis and RNA sequencing were used to identify susceptibility loci and candidate genes. We analysed genomes of patients with gastric cancer to evaluate the relevance of genes identified in the CC mouse model and measured the expression levels of ISG15 by immunohistochemical staining using a gastric adenocarcinoma tissue microarray. Association of gene expression with overall survival (OS) was assessed by Kaplan-Meier analysis.ResultsCC mice displayed a wide range in the incidence and types of spontaneous tumours. More than 40% of CC036 mice developed gastric tumours within 1 year. Genetic association analysis identified Nfκb1 as a candidate susceptibility gene, while RNA sequencing analysis of non-tumour gastric tissues from CC036 mice showed significantly higher expression of inflammatory response genes. In human gastric cancers, the majority of human orthologues of the 166 mouse genes were preferentially altered by amplification or deletion and were significantly associated with OS. Higher expression of the CC036 inflammatory response gene signature is associated with poor OS. Finally, ISG15 protein is elevated in gastric adenocarcinomas and correlated with shortened patient OS.ConclusionsCC strains exhibit tremendous variation in tumour susceptibility, and we present CC036 as a spontaneous laboratory mouse model for studying human gastric tumourigenesis.

Agrobacterium sp. strain 33MFTa1.1 was isolated for functional host-microbe interaction studies from the Thlaspi arvense root-associated microbiome. The complete genome is comprised of a circular chromosome of 2,771,937 bp, a linear chromosome of 2,068,443 bp, and a plasmid of 496,948 bp, with G+C contents of 59%, 59%, and 58%, respectively.

Variations in oral bacterial communities have been linked to oral cancer suggesting that the oral microbiome is an etiological factor that can influence oral cancer development. The 4-nitroquinoline 1-oxide (4-NQO)-induced murine oral and esophageal cancer model is frequently used to assess the effects of preventive and/or therapeutic agents. We used this model to assess the impact of the microbiome on tumorigenesis using axenic (germ-free) and conventionally housed mice. Increased toxicity was observed in germ-free mice, however, no difference in tumor incidence, multiplicity, and size was observed. Transcriptional profiling of liver tissue from germ-free and conventionally housed mice identified 254 differentially expressed genes including ten cytochrome p450 enzymes, the largest family of phase-1 drug metabolizing enzymes in the liver. Gene ontology revealed that differentially expressed genes were enriched for liver steatosis, inflammation, and oxidative stress in livers of germ-free mice. Our observations emphasize the importance of the microbiome in mediating chemical toxicity at least in part by altering host gene expression. Studies on the role of the microbiome in chemical-induced cancer using germ-free animal models should consider the potential difference in dose due to the microbiome-mediated changes in host metabolizing capacity, which might influence the ability to draw conclusions especially for tumor induction models that are dose dependent.

Loss-of-function mutations in the histone demethylases KDM5A, KDM5B, or KDM5C are found in intellectual disability (ID) and autism spectrum disorders (ASD) patients. Here, we use the model organism Drosophila melanogaster to delineate how KDM5 contributes to ID and ASD. We show that reducing KDM5 causes intestinal barrier dysfunction and changes in social behavior that correlates with compositional changes in the gut microbiota. Therapeutic alteration of the dysbiotic microbiota through antibiotic administration or feeding with a probiotic Lactobacillus strain partially rescues the behavioral, lifespan, and cellular phenotypes observed in kdm5-deficient flies. Mechanistically, KDM5 was found to transcriptionally regulate component genes of the immune deficiency (IMD) signaling pathway and subsequent maintenance of host-commensal bacteria homeostasis in a demethylase-dependent manner. Together, our study uses a genetic approach to dissect the role of KDM5 in the gut-microbiome-brain axis and suggests that modifying the gut microbiome may provide therapeutic benefits for ID and ASD patients.

Identifying functional enhancer elements in metazoan systems is a major challenge. Large-scale validation of enhancers predicted by ENCODE reveal false-positive rates of at least 70%. We used the pregrastrula-patterning network of Drosophila melanogaster to demonstrate that loss in accuracy in held-out data results from heterogeneity of functional signatures in enhancer elements. We show that at least two classes of enhancers are active during early Drosophila embryogenesis and that by focusing on a single, relatively homogeneous class of elements, greater than 98% prediction accuracy can be achieved in a balanced, completely held-out test set. The class of well-predicted elements is composed predominantly of enhancers driving multistage segmentation patterns, which we designate segmentation driving enhancers (SDE). Prediction is driven by the DNA occupancy of early developmental transcription factors, with almost no additional power derived from histone modifications. We further show that improved accuracy is not a property of a particular prediction method: after conditioning on the SDE set, naïve Bayes and logistic regression perform as well as more sophisticated tools. Applying this method to a genome-wide scan, we predict 1,640 SDEs that cover 1.6% of the genome. An analysis of 32 SDEs using whole-mount embryonic imaging of stably integrated reporter constructs chosen throughout our prediction rank-list showed >90% drove expression patterns. We achieved 86.7% precision on a genome-wide scan, with an estimated recall of at least 98%, indicating high accuracy and completeness in annotating this class of functional elements.

During terminal erythropoiesis, the splicing machinery in differentiating erythroblasts executes a robust intron retention (IR) program that impacts expression of hundreds of genes. We studied IR mechanisms in the SF3B1 splicing factor gene, which expresses ∼50% of its transcripts in late erythroblasts as a nuclear isoform that retains intron 4. RNA-seq analysis of nonsense-mediated decay (NMD)-inhibited cells revealed previously undescribed splice junctions, rare or not detected in normal cells, that connect constitutive exons 4 and 5 to highly conserved cryptic cassette exons within the intron. Minigene splicing reporter assays showed that these cassettes promote IR. Genome-wide analysis of splice junction reads demonstrated that cryptic noncoding cassettes are much more common in large (>1 kb) retained introns than they are in small retained introns or in nonretained introns. Functional assays showed that heterologous cassettes can promote retention of intron 4 in the SF3B1 splicing reporter. Although many of these cryptic exons were spliced inefficiently, they exhibited substantial binding of U2AF1 and U2AF2 adjacent to their splice acceptor sites. We propose that these exons function as decoys that engage the intron-terminal splice sites, thereby blocking cross-intron interactions required for excision. Developmental regulation of decoy function underlies a major component of the erythroblast IR program.

High-content image acquisition is generally limited to cells grown in culture, requiring complex hardware and preset imaging modalities. Here we report an open source software package, OpenHiCAMM (Open Hi Content Acquisition for μManager), that provides a flexible framework for integration of generic microscope-associated robotics and image processing with sequential work-flows. As an example, we imaged Drosophila embryos, detecting the embryos at low resolution, followed by re-imaging the detected embryos at high resolution, suitable for computational analysis and screening. The OpenHiCAMM package is easy to use and adapt for automating complex microscope image tasks. It expands our abilities for high-throughput image-based screens to a new range of biological samples, such as organoids, and will provide a foundation for bioimaging systems biology.

To develop a catalog of regulatory sites in two major model organisms, Drosophila melanogaster and Caenorhabditis elegans, the modERN (model organism Encyclopedia of Regulatory Networks) consortium has systematically assayed the binding sites of transcription factors (TFs). Combined with data produced by our predecessor, modENCODE (Model Organism ENCyclopedia Of DNA Elements), we now have data for 262 TFs identifying 1.23 M sites in the fly genome and 217 TFs identifying 0.67 M sites in the worm genome. Because sites from different TFs are often overlapping and tightly clustered, they fall into 91,011 and 59,150 regions in the fly and worm, respectively, and these binding sites span as little as 8.7 and 5.8 Mb in the two organisms. Clusters with large numbers of sites (so-called high occupancy target, or HOT regions) predominantly associate with broadly expressed genes, whereas clusters containing sites from just a few factors are associated with genes expressed in tissue-specific patterns. All of the strains expressing GFP-tagged TFs are available at the stock centers, and the chromatin immunoprecipitation sequencing data are available through the ENCODE Data Coordinating Center and also through a simple interface (http://epic.gs.washington.edu/modERN/) that facilitates rapid accessibility of processed data sets. These data will facilitate a vast number of scientific inquiries into the function of individual TFs in key developmental, metabolic, and defense and homeostatic regulatory pathways, as well as provide a broader perspective on how individual TFs work together in local networks and globally across the life spans of these two key model organisms.

Acetobacter pomorum Oregon-R-modENCODE strain BDGP5 was isolated from Drosophila melanogaster for functional host-microbe interaction studies. The complete genome is composed of a single chromosomal circle of 2,848,089 bp, with a G+C content of 53% and three plasmids of 131,455 bp, 19,216 bp, and 9,160 bp.

Acetobacter tropicalis Oregon-R-modENCODE strain BDGP1 was isolated from Drosophila melanogaster for functional host-microbe interaction studies. The complete genome comprises a single chromosomal circle of 3,988,649 bp with a G+C content of 56% and a conjugative plasmid of 151,013 bp.

Lactobacillus plantarum Oregon-R-modENCODE strain BDGP2 was isolated from the gut of Drosophila melanogaster for functional host microbial interaction studies. The complete genome comprised a single circular genome of 3,407,160 bp, with a G+C content of 44%, and four plasmids.

Bacillus kochii Oregon-R-modENCODE strain BDGP4 was isolated from the gut of Drosophila melanogaster for functional host microbial interaction studies. The complete genome comprised a single chromosomal circle of 4,557,232 bp with a G+C content of 37% and a single plasmid of 137,143 bp.

Enterococcus durans Oregon-R-modENCODE strain BDGP3 was isolated from the Drosophila melanogaster gut for functional host-microbe interaction studies. The complete genome is composed of a single circular genome of 2,983,334 bp, with a G+C content of 38%, and a single plasmid of 5,594 bp.

Although the gut microbiome plays important roles in host physiology, health and disease1, we lack understanding of the complex interplay between host genetics and early life environment on the microbial and metabolic composition of the gut. We used the genetically diverse Collaborative Cross mouse system2 to discover that early life history impacts the microbiome composition, whereas dietary changes have only a moderate effect. By contrast, the gut metabolome was shaped mostly by diet, with specific non-dietary metabolites explained by microbial metabolism. Quantitative trait analysis identified mouse genetic trait loci (QTL) that impact the abundances of specific microbes. Human orthologues of genes in the mouse QTL are implicated in gastrointestinal cancer. Additionally, genes located in mouse QTL for Lactobacillales abundance are implicated in arthritis, rheumatic disease and diabetes. Furthermore, Lactobacillales abundance was predictive of higher host T-helper cell counts, suggesting an important link between Lactobacillales and host adaptive immunity.

Although the gut microbiome plays important roles in host physiology, health and disease1, we lack understanding of the complex interplay between host genetics and early life environment on the microbial and metabolic composition of the gut. We used the genetically diverse Collaborative Cross mouse system2 to discover that early life history impacts the microbiome composition, whereas dietary changes have only a moderate effect. By contrast, the gut metabolome was shaped mostly by diet, with specific non-dietary metabolites explained by microbial metabolism. Quantitative trait analysis identified mouse genetic trait loci (QTL) that impact the abundances of specific microbes. Human orthologues of genes in the mouse QTL are implicated in gastrointestinal cancer. Additionally, genes located in mouse QTL for Lactobacillales abundance are implicated in arthritis, rheumatic disease and diabetes. Furthermore, Lactobacillales abundance was predictive of higher host T-helper cell counts, suggesting an important link between Lactobacillales and host adaptive immunity.

Spatial gene expression patterns enable the detection of local covariability and are extremely useful for identifying local gene interactions during normal development. The abundance of spatial expression data in recent years has led to the modeling and analysis of regulatory networks. The inherent complexity of such data makes it a challenge to extract biological information. We developed staNMF, a method that combines a scalable implementation of nonnegative matrix factorization (NMF) with a new stability-driven model selection criterion. When applied to a set ofDrosophilaearly embryonic spatial gene expression images, one of the largest datasets of its kind, staNMF identified 21 principal patterns (PP). Providing a compact yet biologically interpretable representation ofDrosophilaexpression patterns, PP are comparable to a fate map generated experimentally by laser ablation and show exceptional promise as a data-driven alternative to manual annotations. Our analysis mapped genes to cell-fate programs and assigned putative biological roles to uncharacterized genes. Finally, we used the PP to generate local transcription factor regulatory networks. Spatially local correlation networks were constructed for six PP that span along the embryonic anterior-posterior axis. Using a two-tail 5% cutoff on correlation, we reproduced 10 of the 11 links in the well-studied gap gene network. The performance of PP with theDrosophiladata suggests that staNMF provides informative decompositions and constitutes a useful computational lens through which to extract biological insight from complex and often noisy gene expression data.

In this article, we discuss emerging frontiers in RNA biology from a historical perspective. The field is currently undergoing yet another transformative expansion. RNA-seq has revealed that splicing, and, more generally, RNA processing is far more complex than expected, and the mechanisms of regulation are correspondingly sophisticated. Our understanding of the molecular machines involved in RNA metabolism is incomplete and derives from small sample sizes. Even if we manage to complete a catalogue of molecular species, RNA isoforms and the ribonucleoprotein complexes that drive their genesis, the horizons of molecular dynamics and cell-type-specific processing mechanisms await. This is an exciting time to enter into the study of RNA biology; analytical tools, wet and dry, are advancing rapidly, and each new measurement modality brings into view another new function or activity of versatile RNA. Since the dawn of sequence-based RNA biology, we have come a long way.

Steroid hormones induce cascades of gene activation and repression with transformative effects on cell fate . Steroid transduction plays a major role in the development and physiology of nearly all metazoan species, and in the progression of the most common forms of cancer. Despite the paramount importance of steroids in developmental and translational biology, a complete map of transcriptional response has not been developed for any hormone . In the case of 20-hydroxyecdysone (ecdysone) in Drosophila melanogaster, these trajectories range from apoptosis to immortalization. We mapped the ecdysone transduction network in a cohort of 41 cell lines, the largest such atlas yet assembled. We found that the early transcriptional response mirrors the distinctiveness of physiological origins: genes respond in restricted patterns, conditional on the expression levels of dozens of transcription factors. Only a small cohort of genes is constitutively modulated independent of initial cell state. Ecdysone-responsive genes tend to organize into directional same-stranded units, with consecutive genes induced from the same strand. Here, we identify half of the ecdysone receptor heterodimer as the primary rate-limiting step in the response, and find that initial receptor isoform levels modulate the activated cohort of target transcription factors. This atlas of steroid response reveals organizing principles of gene regulation by a model type II nuclear receptor and lays the foundation for comprehensive and predictive understanding of the ecdysone transduction network in the fruit fly.

Evidence has emerged that suggests a link between motor deficits, obesity and many neurological disorders. However, the contributing genetic risk factors are poorly understood. Here we used the Collaborative Cross (CC), a large panel of newly inbred mice that captures 90% of the known variation among laboratory mice, to identify the genetic loci controlling rotarod performance and its relationship with body weight in a cohort of 365 mice across 16 CC strains. Body weight and rotarod performance varied widely across CC strains and were significantly negatively correlated. Genetic linkage analysis identified 14 loci that were associated with body weight. However, 45 loci affected rotarod performance, seven of which were also associated with body weight, suggesting a strong link at the genetic level. Lastly, we show that genes identified in this study overlap significantly with those related to neurological disorders and obesity found in human GWA studies. In conclusion, our results provide a genetic framework for studies of the connection between body weight, the central nervous system and behavior.

Alternative splicing is regulated by RNA binding proteins (RBPs) that recognize pre-mRNA sequence elements and activate or repress adjacent exons. Here, we used RNA interference and RNA-seq to identify splicing events regulated by 56 Drosophila proteins, some previously unknown to regulate splicing. Nearly all proteins affected alternative first exons, suggesting that RBPs play important roles in first exon choice. Half of the splicing events were regulated by multiple proteins, demonstrating extensive combinatorial regulation. We observed that SR and hnRNP proteins tend to act coordinately with each other, not antagonistically. We also identified a cross-regulatory network where splicing regulators affected the splicing of pre-mRNAs encoding other splicing regulators. This large-scale study substantially enhances our understanding of recent models of splicing regulation and provides a resource of thousands of exons that are regulated by 56 diverse RBPs.

In eukaryotic cells, RNAs exist as ribonucleoprotein particles (RNPs). Despite the importance of these complexes in many biological processes, including splicing, polyadenylation, stability, transportation, localization, and translation, their compositions are largely unknown. We affinity-purified 20 distinct RNA-binding proteins (RBPs) from cultured Drosophila melanogaster cells under native conditions and identified both the RNA and protein compositions of these RNP complexes. We identified "high occupancy target" (HOT) RNAs that interact with the majority of the RBPs we surveyed. HOT RNAs encode components of the nonsense-mediated decay and splicing machinery, as well as RNA-binding and translation initiation proteins. The RNP complexes contain proteins and mRNAs involved in RNA binding and post-transcriptional regulation. Genes with the capacity to produce hundreds of mRNA isoforms, ultracomplex genes, interact extensively with heterogeneous nuclear ribonuclear proteins (hnRNPs). Our data are consistent with a model in which subsets of RNPs include mRNA and protein products from the same gene, indicating the widespread existence of auto-regulatory RNPs. From the simultaneous acquisition and integrative analysis of protein and RNA constituents of RNPs, we identify extensive cross-regulatory and hierarchical interactions in post-transcriptional control.

The modENCODE (Model Organism Encyclopedia of DNA Elements) Consortium aimed to map functional elements-including transcripts, chromatin marks, regulatory factor binding sites, and origins of DNA replication-in the model organisms Drosophila melanogaster and Caenorhabditis elegans. During its five-year span, the consortium conducted more than 2,000 genome-wide assays in developmentally staged animals, dissected tissues, and homogeneous cell lines. Analysis of these data sets provided foundational insights into genome, epigenome, and transcriptome structure and the evolutionary turnover of regulatory pathways. These studies facilitated a comparative analysis with similar data types produced by the ENCODE Consortium for human cells. Genome organization differs drastically in these distant species, and yet quantitative relationships among chromatin state, transcription, and cotranscriptional RNA processing are deeply conserved. Of the many biological discoveries of the modENCODE Consortium, we highlight insights that emerged from integrative studies. We focus on operational and scientific lessons that may aid future projects of similar scale or aims in other, emerging model systems.

Recursive splicing is a process in which large introns are removed in multiple steps by re-splicing at ratchet points--5' splice sites recreated after splicing. Recursive splicing was first identified in the Drosophila Ultrabithorax (Ubx) gene and only three additional Drosophila genes have since been experimentally shown to undergo recursive splicing. Here we identify 197 zero nucleotide exon ratchet points in 130 introns of 115 Drosophila genes from total RNA sequencing data generated from developmental time points, dissected tissues and cultured cells. The sequential nature of recursive splicing was confirmed by identification of lariat introns generated by splicing to and from the ratchet points. We also show that recursive splicing is a constitutive process, that depletion of U2AF inhibits recursive splicing, and that the sequence and function of ratchet points are evolutionarily conserved in Drosophila. Finally, we identify four recursively spliced human genes, one of which is also recursively spliced in Drosophila. Together, these results indicate that recursive splicing is commonly used in Drosophila, occurs in humans, and provides insight into the mechanisms by which some large introns are removed.

Drosophila melanogaster plays an important role in molecular, genetic, and genomic studies of heredity, development, metabolism, behavior, and human disease. The initial reference genome sequence reported more than a decade ago had a profound impact on progress in Drosophila research, and improving the accuracy and completeness of this sequence continues to be important to further progress. We previously described improvement of the 117-Mb sequence in the euchromatic portion of the genome and 21 Mb in the heterochromatic portion, using a whole-genome shotgun assembly, BAC physical mapping, and clone-based finishing. Here, we report an improved reference sequence of the single-copy and middle-repetitive regions of the genome, produced using cytogenetic mapping to mitotic and polytene chromosomes, clone-based finishing and BAC fingerprint verification, ordering of scaffolds by alignment to cDNA sequences, incorporation of other map and sequence data, and validation by whole-genome optical restriction mapping. These data substantially improve the accuracy and completeness of the reference sequence and the order and orientation of sequence scaffolds into chromosome arm assemblies. Representation of the Y chromosome and other heterochromatic regions is particularly improved. The new 143.9-Mb reference sequence, designated Release 6, effectively exhausts clone-based technologies for mapping and sequencing. Highly repeat-rich regions, including large satellite blocks and functional elements such as the ribosomal RNA genes and the centromeres, are largely inaccessible to current sequencing and assembly methods and remain poorly represented. Further significant improvements will require sequencing technologies that do not depend on molecular cloning and that produce very long reads.

The transcriptome is the readout of the genome. Identifying common features in it across distant species can reveal fundamental principles. To this end, the ENCODE and modENCODE consortia have generated large amounts of matched RNA-sequencing data for human, worm and fly. Uniform processing and comprehensive annotation of these data allow comparison across metazoan phyla, extending beyond earlier within-phylum transcriptome comparisons and revealing ancient, conserved features. Specifically, we discover co-expression modules shared across animals, many of which are enriched in developmental genes. Moreover, we use expression patterns to align the stages in worm and fly development and find a novel pairing between worm embryo and fly pupae, in addition to the embryo-to-embryo and larvae-to-larvae pairings. Furthermore, we find that the extent of non-canonical, non-coding transcription is similar in each organism, per base pair. Finally, we find in all three organisms that the gene-expression levels, both coding and non-coding, can be quantitatively predicted from chromatin features at the promoter using a 'universal model' based on a single set of organism-independent parameters.

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community.