We describe a new computer program, SnpEff, for rapidly categorizing the effects of variants in genome sequences. Once a genome is sequenced, SnpEff annotates variants based on their genomic locations and predicts coding effects. Annotated genomic locations include intronic, untranslated region, upstream, downstream, splice site, or intergenic regions. Coding effects such as synonymous or non-synonymous amino acid replacement, start codon gains or losses, stop codon gains or losses, or frame shifts can be predicted. Here the use of SnpEff is illustrated by annotating ~356,660 candidate SNPs in ~117 Mb unique sequences, representing a substitution rate of ~1/305 nucleotides, between the Drosophila melanogaster w1118; iso-2; iso-3 strain and the reference y1; cn1 bw1 sp1 strain. We show that ~15,842 SNPs are synonymous and ~4,467 SNPs are non-synonymous (N/S ~0.28). The remaining SNPs are in other categories, such as stop codon gains (38 SNPs), stop codon losses (8 SNPs), and start codon gains (297 SNPs) in...

A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3

https://www.cell.com/cell/pdf/0092-8674(93)90530-4.pdf

Posttranscriptional regulation of the heterochronic gene lin-14 by lin-4 mediates temporal pattern formation in C. elegans

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.

https://science.sciencemag.org/content/sci/298/5595/1039.full.pdf

Role of Histone H3 Lysine 27 Methylation in Polycomb-Group Silencing

The origin of organismal complexity is generally thought to be tightly coupled to the evolution of new gene functions arising subsequent to gene duplication. Under the classical model for the evolution of duplicate genes, one member of the duplicated pair usually degenerates within a few million years by accumulating deleterious mutations, while the other duplicate retains the original function. This model further predicts that on rare occasions, one duplicate may acquire a new adaptive function, resulting in the preservation of both members of the pair, one with the new function and the other retaining the old. However, empirical data suggest that a much greater proportion of gene duplicates is preserved than predicted by the classical model. Here we present a new conceptual framework for understanding the evolution of duplicate genes that may help explain this conundrum. Focusing on the regulatory complexity of eukaryotic genes, we show how complementary degenerative mutations in different regulatory elements of duplicated genes can facilitate the preservation of both duplicates, thereby increasing long-term opportunities for the evolution of new gene functions. The duplication-degeneration-complementation (DDC) model predicts that (1) degenerative mutations in regulatory elements can increase rather than reduce the probability of duplicate gene preservation and (2) the usual mechanism of duplicate gene preservation is the partitioning of ancestral functions rather than the evolution of new functions. We present several examples (including analysis of a new engrailed gene in zebrafish) that appear to be consistent with the DDC model, and we suggest several analytical and experimental approaches for determining whether the complementary loss of gene subfunctions or the acquisition of novel functions are likely to be the primary mechanisms for the preservation of gene duplicates. For a newly duplicated paralog, survival depends on the outcome of the race between entropic decay and chance acquisition of an advantageous regulatory mutation. Sidow (1996, p. 717) On one hand, it may fix an advantageous allele giving it a slightly different, and selectable, function from its original copy. This initial fixation provides substantial protection against future fixation of null mutations, allowing additional mutations to accumulate that refine functional differentiation. Alternatively, a duplicate locus can instead first fix a null allele, becoming a pseudogene. Walsh (1995, p. 426) Duplicated genes persist only if mutations create new and essential protein functions, an event that is predicted to occur rarely. Nadeau and Sankoff (1997, p. 1259) Thus overall, with complex metazoans, the major mechanism for retention of ancient gene duplicates would appear to have been the acquisition of novel expression sites for developmental genes, with its accompanying opportunity for new gene roles underlying the progressive extension of development itself. Cooke et al. (1997, p. 362)

https://www.genetics.org/content/genetics/151/4/1531.full.pdf

Preservation of Duplicate Genes by Complementary, Degenerative Mutations

Background: The recent discoveries of microRNA (miRNA) genes and characterization of the first few target genes regulated by miRNAs in Caenorhabditis elegans and Drosophila melanogaster have set the stage for elucidation of a novel network of regulatory control. We present a computational method for wholegenome prediction of miRNA target genes. The method is validated using known examples. For each miRNA, target genes are selected on the basis of three properties: sequence complementarity using a position-weighted local alignment algorithm, free energies of RNA-RNA duplexes, and conservation of target sites in related genomes. Application to the D. melanogaster, Drosophila pseudoobscura and Anopheles gambiae genomes identifies several hundred target genes potentially regulated by one or more known miRNAs.

/pdf/microrna-targets-in-drosophila-51ilqbrm67.pdf

MicroRNA targets in Drosophila

The metameric organization of the Drosophila embryo is generated in the first 5 h after fertilization. An initially rather simple pattern provides the foundation for subsequent development and diversification of the segmented part of the body. Many of the genes that control the formation of this pattern have been identified and at least twenty have been cloned. By combining the techniques of genetics, molecular biology and experimental embryology, it is becoming possible to unravel the role played by each of these genes. The repeating segment pattern is defined by the persistent expression of engrailed and of other genes of the 'segment polarity' class. The establishment of this pattern is directed by a transient molecular prepattern that is generated in the blastoderm by the activity of the 'pair-rule' genes. Maternal determinants at the poles of the egg coordinate this prepattern and define the anteroposterior sequence of pattern elements. The primary effect of these determinants is not known, but genes required for their production have been identified and the product of one of these, bicoid is known to be localized at the anterior of the egg. One early consequence of their activity is to define domains along the A-P axis within which a series of 'cardinal' genes are transcribed. The activity of the cardinal genes is required both to coordinate the process of segmentation and to define the early domains of homeotic gene expression. Further interactions between the homeotic genes and other classes of segmentation genes refine the initial establishment of segment identities.

The molecular basis for metameric pattern in the Drosophila embryo.

Germ cells play a unique role in gamete production, heredity and evolution. Therefore, to understand the mechanisms that specify germ cells is a central challenge in developmental and evolutionary biology. Data from model organisms show that germ cells can be specified either by maternally inherited determinants (preformation) or by inductive signals (epigenesis). Here we review existing data on 28 metazoan phyla, which indicate that although preformation is seen in most model organisms, it is actually the less prevalent mode of germ cell specification, and that epigenetic germ cell specification may be ancestral to the Metazoa.

/pdf/mechanisms-of-germ-cell-specification-across-the-metazoans-30rvm8f91n.pdf

Mechanisms of germ cell specification across the metazoans: epigenesis and preformation

Oligonucleotide primers have been used to generate a cDNA library covering the entire tobacco mosaic virus (TMV) RNA sequence. Analysis of these clones has enabled us to complete the viral RNA sequence and to study its variability within a viral population. The positive strand coding sequence starts 69 nucleotides from the 5' end with a reading frame for a protein of Mr 125,941 and terminates with UAG. Readthrough of this terminator would give rise to a protein of Mr 183,253. Overlapping the terminal five codons of this readthrough reading frame is a second reading frame coding for a protein of Mr 29,987. This gene terminates two nucleotides before the initiator codon of the coat protein gene. Potential signal sequences responsible for the capping and synthesis of the coat protein and Mr 29,987 protein mRNAs have been identified. Similar sequences within these reading frames may be used in the expression of sets of proteins that share COOH-terminal sequences.

https://www.pnas.org/content/pnas/79/19/5818.full.pdf

Nucleotide sequence of tobacco mosaic virus RNA

The bithorax complex in Drosophila melanogaster is a cluster of homeotic genes that specify developmental pathways for many of the body segments of the fly. The DNA of the bithorax complex has been isolated, and a region of 195,000 base pairs that covers the left half of the complex is described here. The lesions associated with many of the bithorax complex mutants have been identified, and most are due to DNA rearrangements. Most of the spontaneous mutants have insertions of a particular mobile element named "gypsy." This element affects the functions of sequences removed from the site of insertion. Mutant lesions for a given phenotypic class are distributed over large DNA distances of up to 73,000 base pairs.

https://science.sciencemag.org/content/sci/221/4605/23.full.pdf

Molecular Genetics of the Bithorax Complex in Drosophila melanogaster

Understanding the early evolution of animal body plans requires knowledge both of metazoan phylogeny and of the genetic and developmental changes involved in the emergence of particular forms. Recent 18S ribosomal RNA phylogenies suggest a three-branched tree for the Bilateria comprising the deuterostomes and two great protostome clades, the lophotrochozoans1 and ecdysozoans2. Here, we show that the complement of Hox genes in critical protostome phyla reflects these phylogenetic relationships and reveals the early evolution of developmental regulatory potential in bilaterians. We have identified Hox genes that are shared by subsets of protostome phyla. These include a diverged pair of posterior (Abdominal-B -like) genes in both a brachiopod and a polychaete annelid, which supports the lophotrochozoan assemblage, and a distinct posterior Hox gene shared by a priapulid, a nematode and the arthropods, which supports the ecdysozoan clade. The ancestors of each of these two major protostome lineages had a minimum of eight to ten Hox genes. The major period of Hox gene expansion and diversification thus occurred before the radiation of each of the three great bilaterian clades.

Michael Akam

Papers

The molecular basis for metameric pattern in the Drosophila embryo.

Mechanisms of germ cell specification across the metazoans: epigenesis and preformation

Nucleotide sequence of tobacco mosaic virus RNA

Molecular Genetics of the Bithorax Complex in Drosophila melanogaster

Hox genes in brachiopods and priapulids and protostome evolution