Evolutionary analyses of intrinsically disordered regions reveal widespread signals of conservation

Intrinsically disordered regions (IDRs) are segments of proteins without stable three-dimensional structures. As this flexibility allows them to interact with diverse binding partners, IDRs play key roles in cell signaling and gene expression. Despite the prevalence and importance of IDRs in eukaryotic proteomes and various biological processes, associating them with specific molecular functions remains a significant challenge due to their high rates of sequence evolution. However, by comparing the observed values of various IDR-associated properties against those generated under a simulated model of evolution, a recent study found most IDRs across the entire yeast proteome contain conserved features. Furthermore, it showed clusters of IDRs with common "evolutionary signatures," i.e. patterns of conserved features, were associated with specific biological functions. To determine if similar patterns of conservation are found in the IDRs of other systems, in this work we applied a series of phylogenetic models to over 7,500 orthologous IDRs identified in the Drosophila genome to dissect the forces driving their evolution. By comparing models of constrained and unconstrained continuous trait evolution using the Brownian motion and Ornstein-Uhlenbeck models, respectively, we identified signals of widespread constraint, indicating conservation of distributed features is mechanism of IDR evolution common to multiple biological systems. In contrast to the previous study in yeast, however, we observed limited evidence of IDR clusters with specific biological functions, which suggests a more complex relationship between evolutionary constraints and function in the IDRs of multicellular organisms.

Intrinsic protein disorder is insufficient to drive subnuclear clustering in embryonic transcription factors

Modern microscopy has revealed that core nuclear functions, including transcription, replication, and heterochromatin formation, occur in spatially restricted clusters. Previous work from our lab has shown that subnuclear high-concentration clusters of transcription factors may play a role in regulating RNA synthesis in the early Drosophila embryo. A nearly ubiquitous feature of eukaryotic transcription factors is that they contain intrinsically disordered regions (IDRs) that often arise from low complexity amino acid sequences within the protein. It has been proposed that IDRs within transcription factors drive co-localization of transcriptional machinery and target genes into high-concentration clusters within nuclei. Here, we test that hypothesis directly, by conducting a broad survey of the subnuclear localization of IDRs derived from transcription factors. Using a novel algorithm to identify IDRs in the Drosophila proteome, we generated a library of IDRs from transcription factors expressed in the early Drosophila embryo. We used this library to perform a high-throughput imaging screen in Drosophila Schneider-2 (S2) cells. We found that while subnuclear clustering does not occur when the majority of IDRs are expressed alone, it is frequently seen in full-length transcription factors. These results are consistent in live Drosophila embryos, suggesting that IDRs are insufficient to drive the subnuclear clustering behavior of transcription factors. Furthermore, the clustering of transcription factors in living embryos was unaffected by the deletion of IDR sequences. Our results demonstrate that IDRs are unlikely to be the primary molecular drivers of the clustering observed during transcription, suggesting a more complex and nuanced role for these disordered protein sequences.

Single-fly assemblies fill major phylogenomic gaps across the Drosophilidae Tree of Life

Long-read sequencing is driving rapid progress in genome assembly across all major groups of life, including species of the family Drosophilidae, a longtime model system for genetics, genomics, and evolution. We previously developed a cost-effective hybrid Oxford Nanopore (ONT) long-read and Illumina short-read sequencing approach and used it to assemble 101 drosophilid genomes from laboratory cultures, greatly increasing the number of genome assemblies for this taxonomic group. The next major challenge is to address the laboratory culture bias in taxon sampling by sequencing genomes of species that cannot easily be reared in the lab. Here, we build upon our previous methods to perform amplification-free ONT sequencing of single wild flies obtained either directly from the field or from ethanol-preserved specimens in museum collections, greatly improving the representation of lesser studied drosophilid taxa in whole-genome data. Using Illumina Novaseq X Plus and ONT P2 sequencers with R10.4.1 chemistry, we set a new benchmark for inexpensive hybrid genome assembly at US $150 per genome while assembling genomes from as little as 35 ng of genomic DNA from a single fly. We present 183 new genome assemblies for 179 species as a resource for drosophilid systematics, phylogenetics, and comparative genomics. Of these genomes, 62 are from pooled lab strains and 121 from single adult flies. Despite the sample limitations of working with small insects, most single-fly diploid assemblies are comparable in contiguity (>1Mb contig N50), completeness (>98% complete dipteran BUSCOs), and accuracy (>QV40 genome-wide with ONT R10.4.1) to assemblies from inbred lines. We present a well-resolved multi-locus phylogeny for 360 drosophilid and 4 outgroup species encompassing all publicly available (as of August 2023) genomes for this group. Finally, we present a Progressive Cactus whole-genome, reference-free alignment built from a subset of 298 suitably high-quality drosophilid genomes. The new assemblies and alignment, along with updated laboratory protocols and computational pipelines, are released as an open resource and as a tool for studying evolution at the scale of an entire insect family.

Unified bursting strategies in ectopic and endogenous even-skipped expression patterns

Transcription often occurs in bursts as gene promoters switch stochastically between active and inactive states. Enhancers can dictate transcriptional activity in animal development through the modulation of burst frequency, duration, or amplitude. Previous studies observed that different enhancers can achieve a wide range of transcriptional outputs through the same strategies of bursting control. For example, despite responding to different transcription factors, all even-skipped enhancers increase transcription by upregulating burst frequency and amplitude while burst duration remains largely constant. These shared bursting strategies suggest that a unified molecular mechanism constraints how enhancers modulate transcriptional output. Alternatively, different enhancers could have converged on the same bursting control strategy because of natural selection favoring one of these particular strategies. To distinguish between these two scenarios, we compared transcriptional bursting between endogenous and ectopic gene expression patterns. Because enhancers act under different regulatory inputs in ectopic patterns, dissimilar bursting control strategies between endogenous and ectopic patterns would suggest that enhancers adapted their bursting strategies to their trans-regulatory environment. Here, we generated ectopic even-skipped transcription patterns in fruit fly embryos and discovered that bursting strategies remain consistent in endogenous and ectopic even-skipped expression. These results provide evidence for a unified molecular mechanism shaping even-skipped bursting strategies and serve as a starting point to uncover the realm of strategies employed by other enhancers.

Peer review without gatekeeping

eLife is changing its editorial process to emphasize public reviews and assessments of preprints by eliminating accept/reject decisions after peer review.

Impact of essential workers in the context of social distancing for epidemic control

New emerging infectious diseases are identified every year, a subset of which become global pandemics like COVID-19. In the case of COVID-19, many governments have responded to the ongoing pandemic by imposing social policies that restrict contacts outside of the home, resulting in a large fraction of the workforce either working from home or not working. To ensure essential services, however, a substantial number of workers are not subject to these limitations, and maintain many of their pre-intervention contacts. To explore how contacts among such "essential" workers, and between essential workers and the rest of the population, impact disease risk and the effectiveness of pandemic control, we evaluated several mathematical models of essential worker contacts within a standard epidemiology framework. The models were designed to correspond to key characteristics of cashiers, factory employees, and healthcare workers. We find in all three models that essential workers are at substantially elevated risk of infection compared to the rest of the population, as has been documented, and that increasing the numbers of essential workers necessitates the imposition of more stringent controls on contacts among the rest of the population to manage the pandemic. Importantly, however, different archetypes of essential workers differ in both their individual probability of infection and impact on the broader pandemic dynamics, highlighting the need to understand and target intervention for the specific risks faced by different groups of essential workers. These findings, especially in light of the massive human costs of the current COVID-19 pandemic, indicate that contingency plans for future epidemics should account for the impacts of essential workers on disease spread.

Rigorous review and editorial oversight of clinical preprints

Research in many different areas of medicine will benefit from new approaches to peer review and publishing.

A phylogeny for the Drosophila montium species group: A model clade for comparative analyses

The Drosophila montium species group is a clade of 94 named species, closely related to the model species D. melanogaster. The montium species group is distributed over a broad geographic range throughout Asia, Africa, and Australasia. Species of this group possess a wide range of morphologies, mating behaviors, and endosymbiont associations, making this clade useful for comparative analyses. We use genomic data from 42 available species to estimate the phylogeny and relative divergence times within the montium species group, and its relative divergence time from D. melanogaster. To assess the robustness of our phylogenetic inferences, we use 3 non-overlapping sets of 20 single-copy coding sequences and analyze all 60 genes with both Bayesian and maximum likelihood methods. Our analyses support monophyly of the group. Apart from the uncertain placement of a single species, D. baimaii, our analyses also support the monophyly of all seven subgroups proposed within the montium group. Our phylograms and relative chronograms provide a highly resolved species tree, with discordance restricted to estimates of relatively short branches deep in the tree. In contrast, age estimates for the montium crown group, relative to its divergence from D. melanogaster, depend critically on prior assumptions concerning variation in rates of molecular evolution across branches, and hence have not been reliably determined. We discuss methodological issues that limit phylogenetic resolution - even when complete genome sequences are available - as well as the utility of the current phylogeny for understanding the evolutionary and biogeographic history of this clade.

Kinetic sculpting of the seven stripes of the Drosophila even-skipped gene

We used live imaging to visualize the transcriptional dynamics of the Drosophila melanogaster even-skipped gene at single-cell and high-temporal resolution as its seven stripe expression pattern forms, and developed tools to characterize and visualize how transcriptional bursting varies over time and space. We find that despite being created by the independent activity of five enhancers, even-skipped stripes are sculpted by the same kinetic phenomena: a coupled increase of burst frequency and amplitude. By tracking the position and activity of individual nuclei, we show that stripe movement is driven by the exchange of bursting nuclei from the posterior to anterior stripe flanks. Our work provides a conceptual, theoretical and computational framework for dissecting pattern formation in space and time, and reveals how the coordinated transcriptional activity of individual nuclei shapes complex developmental patterns.

Implementing a "publish, then review" model of publishing

From July 2021 eLife will only review manuscripts already published as preprints, and will focus its editorial process on producing public reviews to be posted alongside the preprints.

We need to act now

eLife, like the rest of science, must tackle the many inequalities experienced by Black scientists.

Publishing in the time of COVID-19

eLife is making changes to its policies on peer review in response to the impact of COVID-19 on the scientific community.

Dynamic multifactor hubs interact transiently with sites of active transcription in Drosophila embryos

The regulation of transcription requires the coordination of numerous activities on DNA, yet how transcription factors mediate these activities remains poorly understood. Here, we use lattice light-sheet microscopy to integrate single-molecule and high-speed 4D imaging in developing Drosophila embryos to study the nuclear organization and interactions of the key transcription factors Zelda and Bicoid. In contrast to previous studies suggesting stable, cooperative binding, we show that both factors interact with DNA with surprisingly high off-rates. We find that both factors form dynamic subnuclear hubs, and that Bicoid binding is enriched within Zelda hubs. Remarkably, these hubs are both short lived and interact only transiently with sites of active Bicoid-dependent transcription. Based on our observations, we hypothesize that, beyond simply forming bridges between DNA and the transcription machinery, transcription factors can organize other proteins into hubs that transiently drive multiple activities at their gene targets.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Convergent evolution of gene expression in two high-toothed stickleback populations

Changes in developmental gene regulatory networks enable evolved changes in morphology. These changes can be in cis regulatory elements that act in an allele-specific manner, or changes to the overall trans regulatory environment that interacts with cis regulatory sequences. Here we address several questions about the evolution of gene expression accompanying a convergently evolved constructive morphological trait, increases in tooth number in two independently derived freshwater populations of threespine stickleback fish (Gasterosteus aculeatus). Are convergently evolved cis and/or trans changes in gene expression associated with convergently evolved morphological evolution? Do cis or trans regulatory changes contribute more to gene expression changes accompanying an evolved morphological gain trait? Transcriptome data from dental tissue of ancestral low-toothed and two independently derived high-toothed stickleback populations revealed significantly shared gene expression changes that have convergently evolved in the two high-toothed populations. Comparing cis and trans regulatory changes using phased gene expression data from F1 hybrids, we found that trans regulatory changes were predominant and more likely to be shared among both high-toothed populations. In contrast, while cis regulatory changes have evolved in both high-toothed populations, overall these changes were distinct and not shared among high-toothed populations. Together these data suggest that a convergently evolved trait can occur through genetically distinct regulatory changes that converge on similar trans regulatory environments.

Gene family innovation, conservation and loss on the animal stem lineage

Choanoflagellates, the closest living relatives of animals, can provide unique insights into the changes in gene content that preceded the origin of animals. However, only two choanoflagellate genomes are currently available, providing poor coverage of their diversity. We sequenced transcriptomes of 19 additional choanoflagellate species to produce a comprehensive reconstruction of the gains and losses that shaped the ancestral animal gene repertoire. We identified ~1944 gene families that originated on the animal stem lineage, of which only 39 are conserved across all animals in our study. In addition, ~372 gene families previously thought to be animal-specific, including Notch, Delta, and homologs of the animal Toll-like receptor genes, instead evolved prior to the animal-choanoflagellate divergence. Our findings contribute to an increasingly detailed portrait of the gene families that defined the biology of the Urmetazoan and that may underpin core features of extant animals.

All animals, from sea sponges and reef-building corals to elephants and humans, share a single common ancestor that lived over half a billion years ago. This single-celled predecessor evolved the ability to develop into a creature made up of many cells with specialized jobs. Reconstructing the steps in this evolutionary process has been difficult because the earliest animals were soft-bodied and microscopic and did not leave behind fossils that scientists can study.

Though their bodies have since disintegrated, many of the instructions for building the first animals live on in genes that were passed on to life forms that still exist. Scientists are trying to retrace those genes back to the first animal by comparing the genomes of living animals with their closest relatives, the choanoflagellates. Choanoflagellates are single-celled, colony-forming organisms that live in waters around the world. Comparisons with choanoflagellates may help scientists identify which genes were necessary to help animals evolve and diversify into so many different species. So far, 1,000 animal and two choanoflagellate genomes have been sequenced. But the gene repertoires of most species of choanoflagellates have yet to be analyzed.

Now, Richter et al. have cataloged the genes of 19 more species of choanoflagellates. This added information allowed them to recreate the likely gene set of the first animal and to identify genetic changes that occurred during animal evolution. The analyses showed that modern animals lost about a quarter of the genes present in their last common ancestor with choanoflagellates and gained an equal number of new genes. Richter et al. identified several dozen core animal genes that were gained and subsequently preserved throughout animal evolution. Many of these are necessary so that an embryo can develop properly, but the precise roles of some core genes remain a mystery. Most other genes that emerged in the first animals have been lost in at least one living animal.

The study of Richter et al. also showed that some very important genes in animals, including genes essential for early development and genes that help the immune system detect pathogens, predate animals. These key genes trace back to animals’ last common ancestor with choanoflagellates and may have evolved new roles in animals.

The ecology of the Drosophila-yeast mutualism in wineries

The fruit fly, Drosophila melanogaster, is preferentially found on fermenting fruits. The yeasts that dominate the microbial communities of these substrates are the primary food source for developing D. melanogaster larvae, and adult flies manifest a strong olfactory system-mediated attraction for the volatile compounds produced by these yeasts during fermentation. Although most work on this interaction has focused on the standard laboratory yeast Saccharomyces cerevisiae, a wide variety of other yeasts naturally ferment fallen fruit. Here we address the open question of whether D. melanogaster preferentially associates with distinct yeasts in different, closely-related environments. We characterized the spatial and temporal dynamics of Drosophila-associated fungi in Northern California wineries that use organic grapes and natural fermentation using high-throughput, short-amplicon sequencing. We found that there is nonrandom structure in the fungal communities that are vectored by flies both between and within vineyards. Within wineries, the fungal communities associated with flies in cellars, fermentation tanks, and pomace piles are distinguished by varying abundances of a small number of yeast species. To investigate the origins of this structure, we assayed Drosophila attraction to, oviposition on, larval development in, and longevity when consuming the yeasts that distinguish vineyard microhabitats from each other. We found that wild fly lines did not respond differentially to the yeast species that distinguish winery habitats in habitat specific manner. Instead, this subset of yeast shares traits that make them attractive to and ensure their close association with Drosophila.

Patterns of chromatin accessibility along the anterior-posterior axis in the early Drosophila embryo

As the Drosophila embryo transitions from the use of maternal RNAs to zygotic transcription, domains of open chromatin, with relatively low nucleosome density and specific histone marks, are established at promoters and enhancers involved in patterned embryonic transcription. However it remains unclear how regions of activity are established during early embryogenesis, and if they are the product of spatially restricted or ubiquitous processes. To shed light on this question, we probed chromatin accessibility across the anterior-posterior axis (A-P) of early Drosophila melanogaster embryos by applying a transposon based assay for chromatin accessibility (ATAC-seq) to anterior and posterior halves of hand-dissected, cellular blastoderm embryos. We find that genome-wide chromatin accessibility is highly similar between the two halves, with regions that manifest significant accessibility in one half of the embryo almost always accessible in the other half, even for promoters that are active in exclusively one half of the embryo. These data support previous studies that show that chromatin accessibility is not a direct result of activity, and point to a role for ubiquitous factors or processes in establishing chromatin accessibility at promoters in the early embryo. However, in concordance with similar works, we find that at enhancers active exclusively in one half of the embryo, we observe a significant skew towards greater accessibility in the region of their activity, highlighting the role of patterning factors such as Bicoid in this process.

Rapid Global Spread of wRi-like Wolbachia across Multiple Drosophila

Maternally transmitted Wolbachia, Spiroplasma, and Cardinium bacteria are common in insects [1], but their interspecific spread is poorly understood. Endosymbionts can spread rapidly within host species by manipulating host reproduction, as typified by the global spread of wRi Wolbachia observed in Drosophila simulans [2, 3]. However, because Wolbachia cannot survive outside host cells, spread between distantly related host species requires horizontal transfers that are presumably rare [4-7]. Here, we document spread of wRi-like Wolbachia among eight highly diverged Drosophila hosts (10-50 million years) over only about 14,000 years (5,000-27,000). Comparing 110 wRi-like genomes, we find ≤0.02% divergence from the wRi variant that spread rapidly through California populations of D. simulans. The hosts include both globally invasive species (D. simulans, D. suzukii, and D. ananassae) and narrowly distributed Australian endemics (D. anomalata and D. pandora) [8]. Phylogenetic analyses that include mtDNA genomes indicate introgressive transfer of wRi-like Wolbachia between closely related species D. ananassae, D. anomalata, and D. pandora but no horizontal transmission within species. Our analyses suggest D. ananassae as the Wolbachia source for the recent wRi invasion of D. simulans and D. suzukii as the source of Wolbachia in its sister species D. subpulchrella. Although six of these wRi-like variants cause strong cytoplasmic incompatibility, two cause no detectable reproductive effects, indicating that pervasive mutualistic effects [9, 10] complement the reproductive manipulations for which Wolbachia are best known. "Super spreader" variants like wRi may be particularly useful for controlling insect pests and vector-borne diseases with Wolbachia transinfections [11].

Single Molecule Imaging in Live Embryos Using Lattice Light-Sheet Microscopy

In the past decade, live-cell single molecule imaging studies have provided unique insights on how DNA-binding molecules such as transcription factors explore the nuclear environment to search for and bind to their targets. However, due to technological limitations, single molecule experiments in living specimens have largely been limited to monolayer cell cultures. Lattice light-sheet microscopy overcomes these limitations and has now enabled single molecule imaging within thicker specimens such as embryos. Here we describe a general procedure to perform single molecule imaging in living Drosophila melanogaster embryos using lattice light-sheet microscopy. This protocol allows direct observation of both transcription factor diffusion and binding dynamics. Finally, we illustrate how this Drosophila protocol can be extended to other thick samples using single molecule imaging in live mouse embryos as an example.

Convergence of topological domain boundaries, insulators, and polytene interbands revealed by high-resolution mapping of chromatin contacts in the early Drosophila melanogaster embryo

High-throughput assays of three-dimensional interactions of chromosomes have shed considerable light on the structure of animal chromatin. Despite this progress, the precise physical nature of observed structures and the forces that govern their establishment remain poorly understood. Here we present high resolution Hi-C data from early Drosophila embryos. We demonstrate that boundaries between topological domains of various sizes map to DNA elements that resemble classical insulator elements: short genomic regions sensitive to DNase digestion that are strongly bound by known insulator proteins and are frequently located between divergent promoters. Further, we show a striking correspondence between these elements and the locations of mapped polytene interband regions. We believe it is likely this relationship between insulators, topological boundaries, and polytene interbands extends across the genome, and we therefore propose a model in which decompaction of boundary-insulator-interband regions drives the organization of interphase chromosomes by creating stable physical separation between adjacent domains.

A molecular atlas of the developing ectoderm defines neural, neural crest, placode, and nonneural progenitor identity in vertebrates

During vertebrate neurulation, the embryonic ectoderm is patterned into lineage progenitors for neural plate, neural crest, placodes and epidermis. Here, we use Xenopus laevis embryos to analyze the spatial and temporal transcriptome of distinct ectodermal domains in the course of neurulation, during the establishment of cell lineages. In order to define the transcriptome of small groups of cells from a single germ layer and to retain spatial information, dorsal and ventral ectoderm was subdivided along the anterior-posterior and medial-lateral axes by microdissections. Principal component analysis on the transcriptomes of these ectoderm fragments primarily identifies embryonic axes and temporal dynamics. This provides a genetic code to define positional information of any ectoderm sample along the anterior-posterior and dorsal-ventral axes directly from its transcriptome. In parallel, we use nonnegative matrix factorization to predict enhanced gene expression maps onto early and mid-neurula embryos, and specific signatures for each ectoderm area. The clustering of spatial and temporal datasets allowed detection of multiple biologically relevant groups (e.g., Wnt signaling, neural crest development, sensory placode specification, ciliogenesis, germ layer specification). We provide an interactive network interface, EctoMap, for exploring synexpression relationships among genes expressed in the neurula, and suggest several strategies to use this comprehensive dataset to address questions in developmental biology as well as stem cell or cancer research.

Vertebrate embryo germ layers become progressively regionalized by evolutionarily conserved molecular processes. Catching the early steps of this dynamic spatial cell diversification at the scale of the transcriptome was challenging, even with the advent of efficient RNA sequencing. We have microdissected complementary and defined areas of a single germ layer, the developing ectoderm, and explored how the transcriptome changes over time and space in the ectoderm during the differentiation of frog epidermis, neural plate, and neural crest. We have created EctoMap, a searchable interface using these regional transcriptomes, to predict the expression of the 31 thousand genes expressed in neurulae and their networks of co-expression, predictive of functional relationships. Through several examples, we illustrate how these data provide insights in development, cancer, evolution and stem cell biology.

Dense Bicoid hubs accentuate binding along the morphogen gradient

Morphogen gradients direct the spatial patterning of developing embryos; however, the mechanisms by which these gradients are interpreted remain elusive. Here we used lattice light-sheet microscopy to perform in vivo single-molecule imaging in early Drosophila melanogaster embryos of the transcription factor Bicoid that forms a gradient and initiates patterning along the anteroposterior axis. In contrast to canonical models, we observed that Bicoid binds to DNA with a rapid off rate throughout the embryo such that its average occupancy at target loci is on-rate-dependent. We further observed Bicoid forming transient "hubs" of locally high density that facilitate binding as factor levels drop, including in the posterior, where we observed Bicoid binding despite vanishingly low protein levels. We propose that localized modulation of transcription factor on rates via clustering provides a general mechanism to facilitate binding to low-affinity targets and that this may be a prevalent feature of other developmental transcription factors.

Patterns of Genome-Wide Diversity and Population Structure in the Drosophila athabasca Species Complex

The Drosophila athabasca species complex contains three recently diverged, prezygotically isolated semispecies (Western-Northern, Eastern-A, and Eastern-B) that are distributed across North America and share zones of sympatry. Inferences based on a handful of loci suggest that this complex might be an ideal system for studying the genetics of incipient speciation and the evolution of prezygotic isolating mechanisms, but patterns of differentiation have not been characterized systematically. Here, we assembled a draft genome for D. athabasca and analyze whole-genome re-sequencing data for 28 individuals from across the species range to characterize genome-wide patterns of diversity and population differentiation among semispecies. Patterns of differentiation on the X-chromosome vs. autosomes vary, with the X-chromosome showing better phylogenetic resolution and increased levels of between semispecies divergence. Despite low levels of overall differentiation and a lack of phylogenetic resolution of the autosomes for the most closely related semispecies, individuals do exhibit distinct genetic clustering. Demographic analyses provide some support for a model of isolation with migration within D. athabasca, with divergence times <20 kya. The young divergence times of the semispecies of D. athabasca, together with strong levels of sexual isolation, makes them a promising system for studying the evolution of prezygotic isolation and speciation.

Genome-wide measurement of spatial expression in patterning mutants of Drosophila melanogaster

Patterning in the Drosophila melanogaster embryo is affected by multiple maternal factors, but the effect of these factors on spatial gene expression has not been systematically analyzed. Here we characterize the effect of the maternal factors Zelda, Hunchback and Bicoid by cryosectioning wildtype and mutant blastoderm stage embryos and sequencing mRNA from each slice. The resulting atlas of spatial gene expression highlights the intersecting roles of these factors in regulating spatial patterns, and serves as a resource for researchers studying spatial patterning in the early embryo. We identify a large number of genes with both expected and unexpected patterning changes, and through integrated analysis of transcription factor binding data identify common themes in genes with complex dependence on these transcription factors.

Stable Host Gene Expression in the Gut of Adult Drosophila melanogaster with Different Bacterial Mono-Associations

There is growing evidence that the microbes found in the digestive tracts of animals influence host biology, but we still do not understand how they accomplish this. Here, we evaluated how different microbial species commonly associated with laboratory-reared Drosophila melanogaster impact host biology at the level of gene expression in the dissected adult gut and in the entire adult organism. We observed that guts from animals associated from the embryonic stage with either zero, one or three bacterial species demonstrated indistinguishable transcriptional profiles. Additionally, we found that the gut transcriptional profiles of animals reared in the presence of the yeast Saccharomyces cerevisiae alone or in combination with bacteria could recapitulate those of conventionally-reared animals. In contrast, we found whole body transcriptional profiles of conventionally-reared animals were distinct from all of the treatments tested. Our data suggest that adult flies are insensitive to the ingestion of the bacteria found in their gut, but that prior to adulthood, different microbes impact the host in ways that lead to global transcriptional differences observable across the whole adult body.

The genome of the crustacean Parhyale hawaiensis, a model for animal development, regeneration, immunity and lignocellulose digestion

The amphipod crustacean Parhyale hawaiensis is a blossoming model system for studies of developmental mechanisms and more recently regeneration. We have sequenced the genome allowing annotation of all key signaling pathways, transcription factors, and non-coding RNAs that will enhance ongoing functional studies. Parhyale is a member of the Malacostraca clade, which includes crustacean food crop species. We analysed the immunity related genes of Parhyale as an important comparative system for these species, where immunity related aquaculture problems have increased as farming has intensified. We also find that Parhyale and other species within Multicrustacea contain the enzyme sets necessary to perform lignocellulose digestion ('wood eating'), suggesting this ability may predate the diversification of this lineage. Our data provide an essential resource for further development of Parhyale as an experimental model. The first malacostracan genome will underpin ongoing comparative work in food crop species and research investigating lignocellulose as an energy source.

DOI:http://dx.doi.org/10.7554/eLife.20062.001

The marine crustacean known as Parhyale hawaiensis is related to prawns, shrimps and crabs and is found at tropical coastlines around the world. This species has recently attracted scientific interest as a possible new model to study how animal embryos develop before birth and, because Parhyale can rapidly regrow lost limbs, how tissues and organs regenerate. Indeed, Parhyale has many characteristics that make it a good model organism, being small, fast-growing and easy to keep and care for in the laboratory.

Several research tools have already been developed to make it easier to study Parhyale. This includes the creation of a system for using the popular gene editing technology, CRISPR, in this animal. However, one critical resource that is available for most model organisms was missing; the complete sequence of all the genetic information of this crustacean, also known as its genome, was not available.

Kao, Lai, Stamataki et al. have now compiled the Parhyale genome – which is slightly larger than the human genome – and studied its genetics. Analysis revealed that Parhyale has genes that allow it to fully digest plant material. This is unusual because most animals that do this rely upon the help of bacteria. Kao, Lai, Stamataki et al. also identified genes that provide some of the first insights into the immune system of crustaceans, which protects these creatures from diseases.

Kao, Lai, Stamataki et al. have provided a resource and findings that could help to establish Parhyale as a popular model organism for studying several ideas in biology, including organ regeneration and embryonic development. Understanding how Parhyale digests plant matter, for example, could progress the biofuel industry towards efficient production of greener energy. Insights from its immune system could also be adapted to make farmed shrimp and prawns more resistant to infections, boosting seafood production.

DOI:http://dx.doi.org/10.7554/eLife.20062.002

Oxygen changes drive non-uniform scaling in Drosophila melanogaster embryogenesis

We previously demonstrated that, while changes in temperature produce dramatic shifts in the time elapsed during Drosophila melanogaster embryogenesis, the relative timing of events within embryogenesis does not change. However, it was unclear if this uniform scaling is an intrinsic property of developing embryos, or if it is specific to thermal fluctuations. To investigate this, here we characterize the embryonic response to changes in oxygen concentration, which also impact developmental rate, using time-lapse imaging, and find it fundamentally different from the temperature response. Most notably, changes in oxygen levels drive developmental heterochrony, with the timing of several morphological processes showing distinct scaling behaviors. Gut formation is severely slowed by decreases in oxygen, while head involution and syncytial development are less impacted than the rest of development, and the order of several developmental landmarks is inverted at different oxygen levels. These data reveal that the uniform scaling seen with changes in temperature is not a trivial consequence of adjusting developmental rate. The developmental rate changes produced by changing oxygen concentrations dwarf those induced by temperature, and greatly impact survival. While extreme temperatures increase early embryo mortality, mild hypoxia increases arrest and death during mid-embryogenesis and mild hyperoxia increases survival over normoxia.

Sex Bias and Maternal Contribution to Gene Expression Divergence in Drosophila Blastoderm Embryos

Early embryogenesis is a unique developmental stage where genetic control of development is handed off from mother to zygote. Yet the contribution of this transition to the evolution of gene expression is poorly understood. Here we study two aspects of gene expression specific to early embryogenesis in Drosophila: sex-biased gene expression prior to the onset of canonical X chromosomal dosage compensation, and the contribution of maternally supplied mRNAs. We sequenced mRNAs from individual unfertilized eggs and precisely staged and sexed blastoderm embryos, and compared levels between D. melanogaster, D. yakuba, D. pseudoobscura and D. virilis. First, we find that mRNA content is highly conserved for a given stage and that studies relying on pooled embryos likely systematically overstate the degree of gene expression divergence. Unlike studies done on larvae and adults where most species show a larger proportion of genes with male-biased expression, we find that transcripts in Drosophila embryos are largely female-biased in all species, likely due to incomplete dosage compensation prior to the activation of the canonical dosage compensation mechanism. The divergence of sex-biased gene expression across species is observed to be often due to lineage-specific decrease of expression; the most drastic example of which is the overall reduction of male expression from the neo-X chromosome in D. pseudoobscura, leading to a pervasive female-bias on this chromosome. We see no evidence for a faster evolution of expression on the X chromosome in embryos (no “faster-X” effect), unlike in adults, and contrary to a previous study on pooled non-sexed embryos. Finally, we find that most genes are conserved in regard to their maternal or zygotic origin of transcription, and present evidence that differences in maternal contribution to the blastoderm transcript pool may be due to species-specific divergence of transcript degradation rates.

The early embryonic stage of development is unique as this stage begins under the genetic control of the mother and transitions to genetic control of the zygote. In order to investigate how these developmental specificities have shaped gene expression evolution, we studied the contribution of maternal gene expression to the pool of transcripts in the early embryo and how gene expression differs in early female and male embryos. We find extensive female-bias in gene expression, something not observed in adults. This is largely due to the delay in the establishment of a mechanism that balances gene expression from the differing number of X chromosomes in females and males. The maternal contribution to the transcript pool in the embryo is largely conserved over species, with some exceptions. However, there may also be some species-specific differences in the degradation rate of maternal transcripts, an often-overlooked contributor to variation in transcript level.

Low-cost, low-input RNA-seq protocols perform nearly as well as high-input protocols

Recently, a number of protocols extending RNA-sequencing to the single-cell regime have been published. However, we were concerned that the additional steps to deal with such minute quantities of input sample would introduce serious biases that would make analysis of the data using existing approaches invalid. In this study, we performed a critical evaluation of several of these low-volume RNA-seq protocols, and found that they performed slightly less well in per-gene linearity of response, but with at least two orders of magnitude less sample required. We also explored a simple modification to one of these protocols that, for many samples, reduced the cost of library preparation to approximately $20/sample.

Comparative validation of the D. melanogaster modENCODE transcriptome annotation

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.

Saccharomyces cerevisiae mitochondria are required for optimal attractiveness to Drosophila melanogaster

While screening a large collection of wild and laboratory yeast strains for their ability to attract Drosophila melanogaster adults, we noticed a large difference in fly preference for two nearly isogenic strains of Saccharomyces cerevisiae, BY4741 and BY4742. Using standard genetic analyses, we tracked the preference difference to the lack of mitochondria in the BY4742 strain used in the initial experiment. We used gas chromatography coupled with mass spectroscopy to examine the volatile compounds produced by BY4741 and the mitochondria-deficient BY4742, and found that they differed significantly. We observed that several ethyl esters are present at much higher levels in strains with mitochondria, even in fermentative conditions. We found that nitrogen levels in the substrate affect the production of these compounds, and that they are produced at the highest level by strains with mitochondria when fermenting natural fruit substrates. Collectively these observations demonstrate that core metabolic processes mediate the interaction between yeasts and insect vectors, and highlight the importance mitochondrial functions in yeast ecology.

Characterization of the transcriptome, nucleotide sequence polymorphism, and natural selection in the desert adapted mouse Peromyscus eremicus

As a direct result of intense heat and aridity, deserts are thought to be among the most harsh of environments, particularly for their mammalian inhabitants. Given that osmoregulation can be challenging for these animals, with failure resulting in death, strong selection should be observed on genes related to the maintenance of water and solute balance. One such animal, Peromyscus eremicus, is native to the desert regions of the southwest United States and may live its entire life without oral fluid intake. As a first step toward understanding the genetics that underlie this phenotype, we present a characterization of the P. eremicus transcriptome. We assay four tissues (kidney, liver, brain, testes) from a single individual and supplement this with population level renal transcriptome sequencing from 15 additional animals. We identified a set of transcripts undergoing both purifying and balancing selection based on estimates of Tajima's D. In addition, we used the branch-site test to identify a transcript-Slc2a9, likely related to desert osmoregulation-undergoing enhanced selection in P. eremicus relative to a set of related non-desert rodents.

Establishment of regions of genomic activity during the Drosophila maternal to zygotic transition

We describe the genome-wide distributions and temporal dynamics of nucleosomes and post-translational histone modifications throughout the maternal-to-zygotic transition in embryos of Drosophila melanogaster. At mitotic cycle 8, when few zygotic genes are being transcribed, embryonic chromatin is in a relatively simple state: there are few nucleosome free regions, undetectable levels of the histone methylation marks characteristic of mature chromatin, and low levels of histone acetylation at a relatively small number of loci. Histone acetylation increases by cycle 12, but it is not until cycle 14 that nucleosome free regions and domains of histone methylation become widespread. Early histone acetylation is strongly associated with regions that we have previously shown to be bound in early embryos by the maternally deposited transcription factor Zelda, suggesting that Zelda triggers a cascade of events, including the accumulation of specific histone modifications, that plays a role in the subsequent activation of these sequences.

Structural variation among wild and industrial strains of Penicillium chrysogenum

Strain selection and strain improvement are the first, and arguably most important, steps in the industrial production of biological compounds by microorganisms. While traditional methods of mutagenesis and selection have been effective in improving production of compounds at a commercial scale, the genetic changes underpinning the altered phenotypes have remained largely unclear. We utilized high-throughput Illumina short read sequencing of a wild Penicillium chrysogenum strain in order to make whole genome comparisons to a sequenced improved strain (WIS 54–1255). We developed an assembly-free method of identifying chromosomal rearrangements and validated the in silico predictions with a PCR-based assay and Sanger sequencing. Despite many rounds of mutagen treatment and artificial selection, WIS 54–1255 differs from its wild progenitor at only one of the identified rearrangements. We suggest that natural variants predisposed for high penicillin production were instrumental in the success of WIS 54–1255 as an industrial strain. In addition to finding a previously published inversion in the penicillin biosynthesis cluster, we located several genes related to penicillin production associated with these rearrangements. By comparing the configuration of rearrangement events among several historically important strains known to be high penicillin producers to a collection of recently isolated wild strains, we suggest that wild strains with rearrangements similar to those in known high penicillin producers may be viable candidates for further improvement efforts.

Drosophila embryogenesis scales uniformly across temperature in developmentally diverse species

Temperature affects both the timing and outcome of animal development, but the detailed effects of temperature on the progress of early development have been poorly characterized. To determine the impact of temperature on the order and timing of events during Drosophila melanogaster embryogenesis, we used time-lapse imaging to track the progress of embryos from shortly after egg laying through hatching at seven precisely maintained temperatures between 17.5°C and 32.5°C. We employed a combination of automated and manual annotation to determine when 36 milestones occurred in each embryo. D. melanogaster embryogenesis takes 33 hours at 17.5°C, and accelerates with increasing temperature to a low of 16 hours at 27.5°C, above which embryogenesis slows slightly. Remarkably, while the total time of embryogenesis varies over two fold, the relative timing of events from cellularization through hatching is constant across temperatures. To further explore the relationship between temperature and embryogenesis, we expanded our analysis to cover ten additional Drosophila species of varying climatic origins. Six of these species, like D. melanogaster, are of tropical origin, and embryogenesis time at different temperatures was similar for them all. D. mojavensis, a sub-tropical fly, develops slower than the tropical species at lower temperatures, while D. virilis, a temperate fly, exhibits slower development at all temperatures. The alpine sister species D. persimilis and D. pseudoobscura develop as rapidly as tropical flies at cooler temperatures, but exhibit diminished acceleration above 22.5°C and have drastically slowed development by 30°C. Despite ranging from 13 hours for D. erecta at 30°C to 46 hours for D. virilis at 17.5°C, the relative timing of events from cellularization through hatching is constant across all species and temperatures examined here, suggesting the existence of a previously unrecognized timer controlling the progress of embryogenesis that has been tuned by natural selection as each species diverges.

Temperature profoundly impacts the rate of development of “cold-blooded” animals, which proceeds far faster when it is warm. There is, however, no universal relationship. Closely related species can develop at markedly different speeds at the same temperature. This creates a major challenge when comparing development among species, as it is unclear whether they should be compared at the same temperature or under different conditions to maintain the same developmental rate. Facing this challenge while working with flies (Drosophila species), we found there was little data to inform this decision. So, using time-lapse imaging, precise temperature-control, and computational and manual video-analysis, we tracked the complex process of embryogenesis in 11 species at seven different temperatures. There was over a three-fold difference in developmental rate between the fastest species at its fastest temperature and the slowest species at its slowest temperature. However, our finding that the timing of events within development all scaled uniformly across species and temperatures astonished us. This is good news for developmental biologists, since we can induce species to develop nearly identically by growing them at different temperatures. But it also means flies must possess some unknown clock-like molecular mechanism driving embryogenesis forward.

Sex-specific embryonic gene expression in species with newly evolved sex chromosomes

Sex chromosome dosage differences between females and males are a significant form of natural genetic variation in many species. Like many species with chromosomal sex determination, Drosophila females have two X chromosomes, while males have one X and one Y. Fusions of sex chromosomes with autosomes have occurred along the lineage leading to D. pseudoobscura and D. miranda. The resulting neo-sex chromosomes are gradually evolving the properties of sex chromosomes, and neo-X chromosomes are becoming targets for the molecular mechanisms that compensate for differences in X chromosome dose between sexes. We have previously shown that D. melanogaster possess at least two dosage compensation mechanisms: the well- characterized MSL-mediated dosage compensation active in most somatic tissues, and another system active during early embryogenesis prior to the onset of MSL-mediated dosage compensation. To better understand the developmental constraints on sex chromosome gene expression and evolution, we sequenced mRNA from individual male and female embryos of D. pseudoobscura and D. miranda, from ∼0.5 to 8 hours of development. Autosomal expression levels are highly conserved between these species. But, unlike D. melanogaster, we observe a general lack of dosage compensation in D. pseudoobscura and D. miranda prior to the onset of MSL-mediated dosage compensation. Thus, either there has been a lineage-specific gain or loss in early dosage compensation mechanism(s) or increasing X chromosome dose may strain dosage compensation systems and make them less effective. The extent of female bias on the X chromosomes decreases through developmental time with the establishment of MSL-mediated dosage compensation, but may do so more slowly in D. miranda than D. pseudoobscura. These results also prompt a number of questions about whether species with more sex-linked genes have more sex-specific phenotypes, and how much transcript level variance is tolerable during critical stages of development.

Many animals have sex-specific combinations of chromosomes. In humans, for example, females have two X chromosomes while males have one X and one Y. In most species with XX:XY systems, the Y chromosome is degenerate and gene-poor while the X encodes a large number of functional genes. A variety of systems have evolved to ensure that males with one X chromosome and females with two X chromosomes have the same gene expression level for X-linked genes. The vinegar fly D. melanogaster has at least two dosage compensation systems: one that acts early in development, and another active in later stages. In this paper, we determine expression levels for thousands of genes in male and female embryos at different developmental stages in two species, D. pseudoobscura and D. miranda, that have unusually large fractions of their genomes in X or X-like chromosomes. We show that dosage compensation is established slowly during embryogenesis, and that in these species, dosage compensation appears to be absent in early development. This may be due to a lineage-specific loss or gain of compensation mechanism, or possibly because the machinery of dosage compensation cannot effectively handle the increased demand in these species.

Extensive divergence of transcription factor binding in Drosophila embryos with highly conserved gene expression

To better characterize how variation in regulatory sequences drives divergence in gene expression, we undertook a systematic study of transcription factor binding and gene expression in blastoderm embryos of four species, which sample much of the diversity in the 40 million-year old genus Drosophila: D. melanogaster, D. yakuba, D. pseudoobscura and D. virilis. We compared gene expression, measured by mRNA-seq, to the genome-wide binding, measured by ChIP-seq, of four transcription factors involved in early anterior-posterior patterning. We found that mRNA levels are much better conserved than individual transcription factor binding events, and that changes in a gene's expression were poorly explained by changes in adjacent transcription factor binding. However, highly bound sites, sites in regions bound by multiple factors and sites near genes are conserved more frequently than other binding, suggesting that a considerable amount of transcription factor binding is weakly or non-functional and not subject to purifying selection.

Improving transcriptome assembly through error correction of high-throughput sequence reads

The study of functional genomics, particularly in non-model organisms, has been dramatically improved over the last few years by the use of transcriptomes and RNAseq. While these studies are potentially extremely powerful, a computationally intensive procedure, the de novo construction of a reference transcriptome must be completed as a prerequisite to further analyses. The accurate reference is critically important as all downstream steps, including estimating transcript abundance are critically dependent on the construction of an accurate reference. Though a substantial amount of research has been done on assembly, only recently have the pre-assembly procedures been studied in detail. Specifically, several stand-alone error correction modules have been reported on and, while they have shown to be effective in reducing errors at the level of sequencing reads, how error correction impacts assembly accuracy is largely unknown. Here, we show via use of a simulated and empiric dataset, that applying error correction to sequencing reads has significant positive effects on assembly accuracy, and should be applied to all datasets. A complete collection of commands which will allow for the production of Reptile corrected reads is available at https://github.com/macmanes/error_correction/tree/master/scripts and as File S1.

An essential role for zygotic expression in the pre-cellular Drosophila embryo

The Drosophila embryo proceeds through thirteen mitotic divisions as a syncytium. Its nuclei distribute in the embryo's interior during the first six divisions, dividing synchronously with a cycle time of less than ten minutes. After seven divisions (nuclear cycle 8), the syncytial blastoderm forms as the nuclei approach the embryo surface and slow their cycle time; subsequent divisions proceed in waves that initiate at the poles. Because genetic studies have not identified zygotic mutants that affect the early divisions and because transcription has not been detected before cycle 8, the early, pre-blastoderm embryo has been considered to rely entirely on maternal contributions and to be transcriptionally silent. Our studies identified several abnormal phenotypes in live engrailed (en) mutant embryos prior to cycle 8, as well as a small group of genes that are transcribed in embryos prior to cycle 7. Nuclei in en embryos divide asynchronously, an abnormality that was detected as early as nuclear cycle 2–3. Anti-En antibody detected nuclear En protein in embryos at cycle 2, and expression of an En:GFP fusion protein encoded in the paternal genome was also detected in cycle 2 nuclei. These findings demonstrate that the Drosophila embryo is functionally competent for gene expression prior to the onset of its rapid nuclear divisions and that the embryo requires functions that are expressed in the zygote in order to faithfully prosecute its early, pre-cellularization mitotic cycles.

Genetic studies identified many genes that are required during Drosophila oogenesis to endow the embryo with structures and components it will need to develop; they have also identified many genes that the embryo must express. However, measures of transcription have detected zygotic transcripts only after seven nuclear divisions, and many studies have concluded that zygotic mutants do not affect embryos prior to cellularization. The model that has emerged is that the earliest stages of embryogenesis rely solely on maternal stores and do not receive input from the zygotic genome. The fact that the embryo's nuclei divide rapidly with a cycling time of less than ten minutes has been interpreted to support this model, because it has been assumed that the nuclear cycle is too short for productive gene expression. Using sensitive measures of transcription and histological procedures that detect subtle differences, we found evidence for expression as early as nuclear cycle 2, and we identified a requirement for zygotic gene expression in embryos with just 2–4 nuclei. These findings challenge the idea that the Drosophila embryo is entirely pre-programmed and that its early development is under exclusive maternal control.

Insights into the evolution of Darwin's finches from comparative analysis of the Geospiza magnirostris genome sequence

Background

A classical example of repeated speciation coupled with ecological diversification is the evolution of 14 closely related species of Darwin’s (Galápagos) finches (Thraupidae, Passeriformes). Their adaptive radiation in the Galápagos archipelago took place in the last 2–3 million years and some of the molecular mechanisms that led to their diversification are now being elucidated. Here we report evolutionary analyses of genome of the large ground finch, Geospiza magnirostris.

Results

13,291 protein-coding genes were predicted from a 991.0 Mb G. magnirostris genome assembly. We then defined gene orthology relationships and constructed whole genome alignments between the G. magnirostris and other vertebrate genomes. We estimate that 15% of genomic sequence is functionally constrained between G. magnirostris and zebra finch. Genic evolutionary rate comparisons indicate that similar selective pressures acted along the G. magnirostris and zebra finch lineages suggesting that historical effective population size values have been similar in both lineages. 21 otherwise highly conserved genes were identified that each show evidence for positive selection on amino acid changes in the Darwin's finch lineage. Two of these genes (Igf2r and Pou1f1) have been implicated in beak morphology changes in Darwin’s finches. Five of 47 genes showing evidence of positive selection in early passerine evolution have cilia related functions, and may be examples of adaptively evolving reproductive proteins.

Conclusions

These results provide insights into past evolutionary processes that have shaped G. magnirostris genes and its genome, and provide the necessary foundation upon which to build population genomics resources that will shed light on more contemporaneous adaptive and non-adaptive processes that have contributed to the evolution of the Darwin’s finches.

Dual functions of TAF7L in adipocyte differentiation

The diverse transcriptional mechanisms governing cellular differentiation and development of mammalian tissue remains poorly understood. Here we report that TAF7L, a paralogue of TFIID subunit TAF7, is enriched in adipocytes and white fat tissue (WAT) in mouse. Depletion of TAF7L reduced adipocyte-specific gene expression, compromised adipocyte differentiation, and WAT development as well. Ectopic expression of TAF7L in myoblasts reprograms these muscle precursors into adipocytes upon induction. Genome-wide mRNA-seq expression profiling and ChIP-seq binding studies confirmed that TAF7L is required for activating adipocyte-specific genes via a dual mechanism wherein it interacts with PPARγ at enhancers and TBP/Pol II at core promoters. In vitro binding studies confirmed that TAF7L forms complexes with both TBP and PPARγ. These findings suggest that TAF7L plays an integral role in adipocyte gene expression by targeting enhancers as a cofactor for PPARγ and promoters as a component of the core transcriptional machinery.

DOI:http://dx.doi.org/10.7554/eLife.00170.001

The development of a single fertilized egg into a highly complex animal is determined by its genome, with a process called differential gene regulation exerting exquisite control over gene expression to ensure that various specialized cells are generated and that many types of tissue are produced. However, the mechanisms responsible for controlling gene expression and, therefore mammalian development, are poorly understood.

Researchers have developed a number of in vitro cell culture models to elucidate the details of differential gene regulation, and this approach has been used to characterize adipocytes—cells that store energy in the form of fat—for close to two decades. The formation of adipocytes, a process known as adipogenesis, has been extensively studied, but there remain major gaps in our knowledge: for example, the identities of many of the transcriptional regulators that are responsible for the differentiation of mesenchymal stem cells into adipocytes remain a mystery. This task is complicated by the fact that some of these regulators are involved in the differentiation of multiple cell lines, and that some of them also have multiple roles in the generation of a single cell type. In addition to being of fundamental interest, improving our knowledge of the properties and behavior of adipocytes is essential for tackling the increasing prevalence of obesity in the developed world.

Zhou et al. now report that TAF7L—a gene that was previously thought to be involved only in the production of sperm cells—has two roles in the differentiation of stem cells to form adipocytes. Using a combination of cellular, biochemical, genetic and genomic techniques, they show that TAF7L interacts with PPARγ, an important adipocyte transcriptional regulator at enhancer sites on the genome to increase the transcription of genes that are involved in adipogenesis. They also show that TAF7L interacts with a general transcription factor called TBP (short for TATA-binding protein) at promoter sequences, again to increase the expression of genes involved in adipogenesis. Moreover, they show that the expression of TAF7L in myoblasts—precursor cells that usually become muscle cells—can induce the formation of fat cells rather than muscle cells. Furthermore, mice lacking TAF7L are lean compared to their normal littermates. A clearer understanding of the underlying causes of fat cell formation could lead to the development of new approaches for the treatment of obesity and associated diseases.

DOI:http://dx.doi.org/10.7554/eLife.00170.002

Spatial promoter recognition signatures may enhance transcription factor specificity in yeast

The short length and high degeneracy of sites recognized by DNA-binding transcription factors limit the amount of information they can carry, and individual sites are rarely sufficient to mediate the regulation of specific targets. Computational analysis of microbial genomes has suggested that many factors function optimally when in a particular orientation and position with respect to their target promoters. To investigate this further, we developed and trained spatial models of binding site positioning and applied them to the genome of the yeast Saccharomyces cerevisiae. We found evidence of non-random organization of sites within promoters, differences in binding site density, or both for thirty-eight transcription factors. We show that these signatures allow transcription factors with substantial differences in binding site specificity to share similar promoter specificities. We illustrate how spatial information dictating the positioning and density of binding sites can in principle increase the information available to the organism for differentiating a transcription factor’s true targets, and we indicate how this information could potentially be leveraged for the same purpose in bioinformatic analyses.

A second-generation assembly of the Drosophila simulans genome provides new insights into patterns of lineage-specific divergence

We create a new assembly of the Drosophila simulans genome using 142 million paired short-read sequences and previously published data for strain w⁵⁰¹. Our assembly represents a higher-quality genomic sequence with greater coverage, fewer misassemblies, and, by several indexes, fewer sequence errors. Evolutionary analysis of this genome reference sequence reveals interesting patterns of lineage-specific divergence that are different from those previously reported. Specifically, we find that Drosophila melanogaster evolves faster than D. simulans at all annotated classes of sites, including putatively neutrally evolving sites found in minimal introns. While this may be partly explained by a higher mutation rate in D. melanogaster, we also find significant heterogeneity in rates of evolution across classes of sites, consistent with historical differences in the effective population size for the two species. Also contrary to previous findings, we find that the X chromosome is evolving significantly faster than autosomes for nonsynonymous and most noncoding DNA sites and significantly slower for synonymous sites. The absence of a X/A difference for putatively neutral sites and the robustness of the pattern to Gene Ontology and sex-biased expression suggest that partly recessive beneficial mutations may comprise a substantial fraction of noncoding DNA divergence observed between species. Our results have more general implications for the interpretation of evolutionary analyses of genomes of different quality.

Zelda binding in the early Drosophila melanogaster embryo marks regions subsequently activated at the maternal-to-zygotic transition

The earliest stages of development in most metazoans are driven by maternally deposited proteins and mRNAs, with widespread transcriptional activation of the zygotic genome occurring hours after fertilization, at a period known as the maternal-to-zygotic transition (MZT). In Drosophila, the MZT is preceded by the transcription of a small number of genes that initiate sex determination, patterning, and other early developmental processes; and the zinc-finger protein Zelda (ZLD) plays a key role in their transcriptional activation. To better understand the mechanisms of ZLD activation and the range of its targets, we used chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-Seq) to map regions bound by ZLD before (mitotic cycle 8), during (mitotic cycle 13), and after (late mitotic cycle 14) the MZT. Although only a handful of genes are transcribed prior to mitotic cycle 10, we identified thousands of regions bound by ZLD in cycle 8 embryos, most of which remain bound through mitotic cycle 14. As expected, early ZLD-bound regions include the promoters and enhancers of genes transcribed at this early stage. However, we also observed ZLD bound at cycle 8 to the promoters of roughly a thousand genes whose first transcription does not occur until the MZT and to virtually all of the thousands of known and presumed enhancers bound at cycle 14 by transcription factors that regulate patterned gene activation during the MZT. The association between early ZLD binding and MZT activity is so strong that ZLD binding alone can be used to identify active promoters and regulatory sequences with high specificity and selectivity. This strong early association of ZLD with regions not active until the MZT suggests that ZLD is not only required for the earliest wave of transcription but also plays a major role in activating the genome at the MZT.

The newly fertilized eggs of most animal species begin development with a series of rapid cell divisions. During this time of rapid DNA replication, there is little or no transcription of the embryo's genome, with the synthesis of new proteins being directed by a store of maternally deposited mRNAs. Several hours after fertilization, at a period known as the maternal-to-zygotic transition (MZT), transcription of the embryo's genome begins in earnest, but little is known about how this process is initiated. In this paper we investigate the role of a protein known as Zelda (or ZLD) at the MZT in the laboratory model insect Drosophila melanogaster. ZLD had been previously shown to control the activation of a small number of genes expressed prior to the MZT. Here, using an experimental technique (ChIP-Seq) that allowed us to visualize where on the genome a protein is bound, we show that, approximately an hour prior to the MZT, ZLD is bound to most of the genomic regions active at the MZT. This suggests that ZLD may act as a kind of an “on switch” for the zygotic genome, poising regions where it binds for activation at the MZT, and this raises the possibility that similar master regulators of the MZT exist in other species.

Control of embryonic stem cell lineage commitment by core promoter factor, TAF3

Deciphering the molecular basis of pluripotency is fundamental to our understanding of development and embryonic stem cell function. Here, we report that TAF3, a TBP-associated core promoter factor, is highly enriched in ES cells. In this context, TAF3 is required for endoderm lineage differentiation and prevents premature specification of neuroectoderm and mesoderm. In addition to its role in the core promoter recognition complex TFIID, genome-wide binding studies reveal that TAF3 localizes to a subset of chromosomal regions bound by CTCF/cohesin that are selectively associated with genes upregulated by TAF3. Notably, CTCF directly recruits TAF3 to promoter distal sites and TAF3-dependent DNA looping is observed between the promoter distal sites and core promoters occupied by TAF3/CTCF/cohesin. Together, our findings support a new role of TAF3 in mediating long-range chromatin regulatory interactions that safeguard the finely-balanced transcriptional programs underlying pluripotency.

The role of chromatin accessibility in directing the widespread, overlapping patterns of Drosophila transcription factor binding

Background

In Drosophila embryos, many biochemically and functionally unrelated transcription factors bind quantitatively to highly overlapping sets of genomic regions, with much of the lowest levels of binding being incidental, non-functional interactions on DNA. The primary biochemical mechanisms that drive these genome-wide occupancy patterns have yet to be established.

Results

Here we use data resulting from the DNaseI digestion of isolated embryo nuclei to provide a biophysical measure of the degree to which proteins can access different regions of the genome. We show that the in vivo binding patterns of 21 developmental regulators are quantitatively correlated with DNA accessibility in chromatin. Furthermore, we find that levels of factor occupancy in vivo correlate much more with the degree of chromatin accessibility than with occupancy predicted from in vitro affinity measurements using purified protein and naked DNA. Within accessible regions, however, the intrinsic affinity of the factor for DNA does play a role in determining net occupancy, with even weak affinity recognition sites contributing. Finally, we show that programmed changes in chromatin accessibility between different developmental stages correlate with quantitative alterations in factor binding.

Conclusions

Based on these and other results, we propose a general mechanism to explain the widespread, overlapping DNA binding by animal transcription factors. In this view, transcription factors are expressed at sufficiently high concentrations in cells such that they can occupy their recognition sequences in highly accessible chromatin without the aid of physical cooperative interactions with other proteins, leading to highly overlapping, graded binding of unrelated factors.

Noncanonical compensation of zygotic X transcription in early Drosophila melanogaster development revealed through single-embryo RNA-seq

Mmany genes from the X chromosome are expressed at the same level in female and male embryos during early Drosophila development, prior to the establishment of MSL-mediated dosage compensation, suggesting the existence of a novel mechanism.

When Drosophila melanogaster embryos initiate zygotic transcription around mitotic cycle 10, the dose-sensitive expression of specialized genes on the X chromosome triggers a sex-determination cascade that, among other things, compensates for differences in sex chromosome dose by hypertranscribing the single X chromosome in males. However, there is an approximately 1 hour delay between the onset of zygotic transcription and the establishment of canonical dosage compensation near the end of mitotic cycle 14. During this time, zygotic transcription drives segmentation, cellularization, and other important developmental events. Since many of the genes involved in these processes are on the X chromosome, we wondered whether they are transcribed at higher levels in females and whether this might lead to sex-specific early embryonic patterning. To investigate this possibility, we developed methods to precisely stage, sex, and characterize the transcriptomes of individual embryos. We measured genome-wide mRNA abundance in male and female embryos at eight timepoints, spanning mitotic cycle 10 through late cycle 14, using polymorphisms between parental lines to distinguish maternal and zygotic transcription. We found limited sex-specific zygotic transcription, with a weak tendency for genes on the X to be expressed at higher levels in females. However, transcripts derived from the single X chromosome in males were more abundant that those derived from either X chromosome in females, demonstrating that there is widespread dosage compensation prior to the activation of the canonical MSL-mediated dosage compensation system. Crucially, this new system of early zygotic dosage compensation results in nearly identical transcript levels for key X-linked developmental regulators, including giant (gt), brinker (brk), buttonhead (btd), and short gastrulation (sog), in male and female embryos.

Variation in gene dose can have profound effects on animal development. Yet every generation, animals must cope with differences in sex chromosome numbers. Drosophila compensate for the difference in X chromosome dosage (two in females, one in males) with a mechanism that allows for more transcription of the single X chromosome in males. But this mechanism is not established until over an hour after the embryo begins transcription, during which time a number of important events in development occur such as cellularization and segmentation. Here we use an mRNA sequencing method to characterize gene expression in individual female and male embryos before the onset of the previously characterized dosage compensation system. While we find more transcripts from X chromosomal genes in females, we also find many genes with equal transcript levels in males and females. These results indicate that there is an alternate mechanism to compensate for dosage acting earlier in development, prior to the onset of the previously characterized dosage compensation system.