Expressed sequence tags from the midgut and an epithelial cell line of Chironomus tentans: annotation, bioinformatic classification of unknown transcripts and analysis of expression levels

Expressed sequence tags (ESTs) were generated from two Chironomus tentans cDNA libraries, constructed from an embryo epithelial cell line and from larva midgut tissue. 8584 5'-end ESTs were generated and assembled into 3110 tentative unique transcripts, providing the largest contribution of C. tentans sequences to public databases to date. Annotation using Blast gave 1975 (63.5%) transcripts with a significant match in the major gene/protein databases, 1170 with a best match to Anopheles gambiae and 480 to Drosophila melanogaster. 1091 transcripts (35.1%) had no match to any database. Studies of open reading frames suggest that at least 323 of these contain a coding sequence, indicating that a large proportion of the genes in C. tentans belong to previously unknown gene families.

ESTs from brain and testis of White Leghorn and red junglefowl: annotation, bioinformatic classification of unknown transcripts and analysis of expression levels

We report the generation, assembly and annotation of expressed sequence tags (ESTs) from four chicken cDNA libraries, constructed from brain and testis tissue dissected from red junglefowl and White Leghorn. 21,285 5'-end ESTs were generated and assembled into 2,813 contigs and 9,737 singletons, giving 12,549 tentative unique transcripts. The transcripts were annotated using BLAST by matching to known chicken genes or to putative homologues in other species using the major gene/protein databases. The results for these similarity searches are available on www.sbc.su.se/~arve/chicken. 4,129 (32.9%) of the transcripts remained without a significant match to gene/protein databases, a proportion of unmatched transcripts similar to earlier non-mammalian EST studies. To estimate how many of these transcripts may represent novel genes, they were studied for the presence of coding sequence. It was shown that most of the unique chicken transcripts do not contain coding parts of genes, but it was estimated that at least 400 of the transcripts contain coding sequence, indicating that 3.2% of avian genes belong to previously unknown gene families. Further BLAST search against dbEST left 1,649 (13.1%) of the transcripts unmatched to any library. The number of completely unmatched transcripts containing coding sequence was estimated at 180, giving a measure of the number of putative novel chicken genes identified in this study. 84.3% of the identified transcripts were found only in testis tissue, which has been poorly studied in earlier chicken EST studies. Large differences in expression levels were found between the brain and testis libraries for a large number of transcripts, and among the 525 most frequently represented transcripts, there were at least 20 transcripts with significant difference in expression levels between red junglefowl and White Leghorn.

A novel method for forensic DNA investigations: repeat-type sequence analysis of tandemly repeated mtDNA in domestic dogs

A highly variable and heteroplasmic tandem repeat region situated in the mitochondrial mt DNA control region (CR) in domestic dogs and wolves was studied to evaluate its suitability as a forensic genetic marker for analysis of single hairs. The tandem repeat array is composed of three 10-bp repeat types that are distributed so that a secondary DNA sequence is formed. Thus, the region presents two levels of variation: variation in the number of repeats and variation in the secondary DNA sequence of repeat types. Two analysis methods were therefore tested; fragment length analysis and analysis of the sequence of repeat types. Fragment analysis produced unique profiles that could be used to discriminate between blood samples from maternally closely related individuals. However, different hairs from one individual did not have the same fragment profile, and the method is, therefore, not suitable for analysis of single hairs. In contrast, analysis of the repeat type sequences (array types) is highly informative. When different hairs from one individual were studied, identical array types were found. The repeat-type sequence variation was studied among individuals having identical nonrepetitive CR mtDNA sequence variants. Seven, six, and two individuals, representing three different sequence variants, respectively, were analyzed. All these individuals had different array types, which implies a very high genetic variation between individuals in this region. The analysis method considerably improves the exclusion capacity of mtDNA analysis of domestic dogs compared with sequence analysis of non-repetitive DNA.

mtDNA tandem repeats in domestic dogs and wolves: mutation mechanism studied by analysis of the sequence of imperfect repeats

The mitochondrial (mt) DNA control region (CR) of dogs and wolves contains an array of imperfect 10 bp tandem repeats. This region was studied for 14 domestic dogs representing the four major phylogenetic groups of nonrepetitive CR and for 5 wolves. Three repeat types were found among these individuals, distributed so that different sequences of the repeat types were formed in different molecules. This enabled a detailed study of the arrays and of the mutation events that they undergo. Extensive heteroplasmy was observed in all individuals; 85 different array types were found in one individual, and the total number of types was estimated at 384. Among unrelated individuals, no identical molecules were found, indicating a high rate of evolution of the region. By performing a pedigree analysis, array types which had been inherited from mother to offspring and array types which were the result of somatic mutations, respectively, could be identified, showing that about 20% of the molecules within an individual had somatic mutations. By direct pairwise comparison of the mutated and the original array types, the physiognomy of the inserted or deleted elements (indels) and the approximate positions of the mutations could be determined. All mutations could be explained by replication slippage or point mutations. The majority of the indels were 1-5 repeats long, but deletions of up to 17 repeats were found. Mutations were found in all parts of the arrays, but at a higher frequency in the 5' end. Furthermore, the inherited array types within the mother-offspring pair were aligned and compared so that germ line mutations could be studied. The pattern of the germ line mutations was approximately the same as that of the somatic mutations.

Estimation of reversible substitution matrices from multiple pairs of sequences

We present a method for estimating the most general reversible substitution matrix corresponding to a given collection of pairwise aligned DNA sequences. This matrix can then be used to calculate evolutionary distances between pairs of sequences in the collection. If only two sequences are considered, our method is equivalent to that of Lanave et al. (1984). The main novelty of our approach is in combining data from different sequence pairs. We describe a weighting method for pairs of taxa related by a known tree that results in uniform weights for all branches. Our method for estimating the rate matrix results in fast execution times, even on large data sets, and does not require knowledge of the phylogenetic relationships among sequences. In a test case on a primate pseudogene, the matrix we arrived at resembles one obtained using maximum likelihood, and the resulting distance measure is shown to have better linearity than is obtained in a less general model.