Genome-wide comparative analysis of simple sequence coding repeats among 25 insect species

We present a detailed genome-scale comparative analysis of simple sequence repeats within protein coding regions among 25 insect genomes. The repetitive sequences in the coding regions primarily represented single codon repeats and codon pair repeats. The CAG triplet is highly repetitive in the coding regions of insect genomes. It is frequently paired with the synonymous codon CAA to code for polyglutamine repeats. The codon pairs that are least repetitive code for polyalanine repeats. The frequency of hexanucleotide and dinucleotide motifs of codon pair repeats is significantly (p<0.001) different in the Drosophila species compared to the non-Drosophila species. However, the frequency of synonymous and non-synonymous codon pair repeats varies in a correlated manner (r(2)=0.79) among all the species. Results further show that perfect and imperfect repeats have significant association with the trinucleotide and hexanucleotide coding repeats in most of these insects. However, only select species show significant association between the numbers of perfect/imperfect hexamers and repeat coding for single amino acid/amino acid pair runs. Our data further suggests that genes containing simple sequence coding repeats may be under negative selection as they tend to be poorly conserved across species. The sequences of coding repeats of orthologous genes vary according to the known phylogeny among the species. In conclusion, the study shows that simple sequence coding repeats are important features of genome diversity among insects.     

SAGA: rapid automatic mainchain NMR assignment for large proteins

Here we describe a new algorithm for automatically determining the mainchain sequential assignment of NMR spectra for proteins. Using only the customary triple resonance experiments, assignments can be quickly found for not only small proteins having rather complete data, but also for large proteins, even when only half the residues can be assigned. The result of the calculation is not the single best assignment according to some criterion, but rather a large number of satisfactory assignments that are summarized in such a way as to help the user identify portions of the sequence that are assigned with confidence, vs. other portions where the assignment has some correlated alternatives. Thus very imperfect initial data can be used to suggest future experiments.     

Motif mismatches in microsatellites: insights from genome-wide investigation among 20 insect species

We present a detailed genome-wide comparative study of motif mismatches of microsatellites among 20 insect species representing five taxonomic orders. The results show that varying proportions (∼15–46%) of microsatellites identified in these species are imperfect in motif structure, and that they also vary in chromosomal distribution within genomes. It was observed that the genomic abundance of imperfect repeats is significantly associated with the length and number of motif mismatches of microsatellites. Furthermore, microsatellites with a higher number of mismatches tend to have lower abundance in the genome, suggesting that sequence heterogeneity of repeat motifs is a key determinant of genomic abundance of microsatellites. This relationship seems to be a general feature of microsatellites even in unrelated species such as yeast, roundworm, mouse and human. We provide a mechanistic explanation of the evolutionary link between motif heterogeneity and genomic abundance of microsatellites by examining the patterns of motif mismatches and allele sequences of single-nucleotide polymorphisms identified within microsatellite loci. Using Drosophila Reference Genetic Panel data, we further show that pattern of allelic variation modulates motif heterogeneity of microsatellites, and provide estimates of allele age of specific imperfect microsatellites found within protein-coding genes.     

RNAi-mediated inhibition of HIV-1 by targeting partially complementary viral sequences

Potent antiviral RNAi can be induced by intracellular expression of short hairpin RNAs (shRNAs) and artificial microRNAs (miRNAs). Expression of shRNA and miRNA results in target mRNA degradation (perfect base pairing) or translational repression (partial base pairing). Although efficient inhibition can be obtained, error-prone viruses such as human immunodeficiency virus type 1 (HIV-1) can escape from RNAi-mediated inhibition by mutating the target sequence. Recently, artificial miRNAs have been shown to be potent RNAi inducers due to their efficient processing by the RNAi machinery. Furthermore, miRNAs may be more proficient in suppressing imperfect targets than shRNAs. In this study, we tested the knockdown efficiency of miRNAs and shRNAs against wild-type and RNAi-escape HIV-1 variants with one or two mutations in the target sequence. ShRNAs and miRNAs can significantly inhibit the production of HIV-1 variants with mutated target sequences in the open reading frame. More pronounced mutation-tolerance was measured for targets in the 3′ untranslated region (3′ UTR). Partially complementary sequences within the 3′ UTR of the HIV-1 RNA genome efficiently act as target sites for miRNAs and shRNAs. These data suggest that targeting imperfect target sites by antiviral miRNAs or shRNAs provides an alternative RNAi approach for inhibition of pathogenic viruses.     

Analysis of clustered and longitudinal binary data subject to response misclassification

Misclassified clustered and longitudinal data arise in studies where the response indicates a condition identified through an imperfect diagnostic procedure. Examples include longitudinal studies that use an imperfect diagnostic test to assess whether or not an individual has been infected with a specific virus. This article presents methods to implement both population-averaged and cluster-specific analyses of such data when the misclassification rates are known. The methods exploit the fact that the class of generalized linear models enjoys a closure property in the case of misclassified responses. Data from longitudinal studies of infectious disease will illustrate the findings.     

Searching for Sequence Directed Mutagenesis in Eukaryotes

Sequence directed mutagenesis is a mechanism by which imperfect repeats “repair” each other to become perfect, generating mutations. This process is known to be prevalent in prokaryotes and it has been implicated in several human genetic diseases. Here we test whether sequence directed mutagenesis occurs in the protein coding sequences of eukaryotes using extensive DNA sequence data from humans, mice, Drosophila, nematodes, yeast, and Arabidopsis. Using two tests we find little evidence of sequence directed mutagenesis. We conclude that sequence directed mutagenesis is not prevalent in eukaryotes and that the examples of human diseases, apparently caused by sequence directed mutagenesis, are probably coincidental. [Reviewing Editor: Dr. Richard Kliman]     

Chromatin studies reveal that an ERE is located far upstream of a vitellogenin gene and that a distal tissue-specific hypersensitive site is conserved for two coordinately regulated vitellogenin genes.

Estrogen induces the expression of three vitellogenin genes in chicken hepatocytes. To survey the vitellogenin III (VTGIII) gene region for possible distal regulatory sequences, we identified tissue-specific hypersensitive (HS) sites within a 45 kb chromatin region spanning this gene. Five constitutive HS sites were found to mark the VTGIII gene region in hormone-naive hepatocytes. Strikingly, the constitutive HS site located 5.5 kb upstream of the VTGIII gene and a previously identified HS site located within the coordinately regulated VTGII gene mapped to nearly identical copies of a 72 bp sequence. Moreover, it would appear that there has been evolutionary pressure to retain specifically this 72 bp of VTGII-like sequence near the VTGIII gene subsequent to the VTGIII and VTGII genes becoming unlinked approximately 16 Myr ago. Two additional sets of HS sites were induced in the VTGIII gene region in response to estrogen. One set mapped immediately upstream of the gene in the vicinity of what we show to be a functional estrogen response element (ERE). The other induced HS site mapped 7.5 kb upstream of the gene. This far-upstream region was sequenced and was found to contain two imperfect ERE consensus sequences spaced 88 bp apart. In transient expression assays neither of these individual imperfect ERE sequences was functional, but a fragment spanning both sequences behaved as a strong ERE. In contrast to this synergism between imperfect ERE sequences, the presence of an NF-1 binding site 23 bp away from the more distal imperfect ERE sequence was not sufficient to render the latter a functional ERE in our assays.     

A variable minisatellite sequence in the chloroplast genome of Sorbus L. (Rosaceae: Maloideae).

The chloroplast genome is now known to be more variable than was once thought. Reports of RFLP (restriction fragment length polymorphism) and sequence variation, as well as variation in chloroplast microsatellites, are common. Here, data are presented on the variability of a minisatellite sequence in the chloroplast genome of Sorbus species. RFLP analysis of a PCR product comprising the region between the trnM and rbcL genes of nine Sorbus species identified seven size variants. Sequencing revealed the observed size polymorphism to be due to differences in the number of copies of an imperfect 9-bp motif. A more intensive survey of the variability of the minisatellite was undertaken in populations of Sorbus aucuparia. The potential uses of such regions in chloroplast DNA are discussed and a possible mechanism for the evolution of the minisatellite is presented.     

Patterns of sequence variation in the mitochondrial D-loop region of shrews

Direct sequencing of the mitochondrial displacement loop (D-loop) of shrews (genus Sorex) for the region between the tRNA(Pro) and the conserved sequence block-F revealed variable numbers of 79-bp tandem repeats. These repeats were found in all 19 individuals sequenced, representing three subspecies and one closely related species of the masked shrew group (Sorex cinereus cinereus, S. c. miscix, S. c. acadicus, and S. haydeni) and an outgroup, the pygmy shrew (S. hoyi). Each specimen also possessed an adjacent 76-bp imperfect copy of the tandem repeats. One individual was heteroplasmic for length variants consisting of five and seven copies of the 79-bp tandem repeat. The sequence of the repeats is conducive to the formation of secondary structure. A termination-associated sequence is present in each of the repeats and in a unique sequence region 5' to the tandem array as well. Mean genetic distance between the masked shrew taxa and the pygmy shrew was calculated separately for the unique sequence region, one of the tandem repeats, the imperfect repeat, and these three regions combined. The unique sequence region evolved more rapidly than the tandem repeats or the imperfect repeat. The small genetic distance between pairs of tandem repeats within an individual is consistent with a model of concerted evolution. Repeats are apparently duplicated and lost at a high rate, which tends to homogenize the tandem array. The rate of D- loop sequence divergence between the masked and pygmy shrews is estimated to be 15%-20%/Myr, the highest rate observed in D-loops of mammals. Rapid sequence evolution in shrews may be due either to their high metabolic rate and short generation time or to the presence of variable numbers of tandem repeats.      

Life history and the fitness consequences of imperfect information

The acquisition of information incurs costs in time, energy, exposure to predation, and/or lost opportunity. Without information, however, animals will be unable to assess the costs and benefits of decisions. Obtaining perfect information may be impossible, but how close to perfect do animals need assessments of ecological factors, such as predation risk, before estimation errors affect fitness? A recent article suggested that animals should be tolerant to imperfect information about predation risk, possibly relying on simple rules of thumb. Using a dynamic state variable approach, we examine some of the assumptions underlying this work, and show that tolerance towards imperfect information is dependent on life-history strategy. By changing the relationship between energy and fitness, assumptions about life-history strategies can be modified. Calculations show that tolerance to imperfect information is sensitive to these assumptions with some life histories leading to overestimation, while other life histories result in underestimation. One consistent effect across life histories is that animals with a higher rate of increase in fitness with respect to energy should show greater tolerance to imperfect information.     

Imperfect complementarity and RNA structure

Intra- and intermolecular complementary contacts in RNA are not always perfect: a significant amount of mismatch pairs is frequently found in naturally occurring RNA helices. The state of art in studies on mismatch pairs and examples of imperfect complementarity are reviewed. Two more cases are revealed by nucleotide sequence analysis techniques: imperfect complementary contacts Between ends of intervening sequences in eukaryotic mRNA precursors, and possible “stickiness” of mRNA to the ribosomes. The “stickiness” might arise from specific 3-Base periodicity of protein coding sequences which is found to be as universal as the code itself. The imperfect complementary contacts between mRNA and rRNA which monitor the coding frame provide a structural basis for the explanation of leaky frameshift phenomena.     

The effect of amino acid substitution in the imperfect repeat sequences of α-synuclein on fibrillation

Human α-synuclein is the causative protein of several neurodegenerative diseases, such as Parkinson's disease (PD) and dementia with Lewy Bodies (DLB). The N-terminal half of α-synuclein contains seven imperfect repeat sequences. One of the PD/DLB-causing point mutations, E46K, has been reported in the imperfect repeat sequences of α-synuclein, and is prone to form amyloid fibrils. The presence of seven imperfect repeats in α-synuclein raises the question of whether or not mutations corresponding to E46K in the other imperfect KTKE(Q)GV repeats have similar effects on aggregation and fibrillation, as well as their propensities to form α-helices. To investigate the effect of E(Q)/K mutations in each imperfect repeat sequence, we substituted the amino acid corresponding to E46K in each of the seven repeated sequences with a Lys residue. The mutations in the imperfect KTKE(Q)GV repeat sequences of the N-terminal region were prone to decrease the lag time of fibril formation. In addition, AFM imaging suggested that the Q24K mutant formed twisted fibrils, while the other mutants formed spherical aggregates and short fibrils. These observations indicate that the effect of the mutations on the kinetics of fibril formation and morphology of fibrils varies according to their location.     

Assessing the impacts of imperfect detection on estimates of diversity and community structure through multispecies occupancy modeling

Detecting all species in a given survey is challenging, regardless of sampling effort. This issue, more commonly known as imperfect detection, can have negative impacts on data quality and interpretation, most notably leading to false absences for rare or difficult‐to‐detect species. It is important that this issue be addressed, as estimates of species richness are critical to many areas of ecological research and management. In this study, we set out to determine the impacts of imperfect detection, and decisions about thresholds for inclusion in occupancy, on estimates of species richness and community structure. We collected data from a stream fish assemblage in Algonquin Provincial Park to be used as a representation of ecological communities. We then used multispecies occupancy modeling to estimate species‐specific occurrence probabilities while accounting for imperfect detection, thus creating a more informed dataset. This dataset was then compared to the original to see where differences occurred. In our analyses, we demonstrated that imperfect detection can lead to large changes in estimates of species richness at the site level and summarized differences in the community structure and sampling locations, represented through correspondence analyses.     

Origin and evolution of a major feline satellite DNA

A major satellite DNA has been cloned from the domestic cat (Felis catus) and characterized. The satellite monomer, termed FA-SAT, is 483 base-pairs in size, 64% G + C, and represents about 1 to 2% of the cat genome. A consensus sequence based upon partial sequence data from 21 independently isolated clones demonstrates: (1) FA-SAT is not composed of a series of shorter repeats, although about 25 copies, primarily imperfect, of the hexanucleotide TAACCC appear in the sequence; (2) there are many more CpG dinucleotides present in FA-SAT than expected for a random sequence of its size; and (3) 61% of all base substitutions in FA-SAT involve the replacement of G and C residues by A and T residues, indicating that FA-SAT is rapidly becoming A + T-rich. FA-SAT-related sequences are found in many mammals, where they appear to be scattered throughout the genome and not tandemly arranged as in the cat. An FA-SAT-related sequence was cloned from the domestic dog genome and sequenced, and shown to contain multiple copies of the same TAACCC hexanucleotide found in the cat satellite.     

A large-scale comparison of genomic sequences: one promising approach

We introduce a novel, linguistic-like method of genome analysis. We propose a natural approach to characterizing genomic sequences based on occurrences of fixed length words from a predefined, sufficiently large set of words (strings over the alphabet {A, C, G, T} ). A measure based on this approach is called compositional spectrum and is actually a histogram of imperfect word occurrences. Our results assert that the compositional spectrum is an overall characteristic of a long sequence i.e., a complete genome or an uninterrupted part of a chromosome. This attribute is manifested in the similarity of spectra obtained on different stretches of the same genome, and simultaneously in a broad range of dissimilarities between spectral representations of different genomes. High flexibility characterizes this approach due to imperfect matching and as a result sets of relatively long words can be considered. The proposed approach may have various applications in intra- and intergenomic sequence comparisons.     

Role of -35 sequence and its cooperativity with vir-box for the expression of virE gene

To elucidate the role of the -35 sequence and its cooperativity with vir box in the expression of the virE gene, various mutants were constructed by either site-directed mutation or deletional mutation of the virE promoter. The expression level of pHBAV, a mutant where its putative -35 sequences (CCGAGT) have been substituted with the consensus -35 sequences of the Escherichia coli promoter (TTGACA), was increased by 386%. pECHV, containing the conserved -35 sequence but lacking the vir box and the 5'-half of the imperfect dyad symmetry region (DSR) showed an increase of 286% in its promoter activity. pESHV, containing the conserved -35 sequence but lacking the complete 5'-upstream region from the mid-region of imperfect DSR, exhibited 244% of the native virE promoter activity. pHBCA, containing the conserved -35 sequence but destroying the vir box, was constructed by substitution of A, C, T at the positions -62, -63, and -65 on the vir-box to T, A, C, respectively. These mutations increased promoter activities by 319%. On the other hand, when the vir box was mutated from imperfect DSR to almost perfect DSR with T to A and G to T substitutions at -60 and -61 positions of the virE promoter containing the conserved -35 sequence (pHBNA), a higher activity of 671% was observed. These results demonstrate that when the putative -35 sequence of virE promoter is replaced with the consensus -35 sequence, the virE gene can be expressed independently without the binding of VirG protein to the vir-box and/or the induction of acetosyringone. Moreover, the presence of an almost perfect dyad symmetry of the vir-box can increase the expression of virE synergistically with the consensus -35 sequence.     

Nucleotide sequence analysis of IS427 and its target sites in Agrobacterium tumefaciens T37

We have determined the nucleotide sequence of IS427, an insertion sequence from Agrobacterium tumefaciens T37, IS427 is 1271 bp long, contains 16-bp imperfect terminal inverted repeats, and generates a 2-bp target sequence duplication. It is present at three sites in the pTiT37 plasmid and is absent from the chromosome of A. tumefaciens T37. Each of the IS427 elements sequenced was near a site with sequence homology to integration host factor (IHF)-binding sites which suggested that IHF may be involved in IS427 transposition.     

Entomopathogen ID: a curated sequence resource for entomopathogenic fungi

We report the development of a publicly accessible, curated nucleotide sequence database of hypocrealean entomopathogenic fungi. The goal is to provide a platform for users to easily access sequence data from taxonomic reference strains. The database can be used to accurately identify unknown entomopathogenic fungi based on sequence data for a variety of phylogenetically informative loci. The database provides full multi-locus sequence alignment capabilities. The initial release contains data compiled for 525 strains covering the phylogenetic diversity of three important entomopathogenic families: Clavicipitaceae, Cordycipitaceae, and Ophiocordycipitaceae. Furthermore, Entomopathogen ID can be expanded to other fungal clades of insect pathogens, as sequence data becomes available. The database will allow isolate characterisation and evolutionary analyses. We contend that this freely available, web-accessible database will facilitate the broader community to accurately identify fungal entomopathogens, which will allow users to communicate research results more effectively.