The nuclear autoantigen La/SS-B: mapping and sequencing of the gene and the three retropseudogenes

One target of autoantibodies in sera of patients with systemic lupus erythematosus or primary Sjögren's syndrome is the nuclear autoantigen La/SS-B. Lambda clones and cosmids were isolated, which contained the sequences of the La gene and the three La pseudogenes. They were used for preparation of a physical map. Finally, the La gene and pseudogenes were sequenced. The pseudogenes were characterized as retropseudogenes. Their evolutionary ages were estimated to be approx. 4, 4.5 and 5 million years. Inserts of 4, 16 and 24 nucleotides, which were mostly A-residues, were found in exon 7 of the respective pseudogene. The oldest pseudogene contained the longest insert, the youngest pseudogene contained the smallest insert. The oligonucleotides seem to be the result of repeated inserts of A-residues in a hot spot region of the La genes. Two La cDNAs were isolated which contained either a deletion or an insert of an A-residue at the same position.     

Comparison of complete nuclear receptor sets from the human, Caenorhabditis elegans and Drosophila genomes

Background

The availability of complete genome sequences enables all the members of a gene family to be identified without limitations imposed by temporal, spatial or quantitative aspects of mRNA expression. Using the nearly completed human genome sequence, we combined in silico and experimental approaches to define the complete human nuclear receptor (NR) set. This information was used to carry out a comparative genomic study of the NR superfamily.

Results

Our analysis of the human genome identified two novel NR sequences. Both these contained stop codons within the coding regions, indicating that both are pseudogenes. One (HNF4 γ-related) contained no introns and expressed no detectable mRNA, whereas the other (FXR-related) produced mRNA at relatively high levels in testis. If translated, the latter is predicted to encode a short, non-functional protein. Our analysis indicates that there are fewer than 50 functional human NRs, dramatically fewer than in Caenorhabditis elegans and about twice as many as in Drosophila. Using the complete human NR set we made comparisons with the NR sets of C. elegans and Drosophila. Searches for the >200 NRs unique to C. elegans revealed no human homologs. The comparative analysis also revealed a Drosophila member of NR subfamily NR3, confirming an ancient metazoan origin for this subfamily.

Conclusions

This work provides the basis for new insights into the evolution and functional relationships of NR superfamily members.     

Mouse cytoskeletal gamma-actin: analysis and implications of the structure of cloned cDNA and processed pseudogenes

The nucleotide sequence corresponding to almost the whole of a mouse gamma-cytoskeletal actin mRNA was determined from overlapping cloned DNA copies derived from brain mRNA. Several gamma-actin processed pseudogenes were isolated from a library of cloned DBA mouse genomic DNA, and the nucleotide sequences of these were determined and compared with that of the cDNA. This showed that two of these pseudogenes had arisen from a gene duplication or amplification event, and indicated that they had subsequently undergone partial correction against one another. The relative ages of the pseudogenes were estimated on the basis of their percentage divergence from the cDNA sequence and these were compared with an estimation based on the number of presumed silent mutations in the cDNA since each pseudogene had arisen. Consistent results were obtained, except in the case of one pseudogene which also showed an anomalous regional distribution of differences from the cDNA sequence. One way of accounting for the features of this anomalous pseudogene is by postulating that it is derived from a second functional gene for gamma-actin, different from that represented by the cDNA described here.     

Eleven daughters of NANOG

Nanog is a recently discovered ANTP class homeobox gene. Mouse Nanog is expressed in the inner cell mass and in embryonic stem cells and has roles in self-renewal and maintenance of pluripotency. Here we describe the location, genomic organization, and relative ages of all human NANOG pseudogenes, comprising ten processed pseudogenes and one tandem duplicate. These are compared to the original, intact human NANOG gene. Eleven is an unusually high number of pseudogenes for a homeobox gene and must reflect expression in the human germ line. A pseudogene orthologous to NANOGP4 was found in chimpanzee and an expressed pseudogene in macaque. Examining pseudogenes of differing ages gives insight into pseudogene decay, which involves an excess of deletion mutations over insertions. The mouse genome has two processed pseudogenes, which are not clear orthologues of the primate pseudogenes.     

Comparative genomics search for losses of long-established genes on the human lineage

Taking advantage of the complete genome sequences of several mammals, we developed a novel method to detect losses of well-established genes in the human genome through syntenic mapping of gene structures between the human, mouse, and dog genomes. Unlike most previous genomic methods for pseudogene identification, this analysis is able to differentiate losses of well-established genes from pseudogenes formed shortly after segmental duplication or generated via retrotransposition. Therefore, it enables us to find genes that were inactivated long after their birth, which were likely to have evolved nonredundant biological functions before being inactivated. The method was used to look for gene losses along the human lineage during the approximately 75 million years (My) since the common ancestor of primates and rodents (the euarchontoglire crown group). We identified 26 losses of well-established genes in the human genome that were all lost at least 50 My after their birth. Many of them were previously characterized pseudogenes in the human genome, such as GULO and UOX. Our methodology is highly effective at identifying losses of single-copy genes of ancient origin, allowing us to find a few well-known pseudogenes in the human genome missed by previous high-throughput genome-wide studies. In addition to confirming previously known gene losses, we identified 16 previously uncharacterized human pseudogenes that are definitive losses of long-established genes. Among them is ACYL3, an ancient enzyme present in archaea, bacteria, and eukaryotes, but lost approximately 6 to 8 Mya in the ancestor of humans and chimps. Although losses of well-established genes do not equate to adaptive gene losses, they are a useful proxy to use when searching for such genetic changes. This is especially true for adaptive losses that occurred more than 250,000 years ago, since any genetic evidence of the selective sweep indicative of such an event has been erased.     

More genes in vertebrates?

With the acquisition of complete genome sequences from several animals, there is renewed interest in the pattern of genome evolution on our own lineage. One key question is whether gene number increased during chordate or vertebrate evolution. It is argued here that comparing the total number of genes between a fly, a nematode and human is not appropriate to address this question. Extensive gene loss after duplication is one complication; another is the problem of comparing taxa that are phylogenetically very distant. Amphioxus and tunicates are more appropriate animals for comparison to vertebrates. Comparisons of clustered homeobox genes, where gene loss can be identified, reveals a one to four mode of evolution for Hox and ParaHox genes. Analyses of other gene families in amphioxus and vertebrates confirm that gene duplication was very widespread on the vertebrate lineage. These data confirm that vertebrates have more genes than their closest invertebrate relatives, acquired through gene duplication. abbreviations IHGSC, International Human Genome Sequencing Consortium; TCESC, The C. elegans Sequencing Consortium.     

Segmental duplications in the human genome reveal details of pseudogene formation

Duplicated pseudogenes in the human genome are disabled copies of functioning parent genes. They result from block duplication events occurring throughout evolutionary history. Relatively recent duplications (with sequence similarity ≥90% and length ≥1 kb) are termed segmental duplications (SDs); here, we analyze the interrelationship of SDs and pseudogenes. We present a decision-tree approach to classify pseudogenes based on their (and their parents’) characteristics in relation to SDs. The classification identifies 140 novel pseudogenes and makes possible improved annotation for the 3172 pseudogenes located in SDs. In particular, it reveals that many pseudogenes in SDs likely did not arise directly from parent genes, but are the result of a multi-step process. In these cases, the initial duplication or retrotransposition of a parent gene gives rise to a ‘parent pseudogene’, followed by further duplication creating duplicated–duplicated or duplicated–processed pseudogenes, respectively. Moreover, we can precisely identify these parent pseudogenes by overlap with ancestral SD loci. Finally, a comparison of nucleotide substitutions per site in a pseudogene with its surrounding SD region allows us to estimate the time difference between duplication and disablement events, and this suggests that most duplicated pseudogenes in SDs were likely disabled around the time of the original duplication.     

Evolutionary selection across the nuclear hormone receptor superfamily with a focus on the NR1I subfamily (vitamin D,pregnane X,and constitutive androstane receptors)

Background

The nuclear hormone receptor (NR) superfamily complement in humans is composed of 48 genes with diverse roles in metabolic homeostasis, development, and detoxification. In general, NRs are strongly conserved between vertebrate species, and few examples of molecular adaptation (positive selection) within this superfamily have been demonstrated. Previous studies utilizing two-species comparisons reveal strong purifying (negative) selection of most NR genes, with two possible exceptions being the ligand-binding domains (LBDs) of the pregnane X receptor (PXR, NR1I2) and the constitutive androstane receptor (CAR, NR1I3), two proteins involved in the regulation of toxic compound metabolism and elimination. The aim of this study was to apply detailed phylogenetic analysis using maximum likelihood methods to the entire complement of genes in the vertebrate NR superfamily. Analyses were carried out both across all vertebrates and limited to mammals and also separately for the two major domains of NRs, the DNA-binding domain (DBD) and LBD, in addition to the full-length sequences. Additional functional data is also reported for activation of PXR and the vitamin D receptor (VDR; NR1I1) to gain further insight into the evolution of the NR1I subfamily.

Results

The NR genes appear to be subject to strong purifying selection, particularly in the DBDs. Estimates of the ratio of the non-synonymous to synonymous nucleotide substitution rates (the ω ratio) revealed that only the PXR LBD had a sub-population of codons with an estimated ω ratio greater than 1. CAR was also unusual in showing high relative ω ratios in both the DBD and LBD, a finding that may relate to the recent appearance of the CAR gene (presumably by duplication of a pre-mammalian PXR gene) just prior to the evolution of mammals. Functional analyses of the NR1I subfamily show that human and zebrafish PXRs show similar activation by steroid hormones and early bile salts, properties not shared by sea lamprey, mouse, or human VDRs, or by Xenopus laevis PXRs.

Conclusion

NR genes generally show strong sequence conservation and little evidence for positive selection. The main exceptions are PXR and CAR, genes that may have adapted to cross-species differences in toxic compound exposure.

     

Ancient duplications of the human proglucagon gene

The human proglucagon gene (GCG) is encoded within a finished 576-kb DNA sequence generated by the Human Genome Project. GCG is flanked by 18 kb and 65 kb of DNA, 5' and 3', respectively, that do not encode genes. The genomic sequence that includes GCG was found to have a long history of gene duplication events. Some members of the glucagon-like family of genes, GCG on chromosome 2 and GIP on chromosome 17, may be products of ancient genome duplications on the early vertebrate lineage. A large genomic tandem duplication event that included DPP4-like and GCG genes occurred before the amphibian-mammal divergence, but one of the duplicated copies of GCG has been lost on the human lineage. Recently, a processed pseudogene of the X-chromosome-linked gene TIMM8A was inserted downstream of GCG. Some ancient duplicates of GCG may retain physiological functions in other vertebrates.     

The farnesoid X receptor FXRalpha/NR1H4 acquired ligand specificity for bile salts late in vertebrate evolution

The nuclear receptor FXRalpha (NR1H4) plays a pivotal role in maintaining bile salt and lipid homeostasis by functioning as a bile salt sensor in mammals. In contrast, FXRbeta (NR1H5) from mouse is activated by lanosterol and does not share common ligands with FXRalpha. To further elucidate FXR ligand/receptor and structure/function relationships, we characterized a FXR gene from the marine skate, Leucoraja erinacea, representing a vertebrate lineage that diverged over 400 million years ago. Phylogenetic analysis of sequence data indicated that skate Fxr (sFxr) is a FXRbeta. There is an extra sequence in the middle of the sFxr ligand binding domain (LBD) compared with the LBD of FXRalpha. Luciferase reporter assays demonstrated that sFxr responds weakly to scymnol sulfate, bile salts, and synthetic FXRalpha ligands, in striking difference from human FXRalpha (hFXRalpha). Interestingly, all-trans retinoic acid was capable of transactivating both hFXRalpha and sFxr. When the extra amino acids in the sFxr LBD were deleted and replaced with the corresponding sequence from hFXRalpha, the mutant sFxr gained responsiveness to ursodeoxycholic acid, GW4064, and fexaramine. Surprisingly, chenodeoxycholic acid antagonized this activation. Together, these results indicate that FXR is an ancient nuclear receptor and suggest that FXRalpha may have acquired ligand specificity for bile acids later in evolution by deletion of a sequence from its LBD. Acquisition of this property may be an example of molecular exploitation, where an older molecule is recruited for a new functional role.     

A hotspot of gene order rearrangement by tandem duplication and random loss in the vertebrate mitochondrial genome

Most reported examples of change in vertebrate mitochondrial (mt) gene order could be explained by a tandem duplication followed by random loss of redundant genes (tandem duplication-random loss [TDRL] model). Under this model of evolution, independent loss of genes arising from a single duplication in an ancestral species are predicted, and remnant pseudogenes expected, intermediate states that may remain in rearranged genomes. However, evidence for this is rare and largely scattered across vertebrate lineages. Here, we report new derived mt gene orders in the vertebrate "WANCY" region of four closely related caecilian amphibians. The novel arrangements found in this genomic region (one of them is convergent with the derived arrangement of marsupials), presence of pseudogenes, and positions of intergenic spacers fully satisfy predictions from the TDRL model. Our results, together with comparative data for the available vertebrate complete mt genomes, provide further evidence that the WANCY genomic region is a hotspot for gene order rearrangements and support the view that TDRL is the dominant mechanism of gene order rearrangement in vertebrate mt genomes. Convergent gene rearrangements are not unlikely in hotspots of gene order rearrangement by TDRL.     

Organization and evolution of a gene cluster for human immunoglobulin variable regions of the kappa type

An 80,000 base-pair region from the gene locus encoding the variable regions of the human immunoglobulins of the kappa type (V kappa genes) was cloned and analysed. The region comprises five V kappa sequences of subgroup I and one interspersed V kappa pseudogene of subgroup II. The six genes and pseudogenes are arranged at different distances but in the same orientation. The organization of the cluster can be explained by a series of amplification steps; the existence of a V kappa II pseudogene in a V kappa I gene cluster may have been the result of a transposition event; a final duplication step led to a second closely related copy of the cluster. From sequence data for altogether 16,000 base-pairs it appears that gene conversion-like events and subsequent selection contribute to both homogeneity and diversity of the V kappa repertoire.