Detection of homozygous genotypes for a putatively lethal recessive mutation in the porcine argininosuccinate synthase 1 (ASS1) gene

The sequencing of the pig genome revealed the existence of homozygous individuals for a nonsense mutation in the argininosuccinate synthase 1 (ASS1) gene (rs81212146, c.944T>A, L315X). Paradoxically, an AA homozygous genotype for this polymorphism is expected to abolish the function of the ASS1 enzyme that participates in the urea cycle, leading to citrullinemia, hyperammonemia, coma and death. Sequencing of five Duroc boars that sired a population of 350 Duroc barrows revealed the segregation of the c.944T>A polymorphism, so we aimed to investigate its phenotypic consequences. Genotyping of this mutation in the 350 Duroc barrows revealed the existence of seven individuals homozygous (AA) for the nonsense mutation. These AA pigs had a normal weight despite the fact that mild citrullinemia often involves impaired growth. Sequencing of the region surrounding the mutation in TT, TA and AA individuals revealed that the A substitution in the second position of the codon (c.944T>A) is in complete linkage disequilibrium with a C replacement (c.943T>C) in the first position of the codon. This second mutation would compensate for the potentially damaging effect of the c.944T>A replacement. In fact, this is the most probable reason why pigs with homozygous AA genotypes at the 944 site of the ASS1 coding region are alive. Our results illustrate the complexities of predicting the consequences of nonsense mutations on gene function and phenotypes, not only because of annotation issues but also owing to the existence of genetic mechanisms that sometimes limit the penetrance of highly harmful mutations.     

Evolution of exon 1 of the dopamine D4 receptor (DRD4) gene in primates

The dopamine D4 receptor (DRD4) gene exhibits a large amount of expressed polymorphism in humans. To understand the evolutionary history of the first exon of DRD4-which in humans contains a polymorphic 12bp tandem duplication, a polymorphic 13bp deletion, and other rare variants-we examined the homologous exon in thirteen other primate species. The great apes possess a variable number of tandem repeats in the same region as humans, both within and among species. In this sense, the 12bp tandem repeat of exon 1 is similar to the 48bp VNTR of exon 3 of DRD4, previously shown to be polymorphic in all primate species examined. The Old World monkeys show no variation in length, and a much higher conservation of amino acid sequence than great apes and humans. The New World monkeys show interspecific differences in length in the region of the 12bp polymorphism, but otherwise show the higher conservation seen in Old World monkeys. The different patterns of variation in monkeys compared to apes suggest strong purifying selective pressure on the exon in these monkeys, and somewhat different selection, possibly relaxed selection, in the apes.     

Evolution of the glycophorin gene family in the hominoid primates

Analysis of nucleotide sequences of the human glycophorin A (GPA) and glycophorin B (GPB) genes has indicated that the GPA gene most closely resembles the ancestral gene, whereas the GPB gene likely arose from the GPA gene by homologous recombination. To study the evolution of the glycophorin gene family in the hominoid primates, restricted DNA on Southern blots from man, pygmy chimpanzee, common chimpanzee, gorilla, orangutan, and gibbon was probed with cDNA fragments encoding the human GPA and GPB coding and 3-untranslated regions. This showed the presence in all of the hominoid primates of at least one GPA-like gene. In addition, at least one GPB-like gene was detected in man, both chimpanzee species, and gorilla, strongly suggesting that the event that produced the GPB gene occurred in the common ancestor of man-chimpanzee-gorilla. An unexpected finding in this study was the conservation ofEcoRI restriction sites relative to those of the other four enzymes used; the significance of this observation is unclear, but raises the question of nonrandomness ofEcoRI restriction sites in noncoding regions. Further analysis of the evolution of this multigene family, including nucleotide sequence analysis, will be useful in clarification of the evolutionary relationships of the hominoid primates, in correlation with the structure and function of the glycophorin molecules, and in assessment of the role of evolution in the autogenicity of glycophorin determinants.This work was supported in part by National Institutes of Health Grants AM33463 and CA33000.     

Molecular cloning and nucleotide sequence of the gramicidin S synthetase 1 gene

The entire gene for gramicidin S synthetase 1 (GS 1) was cloned into the plasmid vector pUC18, and the nucleotide sequences of the GS 1 gene and its flanking region were determined. The full-length clone was 4,539 base pairs long and had an open reading frame of 3,294 nucleotides coding for 1,098 amino acids. The calculated molecular weight of 123,474 agreed with the apparent molecular weight of 120,000 found in SDS-PAGE of GS 1 from B. brevis. The nucleotide sequence of GS 1 gene was highly homologous to that of tyrocidine synthetase 1. The overall similarity between the deduced amino acid sequences of the two genes was 57.5%. The gene product of clone GS309 was easily purified to an essentially homogeneous state by ammonium sulfate fractionation followed by DEAE-Sepharose CL-6B, Ultrogel AcA-34, and second DEAE-Sepharose CL-6B column chromatography. The purified protein catalyzed the D-phenylalanine-dependent ATP-32PPi exchange reaction which is specific for GS 1 activity, and the specific activity of the purified product was nearly the same as the purified GS 1 from B. brevis. The product also showed a weak phenylalanine racemase activity.     

Evolution of the AMP-forming acetyl-CoA synthetase gene in the Drosophilidae family

Analysis of the AMP-forming ACS gene was performed in 12 species of the Drosophilidae family. Systematically four introns, aligned at the same positions, were detected, but none of them showed a position similar to those known for species outside the Drosophilidae family. The average length of introns varied from 63 to 75 bp but in two species Drosophila takahashii and D. kikkawai the length of the second intron was 343 and 210 bp, respectively. In coding regions, about 80% of the third codon positions were substituted while first and second positions showed, respectively, 14% and 6% substitutions. Interestingly, the divergence observed at the protein level between species was very low. The phylogenetic tree based on the DNA sequences of the exons was mainly in agreement with taxonomic classification and previous molecular phylogenies except for D. ananassae, which appeared more closely related to D. subobscura and D. funebris than to the species of the melanogaster group.     

Evolution of a pigmentation gene,the melanocortin-1 receptor,in primates

The melanocortin-1 receptor (MC1R) forms a critical switch in the production of orange/red pheomelanin and black/brown eumelanin pigments during hair development in mammals. The molecular evolution of the melanocortin-1 receptor gene was investigated in a broad range of primate species, including several groups with large differences in distribution of orange/red and black hairs. Primate MC1R has been subject to purifying selection throughout most of its evolution, with small changes in selective constraint being detected early in primate evolution. In contrast to the situation in humans and domestic mammals, many intraspecific and intrageneric differences in primate coat color cannot be attributed to changes in the MC1R coding sequence. Nevertheless, important changes in the biochemical function of MC1R are suggested by mutations in sites of known functional importance, particularly in New World monkeys and lemurs. The evolution of the MC1R in lion tamarins is anomalous, with a combination of a high nonsynonymous to synonymous substitution rate (dN/dS) ratio, deletions, and substitutions.     

Evolution of a repeat sequence in the parathyroid hormone-related peptide gene in primates

A polymorphism of the variable number of tandem repeat (VNTR) type is located 97 bp downstream of exon VI of the parathyroid hormone-related peptide (PTHrP) gene in humans. The repeat unit has the general sequence G(TA)_nC, where n equals 4–11. In order to characterize the evolutionary history of this VNTR, we initially tested for its presence in 13 different species representing four main groups of living primates. The sequence is present in the human, great apes, and Old World monkeys, but not in New World monkeys; and this region failed to PCR amplify in the Loris group. Thus, the evolution of the sequence as part of the PTHrP gene started at least 25–35 millions years ago, after divergence of the Old World and New World monkeys, but before divergence of Old World monkeys and great apes and humans. The structural changes occurring during evolution are characterized by a relatively high degree of sequence divergence. In general, the tandem repeat region tends to be longer and more complex in higher primates with the repeat unit motifs all being based on a TA-dinucleotide repeat sequence. Intra-species variability of the locus was demonstrated only in humans and gorilla. The divergence of the TA-dinucleotide repeat sequence and the variable mutation rates observed in different primate species are in contrast to the relative conservation of the flanking sequences during primate evolution. This suggests that the nature of the TA-dinucleotide repeat sequence, rather than its flanking sequences, is responsible for generating variability. Particular features of the sequence may allow it to form stable secondary structures during DNA replication, and this, in turn, could promote slipped-strand mispairing to occur.     

Chromosomal evolution of the PKD1 gene family in primates

Correction to Kirsch S, Pasantes J, Wolf A, Bogdanova N, Münch C, Pennekamp P, Krawczak M, Dworniczak B, Schempp W: Chromosomal evolution of the PKD1 gene family in primates. BMC Evolutionary Biology 2008, 8 :263 (doi:10.1186/1471-2148-8-263)     

High intra-individual sequence variation in the nuclear rDNA LSU-5S intergenic spacer in the Sargassaceae (Fucales, Phaeophyceae)

Of four species of Sargassaceae, representing the genera Carpophyllum, Cystoseira, Landsburgia, and Sargassum, the intergenic spacer of the ribosomal cistron was amplified using a forward primer annealing at the 3′-end of the large subunit (LSU) of the ribosomal cistron and a reverse primer annealing at the 5S rRNA gene. The PCR products were cloned and the DNA sequences of multiple clones were determined. Almost each clone showed a unique DNA sequence. Intra-individual variation of this LSU-5S intergenic spacer was extremely high and was characterized by great length variation and a high number of short tandem repeats. Sequences were unalignable and therefore it was concluded that the LSU-5S intergenic spacer is unsuitable for phylogenetic and phylogeographic studies of sargassacean taxa.     

Molecular phylogenetics of the family Cyprinidae (Actinopterygii: Cypriniformes) as evidenced by sequence variation in the first intron of S7 ribosomal protein-coding gene: further evidence from a nuclear gene of the systematic chaos in the family

The family Cyprinidae is the largest freshwater fish group in the world, including over 200 genera and 2100 species. The phylogenetic relationships of major clades within this family are simply poorly understood, largely because of the overwhelming diversity of the group; however, several investigators have advanced different hypotheses of relationships that pre- and post-date the use of shared-derived characters as advocated through phylogenetic systematics. As expected, most previous investigations used morphological characters. Recently, mitochondrial DNA (mtDNA) sequences and combined morphological and mtDNA investigations have been used to explore and advance our understanding of species relationships and test monophyletic groupings. Limitations of these studies include limited taxon sampling and a strict reliance upon maternally inherited mtDNA variation. The present study is the first endeavor to recover the phylogenetic relationships of the 12 previously recognized monophyletic subfamilies within the Cyprinidae using newly sequenced nuclear DNA (nDNA) for over 50 species representing members of the different previously hypothesized subfamily and family groupings within the Cyprinidae and from other cypriniform families as outgroup taxa. Hypothesized phylogenetic relationships are constructed using maximum parsimony and Basyesian analyses of 1042 sites, of which 971 sites were variable and 790 were phylogenetically informative. Using other appropriate cypriniform taxa of the families Catostomidae (Myxocyprinus asiaticus), Gyrinocheilidae (Gyrinocheilus aymonieri), and Balitoridae (Nemacheilus sp. and Beaufortia kweichowensis) as outgroups, the Cyprinidae is resolved as a monophyletic group. Within the family the genera Raiamas, Barilius, Danio, and Rasbora, representing many of the tropical cyprinids, represent basal members of the family. All other species can be classified into variably supported and resolved monophyletic lineages, depending upon analysis, that are consistent with or correspond to Barbini and Leuciscini. The Barbini includes taxa traditionally aligned with the subfamily Cyprininae sensu previous morphological revisionary studies by Howes (Barbinae, Labeoninae, Cyprininae and Schizothoracinae). The Leuciscini includes six other subfamilies that are mainly divided into three separate lineages. The relationships among genera and subfamilies are discussed as well as the possible origins of major lineages.     

Nucleotide sequence and polymorphism of the pig acyl coenzyme A synthetase long-chain 1 (ACSL1) gene

We have sequenced 3,013 bp of the pig acyl coenzyme A long-chain synthetase 1 (ACSL1) gene. Sequence analysis allowed us to identify three conserved elements in the predicted amino acid sequence, two of them related to the ATP/AMP signature motif and the third involved in enzyme catalysis and fatty acid substrate specificity. In addition, we have identified five C --> T and one G --> A transition polymorphisms located in exon 16 (SNPe16), exon 17 (SNPe17) and the 3' UTR (SNPa to d), which have been genotyped in 143 pigs from the Landrace, Large White, Piétrain, Iberian, and Duroc breeds.     

The human debrisoquine 4-hydroxylase (CYP2D) locus: sequence and identification of the polymorphic CYP2D6 gene, a related gene, and a pseudogene.

The debrisoquine-4-hydroxylase polymorphism is a genetic variation in oxidative drug metabolism characterized by two phenotypes, the extensive metabolizer (EM) and poor metabolizer (PM). Of the Caucasian populations of Europe and North America, 5%-10% are of the PM phenotype and are unable to metabolize debrisoquine and numerous other drugs. The defect is caused by several mutant alleles of the CYP2D6 gene, two of which are detected in about 70% of PMs. We have constructed a genomic library from lymphocyte DNA of an EM positively identified by pedigree analysis to be homozygous for the normal CYP2D6 allele. The normal CYP2D6 gene was isolated; was completely sequenced, including 1,531 and 3,522 bp of 5' and 3' flanking DNA, respectively; and was found to contain nine exons within 4,378 bp. Two other genes, designated CYP2D7 and CYP2D8P, were also cloned and sequenced. CYP2D8P contains several gene-disrupting insertions, deletions, and termination codons within its exons, indicating that this is a pseudogene. CYP2D7, which is just downstream of CYP2D8P, is apparently normal, except for the presence, in the first exon, of an insertion that disrupts the reading frame. A hypothesis is presented that the presence of a pseudogene within the CYP2D subfamily transfers detrimental mutations via gene conversions into the CYP2D6 gene, thus accounting for the high frequency of mutations observed in the CYP2D6 gene in humans.     

Phylogeography of the endangered Cathaya argyrophylla (Pinaceae) inferred from sequence variation of mitochondrial and nuclear DNA

Cathaya argyrophylla is an endangered conifer restricted to subtropical mountains of China. To study phylogeographical pattern and demographic history of C. argyrophylla, species-wide genetic variation was investigated using sequences of maternally inherited mtDNA and biparentally inherited nuclear DNA. Of 15 populations sampled from all four distinct regions, only three mitotypes were detected at two loci, without single region having a mixed composition (G(ST) = 1). Average nucleotide diversity (theta(ws) = 0.0024; pi(s) = 0.0029) across eight nuclear loci is significantly lower than those found for other conifers (theta(ws) = 0.003 approximately 0.015; pi(s) = 0.002 approximately 0.012) based on estimates of multiple loci. Because of its highest diversity among the eight nuclear loci and evolving neutrally, one locus (2009) was further used for phylogeographical studies and eight haplotypes resulting from 12 polymorphic sites were obtained from 98 individuals. All the four distinct regions had at least four haplotypes, with the Dalou region (DL) having the highest diversity and the Bamian region (BM) the lowest, paralleling the result of the eight nuclear loci. An AMOVA revealed significant proportion of diversity attributable to differences among regions (13.4%) and among populations within regions (8.9%). F(ST) analysis also indicated significantly high differentiation among populations (F(ST) = 0.22) and between regions (F(ST) = 0.12-0.38). Non-overlapping distribution of mitotypes and high genetic differentiation among the distinct geographical groups suggest the existence of at least four separate glacial refugia. Based on network and mismatch distribution analyses, we do not find evidence of long distance dispersal and population expansion in C. argyrophylla. Ex situ conservation and artificial crossing are recommended for the management of this endangered species.     

Evolution of the NANOG pseudogene family in the human and chimpanzee genomes

Background  

The NANOG gene is expressed in mammalian embryonic stem cells where it maintains cellular pluripotency. An unusually large family of pseudogenes arose from it with one unprocessed and ten processed pseudogenes in the human genome. This article compares the NANOG gene and its pseudogenes in the human and chimpanzee genomes and derives an evolutionary history of this pseudogene family.     

Characterization of a yeast nuclear gene (MST1) coding for the mitochondrial threonyl-tRNA1 synthetase

The wild-type yeast nuclear gene MST1 complements mutants defective in mitochondrial protein synthesis. The gene has been sequenced and shown to code for a protein of 54,030 kDa. The predicted product of MST1 is 36% identical over its 462 residues to the Escherichia coli threonyl-tRNA synthetase. Amino-acylation of wild-type mitochondrial tRNAs with a mitochondrial extract from mst1 mutants fail to acylate tRNAThr1 (anticodon: 3'-GAU-5') but show normal acylation of tRNAThr2 (anticodon: 3'-UGU-5'). These data suggest the presence of two separate threonyl-tRNA synthetases in yeast mitochondria. Antibodies were prepared against a trpE/MST1 fusion protein containing the 321 residues from the amino-terminal region of the E. coli anthranilate synthetase and 118 residues of the mitochondrial threonyl-tRNA synthetase. Antibodies to the fusion protein detect a 50-55-kDa protein in wild type yeast mitochondria but not in mitochondria of a strain in which the chromosomal MST1 gene was replaced by a copy of the same gene disrupted by insertion of the yeast LEU2 gene. The ability of the mutant with the inactive MST1 gene to charge tRNAThr2 argues strongly for the existence of a second threonyl-tRNA synthetase gene.