Microarray analysis of Drosophila development during metamorphosis

Global gene expression shift during the transition from early neural development to late neuronal differentiation in Drosophila melanogaster

Regulation of transcription is one of the mechanisms involved in animal development, directing changes in patterning and cell fate specification. Large temporal data series, based on microarrays across the life cycle of the fly Drosophila melanogaster, revealed the existence of groups of genes which expression increases or decreases temporally correlated during the life cycle. These groups of genes are enriched in different biological functions. Here, instead of searching for temporal coincidence in gene expression using the entire genome expression data, we searched for temporal coincidence in gene expression only within predefined catalogues of functionally related genes and investigated whether a catalogue's expression profile can be used to generate larger catalogues, enriched in genes necessary for the same function. We analyzed the expression profiles from genes already associated with early neurodevelopment and late neurodifferentiation, at embryonic stages 16 and 17 of Drosophila life cycle. We hypothesized that during this interval we would find global downregulation of genes important for early neuronal development together with global upregulation of genes necessary for the final differentiation of neurons. Our results were consistent with this hypothesis. We then investigated if the expression profile of gene catalogues representing particular processes of neural development matched the temporal sequence along which these processes occur. The profiles of genes involved in patterning, neurogenesis, axogenesis or synaptic transmission matched the prediction, with largest transcript values at the time when the corresponding biological process takes place in the embryo. Furthermore, we obtained catalogues enriched in genes involved in temporally matching functions by performing a genome-wide systematic search for genes with their highest expression levels at the corresponding embryonic intervals. These findings imply the use of gene expression data in combination with known biological information to predict the involvement of functionally uncharacterized genes in particular biological events.

Predicting performance and plasticity in the development of respiratory structures and metabolic systems

The scaling laws governing metabolism suggest that we can predict metabolic rates across taxonomic scales that span large differences in mass. Yet, scaling relationships can vary with development, body region, and environment. Within species, there is variation in metabolic rate that is independent of mass and which may be explained by genetic variation, the environment or their interaction (i.e., metabolic plasticity). Additionally, some structures, such as the insect tracheal respiratory system, change throughout development and in response to the environment to match the changing functional requirements of the organism. We discuss how study of the development of respiratory function meets multiple challenges set forth by the NSF Grand Challenges Workshop. Development of the structure and function of respiratory and metabolic systems (1) is inherently stable and yet can respond dynamically to change, (2) is plastic and exhibits sensitivity to environments, and (3) can be examined across multiple scales in time and space. Predicting respiratory performance and plasticity requires quantitative models that integrate information across scales of function from the expression of metabolic genes and mitochondrial biogenesis to the building of respiratory structures. We present insect models where data are available on the development of the tracheal respiratory system and of metabolic physiology and suggest what is needed to develop predictive models. Incorporating quantitative genetic data will enable mapping of genetic and genetic-by-environment variation onto phenotypes, which is necessary to understand the evolution of respiratory and metabolic systems and their ability to enable respiratory homeostasis as organisms walk the tightrope between stability and change.

DDGni: dynamic delay gene-network inference from high-temporal data using gapped local alignment

MOTIVATION

Inferring gene-regulatory networks is very crucial in decoding various complex mechanisms in biological systems. Synthesis of a fully functional transcriptional factor/protein from DNA involves series of reactions, leading to a delay in gene regulation. The complexity increases with the dynamic delay induced by other small molecules involved in gene regulation, and noisy cellular environment. The dynamic delay in gene regulation is quite evident in high-temporal live cell lineage-imaging data. Although a number of gene-network-inference methods are proposed, most of them ignore the associated dynamic time delay.

RESULTS

Here, we propose DDGni (dynamic delay gene-network inference), a novel gene-network-inference algorithm based on the gapped local alignment of gene-expression profiles. The local alignment can detect short-term gene regulations, that are usually overlooked by traditional correlation and mutual Information based methods. DDGni uses 'gaps' to handle the dynamic delay and non-uniform sampling frequency in high-temporal data, like live cell imaging data. Our algorithm is evaluated on synthetic and yeast cell cycle data, and Caenorhabditis elegans live cell imaging data against other prominent methods. The area under the curve of our method is significantly higher when compared to other methods on all three datasets.

AVAILABILITY

The program, datasets and supplementary files are available at http://www.jjwanglab.org/DDGni/.

The effects of genetic variation on gene expression dynamics during development

The development of a multicellular organism and physiological responses require massive coordinated changes in gene expression across several cell and tissue types. Polymorphic regions of the genome that influence gene expression levels have been identified by expression quantitative trait locus (eQTL) mapping in many species, including loci that have cell-dependent, tissue-dependent and age-dependent effects. However, there has been no comprehensive characterization of how polymorphisms influence the complex dynamic patterns of gene expression that occur during development and in physiology. Here we describe an efficient experimental design to infer gene expression dynamics from single expression profiles in different genotypes, and apply it to characterize the effect of local (cis) and distant (trans) genetic variation on gene expression at high temporal resolution throughout a 12-hour period of the development of Caenorhabditis elegans. Taking dynamic variation into account identifies >50% more cis-eQTLs, including more than 900 that alter the dynamics of expression during this period. Local sequence polymorphisms extensively affect the timing, rate, magnitude and shape of expression changes. Indeed, many local sequence variants both increase and decrease gene expression, depending on the time-point profiled. Expression dynamics during this 12-hour period are also influenced extensively in trans by distal loci. In particular, several trans loci influence genes with quite diverse dynamic expression patterns, but they do so primarily during a common time interval. Trans loci can therefore act as modifiers of expression during a particular period of development. This study provides the first characterization, to our knowledge, of the effect of local and distant genetic variation on the dynamics of gene expression throughout an extensive time period. Moreover, the approach developed here should facilitate the genetic dissection of other dynamic processes, including potentially development, physiology and disease progression in humans.

Metamorphosis revealed: time-lapse three-dimensional imaging inside a living chrysalis

Studies of model insects have greatly increased our understanding of animal development. Yet, they are limited in scope to this small pool of model species: a small number of representatives for a hyperdiverse group with highly varied developmental processes. One factor behind this narrow scope is the challenging nature of traditional methods of study, such as histology and dissection, which can preclude quantitative analysis and do not allow the development of a single individual to be followed. Here, we use high-resolution X-ray computed tomography (CT) to overcome these issues, and three-dimensionally image numerous lepidopteran pupae throughout their development. The resulting models are presented in the electronic supplementary material, as are figures and videos, documenting a single individual throughout development. They provide new insight and details of lepidopteran metamorphosis, and allow the measurement of tracheal and gut volume. Furthermore, this study demonstrates early and rapid development of the tracheae, which become visible in scans just 12 h after pupation. This suggests that there is less remodelling of the tracheal system than previously expected, and is methodologically important because the tracheal system is an often-understudied character system in development. In the future, this form of time-lapse CT-scanning could allow faster and more detailed developmental studies on a wider range of taxa than is presently possible.

Different pattern of Galleria mellonella jhbp gene expression in high five and Sf9 cells

Juvenile hormone binding protein (JHBP) is the key element of the system that transmits hormone signals to target tissues. Recently, we found that the core promoter of the jhbp gene is strongly under the control of the TATA box and the transcription start site. In this report, we have shown that the jhbp promoter contains distal regulatory elements whose functionality clearly depends on the particular cell environment and that the scope of research from one cell line is insufficient to generalize the conclusions of the analysis. Cf1/Usp (where Usp is ultraspiracle protein previously known as Cf1, chorion factor 1) elements suppressed transcription of the reporter gene in the High Five cell line but not in the Sf9 cell line. However, upstream from all three Cf1/Usp elements there is a DNA sequence, containing the Zeste element, which activates jhbp in both systems. We found that juvenile hormone strongly inhibited the activity of the jhbp promoter in the Sf9 cell line, whereas it did not have an effect in the High Five cell line. A second key hormone that controls insect development--20-hydroxyecdysone, was also found to suppress the transcription of jhbp. This is the first report describing how these two hormones affect jhbp gene expression in different cell lines.

Interpreting the regulatory genome: the genomics of transcription factor function in Drosophila melanogaster

Researchers have now had access to the fully sequenced Drosophila melanogaster genome for over a decade, and the sequenced genomes of 11 additional Drosophila species have been available for almost 5 years, with more species' genomes becoming available every year [Adams MD, Celniker SE, Holt RA, et al. The genome sequence of Drosophila melanogaster. Science 2000;287:2185-95; Clark AG, Eisen MB, Smith DR, et al. Evolution of genes and genomes on the Drosophila phylogeny. Nature 2007;450:203-18]. Although the best studied of the D. melanogaster transcription factors (TFs) were cloned before sequencing of the genome, the availability of sequence data promised to transform our understanding of TFs and gene regulatory networks. Sequenced genomes have allowed researchers to generate tools for high-throughput characterization of gene expression levels, genome-wide TF localization and analyses of evolutionary constraints on DNA elements across multiple species. With an estimated 700 DNA-binding proteins in the Drosophila genome, it will be many years before each potential sequence-specific TF is studied in detail, yet the last decade of functional genomics research has already impacted our view of gene regulatory networks and TF DNA recognition.

Genomic analysis of the ecdysone steroid signal at metamorphosis onset using ecdysoneless and EcRnullDrosophila melanogaster mutants

Steroid hormone gene regulation is often depicted as a linear transduction of the signal, from molecule release to the gene level, by activation of a receptor protein after being bound by its steroid ligand. Such an action would require that the hormone be present and bound to the receptor in order to have target gene response. Here, we present data that presents a novel perspective of hormone gene regulation, where the hormone molecule and its receptor have exclusive target gene regulation function, in addition to the traditional direct target genes. Our study is the first genome-wide analysis of conditional mutants simultaneously modeling the steroid and steroid receptor gene expression regulation. We have integrated classical genetic mutant experiments with functional genomics techniques in the Drosophila melanogaster model organism, where we interrogate the 20-hydroxyecdysone signaling response at the onset of metamorphosis. Our novel catalog of ecdysone target genes illustrates the separable transcriptional responses among the hormone, the pre-hormone receptor and the post-hormone receptor. We successfully detected traditional ecdysone target genes as common targets and also identified novel sets of target genes which where exclusive to each mutant condition. Around 12 % of the genome responds to the ecdysone hormone signal at the onset of metamorphosis and over half of these are independent of the receptor. In addition, a significant portion of receptor regulated genes are differentially regulated by the receptor, depending on its ligand state. Gene ontology enrichment analyses confirm known ecdysone regulated biological functions and also validate implicated pathways that have been indirectly associated with ecdysone signaling.

Adaptive evolution and the birth of CTCF binding sites in the Drosophila genome

Comparative ChIP-seq data reveal adaptive evolution of insulator protein CTCF binding in multiple Drosophila species.

Changes in the physical interaction between cis-regulatory DNA sequences and proteins drive the evolution of gene expression. However, it has proven difficult to accurately quantify evolutionary rates of such binding change or to estimate the relative effects of selection and drift in shaping the binding evolution. Here we examine the genome-wide binding of CTCF in four species of Drosophila separated by between ∼2.5 and 25 million years. CTCF is a highly conserved protein known to be associated with insulator sequences in the genomes of human and Drosophila. Although the binding preference for CTCF is highly conserved, we find that CTCF binding itself is highly evolutionarily dynamic and has adaptively evolved. Between species, binding divergence increased linearly with evolutionary distance, and CTCF binding profiles are diverging rapidly at the rate of 2.22% per million years (Myr). At least 89 new CTCF binding sites have originated in the Drosophila melanogaster genome since the most recent common ancestor with Drosophila simulans. Comparing these data to genome sequence data from 37 different strains of Drosophila melanogaster, we detected signatures of selection in both newly gained and evolutionarily conserved binding sites. Newly evolved CTCF binding sites show a significantly stronger signature for positive selection than older sites. Comparative gene expression profiling revealed that expression divergence of genes adjacent to CTCF binding site is significantly associated with the gain and loss of CTCF binding. Further, the birth of new genes is associated with the birth of new CTCF binding sites. Our data indicate that binding of Drosophila CTCF protein has evolved under natural selection, and CTCF binding evolution has shaped both the evolution of gene expression and genome evolution during the birth of new genes.

A large proportion of the diversity of living organisms results from differential regulation of gene transcription. Transcriptional regulation is thought to differ between species because of evolutionary changes in the physical interactions between regulatory DNA elements and DNA-binding proteins; these can generate variation in the spatial and temporal patterns of gene expression. The mechanisms by which these protein–DNA interactions evolve is therefore an important question in evolutionary biology. Does adaptive evolution play a role, or is the process dominated by neutral genetic drift? Insulator proteins are a special group of DNA-binding proteins—instead of directly serving to activate or repress genes, they can function to coordinate the interactions between other regulatory elements (such as enhancers and promoters). Additionally, insulator proteins can limit the spreading of chromatin condensation and help to demarcate the boundaries of regulatory domains in the genome. In spite of their critical role in genome regulation, little is known about the evolution of interactions between insulator proteins and DNA. Here, we use ChIP-seq to examine the distribution of binding sites for CTCF, a highly conserved insulator protein, in four closely related Drosophila species. We find that genome-wide binding profiles of CTCF are highly dynamic across evolutionary time, with frequent births of new CTCF-DNA interactions, and we demonstrate that this evolutionary process is driven by natural selection. By comparing these with RNA-seq data, we find that gain or loss of CTCF binding impacts the expression levels of nearby genes and correlates with structural evolution of the genome. Together these results suggest a potential mechanism of regulatory re-wiring through adaptive evolution of CTCF binding.

Combinatorial control of temporal gene expression in the Drosophila wing by enhancers and core promoters

Background

The transformation of a developing epithelium into an adult structure is a complex process, which often involves coordinated changes in cell proliferation, metabolism, adhesion, and shape. To identify genetic mechanisms that control epithelial differentiation, we analyzed the temporal patterns of gene expression during metamorphosis of the Drosophila wing.

Results

We found that a striking number of genes, approximately 50% of the Drosophila transcriptome, exhibited changes in expression during a time course of wing development. While cis-acting enhancer sequences clearly correlated with these changes, a stronger correlation was discovered between core-promoter types and the dynamic patterns of gene expression within this differentiating tissue. In support of the hypothesis that core-promoter type influences the dynamics of expression, expression levels of several TATA-box binding protein associated factors (TAFs) and other core promoter-associated components changed during this developmental time course, and a testes-specific TAF (tTAF) played a critical role in timing cellular differentiation within the wing.

Conclusions

Our results suggest that the combinatorial control of gene expression via cis-acting enhancer sequences and core-promoter types, determine the complex changes in gene expression that drive morphogenesis and terminal differentiation of the Drosophila wing epithelium.

The ribonuclease Dis3 is an essential regulator of the developmental transcriptome

Background

Dis3 is ribonuclease that acts directly in the processing, turnover, and surveillance of a large number of distinct RNA species. Evolutionarily conserved from eubacteria to eukaryotes and a crucial component of the RNA processing exosome, Dis3 has been shown to be essential in yeast and fly S2 cells. However, it is not known whether Dis3 has essential functions in a metazoan. This study inquires whether Dis3 is required for Drosophila development and viability and how Dis3 regulates the transcriptome in the developing fly.

Results

Using transgenic flies, we show that Dis3 knock down (Dis3KD) retards growth, induces melanotic tumor formation, and ultimately results in 2^nd instar larval lethality. In order to determine whether Dis3KD fly phenotypes were a consequence of disrupting developmentally regulated RNA turnover, we performed RNA deep sequencing analysis on total RNA isolated from developmentally staged animals. Bioinformatic analysis of transcripts from Dis3KD flies reveals substantial transcriptomic changes, most notably down-regulation in early expressed RNAs. Finally, gene ontology analysis of this early stage shows that Dis3 regulates transcripts related to extracellular structure and remodelling, neurogenesis, and nucleotide metabolism.

Conclusions

We conclude that Dis3 is essential for early Drosophila melanogaster development and has specific and important stage-specific roles in regulating RNA metabolism. In showing for the first time that Dis3 is required for the development of a multicellular organism, our work provides mechanistic insight into how Dis3—either independent of or associated with the RNA processing exosome—participates in cell type-specific RNA turnover in metazoan development.

Identification and developmental profiling of conserved and novel microRNAs in Manduca sexta

MicroRNAs (miRNAs) are a group of small RNAs involved in translation inhibition or mRNA degradation. Due to its large size, Manduca sexta has long been used as a model to study insect physiology and biochemistry. While transcriptome studies have greatly enriched our knowledge on M. sexta structural genes, little is known about posttranscriptional regulation by miRNAs in this lepidopteran species. We constructed four small RNA libraries from embryos, 4th instar feeding larvae, pupae, and adults, obtained 21 million reads of 18-31 nucleotides by Illumina sequencing, and found 163 conserved and 13 novel miRNAs. By searching the M. sexta genome assembly, we identified precursors of 82 conserved miRNAs, 76 of which had mapped reads in one or more of these libraries. After normalization, we compared numbers of miRNA and miRNA-star reads in these libraries and observed abundance changes during development. Interestingly, mse-miR-281-star, mse-miR-31-star, mse-miR-965-star, mse-miR-9a-star, mse-miR-9b-star, mse-miR-2a-star, mse-miR-92b-star and mse-miR-279c-star are either more abundant or maintained at similar levels compared to respective mature miRNA strand. Expression profiling of the first set of miRNAs provided insights to their possible involvement in developmental regulation. This study will aid in the annotation of miRNA genes in the genome.

Transcriptome profiling of the demosponge Amphimedon queenslandica reveals genome-wide events that accompany major life cycle transitions

Background

The biphasic life cycle with pelagic larva and benthic adult stages is widely observed in the animal kingdom, including the Porifera (sponges), which are the earliest branching metazoans. The demosponge, Amphimedon queenslandica, undergoes metamorphosis from a free-swimming larva into a sessile adult that bears no morphological resemblance to other animals. While the genome of A. queenslandica contains an extensive repertoire of genes very similar to that of complex bilaterians, it is as yet unclear how this is drawn upon to coordinate changing morphological features and ecological demands throughout the sponge life cycle.

Results

To identify genome-wide events that accompany the pelagobenthic transition in A. queenslandica, we compared global gene expression profiles at four key developmental stages by sequencing the poly(A) transcriptome using SOLiD technology. Large-scale changes in transcription were observed as sponge larvae settled on the benthos and began metamorphosis. Although previous systematics suggest that the only clear homology between Porifera and other animals is in the embryonic and larval stages, we observed extensive use of genes involved in metazoan-associated cellular processes throughout the sponge life cycle. Sponge-specific transcripts are not over-represented in the morphologically distinct adult; rather, many genes that encode typical metazoan features, such as cell adhesion and immunity, are upregulated. Our analysis further revealed gene families with candidate roles in competence, settlement, and metamorphosis in the sponge, including transcription factors, G-protein coupled receptors and other signaling molecules.

Conclusions

This first genome-wide study of the developmental transcriptome in an early branching metazoan highlights major transcriptional events that accompany the pelagobenthic transition and point to a network of regulatory mechanisms that coordinate changes in morphology with shifting environmental demands. Metazoan developmental and structural gene orthologs are well-integrated into the expression profiles at every stage of sponge development, including the adult. The utilization of genes involved in metazoan-associated processes throughout sponge development emphasizes the potential of the genome of the last common ancestor of animals to generate phenotypic complexity.

Transcription profiling associated with life cycle of Anopheles gambiae

Complex biological events occur during the developmental process of the mosquito Anopheles gambiae (Giles). Using cDNA expression microarrays, the expression patterns of 13,440 clones representing 8,664 unique transcripts were revealed from six different developmental stages: early larvae (late third instar/early fourth instar), late larvae (late fourth instar), early pupae (< 30 min after pupation), late pupae (after tanning), and adult female and male mosquitoes (24 h postemergence). After microarray analysis, 560 unique transcripts were identified to show at least a fourfold up- or down-regulation in at least one developmental stage. Based on the expression patterns, these gene products were clustered into 13 groups. In total, eight genes were analyzed by quantitative real-time polymerase chain reaction to validate microarray results. Among 560 unique transcripts, 446 contigs were assigned to respective genes from the An. gambiae genome. The expression patterns and annotations of the genes in the 13 groups are discussed in the context of development including metabolism, transport, protein synthesis and degradation, cellular processes, cellular communication, intra- or extra-cellular architecture maintenance, response to stress or immune-related defense, and spermatogenesis.

Coordinating growth and maturation - insights from Drosophila

Adult body size in higher animals is dependent on the amount of growth that occurs during the juvenile stage. The duration of juvenile development, therefore, must be flexible and responsive to environmental conditions. When immature animals experience environmental stresses such as malnutrition or disease, maturation can be delayed until conditions improve and normal growth can resume. In contrast, when animals are raised under ideal conditions that promote rapid growth, internal checkpoints ensure that maturation does not occur until juvenile development is complete. Although the mechanisms that regulate growth and gate the onset of maturation have been investigated for decades, the emerging links between childhood obesity, early onset puberty, and adult metabolic disease have placed a new emphasis on this field. Remarkably, genetic studies in the fruit fly Drosophila melanogaster have shown that the central regulatory pathways that control growth and the timing of sexual maturation are conserved through evolution, and suggest that this aspect of animal life history is regulated by a common genetic architecture. This review focuses on these conserved mechanisms and highlights recent studies that explore how Drosophila coordinates developmental growth with environmental conditions.

Global gene expression analysis of murine limb development

Detailed information about stage-specific changes in gene expression is crucial for understanding the gene regulatory networks underlying development and the various signal transduction pathways contributing to morphogenesis. Here we describe the global gene expression dynamics during early murine limb development, when cartilage, tendons, muscle, joints, vasculature and nerves are specified and the musculoskeletal system of limbs is established. We used whole-genome microarrays to identify genes with differential expression at 5 stages of limb development (E9.5 to 13.5), during fore- and hind-limb patterning. We found that the onset of limb formation is characterized by an up-regulation of transcription factors, which is followed by a massive activation of genes during E10.5 and E11.5 which levels off at later time points. Among the 3520 genes identified as significantly up-regulated in the limb, we find ∼30% to be novel, dramatically expanding the repertoire of candidate genes likely to function in the limb. Hierarchical and stage-specific clustering identified expression profiles that are likely to correlate with functional programs during limb development and further characterization of these transcripts will provide new insights into specific tissue patterning processes. Here, we provide for the first time a comprehensive analysis of developmentally regulated genes during murine limb development, and provide some novel insights into the expression dynamics governing limb morphogenesis.

Generalized shrinkage F-like statistics for testing an interaction term in gene expression analysis in the presence of heteroscedasticity

Background

Many analyses of gene expression data involve hypothesis tests of an interaction term between two fixed effects, typically tested using a residual variance. In expression studies, the issue of variance heteroscedasticity has received much attention, and previous work has focused on either between-gene or within-gene heteroscedasticity. However, in a single experiment, heteroscedasticity may exist both within and between genes. Here we develop flexible shrinkage error estimators considering both between-gene and within-gene heteroscedasticity and use them to construct F-like test statistics for testing interactions, with cutoff values obtained by permutation. These permutation tests are complicated, and several permutation tests are investigated here.

Results

Our proposed test statistics are compared with other existing shrinkage-type test statistics through extensive simulation studies and a real data example. The results show that the choice of permutation procedures has dramatically more influence on detection power than the choice of F or F-like test statistics. When both types of gene heteroscedasticity exist, our proposed test statistics can control preselected type-I errors and are more powerful. Raw data permutation is not valid in this setting. Whether unrestricted or restricted residual permutation should be used depends on the specific type of test statistic.

Conclusions

The F-like test statistic that uses the proposed flexible shrinkage error estimator considering both types of gene heteroscedasticity and unrestricted residual permutation can provide a statistically valid and powerful test. Therefore, we recommended that it should always applied in the analysis of real gene expression data analysis to test an interaction term.

Quantitative trait locus analysis of stage-specific inbreeding depression in the Pacific oyster Crassostrea gigas

Inbreeding depression and genetic load have been widely observed, but their genetic basis and effects on fitness during the life cycle remain poorly understood, especially for marine animals with high fecundity and high, early mortality (type-III survivorship). A high load of recessive mutations was previously inferred for the Pacific oyster Crassostrea gigas, from massive distortions of zygotic, marker segregation ratios in F(2) families. However, the number, genomic location, and stage-specific onset of mutations affecting viability have not been thoroughly investigated. Here, we again report massive distortions of microsatellite-marker segregation ratios in two F(2) hybrid families, but we now locate the causative deleterious mutations, using a quantitative trait locus (QTL) interval-mapping model, and we characterize their mode of gene action. We find 14-15 viability QTL (vQTL) in the two families. Genotypic frequencies at vQTL generally suggest selection against recessive or partially recessive alleles, supporting the dominance theory of inbreeding depression. No epistasis was detected among vQTL, so unlinked vQTL presumably have independent effects on survival. For the first time, we track segregation ratios of vQTL-linked markers through the life cycle, to determine their stage-specific expression. Almost all vQTL are absent in the earliest life stages examined, confirming zygotic viability selection; vQTL are predominantly expressed before the juvenile stage (90%), mostly at metamorphosis (50%). We estimate that, altogether, selection on vQTL caused 96% mortality in these families, accounting for nearly all of the actual mortality. Thus, genetic load causes substantial mortality in inbred Pacific oysters, particularly during metamorphosis, a critical developmental transition warranting further investigation.

Steroidal regulation of hydrolyzing activity of the dietary carbohydrates in the silkworm, Bombyx mori

Blood sugar is an essential energy source for growth and development and is maintained at a constant level through precise regulation of formation and utilization. Sugars are produced from dietary carbohydrates by enzymatic hydrolysis in the digestive tract, which are under the homeostatic control of paracrine and prandial mechanisms in mammals. Here, we show that dietary carbohydrates hydrolyzing activity of the digestive tract is developmentally regulated by the steroid hormone ecdysone in the silkworm, Bombyx mori. The dietary carbohydrates hydrolyzing activity remained high throughout the last larval period and then decreased to negligible levels until the pupal period. However, dietary carbohydrates digestive activities were constitutively high when the steroidogenic organ, prothoracic glands were ablated. The prothoracic glands produced and released a large amount of ecdysone at the end of the larval period, suggesting that ecdysone is responsible for the decrease in dietary carbohydrates hydrolyzing activity. In fact, ecdysone decreased the activity to negligible levels in silkworms lacking the prothoracic glands. The present results indicate that the dietary carbohydrates hydrolyzing activity is regulated by ecdysone and that an increase in ecdysone titer decreases that activity at the end of the larval period, suggesting that ecdysone is essential for metabolic coordination during development.

Development and application of a modified dynamic time warping algorithm (DTW-S) to analyses of primate brain expression time series

Background

Comparing biological time series data across different conditions, or different specimens, is a common but still challenging task. Algorithms aligning two time series represent a valuable tool for such comparisons. While many powerful computation tools for time series alignment have been developed, they do not provide significance estimates for time shift measurements.

Results

Here, we present an extended version of the original DTW algorithm that allows us to determine the significance of time shift estimates in time series alignments, the DTW-Significance (DTW-S) algorithm. The DTW-S combines important properties of the original algorithm and other published time series alignment tools: DTW-S calculates the optimal alignment for each time point of each gene, it uses interpolated time points for time shift estimation, and it does not require alignment of the time-series end points. As a new feature, we implement a simulation procedure based on parameters estimated from real time series data, on a series-by-series basis, allowing us to determine the false positive rate (FPR) and the significance of the estimated time shift values. We assess the performance of our method using simulation data and real expression time series from two published primate brain expression datasets. Our results show that this method can provide accurate and robust time shift estimates for each time point on a gene-by-gene basis. Using these estimates, we are able to uncover novel features of the biological processes underlying human brain development and maturation.

Conclusions

The DTW-S provides a convenient tool for calculating accurate and robust time shift estimates at each time point for each gene, based on time series data. The estimates can be used to uncover novel biological features of the system being studied. The DTW-S is freely available as an R package TimeShift at http://www.picb.ac.cn/Comparative/data.html.

Linking fold, function and phylogeny: a comparative genomics view on protein (domain) evolution

Domains are the building blocks of all globular proteins and present one of the most useful levels at which protein function can be understood. Through recombination and duplication of a limited set of domains, proteomes evolved and the collection of protein superfamilies in an organism formed. As such, the presence of a shared domain can be regarded as an indicator of similar function and evolutionary history, but it does not necessarily imply it since convergent evolution may give rise to similar gene functions as well as architectures.Through the wealth of sequences and annotation data brought about by genomics, evolutionary links can be sought for via homology relationships and comparative genomics, structural modeling and phylogenetics. The goal hereby is not only to predict the function of newly discovered proteins, but also to spell out their pathway of evolution and, possibly, identify their most likely origin. This can ultimately help to understand protein function and functional relationships of protein families. Additionally, through comparison with transcriptional data, evolutionary data can be linked to gene (and genome) activity and thus allow for the identification of common principles behind fast evolving proteins and relatively stable ones.In this review, we describe the basic principles of studying protein (domain) evolution and illustrate recent developments in molecular evolution and give valuable new insights in the field of comparative genomics. As an example, we include here molecular models of the multiple PDZ domain protein MUPP-1 and present a simple comparative genomic view on its structural course of evolution.

High-dimensional variance partitioning reveals the modular genetic basis of adaptive divergence in gene expression during reproductive character displacement

Although adaptive change is usually associated with complex changes in phenotype, few genetic investigations have been conducted on adaptations that involve sets of high-dimensional traits. Microarrays have supplied high-dimensional descriptions of gene expression, and phenotypic change resulting from adaptation often results in large-scale changes in gene expression. We demonstrate how genetic analysis of large-scale changes in gene expression generated during adaptation can be accomplished by determining high-dimensional variance partitioning within classical genetic experimental designs. A microarray experiment conducted on a panel of recombinant inbred lines (RILs) generated from two populations of Drosophila serrata that have diverged in response to natural selection, revealed genetic divergence in 10.6% of 3762 gene products examined. Over 97% of the genetic divergence in transcript abundance was explained by only 12 genetic modules. The two most important modules, explaining 50% of the genetic variance in transcript abundance, were genetically correlated with the morphological traits that are known to be under selection. The expression of three candidate genes from these two important genetic modules was assessed in an independent experiment using qRT-PCR on 430 individuals from the panel of RILs, and confirmed the genetic association between transcript abundance and morphological traits under selection.

Developmental transcriptome of Aplysia californica

Genome-wide transcriptional changes in development provide important insight into mechanisms underlying growth, differentiation, and patterning. However, such large-scale developmental studies have been limited to a few representatives of Ecdysozoans and Chordates. Here, we characterize transcriptomes of embryonic, larval, and metamorphic development in the marine mollusc Aplysia californica and reveal novel molecular components associated with life history transitions. Specifically, we identify more than 20 signal peptides, putative hormones, and transcription factors in association with early development and metamorphic stages-many of which seem to be evolutionarily conserved elements of signal transduction pathways. We also characterize genes related to biomineralization-a critical process of molluscan development. In summary, our experiment provides the first large-scale survey of gene expression in mollusc development, and complements previous studies on the regulatory mechanisms underlying body plan patterning and the formation of larval and juvenile structures. This study serves as a resource for further functional annotation of transcripts and genes in Aplysia, specifically and molluscs in general. A comparison of the Aplysia developmental transcriptome with similar studies in the zebra fish Danio rerio, the fruit fly Drosophila melanogaster, the nematode Caenorhabditis elegans, and other studies on molluscs suggests an overall highly divergent pattern of gene regulatory mechanisms that are likely a consequence of the different developmental modes of these organisms.

Transcriptomics and proteomics in human African trypanosomiasis: current status and perspectives

Human African trypanosomiasis, or sleeping sickness, is a neglected vector-borne parasitic disease caused by protozoa of the species Trypanosoma brucei sensu lato. Within this complex species, T. b. gambiense is responsible for the chronic form of sleeping sickness in Western and Central Africa, whereas T. b. rhodesiense causes the acute form of the disease in East Africa. Presently, 1.5 million disability-adjusted life years (DALYs) per year are lost due to sleeping sickness. In addition, on the basis of the mortality, the disease is ranked ninth out of 25 human infectious and parasitic diseases in Africa. Diagnosis is complex and needs the intervention of a specialized skilled staff; treatment is difficult and expensive and has potentially life-threatening side effects. The use of transcriptomic and proteomic technologies, currently in rapid development and increasing in sensitivity and discriminating power, is already generating a large panel of promising results. The objective of these technologies is to significantly increase our knowledge of the molecular mechanisms governing the parasite establishment in its vector, the development cycle of the parasite during the parasite's intra-vector life, its interactions with the fly and the other microbial inhabitants of the gut, and finally human host-trypanosome interactions. Such fundamental investigations are expected to provide opportunities to identify key molecular events that would constitute accurate targets for further development of tools dedicated to field work for early, sensitive, and stage-discriminant diagnosis, epidemiology, new chemotherapy, and potentially vaccine development, all of which will contribute to fighting the disease. The present review highlights the contributions of the transcriptomic and proteomic analyses developed thus far in order to identify potential targets (genes or proteins) and biological pathways that may constitute a critical step in the identification of new targets for the development of new tools for diagnostic and therapeutic purposes.

Changed gene expression for candidate ageing genes in long-lived Bicyclus anynana butterflies

Candidate genes for the regulation of lifespan have emerged from studies that use mutants and genetically manipulated model organisms. However, it is rarely addressed whether these genes contribute to lifespan variation in populations of these species that capture natural standing genetic variation. Here, we explore expression variation in three candidate ageing genes, Indy, sod2, and catalase, in Bicyclus anynana, a butterfly with well understood ecology. We used lines established from natural populations and artificially selected for increased adult starvation resistance. They show a considerable increase in adult lifespan under both starvation and optimal food conditions. We measured adult butterflies of various ages, under a range of optimal and starvation diets, from two selected populations and one unselected control population. In all lines, Indy and catalase are up-regulated in response to starvation while this is not evident for sod2. Under starvation, Indy and catalase are up-regulated in, while this is not evident for sod2. Under optimal food conditions, Indy is down-regulated at a later age, with Indy expression showing relatively high inter-individual variation. We find differences between the selected lines and the unselected line. Under starvation conditions, expression is higher for catalase in one, and for sod2 in both selected lines. Importantly, sod2 expression is also higher in the selected populations under optimal food conditions. We conclude that sod2, but not Indy, is involved in the response to artificial selection for increased starvation resistance. The role of catalase is less clear because of the differences between the two selected lines. Moreover, sod2 appears to be a candidate gene that underpins the genetic correlation between starvation resistance and longevity. Our study indicates that some, but not all, genes identified through mutant screens in other organisms may underpin standing genetic variation for ageing-related traits in stocks of Bicyclus butterflies established from natural populations. Clearly, this needs to be investigated in other organisms as well, especially in the organisms to which mutants screens were applied. This information will narrow down the list of genes that underpin variation in lifespan and ageing in extant populations of organisms, and which may serve as candidate genes in humans.

Rift valley fever virus infection of human cells and insect hosts is promoted by protein kinase C epsilon

As an arthropod-borne human pathogen, Rift Valley fever virus (RVFV) cycles between an insect vector and mammalian hosts. Little is known about the cellular requirements for infection in either host. Here we developed a tissue culture model for RVFV infection of human and insect cells that is amenable to high-throughput screening. Using this approach we screened a library of 1280 small molecules with pharmacologically defined activities and identified 59 drugs that inhibited RVFV infection with 15 inhibiting RVFV replication in both human and insect cells. Amongst the 15 inhibitors that blocked infection in both hosts was a subset that inhibits protein kinase C. Further studies found that infection is dependent upon the novel protein kinase C isozyme epsilon (PKCε) in both human and insect cells as well as in adult flies. Altogether, these data show that inhibition of cellular factors required for early steps in the infection cycle including PKCε can block RVFV infection, and may represent a starting point for the development of anti-RVFV therapeutics.

Temporal waves of coherent gene expression during Drosophila embryogenesis

MOTIVATION

Animal development depends on localized patterns of gene expression. Whole-genome methods permit the global identification of differential expression patterns. However, most gene-expression-clustering methods focus on the analysis of entire expression profiles, rather than temporal segments or time windows.

RESULTS

In the current study, local clustering of temporal time windows was applied to developing embryos of the fruitfly, Drosophila melanogaster. Large-scale developmental events, involving temporal activation of hundreds of genes, were identified as discrete gene clusters. The time-duration analysis revealed six temporal waves of coherent gene expression during Drosophila embryogenesis. The most powerful expression waves preceded major morphogenetic movements, such as germ band elongation and dorsal closure. These waves of gene expression coincide with the inhibition of maternal transcripts during early development, the specification of ectoderm, differentiation of the nervous system, differentiation of the digestive tract, deposition of the larval cuticle and the reorganization of the cytoskeleton during global morphogenetic events. We discuss the implications of these findings with respect to the gene regulatory networks governing Drosophila development.

AVAILABILITY

Data and software are available from the UC Berkeley web resource http://flydev.berkeley.edu/cgi-bin/GTEM/dmap_dm-ag/index_dmap.htm

DroID 2011: a comprehensive, integrated resource for protein, transcription factor, RNA and gene interactions for Drosophila

DroID (http://droidb.org/), the Drosophila Interactions Database, is a comprehensive public resource for Drosophila gene and protein interactions. DroID contains genetic interactions and experimentally detected protein–protein interactions curated from the literature and from external databases, and predicted protein interactions based on experiments in other species. Protein interactions are annotated with experimental details and periodically updated confidence scores. Data in DroID is accessible through user-friendly, intuitive interfaces that allow simple or advanced searches and graphical visualization of interaction networks. DroID has been expanded to include interaction types that enable more complete analyses of the genetic networks that underlie biological processes. In addition to protein–protein and genetic interactions, the database now includes transcription factor–gene and regulatory RNA–gene interactions. In addition, DroID now has more gene expression data that can be used to search and filter interaction networks. Orthologous gene mappings of Drosophila genes to other organisms are also available to facilitate finding interactions based on gene names and identifiers for a number of common model organisms and humans. Improvements have been made to the web and graphical interfaces to help biologists gain a comprehensive view of the interaction networks relevant to the genes and systems that they study.

Developmental regulation of glycolysis by 20-hydroxyecdysone and juvenile hormone in fat body tissues of the silkworm, Bombyx mori

20-Hydroxyecdysone (20E) and juvenile hormone (JH) control a variety of physiological events during insect development and metamorphosis. To understand how 20E and JH developmentally regulate energy metabolism in insects, we performed a genome-wide microarray analysis of fat body tissues isolated from the silkworm, Bombyx mori. Many genes involved in energy metabolism, including genes in the glycolytic pathway, were down-regulated during molting and pupation, when 20E levels are high. Notably, 20E treatment exhibited inhibitory effects on key glycolytic enzyme mRNA levels and activities, and RNA interference of the 20E receptor EcR-USP had the opposite effects to 20E treatment. Meanwhile, JH treatment stimulated both mRNA levels and activities of the key glycolytic enzymes, presumably via antagonizing the 20E action. Taken together, we conclude that 20E acts as a general blocker for glycolysis in the Bombyx fat body during molting and pupation, whereas the physiological role of JH is contrast with 20E during molting.

The nature of protein domain evolution: shaping the interaction network

The proteomes that make up the collection of proteins in contemporary organisms evolved through recombination and duplication of a limited set of domains. These protein domains are essentially the main components of globular proteins and are the most principal level at which protein function and protein interactions can be understood. An important aspect of domain evolution is their atomic structure and biochemical function, which are both specified by the information in the amino acid sequence. Changes in this information may bring about new folds, functions and protein architectures. With the present and still increasing wealth of sequences and annotation data brought about by genomics, new evolutionary relationships are constantly being revealed, unknown structures modeled and phylogenies inferred. Such investigations not only help predict the function of newly discovered proteins, but also assist in mapping unforeseen pathways of evolution and reveal crucial, co-evolving inter- and intra-molecular interactions. In turn this will help us describe how protein domains shaped cellular interaction networks and the dynamics with which they are regulated in the cell. Additionally, these studies can be used for the design of new and optimized protein domains for therapy. In this review, we aim to describe the basic concepts of protein domain evolution and illustrate recent developments in molecular evolution that have provided valuable new insights in the field of comparative genomics and protein interaction networks.

A mathematical model for suppression subtractive hybridization

Suppression subtractive hybridization (SSH) is frequently used to unearth differentially expressed genes on a whole-genome scale. Its versatility is based on combining cDNA library subtraction and normalization, which allows the isolation of sequences of varying degrees of abundance and differential expression. SSH is a complex process with many adjustable parameters that affect the outcome of gene isolation.We present a mathematical model of SSH based on DNA hybridization kinetics for assessing the effect of various parameters to facilitate its optimization. We derive an equation for the probability that a particular differentially expressed species is successfully isolated and use this to quantify the effect of the following parameters related to the cDNA sample: (a) mRNA abundance; (b) partial sequence complementarity to other species; and (3) degree of differential expression. We also evaluate the effect of parameters related to the process, including: (a) reaction times; and (b) extent of driver excess used in the two hybridization reactions. The optimum set of process parameters for successful isolation of differentially expressed species depends on transcript abundance. We show that the reaction conditions have a significant effect on the occurrence of false-positives and formulate strategies to isolate specific subsets of differentially expressed genes. We also quantify the effect of non-specific hybridization on the false-positive results and present strategies for spiking cDNA sequences to address this problem.

Bowman-Birk inhibitor affects pathways associated with energy metabolism in Drosophila melanogaster

Bowman-Birk inhibitor (BBI) is toxic when fed to certain insects, including the fruit fly, Drosophila melanogaster. Dietary BBI has been demonstrated to slow growth and increase insect mortality by inhibiting the digestive enzymes trypsin and chymotrypsin, resulting in a reduced supply of amino acids. In mammals, BBI influences cellular energy metabolism. Therefore, we tested the hypothesis that dietary BBI affects energy-associated pathways in the D. melanogaster midgut. Through microarray and metabolomic analyses, we show that dietary BBI affects energy utilization pathways in the midgut cells of D. melanogaster. In addition, ultrastructure studies indicate that microvilli are significantly shortened in BBI-fed larvae. These data provide further insights into the complex cellular response of insects to dietary protease inhibitors.

Gene expression atlas for human embryogenesis

Human embryogenesis is believed to involve an integrated set of complex yet coordinated development of different organs and tissues mediated by the changes in the spatiotemporal expression of many genes. Here, we report a genome-wide expression analysis during wk 4-9 of human embryogenesis, a critical period when most organs develop. About half of all human genes are expressed, and 18.6% of the expressed genes were significantly regulated during this important period. We further identified >5000 regulated genes, most of which previously were not known to be associated with animal development. Our study fills an important gap in mammalian developmental studies by identifying functional pathways involved in this critical but previously not studied period. Our study also revealed that the genes involved here are distinct from those during early embryogenesis, which include three groups of maternal genes. Furthermore, we discovered that genes in a given developmental process are regulated coordinately. This led us to develop an easily searchable database of this entire collection of gene expression profiles, allowing for the identification new genes important for a particular developmental process/pathway and deducing the potential function of a novel gene. The validity of the predictions from the database was demonstrated with two examples through spatiotemporal analyses of the two novel genes. Such a database should serve as a highly valuable resource for the molecular analysis of human development and pathogenesis.

Intensity-based analysis of dual-color gene expression data as an alternative to ratio-based analysis to enhance reproducibility

BACKGROUND

Ratio-based analysis is the current standard for the analysis of dual-color microarray data. Indeed, this method provides a powerful means to account for potential technical variations such as differences in background signal, spot size and spot concentration. However, current high density dual-color array platforms are of very high quality, and inter-array variance has become much less pronounced. We therefore raised the question whether it is feasible to use an intensity-based analysis rather than ratio-based analysis of dual-color microarray datasets. Furthermore, we compared performance of both ratio- and intensity-based analyses in terms of reproducibility and sensitivity for differential gene expression.

RESULTS

By analyzing three distinct and technically replicated datasets with either ratio- or intensity-based models, we determined that, when applied to the same dataset, intensity-based analysis of dual-color gene expression experiments yields 1) more reproducible results, and 2) is more sensitive in the detection of differentially expressed genes. These effects were most pronounced in experiments with large biological variation and complex hybridization designs. Furthermore, a power analysis revealed that for direct two-group comparisons above a certain sample size, ratio-based models have higher power, although the difference with intensity-based models is very small.

CONCLUSIONS

Intensity-based analysis of dual-color datasets results in more reproducible results and increased sensitivity in the detection of differential gene expression than the analysis of the same dataset with ratio-based analysis. Complex dual-color setups such as interwoven loop designs benefit most from ignoring the array factor. The applicability of our approach to array platforms other than dual-color needs to be further investigated.

Deep sequencing of small RNA libraries reveals dynamic regulation of conserved and novel microRNAs and microRNA-stars during silkworm development

BACKGROUND

In eukaryotes, microRNAs (miRNAs) have emerged as critical regulators of gene expression. The Silkworm (Bombyx mori L.) is one of the most suitable lepidopteran insects for studying the molecular aspects of metamorphosis because of its large size, availability of mutants and genome sequence. Besides, this insect also has been amply studied from a physiological and biochemical perspective. Deep sequencing of small RNAs isolated from different stages of silkworm is a powerful tool not only for measuring the changes in miRNA profile but also for discovering novel miRNAs.

RESULTS

We generated small RNA libraries from feeding larvae, spinning larvae, pupae and adults of B. mori and obtained approximately 2.5 million reads of 18-30 nt. Sequence analysis identified 14 novel and 101 conserved miRNAs. Most novel miRNAs are preferentially expressed in pupae, whereas more than 95% of the conserved miRNAs are dynamically regulated during different developmental stages. Remarkably, the miRNA-star (miR*) of four miRNAs are expressed at much higher levels than their corresponding miRNAs, and their expression profiles are distinct from their corresponding miRNA profiles during different developmental stages. Additionally, we detected two antisense miRNA loci (miR-263-S and miR-263-AS; miR-306-S and miR-306-AS) that are expressed in sense and antisense directions. Interestingly, miR-263 and miR-306 are preferentially and abundantly expressed in pupae and adults, respectively.

CONCLUSIONS

We identified 101 homologs of conserved miRNAs, 14 species-specific and two antisense miRNAs in the silkworm. Our results provided deeper insights into changes in conserved and novel miRNA and miRNA* accumulation during development.

Deep sequencing of small RNA libraries reveals dynamic regulation of conserved and novel microRNAs and microRNA-stars during silkworm development

BACKGROUND

In eukaryotes, microRNAs (miRNAs) have emerged as critical regulators of gene expression. The Silkworm (Bombyx mori L.) is one of the most suitable lepidopteran insects for studying the molecular aspects of metamorphosis because of its large size, availability of mutants and genome sequence. Besides, this insect also has been amply studied from a physiological and biochemical perspective. Deep sequencing of small RNAs isolated from different stages of silkworm is a powerful tool not only for measuring the changes in miRNA profile but also for discovering novel miRNAs.

RESULTS

We generated small RNA libraries from feeding larvae, spinning larvae, pupae and adults of B. mori and obtained approximately 2.5 million reads of 18-30 nt. Sequence analysis identified 14 novel and 101 conserved miRNAs. Most novel miRNAs are preferentially expressed in pupae, whereas more than 95% of the conserved miRNAs are dynamically regulated during different developmental stages. Remarkably, the miRNA-star (miR*) of four miRNAs are expressed at much higher levels than their corresponding miRNAs, and their expression profiles are distinct from their corresponding miRNA profiles during different developmental stages. Additionally, we detected two antisense miRNA loci (miR-263-S and miR-263-AS; miR-306-S and miR-306-AS) that are expressed in sense and antisense directions. Interestingly, miR-263 and miR-306 are preferentially and abundantly expressed in pupae and adults, respectively.

CONCLUSIONS

We identified 101 homologs of conserved miRNAs, 14 species-specific and two antisense miRNAs in the silkworm. Our results provided deeper insights into changes in conserved and novel miRNA and miRNA* accumulation during development.