Nucleotide diversity and molecular evolution of the WAG-2 gene in common wheat (Triticum aestivum L) and its relatives

In this work, we examined the genetic diversity and evolution of the WAG-2 gene based on new WAG-2 alleles isolated from wheat and its relatives. Only single nucleotide polymorphisms (SNP) and no insertions and deletions (indels) were found in exon sequences of WAG-2 from different species. More SNPs and indels occurred in introns than in exons. For exons, exons+introns and introns, the nucleotide polymorphism π decreased from diploid and tetraploid genotypes to hexaploid genotypes. This finding indicated that the diversity of WAG-2 in diploids was greater than in hexaploids because of the strong selection pressure on the latter. All dn/ds ratios were < 1.0, indicating that WAG-2 belongs to a conserved gene affected by negative selection. Thirty-nine of the 57 particular SNPs and eight of the 10 indels were detected in diploid species. The degree of divergence in intron length among WAG-2 clones and phylogenetic tree topology suggested the existence of three homoeologs in the A, B or D genome of common wheat. Wheat AG-like genes were divided into WAG-1 and WAG-2 clades. The latter clade contained WAG-2, OsMADS3 and ZMM2 genes, indicating functional homoeology among them.     

Structure and evolution of the human prosaposin chromosomal gene.

The gene for prosaposin was characterized by sequence analysis of chromosomal DNA to gain insight into the evolution of this locus that encodes four highly conserved sphingolipid activator proteins or saposins. The 13 exons ranged in size from 57 to 1200 bp, while the introns were from 91 to 3812 bp in length. The regions encoding saposins A, B, and D each had three exons, while that for saposin C had only two. This sequence included the regions that encode the carboxy terminus of the signal peptide, the four mature prosaposin proteins, and the 3' untranslated region. Primer extension studies indicated that over 99% of the coding sequence was contained in these 19,985 bp. Use of PCR and reverse PCR techniques indicated that the most 5' coding approximately 140 bp contained large introns and at least two small exons. Analyses of the intronic positions in the saposin regions indicated that this gene evolved from an ancestral gene by two duplication events and at least one gene rearrangement involving a double crossover after introns had been inserted into the gene.     

Asymmetrical distribution of CpG in an 'average' mammalian gene.

The frequency and distribution of the rare dinucleotide CpG was examined in 15 mammalian genes. CpG is highly methylated at cytosine in mammalian DNA (1,2) and 5-methylcytosine (5mC) is thought to undergo a transition mutation via deamination to produce thymine (3). This would result in the accumulation of TpG and CpA and depletion of CpG during evolution (4). Consistent with this hypothesis, the gene sample of 26,541 dinucleotides contained CpG at 40% the frequency expected by base composition and the CpG transition products, TpG+CpA, were significantly elevated at 124% of expected random frequency. However, because CpG occurs at only 25% of expected random frequency in the genome, the sampled genes were considerably enriched in this dinucleotide. CpGs were asymmetrically distributed in sequences flanking the genes. 5'-flanking sequences were enriched in CpG at 135% of the frequency expected assuming a symmetrical distribution of all the CpGs in the sampled genes (p less than 0.01), while 3'-flanking regions were depleted in CpG at 40% of expected values (p less than 0.0001). This asymmetry may reflect the role of 5-methylcytosine in gene expression. In contrast the frequencies of GpC and GpT+ ApC did not differ significantly from that predicted by base composition and these dinucleotides were not asymmetrically distributed.     

Concerted and divergent evolution within the rat gamma-crystallin gene family

The nucleotide sequences of six rat gamma-crystallin genes have been determined. All genes have the same mosaic structure: the first exons contain a relatively short (25 to 44 base-pair) 5' non-coding region and the first nine base-pairs of the coding sequence, the second exons encode protein motifs I and II, while protein motifs III and IV are encoded by the third exons. The third exons also contain a 60 to 67-base-pair long 3' non-coding region. In the gamma 1-2 gene, the splice acceptor site of the third exon has been shifted three base-pairs upstream. Hence, the protein product of this gene is one amino acid residue longer. The first introns, though varying in length from 85 to 100 base-pairs, are conserved in sequence. The second introns vary considerably in length (0.9 X 10(3) to 1.9 X 10(3) base-pairs) and sequence. The second exons of the genes show concerted evolution and have undergone multiple gene conversions. In contrast, the third exons show divergent evolution. From the sequences of the third exons, an evolutionary tree of the gene family was constructed. This tree suggests that three of the present genes derive directly from the genes that originated from a tandem duplication of a two-gene cluster. Two duplications of the last gene of the four-gene cluster then yielded the other three genes. Region a' of the third exon, encoding protein motif III, is variable, while the region encoding protein motif IV (b') is constant. We postulate that this variability in region a' is due to a period of radiation after each gene duplication. A comparison of the rat sequences with those of orthologous sequences from other species shows that the variation in region a' is now preserved. Hence, it might specify the specific functional property of each gamma-crystallin protein within the lens.     

Utility of NADP-dependent isocitrate dehydrogenase for species-level evolutionary inference in angiosperm phylogeny: a case study in Saltugilia

NADP-dependent isocitrate dehydrogenase is a low-copy nuclear gene family. We have sequenced two regions from an idh gene (idhB) near the 3' terminal end. The first fragment encodes 4 exons and 3 introns and is between approximately 600 and 950 bp in length. The second fragment includes three additional exons and introns and is between approximately 1200 and 1500 bp in length. The phylogenetic utility of the two sequence regions was evaluated in Polemoniaceae with a focus on Saltugilia, an incipient species complex that lacks phylogenetic resolution among these same taxa based on nuclear ribosomal ITS and chloroplast trnL. Multiple sequences from several individuals, multiple individuals from several populations, and multiple populations from all Saltugilia species were sampled to evaluate the taxonomic level at which idhB was useful as a phylogenetic marker in this clade. Phylogenies based on idhB sequences were compared with topological resolution and clade composition in ITS and trnL phylogenies. Phylogenies based on idhB and idhB in combination with ITS and trnL are better resolved than any other phylogenies for Saltugilia published to date, and character evolution within Saltugilia is explored.     

Structure and evolution of the bovine prothrombin gene

The cloned bovine prothrombin gene has been characterized by partial DNA sequence analysis, including the 5' and 3' flanking sequences and all the intron-exon junctions. The gene is approximately 15.4 x 10(3) base-pairs in length and comprises 14 exons interrupted by 13 introns. The exons coding for the prepro-leader peptide and the gamma-carboxyglutamic acid-containing region are similar in organization to the corresponding exons in the factor IX and protein C genes. This region has probably evolved as a result of recent gene duplication and exon shuffling events. The exons coding for the kringles and the serine protease region of the prothrombin gene are different in organization from the homologous regions in other genes, suggesting that introns have been inserted into these regions after the initial gene duplication events.     

Sequence of EndoA gene encoding mouse cytokeratin and its methylation state in the CpG-rich region.

A genomic clone obtained from mouse liver DNA using a mouse cytokeratin EndoA cDNA probe revealed the complete sequence of the EndoA gene. The gene is divided into nine exons and the exon-intron pattern has been conserved compared to that of other type-II cytokeratin-encoding genes. The 5' upstream, 3' downstream and first and third introns contain potential regulatory sequences, including polyoma virus enhancer motifs (PEA1 and PEA3) and AP-1 elements. The 5' regions upstream of the EndoA, EndoB and Ck8 genes contain homologous sequences surrounding the TATA boxes. In addition, a CpG dinucleotide cluster region was located around the first exon. This CpG cluster region was found to be hypomethylated in endodermal PYS-2 cells, retinoic acid-treated F9 cells, and F9 embryonal carcinoma cells, but hypermethylated in BALB/C 3T3 fibroblast cells that do not express EndoA. These findings may provide a clue to understanding the molecular mechanisms of EndoA gene expression.     

Diversity and diversification of HLA-A,B,C alleles

The nucleotide sequences encoding 14 HLA-A,B,C and 5 ChLA-A,B,C molecules have been determined. Combining these sequences with published data has enabled the polymorphism in 40 HLA-A,B,C and 9 ChLA-A,B,C alleles to be analyzed. Diversity is generated through assortment of point mutations by recombinational mechanisms including gene and allelic conversions. The distribution and frequency of silent and replacement substitutions indicate that there has been positive selection for allelic diversity in the 5' part of the gene (exons 1 to 3) and for allelic homogenization and locus specificity in the 3' part of the gene (exons 4 to 8). These differences may correlate with the lengths of converted sequences in the two parts of the gene and frequency of the CpG dinucleotide. Locus-specific divergence of HLA-A,B, and C demonstrates that recombinational events involving alleles of a locus have been more important than conversion between loci. This contrasts with the predominance of gene conversion events in the evolution of mutants of the H-2Kb gene. However, a striking example of gene conversion involving HLA-B and C alleles of an oriental haplotype has been found. Comparison of human and chimpanzee alleles reveals extensive sharing of polymorphisms, confirming that diversification is a slow process, and that much of contemporary polymorphism originated in ancestral primate species before the emergence of Homo sapiens. There is less polymorphism at the HLA-A locus compared to HLA-B, with greater similarity also being seen between HLA-A and ChLA-A alleles than between HLA-B and ChLA-B alleles. Although greater diversity is seen in the 5' "variable" exons of HLA-B compared to HLA-A, there is increased heterogeneity in the 3' "conserved" exons of HLA-A compared to HLA-B.     

Natural selection affects frequencies of AG and GT dinucleotides at the 5' and 3' ends of exons

GT and AG, located at the 5' and 3' ends of introns, are important for correct splicing. It is anticipated that natural selection decreases frequency of AG and GT near the 5' and 3' ends of exons, preventing appearance of cryptic splicing sites. The data presented in this article support the expectation.     

Alternative splicing of RNA triplets is often regulated and accelerates proteome evolution

Thousands of human genes contain introns ending in NAGNAG (N any nucleotide), where both NAGs can function as 3' splice sites, yielding isoforms that differ by inclusion/exclusion of three bases. However, few models exist for how such splicing might be regulated, and some studies have concluded that NAGNAG splicing is purely stochastic and nonfunctional. Here, we used deep RNA-Seq data from 16 human and eight mouse tissues to analyze the regulation and evolution of NAGNAG splicing. Using both biological and technical replicates to estimate false discovery rates, we estimate that at least 25% of alternatively spliced NAGNAGs undergo tissue-specific regulation in mammals, and alternative splicing of strongly tissue-specific NAGNAGs was 10 times as likely to be conserved between species as was splicing of non-tissue-specific events, implying selective maintenance. Preferential use of the distal NAG was associated with distinct sequence features, including a more distal location of the branch point and presence of a pyrimidine immediately before the first NAG, and alteration of these features in a splicing reporter shifted splicing away from the distal site. Strikingly, alignments of orthologous exons revealed a ～15-fold increase in the frequency of three base pair gaps at 3' splice sites relative to nearby exon positions in both mammals and in Drosophila. Alternative splicing of NAGNAGs in human was associated with dramatically increased frequency of exon length changes at orthologous exon boundaries in rodents, and a model involving point mutations that create, destroy, or alter NAGNAGs can explain both the increased frequency and biased codon composition of gained/lost sequence observed at the beginnings of exons. This study shows that NAGNAG alternative splicing generates widespread differences between the proteomes of mammalian tissues, and suggests that the evolutionary trajectories of mammalian proteins are strongly biased by the locations and phases of the introns that interrupt coding sequences.     

Isolation and nucleotide sequence analysis of the beta-type globin pseudogene from human, gorilla and chimpanzee

The beta-globin gene cluster of human, gorilla and chimpanzee contain the same number and organization of beta-type globin genes: 5'-epsilon (embryonic)-G gamma and A gamma (fetal)-psi beta (inactive)-delta and beta (adult)-3'. We have isolated the psi beta-globin gene regions from the three species and determined their nucleotide sequences. These three pseudogenes each share the same substitutions in the initiator codon (ATG----GTA), a substitution in codon 15 which generates a termination signal TGG----TGA, nucleotide deletion in codon 20 and the resulting frame shift which yields many termination signals in exons 2 and 3. The basic structure of these psi beta-globin genes, however, remains consistent with that found for functional beta-globin genes: their coding regions are split by two introns, IVS 1 (which splits codon 30, 121 base-pairs in length) and IVS 2 (which splits codon 104, 840 to 844 base-pairs in length). These introns retain the normal splice junctions found in other eukaryotic split genes. The three hominoid psi beta-globin genes show a high degree of sequence correspondence, with the number of differences found among them being only about one-third of that predicted for DNA sites evolving at the neutral rate (i.e. for sites evolving in the absence of purifying selection). Thus, there appears to be a deceleration in the rate of evolution of the psi beta-globin locus in higher primates.