Loss of phylogenetic information in chorion gene families of Bombyx mori gene conversion

The silkmoth chorion has provided a stimulating model for the study of evolution and developmental regulation of gene families. Previous attempts at inferring relationships among chorion sequences have been based on pairwise comparisons of overall similarity, a potentially problematic approach. To remedy this, we identified the alignable regions of low sequence variability and then analyzed this restricted database by parsimony and neighbor-joining methods. At the deepest level, the chorion sequence tree is split into two branches, called "alpha" and "beta." Within each branch, early- and late-expressing genes each constitute monophyletic groups, while the situation with middle-expressing genes remains uncertain. The HcB gene family appears to be the most basal beta-branch group, but this conclusion is qualified because the effect of gene conversion on branching order is unknown. Previous studies by Eickbush and colleagues have strongly suggested that ErA, HcA, and HcB families undergo gene conversion within a gene family, whereas the ErB family does not. The occurrence of conversion correlates with a particular tree structure; namely, branch lengths are much greater at the base of the family than at higher internodes and terminal branches. These observations raise the possibility that chorion gene families are defined by gene conversion events (reticulate evolution) rather than by descent with modification (synapomorphy).      

Evolution of two major chorion multigene families as inferred from cloned cDNA and protein sequences

Complete or partial sequences are reported from six chorion cDNA clones of the silkmoth Antheraea polyphemus. The proteins encoded belong to the two major chorion protein classes, A and B, each of which is encoded by a multigene family. The sequence comparisons define some major features of the families and suggest how these genes may be evolving. Deletions and insertions might be involved in expanding or contracting internally repetitive regions. Sequence divergence is localized, thus defining sequence domains of distinct evolutionary properties and presumably distinct functions.     

The silkmoth late chorion locus. I. Variation within two paired multigene families

The 140 X 10(3) base late chorion locus of Bombyx mori contains two 15-member multigene families arranged in tightly linked pairs, which are divergently transcribed (the high-cysteine A (HcA) and the high-cysteine B (HcB) families). Previous DNA hybridization experiments have indicated that all members of these gene families contain a complex pattern of shared sequence variation. The sequence analysis in this paper involving all 15 gene pairs allows a comprehensive examination of the nature of this variation. Average sequence homology between gene pairs is: 95% for the protein-encoding regions; 93% for the common 272 base-pair 5' flanking region; 87% for the introns; and 88% for the 3' untranslated regions. Considering the great degree of sequence homology in the coding regions, an unexpectedly high level of variation is found in the deduced protein sequences. Over 50% of the nucleotide substitutions in the protein-encoding regions lead to amino acid replacements, most of which involve a change in charge or effect the secondary structure of the protein. In addition, significant differences in length between the proteins occur in the carboxyl-terminal arm. In both families, the major portion of this arm is composed of Cys-Gly-Gly and Cys-Gly subrepeats forming a (Cys-Gly-Gly)2-(Cys-Gly)2 major repeat. Differences in the number of complete and partial repeats results in deduced protein sequences that contain arms varying from 32 to 54 amino acid residues for members of the HcA family and 14 to 88 residues for the HcB family. The high level of variation in protein composition indicates a lack of strong selective pressure. We suggest the high level of DNA sequence homology maintained by these genes in the coding as well as in the non-coding regions is the result of sequence exchange between family members.     

Gene conversions within the skeletal myosin multigene family.

Comparisons of the nucleotide sequences of the light meromyosin (LMM) region of developmentally regulated fast chicken myosin heavy chain (MHC) isoforms indicates that chicken MHC isoforms are more similar to each other than to MHC isoforms in other species. The sequence data provide evidence that gene conversion events have occurred recently among the isoforms. An embryonic (Cemb1) isoform and neonatal isoform have the most extensive regions of sequence identity. Similar gene conversion events are present in the rat alpha- and beta-cardiac MHCs, but were not obvious in the LMM of developmentally regulated fast human MHC isoforms. The data suggest that gene conversion events can play a significant role in the evolution of the MHC multigene families and that concerted evolution of the chicken multigene family occurred after the divergence of mammals and avians.     

The biolistic method as a tool for testing the differential activity of putative silkmoth chorion gene promoters

Bombyx mori unpaired early chorion gene copies 6F6.1,.2 and.3 are exceptions to the typical organization and distribution pattern of known early ErA/ErB, middle A/B and late HcA/HcB divergently transcribed gene pairs. Contrary to such pairs, the boundaries of the 6F6 regulatory sequences are not easily defined; moreover, they share common sequence elements with the regulatory sequences of middle and late genes. In order to perform a functional study of the tissue and temporal specificity of the 6F6 putative promoter region, we decided to apply biolistics. In the present work, use of a region from the 6F6.2 5' untranslated sequence, spanning nucleotides -138 to the cap site, gave an expected expression pattern of a lacZ reporter gene. Temporal specificity was further verified by control experiments using the cloned intergenic sequence of the late gene pair HcA/B.12, which resulted in lacZ expression in late choriogenic follicles. At present, despite the recent successful germinal transgenesis of Bombyx mori, the biolistic transient expression system seems to be the most rapid technique to pursue the functional study of the promoter region of early chorion genes, including the three unconventional early 6F6 genes.     

Diversity in a chorion multigene family created by tandem duplications and a putative gene-conversion event

Summary Two families of high-cysteine chorion proteins inBombyx mori are encoded in 15 tandemly arranged nonidentical gene pairs. It is assumed that this locus arose by duplication with subsequent sequence divergence. We have compared DNA sequences from two such neighboring pairs of genes in an attempt to understand the manner in which diversity has been generated and/or removed. A high level of sequence identity (91%–99%) was found between the repeats throughout the transcribed and flanking regions, with two significant exceptions. First, in the DNA segment encoding a conserved region of the chorion proteins, ten substitutions were detected in a 39-base-pair region. This localized region of high variability would suggest an intergene conversion-like event. Second, a length difference of 141 base pairs was detected in a region encoding the carboxy-terminal arm of the protein. This difference can be explained by three separate reiterations of single codons (3 base pairs) separated in time by duplication or triplication events.     

Population genetics theory of concerted evolution and its application to the immunoglobulin V gene tree

Summary The previous simple model for treating concerted evolution of multigene families has been revised to be compatible with various new observations on the immunoglobulin variable region family and other families. In the previous model, gene conversion and unequal crossing-over were considered, and it was assumed that genes are randomly arranged on the chromosome; neither subdivision nor correlation of gene identity and chromosomal distance were considered. Although this model satisfactorily explains the observed amino acid diversity within and between species, it fails to predict the very ancient branching of the mouse immunoglobulin heavy chain V-gene family. By incorporating subdivided structure and genetic correlation with chromosomal distance into the simple model, the data of divergence may be satisfactorily explained, as well as the rate of nucleotide substitution and the amino acid diversity. The rate at which a V-gene is duplicated or deleted by conversion or by unequal crossing-over is estimated by the new model to be on the order of 10^–6 per year. The model may be applicable to other multigene families, such as those coding for silkmoth chorion or mammalian kallikrein.Contribution no. 1560 from the National Institute of Genetics, Mishima, 411 Japan     

Polymorphism and Gene Conversion in Mouse α-Globin Haplotypes

We have cloned and characterized three distinct alpha-globin haplotypes obtained from inbred strains of the mouse, Mus domesticus. We report here the complete nucleotide sequence of the six alpha-globin genes that the haplotypes contain. Our analysis of these genes and those from one other previously described haplotype indicates that recurrent gene conversion events have played a major role in their history. The pattern of nucleotide substitutions suggests that conversions have occurred both within and between haplotypes. Limited segments of coding and noncoding DNA have been involved in these gene conversion events. In two of the haplotypes, the nonallelic genes of each maintain DNA sequence identity over discrete intervals and encode the same alpha-globin polypeptide. On the other hand, the coding regions of some genes have accumulated replacement changes that result in distinct alpha-globins. In one instance, these changes appear to reflect positive selection of advantageous mutations.     

Comparison of Bombyx mori and Helicoverpa armigera cytoplasmic actin genes provides clues to the evolution of actin genes in insects

The cytoplasmic actin genes BmA3 and BmA4 of Bombyx mori were found clustered in a single genomic clone in the same orientation. As a similar clustering of the two cytoplasmic actin genes Ha3a and Ha3b also occurs in another lepidopteran, Helicoverpa armigera, we analyzed the sequence of the pair of genes from each species. Due to the high conservation of cytoplasmic actins, the coding sequence of the four genes was easily aligned, allowing the detection of similarities in noncoding exon and intron sequences as well as in flanking sequences. All four genes exhibited a conserved intron inserted in codon 117, an original position not encountered in other species. It can thus be postulated that all of these genes derived from a common ancestral gene carrying this intron after a single event of insertion. The comparison of the four genes revealed that the genes of B. mori and H. armigera are related in two different ways: the coding sequence and the intron that interrupts it are more similar between paralogous genes within each species than between orthologous genes of the two species. In contrast, the other (noncoding) regions exhibited the greatest similarity between a gene of one species and a gene of the other species, defining two pairs of orthologous genes, BmA3 and HaA3a on one hand and BmA4 and HaA3b on the other. However, in each species, the very high similarities of the coding sequence and of the single intron that interrupts it strongly suggest that gene conversion events have homogenized this part of the sequence. As the divergence of the B. mori genes was higher than that of the H. armigera genes, we postulated that the gene conversion occurred earlier in the B. mori lineage. This leads us to hypothesize that gene conversion could also be responsible for the original transfer of the common intron to the second gene copy before the divergence of the B. mori and H. armigera lineages.     

Molecular Evolution of Drosophila Cuticular Protein Genes

Several multigene families have been described that together encode scores of structural cuticular proteins in Drosophila, although the functional significance of this diversity remains to be explored. Here I investigate the evolutionary histories of several multigene families (CPR, Tweedle, CPLCG, and CPF/CPFL) that vary in age, size, and sequence complexity, using sequenced Drosophila genomes and mosquito outgroups. My objective is to describe the rates and mechanisms of ‘cuticle-ome’ divergence, in order to identify conserved and rapidly evolving elements. I also investigate potential examples of interlocus gene conversion and concerted evolution within these families during Drosophila evolution. The absolute rate of change in gene number (per million years) is an order of magnitude lower for cuticular protein families within Drosophila than it is among Drosophila and the two mosquito taxa, implying that major transitions in the cuticle proteome have occurred at higher taxonomic levels. Several hotspots of intergenic conversion and/or gene turnover were identified, e.g. some gene pairs have independently undergone intergenic conversion within different lineages. Some gene conversion hotspots were characterized by conversion tracts initiating near nucleotide repeats within coding regions, and similar repeats were found within concertedly evolving cuticular protein genes in Anopheles gambiae. Rates of amino-acid substitution were generally severalfold higher along the branch connecting the Sophophora and Drosophila species groups, and 13 genes have Ka/Ks significantly greater than one along this branch, indicating adaptive divergence. Insect cuticular proteins appear to be a source of adaptive evolution within genera and, at higher taxonomic levels, subject to periods of gene-family expansion and contraction followed by quiescence. However, this relative stasis is belied by hotspots of molecular evolution, particularly concerted evolution, during the diversification of Drosophila. The prominent association between interlocus gene conversion and repeats within the coding sequence of interacting genes suggests that the latter promote strand exchange.     

Patterns of Sequence Divergence in Daphniid Hemoglobin Genes

In common with other multigene families, sequence diversity in the hemoglobin genes of cladoceran crustaceans has been heavily impacted by gene conversion events. Because of their structural complexity (six exons, five introns), these genes provide a good opportunity to study the influence of intron length and position on the conversion process. This study surveys the patterns of divergence in variants of one hemoglobin gene (H1) from two closely related species of Daphnia using a PCR-based approach. Although its effects were most pronounced at their 5' ends, intron and exon regions of these genes showed similar exposure to gene conversion, excepting intron 2. This intron, which was the only one with a marked length difference among variants, showed substantial sequence divergence, suggesting that gene conversion was disrupted. These results, together with those on hemoglobin gene families in other organisms, indicate that sequence tracts showing gene conversion are often distributed in a mosaic fashion. The reactivation of gene conversion downstream of a block protected from its effects suggests that there are multiple initiation points, and the distribution of conversion tracts suggests that exon/intron splice sites are important in this regard.     

Gene Conversions Can Generate Sequence Variants in the Late Chorion Multigene Families of Bombyx Mori

The 140-kbp late chorion locus of Bombyx mori strain 703 contains 15 divergently oriented gene pairs encoding the high cysteine (Hc) eggshell proteins. Sequence homology is approximately 91% for the 2-kb region of each gene pair, including the 5' flanking region, intron and exons. The homology rapidly disappears within a few hundred basepairs of the 3' end of most genes. Here we present the results of the nucleotide sequence and genomic blot comparison of Hc genes from different races of B. mori. Comparison of the nucleotide sequences of the same gene pair in two different races reveals that most of the nucleotide differences occur in clusters or patches and correspond to sequences present in other Hc genes in the locus. The number of nucleotide differences that have accumulated in the highly conserved regions of the gene pair (2.3/100 bp), most of which are attributable to patchwork exchanges, is significantly higher than the number of differences in the poorly conserved 3' flanking regions (0.6/100 bp), due primarily to new mutations. These data are consistent with a gene conversion process, which in the short-term generates new combinations of sequence variants, but in the long-term results in concerted evolution. Genomic blot analyses of different geographical races of B. mori reveal that there is variation in the number of Hc gene pairs (14-19 gene pairs), indicating that unequal crossovers also occur in the locus.     

Independent gene evolution in the potato actin gene family demonstrated by phylogenetic procedures for resolving gene conversions and the phylogeny of angiosperm actin genes

Summary Nine different actin DNA sequences were isolated from the common potato,Solanum tuberosum, and the nucleotide sequence of five actin loci and of two allelic variants are presented. Unlike the wide variation in intron position among animal actin genes, the potato actin genes have three introns situated in the same positions as reported for all other angiosperm actin genes. Using a novel combination of analytical procedures (G-test and compatibility analysis), we could not find evidence of frequent large or small nonreciprocal exchanges of genetic material between the sequenced loci, although there were a few candidates. Resolution of such gene conversion events and the quantification of independence of gene evolution in multigene families is critical to the inference of phylogenetic relationships. Comparison with actin genes in other angiosperm species suggests that the actin multigene family can be divided into a number of subfamilies, evolved by descent rather than gene conversion, which are of possible functional origin, with one major subfamily diversification occurring before the divergence of monocots and dicots. The silent rate of nucleotide substitution was estimated to be similar to that suggested for a number of other plant nuclear genes, whereas the replacement rate was extremely slow, suggestive of selective constraints.     

The keratin BIIIB gene family: isolation of cDNA clones and structure of a gene and a related pseudogene

The nucleotide sequence of cDNA clones encoding the three major BIIIB high-sulfur wool keratin proteins (BIIIB2, 3, and 4) and the structure of a BIIIB4 gene and a BIIIB3 pseudogene are reported. Although Southern blot analysis indicates that the BIIIB genes comprise a multigene family in the sheep genome, they are poorly represented in genomic DNA libraries. The family sequence homology of the coding region extends into the 5' and 3' untranslated regions and the near 5' flanking region of the BIIIB3 and 4 genes. These homologies suggest that the BIIIB3 and 4 genes represent the latest gene duplication event in the evolution of the BIIIB multigene family. Like the genes coding for other wool keratin matrix protein components, the BIIIB genes have the conserved 18-bp sequence immediately 5' to the initiation codon and also appear to lack introns.     

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.     

A walk in the chorion locus of bombyx mori

We have recovered overlapping clones that represent in the aggregate a contiguous segment of chromosomal DNA 270 kb in length, or probably one third of the chorion locus of Bombyx mori. Approximately 70 genes have been identified, the majority of which are arranged in coordinately expressed pairs. The nonidentical genes expressed in the late period of choriogenesis are clustered within a single, 130 kb region, which is flanked by regions containing genes that are active during the middle developmental period. The late genes encode two families of high-cysteine proteins; the evolutionarily persistent clustering of these families contrasts sharply with the extensive sequence diversification of the structural genes and their flanking DNA elements. We discuss the possible regulatory significance of the clustered arrangement, as well as certain features of multigene family evolution.