Messenger RNA sequence complexity and homology in developmental stages of Drosophila

Total polysomal RNA and polyadenylated mRNA from third instar larvae, pupae, and adults of D. melanogaster were hybridized in vast excess to labeled single-copy DNA in order to measure the sequence complexity of each RNA population. Then, to measure the sequence homology between the populations, each was hybridized to DNA enriched for messenger coding sequences in third instar larvae and to DNA depleted of these sequences. Our results show that a similar number of genes, approximately 16,000, is expressed in larvae, pupae, and adults, and that only one-third of these is expressed as polyadenylated mRNA. Further, the composition of both polyadenylated and nonpolyadenylated mRNA classes is shown to change very little between these three stages of development. Finally, the head of adult Drosophila is shown to contain 11,000 RNA species, approximately 70% of the number contained in the entire adult.     

Comparative analysis of cloned larval and adult globin cDNA sequences of Xenopus laevis

cDNA clones complementary to 9 S poly(A)⁺ RNA from erythroblasts of anemic larvae and adults of Xenopus laevis have been prepared. Clones, containing at least 400 bp of cDNA, have been analyzed by cross-hybridization and restriction mapping. They were found to comprise four unrelated main groups of sequences (two larval and two adult) and each main group contained two related subgroups. Partial sequence analysis and comparison of restriction data to previously published maps allowed the four main groups to be identified as α- or β-globin sequences. The sequence divergence between the subgroups was determined by melting curves of homo- and heteroduplexes. We found that the larval sequences have diverged twice as far as the adult ones. To account for this result, different hypotheses on globin gene evolution are proposed.     

Application of the Deoxyribonucleic Acid/Ribonucleic Acid Hybridization Technique in Bdellovibrio as a Model for Studying Ribonucleic Acid Turnover in Host-Parasite Systems

The kinetics of host ribonucleic acid (RNA) degradation and its resynthesis into Bdellovibrio-specific polyribonucleotides has been studied. The kinetics of RNA turnover was followed during a one-step synchronous growth cycle of Bdellovibrio growing within ³²PO₄-labeled Escherichia coli host cells. The species of labeled RNA present at any given time was ascertained through the specificity of the deoxyribonucleic acid (DNA)/RNA hybridization technique. At nearsaturating levels of RNA and at zero time, 7% of the host DNA sequences and only 0.04% of the Bdellovibrio DNA became hybridized with ³²P-labeled host cell RNA (greater than 99% host specific). At the end of the burst, 98% of the labeled RNA sequences were specific for Bdellovibrio DNA. About 74% of the initial labeled host cell RNA became turned over into Bdellovibrio-specific sequences. We provide data indicating that host cell ribosomal RNA is assimilated by Bdellovibrio. Degradation of host cell RNA occurs in a gradual fashion over most of the Bdellovibrio developmental growth cycle. This application of the DNA/RNA hybridization technique and its general concept should be of value in elucidating the kinetics of nucleic acid turnover in other types of host-parasite systems.     

A small nuclear precursor of messenger RNA in the cellular slime mold Dictyostelium discoideum

Previous work (Firtel et al., 1972) showed that messenger RNA from the cellular slime mold Dictyostelium discoideum, like that from mammalian cells, contains a sequence of about 100 adenylic acid residues at the 3′ end. We show here that Dictyostelium nuclei, labeled under a variety of conditions, do not contain material analogous to the large nuclear heterogeneous RNA found in mammalian cells. Rather, the majority of pulse-labeled nuclear RNA that is not a precursor of ribosomal RNA does contain at least one sequence of polyadenylic acid; this RNA, with an average molecular weight of 500,000, appears to be only 20% larger than cytoplasmic messenger RNA.Pulse-labeling experiments show that the nuclear poly(A)-containing RNA is a material precursor of messenger RNA. Whereas previous work showed that over 90% of messenger RNA sequences are transcribed from non-reiterated DNA, we show here that about 25% of nuclear poly (A)-containing RNA is transcribed from reiterated DNA sequences and only 75% from single-copy DNA. We present evidence that a large fraction of the nuclear poly(A)-containing RNA contains, at the 5′ end, a sequence of about 300 nucleotides that is transcribed from repetitive DNA, and which is lost before transport of messenger RNA into the cytoplasm.Based on these and other results, we present a model of arrangement of repetitive and single-copy DNA sequences in the Dictyostelium chromosome.     

Characterization of Sus scrofa Small Non-Coding RNAs Present in Both Female and Male Gonads

Small non-coding RNAs (sncRNAs) are indispensable for proper germ cell development, emphasizing the need for greater elucidation of the mechanisms of germline development and regulation of this process by sncRNAs. We used deep sequencing to characterize three families of small non-coding RNAs (piRNAs, miRNAs, and tRFs) present in Sus scrofa gonads and focused on the small RNA fraction present in both male and female gonads. Although similar numbers of reads were obtained from both types of gonads, the number of unique RNA sequences in the ovaries was several times lower. Of the sequences detected in the testes, 2.6% of piRNAs, 9% of miRNAs, and 10% of tRFs were also present in the ovaries. Notably, the majority of the shared piRNAs mapped to ribosomal RNAs and were derived from clustered loci. In addition, the most abundant miRNAs present in the ovaries and testes are conserved and are involved in many biological processes such as the regulation of homeobox genes, the control of cell proliferation, and carcinogenesis. Unexpectedly, we detected a novel sncRNA type, the tRFs, which are 30–36-nt RNA fragments derived from tRNA molecules, in gonads. Analysis of S. scrofa piRNAs show that testes specific piRNAs are biased for 5′ uracil but both testes and ovaries specific piRNAs are not biased for adenine at the 10^th nucleotide position. These observations indicate that adult porcine piRNAs are predominantly produced by a primary processing pathway or other mechanisms and secondary piRNAs generated by ping-pong mechanism are absent.     

In silico selection of RNA aptamers

In vitro selection of RNA aptamers that bind to a specific ligand usually begins with a random pool of RNA sequences. We propose a computational approach for designing a starting pool of RNA sequences for the selection of RNA aptamers for specific analyte binding. Our approach consists of three steps: (i) selection of RNA sequences based on their secondary structure, (ii) generating a library of three-dimensional (3D) structures of RNA molecules and (iii) high-throughput virtual screening of this library to select aptamers with binding affinity to a desired small molecule. We developed a set of criteria that allows one to select a sequence with potential binding affinity from a pool of random sequences and developed a protocol for RNA 3D structure prediction. As verification, we tested the performance of in silico selection on a set of six known aptamer–ligand complexes. The structures of the native sequences for the ligands in the testing set were among the top 5% of the selected structures. The proposed approach reduces the RNA sequences search space by four to five orders of magnitude—significantly accelerating the experimental screening and selection of high-affinity aptamers.     

Novel Virus Discovery and Genome Reconstruction from Field RNA Samples Reveals Highly Divergent Viruses in Dipteran Hosts

We investigated whether small RNA (sRNA) sequenced from field-collected mosquitoes and chironomids (Diptera) can be used as a proxy signature of viral prevalence within a range of species and viral groups, using sRNAs sequenced from wild-caught specimens, to inform total RNA deep sequencing of samples of particular interest. Using this strategy, we sequenced from adult Anopheles maculipennis s.l. mosquitoes the apparently nearly complete genome of one previously undescribed virus related to chronic bee paralysis virus, and, from a pool of Ochlerotatus caspius and Oc. detritus mosquitoes, a nearly complete entomobirnavirus genome. We also reconstructed long sequences (1503-6557 nt) related to at least nine other viruses. Crucially, several of the sequences detected were reconstructed from host organisms highly divergent from those in which related viruses have been previously isolated or discovered. It is clear that viral transmission and maintenance cycles in nature are likely to be significantly more complex and taxonomically diverse than previously expected.     

Comparison of an Avian Osteopetrosis Virus with an Avian Lymphomatosis Virus by RNA-DNA Hybridization

Myeloblastosis-associated virus (MAV)-2(0), a virus which was derived from avian myeloblastosis virus and induced a high incidence of osteopetrosis, was compared with avian lymphomatosis virus 5938, a recent field isolate which induced a high incidence of lymphomatosis. The following information was obtained. (i) MAV-2(0) induced osteopetrosis, nephroblastoma, and a very low incidence of hepatocellular carcinoma. No difference was seen in the oncogenic spectrum of end point and plaque-purified MAV-2(0). (ii) ¹²⁵I-labeled RNA sequences from MAV-2(0) formed hybrids with DNA extracted from osteopetrotic bone at a rate suggesting five proviral copies per haploid cell genome. The extent of hybridization of MAV-2(0) RNA with DNA from osteopetrotic tissue was more extensive (87%) than was observed in reactions with DNA from uninfected chicken embryos (52%). (iii) Competition of unlabeled viral RNA in hybridization reactions between the radioactive RNA from the two viruses and their respective proviral sequences present in tumor tissues showed that 15 to 20% of the viral sequences detected in these reactions were unshared. In contrast, no differences were detected in competition analyses of RNA sequences from the two viruses detected in DNA of normal chicken cells. (iv) MAV-2(0) 35S RNA was indistinguishable in size from avian lymphomatosis virus 5938 35S RNA by polyacrylamide gel electrophoresis.     

Ribonucleotide Sequence Homology Among Avian Oncornaviruses

RNA sequence relatedness among avian RNA tumor virus genomes was analyzed by inhibition of DNA-RNA hybrid formation between ³H-labeled 35S viral RNA and an excess of leukemic or normal chicken cell DNA with increasing concentrations of unlabeled 35S viral RNA. The avian viruses tested were Rous associated virus (RAV)-0, avian myeloblastosis virus (AMV), RAV-60, RAV-61, and B-77 sarcoma virus. Hybridization of ³H-labeled 35S AMV RNA with DNA from normal chicken cells was inhibited by unlabeled 35S RAV-0 RNA as efficiently (100%) as by unlabeled AMV RNA. Hybridization between ³H-labeled 35S AMV RNA and DNA from leukemic chicken myeloblasts induced by AMV was suppressed 100 and 68% by unlabeled 35S RNA from AMV and RAV-0, respectively. Hybridization between ³H-labeled RAV-0 and leukemic chicken myeloblast DNA was inhibited 100 and 67% by unlabeled 35S RNA from RAV-0 and AMV, respectively. It appears therefore that the AMV and RAV-0 genomes are 67 to 70% homologous and that AMV hybridizes to RAV-0 like sequences in normal chicken DNA. Hybridization between AMV RNA and leukemic chicken DNA was inhibited 40% by RNA from RAV-60 or RAV-61 and 50% by B-77 RNA. Hybridization between RAV-0 RNA and leukemic chicken DNA was inhibited 80% by RAV-60 or RAV-61 and 70% by B-77 RNA. Hybridization between ³H-labeled 35S RNA from RAV-60 or RAV-61 and leukemic chicken myeloblast DNA was reduced equally by RNA from RAV-60, RAV-61, AMV or RAV-0; this suggests that RNA from RAV-60 and RAV-61 hybridizes with virus-specific sequences in leukemic DNA which are shared by AMV, RAV-0, RAV-60, and RAV-61 RNAs. Hybridization between ³H-labeled 35S RNA from RAV-61 and normal pheasant DNA was inhibited 100% by homologous viral RNA, 22 to 26% by RNA from AMV or RAV-0, and 30 to 33% by RNA from RAV-60 or B-77. Nearly complete inhibition of hybridization between RAV-0 RNA and leukemic chicken DNA by a mixture of AMV and B-77 35S RNAs indicates that the RNA sequences shared by B-77 virus and RAV-0 are different from the sequences shared by AMV and RAV-0. It appears that different avian RNA tumor virus genomes have from 50 to 80% homology in nucleotide sequences and that the degree of hybridization between normal chicken cell DNA and a given viral RNA can be predicted from the homology that exists between the viral RNA tested and RAV-0 RNA.     

Next generation sequencing from Hepatozoon canis (Apicomplexa: Coccidia: Adeleorina): Complete apicoplast genome and multiple mitochondrion-associated sequences

Extrachromosomal genomes of the adeleorinid parasite Hepatozoon canis infecting an Israeli dog were investigated using next-generation and standard sequencing technologies. A complete apicoplast genome and several mitochondrion-associated sequences were generated. The apicoplast genome (31,869?bp) possessed two copies of both large subunit (23S) and small subunit (16S) ribosomal RNA genes (rDNA) within an inverted repeat region, as well as 22 protein-coding sequences, 25 transfer RNA genes (tDNA) and seven open reading frames of unknown function. Although circular-mapping, the apicoplast genome was physically linear according to next-generation data. Unlike other apicoplast genomes, genes encoding ribosomal protein S19 and tDNAs for alanine, aspartic acid, histidine, threonine and valine were not identified. No complete mitochondrial genome was recovered using next-generation data or directed PCR amplifications. Eight mitochondrion-associated (215–3523?bp) contigs assembled from next-generation data encoded a complete cytochrome c oxidase subunit I coding sequence, a complete cytochrome c oxidase subunit III coding sequence, two complete cytochrome B coding sequences, a non-coding, pseudogene for cytochrome B and multiple fragmented mitochondrial rDNA genes (SSUA, SSUB, SSUD, LSUC, LSUG, RNA6, RNA10, RNA14, RNA18). The paucity of NGS reads generating each of the mitochondrion-like sequences suggested that a complete mitochondrial genome at typically high copy number was absent in H. canis. In contrast, the complete nuclear rDNA unit sequence of H. canis (18S rDNA to 28S rDNA, 6977?bp) had >1000-fold next-generation coverage. Multiple divergent (from 93.6% to 99.9% pairwise identities) nuclear 18S rDNA contigs were generated (three types with 10 subtypes total). To our knowledge this is the first apicoplast genome sequenced from any adeleorinid coccidium and the first mitochondrion-associated sequences from this serious pathogen of wild and domestic canids. These newly generated sequences may provide useful genetic loci for high-resolution species-level genotyping that is currently impossible using existing nuclear rDNA targets.     

Determination of nucleotide sequences from double-stranded regions of HeLa cell nuclear RNA

Double-stranded regions which comprise about 4% of isolated HeLa cell heterogeneous nuclear RNA have been characterized by RNA fingerprinting and sequencing analysis. The simplicity of the pattern in two-dimensional RNA fingerprints suggests a sequence complexity of about 1000 nucleotides. The nucleotide sequences of six prominent RNase T₁-resistant oligonucleotides (ranging in size from 7 to 9 bases) have been determined using isolated double-stranded nuclear RNA labeled in vivo with ³²P-labeled inorganic phosphate. We conclude that (here exists a substantial subpopulation of simple, potentially complementary sequences common to much of the heterogeneous nuclear RNA population and interspersed with other kinds of sequences.     

Similar tissue-specific expression of the Adh genes from different Drosophila species is mediated by distinct arrangements of cis-acting sequences

In Drosophila melanogaster transformants, the alcohol dehydrogenase (Adh) genes from D. affinidisjuncta and D. grimshawi show similar levels of expression except in the adult midgut where the D. affinidisjuncta gene is expressed about 10- to 20-fold more strongly. To study the arrangement of cis-acting sequences responsible for this regulatory difference, homologous restriction sites were used to create a series of chimeric genes that switched fragments from the 5′ and 3′ flanking regions of these two genes. Chimeric genes were introduced into the germ-line of D. melanogaster, and Adh gene expression was analyzed by measuring RNA levels. Various gene fragments in the promoter region and elsewhere influence expression in the adult midgut and in whole larvae and adults. Comparison of these results with earlier studies involving chimeras between the D. affinidisjuncta and D. hawaiiensis genes indicates that expression in the adult midgut is influenced by multiple regulatory sequences and that distinct arrangements of regulatory sequences can result in similar levels of expression both in the adult midgut and in the whole organism.     

Developmental gene regulation in Aspergillus nidulans

The Ascomycete fungus Aspergillus nidulans reproduces asexually by differentiating conidiophores and conidia. Gene regulation during asexual reproduction was investigated by comparing poly(A) RNA populations derived from somatic hyphae, conidiating cultures and purified conidia. Single-copy and complementary DNA hybridization experiments showed that vegetative cells contained 5600–6000 diverse, average-sized poly(A) RNA sequences distributed into three prevalence classes. cDNA hybridization experiments indicated that a significant proportion of the poly(A) RNA derived from either conidiating cultures or spores consisted of sequences absent from somatic hyphae. To assess accurately the degree to which the poly(A) RNA populations differed, cDNA preparations were isolated which were complementary to sequences present only in conidia or in conidiating cultures. Hybridization of these cDNAs with poly(A) RNA from conidiating cultures showed that approximately 18.5% of the poly(A) RNA mass comprised 1300 diverse sequences not present in somatic cells. Of these, about 300 were present only in conidia. The remainder were accumulated specifically during sporulation, but were absent from spores. Analogous experiments showed that the great majority of the poly(A) RNA sequences accumulated by vegetative hyphae were also present in conidiating cultures. Thus, cell differentiation during A. nidulans asexual reproduction involves the accumulation of many new poly(A) RNA sequences, but not the loss of preexisting ones.     

Identification of development and tissue-specific gene expression in the fathead minnow Pimephales promelas, Rafinesque using computational and DNA microarray methods

The fathead minnow Pimephales promelas serves as a model organism for assessing the effects of environmental contaminants on early life stage growth and development. Yet, the utilization of genomic tools has been hindered by the lack of genome sequence and genomic information known from this model species. Utilizing published cDNA library sequences, the authors used sequence similarity to compare 4105 cDNAs isolated from fathead minnow fry (<14 days old) with over 250 000 adult cDNA sequences derived from whole body and various tissue types. The objectives of the computational subtraction were to (1) assess the extent of sequence similarity between developing and adult cDNA libraries and (2) predict which cDNA clones are expressed only in developing organisms. The results of the computational predictions were assessed through the construction of a development‐specific DNA microarray targeting all 4105 sequences in the fry cDNA library as well as 56 known mRNAs in P. promelas. Gene expression was determined by comparing total RNA isolated from fry with total RNA isolated from adult samples (whole animal, kidney, liver, brain, ovary and testes). The results showed that 1381 of the targeted fry cDNA sequences (34%) displayed expression across all sample comparisons, and of these, only 166 genes were found to harbour fry‐specific expression (i.e. no expression in adult samples). Of note, 69% of the genes computationally predicted to be fry specific were found across all experimental results; yet, only 27% of the computationally predicted fry‐specific sequences were experimentally confirmed to be fry specific. An important result was the identification of many novel mRNA sequences specific to the developing minnow, which lack homology with any other known sequence. In addition, the study results included tissue‐specific expression in adult samples. These results demonstrate the capabilities and limitations of inter‐library sequence comparisons as a predictor of gene activity in non‐sequenced organisms and tissues, as well as DNA microarray gene expression studies in non‐sequenced organisms.