A large-scale binding and functional map of human RNA-binding proteins.

Many proteins regulate the expression of genes by binding to specific regions encoded in the genome1. Here we introduce a new data set of RNA elements in the human genome that are recognized by RNA-binding proteins (RBPs), generated as part of the Encyclopedia of DNA Elements (ENCODE) project phase III. This class of regulatory elements functions only when transcribed into RNA, as they serve as the binding sites for RBPs that control post-transcriptional processes such as splicing, cleavage and polyadenylation, and the editing, localization, stability and translation of mRNAs. We describe the mapping and characterization of RNA elements recognized by a large collection of human RBPs in K562 and HepG2 cells. Integrative analyses using five assays identify RBP binding sites on RNA and chromatin in vivo, the in vitro binding preferences of RBPs, the function of RBP binding sites and the subcellular localization of RBPs, producing 1,223 replicated data sets for 356 RBPs. We describe the spectrum of RBP binding throughout the transcriptome and the connections between these interactions and various aspects of RNA biology, including RNA stability, splicing regulation and RNA localization. These data expand the catalogue of functional elements encoded in the human genome by the addition of a large set of elements that function at the RNA level by interacting with RBPs.

RNA-guided RNA cleavage by a CRISPR RNA-Cas protein complex.

Compelling evidence indicates that the CRISPR-Cas system protects prokaryotes from viruses and other potential genome invaders. This adaptive prokaryotic immune system arises from the clustered regularly interspaced short palindromic repeats (CRISPRs) found in prokaryotic genomes, which harbor short invader-derived sequences, and the CRISPR-associated (Cas) protein-coding genes. Here, we have identified a CRISPR-Cas effector complex that is comprised of small invader-targeting RNAs from the CRISPR loci (termed prokaryotic silencing (psi)RNAs) and the RAMP module (or Cmr) Cas proteins. The psiRNA-Cmr protein complexes cleave complementary target RNAs at a fixed distance from the 3' end of the integral psiRNAs. In Pyrococcus furiosus, psiRNAs occur in two size forms that share a common 5' sequence tag but have distinct 3' ends that direct cleavage of a given target RNA at two distinct sites. Our results indicate that prokaryotes possess a unique RNA silencing system that functions by homology-dependent cleavage of invader RNAs.

Genome-wide identification of zero nucleotide recursive splicing in Drosophila.

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.

Pelvic pain: What are the symptoms and predictors for surgery, endometriosis and endometriosis severity.

BACKGROUND: Chronic pelvic pain (CPP) is a common condition which significantly impacts the quality of life and wellbeing of many women. Laparoscopy with histopathology is recommended for investigation of pelvic pain and identification of endometriosis with concurrent removal. Never-the-less, the association between endometriosis and pelvic pain is challenging, with endometriosis identified in only 30-50% of women with pain. AIMS: To explore the predictors for undergoing surgery, for identifying endometriosis and endometriosis severity in a cohort of women with CPP. MATERIALS AND METHODS: This study forms part of the Persistent Pelvic Pain project, a prospective observational cohort study (ANZCTR:ACTRN12616000150448). Women referred to a public gynaecology clinic with pain were randomised to one of two gynaecology units for routine care and followed for 36 months with 6-monthly surveys assessing demographics, medical history, quality of life, and pain symptoms measured on a Likert scale. Operative notes were reviewed and endometriosis staged. RESULTS: Of 471 women recruited, 102 women underwent laparoscopy or laparotomy, of whom 52 had endometriosis (n = 37 stage I-II; n = 15 Stage III-IV). Gynaecology unit, pelvic pain intensity and lower parity were all predictors of surgery (odds ratio (OR) 0.342; 95% CI 0.209-0.561; OR 1.303; 95% CI: 1.079-1.573; OR 0.767; 95% CI: 0.620-0.949, respectively). There were no predictors identified for endometriosis diagnosis and the only predictor of severity was increasing age (OR 1.155; 95% CI: 1.047-1.310). CONCLUSIONS: Gynaecology unit and pain intensity were key predictors of undergoing laparoscopy; however, pain severity did not predict endometriosis diagnosis or staging. These findings indicate the need to review current frameworks guiding practice toward surgery for pelvic pain.

Diversity and dynamics of the Drosophila transcriptome.

Animal transcriptomes are dynamic, with each cell type, tissue and organ system expressing an ensemble of transcript isoforms that give rise to substantial diversity. Here we have identified new genes, transcripts and proteins using poly(A)+ RNA sequencing from Drosophila melanogaster in cultured cell lines, dissected organ systems and under environmental perturbations. We found that a small set of mostly neural-specific genes has the potential to encode thousands of transcripts each through extensive alternative promoter usage and RNA splicing. The magnitudes of splicing changes are larger between tissues than between developmental stages, and most sex-specific splicing is gonad-specific. Gonads express hundreds of previously unknown coding and long non-coding RNAs (lncRNAs), some of which are antisense to protein-coding genes and produce short regulatory RNAs. Furthermore, previously identified pervasive intergenic transcription occurs primarily within newly identified introns. The fly transcriptome is substantially more complex than previously recognized, with this complexity arising from combinatorial usage of promoters, splice sites and polyadenylation sites.

Comparative analysis of the transcriptome across distant species.

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.

Crowders Steal Dihydrofolate Reductase Ligands through Quinary Interactions.

Dihydrofolate reductase (DHFR) reduces dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. Due to its role in one carbon metabolism, chromosomal DHFR is the target of the antibacterial drug, trimethoprim. Resistance to trimethoprim has resulted in a type II DHFR that is not structurally related to the chromosomal enzyme target. Because of its metabolic significance, understanding DHFR kinetics and ligand binding behavior in more cell-like conditions, where the total macromolecule concentration can be as great as 300 mg/mL, is important. The progress-curve kinetics and ligand binding properties of the drug target (chromosomal E. coli DHFR) and the drug resistant (R67 DHFR) enzymes were studied in the presence of macromolecular cosolutes. There were varied effects on NADPH oxidation and binding to the two DHFRs, with some cosolutes increasing affinity and others weakening binding. However, DHF binding and reduction in both DHFRs decreased in the presence of all cosolutes. The decreased binding of ligands is mostly attributed to weak associations with the macromolecules, as opposed to crowder effects on the DHFRs. Computer simulations found weak, transient interactions for both ligands with several proteins. The net charge of protein cosolutes correlated with effects on NADP+ binding, with near neutral and positively charged proteins having more detrimental effects on binding. For DHF binding, effects correlated more with the size of binding pockets on the protein crowders. These nonspecific interactions between DHFR ligands and proteins predict that the in vivo efficiency of DHFRs may be much lower than expected from their in vitro rates.

Five-year oncological outcomes after selective neoadjuvant radiotherapy for resectable rectal cancer.

Introduction: There is considerable variation in selection of patients for and type of neoadjuvant radiotherapy administered in the treatment of resectable rectal cancer. The aim of this study was to report outcomes for patients with resected rectal cancer from a unit with step-wise selection for surgery alone, short course radiotherapy (SCRT) or downstaging long course chemoradiotherapy (LCCRT). Material and methods: Cohort analysis of patients with rectal adenocarcinoma resected with curative intent between 2008 and 2012 at a specialist regional colorectal surgery center. The primary endpoints were local recurrence, metastatic recurrence, disease-free survival and overall survival. Exploratory uni- and multi-variable regression analyses were performed to identify predictive factors. Results: About 240 patients were treated by surgery alone, 90 patients received SCRT and 91 patients received LCCRT. Five-year local recurrence was 10.8% in the surgery alone group, 3.3% with SCRT and 18.7% with LCCRT. Metachronous distant metastasis was highest in the SCRT group (13.8% surgery alone, 25.6% SCRT, 15.4% LCCRT). Uni- and multi-variable regression analysis found that local and distant recurrence was attributable predominantly to adverse tumor biology. Conclusions: Patients selected for SCRT had a lower rate of local recurrence than patients selected for surgery alone, but were more likely to develop distant metastasis. There was no difference in overall survival. With low local recurrence rates, distant metastasis is the predominant risk for patients with resectable rectal cancer.

Modulating Enzyme Activity by Altering Protein Dynamics with Solvent.

Optimal enzyme activity depends on a number of factors, including structure and dynamics. The role of enzyme structure is well recognized; however, the linkage between protein dynamics and enzyme activity has given rise to a contentious debate. We have developed an approach that uses an aqueous mixture of organic solvent to control the functionally relevant enzyme dynamics (without changing the structure), which in turn modulates the enzyme activity. Using this approach, we predicted that the hydride transfer reaction catalyzed by the enzyme dihydrofolate reductase (DHFR) from Escherichia coli in aqueous mixtures of isopropanol (IPA) with water will decrease by ∼3 fold at 20% (v/v) IPA concentration. Stopped-flow kinetic measurements find that the pH-independent khydride rate decreases by 2.2 fold. X-ray crystallographic enzyme structures show no noticeable differences, while computational studies indicate that the transition state and electrostatic effects were identical for water and mixed solvent conditions; quasi-elastic neutron scattering studies show that the dynamical enzyme motions are suppressed. Our approach provides a unique avenue to modulating enzyme activity through changes in enzyme dynamics. Further it provides vital insights that show the altered motions of DHFR cause significant changes in the enzyme's ability to access its functionally relevant conformational substates, explaining the decreased khydride rate. This approach has important implications for obtaining fundamental insights into the role of rate-limiting dynamics in catalysis and as well as for enzyme engineering.

Extensive cross-regulation of post-transcriptional regulatory networks in Drosophila.

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.

Regulation of alternative splicing in Drosophila by 56 RNA binding proteins.

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 Vivo Titration of Folate Pathway Enzymes.

How enzymes behave in cells is likely different from how they behave in the test tube. Previous in vitro studies find that osmolytes interact weakly with folate. Removal of the osmolyte from the solvation shell of folate is more difficult than removal of water, which weakens binding of folate to its enzyme partners. To examine if this phenomenon occurs in vivo, osmotic stress titrations were performed with Escherichia coli Two strategies were employed: resistance to an antibacterial drug and complementation of a knockout strain by the appropriate gene cloned into a plasmid that allows tight control of expression levels as well as labeling by a degradation tag. The abilities of the knockout and complemented strains to grow under osmotic stress were compared. Typically, the knockout strain could grow to high osmolalities on supplemented medium, while the complemented strain stopped growing at lower osmolalities on minimal medium. This pattern was observed for an R67 dihydrofolate reductase clone rescuing a ΔfolA strain, for a methylenetetrahydrofolate reductase clone rescuing a ΔmetF strain, and for a serine hydroxymethyltransferase clone rescuing a ΔglyA strain. Additionally, an R67 dihydrofolate reductase clone allowed E. coli DH5α to grow in the presence of trimethoprim until an osmolality of ∼0.81 is reached, while cells in a control titration lacking antibiotic could grow to 1.90 osmol.IMPORTANCEE. coli can survive in drought and flooding conditions and can tolerate large changes in osmolality. However, the cell processes that limit bacterial growth under high osmotic stress conditions are not known. In this study, the dose of four different enzymes in E. coli was decreased by using deletion strains complemented by the gene carried in a tunable plasmid. Under conditions of limiting enzyme concentration (lower than that achieved by chromosomal gene expression), cell growth can be blocked by osmotic stress conditions that are normally tolerated. These observations indicate that E. coli has evolved to deal with variations in its osmotic environment and that normal protein levels are sufficient to buffer the cell from environmental changes. Additional factors involved in the osmotic pressure response may include altered protein concentration/activity levels, weak solute interactions with ligands which can make it more difficult for proteins to bind their substrates/inhibitors/cofactors in vivo, and/or viscosity effects.

Small Angle Neutron Scattering Studies of R67 Dihydrofolate Reductase, a Tetrameric Protein with Intrinsically Disordered N-Termini.

R67 dihydrofolate reductase (DHFR) is a homotetramer with a single active site pore and no sequence or structural homology with chromosomal DHFRs. The R67 enzyme provides resistance to trimethoprim, an active site-directed inhibitor of Escherichia coli DHFR. Sixteen to twenty N-terminal amino acids are intrinsically disordered in the R67 dimer crystal structure. Chymotrypsin cleavage of 16 N-terminal residues results in an active enzyme with a decreased stability. The space sampled by the disordered N-termini of R67 DHFR was investigated using small angle neutron scattering. From a combined analysis using molecular dynamics and the program SASSIE ( http://www.smallangles.net/sassie/SASSIE_HOME.html ), the apoenzyme displays a radius of gyration (Rg) of 21.46 ± 0.50 Å. Addition of glycine betaine, an osmolyte, does not result in folding of the termini as the Rg increases slightly to 22.78 ± 0.87 Å. SASSIE fits of the latter SANS data indicate that the disordered N-termini sample larger regions of space and remain disordered, suggesting they might function as entropic bristles. Pressure perturbation calorimetry also indicated that the volume of R67 DHFR increases upon addition of 10% betaine and decreased at 20% betaine because of the dehydration of the protein. Studies of the hydration of full-length R67 DHFR in the presence of the osmolytes betaine and dimethyl sulfoxide find around 1250 water molecules hydrating the protein. Similar studies with truncated R67 DHFR yield around 400 water molecules hydrating the protein in the presence of betaine. The difference of ∼900 waters indicates the N-termini are well-hydrated.

A Structural Basis for Biguanide Activity.

Metformin is the most commonly prescribed treatment for type II diabetes and related disorders; however, molecular insights into its mode(s) of action have been limited by an absence of structural data. Structural considerations along with a growing body of literature demonstrating its effects on one-carbon metabolism suggest the possibility of folate mimicry and anti-folate activity. Motivated by the growing recognition that anti-diabetic biguanides may act directly upon the gut microbiome, we have determined structures of the complexes formed between the anti-diabetic biguanides (phenformin, buformin, and metformin) and Escherichia coli dihydrofolate reductase (ecDHFR) based on nuclear magnetic resonance, crystallographic, and molecular modeling studies. Interligand Overhauser effects indicate that metformin can form ternary complexes with p-aminobenzoyl-l-glutamate (pABG) as well as other ligands that occupy the region of the folate-binding site that interacts with pABG; however, DHFR inhibition is not cooperative. The biguanides competitively inhibit the activity of ecDHFR, with the phenformin inhibition constant being 100-fold lower than that of metformin. This inhibition may be significant at concentrations present in the gut of treated individuals, and inhibition of DHFR in intestinal mucosal cells may also occur if accumulation levels are sufficient. Perturbation of folate homeostasis can alter the pyridine nucleotide redox ratios that are important regulators of cellular metabolism.

Pathogenicity of modified bat influenza virus with different M genes and its reassortment potential with swine influenza A virus.

In our previous studies, the reassortant virus containing only the PR8 H1N1 matrix (M) gene in the background of the modified bat influenza Bat09 : mH1mN1 virus could be generated. However, whether M genes from other origins can be rescued in the background of the Bat09 : mH1mN1 virus and whether the resulting novel reassortant virus is virulent remain unknown. Herein, two reassortant viruses were generated in the background of the Bat09 : mH1mN1 virus containing either a North American or a Eurasian swine influenza virus M gene. These two reassortant viruses and the reassortant virus with PR8 M as well as the control Bat09 : mH1mN1 virus replicated efficiently in cultured cells, while the reassortant virus with PR8 M grew to a higher titre than the other three viruses in tested cells. Mouse studies showed that reassortant viruses with either North American or Eurasian swine influenza virus M gene did not enhance virulence, whereas the reassortant virus with PR8 M gene displayed higher pathogenicity when compared to the Bat09 : mH1mN1 virus. This is most likely due to the fact that the PR8 H1N1 virus is a mouse-adapted virus. Furthermore, reassortment potential between the Bat09 : mH1mN1 virus and an H3N2 swine influenza virus (A/swine/Texas/4199-2/1998) was investigated using co-infection of Madin-Darby canine kidney cells, but no reassortant viruses were detected. Taken together, our results indicate that the modified bat influenza virus is most likely incapable of reassortment with influenza A viruses with in vitro co-infection experiments, although reassortant viruses with different M genes can be generated by reverse genetics.

Effects of PB1-F2 on the pathogenicity of H1N1 swine influenza virus in mice and pigs.

Although several studies have exploited the effects of PB1-F2 in swine influenza viruses, its contribution to the pathogenicity of swine influenza viruses remains unclear. Herein, we investigated the effects of PB1-F2 on the pathogenicity of influenza virus using a virulent H1N1 A/swine/Kansas/77778/2007 (KS07) virus, which expresses a full-length PB1-F2, in mice and pigs. Using reverse genetics, we generated the wild-type KS07 (KS07_WT), a PB1-F2 knockout mutant (KS07_K/O) and its N66S variant (KS07_N66S). KS07_K/O showed similar pathogenicity in mice to the KS07_WT, whereas KS07_N66S displayed enhanced virulence when compared to the other two viruses. KS07_WT exhibited more efficient replication in lungs and nasal shedding in infected pigs than the other two viruses. Pigs infected with the KS07_WT had higher pulmonary levels of granulocyte-macrophage colony-stimulating factor, IFN-γ, IL-6 and IL-8 at 3 and 5 days post-infection, as well as lower levels of IL-2, IL-4 and IL-12 at 1 day post-infection compared to those infected with the KS07_K/O. These results indicate that PB1-F2 modulates KS07 H1N1 virus replication, pathogenicity and innate immune responses in pigs and the single substitution at position 66 (N/S) in the PB1-F2 plays a critical role in virulence in mice. Taken together, our results provide new insights into the effects of PB1-F2 on the virulence of influenza virus in swine and support PB1-F2 as a virulence factor of influenza A virus in a strain- and host-dependent manner.

Global Analysis of Mouse Polyomavirus Infection Reveals Dynamic Regulation of Viral and Host Gene Expression and Promiscuous Viral RNA Editing.

Mouse polyomavirus (MPyV) lytically infects mouse cells, transforms rat cells in culture, and is highly oncogenic in rodents. We have used deep sequencing to follow MPyV infection of mouse NIH3T6 cells at various times after infection and analyzed both the viral and cellular transcriptomes. Alignment of sequencing reads to the viral genome illustrated the transcriptional profile of the early-to-late switch with both early-strand and late-strand RNAs being transcribed at all time points. A number of novel insights into viral gene expression emerged from these studies, including the demonstration of widespread RNA editing of viral transcripts at late times in infection. By late times in infection, 359 host genes were seen to be significantly upregulated and 857 were downregulated. Gene ontology analysis indicated transcripts involved in translation, metabolism, RNA processing, DNA methylation, and protein turnover were upregulated while transcripts involved in extracellular adhesion, cytoskeleton, zinc finger binding, SH3 domain, and GTPase activation were downregulated. The levels of a number of long noncoding RNAs were also altered. The long noncoding RNA MALAT1, which is involved in splicing speckles and used as a marker in many late-stage cancers, was noticeably downregulated, while several other abundant noncoding RNAs were strongly upregulated. We discuss these results in light of what is currently known about the MPyV life cycle and its effects on host cell growth and metabolism.

Glucose oxidase stabilization against thermal inactivation using high hydrostatic pressure and hydrophobic modification.

High hydrostatic pressure (HHP) stabilized glucose oxidase (GOx) against thermal inactivation. The apparent first-order kinetics of inactivation of GOx were investigated at 0.1-300 MPa and 58.8-80.0°C. At 240 MPa and 74.5°C, GOx inactivated at a rate 50 times slower than at atmospheric pressure at the same temperature. The apparent activation energy of inactivation at 300 MPa was 281.0 ± 17.4 kJ mol-1 or 1.3-fold smaller than for the inactivation at atmospheric pressure (378.1 ± 25.6 kJ mol-1 ). The stabilizing effect of HHP was greatest at 74.5°C, where the activation volume of 57.0 ± 12.0 cm3 mol-1 was highest compared to all other studied temperatures. Positive apparent activation volumes for all the treatment temperatures confirmed that HHP favors GOx stabilization. A second approach to increase GOx stability involved crosslinking with N-(3-dimethylaminopropyl)-N'-ethylcarbodiimide hydrochloride (EDC) and either aniline or benzoate. The modified enzyme remained fully active with only slight increases in KM (1.3-1.9-fold increases for aniline and benzoate modification, respectively). The thermal stability of GOx increased by 8°C with aniline modification, while it decreased by 0.9°C upon modification with benzoate. Biotechnol. Bioeng. 2017;114: 516-525. © 2016 Wiley Periodicals, Inc.

The developmental transcriptome of Drosophila melanogaster.

Drosophila melanogaster is one of the most well studied genetic model organisms; nonetheless, its genome still contains unannotated coding and non-coding genes, transcripts, exons and RNA editing sites. Full discovery and annotation are pre-requisites for understanding how the regulation of transcription, splicing and RNA editing directs the development of this complex organism. Here we used RNA-Seq, tiling microarrays and cDNA sequencing to explore the transcriptome in 30 distinct developmental stages. We identified 111,195 new elements, including thousands of genes, coding and non-coding transcripts, exons, splicing and editing events, and inferred protein isoforms that previously eluded discovery using established experimental, prediction and conservation-based approaches. These data substantially expand the number of known transcribed elements in the Drosophila genome and provide a high-resolution view of transcriptome dynamics throughout development.