Further characterization of the cathepsin L-associated protein and its gene in two species of the brine shrimp, Artemia

The major cysteine protease in embryos and larvae of the brine shrimp Artemia franciscana is a heterodimer composed of a cathepsin L-like polypeptide of 28.5 kDa and a 31.5 kDa polypeptide called the cathepsin L-associated protein or CLAP. In a previous study, CLAP was shown to be a cell adhesion protein containing two Fas I domains and two GTP/ATP binding sites known as Walker A and B motifs. Here, we have characterized CLAP and its genes to better understand the role of this protein in Artemia development. The polymerase chain reaction was used to investigate the structure of the CLAP gene in two species of Artemia, the New World bisexual diploid A. franciscana and the Old World parthenogenetic tetraploid Artemia parthenogenetica. The protein coding region of the CLAP gene from each species was 99.5% identical for a protein of 332 amino acids, while the 3' non-coding region, representing nearly 45% of the gene, was only 86% identical between the two related species. However, while the CLAP gene is intronless in A. franciscana, in A. parthenogenetica the gene contained a mini-intron of 30 base pairs in the 3' non-coding region. The sequences representing the CLAP gene in A. franciscana and A. parthenogenetica have been entered into the NCBI database as AY757920 and DQ100385, respectively. Northern blot analysis showed that while the cathepsin L gene is expressed constitutively in Artemia franciscana embryos and young larvae, the CLAP gene is not expressed in late embryos and young larvae. In contrast, Western blots indicated that CLAP is present in developing embryos and young larvae, at least to the first larval molt, supporting results obtained previously showing CLAP's resistance to degradation by its dimeric partner, cathepsin L. At the protein level we showed that the GTP/ATP binding sites in CLAP are functional with rate constants of 0.024 and 0.022 for GTP and ATP hydrolase activity, respectively. GTP but not ATP also had a slight stimulatory effect on cathepsin L activity of the heterodimeric protease containing CLAP. Our results support the hypothesis that CLAP plays an important role in targeting and expression regulation of cathepsin L activity during early development of Artemia.     

Molecular evolution of cadherin-related neuronal receptor/protocadherin(alpha) (CNR/Pcdh(alpha)) gene cluster in Mus musculus subspecies

The mouse cadherin-related neuronal receptor/protocadherin (CNR/Pcdh) gene clusters are located on chromosome 18. We sequenced single-nucleotide polymorphisms (SNPs) of the CNR/Pcdh(alpha)-coding region among 12 wild-derived and four laboratory strains; these included the four major subspecies groups of Mus musculus: domesticus, musculus, castaneus, and bactrianus. We detected 883 coding SNPs (cSNPs) in the CNR/Pcdh(alpha) variable exons and three in the constant exons. Among all the cSNPs, 586 synonymous (silent) and 297 nonsynonymous (amino acid exchanged) substitutions were found; therefore, the K(a)/K(s) ratio (nonsynonymous substitutions per synonymous substitution) was 0.51. The synonymous cSNPs were relatively concentrated in the first and fifth extracellular cadherin domain-encoding regions (ECs) of CNR/Pcdh(alpha). These regions have high nucleotide homology among the CNR/Pcdh(alpha) paralogs, suggesting that gene conversion events in synonymous and homologous regions of the CNR/Pcdh(alpha) cluster are related to the generation of cSNPs. A phylogenetic analysis revealed gene conversion events in the EC1 and EC5 regions. Assuming that the common sequences between rat and mouse are ancestral, the GC content of the third codon position has increased in the EC1 and EC5 regions, although biased substitutions from GC to AT were detected in all the codon positions. In addition, nonsynonymous substitutions were extremely high (11 of 13, K(a)/K(s) ratio 5.5) in the laboratory mouse strains. The artificial environment of laboratory mice may allow positive selection for nonsynonymous amino acid variations in CNR/Pcdh(alpha) during inbreeding. In this study, we analyzed the direction of cSNP generation, and concluded that subspecies-specific nucleotide substitutions and region-restricted gene conversion events may have contributed to the generation of genetic variations in the CNR/Pcdh genes within and between species.     

Hypervariable noncoding sequences in Saccharomyces cerevisiae

Compared to protein-coding sequences, the evolution of noncoding sequences and the selective constraints placed on these sequences is not well characterized. To compare the evolution of coding and noncoding sequences, we have conducted a survey for DNA polymorphism at five randomly chosen loci among a diverse collection of 81 strains of Saccharomyces cerevisiae. Average rates of both polymorphism and divergence are 40% lower at noncoding sites and 90% lower at nonsynonymous sites in comparison to synonymous sites. Although noncoding and coding sequences show substantial variability in ratios of polymorphism to divergence, two of the loci, MLS1 and PDR10, show a higher rate of polymorphism at noncoding compared to synonymous sites. The high rate of polymorphism is not accompanied by a high rate of divergence and is limited to a few small regions. These hypervariable regions include sites with three segregating bases at a single site and adjacent polymorphic sites. We show that this clustering of polymorphic sites is significantly greater than one would expect on the basis of the spacing between polymorphic fourfold degenerate sites. Although hypervariable noncoding sequences could result from selection on regulatory mutations, they could also result from transient mutational hotspots.     

Site-directed point mutations in embryonic stem cells: a gene-targeting tag-and-exchange strategy.

Sequential gene targeting was used to introduce point mutations into one alpha 2 isoform Na,K-ATPase homolog in mouse embryonic stem (ES) cells. In the first round of targeted replacement, the gene was tagged with selectable markers by insertion of a Neor/HSV-tk gene cassette, and this event was selected for by gain of neomycin (G418) resistance. In the second targeted replacement event, the tagged genomic sequence was exchanged with a vector consisting of homologous genomic sequences carrying five site-directed nucleotide substitutions. Embryonic stem cell clones modified by exchange with the mutation vector were selected for loss of the HSV-tk gene by resistance to ganciclovir. Candidate clones were further screened and identified by polymerase chain reaction and Southern blot analysis. By this strategy, the endogenous alpha 2 isoform Na,K-ATPase gene was altered to encode two other amino acids so that the enzyme is resistant to inhibition by cardiac glycosides while maintaining its transmembrane ion-pumping function. Since the initial tagging event and the subsequent mutation-exchange event are independent of one another, a tagged cell line can be used to generate a variety of mutant lines by exchange with various mutation vectors at the tagged locus. This method should be useful for testing specific mutations introduced into the genomes of tissue culture cells and animals and for developing animal models encompassing the mutational variability of known genetic disorders.     

DNA Sequence Evolution of the Amylase Multigene Family in Drosophila Pseudoobscura

The alpha-Amylase locus in Drosophila pseudoobscura is a multigene family of one, two or three copies on the third chromosome. The nucleotide sequences of the three Amylase genes from a single chromosome of D. pseudoobscura are presented. The three Amylase genes differ at about 0.5% of their nucleotides. Each gene has a putative intron of 71 (Amy1) or 81 (Amy2 and Amy3) bp. In contrast, Drosophila melanogaster Amylase genes do not have an intron. The functional Amy1 gene of D. pseudoobscura differs from the Amy-p1 gene of D. melanogaster at an estimated 13.3% of the 1482 nucleotides in the coding region. The estimated rate of synonymous substitutions is 0.398 +/- 0.043, and the estimated rate of nonsynonymous substitutions is 0.068 +/- 0.008. From the sequence data we infer that Amy2 and Amy3 are more closely related to each other than either is to Amy1. From the pattern of nucleotide substitutions we reason that there is selection against synonymous substitutions within the Amy1 sequence; that there is selection against nonsynonymous substitutions within the Amy2 sequence, or that Amy2 has recently undergone a gene conversion with Amy1; and that Amy3 is nonfunctional and subject to random genetic drift.     

Cloning of a cDNA encoding an Artemia franciscana Na/K ATPase alpha-subunit.

Clones of cDNA that code for an isoform of the Artemia franciscana Na/K ATPase alpha subunit (NaKA alpha) have been isolated. The sequence of the longest of these clones (pArATNa136) is 3595 nucleotides; it codes for a 1004-amino acid protein whose sequence is identical to that of two previously sequenced Artemia NaKA alpha peptides. The encoded protein is over 73% identical to Drosophila melanogaster and vertebrate NaKA alpha s, and 73.8% identical to another Artemia NaKA alpha isoform previously described (named alpha 2850 in this article). The two Artemia cDNA clones code for mRNAs of different size; the clone pArATNa136 codes for a 4.5-kb mRNA while the alpha 2850 clone codes for a 3.6-kb mRNA. The degree of homology and the different size of the mRNAs encoded by both cDNAs suggest that they code for two different isoforms of the protein.     

Micro-scale signature of purifying selection in Marburg virus genomes

In the seven protein-coding genes in the Marburg virus (MARV) genome, the synonymous nucleotide diversity substantially exceeded the nonsynonymous nucleotide diversity, indicating strong purifying selection. Likewise, there was evidence of purifying selection on 5'UTR and 3'UTR, where nucleotide diversity (pi) was significantly less than piS in the coding regions. Nonsynonymous polymorphic sites showed significantly reduced mean gene diversity in comparison to other polymorphic sites, indicating that purifying selection at certain slightly deleterious nonsynonymous polymorphisms is ongoing. Moreover, nonsynonymous polymorphic sites showed significantly reduced gene diversity in comparison to adjacent synonymous sites, even though the vast majority of such adjacent synonymous sites were in the same codon or an adjacent codon. Thus purifying selection, in conjunction with recombination and/or backward mutation, can act to break up linkage relationships at a micro-scale in the MARV genome. The ability of purifying selection to break up linkage between synonymous and nonsynonymous polymorphisms on such a fine scale has not been reported in any other genome.     

Expression of an ouabain-resistant Na,K-ATPase in CV-1 cells after transfection with a cDNA encoding the rat Na,K-ATPase alpha 1 subunit

We have used a gene transfer system to investigate the relationship between expression of the rat Na,K-ATPase alpha 1 subunit gene and ouabain-resistant Na,K-ATPase activity. A cDNA clone encoding the entire rat Na,K-ATPase alpha 1 subunit was inserted into the expression vector pSV2neo. This construct (pSV2 alpha 1) conferred resistance to 100 microM ouabain to ouabain-sensitive CV-1 cells. Hybridization analysis of transfected clones revealed the presence of both rat-specific and endogenous Na,K-ATPase alpha 1 subunit DNA and mRNA sequences. A single form of highly ouabain-sensitive 86Rb+ uptake was detected in CV-1 cells, whereas two distinct classes of ouabain-inhibitable uptake were observed in transfectants. One class exhibited the high ouabain sensitivity of the endogenous monkey Na,K-ATPase, while the second class showed the reduced ouabain sensitivity characteristic of the rodent renal Na,K-ATPase. Examination of the ouabain-sensitive, sodium-dependent ATPase activity of the transfectants also revealed a low affinity component of Na,K-ATPase activity characteristic of the rodent kidney enzyme. These results suggest that expression of the rat alpha 1 subunit gene is directly responsible for ouabain-resistant Na,K-ATPase activity in transfected CV-1 cells.     

Evolution of the Sry genes

Existing DNA sequence data on the Sry gene, the mammalian sex- determining locus in the Y chromosome, were analyzed for primates, rodents, and bovids. In all three taxonomic groups, the terminal sequences evolved faster than the HMG (high mobility group) boxes, and this applies both to synonymous (Ks) and nonsynonymous (Ka) nucleotide substitutions. Similar intragenic correlation between synonymous and nonsynonymous substitution rates was not found either in other mammalian genes that contain a conservative box (Sox, Msx) or in the MADS-box genes of plants. The rate of nonsynonymous substitutions exceeds significantly that of synonymous substitutions in the terminal Sry sequences of apes. We did not find good support for the hypothesis that the high evolutionary rate of Sry would be associated with a promiscuous mating system.      

Somatic variation plays a key role in the evolution of the Vf gene family residing in the Vf locus that confers resistance to apple scab disease

A cluster of four receptor-like genes has been previously identified in the Vf locus of the crabapple Malus floribunda clone 821 that confers resistance to five races of the fungal pathogen Venturia inaequalis, the casual agent of apple scab disease. Pairwise comparisons of the four Vf paralogs in both promoter and coding regions reveal their timeline evolutionary history. The four Vf paralogs have evolved from four ancient Vf members resulting from two sequential duplication events of a single Vf progenitor initially present in the Malus genome. The coding sequences of the four Vf paralogs are characterized with high numbers of unique polymorphic nucleotides, a number of short duplications/deletions, various deletions of complete LRR copy units, and a casual insert of a transposon-like element. Significant high ratios of nonsynonymous to synonymous substitutions, Ka/Ks, are observed in the putative ligand binding residues in the LRR domains. No sequence exchange between the four Vf paralogs is observed. Compared with promoter regions, only nucleotide substitutions are dramatically elevated in the coding regions. The results presented in this study strongly indicate that the Vf locus is under strong and steady horizontal selective pressures imposed by the fungal pathogen V. inaequalis, and divergent selection on somatic variations plays a key role in shaping the resistance specificity.     

Pattern of synonymous and nonsynonymous substitutions: An indicator of mechanisms of molecular evolution

Comparison of numbers of synonymous and nonsynonymous substitutions is useful for understanding mechanisms of molecular evolution. In this paper, I examine the statistical properties of six methods of estimating numbers of synonymous and nonsynonymous substitutions. The six methods are Miyata and Yasunaga’s (MY) method; Nei and Gojobori’s (NG) method; Li, Wu and Luo’s (LWL) method; Pamilo, Bianchi and Li’s (PBL) method; and Ina’s (Ina) two methods. When the transition/transversion bias at the mutation level is strong, the numbers of synonymous and nonsynonymous substitutions are estimated more accurately by the PBL and Ina methods than by the NG, MY and LWL methods. When the nucleotide-frequency bias is strong and distantly related sequences are compared, all the six methods give underestimates of the number of synonymous substitutions. The concept of synonymous and nonsynonymous categories is also useful for analysis of DNA polymorphism data.     

Gene conversion as a source of nucleotide diversity in Plasmodium falciparum

Examination of polymorphisms in the Plasmodium falciparum gene for falcipain 2 revealed that this gene is one of two paralogs separated by 10.8 kb in chromosome 11. We designate the annotated gene denoted chr11.gen_424 as encoding falcipain 2A and the annotated gene denoted chr11.gen_427 as encoding falcipain 2B. The paralogs are 96% identical at the nucleotide level and 93% identical at the amino acid level. The consensus sequences differ in 31/309 synonymous sites and 45/1140 nonsynonymous sites, including three amino acid replacements (V393I, A400P, and Q414E) that are near the catalytic site and that may affect substrate affinity or specificity. In six reference isolates, among 36 synonymous sites and 46 nonsynonymous sites that are polymorphic in the gene for falcipain 2A, falcipain 2B, or both, significant spatial clustering is observed. All but one of the polymorphisms appear to result from gene conversion between the paralogs. The estimated rate of gene conversion between the paralogs may be as many as 1,400 to 1,700 times greater than the rate of mutation. Owing to gene conversion, one of the falcipain 2A alleles is more similar to the falcipain 2B alleles than it is to other falcipain 2A alleles. Divergence among the synonymous sites suggests that the paralogous genes last shared a common ancestor 15.2 MYA, with a range of 8.8 to 20.6 MYA. During this period, the paralogs have acquired 0.10 synonymous substitutions per synonymous site in the coding region. The 5' and 3' flanking regions differ in 47.7% and 39.8% of the nucleotide sites, respectively. Hence synonymous sites and flanking regions are not conserved in sequence in spite of their high AT content and T skew.     

The Na,K-ATPase alpha4 gene (Atp1a4) encodes a ouabain-resistant alpha subunit and is tightly linked to the alpha2 gene (Atp1a2) on mouse chromosome 1.

We have isolated and characterized cDNA clones encoding the murine homologue of a putative fourth Na,K-ATPase alpha subunit isoform (alpha4). The predicted polypeptide is 1032 amino acids in length and exhibits 75% amino acid sequence identity to the rat alpha1, alpha2, and alpha3 subunits. Within the first extracellular loop, the alpha4 subunit is highly divergent from other Na,K-ATPase alpha subunits. Because this region of Na,K-ATPase is a major determinant of ouabain sensitivity, we tested the ability of the rodent alpha4 subunit to transfer ouabain resistance in a transfection protocol. We find that a cDNA containing the complete rodent alpha4 ORF is capable of conferring low levels of ouabain resistance upon HEK 293 cells, an indication that the alpha4 subunit can substitute for the endogenous ouabain-sensitive alpha subunit of human cells. Nucleotide sequences specific for the murine alpha4 subunit were used to identify the chromosomal position of the alpha4 subunit gene. By hybridizing an alpha4 probe with a series of BACs, we localized the alpha4 subunit gene (Atp1a4) to the distal portion of mouse chromosome 1, in very close proximity to the murine Na,K-ATPase alpha2 subunit gene. In adult mouse tissues, we detected expression of the alpha4 subunit gene almost exclusively in testis, with low levels of expression in epididymis. The close similarities in the organization and expression pattern of the murine and human alpha4 subunit genes suggest that these two genes are orthologous. Together, our studies indicate that the alpha4 subunit represents a functional Na,K-ATPase alpha subunit isoform.     

中国地方鸡种MHC B-LBⅡ新等位基因的遗传多态性研究

为研究鸡MHC B-LBⅡ基因的遗传多态性,首先在8个中国地方鸡种(藏鸡、仙居鸡、北京油鸡、固始鸡、斗鸡、丝羽乌骨鸡、白耳鸡和狼山鸡)B-LBⅡ基因第二外显子扩增了一长度为 175 bp 的 DNA 片段并进行 SSCP 基因型分析;在8 个地方鸡种共 467 个个体中检测到 37 个 PCR-SSCP 基因型;从被检样品中筛选出不同基因型的个体,并在其 B-LBⅡ基因组中扩增了一个包括其第二外显子和第二内含子在内长度为374 bp的片段,通过克隆和测序获得了该片段的核苷酸序列。经序列分析,在前述地方鸡种被筛选出的 30 个无血缘关系的个体中发现了 31 个 B-LBⅡ新等位基因,并参照哺乳动物 MHC II 类 B 等位基因命名规则进行了命名。对这 31 个 B-LBⅡ新等位基因长度为 374 bp 的 DNA 片段进行比对表明,在其第二外显子序列上共有 68 个多态性变异位点,其中简约性信息位点 51 个,单变异位点 17 个,具有丰富的遗传多态性。在这些多态性变异位点中,出现在遗传密码子第一和第二位上的碱基替换率分别为 36.76% 和 35.29%。等位基因序列间的相似性估测为 90.6%-99.5%;B-LBⅡ基因第二外显子的错义替换率和同义替换率分别为 14.64±2.67%和 2.92±0.94%。结果表明,B-LBⅡ基因的丰富遗传多态性主要是由基因重组和平衡选择效应所引起的。对 B-LBⅡ等位基因第二外显子所编码的 B-LBⅡ分子β1 结构域氨基酸序列比对发现,31 个 B-LBⅡ新等位基因属于 26 个等位基因主型;在β1结构域氨基酸序列的 33个变异位点上,存在 6 个同义替换和 27 个错义替换。分析认为,那些发生在多肽结合位点上的氨基酸错义替换与鸡 MHC B-LBⅡ分子的免疫特异性有关。该结果可为鸡的抗病育种研究提供分子生物学依据。     

Molecular Population Genetics of Sex Determination Genes: The Transformer Gene of Drosophila Melanogaster

The transformer locus (tra) produces an RNA processing protein that alternatively splices the doublesex pre-mRNA in the sex determination hierarchy of Drosophila melanogaster. Comparisons of the tra coding region among Drosophila species have revealed an unusually high degree of divergence in synonymous and nonsynonymous sites. In this study, we tested the hypothesis that the tra gene will be polymorphic in synonymous and nonsynonymous sites within species by investigating nucleotide sequence variation in eleven tra alleles within D. melanogaster. Of the 1063 nucleotides examined, two synonymous sites were polymorphic and no amino acid variation was detected. Three statistical tests were used to detect departures from an equilibrium neutral model. Two tests failed to reject a neutral model of molecular evolution because of low statisitical power associated with low levels of genetic variation (Tajima/Fu and Li). The Hudson, Kreitman, and Aguade test rejected a neutral model when the tra region was compared to the 5'-flanking region of alcohol dehydrogenase (Adh). The lack of variability in the tra gene is consistent with a recent selective sweep of a beneficial allele in or near the tra locus.     

Transpecific polymorphisms in an inversion linked esterase locus in Drosophila buzzatii

Nucleotide variation was studied in a 1.1 kb section of the coding region of an Esterase gene (Est-A) that maps in the center of the segments rearranged by polymorphic inversions in the cactophilic Drosophila buzzatii. We examine 30 homozygous second-chromosome lines differing in gene arrangement and three D. koepferae isofemale lines as outgroups. Our data show that Est-A is a highly polymorphic gene at both synonymous and replacement sites. Significant departures from homogeneity in the distribution of the ratio of silent polymorphism to divergence predicted by the neutral theory reveals a local excess of silent polymorphism. This is consistent with the presence of two apparent narrow peaks of elevated silent polymorphism surrounding nonconservative amino acid substitutions. These polymorphisms as well as others at synonymous and nonsynonymous sites are shared with D. koepferae. We suggest that the presence of shared nucleotide polymorphisms is probably due to interspecific gene flow and/or balancing selection acting on replacement variants and/or to a decreased probability of loss of ancestral polymorphisms caused by linkage to an adaptive inversion polymorphism. Recurrent mutation and persistence of neutral ancestral polymorphisms cannot, however, be ruled out. The analysis of the distribution of nucleotide variation among the three chromosomal arrangements sampled reveals that derived arrangements (J and JZ(3)) are less polymorphic than the ancestral ST, and that the widely distributed ST and J arrangements are genetically differentiated. However, a significant number of polymorphisms are shared between arrangements, suggesting frequent exchange either from gene conversion or from double crossovers in heterokaryotypes. Finally, our present results in combination with data of sequence variation at the breakpoints of inversion J suggest that this old gene arrangement has risen in frequency in relatively recent times.     

A statistical examination of hypervariability in complementarity-determining regions of immunoglobulins

To determine the relative importance of gene conversion followed by natural selection and of natural selection for point mutation in generating variability in immunoglobulins, the numbers of synonymous and nonsynonymous substitutions in immunoglobulin sequences of various subgroups were estimated for complementarity-determining regions (CDRs) and for framework regions (FRs). Both the number of synonymous substitutions and the number of nonsynonymous substitutions in the CDR were found to exceed the corresponding numbers in the FR. Therefore, gene conversion is likely to be an important mechanism for providing variability in the CDR of immunoglobulins. The correlation coefficients between the number of synonymous substitutions and the number of nonsynonymous substitutions and between the substitution number in the CDR and that in the FR were found to be very low. Again, gene conversion is thought to be responsible for this finding.