Duplication of the gamma-globin gene mediated by L1 long interspersed repetitive elements in an early ancestor of simian primates

Varying strength of selection contributes to the intragenomic diversity of rRNA genes

Ribosome biogenesis in eukaryotes is supported by hundreds of ribosomal RNA (rRNA) gene copies that are encoded in the ribosomal DNA (rDNA). The multiple copies of rRNA genes are thought to have low sequence diversity within one species. Here, we present species-wide rDNA sequence analysis in Saccharomyces cerevisiae that challenges this view. We show that rDNA copies in this yeast are heterogeneous, both among and within isolates, and that many variants avoided fixation or elimination over evolutionary time. The sequence diversity landscape across the rDNA shows clear functional stratification, suggesting different copy-number thresholds for selection that contribute to rDNA diversity. Notably, nucleotide variants in the most conserved rDNA regions are sufficiently deleterious to exhibit signatures of purifying selection even when present in only a small fraction of rRNA gene copies. Our results portray a complex evolutionary landscape that shapes rDNA sequence diversity within a single species and reveal unexpectedly strong purifying selection of multi-copy genes.

Gene fusion of heterophyletic gamma-globin genes in platyrrhine primates

We performed phylogenetic analyses of HBG genes to assess its origin and interspecific variation among primates. Our analyses showed variation in HBG genes copy number ranging from one to three, some of them pseudogenes. For platyrrhines HBG genes, phylogenetic reconstructions of flanking regions recovered orthologous clades with distinct topologies for 5' and 3' flanking regions. The 5' region originated in the common ancestor of platyrrhines but the 3' region had an anthropoid origin. We hypothesize that the platyrrhine HBG genes of 5' and 3' heterophyletic origins arose from subsequent fusions of the (earlier) platyrrhine 5' portion and the (later) anthropoid 3' portion.

β2-microglobulin gene duplication in cetartiodactyla remains intact only in pigs and possibly confers selective advantage to the species

Several β2-microglobulin (B2M) -bound protein complexes undertake key roles in various immune system pathways, including the neonatal Fc receptor (FcRn), cluster of differentiation 1 (CD1) protein, non-classical major histocompatibility complex (MHC), and well-known MHC class I molecules. Therefore, the duplication of B2M may lead to an increase in the biological competence of organisms to the environment. Based on the pig genome assembly SSC10.2, a segmental duplication of ~45.5 kb, encoding the entire B2M protein, was identified in pig chromosome 1. Through experimental validation, we confirmed the functional duplication of the B2M gene with a completely identical coding sequence between two copies in pigs. Considering the importance of B2M in the immune system, we performed the phylogenetic analysis of B2M duplication in ten mammalian species, confirming the presence of B2M duplication in cetartioldactyls, like cattle, sheep, goats, pigs and whales, but non-cetartiodactyl species, like mice, cats, dogs, horses, and humans. The density of long interspersed nuclear element (LINE) at the edges of duplicated blocks (39 to 66%) was found to be 2 to 3-fold higher than the average (20.12%) of the pig genome, suggesting its role in the duplication event. The B2M mRNA expression level in pigs was 12.71 and 7.57 times (2-ΔΔCt values) higher than humans and mice, respectively. However, we were unable to experimentally demonstrate the difference in the level of B2M protein because species specific anti-B2M antibodies are not available. We reported, for the first time, the functional duplication of the B2M gene in animals. The identification of partially remaining duplicated B2M sequences in the genomes of only cetartiodactyls indicates that the event was lineage specific. B2M duplication could be beneficial to the immune system of pigs by increasing the availability of MHC class I light chain protein, B2M, to complex with the proteins encoded by the relatively large number of MHC class I heavy chain genes in pigs. Further studies are necessary to address the biological meaning of increased expression of B2M.

The Role of Retrotransposons in Gene Family Expansions in the Human and Mouse Genomes

Retrotransposons comprise a large portion of mammalian genomes. They contribute to structural changes and more importantly to gene regulation. The expansion and diversification of gene families have been implicated as sources of evolutionary novelties. Given the roles retrotransposons play in genomes, their contribution to the evolution of gene families warrants further exploration. In this study, we found a significant association between two major retrotransposon classes, LINEs and LTRs, and lineage-specific gene family expansions in both the human and mouse genomes. The distribution and diversity differ between LINEs and LTRs, suggesting that each has a distinct involvement in gene family expansion. LTRs are associated with open chromatin sites surrounding the gene families, supporting their involvement in gene regulation, whereas LINEs may play a structural role promoting gene duplication. Our findings also suggest that gene family expansions, especially in the mouse genome, undergo two phases. The first phase is characterized by elevated deposition of LTRs and their utilization in reshaping gene regulatory networks. The second phase is characterized by rapid gene family expansion due to continuous accumulation of LINEs and it appears that, in some instances at least, this could become a runaway process. We provide an example in which this has happened and we present a simulation supporting the possibility of the runaway process. Altogether we provide evidence of the contribution of retrotransposons to the expansion and evolution of gene families. Our findings emphasize the putative importance of these elements in diversification and adaptation in the human and mouse lineages.

Drawing a fine line on endogenous retroelement activity

Endogenous retroelements (EREs) are essential motors of evolution yet require careful control to prevent genomic catastrophes, notably during the vulnerable phases of epigenetic reprogramming that occur immediately after fertilization and in germ cells. Accordingly, a variety of mechanisms restrict these mobile genetic units. Previous studies have revealed the importance of KRAB-containing zinc finger proteins (KRAB-ZFPs) and their cofactor, KAP1, in the early embryonic silencing of endogenous retroviruses and so-called SVAs, but the implication of this transcriptional repression system in the control of LINE-1, the only known active autonomous retrotransposon in the human genome, was thought to be marginal. Two recent studies straighten the record by revealing that the KRAB/KAP system is key to the control of L1 in embryonic stem (ES) cells, and go further in demonstrating that DNA methylation and KRAB/KAP1-induced repression contribute to this process in an evolutionally dynamic fashion. These results shed light on the delicate equilibrium between higher vertebrates and endogenous retroelements, which are not just genetic invaders calling for strict control but rather a constantly renewed and nicely exploitable source of evolutionary potential.

The role of retrotransposons in gene family expansions: insights from the mouse Abp gene family

Background

Retrotransposons have been suggested to provide a substrate for non-allelic homologous recombination (NAHR) and thereby promote gene family expansion. Their precise role, however, is controversial. Here we ask whether retrotransposons contributed to the recent expansions of the Androgen-binding protein (Abp) gene families that occurred independently in the mouse and rat genomes.

Results

Using dot plot analysis, we found that the most recent duplication in the Abp region of the mouse genome is flanked by L1Md_T elements. Analysis of the sequence of these elements revealed breakpoints that are the relicts of the recombination that caused the duplication, confirming that the duplication arose as a result of NAHR using L1 elements as substrates. L1 and ERVII retrotransposons are considerably denser in the Abp regions than in one Mb flanking regions, while other repeat types are depleted in the Abp regions compared to flanking regions. L1 retrotransposons preferentially accumulated in the Abp gene regions after lineage separation and roughly followed the pattern of Abp gene expansion. By contrast, the proportion of shared vs. lineage-specific ERVII repeats in the Abp region resembles the rest of the genome.

Conclusions

We confirmed the role of L1 repeats in Abp gene duplication with the identification of recombinant L1Md_T elements at the edges of the most recent mouse Abp gene duplication. High densities of L1 and ERVII repeats were found in the Abp gene region with abrupt transitions at the region boundaries, suggesting that their higher densities are tightly associated with Abp gene duplication. We observed that the major accumulation of L1 elements occurred after the split of the mouse and rat lineages and that there is a striking overlap between the timing of L1 accumulation and expansion of the Abp gene family in the mouse genome. Establishing a link between the accumulation of L1 elements and the expansion of the Abp gene family and identification of an NAHR-related breakpoint in the most recent duplication are the main contributions of our study.

Molecular cytogenetic characterization of the Amazon River dolphin Inia geoffrensis

Classical and molecular cytogenetic (18S rDNA, telomeric sequence, and LINE-1 retrotransposon probes) studies were carried out to contribute to an understanding of the organization of repeated DNA elements in the Amazon River dolphin (boto, Inia geoffrensis). Twenty-seven specimens were examined, each presenting 2n = 44 chromosomes, the karyotype formula 12m + 14sm + 6st + 10t + XX/XY, and fundamental number (FN) = 74. C-positive heterochromatin was observed in terminal and interstitial positions, with the occurrence of polymorphism. Interstitial telomeric sequences were not observed. The nucleolar organizer region (NOR) was located at a single site on a smallest autosomal pair. LINE-1 was preferentially distributed in the euchromatin regions, with the greatest accumulation on the X chromosome. Although the karyotype structure in cetaceans is considered to be conserved, the boto karyotype demonstrated significant variations in its formula, heterochromatin distribution, and the location of the NOR compared to other cetacean species. These results contribute to knowledge of the chromosome organization in boto and to a better understanding of karyoevolution in cetaceans.

Revealing mammalian evolutionary relationships by comparative analysis of gene clusters

Many software tools for comparative analysis of genomic sequence data have been released in recent decades. Despite this, it remains challenging to determine evolutionary relationships in gene clusters due to their complex histories involving duplications, deletions, inversions, and conversions. One concept describing these relationships is orthology. Orthologs derive from a common ancestor by speciation, in contrast to paralogs, which derive from duplication. Discriminating orthologs from paralogs is a necessary step in most multispecies sequence analyses, but doing so accurately is impeded by the occurrence of gene conversion events. We propose a refined method of orthology assignment based on two paradigms for interpreting its definition: by genomic context or by sequence content. X-orthology (based on context) traces orthology resulting from speciation and duplication only, while N-orthology (based on content) includes the influence of conversion events. We developed a computational method for automatically mapping both types of orthology on a per-nucleotide basis in gene cluster regions studied by comparative sequencing, and we make this mapping accessible by visualizing the output. All of these steps are incorporated into our newly extended CHAP 2 package. We evaluate our method using both simulated data and real gene clusters (including the well-characterized α-globin and β-globin clusters). We also illustrate use of CHAP 2 by analyzing four more loci: CCL (chemokine ligand), IFN (interferon), CYP2abf (part of cytochrome P450 family 2), and KIR (killer cell immunoglobulin-like receptors). These new methods facilitate and extend our understanding of evolution at these and other loci by adding automated accurate evolutionary inference to the biologist's toolkit. The CHAP 2 package is freely available from http://www.bx.psu.edu/miller_lab.

Transposable elements as tools for reshaping the genome: it is a huge world after all!

Transposable elements (TEs) are discrete pieces of DNA that can move from one site to another within genomes and sometime between genomes. They are found in all major branches of life. Because of their wide distribution and considerable diversity, they are a considerable source of genomic variation and as such, they constitute powerful drivers of genome evolution. Moreover, it is becoming clear that the epigenetic regulation of certain genes is derived from defense mechanisms against the activity of ancestral transposable elements. TEs now tend to be viewed as natural molecular tools that can reshape the genome, which challenges the idea that TEs are natural tools used to answer biological questions. In the first part of this chapter, we review the classification and distribution of TEs, and look at how they have contributed to the structural and transcriptional reshaping of genomes. In the second part, we describe methodological innovations that have modified their contribution as molecular tools.

Conversion events in gene clusters

BACKGROUND

Gene clusters containing multiple similar genomic regions in close proximity are of great interest for biomedical studies because of their associations with inherited diseases. However, such regions are difficult to analyze due to their structural complexity and their complicated evolutionary histories, reflecting a variety of large-scale mutational events. In particular, conversion events can mislead inferences about the relationships among these regions, as traced by traditional methods such as construction of phylogenetic trees or multi-species alignments.

RESULTS

To correct the distorted information generated by such methods, we have developed an automated pipeline called CHAP (Cluster History Analysis Package) for detecting conversion events. We used this pipeline to analyze the conversion events that affected two well-studied gene clusters (α-globin and β-globin) and three gene clusters for which comparative sequence data were generated from seven primate species: CCL (chemokine ligand), IFN (interferon), and CYP2abf (part of cytochrome P450 family 2). CHAP is freely available at http://www.bx.psu.edu/miller_lab.

CONCLUSIONS

These studies reveal the value of characterizing conversion events in the context of studying gene clusters in complex genomes.

Mobile DNA and the TE-Thrust hypothesis: supporting evidence from the primates

Transposable elements (TEs) are increasingly being recognized as powerful facilitators of evolution. We propose the TE-Thrust hypothesis to encompass TE-facilitated processes by which genomes self-engineer coding, regulatory, karyotypic or other genetic changes. Although TEs are occasionally harmful to some individuals, genomic dynamism caused by TEs can be very beneficial to lineages. This can result in differential survival and differential fecundity of lineages. Lineages with an abundant and suitable repertoire of TEs have enhanced evolutionary potential and, if all else is equal, tend to be fecund, resulting in species-rich adaptive radiations, and/or they tend to undergo major evolutionary transitions. Many other mechanisms of genomic change are also important in evolution, and whether the evolutionary potential of TE-Thrust is realized is heavily dependent on environmental and ecological factors. The large contribution of TEs to evolutionary innovation is particularly well documented in the primate lineage. In this paper, we review numerous cases of beneficial TE-caused modifications to the genomes of higher primates, which strongly support our TE-Thrust hypothesis.

A combined top-down and bottom-up MS approach for the characterization of hemoglobin variants in Rhesus monkeys

Sickle cell disease is caused by one of the 1200 known hemoglobin variations. A single-point mutation β6(A3)Glu→Val leads to sickling of red blood cells, which in turn causes a lack of oxygen supply to tissue and organs. Although sickle cell disease is well understood, treatment options are currently underdeveloped. The only Food and Drug Administration-approved drug is hydroxyurea, an inducer of fetal γ-hemoglobin, which is known to have a higher oxygen affinity than adult hemoglobins and thus alleviates symptoms. In the search for better cures, Rhesus monkeys (Macaca mulatta) serve as models for monitoring success of induction of fetal γ-hemoglobins and with recent advances in proteomics, MS has become the leading technique to determine globin expression. Similar to humans, Rhesus monkeys possess hemoglobin variants that have not been sufficiently characterized to initiate such a study. Therefore, we developed a combined bottom-up and top-down approach to identify and characterize novel hemoglobin variants of the umbilical cord blood of Rhesus monkeys. A total of four different variants were studied: α, β, γ1 and γ2. A new α- and β-hemoglobin variant was identified, and the two previously hypothesized γ-hemoglobins were identified. In addition, glutathionylation of both γ-hemoglobin variants at their cysteines has been characterized. The combined approach outperformed either bottom-up or top-down alone and can be used for characterization of unknown hemoglobin variants and their PTMs.

Genes devoid of full-length transposable element insertions are involved in development and in the regulation of transcription in human and closely related species

Transposable elements (TEs) are major components of mammalian genomes, and their impact on genome evolution is now well established. In recent years several findings have shown that they are associated with the expression level and function of genes. In this study, we analyze the relationships between human genes and full-length TE copies in terms of three factors (gene function, expression level, and selective pressure). We classified human genes according to their TE density, and found that TE-free genes are involved in important functions such as development, transcription, and the regulation of transcription, whereas TE-rich genes are involved in functions such as transport and metabolism. This trend is conserved through evolution. We show that this could be explained by a stronger selection pressure acting on both the coding and non-coding regions of TE-free genes than on those of TE-rich genes. The higher level of expression found for TE-rich genes in tumor and immune system tissues suggests that TEs play an important role in gene regulation.

A LINE-1 component to human aging: do LINE elements exact a longevity cost for evolutionary advantage?

Advancing age remains the largest risk factor for devastating diseases, such as heart disease, stroke, and cancer. The mechanisms by which advancing age predisposes to disease are now beginning to unfold, due in part, to genetic and environmental manipulations of longevity in lower organisms. Converging lines of evidence suggest that DNA damage may be a final common pathway linking several proposed mechanisms of aging. The present review forwards a theory for an additional aging pathway that involves modes of inherent genetic instability. Long interspersed nuclear elements (LINEs) are endogenous non-LTR retrotransposons that compose about 20% of the human genome. The LINE-1 (L1) gene products, ORF1p and ORF2p, possess mRNA binding, endonuclease, and reverse transcriptase activity that enable retrotransposition. While principally active only during embryogenesis, L1 transcripts are detected in adult somatic cells under certain conditions. The present hypothesis proposes that L1s act as an 'endogenous clock', slowly eroding genomic integrity by competing with the organism's double-strand break repair mechanism. Thus, while L1s are an accepted mechanism of genetic variation fueling evolution, it is proposed that longevity is negatively impacted by somatic L1 activity. The theory predicts testable hypotheses about the relationship between L1 activity, DNA repair, healthy aging, and longevity.

High divergence in primate-specific duplicated regions: human and chimpanzee chorionic gonadotropin beta genes

Background

Low nucleotide divergence between human and chimpanzee does not sufficiently explain the species-specific morphological, physiological and behavioral traits. As gene duplication is a major prerequisite for the emergence of new genes and novel biological processes, comparative studies of human and chimpanzee duplicated genes may assist in understanding the mechanisms behind primate evolution. We addressed the divergence between human and chimpanzee duplicated genomic regions by using Luteinizing Hormone Beta (LHB)/Chorionic Gonadotropin Beta (CGB) gene cluster as a model. The placental CGB genes that are essential for implantation have evolved from an ancestral pituitary LHB gene by duplications in the primate lineage.

Results

We shotgun sequenced and compared the human (45,165 bp) and chimpanzee (39,876 bp) LHB/CGB regions and hereby present evidence for structural variation resulting in discordant number of CGB genes (6 in human, 5 in chimpanzee). The scenario of species-specific parallel duplications was supported (i) as the most parsimonious solution requiring the least rearrangement events to explain the interspecies structural differences; (ii) by the phylogenetic trees constructed with fragments of intergenic regions; (iii) by the sequence similarity calculations. Across the orthologous regions of LHB/CGB cluster, substitutions and indels contributed approximately equally to the interspecies divergence and the distribution of nucleotide identity was correlated with the regional repeat content. Intraspecies gene conversion may have shaped the LHB/CGB gene cluster. The substitution divergence (1.8–2.59%) exceeded two-three fold the estimates for single-copy loci and the fraction of transversional mutations was increased compared to the unique sequences (43% versus ~30%). Despite the high sequence identity among LHB/CGB genes, there are signs of functional differentiation among the gene copies. Estimates for d_n/d_s rate ratio suggested a purifying selection on LHB and CGB8, and a positive evolution of CGB1.

Conclusion

If generalized, our data suggests that in addition to species-specific deletions and duplications, parallel duplication events may have contributed to genetic differences separating humans from their closest relatives. Compared to unique genomic segments, duplicated regions are characterized by high divergence promoted by intraspecies gene conversion and species-specific chromosomal rearrangements, including the alterations in gene copy number.

Characterization of the −16C>G sequence variation in the promoters of both HBG1 and HBG2: Convergent evolution of the human γ-globin genes

We encountered a homozygous -16C>G mutation in cis at identical positions in the promoters of both human gamma-globin genes in a subject who was also homozygous for Hemoglobin C (HbC). Subsequent analysis of normal control individuals of African American ancestry revealed that both mutations were always present in cis with an allelic frequency of 3%. Furthermore, 10 out of 11 HbC subjects carried the -16C>G sequence variations, suggesting an association with HbC. The -16C>G mutation disrupts a putative CACCC box positioned between the TATA box and the transcriptional start site. However, the absence of high levels of HbF in HbC subjects homozygous and heterozygous for the -16C>G sequence variation suggested no effect of this mutation on gamma-globin gene expression in the adult stage of development. Further functional characterization by means of transient transfections in human erythroleukemic K562 cells showed that the -16C>G promoter sequence variation did not have an effect on gamma-globin expression in the fetal stage of development either. We therefore conclude that the -16C>G gamma-globin sequence variations are not beneficial to the clinical phenotype of HbC. The unique concurrent presence of this non-functional sequence variation is likely the result of a gene conversion event, and supports the concept of sequence homogenization between the two human gamma-globin genes.

Phylogenetic comparisons suggest that distance from the locus control region guides developmental expression of primate beta-type globin genes

Phylogenetic inferences drawn from comparative data on mammalian beta-globin gene clusters indicate that the ancestral primate cluster contained a locus control region (LCR) and five paralogously related beta-type globin loci (5'-LCR-epsilon-gamma-psieta-delta-beta-3'), with epsilon and gamma expressed solely during embryonic life. A gamma locus tandem duplication (5'-gamma(1)-gamma(2)-3') triggered gamma's evolution toward fetal expression but by a different trajectory in platyrrhines (New World monkeys) than in catarrhines (Old World monkeys and apes, including humans). In platyrrhine (e.g., Cebus) fetuses, gamma(1) at the ancestral distance from epsilon is down-regulated, whereas gamma(2) at increased distance is up-regulated. Catarrhine gamma(1) and gamma(2) acquired longer distances from epsilon (14 and 19 kb, respectively), and both are up-regulated throughout fetal life with gamma(1)'s expression predominating over gamma(2)'s. On enlarging the platyrrhine expression data, we find Aotus gamma is embryonic, Alouatta gamma is inactive at term, and in Callithrix, gamma(1) is down-regulated fetally, whereas gamma(2) is up-regulated. Of eight mammalian taxa now represented per taxon by embryonic, fetal, and postnatal beta-type globin gene expression data, four taxa are primates, and data for three of these primates are from this laboratory. Our results support a model in which a short distance (<10 kb) between epsilon and the adjacent gamma is a plesiomorphic character that allows the LCR to drive embryonic expression of both genes, whereas a longer distance (>10 kb) impedes embryonic activation of the downstream gene.

Repositioning-dependent fate of duplicate genes

Gene duplication is the main source of evolutionary novelties. However, the problem with duplicates is that the purifying selection overlooks deleterious mutations in the redundant sequence, which therefore, instead of gaining a new function, often degrades into a functionless pseudogene. This risk of functional loss instead of gain is much higher for small populations of higher organisms with a slow and complex development. We propose that it is the epigenetic tissue/stage-complementary silencing of duplicates that makes them exposable to the purifying selection, thus saving them from pseudogenization and opening the way towards new function(s). Our genome-wide analyses of gene duplicates in several eukaryotic species combined with the phylogenetic comparison of vertebrate alpha- and beta-globin gene clusters strongly support this epigenetic complementation (EC) model. The distinctive condition for a new duplicate to survive by the EC mechanism seems to be its repositioning to an ectopic site, which is accompanied by changes in the rate and direction of mutagenesis. The most distinguished in this respect is the human genome. In this review, we extend and discuss the data on the EC- and repositioning-dependent fate of gene duplicates with the special emphasis on the problem of detecting brief postduplication period of adaptive evolution driven by positive selection. Accordingly, we propose a new CpG-focused measure of selection that is insensitive to translocation-caused biases in mutagenesis.

The genomic distribution of L1 elements: the role of insertion bias and natural selection

LINE-1 (L1) retrotransposons constitute the most successful family of retroelements in mammals and account for as much as 20% of mammalian DNA. L1 elements can be found in all genomic regions but they are far more abundant in AT-rich, gene-poor, and low-recombining regions of the genome. In addition, the sex chromosomes and some genes seem disproportionately enriched in L1 elements. Insertion bias and selective processes can both account for this biased distribution of L1 elements. L1 elements do not appear to insert randomly in the genome and this insertion bias can at least partially explain the genomic distribution of L1. The contrasted distribution of L1 and Alu elements suggests that postinsertional processes play a major role in shaping L1 distribution. The most likely mechanism is the loss of recently integrated L1 elements that are deleterious (negative selection) either because of disruption of gene function or their ability to mediate ectopic recombination. By comparison, the retention of L1 elements because of some positive effect is limited to a small fraction of the genome. Understanding the respective importance of insertion bias and selection will require a better knowledge of insertion mechanisms and the dynamics of L1 inserts in populations.

Gene conversion and functional divergence in the beta-globin gene family

Different models of gene family evolution have been proposed to explain the mechanism whereby gene copies created by gene duplications are maintained and diverge in function. Ohta proposed a model which predicts a burst of nonsynonymous substitutions following gene duplication and the preservation of duplicates through positive selection. An alternative model, the duplication-degeneration-complementation (DDC) model, does not explicitly require the action of positive Darwinian selection for the maintenance of duplicated gene copies, although purifying selection is assumed to continue to act on both copies. A potential outcome of the DDC model is heterogeneity in purifying selection among the gene copies, due to partitioning of subfunctions which complement each other. By using the d(N)/ d(S) (omega) rate ratio to measure selection pressure, we can distinguish between these two very different evolutionary scenarios. In this study we investigated these scenarios in the beta-globin family of genes, a textbook example of evolution by gene duplication. We assembled a comprehensive dataset of 72 vertebrate beta-globin sequences. The estimated phylogeny suggested multiple gene duplication and gene conversion events. By using different programs to detect recombination, we confirmed several cases of gene conversion and detected two new cases. We tested evolutionary scenarios derived from Ohta's model and the DDC model by examining selective pressures along lineages in a phylogeny of beta-globin genes in eutherian mammals. We did not find significant evidence for an increase in the omega ratio following major duplication events in this family. However, one exception to this pattern was the duplication of gamma-globin in simian primates, after which a few sites were identified to be under positive selection. Overall, our results suggest that following gene duplications, paralogous copies of beta-globin genes evolved under a nonepisodic process of functional divergence.

Characterization of in vivo recombination activities in the mouse embryo

Homologous recombination makes use of sequence homology to repair DNA and to rearrange genetic material. In mammals, these processes have mainly been characterized using cultured cell systems. We have developed an assay that allows us to quantitatively analyze homologous recombination in vivo in the mouse embryo. Transgenic mouse lines were generated by microinjection into a fertilized mouse ovum of a vector containing two homologous LINE-1 (L1) sequences arranged as a direct repeat: these sequences can recombine with each other and with endogenous L1 sequences before, during or after integration of the vector into the genome. Using a plasmid rescue procedure, we determined the composition of the integrated vector array in several transgenic mice and their descendants. Homologous recombination frequencies were found to be strikingly high, involving 70% of integrated vectors in some arrays, with homologous deletions being five times more frequent than gene conversion without crossing-over. Interestingly, non-homologous recombination was found to be much less frequent. We also found that endogenous L1 sequences could be involved in homologous recombination events in the mouse embryo, and that the integrated arrays could be modified from generation to generation by homologous recombination between the integrated L1 sequences.

Co-duplication of olfactory receptor and MHC class I genes in the mouse major histocompatibility complex

We report the 897 kb sequence of a cluster of olfactory receptor (OR) genes located at the distal end of the major histocompatibility complex (MHC) class I region on mouse chromosome 17 of strain 129/SvJ (H2bc). With additional information from the mouse genome draft sequence, we identified 59 OR loci (approximately 20% pseudogenes) in contrast to only 25 OR loci (approximately 50% pseudogenes) in the corresponding centromeric OR cluster that is part of the 'extended MHC class I region' on human chromosome 6. Comparative analysis leads to three major observations: (i) most of the OR subfamilies have evolved independently in the two species, expanding more in the mouse, and resulting in co-orthologs--subfamilies of highly similar paralogs that keep orthologous relationships with their human counterparts; (ii) three of the mouse OR subfamilies have no orthologs in humans; and (iii) MHC class I loci are interspersed in the OR cluster in mouse but not in human, and were subjected to co-duplication with OR genes. Screening of our sequence against the available sequences of other strains/haplotypes revealed that most of the OR loci are polymorphic and that the number of OR loci may vary among strains/haplotypes. Our findings that MHC-linked OR loci share duplication with MHC class I loci, have duplicated extensively and are polymorphic revives questions about potential reciprocal influences acting on the dynamics and evolution of the H2 region and the H2-linked OR loci.

An overview of the kallikrein gene families in humans and other species: emerging candidate tumour markers

Kallikreins are serine proteases with diverse physiologic functions. They are represented by multigene families in many animal species, especially in rat and mouse. Recently, the human kallikrein gene family has been fully characterized and includes 15 members, tandemly localized on chromosome 19q13.4. A new definition has now been proposed for kallikreins, which is not based on function but, rather, on close proximity and structural similarities. In this review, we summarize available information about kallikreins in many animal species with special emphasis on human kallikreins. We discuss the common structural features of kallikreins at the DNA, mRNA and protein levels and overview their evolutionary history. Kallikreins are expressed in a wide range of tissues including the salivary gland, endocrine or endocrine-related tissues such as testis, prostate, breast and endometrium and in the central nervous system. Most, if not all, genes are under steroid hormone regulation. Accumulating evidence indicates that kallikreins are involved in many pathologic conditions. Of special interest is the potential role of kallikreins in the central nervous system. In addition, many kallikreins seem to be candidate tumor markers for many malignancies, especially those of endocrine-related organs.

Atypical expansion in mice of the sensory neuron-specific Mrg G protein-coupled receptor family

The Mas-related genes (Mrgs) comprise a family of >50 G protein-coupled receptors (GPCRs), many of which are expressed in specific subsets of nociceptive sensory neurons in mice. In contrast, humans contain a related but nonorthologous family of genes, called MrgXs or sensory neuron-specific receptors, of which many fewer appear to be expressed in sensory neurons. To determine whether the diversity of murine Mrgs is generic to rodents or is an atypical feature of mice, we characterized MrgA, MrgB, MrgC, and MrgD subfamilies in rat and gerbil. Surprisingly, although mice have approximately 22 MrgA and approximately 14 MrgC genes, rats and gerbils have just a single MrgA and MrgC gene. This murine-specific expansion likely reflects recent retrotransposon-mediated unequal crossover events. The expression of Mrgs in rat sensory ganglia suggests that the extensive cellular diversity in mice can be simplified to a core subset of approximately four different genes (MrgA, MrgB, MrgC, and MrgD), defining a similar number of neuronal subpopulations. Our results suggest more generally that mouse-human genomic comparisons may sometimes reveal differences atypical of rodents.

Type III intermediate filament proteins interact with four-way junction DNA and facilitate its cleavage by the junction-resolving enzyme T7 endonuclease I

The isolation from proliferating mouse and human embryo fibroblasts of SDS-stable crosslinkage products of vimentin with DNA fragments containing inverted repeats capable of cruciform formation under superhelical stress and the competitive effect of a synthetic Holliday junction on the binding of cytoplasmic intermediate filament (cIF) proteins to supercoiled DNA prompted a detailed investigation of the proteins' capacity to associate with four-way junction DNA and to influence its processing by junction-resolving endonucleases. Electrophoretic mobility shift analysis of reaction products obtained from vimentin and Holliday junctions under varying ionic conditions revealed efficient complex formation of the filament protein not only with the unstacked, square-planar configuration of the junctions but also with their coaxially stacked X-conformation. Glial fibrillary acidic protein (GFAP) was less efficient and desmin virtually inactive in complex formation. Electron microscopy showed binding of vimentin tetramers or octamers almost exclusively to the branchpoint of the Holliday junctions under physiological ionic conditions. Even at several hundredfold molar excess, sequence-related single- and double-stranded DNAs were unable to chase Holliday junctions from their complexes with vimentin. Vimentin also stimulated bacteriophage T7 endonuclease I in introducing single-strand cuts diametrically across the branchpoint and thus in the resolution of the Holliday junctions. This effect is very likely due to vimentin-induced structural distortion of the branchpoint, as suggested by the results of hydroxyl radical footprinting of Holliday junctions in the absence and the presence of vimentin. Moreover, vimentin, and to a lesser extent GFAP and desmin, interacted with the cruciform structures of inverted repeats inserted into a supercoiled vector plasmid, thereby changing their configuration via branch migration and sensibilizing them to processing by T7 endonuclease I. This refers to both plasmid relaxation caused by unilateral scission and, particularly, linearization via bilateral scission at primary and cIF protein-induced secondary cruciform branchpoints that were identified by T7 endonuclease I footprinting. cIF proteins share these activities with a variety of other architectural proteins interacting with and structurally modulating four-way DNA junctions. In view of the known and hypothetical functions of four-way DNA junctions and associated protein factors in DNA metabolism, cIF proteins as complementary nuclear matrix proteins may play important roles in such nuclear matrix-associated processes as DNA replication, recombination, repair, and transcription, with special emphasis on both the preservation and evolution of the genome.

LINE-1 preTa elements in the human genome

The preTa subfamily of long interspersed elements (LINEs) is characterized by a three base-pair "ACG" sequence in the 3' untranslated region, contains approximately 400 members in the human genome, and has low level of nucleotide divergence with an estimated average age of 2.34 million years old suggesting that expansion of the L1 preTa subfamily occurred just after the divergence of humans and African apes. We have identified 362 preTa L1 elements from the draft human genomic sequence, investigated the genomic characteristics of preTa L1 insertions, and screened individual elements across diverse human populations and various non-human primate species using polymerase chain reaction (PCR) assays to determine the phylogenetic origin and levels of human genomic diversity associated with the L1 elements. All of the preTa L1 elements analyzed by PCR were absent from the orthologous positions in non-human primate genomes with 33 (14%) of the L1 elements being polymorphic with respect to insertion presence or absence in the human genome. The newly identified L1 insertion polymorphisms will prove useful as identical by descent genetic markers for the study of human population genetics. We provide evidence that preTa L1 elements show an integration site preference for genomic regions with low GC content. Computational analysis of the preTa L1 elements revealed that 29% of the elements amenable to complete sequence analysis have apparently escaped 5' truncation and are essentially full-length (approximately 6kb). In all, 29 have two intact open reading frames and may be capable of retrotransposition.

Globin gene switching in primates

Evolutionary approaches to the identification of DNA sequences required for transcription of the genes of the beta-globin cluster are reviewed. Sequence alignments of non-coding regions from widely divergent species revealed many conserved motifs (phylogenetic footprints) that were putative transcription factor binding sites and candidate binding proteins were identified. The differential timing of the prosimian and simian gamma-globin genes was analyzed by identifying base changes in the vicinity of the phylogenetic footprints. These differential phylogenetic footprints were shown to bind different nuclear factors, and the behavior of constructs with human or galago gamma-promoters in transgenic mice indicated that DNA motifs near the gamma-globin genes are sufficient to determine the developmental stage of expression. Locus control region alignments have identified many conserved sequence differences outside of the hypersensitive sites. Globin protein and mRNA expression profiles during embryological development in a series of catarrhine (Old World monkeys and apes) and platyrrhine (New World monkeys) primates have been determined. While all catarrhines examined to date have globin expression patterns that are highly similar to the well-established human switching behavior, platyrrhines have inactivated their gamma 1 genes by a variety of mechanisms, and have an earlier gamma-beta switch.

A comprehensive analysis of recently integrated human Ta L1 elements

The Ta (transcribed, subset a) subfamily of L1 LINEs (long interspersed elements) is characterized by a 3-bp ACA sequence in the 3' untranslated region and contains approximately 520 members in the human genome. Here, we have extracted 468 Ta L1Hs (L1 human specific) elements from the draft human genomic sequence and screened individual elements using polymerase-chain-reaction (PCR) assays to determine their phylogenetic origin and levels of human genomic diversity. One hundred twenty-four of the elements amenable to complete sequence analysis were full length ( approximately 6 kb) and have apparently escaped any 5' truncation. Forty-four of these full-length elements have two intact open reading frames and may be capable of retrotransposition. Sequence analysis of the Ta L1 elements showed a low level of nucleotide divergence with an estimated age of 1.99 million years, suggesting that expansion of the L1 Ta subfamily occurred after the divergence of humans and African apes. A total of 262 Ta L1 elements were screened with PCR-based assays to determine their phylogenetic origin and the level of human genomic variation associated with each element. All of the Ta L1 elements analyzed by PCR were absent from the orthologous positions in nonhuman primate genomes, except for a single element (L1HS72) that was also present in the common (Pan troglodytes) and pygmy (P. paniscus) chimpanzee genomes. Sequence analysis revealed that this single exception is the product of a gene conversion event involving an older preexisting L1 element. One hundred fifteen (45%) of the Ta L1 elements were polymorphic with respect to insertion presence or absence and will serve as identical-by-descent markers for the study of human evolution.

Genomic organization of the siglec gene locus on chromosome 19q13.4 and cloning of two new siglec pseudogenes

The sialic acid binding immunoglobulin-like lectin (Siglec) family of genes is a recently described member of the immunoglobulin superfamily. Within this Siglec family there is a subgroup of genes which bear a high degree of homology with Siglec-3 (CD33), thus designated the Siglec-3-like subgroup of Siglecs. While their mRNA structure has been reported, the full genomic organization of these genes, is not known. Genes of this subgroup have been mapped to chromosome 19q13.4, primarily through in situ hybridization. Through analysis of several bacterial artificial chromosome (BAC) clones, we studied an approximate 700 kb region that encompasses the putative Siglec gene locus on chromosome 19q13.4. We established the first detailed map of the locus, which contains 8 Siglec and Siglec-like genes. Our map shows the relative position of all genes and the precise distances between them, along with the direction of transcription of each gene. To our knowledge, this is the first report that describes the full genomic organization of all members of the CD33-like subgroup of Siglecs, including the promoter sequences of all genes. Members of this subfamily exhibit two patterns of organization of the signal peptide, which is followed by one V-set domain (except for the long form of the siglecL1 gene). Exons containing the C2-set domains are all comparable in size and are separated by linker exons. The transmembrane domain is encoded for by a separate exon of almost the same size in all genes. The total number of exons differs according to the number of C2-set Ig domains, but intron phases are identical. The cytoplasmic domain is always encoded by two exons. We further identified two new Siglec pseudogenes in this locus, and analyzed their tissue expression pattern and their structural features.

Genomic structure of newly identified paralogue of RNA helicase II/Gu: detection of pseudogenes and multiple alternatively spliced mRNAs

RNA helicase II/Gu (RH-II/Gu or DDX21) is a DEAD-box enzyme that localizes to the nucleoli and may be involved in ribosomal RNA synthesis or processing. It has two paralogues, RH-II/Gualpha and RH-II/Gubeta, both genes of which are on chromosome 10. Their similar genomic structures suggest the two genes arose by gene duplication. Both genes are expressed at higher levels in some normal human tissues compared to matching tumor tissues. Pseudogenes for RH-II/Gubeta exist on chromosomes 2, 3 and 4. No pseudogene was identified for RH-II/Gualpha. Both exon inclusion and exon skipping were found to post-transcriptionally regulate RH-II/Gubeta gene expression. No alternative splicing was identified for RH-II/Gualpha. Overall, the results suggest that the two paralogues of RH-II/Gu arose by gene duplication but the resulting genes are differentially regulated.

Biology of mammalian L1 retrotransposons

L1 retrotransposons comprise 17% of the human genome. Although most L1s are inactive, some elements remain capable of retrotransposition. L1 elements have a long evolutionary history dating to the beginnings of eukaryotic existence. Although many aspects of their retrotransposition mechanism remain poorly understood, they likely integrate into genomic DNA by a process called target primed reverse transcription. L1s have shaped mammalian genomes through a number of mechanisms. First, they have greatly expanded the genome both by their own retrotransposition and by providing the machinery necessary for the retrotransposition of other mobile elements, such as Alus. Second, they have shuffled non-L1 sequence throughout the genome by a process termed transduction. Third, they have affected gene expression by a number of mechanisms. For instance, they occasionally insert into genes and cause disease both in humans and in mice. L1 elements have proven useful as phylogenetic markers and may find other practical applications in gene discovery following insertional mutagenesis in mice and in the delivery of therapeutic genes.

Selection against deleterious LINE-1-containing loci in the human lineage

We compared sex chromosomal and autosomal regions of similar GC contents and found that the human Y chromosome contains nine times as many full-length (FL) ancestral LINE-1 (L1) elements per megabase as do autosomes and that the X chromosome contains three times as many. In addition, both sex chromosomes contain a ca. twofold excess of elements that are >500 bp but not long enough to be capable of autonomous replication. In contrast, the autosomes are not deficient in short (<500 bp) L1 elements or SINE elements relative to the sex chromosomes. Since neither the Y nor the X chromosome, when present in males, can be cleared of deleterious genetic loci by recombination, we conclude that most FL L1s were deleterious and thus subject to purifying selection. Comparison between nonrecombining and recombining regions of autosome 21 supported this conclusion. We were able to identify a subset of loci in the human DNA database that once contained active L1 elements, and we found by using the polymerase chain reaction that 72% of them no longer contain L1 elements in a representative of each of eight different ethnic groups. Genetic damage produced by both L1 retrotransposition and ectopic (nonallelic) recombination between L1 elements could provide the basis for their negative selection.

Characterization of clustered MHC-linked olfactory receptor genes in human and mouse

Olfactory receptor (OR) loci frequently cluster and are present on most human chromosomes. They are members of the seven transmembrane receptor (7-TM) superfamily and, as such, are part of one of the largest mammalian multigene families, with an estimated copy number of up to 1000 ORs per haploid genome. As their name implies, ORs are known to be involved in the perception of odors and possibly also in other, nonolfaction-related, functions. Here, we report the characterization of ORs that are part of the MHC-linked OR clusters in human and mouse (partial sequence only). These clusters are of particular interest because of their possible involvement in olfaction-driven mate selection. In total, we describe 50 novel OR loci (36 human, 14 murine), making the human MHC-linked cluster the largest sequenced OR cluster in any organism so far. Comparative and phylogenetic analyses confirm the cluster to be MHC-linked but divergent in both species and allow the identification of at least one ortholog that will be useful for future regulatory and functional studies. Quantitative feature analysis shows clear evidence of duplications of blocks of OR genes and reveals the entire cluster to have a genomic environment that is very different from its neighboring regions. Based on in silico transcript analysis, we also present evidence of extensive long-distance splicing in the 5'-untranslated regions and, for the first time, of alternative splicing within the single coding exon of ORs. Taken together with our previous finding that ORs are also polymorphic, the presented data indicate that the expression, function, and evolution of these interesting genes might be more complex than previously thought.

Characterization of Clustered MHC-Linked Olfactory Receptor Genes in Human and Mouse

Olfactory receptor (OR) loci frequently cluster and are present on most human chromosomes. They are members of the seven transmembrane receptor (7-TM) superfamily and, as such, are part of one of the largest mammalian multigene families, with an estimated copy number of up to 1000 ORs per haploid genome. As their name implies, ORs are known to be involved in the perception of odors and possibly also in other, nonolfaction-related, functions. Here, we report the characterization of ORs that are part of the MHC-linked OR clusters in human and mouse (partial sequence only). These clusters are of particular interest because of their possible involvement in olfaction-driven mate selection. In total, we describe 50 novel OR loci (36 human, 14 murine), making the human MHC-linked cluster the largest sequenced OR cluster in any organism so far. Comparative and phylogenetic analyses confirm the cluster to be MHC-linked but divergent in both species and allow the identification of at least one ortholog that will be useful for future regulatory and functional studies. Quantitative feature analysis shows clear evidence of duplications of blocks of OR genes and reveals the entire cluster to have a genomic environment that is very different from its neighboring regions. Based on in silico transcript analysis, we also present evidence of extensive long-distance splicing in the 5′-untranslated regions and, for the first time, of alternative splicing within the single coding exon of ORs. Taken together with our previous finding that ORs are also polymorphic, the presented data indicate that the expression, function, and evolution of these interesting genes might be more complex than previously thought.

[The sequence data described in this paper have been submitted to the EMBL nucleotide data library under accession nos.Z84475, Z98744, Z98745, AL021807, AL021808, AL022723, AL022727,AL031893, AL035402, AL035542, AL050328, AL050339, AL078630, AL096770,AL121944, AL133160, and AL133267.]

Gene conversion may aid adaptive peak shifts

Gene conversion is often viewed as a homogenizing force that opposes adaptive evolution. The objective of this study is to suggest a potential role for gene conversion in adaptive evolution of proteins through aiding the transfer of a population from one adaptive peak to another. Our hypothesis starts with the observation that a tandem gene duplication may result in an extra gene copy that is released from selective constraints. In such cases, individually deleterious mutations may accumulate on the extra copy of the gene, and through gene conversion these mutations may subsequently be presented to the functioning gene for selection en masse. Thus, groups of mutations that jointly confer a selective advantage may regularly be made available for selection. We present a mathematical model of this process and identify the range of rates of gene conversion, gene duplication and mutation under which it may operate. The results indicate that the process may be biologically feasible if the rate of appearance of the potentially beneficial mutations is not too small in relation to the rates of null mutation and of gene conversion. This process appears to be a possible mechanism for effecting adaptive peak shifts in large populations. We show that all the evolutionary steps in the proposed model may have occurred in the evolution of primate gamma -globin genes. We suggest that hide-and-release mechanisms for genetic variation may constitute a more general principal of evolvability.

Genetics. L1 retrotransposons shape the mammalian genome

What are all of those retrotransposons doing buzzing about in our genome? According to Kazazian, in his Perspective this week, these mobile pieces of DNA are busy reshaping our genome, making it more diverse and enabling us to survive and thrive through the vagaries of evolution. And just how do they do this?...zip to page 1152 to find out.

Cloning, characterization, and phylogenetic analysis of siglec-9, a new member of the CD33-related group of siglecs. Evidence for co-evolution with sialic acid synthesis pathways

The Siglecs are a subfamily of I-type lectins (immunoglobulin superfamily proteins that bind sugars) that specifically recognize sialic acids. We report the cloning and characterization of human Siglec-9. The cDNA encodes a type 1 transmembrane protein with three extracellular immunoglobulin-like domains and a cytosolic tail containing two tyrosines, one within a typical immunoreceptor tyrosine-based inhibitory motif (ITIM). The N-terminal V-set Ig domain has most amino acid residues typical of Siglecs. Siglec-9 is expressed on granulocytes and monocytes. Expression of the full-length cDNA in COS cells induces sialic-acid dependent erythrocyte binding. A recombinant soluble form of the extracellular domain binds to alpha2-3 and alpha2-6-linked sialic acids. Typical of Siglecs, the carboxyl group and side chain of sialic acid are essential for recognition, and mutation of a critical arginine residue in domain 1 abrogates binding. The underlying glycan structure also affects binding, with Galbeta1-4Glc[NAc] being preferred. Siglec-9 shows closest homology to Siglec-7 and both belong to a Siglec-3/CD33-related subset of Siglecs (with Siglecs-5, -6, and -8). The Siglec-9 gene is on chromosome 19q13.3-13.4, in a cluster with all Siglec-3/CD33-related Siglec genes, suggesting their origin by gene duplications. A homology search of the Drosophila melanogaster and Caenorhabditis elegans genomes suggests that Siglec expression may be limited to animals of deuterostome lineage, coincident with the appearance of the genes of the sialic acid biosynthetic pathway.

DNA archives and our nearest relative: the trichotomy problem revisited

Ever since Thomas H. Huxley correctly identified the chimpanzee and the gorilla as the two closest relatives of the human, the problem of the relationship among the three species ("the trichotomy problem") has remained unresolved. Comparative morphology and other classical methods of biological investigation have failed to answer definitively whether the chimpanzee or the gorilla is the closest relative of the human species. DNA sequences, both mitochondrial and nuclear, too, have provided equivocal solutions, depending on the region of the genome analyzed. Random sorting of ancestral allelic lineages, sequence convergence, and sequence exchanges between alleles or duplicated loci have been identified as likely factors confounding the interpretation of the interrelationships among the three species. In the present study most of these difficulties are overcome by identifying evolutionary causes that might potentially provide misleading information. Altogether, 45 loci consisting of 46, 855 bp are analyzed. About 60% of the loci and approximately the same proportion of phylogenetically informative sites support the human-chimpanzee clade. The remaining 40% of loci and sites support the two alternatives equally. It is demonstrated that, while incompatibility between loci can be explained by random sorting of allelic lineages, incompatibility within loci must be attributed largely to the joint effect of recombination and genetic drift. The trichotomy problem can be properly addressed only within this framework.

Detection of protein coding sequences using a mixture model for local protein amino acid sequence

Locating protein coding regions in genomic DNA is a critical step in accessing the information generated by large scale sequencing projects. Current methods for gene detection depend on statistical measures of content differences between coding and noncoding DNA in addition to the recognition of promoters, splice sites, and other regulatory sites. Here we explore the potential value of recurrent amino acid sequence patterns 3-19 amino acids in length as a content statistic for use in gene finding approaches. A finite mixture model incorporating these patterns can partially discriminate protein sequences which have no (detectable) known homologs from randomized versions of these sequences, and from short (< or = 50 amino acids) non-coding segments extracted from the S. cerevisiea genome. The mixture model derived scores for a collection of human exons were not correlated with the GENSCAN scores, suggesting that the addition of our protein pattern recognition module to current gene recognition programs may improve their performance.

Characterization of a new ATP-binding cassette transporter in Trypanosoma cruzi associated to a L1Tc retrotransposon

We have characterized the tcpgp1-like gene of Trypanosoma cruzi, a new ATP-binding cassette (ABC) transporter. tcpgp1 codes for a 1035 amino acid protein with a considerable homology to LtpgpA of Leishmania. Tcpgp1 lacks the conserved sequences corresponding to the second nucleotide-binding domain of other ABC transporters due to the insertion of the L1Tc non-LTR retrotransposon.