Structural conservation and variation in the D-loop-containing region of vertebrate mitochondrial DNA

The nucleotide sequences of the D-loop-containing regions of three rat mitochondrial DNAs (mtDNAs), two from the species Rattus norvegicus and one from R. rattus, were determined. Comparisons made among these sequences and with the mouse sequence showed that, on the basis of both base composition and frequency of nucleotide alterations, three domains could be defined within the D-loop-containing region: a central conserved segment, poor in L-strand adenine, flanked by two divergent, adenine-rich regions. Deletions and insertions were found to occur at an unexpectedly high frequency in these sequences and the conserved sequence block called CSB-1 was found not to be intact in the R. rattus sequence. Although in comparisons of more distantly related mtDNAs the D-loop region is the most divergent on the molecule, it does not diverge more than typical protein genes between R. norvegicus and R. rattus, and its central conserved domain appears to be one of the molecule's most conserved regions. The most variable domain borders the tRNAPhe gene and contains the L and H-strand promoters and the 5' terminus for H-strand DNA synthesis. Within this region we have found sequences in all the mtDNAs we have examined, including those of human, two artiodactyls and Xenopus, that are capable of folding into cloverleaf structures. In the other divergent domain of the same mtDNAs, we find sequences capable of assuming similar secondary structural configurations at or near the sites for the termination of D-loop DNA synthesis. The evolutionary preservation of the potential to form such structures despite the high primary-structural divergence of the regions they occur in, suggests the structures are of principal importance for some processes occurring in the D-loop-containing region.     

Establishment and maintenance of nuclear position during zoospore formation in allomyces macrogynus: roles of the cytoskeleton

The involvement of the microtubule (MT) and actin microfilament (MF) cytoskeletons in establishing nuclear positions during zoosporogenesis in Allomyces macrogynus was assessed using selective cytoskeletal disrupting treatments and documented with light microscopy. These experiments were coupled with low-speed centrifugation studies to determine the degree to which cytoskeletal elements anchor nuclear position. At the onset of zoospore formation, nuclei were positioned only in cortical cytoplasmic regions of the zoosporangia (ZS). Immunofluorescence microscopy revealed that MTs primarily emanated from centrosomal regions into the surrounding cytoplasm at this stage. During delimitation of the cytoplasm into individual uninucleate zoospores, nuclei migrated from cortical regions to become distributed throughout the cytoplasm. Coincident with nuclear migrations, MTs were primarily organized at and emanated from nuclear surfaces, forming extensive perinuclear arrays. Nuclear migrations were suppressed in ZS induced to sporulate in the presence of cytochalasin D, an actin MF inhibiting compound. Disruption of MTs with nocodazole did not block nuclear migrations, although resultant nuclear spacing was irregular. Centrifugation treatments of control and drug-treated ZS demonstrated that nuclear positions were stabilized by perinuclear MT arrays. The results indicate that nuclear motility in ZS of A. macrogynus is the result of an actin-based system while perinuclear MTs arrays function to establish and fix nuclear position during zoospore formation. Copyright 1998 Academic Press.     

The ability to form intrastrand tetraplexes is an evolutionarily conserved feature of the 3' end of L1 retrotransposons

Mammalian genomes contain many thousands of members of a family of retrotransponsons known as L1 (or LINE-1) elements. These elements lack long terminal repeats (LTRs), and are thought to use a retroposition mechanism than differs from that of retroviruses and other LTR- containing retroelements. In order to define those regions of the L1 element that may be important for L1 retroposition, we examined the 3' untranslated regions (UTRs) of L1 elements from a diverse group of mammals. We show that while the 3' UTRs of L1 elements from different species share little if any sequence homology, they all contain a G- rich polypurine tract of variable length and sequence which can form one or more intrastrand tetraplexes. This conservation over the 100 Myr since the mammalian radiation suggests that either the G-rich motif itself or a conserved structure such as the tetraplex that can be formed by this motif is a significant structural feature of L1 elements and may play a role in their propagation.      

Retroviral and pseudogene insertion sites reveal the lineage of human salivary and pancreatic amylase genes from a single gene during primate evolution.

We have analyzed the junction regions of inserted elements within the human amylase gene complex. This complex contains five genes which are expressed at high levels either in the pancreas or in the parotid gland. The proximal 5'-flanking regions of these genes contain two inserted elements. A gamma-actin pseudogene is located at a position 200 base pairs upstream of the first coding exon. All of the amylase genes contain this insert. The subsequent insertion of an endogenous retrovirus interrupted the gamma-actin pseudogene within its 3'-untranslated region. Nucleotide sequence analysis of the inserted elements associated with each of the five human amylase genes has revealed a series of molecular events during the recent history of this gene family. The data indicate that the entire gene family was generated during primate evolution from one ancestral gene copy and that the retroviral insertion activated a cryptic promoter.     

The 28S ribosomal RNA-encoding gene of Hymenoptera: inserted sequences in the retrotransposon-rich regions.

The genomes of two parasitoid wasps, Diadromus pulchellus and Eupelmus vuilleti, and the honey bee, Apis mellifera, contain few interspersed repeated sequences corresponding to transposons (Tn). This suggests that the genomic organisation of Hymenoptera could be due to the elimination of deleterious Tn in haploid males. We have used restriction-fragment length polymorphism analysis to show that nondeleterious Tn are present in the DNA (rDNA) encoding ribosomal RNA of twelve species of Hymenoptera. Sequence analysis of the 28S rDNA type-I and type-II insertion-rich regions of 80 species showed that this region is very highly conserved (95.8%). A consensus sequence and restriction map of the rDNA region were established. These sequence data were used to develop a strategy for detecting inserted elements in the rDNA fragments containing type-I or type-II insertion sites, and this strategy was used to screen twelve hymenopteran species and four non-Hymenoptera control species. The rDNA fragments from the Hymenoptera and control species contained inserted sequences in the area where type-I and type-II elements are inserted in the 28S rDNA retrotransposon-rich region of Diptera and Lepidoptera. The hymenopteran genomes therefore appear to contain repeated elements, the mobility and nature of which remain to be determined.     

Evaluation and improvements in the automatic alignment of protein sequences

The accuracy of protein sequence alignment obtained by applying a commonly used global sequence comparison algorithm is assessed. Alignments based on the superposition of the three-dimensional structures are used as a standard for testing the automatic, sequence-based methods. Alignments obtained from the global comparison of five pairs of homologous protein sequences studied gave 54% agreement overall for residues in secondary structures. The inclusion of information about the secondary structure of one of the proteins in order to limit the number of gaps inserted in regions of secondary structure, improved this figure to 68%. A similarity score of greater than six standard deviation units suggests that an alignment which is greater than 75% correct within secondary structural regions can be obtained automatically for the pair of sequences.     

Improving protein secondary structure prediction with aligned homologous sequences.

Most recent protein secondary structure prediction methods use sequence alignments to improve the prediction quality. We investigate the relationship between the location of secondary structural elements, gaps, and variable residue positions in multiple sequence alignments. We further investigate how these relationships compare with those found in structurally aligned protein families. We show how such associations may be used to improve the quality of prediction of the secondary structure elements, using the Quadratic-Logistic method with profiles. Furthermore, we analyze the extent to which the number of homologous sequences influences the quality of prediction. The analysis of variable residue positions shows that surprisingly, helical regions exhibit greater variability than do coil regions, which are generally thought to be the most common secondary structure elements in loops. However, the correlation between variability and the presence of helices does not significantly improve prediction quality. Gaps are a distinct signal for coil regions. Increasing the coil propensity for those residues occurring in gap regions enhances the overall prediction quality. Prediction accuracy increases initially with the number of homologues, but changes negligibly as the number of homologues exceeds about 14. The alignment quality affects the prediction more than other factors, hence a careful selection and alignment of even a small number of homologues can lead to significant improvements in prediction accuracy.     

Functional necessity of the cytoskeleton during cleavage membrane development and zoosporogenesis in Allomyces macrogynus

Cleavage membrane development and cytokinesis were examined in zoosporangia of Allomyces macrogynus treated with cytoskeletal inhibitors and compared to zoosporogenesis under control conditions. Developing membranes were visualized in living zoosporangia with laser-scanning confocal microscopy using the lipophilic membrane dye FM4-64. Under control conditions, cleavage membranes developed in four discrete stages, ultimately interconnecting to delimit the cytoplasm into polygonal uninucleate domains of near uniform size. Disruption of microtubules did not impede the normal four-stage development of cleavage membranes, and cytokinesis occurred with only minor detectable anomalies, although zoospores lacked flagella. Disruption of actin microfilaments did not inhibit membrane formation but blocked nuclear migration and significantly disrupted membrane alignment and cytoplasmic delimitation. This resulted in masses of membrane that remained primarily in cortical regions of the zoosporangia, as did nuclei, throughout zoosporogenesis. Zoospores formed in the absence of microtubules had only a slightly larger mean diameter than control zoospores, although nearly 50% of spores contained two or more nuclei. Microfilament inhibitor treatments produced spores with substantially larger mean diameters and correspondingly larger numbers of nuclei per spore, with greater than 85% containing three or more nuclei. These results showed that a functional actin microfilament cytoskeleton was required for proper alignment of cleavage elements and cytokinesis in Allomyces zoosporangia while microtubules played a less significant role.     

A proposal for the secondary structure of a variable area of eukaryotic small ribosomal subunit RNA involving the existence of a pseudoknot.

Eukaryotic small ribosomal subunit RNAs contain an area of variable structure, V4, which comprises about 250 nucleotides in most species, whereas the corresponding area in bacterial small ribosomal subunit RNAs consists of about 64 nucleotides folded into a single hairpin. There is no consensus on the secondary structure of area V4 in eukaryotes, about 10 different models having been proposed. The prediction of a model on a comparative basis poses special problems because, due to the variability of the area in length as well as sequence, a dependable alignment is very difficult to achieve. A new model was derived by systematic examination of all combinations of helices that have been hitherto proposed, plus some new ones. The following properties of the helices were examined: transposability to all presently known sequences, presence of compensating substitutions, and thermodynamic stability. A model was selected by ranking all possible combinations of transposable helices according to the number of compensating substitutions scored. The optimal model comprises a pseudoknot and four hairpin structures. Certain species contain additional hairpins inserted between these structural elements, while in others the structure is partially or entirely deleted.     

Intron/exon structure of the chicken pyruvate kinase gene

The chicken pyruvate kinase gene is interrupted by at least ten introns, including nine introns within the coding region. We compare the structure of this gene with the three-dimensional protein structure of the homologous cat muscle enzyme. The introns are not randomly placed--they divide the coding sequence into fairly uniformly sized pieces encoding discrete elements of secondary structure. The introns tend to fall at interruptions between stretches of alpha-helix or beta-sheet residues, and each of the six exons that contribute to the barrel-shaped central domain include one or two repeats of a simple unit, an alpha-helix plus a beta strand. This structure suggests that introns were not inserted into a previously uninterrupted coding sequence, but instead are products of the evolution of the first pyruvate kinase gene. We have found some sequence homology between a segment of pyruvate kinase and the structurally homologous mononucleotide binding fold of alcohol dehydrogenase. The superposition of these two regions aligns an intron from the maize alcohol dehydrogenase gene four nucleotides from an intron in the chicken pyruvate kinase gene.     

Gene rearrangements in snake mitochondrial genomes: highly concerted evolution of control-region-like sequences duplicated and inserted into a tRNA gene cluster

Mitochondrial DNA (mtDNA) regions corresponding to two major tRNA gene clusters were amplified and sequenced for the Japanese pit viper, himehabu. In one of these clusters, which in most vertebrates characterized to date contains three tightly connected genes for tRNA(Ile), and tRNA(Gln), and tRNA(Met), a sequence of approximately 1.3 kb was found to be inserted between the genes for tRNA(Ile) and tRNA(Gln). The insert consists of a control-region-like sequence possessing some conserved sequence blocks, and short flanking sequences which may be folded into tRNA(Pro), tRNA(Phe), and tRNA(Leu) genes. Several other snakes belonging to different families were also found to possess a control-region-like sequence and tRNA(Leu) gene between the tRNA(Ile)and tRNA(Gln) genes. We also sequenced a region surrounded by genes for cytochrome b and 12S rRNA, where the control region and genes for tRNA(Pro) and tRNA(Phe) are normally located in the mtDNAs of most vertebrates. In this region of three examined snakes, a control-region- like sequence exists that is almost completely identical to the one found between the tRNA(Ile) and tRNA(Gln) genes. The mtDNAs of these snakes thus possess two nearly identical control-region-like sequences which are otherwise divergent to a large extent between the species. These results suggest that the duplicate state of the control-region- like sequences has long persisted in snake mtDNAs, possibly since the original insertion of the control-region-like sequence and tRNA(Leu) gene into the tRNA gene cluster, which occurred in the early stage of the divergence of snakes. It is also suggested that the duplicated control-region-like sequences at two distant locations of mtDNA have evolved concertedly by a mechanism such as frequent gene conversion. The secondary structures of the determined tRNA genes point to the operation of simplification pressure on the T psi C arm of snake mitochondrial tRNAs.      

Cell wall synthesis in Allomyces macrogynus

Scanning electron microscopy and freeze-etching/cleaving have been employed to examine events in the synchronized development of gametophytic germlings of the aquatic Phycomycete Allomyces macrogynus. Motile spores were induced to start synchronized development and the sequence of surface changes associated with the encystment process was studied. Time course studies show that small vesicles (apparently blebbed off from the gamma-particles) start to accumulate on the surface of the plasma membrane after 6 min of synchronized growth at the same time as the first cell wall material can be detected. The vesicles increase in number during encystment. After 15 min of synchronized growth the number of vesicles decrease and after 20 min of growth no vesicle can be observed on the cell surface. During this period the cell surface appears increasingly smooth, probably due to cell wall formation. In freeze-etching/cleaving electron micrographs from this period, both intact and what appear to be ruptured vesicles outside the cell surface, can be observed. The intact vesicle has a characteristic surface pattern presumably of membrane particles. This surface view of the encystment processes supports the hypothesis that the gamma-particles through gamma vesicle formation participate in the cell wall synthesis during encystment in Allomyces.     

Conserved secondary structures in Aspergillus

Background

Recent evidence suggests that the number and variety of functional RNAs (ncRNAs as well as cis-acting RNA elements within mRNAs ) is much higher than previously thought; thus, the ability to computationally predict and analyze RNAs has taken on new importance. We have computationally studied the secondary structures in an alignment of six Aspergillus genomes. Little is known about the RNAs present in this set of fungi, and this diverse set of genomes has an optimal level of sequence conservation for observing the correlated evolution of base-pairs seen in RNAs.

Methodology/Principal Findings

We report the results of a whole-genome search for evolutionarily conserved secondary structures, as well as the results of clustering these predicted secondary structures by structural similarity. We find a total of 7450 predicted secondary structures, including a new predicted ∼60 bp long hairpin motif found primarily inside introns. We find no evidence for microRNAs. Different types of genomic regions are over-represented in different classes of predicted secondary structures. Exons contain the longest motifs (primarily long, branched hairpins), 5′ UTRs primarily contain groupings of short hairpins located near the start codon, and 3′ UTRs contain very little secondary structure compared to other regions. There is a large concentration of short hairpins just inside the boundaries of exons. The density of predicted intronic RNAs increases with the length of introns, and the density of predicted secondary structures within mRNA coding regions increases with the number of introns in a gene.

Conclusions/Sigificance

There are many conserved, high-confidence RNAs of unknown function in these Aspergillus genomes, as well as interesting spatial distributions of predicted secondary structures. This study increases our knowledge of secondary structure in these aspergillus organisms.     

Upstream elements present in the 3'-untranslated region of collagen genes influence the processing efficiency of overlapping polyadenylation signals

3'-Untranslated regions (UTRs) of genes often contain key regulatory elements involved in gene expression control. A high degree of evolutionary conservation in regions of the 3'-UTR suggests important, conserved elements. In particular, we are interested in those elements involved in regulation of 3' end formation. In addition to canonical sequence elements, auxiliary sequences likely play an important role in determining the polyadenylation efficiency of mammalian pre-mRNAs. We identified highly conserved sequence elements upstream of the AAUAAA in three human collagen genes, COL1A1, COL1A2, and COL2A1, and demonstrate that these upstream sequence elements (USEs) influence polyadenylation efficiency. Mutation of the USEs decreases polyadenylation efficiency both in vitro and in vivo, and inclusion of competitor oligoribonucleotides representing the USEs specifically inhibit polyadenylation. We have also shown that insertion of a USE into a weak polyadenylation signal can enhance 3' end formation. Close inspection of the COL1A2 3'-UTR reveals an unusual feature of two closely spaced, competing polyadenylation signals. Taken together, these data demonstrate that USEs are important auxiliary polyadenylation elements in mammalian genes.     

Nucleotide composition and kinetic complexity of the genomic DNA of an intracellular symbiont in the pea aphidAcyrthosiphon pisum

Summary A detailed comparative study of the regions surrounding the origin of replication in vertebrate mitochondrial DNA (mtDNA) has revealed a number of interesting properties. This region, called the D-loop-containing region, can be divided into three domains. The left (L) and right (R) domains, which have a low G content and contain the 5 and the 3 D-loop ends, respectively, are highly variable for both base sequence and length. They, however, contain thermodynamically stable secondary structures which include the conserved sequence blocks called CSB-1 and TAS which are associated with the start and stop sites, respectively, for D-loop strand synthesis. We have found that a mirror symmetry exists between the CSB-1 and TAS elements, which suggests that they can act as specific recognition sites for regulatory, probably dimeric, proteins. Long, statistically significant repeats are found in the L and R domains.Between the L and R domains we observed in all mtDNA sequences a region with a higher G content which was apparently free of complex secondary structure. This central domain, well preserved in mammals, contains an open reading frame of variable length in the organisms considered.The identification of common features well preserved in evolution despite the high primary structural divergence of the D-loop-containing region of vertebrate mtDNA suggests that these properties are of prime importance for the mitochondrial processes that occur in this region and may be useful for singling out the sites on which one should operate experimentally in order to discover functionally important elements.Presented at the FEBS Symposium on Genome Organization and Evolution, held in Crete, Greece, September 1–5, 1986     

A CT promoter element binding protein: definition of a double-strand and a novel single-strand DNA binding motif.

Numerous genes contain promoter elements that are nuclease hypersensitive. These elements frequently possess polypurine/polypyrimidine stretches and are usually associated with altered chromatin structure. We have previously isolated a clone that binds a class of CT-rich promoter elements. We have further characterized this clone, termed the nuclease-sensitive element protein-1, or NSEP-1. NSEP-1 binds both duplex CT elements and the CT-rich strand of these elements in a 'generic' sequence specific manner and has overlapping but distinct single-and double-strand DNA binding domains. The minimal peptide region sufficient for both duplex and single-strand DNA binding includes two regions rich in basic amino acids flanking an RNP-CS-1 like octapeptide motif. Deletion analysis shows that the single-strand DNA binding activity is mediated by the RNP-CS-1 like octapeptide motif and is the key peptide region necessary for single-strand binding. NSEP-1's affinity for CT rich promoter elements with strand asymmetry in addition to its double- and single-strand DNA binding properties suggests that it may be a member of a class of DNA binding proteins that modulate gene expression by their ability to recognize DNA with unusual secondary structure.