MicroRNAs (miRNAs) are endogenous approximately 22 nt RNAs that can play important regulatory roles in animals and plants by targeting mRNAs for cleavage or translational repression. Although they escaped notice until relatively recently, miRNAs comprise one of the more abundant classes of gene regulatory molecules in multicellular organisms and likely influence the output of many protein-coding genes.

High-throughput sequencing assays such as RNA-Seq, ChIP-Seq or barcode counting provide quantitative readouts in the form of count data. To infer differential signal in such data correctly and with good statistical power, estimation of data variability throughout the dynamic range and a suitable error model are required. We propose a method based on the negative binomial distribution, with variance and mean linked by local regression and present an implementation, DESeq, as an R/Bioconductor package.

We have designed a system for targeted gene expression that allows the selective activation of any cloned gene in a wide variety of tissue- and cell-specific patterns. The gene encoding the yeast transcriptional activator GAL4 is inserted randomly into the Drosophila genome to drive GAL4 expression from one of a diverse array of genomic enhancers. It is then possible to introduce a gene containing GAL4 binding sites within its promoter, to activate it in those cells where GAL4 is expressed, and to observe the effect of this directed misexpression on development. We have used GAL4-directed transcription to expand the domain of embryonic expression of the homeobox protein even-skipped. We show that even-skipped represses wingless and transforms cells that would normally secrete naked cuticle into denticle secreting cells. The GAL4 system can thus be used to study regulatory interactions during embryonic development. In adults, targeted expression can be used to generate dominant phenotypes for use in genetic screens. We have directed expression of an activated form of the Dras2 protein, resulting in dominant eye and wing defects that can be used in screens to identify other members of the Dras2 signal transduction pathway.

Microorganisms that invade a vertebrate host are initially recognized by the innate immune system through germline-encoded pattern-recognition receptors (PRRs). Several classes of PRRs, including Toll-like receptors and cytoplasmic receptors, recognize distinct microbial components and directly activate immune cells. Exposure of immune cells to the ligands of these receptors activates intracellular signaling cascades that rapidly induce the expression of a variety of overlapping and unique genes involved in the inflammatory and immune responses. New insights into innate immunity are changing the way we think about pathogenesis and the treatment of infectious diseases, allergy, and autoimmunity.

RNA interference (RNAi) is the process of sequence-specific, post-transcriptional gene silencing in animals and plants, initiated by double-stranded RNA (dsRNA) that is homologous in sequence to the silenced gene. The mediators of sequence-specific messenger RNA degradation are 21- and 22-nucleotide small interfering RNAs (siRNAs) generated by ribonuclease III cleavage from longer dsRNAs. Here we show that 21-nucleotide siRNA duplexes specifically suppress expression of endogenous and heterologous genes in different mammalian cell lines, including human embryonic kidney (293) and HeLa cells. Therefore, 21-nucleotide siRNA duplexes provide a new tool for studying gene function in mammalian cells and may eventually be used as gene-specific therapeutics.

The newest version of MUMmer easily handles comparisons of large eukaryotic genomes at varying evolutionary distances, as demonstrated by applications to multiple genomes. Two new graphical viewing tools provide alternative ways to analyze genome alignments. The new system is the first version of MUMmer to be released as open-source software. This allows other developers to contribute to the code base and freely redistribute the code. The MUMmer sources are available at http://www.tigr.org/software/mummer.

Methods for alignment of protein sequences typically measure similarity by using a substitution matrix with scores for all possible exchanges of one amino acid with another. The most widely used matrices are based on the Dayhoff model of evolutionary rates. Using a different approach, we have derived substitution matrices from about 2000 blocks of aligned sequence segments characterizing more than 500 groups of related proteins. This led to marked improvements in alignments and in searches using queries from each of the groups.

The innate immune system in drosophila and mammals senses the invasion of microorganisms using the family of Toll receptors, stimulation of which initiates a range of host defense mechanisms. In drosophila antimicrobial responses rely on two signaling pathways: the Toll pathway and the IMD pathway. In mammals there are at least 10 members of the Toll-like receptor (TLR) family that recognize specific components conserved among microorganisms. Activation of the TLRs leads not only to the induction of inflammatory responses but also to the development of antigen-specific adaptive immunity. The TLR-induced inflammatory response is dependent on a common signaling pathway that is mediated by the adaptor molecule MyD88. However, there is evidence for additional pathways that mediate TLR ligand-specific biological responses.

RNA interference (RNAi) is the mechanism through which double-stranded RNAs silence cognate genes. In plants, this can occur at both the transcriptional and the post-transcriptional levels; however, in animals, only post-transcriptional RNAi has been reported to date. In both plants and animals, RNAi is characterized by the presence of RNAs of about 22 nucleotides in length that are homologous to the gene that is being suppressed. These 22-nucleotide sequences serve as guide sequences that instruct a multicomponent nuclease, RISC, to destroy specific messenger RNAs. Here we identify an enzyme, Dicer, which can produce putative guide RNAs. Dicer is a member of the RNase III family of nucleases that specifically cleave double-stranded RNAs, and is evolutionarily conserved in worms, flies, plants, fungi and mammals. The enzyme has a distinctive structure, which includes a helicase domain and dual RNase III motifs. Dicer also contains a region of homology to the RDE1/QDE2/ARGONAUTE family that has been genetically linked to RNAi.

Programmed cell death is a distinct genetic and biochemical pathway essential to metazoans. An intact death pathway is required for successful embryonic development and the maintenance of normal tissue homeostasis. Apoptosis has proven to be tightly interwoven with other essential cell pathways. The identification of critical control points in the cell death pathway has yielded fundamental insights for basic biology, as well as provided rational targets for new therapeutics.

Induction of the adaptive immune response depends on the expression of co-stimulatory molecules and cytokines by antigen-presenting cells. The mechanisms that control the initial induction of these signals upon infection are poorly understood. It has been proposed that their expression is controlled by the non-clonal, or innate, component of immunity that preceded in evolution the development of an adaptive immune system in vertebrates. We report here the cloning and characterization of a human homologue of the Drosophila toll protein (Toll) which has been shown to induce the innate immune response in adult Drosophila. Like Drosophila Toll, human Toll is a type I transmembrane protein with an extracellular domain consisting of a leucine-rich repeat (LRR) domain, and a cytoplasmic domain homologous to the cytoplasmic domain of the human interleukin (IL)-1 receptor. Both Drosophila Toll and the IL-1 receptor are known to signal through the NF-kappaB pathway. We show that a constitutively active mutant of human Toll transfected into human cell lines can induce the activation of NF-kappaB and the expression of NF-kappaB-controlled genes for the inflammatory cytokines IL-1, IL-6 and IL-8, as well as the expression of the co-stimulatory molecule B7.1, which is required for the activation of naive T cells.

Polycomb group (PcG) proteins play important roles in maintaining the silent state of HOX genes. Recent studies have implicated histone methylation in long-term gene silencing. However, a connection between PcG-mediated gene silencing and histone methylation has not been established. Here we report the purification and characterization of an EED-EZH2 complex, the human counterpart of the Drosophila ESC-E(Z) complex. We demonstrate that the complex specifically methylates nucleosomal histone H3 at lysine 27 (H3-K27). Using chromatin immunoprecipitation assays, we show that H3-K27 methylation colocalizes with, and is dependent on, E(Z) binding at an Ultrabithorax (Ubx) Polycomb response element (PRE), and that this methylation correlates with Ubx repression. Methylation on H3-K27 facilitates binding of Polycomb (PC), a component of the PRC1 complex, to histone H3 amino-terminal tail. Thus, these studies establish a link between histone methylation and PcG-mediated gene silencing.

Proteins often differ in amino-acid sequence across species. This difference has evolved by the accumulation of neutral mutations by random drift, the fixation of adaptive mutations by selection, or a mixture of the two. Here we propose a simple statistical test of the neutral protein evolution hypothesis based on a comparison of the number of amino-acid replacement substitutions to synonymous substitutions in the coding region of a locus. If the observed substitutions are neutral, the ratio of replacement to synonymous fixed differences between species should be the same as the ratio of replacement to synonymous polymorphisms within species. DNA sequence data on the Adh locus (encoding alcohol dehydrogenase, EC 1.1.1.1) in three species in the Drosophila melanogaster species subgroup do not fit this expectation; instead, there are more fixed replacement differences between species than expected. We suggest that these excess replacement substitutions result from adaptive fixation of selectively advantageous mutations.

A number of different intracellular signaling pathways have been shown to be activated by receptor tyrosine kinases. These activation events include the phosphoinositide 3-kinase, 70 kDa S6 kinase, mitogen-activated protein kinase (MAPK), phospholipase C-gamma, and the Jak/STAT pathways. The precise role of each of these pathways in cell signaling remains to be resolved, but studies on the differentiation of mammalian PC12 cells in tissue culture and the genetics of cell fate determination in Drosophila and Caenorhabditis suggest that the extracellular signal-regulated kinase (ERK-regulated) MAPK pathway may be sufficient for these cellular responses. Experiments with PC12 cells also suggest that the duration of ERK activation is critical for cell signaling decisions.

Exogenous DNA sequences were introduced into the Drosophila germ line. A rosy transposon (ry1), constructed by inserting a chromosomal DNA fragment containing the wild-type rosy gene into a P transposable element, transformed germ line cells in 20 to 50 percent of the injected rosy mutant embryos. Transformants contained one or two copies of chromosomally integrated, intact ry1 that were stably inherited in subsequent generations. These transformed flies had wild-type eye color indicating that the visible genetic defect in the host strain could be fully and permanently corrected by the transferred gene. To demonstrate the generality of this approach, a DNA segment that does not confer a recognizable phenotype on recipients was also transferred into germ line chromosomes.

Toll-like receptors have a crucial role in the detection of microbial infection in mammals and insects. In mammals, these receptors have evolved to recognize conserved products unique to microbial metabolism. This specificity allows the Toll proteins to detect the presence of infection and to induce activation of inflammatory and antimicrobial innate immune responses. Recognition of microbial products by Toll-like receptors expressed on dendritic cells triggers functional maturation of dendritic cells and leads to initiation of antigen-specific adaptive immune responses.

The mammalian TOR (mTOR) pathway integrates nutrient- and growth factor-derived signals to regulate growth, the process whereby cells accumulate mass and increase in size. mTOR is a large protein kinase and the target of rapamycin, an immunosuppressant that also blocks vessel restenosis and has potential anticancer applications. mTOR interacts with the raptor and GbetaL proteins to form a complex that is the target of rapamycin. Here, we demonstrate that mTOR is also part of a distinct complex defined by the novel protein rictor (rapamycin-insensitive companion of mTOR). Rictor shares homology with the previously described pianissimo from D. discoidieum, STE20p from S. pombe, and AVO3p from S. cerevisiae. Interestingly, AVO3p is part of a rapamycin-insensitive TOR complex that does not contain the yeast homolog of raptor and signals to the actin cytoskeleton through PKC1. Consistent with this finding, the rictor-containing mTOR complex contains GbetaL but not raptor and it neither regulates the mTOR effector S6K1 nor is it bound by FKBP12-rapamycin. We find that the rictor-mTOR complex modulates the phosphorylation of Protein Kinase C alpha (PKCalpha) and the actin cytoskeleton, suggesting that this aspect of TOR signaling is conserved between yeast and mammals.

Alternative pre-mRNA splicing is a central mode of genetic regulation in higher eukaryotes. Variability in splicing patterns is a major source of protein diversity from the genome. In this review, I describe what is currently known of the molecular mechanisms that control changes in splice site choice. I start with the best-characterized systems from the Drosophila sex determination pathway, and then describe the regulators of other systems about whose mechanisms there is some data. How these regulators are combined into complex systems of tissue-specific splicing is discussed. In conclusion, very recent studies are presented that point to new directions for understanding alternative splicing and its mechanisms.

MicroRNAs are key regulators of gene expression, but the precise mechanisms underlying their interaction with their mRNA targets are still poorly understood. Here, we systematically investigate the role of target-site accessibility, as determined by base-pairing interactions within the mRNA, in microRNA target recognition. We experimentally show that mutations diminishing target accessibility substantially reduce microRNA-mediated translational repression, with effects comparable to those of mutations that disrupt sequence complementarity. We devise a parameter-free model for microRNA-target interaction that computes the difference between the free energy gained from the formation of the microRNA-target duplex and the energetic cost of unpairing the target to make it accessible to the microRNA. This model explains the variability in our experiments, predicts validated targets more accurately than existing algorithms, and shows that genomes accommodate site accessibility by preferentially positioning targets in highly accessible regions. Our study thus demonstrates that target accessibility is a critical factor in microRNA function.

We describe a genetic mosaic system in Drosophila, in which a dominant repressor of a cell marker is placed in trans to a mutant gene of interest. Mitotic recombination events between homologous chromosomes generate homozygous mutant cells, which are exclusively labeled due to loss of the repressor. Using this system, we are able to visualize axonal projections and dendritic elaboration in large neuroblast clones and single neuron clones with a membrane-targeted GFP marker. This new method allows for the study of gene functions in neuroblast proliferation, axon guidance, and dendritic elaboration in the complex central nervous system. As an example, we show that the short stop gene is required in mushroom body neurons for the extension and guidance of their axons.

Double-stranded RNA (dsRNA) induces sequence-specific posttranscriptional gene silencing in many organisms by a process known as RNA interference (RNAi). Using a Drosophila in vitro system, we demonstrate that 21- and 22-nt RNA fragments are the sequence-specific mediators of RNAi. The short interfering RNAs (siRNAs) are generated by an RNase III-like processing reaction from long dsRNA. Chemically synthesized siRNA duplexes with overhanging 3' ends mediate efficient target RNA cleavage in the lysate, and the cleavage site is located near the center of the region spanned by the guiding siRNA. Furthermore, we provide evidence that the direction of dsRNA processing determines whether sense or antisense target RNA can be cleaved by the siRNA-protein complex.

Comparison of relative fixation rates of synonymous (silent) and nonsynonymous (amino acid-altering) mutations provides a means for understanding the mechanisms of molecular sequence evolution. The nonsynonymous/synonymous rate ratio (omega = d(N)d(S)) is an important indicator of selective pressure at the protein level, with omega = 1 meaning neutral mutations, omega < 1 purifying selection, and omega>> 1 diversifying positive selection. Amino acid sites in a protein are expected to be under different selective pressures and have different underlying omega ratios. We develop models that account for heterogeneous omega ratios among amino acid sites and apply them to phylogenetic analyses of protein-coding DNA sequences. These models are useful for testing for adaptive molecular evolution and identifying amino acid sites under diversifying selection. Ten data sets of genes from nuclear, mitochondrial, and viral genomes are analyzed to estimate the distributions of omega among sites. In all data sets analyzed, the selective pressure indicated by the omega ratio is found to be highly heterogeneous among sites. Previously unsuspected Darwinian selection is detected in several genes in which the average omega ratio across sites is <1, but in which some sites are clearly under diversifying selection with omega>> 1. Genes undergoing positive selection include the beta-globin gene from vertebrates, mitochondrial protein-coding genes from hominoids, the hemagglutinin (HA) gene from human influenza virus A, and HIV-1 env, vif, and pol genes. Tests for the presence of positively selected sites and their subsequent identification appear quite robust to the specific distributional form assumed for omega and can be achieved using any of several models we implement. However, we encountered difficulties in estimating the precise distribution of omega among sites from real data sets.

Comparative analysis of multiple genomes in a phylogenetic framework dramatically improves the precision and sensitivity of evolutionary inference, producing more robust results than single-genome analyses can provide. The genomes of 12 Drosophila species, ten of which are presented here for the first time (sechellia, simulans, yakuba, erecta, ananassae, persimilis, willistoni, mojavensis, virilis and grimshawi), illustrate how rates and patterns of sequence divergence across taxa can illuminate evolutionary processes on a genomic scale. These genome sequences augment the formidable genetic tools that have made Drosophila melanogaster a pre-eminent model for animal genetics, and will further catalyse fundamental research on mechanisms of development, cell biology, genetics, disease, neurobiology, behaviour, physiology and evolution. Despite remarkable similarities among these Drosophila species, we identified many putatively non-neutral changes in protein-coding genes, non-coding RNA genes, and cis-regulatory regions. These may prove to underlie differences in the ecology and behaviour of these diverse species.