Distribution of haplotypes from a chromosome 21 region distinguishes multiple prehistoric human migrations

Despite mounting genetic evidence implicating a recent origin of modern humans, the elucidation of early migratory gene-flow episodes remains incomplete. Geographic distribution of haplotypes may show traces of ancestral migrations. However, such evolutionary signatures can be erased easily by recombination and mutational perturbations. A 565-bp chromosome 21 region near the MX1 gene, which contains nine sites frequently polymorphic in human populations, has been found. It is unaffected by recombination and recurrent mutation and thus reflects only migratory history, genetic drift, and possibly selection. Geographic distribution of contemporary haplotypes implies distinctive prehistoric human migrations: one to Oceania, one to Asia and subsequently to America, and a third one predominantly to Europe. The findings with chromosome 21 are confirmed by independent evidence from a Y chromosome phylogeny. Loci of this type will help to decipher the evolutionary history of modern humans.

Complex promoter and coding region β2-adrenergic receptor haplotypes alter receptor expression and predict in vivo responsiveness

The human β2-adrenergic receptor gene has multiple single-nucleotide polymorphisms (SNPs), but the relevance of chromosomally phased SNPs (haplotypes) is not known. The phylogeny and the in vitro and in vivo consequences of variations in the 5′ upstream and ORF were delineated in a multiethnic reference population and an asthmatic cohort. Thirteen SNPs were found organized into 12 haplotypes out of the theoretically possible 8,192 combinations. Deep divergence in the distribution of some haplotypes was noted in Caucasian, African-American, Asian, and Hispanic-Latino ethnic groups with >20-fold differences among the frequencies of the four major haplotypes. The relevance of the five most common β2-adrenergic receptor haplotype pairs was determined in vivo by assessing the bronchodilator response to β agonist in asthmatics. Mean responses by haplotype pair varied by >2-fold, and response was significantly related to the haplotype pair (P = 0.007) but not to individual SNPs. Expression vectors representing two of the haplotypes differing at eight of the SNP loci and associated with divergent in vivo responsiveness to agonist were used to transfect HEK293 cells. β2-adrenergic receptor mRNA levels and receptor density in cells transfected with the haplotype associated with the greater physiologic response were ≈50% greater than those transfected with the lower response haplotype. The results indicate that the unique interactions of multiple SNPs within a haplotype ultimately can affect biologic and therapeutic phenotype and that individual SNPs may have poor predictive power as pharmacogenetic loci.

Determination of gene organization in individual haplotypes by analyzing single DNA fragments from single spermatozoa

To determine human Ig heavy chain variable region (VH) gene segment organization on individual homologous chromosomes, an efficient approach has been developed. Single spermatozoa were used as subjects for the study. Upon sperm lysis, VH regions in each sperm were randomly sheared into fragments by the random Brownian force. The fragments were separated from each other by aliquoting the lysate into a certain number of tubes. The gene segments in the VH1 and VH4 families in each tube were identified by denaturing gradient gel electrophoresis after PCR amplification. The polymorphic VH sequences were used to determine the parental origins of the analyzed sperm. VH segment organization in the parental haplotypes was determined by aligning the overlapping fragments from the spermatozoa with the corresponding haplotypes. Based on this comparison between the resulting haplotype maps and the composite map reported previously, the VH region on chromosome 14 could be subdivided into four portions. The numbers and compositions of the VH gene segments differ considerably among the maps in two portions, but are highly conserved in the other two. The data also indicate that the VH region on chromosome 15 may contain a large duplicated block with copy number varying among haplotypes. The approach used in the present study may be used to construct high-resolution haplotype maps without molecular cloning.

Xce haplotypes show modified methylation in a region of the active X chromosome lying 3' to Xist.

During early mammalian embryogenesis, one of the two X chromosomes in somatic cells of the female becomes inactivated through a process that is thought to depend on a unique initiator region, the X-chromosome inactivation center (Xic). The recently characterized Xist sequence (X-inactive-specific transcript) is thought to be a possible candidate for Xic. In mice a further genetic element, the X chromosome-controlling element (Xce), is also known to influence the choice of which of the two X chromosomes is inactivated. We report that a region of the mouse X chromosome lying 15 kb distal to Xist contains several sites that show hypermethylation specifically associated with the active X chromosome. Analysis of this region in various Xce strains has revealed a correlation between the strength of the Xce allele carried and the methylation status of this region. We propose that such a region could be involved in the initial stages of the inactivation process and in particular in the choice of which of the two X chromosomes present in a female cell will be inactivated.

Cladistic association analysis of Y chromosome effects on alcohol dependence and related personality traits

Association between Y chromosome haplotype variation and alcohol dependence and related personality traits was investigated in a large sample of psychiatrically diagnosed Finnish males. Haplotypes were constructed for 359 individuals using alleles at eight loci (seven microsatellite loci and a nucleotide substitution in the DYZ3 alphoid satellite locus). A cladogram linking the 102 observed haplotype configurations was constructed by using parsimony with a single-step mutation model. Then, a series of contingency tables nested according to the cladogram hierarchy were used to test for association between Y haplotype and alcohol dependence. Finally, using only alcohol-dependent subjects, we tested for association between Y haplotype and personality variables postulated to define subtypes of alcoholism—antisocial personality disorder, novelty seeking, harm avoidance, and reward dependence. Significant association with alcohol dependence was observed at three Y haplotype clades, with significance levels of P = 0.002, P = 0.020, and P = 0.010. Within alcohol-dependent subjects, no relationship was revealed between Y haplotype and antisocial personality disorder, novelty seeking, harm avoidance, or reward dependence. These results demonstrate, by using a fully objective association design, that differences among Y chromosomes contribute to variation in vulnerability to alcohol dependence. However, they do not demonstrate an association between Y haplotype and the personality variables thought to underlie the subtypes of alcoholism.

Differential expression of major histocompatibility complex class II genes on murine macrophages associated with T cell cytokine profile and protective/suppressive effects

Protective/suppressive major histocompatibility complex (MHC) class II alleles have been identified in humans and mice where they exert a disease-protective and immunosuppressive effect. Various modes of action have been proposed, among them differential expression of MHC class II genes in different types of antigen-presenting cells impacting on the T helper type 1 (Th1)–Th2 balance. To test this possibility, the expression of H-2 molecules from the four haplotypes H-2b, H-2d, H-2k, and H-2q was determined on bone marrow-derived macrophages (BMDMs) and splenic B cells. The I-Ab and I-Ek molecules, both well characterized as protective/suppressive, are expressed at a high level on almost all CD11b+ BMDMs for 5–8 days, after which expression slowly declines. In contrast, I-Ad, I-Ak, and I-Aq expression is lower, peaks over a shorter period, and declines more rapidly. No differential expression could be detected on B cells. In addition, the differential MHC class II expression found on macrophages skews the cytokine response of T cells as shown by an in vitro restimulation assay with BMDMs as antigen-presenting cells. The results indicate that macrophages of the protective/suppressive haplotypes express MHC class II molecules at a high level and exert Th1 bias, whereas low-level expression favors a Th2 response. We suggest that the extent of expression of the class II gene gates the back signal from T cells and in this way controls the activity of macrophages. This effect mediated by polymorphic nonexon segments of MHC class II genes may play a role in determining disease susceptibility in humans and mice.

Philopatry of male marine turtles inferred from mitochondrial DNA markers

Recent studies of mitochondrial DNA (mtDNA) variation among marine turtle populations are consistent with the hypothesis that females return to beaches in their natal region to nest as adults. In contrast, less is known about breeding migrations of male marine turtles and whether they too are philopatric to natal regions. Studies of geographic structuring of restriction fragment and microsatellite polymorphisms at anonymous nuclear loci in green turtle (Chelonia mydas) populations indicate that nuclear gene flow is higher than estimates from mtDNA analyses. Regional populations from the northern and southern Great Barrier Reef were distinct for mtDNA but indistinguishable at nuclear loci, whereas the Gulf of Carpentaria (northern Australia) population was distinct for both types of marker. To assess whether this result was due to reduced philopatry of males across the Great Barrier Reef, we determined the mtDNA haplotypes of breeding males at courtship areas for comparison with breeding females from the same three locations. We used a PCR-restriction fragment length polymorphism approach to determine control region haplotypes and designed mismatch primers for the identification of specific haplotypes. The mtDNA haplotype frequencies were not significantly different between males and females at any of the three areas and estimates of Fst among the regions were similar for males and females (Fst = 0.78 and 0.73, respectively). We conclude that breeding males, like females, are philopatric to courtship areas within their natal region. Nuclear gene flow between populations is most likely occurring through matings during migrations of both males and females through nonnatal courtship areas.

Phylogenetic assessment of length variation at a microsatellite locus

Sixty-six haplotypes at a locus containing a simple dinucleotide (CA)n microsatellite repeat were isolated by PCR–single-strand conformational polymorphism from populations of the horseshoe crab Limulus polyphemus. These haplotypes were sequenced to assess nucleotide variation directly. Thirty-four distinct sequences (alleles) were identified in a region 570 bp long that included the microsatellite motif. In the repeat region itself, CA-number varied in integer values from 5 to 11 across alleles, except that a (CA)8 class was not observed. Differences among alleles were due also to polymorphisms at 22 sites in regions immediately flanking the microsatellite repeats. Nucleotide substitutions in these regions were used to estimate phylogenetic relationships among alleles, and the gene phylogeny was used to trace the evolution of length variation and CA repeat numbers. A low correlation between size variation and genealogical relationships among alleles suggests that absolute fragment size (as normally scored in microsatellite assays) is an unreliable indicator of historical affinities among alleles. This finding on the molecular fine structure of microsatellite variation suggests the need for caution in the use of repeat counts at microsatellite loci as secure indicators of allelic relationships.

Evolution of the Friedreich’s ataxia trinucleotide repeat expansion: Founder effect and premutations

Friedreich’s ataxia, the most frequent inherited ataxia, is caused, in the vast majority of cases, by large GAA repeat expansions in the first intron of the frataxin gene. The normal sequence corresponds to a moderately polymorphic trinucleotide repeat with bimodal size distribution. Small normal alleles have approximately eight to nine repeats whereas a more heterogeneous mode of large normal alleles ranges from 16 to 34 GAA. The latter class accounts for ≈17% of normal alleles. To identify the origin of the expansion mutation, we analyzed linkage disequilibrium between expansion mutations or normal alleles and a haplotype of five polymorphic markers within or close to the frataxin gene; 51% of the expansions were associated with a single haplotype, and the other expansions were associated with haplotypes that could be related to the major one by mutation at a polymorphic marker or by ancient recombination. Of interest, the major haplotype associated with expansion is also the major haplotype associated with the larger alleles in the normal size range and was almost never found associated with the smaller normal alleles. The results indicate that most if not all large normal alleles derive from a single founder chromosome and that they represent a reservoir for larger expansion events, possibly through “premutation” intermediates. Indeed, we found two such alleles (42 and 60 GAA) that underwent cataclysmic expansion to pathological range in a single generation. This stepwise evolution to large trinucleotide expansions already was suggested for myotonic dystrophy and fragile X syndrome and may relate to a common mutational mechanism, despite sequence motif differences.

Plasmodium falciparum antigenic diversity: Evidence of clonal population structure

Plasmodium falciparum, the agent of malignant malaria, is one of mankind’s most severe scourges. Efforts to develop preventive vaccines or remedial drugs are handicapped by the parasite’s rapid evolution of drug resistance and protective antigens. We examine 25 DNA sequences of the gene coding for the highly polymorphic antigenic circumsporozoite protein. We observe total absence of silent nucleotide variation in the two nonrepeated regions of the gene. We propose that this absence reflects a recent origin (within several thousand years) of the world populations of P. falciparum from a single individual; the amino acid polymorphisms observed in these nonrepeat regions would result from strong natural selection. Analysis of these polymorphisms indicates that: (i) the incidence of recombination events does not increase with nucleotide distance; (ii) the strength of linkage disequilibrium between nucleotides is also independent of distance; and (iii) haplotypes in the two nonrepeat regions are correlated with one another, but not with the central repeat region they span. We propose two hypotheses: (i) variation in the highly polymorphic central repeat region arises by mitotic intragenic recombination, and (ii) the population structure of P. falciparum is clonal—a state of affairs that persists in spite of the necessary stage of physiological sexuality that the parasite must sustain in the mosquito vector to complete its life cycle.

Sequence variation in the Fanconi anemia gene FAA

Fanconi anemia (FA) is a genetically heterogeneous autosomal recessive syndrome associated with chromosomal instability, hypersensitivity to DNA crosslinking agents, and predisposition to malignancy. The gene for FA complementation group A (FAA) recently has been cloned. The cDNA is predicted to encode a polypeptide of 1,455 amino acids, with no homologies to any known protein that might suggest a function for FAA. We have used single-strand conformational polymorphism analysis to screen genomic DNA from a panel of 97 racially and ethnically diverse FA patients from the International Fanconi Anemia Registry for mutations in the FAA gene. A total of 85 variant bands were detected. Forty-five of the variants are probably benign polymorphisms, of which nine are common and can be used for various applications, including mapping studies for other genes in this region of chromosome 16q. Amplification refractory mutation system assays were developed to simplify their detection. Forty variants are likely to be pathogenic mutations. Seventeen of these are microdeletions/microinsertions associated with short direct repeats or homonucleotide tracts, a type of mutation thought to be generated by a mechanism of slipped-strand mispairing during DNA replication. A screening of 350 FA probands from the International Fanconi Anemia Registry for two of these deletions (1115–1118del and 3788–3790del) revealed that they are carried on about 2% and 5% of the FA alleles, respectively. 3788–3790del appears in a variety of ethnic groups and is found on at least two different haplotypes. We suggest that FAA is hypermutable, and that slipped-strand mispairing, a mutational mechanism recognized as important for the generation of germ-line and somatic mutations in a variety of cancer-related genes, including p53, APC, RB1, WT1, and BRCA1, may be a major mechanism for FAA mutagenesis.

Genetic epidemiology of single-nucleotide polymorphisms

On the causal hypothesis, most genetic determinants of disease are single-nucleotide polymorphisms (SNPs) that are likely to be selected as markers for positional cloning. On the proximity hypothesis, most disease determinants will not be included among markers but may be detected through linkage disequilibrium with other SNPs. In that event, allelic association among SNPs is an essential factor in positional cloning. Recent simulation based on monotonic population expansion suggests that useful association does not usually extend beyond 3 kb. This is contradicted by significant disequilibrium at much greater distances, with corresponding reduction in the number of SNPs required for a cost-effective genome scan. A plausible explanation is that cyclical expansions follow population bottlenecks that establish new disequilibria. Data on more than 1,000 locus pairs indicate that most disequilibria trace to the Neolithic, with no apparent difference between haplotypes that are random or selected through a major disease gene. Short duration may be characteristic of alleles contributing to disease susceptibility and haplotypes characteristic of particular ethnic groups. Alleles that are highly polymorphic in all ethnic groups may be older, neutral, or advantageous, in weak disequilibrium with nearby markers, and therefore less useful for positional cloning of disease genes. Significant disequilibrium at large distance makes the number of suitably chosen SNPs required for genome screening as small as 30,000, or 1 per 100 kb, with greater density (including less common SNPs) reserved for candidate regions.

Role of the major histocompatibility complex class II Ea gene in lupus susceptibility in mice

The gene(s) encoded within major histocompatibility complex (MHC) act as one of the major genetic elements contributing to the susceptibility of murine systemic lupus erythematosus (SLE). We have recently demonstrated that lupus susceptibility is more closely linked to the I-E− H-2b haplotype than to the I-E+ H-2d haplotype in lupus-prone BXSB and (NZB × BXSB)F1 hybrid mice. To investigate whether the reduced susceptibility to SLE in H-2d mice is related to the expression of the MHC class II Ea gene (absent in H-2b mice), we determined the possible role of the Ea gene as a lupus protective gene in mice. Our results showed that (i) the development of SLE was almost completely prevented in BXSB (H-2b) mice expressing two copies of the Ead transgene at the homozygous level as well as in BXSB H-2k (I-E+) congenic mice as for H-2d BXSB mice, and (ii) the expression of two functional Ea (transgenic and endogenous) genes in either H-2d/b (NZB × BXSB)F1 or H-2k/b (MRL × BXSB)F1 mice provided protection from SLE at levels comparable to those conferred by the H-2d/d or H-2k/k haplotype. In addition, the level of the Ea gene-mediated protection appeared to be dependent on the genetic susceptibility to SLE in individual lupus-prone mice. Our results indicate that the reduced susceptibility associated with the I-E+ H-2d and H-2k haplotypes (versus the I-E− H-2b haplotype) is largely, if not all, contributed by the apparent autoimmune suppressive effect of the Ea gene, independently of the expression of the I-A or other MHC-linked genes.

Protective anti-tumor immunity induced by vaccination with recombinant adenoviruses encoding multiple tumor-associated cytotoxic T lymphocyte epitopes in a string-of-beads fashion

Vaccines harboring genes that encode functional oncoproteins are intrinsically hazardous, as their application may lead to introduction of these genes into normal cells and thereby to tumorigenesis. On the other hand, oncoproteins are especially attractive targets for immunotherapy of cancer, as their expression is generally required for tumor growth, making the arisal of tumor variants lacking these antigens unlikely. Using murine tumor models, we investigated the efficacy of polyepitope recombinant adenovirus (rAd) vaccines, which encode only the immunogenic T cell epitopes derived from several oncogenes, for the induction of protective anti-tumor immunity. We chose to employ rAd, as these are safe vectors that do not induce the side effects associated with, for example, vaccinia virus vaccines. A single polyepitope rAd was shown to give rise to presentation of both H-2 and human leukocyte antigen-restricted cytotoxic T lymphocyte (CTL) epitopes. Moreover, vaccination with a rAd encoding H-2-restricted CTL epitopes, derived from human adenovirus type 5 early region 1 and human papilloma virus type 16-induced tumors, elicited strong tumor-reactive CTL and protected the vaccinated animals against an otherwise lethal challenge with either of these tumors. The protection induced was superior compared with that obtained by vaccination with irradiated tumor cells. Thus, vaccination with polyepitope rAd is a powerful approach for the induction of protective anti-tumor immunity that allows simultaneous immunization against multiple tumor-associated T cell epitopes, restricted by various major histocompatibility complex haplotypes.