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.

Multiple independent transpositions of mitochondrial DNA control region sequences to the nucleus

Transpositions of mtDNA sequences to the nuclear genome have been documented in a wide variety of individual taxa, but little is known about their taxonomic frequency or patterns of variation. We provide evidence of nuclear sequences homologous to the mtDNA control region in seven species of diving ducks (tribe Aythyini). Phylogenetic analysis places each nuclear sequence as a close relative of the mtDNA haplotypes of the specie(s) in which it occurs, indicating that they derive from six independent transposition events, all occurring within the last ≈1.5 million years. Relative-rate tests and comparison of intraspecific variation in nuclear and mtDNA sequences confirm the expectation of a greatly reduced rate of evolution in the nuclear copies. By representing mtDNA haplotypes from ancestral populations, nuclear insertions may be valuable in some phylogenetic analyses, but they also confound the accurate determination of mtDNA sequences. In particular, our data suggest that the presumably nonfunctional but more slowly evolving nuclear sequences often will not be identifiable by changes incompatible with function and may be preferentially amplified by PCR primers based on mtDNA sequences from related taxa.

The pollen determinant of self-incompatibility in Brassica campestris

Many flowering plants possess self-incompatibility (SI) systems that prevent inbreeding. In Brassica, SI is controlled by a single polymorphic locus, the S locus. Two highly polymorphic S locus genes, SLG (S locus glycoprotein) and SRK (S receptor kinase), have been identified, both of which are expressed predominantly in the stigmatic papillar cell. We have shown recently that SRK is the determinant of the S haplotype specificity of the stigma. SRK is thought to serve as a receptor for a pollen ligand, which presumably is encoded by another polymorphic gene at the S locus. We previously have identified an S locus gene, SP11 (S locus protein 11), of the S9 haplotype of Brassica campestris and proposed that it potentially encodes the pollen ligand. SP11 is a novel member of the PCP (pollen coat protein) family of proteins, some members of which have been shown to interact with SLG. In this work, we identified the SP11 gene from three additional S haplotypes and further characterized the gene. We found that (i) SP11 showed an S haplotype-specific sequence polymorphism; (ii) SP11 was located in the immediate flanking region of the SRK gene of the four S haplotypes examined; (iii) SP11 was expressed in the tapetum of the anther, a site consistent with sporophytic control of Brassica SI; and (iv) recombinant SP11 of the S9 haplotype applied to papillar cells of S9 stigmas, but not of S8 stigmas, elicited SI response, resulting in inhibition of hydration of cross-pollen. All these results taken together strongly suggest that SP11 is the pollen S determinant in SI.

Mutations in the nebulin gene associated with autosomal recessive nemaline myopathy

The congenital nemaline myopathies are rare hereditary muscle disorders characterized by the presence in the muscle fibers of nemaline bodies consisting of proteins derived from the Z disc and thin filament. In a single large Australian family with an autosomal dominant form of nemaline myopathy, the disease is caused by a mutation in the α-tropomyosin gene TPM3. The typical form of nemaline myopathy is inherited as an autosomal recessive trait, the locus of which we previously assigned to chromosome 2q21.2-q22. We show here that mutations in the nebulin gene located within this region are associated with the disease. The nebulin protein is a giant protein found in the thin filaments of striated muscle. A variety of nebulin isoforms are thought to contribute to the molecular diversity of Z discs. We have studied the 3′ end of the 20.8-kb cDNA encoding the Z disc part of the 800-kDa protein and describe six disease-associated mutations in patients from five families of different ethnic origins. In two families with consanguineous parents, the patients were homozygous for point mutations. In one family with nonconsanguineous parents, the affected siblings were compound heterozygotes for two different mutations, and in two further families with one detected mutation each, haplotypes are compatible with compound heterozygosity. Immunofluorescence studies with antibodies specific to the C-terminal region of nebulin indicate that the mutations may cause protein truncation possibly associated with loss of fiber-type diversity, which may be relevant to disease pathogenesis.

Microsatellites provide evidence for Y chromosome diversity among the founders of the New World

Recently, Y chromosome markers have begun to be used to study Native American origins. Available data have been interpreted as indicating that the colonizers of the New World carried a single founder haplotype. However, these early studies have been based on a few, mostly complex polymorphisms of insufficient resolution to determine whether observed diversity stems from admixture or diversity among the colonizers. Because the interpretation of Y chromosomal variation in the New World depends on founding diversity, it is important to develop marker systems with finer resolution. Here we evaluate the hypothesis of a single-founder Y haplotype for Amerinds by using 11 Y-specific markers in five Colombian Amerind populations. Two of these markers (DYS271, DYS287) are reliable indicators of admixture and detected three non-Amerind chromosomes in our sample. Two other markers (DYS199, M19) are single-nucleotide polymorphisms mostly restricted to Native Americans. The relatedness of chromosomes defined by these two markers was evaluated by constructing haplotypes with seven microsatellite loci (DYS388 to 394). The microsatellite backgrounds found on the two haplogroups defined by marker DYS199 demonstrate the existence of at least two Amerind founder haplotypes, one of them (carrying allele DYS199 T) largely restricted to Native Americans. The estimated age and distribution of these haplogroups places them among the founders of the New World.

Polynesian origins: Insights from the Y chromosome

The question surrounding the colonization of Polynesia has remained controversial. Two hypotheses, one postulating Taiwan as the putative homeland and the other asserting a Melanesian origin of the Polynesian people, have received considerable attention. In this work, we present haplotype data based on the distribution of 19 biallelic polymorphisms on the Y chromosome in a sample of 551 male individuals from 36 populations living in Southeast Asia, Taiwan, Micronesia, Melanesia, and Polynesia. Surprisingly, nearly none of the Taiwanese Y haplotypes were found in Micronesia and Polynesia. Likewise, a Melanesian-specific haplotype was not found among the Polynesians. However, all of the Polynesian, Micronesian, and Taiwanese haplotypes are present in the extant Southeast Asian populations. Evidently, the Y-chromosome data do not lend support to either of the prevailing hypotheses. Rather, we postulate that Southeast Asia provided a genetic source for two independent migrations, one toward Taiwan and the other toward Polynesia through island Southeast Asia.

Global survey of genetic variation in CCR5, RANTES, and MIP-1α: Impact on the epidemiology of the HIV-1 pandemic

Expression of CC chemokine receptor 5 (CCR5), the major coreceptor for HIV-1 cell entry, and its ligands (e.g., RANTES and MIP-1α) is widely regarded as central to the pathogenesis of HIV-1 infection. By surveying nearly 3,000 HIV+ and HIV− individuals from worldwide populations for polymorphisms in the genes encoding RANTES, MIP-1α, and CCR5, we show that the evolutionary histories of human populations have had a significant impact on the distribution of variation in these genes, and that this may be responsible, in part, for the heterogeneous nature of the epidemiology of the HIV-1 pandemic. The varied distribution of RANTES haplotypes (AC, GC, and AG) associated with population-specific HIV-1 transmission- and disease-modifying effects is a striking example. Homozygosity for the AC haplotype was associated with an increased risk of acquiring HIV-1 as well as accelerated disease progression in European Americans, but not in African Americans. Yet, the prevalence of the ancestral AC haplotype is high in individuals of African origin, but substantially lower in non-Africans. In a Japanese cohort, AG-containing RANTES haplotype pairs were associated with a delay in disease progression; however, we now show that their contribution to HIV-1 pathogenesis and epidemiology in other parts of the world is negligible because the AG haplotype is infrequent in non-Far East Asians. Thus, the varied distribution of RANTES, MIP-1α, and CCR5 haplotype pairs and their population-specific phenotypic effects on HIV-1 susceptibility and disease progression results in a complex pattern of biological determinants of HIV-1 epidemiology. These findings have important implications for the design, assessment, and implementation of effective HIV-1 intervention and prevention strategies.

Molecular and genealogical evidence for a founder effect in Fanconi anemia families of the Afrikaner population of South Africa

Fanconi anemia (FA) is a rare, genetically heterogeneous autosomal recessive disorder associated with progressive aplastic anemia, congenital abnormalities, and cancer. FA has a very high incidence in the Afrikaner population of South Africa, possibly due to a founder effect. Previously we observed allelic association between polymorphic markers flanking the FA group A gene (FANCA) and disease chromosomes in Afrikaners. We genotyped 26 FA families with microsatellite and single nucleotide polymorphic markers and detected five FANCA haplotypes. Mutation scanning of the FANCA gene revealed association of these haplotypes with four different mutations. The most common was an intragenic deletion of exons 12–31, accounting for 60% of FA chromosomes in 46 unrelated Afrikaner FA patients, while two other mutations accounted for an additional 20%. Screening for these mutations in the European populations ancestral to the Afrikaners detected one patient from the Western Ruhr region of Germany who was heterozygous for the major deletion. The mutation was associated with the same unique FANCA haplotype as in Afrikaner patients. Genealogical investigation of 12 Afrikaner families with FA revealed that all were descended from a French Huguenot couple who arrived at the Cape on June 5, 1688, whereas mutation analysis showed that the carriers of the major mutation were descendants of this same couple. The molecular and genealogical evidence is consistent with transmission of the major mutation to Western Germany and the Cape near the end of the 17th century, confirming the existence of a founder effect for FA in South Africa.

Genes, peoples, and languages

The genetic history of a group of populations is usually analyzed by reconstructing a tree of their origins. Reliability of the reconstruction depends on the validity of the hypothesis that genetic differentiation of the populations is mostly due to population fissions followed by independent evolution. If necessary, adjustment for major population admixtures can be made. Dating the fissions requires comparisons with paleoanthropological and paleontological dates, which are few and uncertain. A method of absolute genetic dating recently introduced uses mutation rates as molecular clocks; it was applied to human evolution using microsatellites, which have a sufficiently high mutation rate. Results are comparable with those of other methods and agree with a recent expansion of modern humans from Africa. An alternative method of analysis, useful when there is adequate geographic coverage of regions, is the geographic study of frequencies of alleles or haplotypes. As in the case of trees, it is necessary to summarize data from many loci for conclusions to be acceptable. Results must be independent from the loci used. Multivariate analyses like principal components or multidimensional scaling reveal a number of hidden patterns and evaluate their relative importance. Most patterns found in the analysis of human living populations are likely to be consequences of demographic expansions, determined by technological developments affecting food availability, transportation, or military power. During such expansions, both genes and languages are spread to potentially vast areas. In principle, this tends to create a correlation between the respective evolutionary trees. The correlation is usually positive and often remarkably high. It can be decreased or hidden by phenomena of language replacement and also of gene replacement, usually partial, due to gene flow.

Gene genealogies and population variation in plants

Early in the development of plant evolutionary biology, genetic drift, fluctuations in population size, and isolation were identified as critical processes that affect the course of evolution in plant species. Attempts to assess these processes in natural populations became possible only with the development of neutral genetic markers in the 1960s. More recently, the application of historically ordered neutral molecular variation (within the conceptual framework of coalescent theory) has allowed a reevaluation of these microevolutionary processes. Gene genealogies trace the evolutionary relationships among haplotypes (alleles) with populations. Processes such as selection, fluctuation in population size, and population substructuring affect the geographical and genealogical relationships among these alleles. Therefore, examination of these genealogical data can provide insights into the evolutionary history of a species. For example, studies of Arabidopsis thaliana have suggested that this species underwent rapid expansion, with populations showing little genetic differentiation. The new discipline of phylogeography examines the distribution of allele genealogies in an explicit geographical context. Phylogeographic studies of plants have documented the recolonization of European tree species from refugia subsequent to Pleistocene glaciation, and such studies have been instructive in understanding the origin and domestication of the crop cassava. Currently, several technical limitations hinder the widespread application of a genealogical approach to plant evolutionary studies. However, as these technical issues are solved, a genealogical approach holds great promise for understanding these previously elusive processes in plant evolution.