Gene regulation and evolution in the chorion locus of Bombyx mori. Structural and developmental characterization of four eggshell genes and their flanking DNA regions

We report on the detailed structural and developmental characterization of four chorion genes and a truncated pseudogene located within a 9.5 X 10(3) base chromosomal segment. These genes belong to the A and B multigene families and, like previously characterized moth chorion genes, are arranged in tightly linked pairs, which are divergently transcribed (A/B.L11 and A/B.L12). On the basis of their high degree of sequence divergence, the A genes define two distinct subfamilies, while the more homologous B genes represent different copies of the same gene type. The A.L11 and B.L11 introns are much longer, in each case because of a single inserted DNA segment that is missing from A.L12 or B.L12. The 2.1 X 10(3) base insertion in A.L11 is the first retrovirus-like transposable element characterized in Bombyx mori. The very short 5' flanking sequences of A/B.L11 and A/B.L12 (277 and 276 base-pairs) are distinct as shown by hybridization but both recur in additional chorion gene pairs, forming two respective classes that are expressed during distinctly different developmental periods. The divergently transcribed genes of each pair, which border the same 5' flanking sequence, are expressed co-ordinately, during the same developmental period. Detailed comparisons of the 5' flanking regions, and of the corresponding region of the Drosophila s15-1 chorion gene, revealed numerous, very short sequence elements that are shared. One such element, T-C-A-C-G-T, is also associated with all five sequenced Drosophila chorion genes. Some elements are repeated in a dyad symmetrical pattern, i.e. are associated with each of the two genes in a pair, while others, including T-C-A-C-G-T, occur only once per 5' flanking region, and, if functionally important, would presumably act bi-directionally on both genes of the pair.     

Gene evolution and regulation in the chorion complex of Bombyx mori. Hybridization and sequence analysis of multiple developmentally middle A/B chorion gene pairs

Twenty-two pairs of chorion genes belonging to the A and B multigene families have been characterized and mapped within two segments of a 320 kb (1 kb = 10(3) bases or base-pairs) chromosomal walk in the domesticated silkmoth Bombyx mori. Eighteen of the gene pairs belong to two groups that are typified by the previously characterized A/B.L12 and A/B.L11 chorion gene pairs, and are defined by two respective types of short (approx. 280 base-pairs) bidirectional promoter sequences. In the chromosome, the L12-like and L11-like pairs are interspersed with each other and with the remaining four gene pairs, which have unrelated promoter sequences. We have sequenced the promoter regions and adjacent small exons of all L12-like and L11-like A and B genes in the walk. The L12-like promoters are highly conserved, whereas L11-like promoters are somewhat more variable. Reconsideration of previous data on RNA accumulation and disappearance during choriogenesis, in the light of the sequences, indicates that L12-like genes are developmentally early-middle, while L11-like genes correspond to two developmental subgroups, middle I and middle II. Sequence comparisons of all these promoters, as well as the previously characterized promoters of the developmentally late HcA and HcB genes, identify short elements of possible regulatory significance. The sequences, as well as extensive cross-hybridization analysis with short probes derived from the reference A/B.L12 gene pair, under carefully controlled conditions of stringency, indicate the occurrence of sequence transfers among A or B genes. These sequence transfers, which could result from gene conversions or unequal crossovers, are less abundant than in the HcA and HcB families, but do result in a patchwork of similarities and differences in the A and B genes. The transfers appear to be least frequent between the moderately divergent A genes that belong to different temporal classes, while the L12-like and L11-like B genes appear to be extensively homogenized in sequence.     

On the statistical significance of homologous structures among the Escherichia coli ribosomal proteins

Summary Completion of the sequence determination of all 52 Escherichia coli ribosomal proteins enabled a final comparison of their sequences. Similarities in amino acid compositions were compared to the relatedness of the sequences, which was analyzed statistically with the aid of the computer programs RELATE and ALIGN.Among the examined 52×52 possible protein pairs at least 40 pairs were found that can be regarded as distantly related (showing segment comparison score values slightly above 3.0 S.D. units). These protein pairs were further examined with the programs ALIGN and SEEK to locate homologous sequence stretches. In no case were two complete homologous sequences found (with the exception of the known identical pairs L7/L12 and S20/L26). However, short homologous sequence regions were observed. Beside those protein pairs that show significant although distant relatedness, other pairs were slightly below the threshold value of 3.0 S.D. units.Those pairs observed to be distantly related consisted either of two proteins from the same subunit or of one protein from each of the different subunits. A further analysis of these pairs revealed a correlation between their relatedness and their time of incorporation into the ribosome during assembly.     

Mitochondrial DNA Sequence Divergence among Lycopersicon and Related Solanum Species

Sequence divergence among the mitochondrial (mt) DNA of nine Lycopersicon and two closely related Solanum species was estimated using the shared fragment method. A portion of each mt genome was highlighted by probing total DNA with a series of plasmid clones containing mt-specific DNA fragments from Lycopersicon pennellii. A total of 660 fragments were compared. As calculated by the shared fragment method, sequence divergence among the mtDNAs ranged from 0.4% for the L. esculentum-L. esculentum var. cerasiforme pair to 2.7% for the Solanum rickii-L. pimpinellifolium and L. cheesmanii-L. chilense pairs. The mtDNA divergence is higher than that reported for Lycopersicon chloroplast (cp) DNA, which indicates that the DNAs of the two plant organelles are evolving at different rates. The percentages of shared fragments were used to construct a phenogram that illustrates the present-day relationships of the mtDNAs. The mtDNA-derived phenogram places L. hirsutum closer to L. esculentum than taxonomic and cpDNA comparisons. Further, the recent assignment of L. pennellii to the genus Lycopersicon is supported by the mtDNA analysis.     

Molecular analysis of the D-genome of the Triticeae

Summary The chromosome of three tetraploid Aegilops L. species containing the D-genome were analyzed by in situ hybridization with a repeated DNA sequence clone pAS1 isolated from Aegilops squarrosa and observed to be D-genome specific. This sequence is found on all seven D-genome chromosome pairs of A. squarrosa and hexaploid wheat. Two distinct D-genome patterns were observed in the tetraploid species. The D-genome of A. cylindrica was similar to hexaploid wheat. Seven pairs of chromosomes having large amounts and numerous sites of the sequence were observed. Five chromosome pairs with fewer and smaller sites of the repetitive sequence were observed in the D-genomes of A. crassa and A. ventricosa. In addition to these major repeated sequence differences, chromosomal modifications appear to have occurred between T. aestivum and A. cylindrica and between A. crassa and A. ventricosa resulting in changes with respect to location of the sequence between the respective species. D-genome divergence with respect to pAS1 sequence appears to have occurred at least in two forms, one characterized by the changes in amount of repetitive sequence and the second by changes in location of the sequence.     

Types, levels and patterns of low-copy DNA sequence divergence, and phylogenetic implications, for Gossypium genome types

To explore types, levels and patterns of genetic divergence among diploid Gossypium (cotton) genomes, 780 cDNA, genomic DNA and simple sequence repeat (SSR) loci were re-sequenced in Gossypium herbaceum (A1 genome), G. arboreum (A2), G. raimondii (D5), G. trilobum (D8), G. sturtianum (C1) and an outgroup, Gossypioides kirkii. Divergence among these genomes ranged from 7.32 polymorphic base pairs per 100 between G. kirkii and G. herbaceum (A1) to only 1.44 between G. herbaceum (A1) and G. arboreum (A2). SSR loci are least conserved with 12.71 polymorphic base pairs and 3.77 polymorphic sites per 100 base pairs, whereas expressed sequence tags are most conserved with 3.96 polymorphic base pairs and 2.06 sites. SSR loci also exhibit the highest percentage of 'extended polymorphisms' (spanning multiple consecutive nucleotides). The A genome lineage was particularly rapidly evolving, with the D genome also showing accelerated evolution relative to the C genome. Unexpected asymmetry in mutation rates was found, with much more transition than transversion mutation in the D genome after its divergence from a common ancestor shared with the A genome. This large quantity of orthologous DNA sequence strongly supports a phylogeny in which A-C divergence is more recent than A-D divergence, a subject that is of much importance in view of A-D polyploid formation being key to the evolution of the most productive and finest-quality cottons. Loci that are monomorphic within A or D genome types, but polymorphic between genome types, may be of practical importance for identifying locus-specific DNA markers in tetraploid cottons including leading cultivars.     

Mitochondrial <Emphasis Type="Italic">cox</Emphasis>1 and plastid <Emphasis Type="Italic">rbc</Emphasis>L genes of <Emphasis Type="Italic">Gracilaria vermiculophylla</Emphasis> (Gracilariaceae,Rhodophyta)

The mitochondrial cytochrome c oxidase subunit I gene sequence was recently developed for DNA barcoding of red algal species. We determined the 1245 base pairs of the gene from 27 taxa of an agar-producing species, Gracilaria vermiculophylla, and putative relatives and compared the results with rbcL data from the same species. A total of 392 positions (31.5%) were variable, 282 positions (22.6%) were parsimoniously informative, and average sequence divergence was 13% in an ingroup. Within G. vermiculophylla, pairwise divergence of the gene was variable up to 11 bp (0.9%). Seven recognized haplotypes of cox1 tended to be geographically related. In the aligned 1386 bp of rbcL, three haplotypes were recognized. These results suggest that cox1 is a valuable molecular marker within species and will be very useful in haplotype analyses.     

Tempo and mode of concerted evolution in the L1 repeat family of mice

A 300-bp DNA sequence has been determined for 30 (10 from each of three species of mice) random isolates of a subset of the long interspersed repeat family L1. From these data we conclude that members of the L1 family are evolving in concert at the DNA sequence level in Mus domesticus, Mus caroli, and Mus platythrix. The mechanism responsible for this phenomenon may be either duplicative transposition, gene conversion, or a combination of the two. The amount of intraspecies divergence averages 4.4%, although between species base substitutions accumulate at the rate of approximately 0.85%/Myr to a maximum divergence of 9.1% between M. platythrix and both M. domesticus and M. caroli. Parsimony analysis reveals that the M. platythrix L1 family has evolved into a distinct clade in the 10-12 Myr since M. platythrix last shared a common ancestor with M. domesticus and M. caroli. The parsimony tree also provides a means to derive the average half-life of L1 sequences in the genome. The rates of gain and loss of individual copies of L1 were estimated to be approximately equal, such that approximately one-half of them turn over every 3.3 Myr.      

A chloroplast DNA phylogeny of lilacs (Syringa, Oleaceae): plastome groups show a strong correlation with crossing groups

Phylogenetic relationships and genomic compatibility were compared for 60 accessions of Syringa using chloroplast DNA (cpDNA) and nuclear ribosomal DNA (rDNA) markers. A total of 669 cpDNA variants, 653 of which were potentially phylogenetically informative, was detected using 22 restriction enzymes. Phylogenetic analyses reveal four strongly supported plastome groups that correspond to four genetically incompatible crossing groups. Relationships of the four plastome groups (I(II(III,IV))) correlate well with the infrageneric classification except for ser. Syringa and Pinnatifoliae. Group I, which includes subg. Ligustrina, forms a basal lineage within Syringa. Group II includes ser. Syringa and Pinnatifoliae and the two series have high compatibility and low sequence divergence. Group III consists of three well-defined species groups of ser. Pubescentes. Group IV comprises all members of ser. Villosae and has the lowest interspecific cpDNA sequence divergences. Comparison of cpDNA sequence divergence with crossability data indicates that hybrids have not been successfully generated between species with divergence greater than 0.7%. Hybrid barriers are strong among the four major plastome groups, which have sequence divergence estimates ranging from 1.096 to 1.962%. In contrast, fully fertile hybrids occur between species pairs with sequence divergence below 0.4%. Three regions of the plastome have length variants of greater than 100 bp, and these indels identify 12 different plastome types that correlate with phylogenetic trees produced from cpDNA restriction site data. Biparentally inherited nuclear rDNA and maternally inherited cpDNA length variants enable the identification of the specific parentage of several lilac hybrids.     

Molecular evolution of theSaccharomyces cerevisiae histone gene loci

Summary The core histone genes ofSaccharomyces cerevisiae are arranged as duplicate nonallelic sets of specifically paired genes. The identity of structural organization between the duplicated gene pairs would have its simplest evolutionary origin in the duplication of a complete locus in a single event. In such a case, the time since the duplication of one of the genes should be identical to that since duplication of the gene adjacent to it on the chromosome. A calculation of the evolutionary distances between the coding DNA sequences of the histone genes leads to a duplication paradox: The extents of sequence divergence in the silent component of third-base positions for adjacent pairs of genes are not identical. Estimates of the evolutionary distance between the two H3-H4 noncoding intergene DNA sequences are large; the divergence between the two separate sequences is indistinguishable from the divergence between either of the regions and a randomly generated permutation of itself. These results suggest that the duplication event may have occurred much earlier than previously estimated. The potential age of the duplication, and the attractive simplicity of the duplication of both the H3-H4 and the H2A-H2B gene pairs having taken place in a single event, leads to the hypothesis that modern haploidS. cerevisiae may have evolved by diploidization or fusion of two ancient fungi.     

Genetic diversity within human erythroviruses: identification of three genotypes

B19 virus is a human virus belonging to the genus Erythrovirus: The genetic diversity among B19 virus isolates has been reported to be very low, with less than 2% nucleotide divergence in the whole genome sequence. We have previously reported the isolation of a human erythrovirus isolate, termed V9, whose sequence was markedly distinct (>11% nucleotide divergence) from that of B19 virus. To date, the V9 isolate remains the unique representative of a new variant in the genus Erythrovirus, and its taxonomic position is unclear. We report here the isolation of 11 V9-related viruses. A prospective study conducted in France between 1999 and 2001 indicates that V9-related viruses actually circulate at a significant frequency (11.4%) along with B19 viruses. Analysis of the nearly full-length genome sequence of one V9-related isolate (D91.1) indicates that the D91.1 sequence clusters together with but is notably distant from the V9 sequence (5.3% divergence) and is distantly related to B19 virus sequences (13.8 to 14.2% divergence). Additional phylogenetic analysis of partial sequences from the V9-related isolates combined with erythrovirus sequences available in GenBank indicates that the erythrovirus group is more diverse than thought previously and can be divided into three well-individualized genotypes, with B19 viruses corresponding to genotype 1 and V9-related viruses being distributed into genotypes 2 and 3.     

Recognition specificity of self-incompatibility maintained after the divergence of Brassica oleracea and Brassica rapa

The determinants of recognition specificity of self-incompatibility in Brassica are SRK in the stigma and SP11/SCR in the pollen, respectively. In the pair of S haplotypes BrS46 (S46 in B. rapa) and BoS7 (S7 in B. oleracea), which have highly similar SRK alleles, the SP11 alleles were found to be similar, with 96.1% identity in the deduced amino acid sequence. Two other pairs of S haplotypes, BrS47 and BoS12, and BrS8 and BoS32, having highly similar SRK and SP11 alleles between the two species were also found. The haplotypes in each pair are considered to have been derived from a single S haplotype in the ancestral species. The allotetraploid produced by interspecific hybridization between homozygotes of BrS46 and BoS15 showed incompatibility with a BoS7 homozygote and compatibility with other B. oleracea S haplotypes in reciprocal crossings. This result indicates that BrS46 and BoS7 have maintained the same recognition specificity after the divergence of the two species and that amino acid substitutions found in such cases in both SRK alleles and SP11 alleles do not alter the recognition specificity. DNA blot analysis of SRK, SP11, SLG and other S-locus genes showed different DNA fragment sizes between the interspecific pairs of S haplotypes. A much lower level of sequence similarity was observed outside the genes of SRK and SP11 between BrS46 and BoS7. These results suggest that the DNA sequences of the regions intervening between the S-locus genes were diversified after or at the time of speciation. This is the first report demonstrating the presence of common S haplotypes in different plant species and presenting definite evidence of the trans-specific evolution of self-incompatibility genes.     

Evolutionary divergence of mitochondrial DNA from Paramecium aurelia

Summary Mitochondrial (mt) DNA from four sibling species within the Paramecium aurelia complex, including stocks of different geographic origin and mutants, were analyzed using four 6-bp recognition site and one 4-bp recognition site endonucleases and the sequence divergence was estimated using Upholt's (1977) statistical procedure. All four species were readily distinguishable regardless of the restriction endonuclease employed. With intraspecies comparisons, no differences were observed which could be accounted for on the basis of geographic origin. Except for species 4, each stock (and mutant) gave a species-specific fragment pattern. For species 4, while the patterns were distinct from the other species, two species-specific type of patterns were found, designated A and B. The sequence divergence between these was estimated to be between 1 and 2 percent. With interspecies comparisons, the sequence divergence ranged from 3.9 to 10.3% with the greatest divergence being between species 1 and 4, and the least between species 1 and 5. The similarity between species 1 and 5 is in accord with other criteria for interspecies comparisons. The degree of sequence divergence measured here in Paramecium mt DNA is well within the range reported for rodents and primates. All four species mt DNA were cleaved to many DNA fragments by DPN II, an enzyme which recognizes non-methylated sites, and not by DPNI, the methyl-site specific counterpart of DPN II, suggesting that mt DNA from Paramecium aurelia is not appreciably methylated, if at all.     

Gene Duplication and Gene Conversion in the Caenorhabditis elegans Genome

A comprehensive analysis of duplication and gene conversion for 7394 Caenorhabditis elegans genes (about half the expected total for the genome) is presented. Of the genes examined, 40% are involved in duplicated gene pairs. Intrachromosomal or cis gene duplications occur approximately two times more often than expected. In general the closer the members of duplicated gene pairs are, the more likely it is that gene orientation is conserved. Gene conversion events are detectable between only 2% of the duplicated pairs. Even given the excesses of cis duplications, there is an excess of gene conversion events between cis duplicated pairs on every chromosome except the X chromosome. The relative rates of cis and trans gene conversion and the negative correlation between conversion frequency and DNA sequence divergence for unconverted regions of converted pairs are consistent with previous experimental studies in yeast. Three recent, regional duplications, each spanning three genes are described. All three have already undergone substantial deletions spanning hundreds of base pairs. The relative rates of duplication and deletion may contribute to the compactness of the C. elegans genome. Received: 30 July 1998 / Accepted: 12 October 1998     

Single copy DNA and structural gene sequence relationships among four sea urchin species

Measurements of the divergence of single copy DNA sequences among four sea urchin species are presented. At a standard criterion for reassociation (0.12 M phosphate buffer, 60° C, hydroxyapatite binding) we observe the following extents of reaction and reductions in thermal stability for single copy DNA reassociation between Strongylocentrotus purpuratus tracer and heterologous driver DNA: S. dröbachiensis 68% and 2.5°C; S. franciscanus 51% and 3.5° C; Lytechinus pictus 12% and 7.5° C. The implied extents of sequence relatedness are consistent with the phylogenetic relationships of these species. The rate of single copy sequence divergence in the evolutionary lines leading to the Strongylocentrotus species is estimated to be 0.06–0.35% per million years. The rate of divergence of total single copy sequence has been compared to that of structural gene sequences represented in S. purpuratus gastrula polysomal messenger RNA. When closely related species, S. purpuratus and S. franciscanus, are compared, these polysomal sequences are found to diverge at a lower rate than does the total single copy sequence. For two very distantly related species, S. purpuratus and L. pictus, a small fraction of the single copy DNA sequence is probably conserved. These conserved sequences are not enriched in their content of structural gene sequences.Also staff member, Carnegie Institution of Washington, Washington, D.C. 20015     

ITS-1 ribosomal DNA sequence variants are maintained in different species and strains of Echinococcus

This study investigated sequence heterogeneity in the first internal transcribed spacer (ITS-1) of ribosomal DNA within and among species and strains of Echinococcus. Different ITS-1 sequence variants exist in Echinococcus granulosus and Echinococcus multilocularis, which represent at least four evolutionary lineages: (1) a sheep strain-lineage of E. granulosus, (2) a sister lineage of a cervid and camel E. granulosus ITS-1 variants, (3) a lineage including the ITS-1 variants representing horse, bovine and camel strains of E. granulosus, as well as variants from E. multilocularis, Echinococcus oligarthrus and Echinococcus vogeli and (4) a distinct lineage of ITS-1 variants including E. granulosus strains from sheep and cervid, and E. multilocularis. At least two of the species (E. granulosus and E. multilocularis) were paraphyletic for ITS-1. Divergent ITS-1 variants from these two species shared distinct evolutionary lineages. The sequence data provided evidence that at least two turnover mechanisms, namely slippage and unequal crossing over/transposition, have led to the divergence and maintenance of sequence variants in Echinococcus species and strains.     

Features of two hepatitis B virus (HBV) DNA integrations suggest mechanisms of HBV integration.

Two integrated hepatitis B virus (HBV) DNA molecules were cloned from two primary hepatocellular carcinomas each containing only a single integration. One integration (C3) contained a single linear segment of HBV DNA, and the other integration (C4) contained a large inverted duplication of viral DNA at the site of a chromosome translocation (O. Hino, T.B. Shows, and C.E. Rogler, Proc. Natl. Acad. Sci. USA 83:8338-8342, 1986). Sequence analysis of the virus-cell junctions of C3 placed the left virus-cell junction at nucleotide 1824, which is at the 5' end of the directly repeated DR1 sequence and is 6 base pairs from the 3' end of the long (L) negative strand. The right virus-cell junction was at nucleotide 1762 in a region of viral DNA (within the cohesive overlap) which shared 5-base-pair homology with cellular DNA. Sequence analysis of the normal cellular DNA across the integration site showed that 11 base pairs of cellular DNA were deleted at the site of integration. On the basis of this analysis, we suggest a mechanism for integration of the viral DNA molecule which involves strand invasion of the 3' end of the L negative strand of an open circular or linear HBV DNA molecule (at the DR1 sequence) and base pairing of the opposite end of the molecule with cellular DNA, accompanied by the deletion of 11 base pairs of cellular DNA during the double recombination event. Sequencing across the inverted duplication of HBV DNA in clone C4 located one side of the inversion at nucleotide 1820, which is 2 base pairs from the 3' end of the L negative strand. Both this sequence and the left virus-cell junction of C3 are within the 9-nucleotide terminally redundant region of the HBV L negative strand DNA. We suggest that the terminal redundancy is a preferred topoisomerase I nicking region because of both its base sequence and forked structure. Such nicking would lead to integration and rearrangement of HBV molecules within the terminal redundancy, as we have observed in both our clones.     

BIOGEOGRAPHYOF BATRACHOSPERMUM GELATINOSUM (BATRACHOSPERMALES,RHODOPHYTA) IN NORTH AMERICA BASED ON MOLECULAR AND MORPHOLOGICAL DATA1

Fifteen populations of the widespread fieshwater red alga Batrachospermum gelatinosum (L.) De Candolle were sampled throughout the geographic range in North America from central Alabama, U.S.A. (33° N), to Ellesmere Island, Northwest Territories (NWT), Canada (80° N). Analysis of ribosomal DNA internal transcribed spacer (ITS) 1 and 2 sequences yielded a parsimony tree with a large polytomy consisting of most populations plus a branch with one Nova Scotia and two NWT populations. The nucleotide variation, both within the polytomy and within the branch, was small (< 1%). The sequence divergence between the branch and polytomy was 3%. The lengths of the ITS 1 and 2 sequences of B. gelatinosum, 216–229 and 448–458 base pairs, respectively, fall within the very broad ranges reported for other red algae. The cluster analysis of 11 morphometric characteristics revealed three groupings of populations, partly based on geographic distribution. All tundra, eastern boreal forest, and mid-western hemlock-hardwood populations were in one grouping, whereas the deciduous forest, coastal plain, and eastern hemlock-hardwood populations were in a second. How ever, one deciduous forest population from Rhode Island, U.S.A. was unassociated. There was considerable overlap in morphometric characteristics among the three groupings. Based on this fact and the relatively small nucleotade variation in ITS sequences, we conclude that B. gelatinosum is a morphologically variable and geographically widespread species that is a valid taxonomic entity.