Messenger RNA sequence complexity and homology in developmental stages of Drosophila

Ultrasensitive proteome analysis using paramagnetic bead technology

In order to obtain a systems-level understanding of a complex biological system, detailed proteome information is essential. Despite great progress in proteomics technologies, thorough interrogation of the proteome from quantity-limited biological samples is hampered by inefficiencies during processing. To address these challenges, here we introduce a novel protocol using paramagnetic beads, termed Single-Pot Solid-Phase-enhanced Sample Preparation (SP3). SP3 provides a rapid and unbiased means of proteomic sample preparation in a single tube that facilitates ultrasensitive analysis by outperforming existing protocols in terms of efficiency, scalability, speed, throughput, and flexibility. To illustrate these benefits, characterization of 1,000 HeLa cells and single Drosophila embryos is used to establish that SP3 provides an enhanced platform for profiling proteomes derived from sub-microgram amounts of material. These data present a first view of developmental stage-specific proteome dynamics in Drosophila at a single-embryo resolution, permitting characterization of inter-individual expression variation. Together, the findings of this work position SP3 as a superior protocol that facilitates exciting new directions in multiple areas of proteomics ranging from developmental biology to clinical applications.

Genomic basis of ecological niche divergence among cryptic sister species of non-biting midges

BACKGROUND

There is a lack of understanding the evolutionary forces driving niche segregation of closely related organisms. In addition, pinpointing the genes driving ecological divergence is a key goal in molecular ecology. Here, larval transcriptome sequences obtained by next-generation-sequencing are used to address these issues in a morphologically cryptic sister species pair of non-biting midges (Chironomus riparius and C. piger).

RESULTS

More than eight thousand orthologous open reading frames were screened for interspecific divergence and intraspecific polymorphisms. Despite a small mean sequence divergence of 1.53% between the sister species, 25.1% of 18,115 observed amino acid substitutions were inferred by α statistics to be driven by positive selection. Applying McDonald-Kreitman tests to 715 alignments of gene orthologues identified eleven (1.5%) genes driven by positive selection.

CONCLUSIONS

Three candidate genes were identified as potentially responsible for the observed niche segregation concerning nitrite concentration, habitat temperature and water conductivity. Additionally, signs of positive selection in the hydrogen sulfide detoxification pathway were detected, providing a new plausible hypothesis for the species' ecological differentiation. Finally, a divergently selected, nuclear encoded mitochondrial ribosomal protein may contribute to reproductive isolation due to cytonuclear coevolution.

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.

Protein evolution in the context of Drosophila development

The tempo at which a protein evolves depends not only on the rate at which mutations arise but also on the selective effects that those mutations have at the organismal level. It is intuitive that proteins functioning during different stages of development may be predisposed to having mutations of different selective effects. For example, it has been hypothesized that changes to proteins expressed during early development should have larger phenotypic consequences because later stages depend on them. Conversely, changes to proteins expressed much later in development should have smaller consequences at the organismal level. Here we assess whether proteins expressed at different times during Drosophila development vary systematically in their rates of evolution. We find that proteins expressed early in development and particularly during mid-late embryonic development evolve unusually slowly. In addition, proteins expressed in adult males show an elevated evolutionary rate. These two trends are independent of each other and cannot be explained by peculiar rates of mutation or levels of codon bias. Moreover, the observed patterns appear to hold across several functional classes of genes, although the exact developmental time of the slowest protein evolution differs among each class. We discuss our results in connection with data on the evolution of development.

Gene expression during the life cycle of Drosophila melanogaster

Molecular genetic studies of Drosophila melanogaster have led to profound advances in understanding the regulation of development. Here we report gene expression patterns for nearly one-third of all Drosophila genes during a complete time course of development. Mutations that eliminate eye or germline tissue were used to further analyze tissue-specific gene expression programs. These studies define major characteristics of the transcriptional programs that underlie the life cycle, compare development in males and females, and show that large-scale gene expression data collected from whole animals can be used to identify genes expressed in particular tissues and organs or genes involved in specific biological and biochemical processes.

Courtship and visual defects of cacophony mutants reveal functional complexity of a calcium-channel alpha1 subunit in Drosophila

We show by molecular analysis of behavioral and physiological mutants that the Drosophila Dmca1A calcium-channel alpha1 subunit is encoded by the cacophony (cac) gene and that nightblind-A and lethal(1)L13 mutations are allelic to cac with respect to an expanded array of behavioral and physiological phenotypes associated with this gene. The cacS mutant, which exhibits defects in the patterning of courtship lovesong and a newly revealed but subtle abnormality in visual physiology, is mutated such that a highly conserved phenylalanine (in one of the quasi-homologous intrapolypeptide regions called IIIS6) is replaced by isoleucine. The cacH18 mutant exhibits defects in visual physiology (including complete unresponsiveness to light in certain genetic combinations) and visually mediated behaviors; this mutant (originally nbAH18) has a stop codon in an alternative exon (within the cac ORF), which is differentially expressed in the eye. Analysis of the various courtship and visual phenotypes associated with this array of cac mutants demonstrates that Dmca1A calcium channels mediate multiple, separable biological functions; these correlate in part with transcript diversity generated via alternative splicing.

Molecular and behavioral analysis of four period mutants in Drosophila melanogaster encompassing extreme short, novel long, and unorthodox arrhythmic types

Of the mutationally defined rhythm genes in Drosophila melanogaster, period (per) has been studied the most. We have molecularly characterized three older per mutants-perT, perClk, and per04-along with a novel long-period one (perSLIH). Each mutant is the result of a single nucleotide change. perT, perClk, and perSLIH are accounted for by amino acid substitutions; per04 is altered at a splice site acceptor and causes aberrant splicing. perSLIH exhibits a long period of 27 hr in constant darkness and entrains to light/dark (L/D) cycles with a later-than-normal evening peak of locomotion. perSLIH males are more rhythmic than females. perSLIH's clock runs faster at higher temperatures and slower at lower ones, exhibiting a temperature-compensation defect opposite to that of perLong. The per-encoded protein (PER) in the perT mutant cycles in L/D with an earlier-than-normal peak; this peak in perSLIH is later than normal, and there was a slight difference in the PER timecourse of males vs. females. PER in per04 was undetectable. Two of these mutations, perSLIH and perClk, lie within regions of PER that have not been studied previously and may define important functional domains of this clock protein.

Defective learning in mutants of the Drosophila gene for a regulatory subunit of cAMP-dependent protein kinase

Disruptions of a Drosophila gene encoding a regulatory subunit of cAMP-dependent protein kinase homologous to mammalian RIbeta (dPKA-RI) were targeted to the first (noncoding) exon of dPKA-RI via site-selected P element mutagenesis. Flies homozygous for either of two mutant alleles showed specific defects in olfactory learning but not in subsequent memory decay. In contrast, olfactory acuity and shock reactivity, component behaviors required for normal odor avoidance learning, were normal in these mutants. Northern and Western blot analyses of mRNA and protein extracted from adult heads have revealed a complex lesion of the PKA-RI locus, including expression of a novel product and over- or underexpression of wild-type products in mutants. Western blot analysis revealed reductions in RI protein in mutants. PKA activity in the absence of exogenous cAMP also was significantly higher than normal in homogenates from mutant adult heads. These two mutant alleles failed to complement each other for each of these phenotypic defects, eliminating second-site mutations as a possible explanation. These results establish a role for an RI regulatory subunit of PKA in Pavlovian olfactory conditioning.

An on-line two-dimensional polyacrylamide gel electrophoresis protein database of adult Drosophila melanogaster

An annotated two-dimensional polyacrylamide gel electrophoresis (2-D PAGE) protein database of adult Drosophila melanogaster has been constructed, based on the protein patterns of heads, thoraces and abdomens of adult male and female Drosophila melanogaster. About 1200 major protein spots are catalogued. Common proteins, found in all body parts, as well as bodypart- and sex-specifically expressed proteins are reported. Of the major proteins, 91, or 7.5%, are differentially expressed in the two sexes or in different body parts, at least in part reflecting specific functional requirements. At the present time 43 proteins, or about 3.5% of the detected proteins, have been identified. These data can be accessed interactively from our World Wide Web (WWW) server through clickable inline gel images and hypertext links. Identified protein spots are cross-referenced, through hypertext links, to the SWISS-PROT annotated database of protein primary sequences and the Fly-Base database of Drosophila genomic data. Our reference gels can be used to gain immediate access to protein spot identify and to the pattern of differentially expressed proteins in Drosophila melanogaster. The work presented in this article ties together information from protein 2-D PAGE, molecular biology and genetics and offers a uniform way to access this large volume of data.

Gene number, noise reduction and biological complexity

Preliminary estimates suggest that gene number, and hence biological complexity, increased suddenly at two periods of macroevolutionary change (the origin of eukaryotes and the origin of vertebrates), but otherwise remained relatively constant. As the genome is in constant flux, what normally constrains the number of different genes that an organism can retain? Here, I suggest that an important limitation on gene number is the efficiency of mechanisms that reduce transcriptional background noise. The appearance of both eukaryotes and vertebrates coincided with novel mechanisms of noise reduction.

The non-dosage compensated LSP1-alpha gene of Drosophila melanogaster lies immediately downstream of the dosage compensated L12 gene

The X-linked gene LSP1-alpha of Drosophila melanogaster, expressed in the third larval instar, does not exhibit dosage compensation at its normal locus but does compensate when it is relocated to ectopic sites on the X chromosome. A transcription unit designated L12, which is active in the second larval instar and capable of encoding a putative protein of 28.5 kDa, lies immediately downstream from LSP1-alpha. We have determined that L12 is dosage compensated by measuring the steady-state level of its transcript in male and female larvae. The difference in response of these two adjacent genes should be taken into consideration when models of the mechanism of dosage compensation are formulated.

Transforming potential of a myc-containing variant of feline leukemia virus in vitro in early-passage feline cells

We studied a naturally occurring variant of feline leukemia virus (FeLV) in which the oncogene myc has substituted for a portion of the viral structural genes (myc-FeLV). myc-FeLV was rescued by replication in the presence of FeLV as helper, and its biological activity was examined in early-passage feline cells in vitro. Infection of leukocytes from peripheral blood, spleen, or thymus, or of kitten fibroblasts did not immortalize these cells or alter them morphologically. Northern blot (RNA blot) analysis of virion RNA prepared from the supernatant of infected cells demonstrated the 8.2-kilobase genome of FeLV, but did not demonstrate the 5.0-kilobase genome of myc-FeLV. Apparently, the myc-FeLV genome was lost in the absence of the selective pressure of transformation. In contrast, infection of embryonic fibroblasts with myc-FeLV(FeLV) rendered these cells capable of greatly increased, if not infinite, proliferative potential. The cells were morphologically altered compared with controls and were only loosely adherent to the substrate. The cells failed to proliferate in semisolid medium and did not form tumors when inoculated subcutaneously into athymic mice. Blot analyses demonstrated the presence and expression of integrated proviral DNAs of both FeLV and myc-FeLV in these cells. They appear, then, to represent cells partially transformed by infection with myc-FeLV(FeLV). The action of feline v-myc in early-passage cells in vitro was compared to that of avian v-myc.

Molecular cloning of chromosome I DNA from Saccharomyces cerevisiae: isolation and characterization of the CDC24 gene and adjacent regions of the chromosome

Molecular cloning techniques were used to isolate and characterize the DNA including and surrounding the CDC24 and PYK1 genes on the left arm of chromosome I of the yeast Saccharomyces cerevisiae. A plasmid that complemented a temperature-sensitive cdc24 mutation was isolated from a yeast genomic DNA library in a shuttle vector. Plasmids containing pyk1-complementing DNA were obtained from other investigators. Several lines of evidence (including one-step gene replacement experiments) demonstrated that the complementing plasmids contained the bona fide CDC24 and PYK1 genes. These sequences were then used to isolate additional DNA from chromosome I by probing a yeast genomic DNA library in a lambda vector. A total of 28 kilobases (kb) of contiguous DNA surrounding the CDC24 and PYK1 genes was isolated, and a restriction map was determined. Electron microscopy of R-loop-containing DNA and RNA blot hybridization analyses indicated that an 18-kb segment contained at least seven transcribed regions, only three of which corresponded to previously known genes (CDC24, PYK1, and CYC3). Southern blot hybridization experiments suggested that none of the genes in this region was duplicated elsewhere in the yeast genome. The centers of CDC24 and PYK1 were only approximately 7.5 kb apart, although the genetic map distance between them is approximately 13 centimorgans. As previous studies with S. cerevisiae have indicated that 1 centimorgan generally corresponds to approximately 3 kb, the region between CDC24 and PYK1 appears to undergo meiotic recombination at an unusually high frequency.

Molecular cloning of chromosome I DNA from Saccharomyces cerevisiae: isolation and characterization of the CDC24 gene and adjacent regions of the chromosome

Molecular cloning techniques were used to isolate and characterize the DNA including and surrounding the CDC24 and PYK1 genes on the left arm of chromosome I of the yeast Saccharomyces cerevisiae. A plasmid that complemented a temperature-sensitive cdc24 mutation was isolated from a yeast genomic DNA library in a shuttle vector. Plasmids containing pyk1-complementing DNA were obtained from other investigators. Several lines of evidence (including one-step gene replacement experiments) demonstrated that the complementing plasmids contained the bona fide CDC24 and PYK1 genes. These sequences were then used to isolate additional DNA from chromosome I by probing a yeast genomic DNA library in a lambda vector. A total of 28 kilobases (kb) of contiguous DNA surrounding the CDC24 and PYK1 genes was isolated, and a restriction map was determined. Electron microscopy of R-loop-containing DNA and RNA blot hybridization analyses indicated that an 18-kb segment contained at least seven transcribed regions, only three of which corresponded to previously known genes (CDC24, PYK1, and CYC3). Southern blot hybridization experiments suggested that none of the genes in this region was duplicated elsewhere in the yeast genome. The centers of CDC24 and PYK1 were only approximately 7.5 kb apart, although the genetic map distance between them is approximately 13 centimorgans. As previous studies with S. cerevisiae have indicated that 1 centimorgan generally corresponds to approximately 3 kb, the region between CDC24 and PYK1 appears to undergo meiotic recombination at an unusually high frequency.

Cloning of a major surface-antigen gene of Trypanosoma cruzi and identification of a nonapeptide repeat

The parasitic protozoan Trypanosoma cruzi can establish infection in humans and other vertebrate hosts through direct penetration of host cells by trypomastigotes transmitted by the insect vector. Although the molecular processes involved in trypomastigote interiorization of vertebrate cells are unknown, several studies suggest that surface glycoproteins are involved. It is likely that the proteins involved are specific to the trypomastigote stage of the parasite, since only trypomastigotes found in both the insect vector and the vertebrate host bloodstream are capable of invading vertebrate cells. In contrast, the epimastigote stage, found exclusively in the vector, and the amastigote stage, an intracellular stage in the vertebrate host, cannot penetrate the cell directly. We have therefore concentrated our efforts on trypomastigote surface proteins and, along with others, have identified two trypomastigote-specific surface glycoproteins of relative molecular mass (Mr) 90,000 (90K) and 85,000 (85K). Antibody neutralization experiments indicate that the 85K glycoprotein is necessary for efficient interiorization of trypomastigotes in mammalian cells. Here we describe the molecular cloning of a genomic DNA fragment that encodes antigenic determinants present in the 85K trypomastigote surface antigen. The polypeptide fragment encoded by the cloned DNA is recognized by serum from a T. cruzi-infected host and is inferred by DNA sequence analysis to contain a nonapeptide unit that is tandemly repeated five times. Also, the messenger complementary to the cloned DNA fragment is present only in the trypomastigote stage of the parasite.

Polyadenylation state of abundant mRNAs during Drosophila development

We have used a two-dimensional gel analysis of cell-free translation products to determine whether individual mRNAs present in Drosophila melanogaster embryos, larvae, pupae, and adults are predominantly polyadenylated or nonadenylated. While the majority of the embryonic mRNAs we detected exist mainly in the polyadenylated form, these mRNAs become more evenly distributed between the poly(A)+ and poly(A)- RNA fractions during postembryonic development. Although DNA:RNA hybridization experiments have indicated that Drosophila RNA populations contain a large group of rare class mRNAs restricted to the poly(A)- RNA compartment, this is not true for the 150 more abundant mRNA species analyzed by our methods. The histone mRNAs are the only abundant mRNA species which appear to be exclusively in the poly(A)- RNA class.

Isolation and characterization of the glucose-6-phosphate dehydrogenase gene of Drosophila melanogaster

To investigate the molecular basis of dosage compensation in Drosophila, a recombinant lambda phage containing the Drosophila melanogaster glucose-6-phosphatase dehydrogenase (G6PD) gene was isolated by differential screening of a Drosophila genomic lambda library with poly(A)+RNA obtained from polyribosomes enriched for or depleted of G6PD mRNA sequences. Of 44 000 plaques screened, a single phage, lambda DmG21, showed hybridization with the enriched poly(A)+RNA but not the depleted one. Confirmation that the Drosophila DNA fragment cloned in lambda DmG21 contains the G6PD gene sequence is based on the following observations. lambda DmG21 DNA shows hybridization only to the 18D region of the salivary gland X-chromosome, which is the known cytological locus for the G6PD gene. In vitro translation of the poly(A)+mRNA selected by hybridization to lambda DmG21 DNA sequences shows a polypeptide product of apparent Mr 55 000, identical to that of the monomeric unit of G6PD. When the putative coding sequence of G6PD is cloned into the expression vector lambda gt11, recombinant plaques are recognized by anti-G6PD immunoglobulin. A transcriptional map of the G6PD gene shows that it is divided into two exons, 0.9 kb (exon I) and 1.8 kb (exon II) long, which are separated by a 2.4-kb intron. The G6PD mRNA is 2.0 kb in length and the steady-state level of the mRNA is similar in both sexes. Measurement of the copy number of the G6PD gene in males and females shows the gene to be present once per X-chromosome in both sexes. No amplification of the gene sequence was observed in males. These results are, therefore, in agreement with the previous suggestion that dosage compensation is the result of enhanced transcription of X-linked genes in males.

Genetic activity along 315 kb of the Drosophila chromosome

Transcripts from different tissues were mapped along a 315 kb segment of the Drosophila chromosome, a region which includes the rosy and Ace loci. Forty-three distinct RNA species were detected, though only 12 recessive lethal complementation groups had been mapped in the interval. The sum of the sizes of the transcripts covers 33% of the genomic DNA. The distribution of transcription units along the walk is very uneven. Sixty-three kb of genomic DNA at the proximal end of the walk encode 18 transcripts, while only seven are found in the next 153 kb. Each tissue exhibits a specific spectrum of transcripts. No clustering was seen among genes expressed coordinately. In salivary glands, the number of transcripts detected corresponds to the number of chromomeric units in the polytene chromosomes of this tissue. Moreover, the density distribution of transcripts along the DNA walk is parallel to the density distribution of chromomeric units.

Molecular analysis of the period locus in Drosophila melanogaster and identification of a transcript involved in biological rhythms

We have isolated and analyzed DNA sequences encompassing the period (per) locus of Drosophila melanogaster. The location of this clock gene was delimited by the molecular mapping of chromosome aberrations at or very near the per locus. At least five RNAs are transcribed from this region. One of these transcripts, a 0.9 kb species, is strongly implicated in per's control of biological rhythms. Two independently isolated arrhythmic mutations at the per locus dramatically reduce the level of this transcript. Furthermore, the level of the 0.9 kb transcript is strongly modulated during a light/dark cycle. We discuss evidence, from previously reported genetic and phenotypic analysis of per's function, suggesting that this region may be complex and that several gene products from the per region, including this 0.9 kb transcript, may be involved in the different aspects of normal rhythmicity influenced by this clock gene.

Molecular and genetic studies on the euchromatin-heterochromatin transition region of the X chromosome of Drosophila melanogaster. 1. A cloned entry point near to the uncoordinated (unc) locus

A recombinant Charon 4 bacteriophage has been isolated on the basis of RNAs which are enriched in the head of the adult Drosophila melanogaster and hence are likely to be of neural origin. The cloned insert maps to the near vicinity of the uncoordinated locus in polytene chromosome band 19E8. This band is within the transition zone between the euchromatic and heterochromatic regions of the X chromosome, a region which has been well characterized cytogenetically. The insert contains both repetitious and low copy number sequences, some of which vary extensively in both frequency and restriction fragment size between different laboratory strains. One particular family of moderately repeated sequences occurs predominantly in divisions 19 and 20 of the X chromosome and perhaps the distally located X heterochromatin. The molecular landscape surrounding the initial entry point contains many repeated sequences and is thus unlike those observed in most published chromosomal walks. The possible significance of the presence of repeated families in the distinct properties of this region are discussed.

Genetic expression in the developing brain

The adult mouse brain contains complex populations of polyadenylated [poly(A)+] and nonpolyadenylated [poly(A)-] messenger RNA's (mRNA's). These mRNA's are separate sequence populations, similar in complexity, and in combination are equivalent to approximately 150,000 different mRNA sequences, of average length. Essentially all of the "adult" poly(A)+ mRNA's are present in the brain at birth. In contrast, most of the poly(A)- mRNA's are absent. Brain poly(A)- mRNA's begin to appear soon after birth, but the full adult complement is not reached until young adulthood. This suggests that these poly(A)- mRNA's specify proteins required for the biological capabilities of the brain that emerge during the course of postnatal development.

Nonadenylylated mRNA is present as polyadenylylated RNA in nuclei of Drosophila

The sequence complexity of nuclear total RNA and nuclear poly(A)+RNA from Drosophila third-instar larvae was determined by hybridization of these RNAs to labeled single-copy DNA. At saturation, the nuclear poly(A)+ - and total RNA hybridized to 11% and 22.5% of the single-copy DNA, respectively. The increase in complexity of nuclear total RNA over that observed for nuclear poly(A)+RNA indicates the presence of a discrete class of nonoadenylylated nuclear RNA molecules. The relationship between DNA sequences coding for nuclear RNA and mRNA was then determined by hybridization of nuclear total and poly(A)+RNA to DNA enriched for mRNA coding sequences. The results of these studies show that those single-copy DNA sequences that are represented in either the poly(A)+ - or poly(A)- mRNA population are transcribed into RNA molecules that appear in the nuclear poly(A)+RNA population.