An anomalous haplotype distribution of the arrestin domain-containing 4 gene (ARRDC4) haplotypes in Caucasians

Little was known about the sequence variability of the human Arrestin domain-containing 4 gene (ARRDC4). We sequenced its DNA from exon 2 to exon 8 in a sample of 92 Russians. Seven variants were identified; one of them has not been described yet. It causes an amino acid change from Thr to Met. Identified variants were genotyped in the complete sample of 253 unrelated men and women to analyze haplotype distribution. Fifteen haplotypes were inferred. Nine haplotypes had estimated frequencies > 1%. Ninety-five percent of all haplotypes were determined by five haplotype-tagging single nucleotide polymorphisms. Haplotypes form two clades. The two most common haplotypes cover 76% of all haplotypes. The certainty of the haplotype reconstruction does not depend on the haplotype-inferring algorithms, but is a result of the anomalous haplotype distribution of ARRDC4, which makes this gene a suitable candidate gene for haplotype association studies. Interestingly, there is a great evolutionary distance between the two most common haplotypes, which could suggest a more complicated coalescent process with either past gene flow, selections, or bottlenecks.     

Haplotypes at ATM identify coding-sequence variation and indicate a region of extensive linkage disequilibrium

Genetic variation in the human population may lead to functional variants of genes that contribute to risk for common chronic diseases such as cancer. In an effort to detect such possible predisposing variants, we constructed haplotypes for a candidate gene and tested their efficacy in association studies. We developed haplotypes consisting of 14 biallelic neutral-sequence variants that span 142 kb of the ATM locus. ATM is the gene responsible for the autosomal recessive disease ataxia-telangiectasia (AT). These ATM noncoding single-nucleotide polymorphisms (SNPs) were genotyped in nine CEPH families (89 individuals) and in 260 DNA samples from four different ethnic origins. Analysis of these data with an expectation-maximization algorithm revealed 22 haplotypes at this locus, with three major haplotypes having frequencies > or = .10. Tests for recombination and linkage disequilibrium (LD) show reduced recombination and extensive LD at the ATM locus, in all four ethnic groups studied. The most striking example was found in the study population of European ancestry, in which no evidence for recombination could be discerned. The potential of ATM haplotypes for detection of genetic variants through association studies was tested by analysis of 84 individuals carrying one of three ATM coding SNPs. Each coding SNP was detected by association with an ATM haplotype. We demonstrate that association studies with haplotypes for candidate genes have significant potential for the detection of genetic backgrounds that contribute to disease.     

Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium

Common genetic polymorphisms may explain a portion of the heritable risk for common diseases. Within candidate genes, the number of common polymorphisms is finite, but direct assay of all existing common polymorphism is inefficient, because genotypes at many of these sites are strongly correlated. Thus, it is not necessary to assay all common variants if the patterns of allelic association between common variants can be described. We have developed an algorithm to select the maximally informative set of common single-nucleotide polymorphisms (tagSNPs) to assay in candidate-gene association studies, such that all known common polymorphisms either are directly assayed or exceed a threshold level of association with a tagSNP. The algorithm is based on the r(2) linkage disequilibrium (LD) statistic, because r(2) is directly related to statistical power to detect disease associations with unassayed sites. We show that, at a relatively stringent r(2) threshold (r2>0.8), the LD-selected tagSNPs resolve >80% of all haplotypes across a set of 100 candidate genes, regardless of recombination, and tag specific haplotypes and clades of related haplotypes in nonrecombinant regions. Thus, if the patterns of common variation are described for a candidate gene, analysis of the tagSNP set can comprehensively interrogate for main effects from common functional variation. We demonstrate that, although common variation tends to be shared between populations, tagSNPs should be selected separately for populations with different ancestries.     

Length Variation in 18S rRNA Expansion Segment 43/e4 of <Emphasis Type="Italic">Daphnia obtusa</Emphasis>: Ancient or Recurring Polymorphism?

Expansion segments in ribosomal DNA (rDNA) can show length variation at the level of the individual, yet our understanding of the evolutionary forces shaping this variation is incomplete. Previous studies of expansion segment 43/e4 of the 18S rRNA gene in Daphnia obtusa have examined this variation in six individuals; however, it is not known if the variation documented at this locus is representative of variation across the species’ geographic range. Furthermore, it is unclear whether length variants found in multiple individuals share common ancestry, or were generated de novo through recombination. We quantified expansion segment length variant frequencies in 134 individual D. obtusa from 33 populations at 15 sites across the species range in the US, and used a phylogeographic approach to determine whether recombination continues to add to the standing crop of variation at this locus. We identified seven length variants across the sampling range, which spans almost 3000 km. Based on the phylogeographic distribution of length variants in the expansion segment, we conclude that they are shared ancient polymorphisms that have persisted despite the operation of molecular mechanisms that cause the concerted evolution of multigene families such as rDNA. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Contrasting patterns of evolution between allelic groups at a single locus in <Emphasis Type="Italic">Arabidopsis</Emphasis>

Heterogeneities in evolutionary pattern among different loci are commonly observed. To see whether the heterogeneity can also be observed among allelic groups in a single locus, we investigated the coding sequence and the flanking regions of Rpp13, a disease resistance gene in up to 60 accession lines from worldwide populations in Arabidopsis thaliana. An extraordinarily high level of polymorphism (π=0.098) and four distinct clades were found in the leucine-rich repeat (LRR) region in this gene. No obvious geographic relationship with the clades was observed, and such clades were not observed in the other regions in and around this gene. The average genetic diversity among the clades ranged from 10 to 14.6% in the LRR. The levels of polymorphism within each clade varied largely, and significant heterogeneity in evolutionary rates among clades was detected. A statistically significant departure from neutrality was also detected by Fu & Li’s tests. These results suggest that both directional and diversifying selection are working on this locus, and that natural selection can cause heterogeneity in evolutionary rate, even among allele groups in a locus. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Corticotropin releasing hormone (CRH) gene variation: comprehensive resequencing for variant and molecular haplotype discovery in monosomic hybrid cell lines.

Candidate gene association studies have met with mixed success due to many reasons including incomplete surveys of genetic variation and differences in patterns of genetic variation among study populations. We present the results of comprehensive variant discovery for the corticotropin releasing hormone gene (CRH on chromosome 8) encoding a neuropeptide that is central to many physiologic pathways. Mouse-human hybrid cell lines were constructed that are monosomic for human chromosome 8 for resequencing of separated CRH alleles to identify variants and directly determine their chromosomal phase for three major ethnic groups including African Americans (AA), Mexican Americans (MA) and European Americans (EA). We also resequenced diploid individuals to evaluate single nucleotide polymorphism (SNP) discovery in the limited numbers of monosomic hybrid cell lines. Our results show that CRH variation is very different in AA, yielding larger numbers of variants and haplotypes compared to MA and EA. Analysis of LD structure found three haplotype blocks in AA and two blocks in EA. Comparisons between AA and EA groups yielded extremely high measures of genetic differentiation (Wright's F(ST)>0.6), likely reflecting disruptive selection in CRH evolution. Network analysis showed that AA have retained an ancestral CRH haplotype, while the most common EA haplotype is derived from a single recombination event.     

Molecular phylogenetics of Candida albicans

We analyzed data on multilocus sequence typing (MLST), ABC typing, mating type-like locus (MAT) status, and antifungal susceptibility for a panel of 1,391 Candida albicans isolates. Almost all (96.7%) of the isolates could be assigned by MLST to one of 17 clades. eBURST analysis revealed 53 clonal clusters. Diploid sequence type 69 was the most common MLST strain type and the founder of the largest clonal cluster, and examples were found among isolates from all parts of the world. ABC types and geographical origins showed statistically significant variations among clades by univariate analysis of variance, but anatomical source and antifungal susceptibility data were not significantly associated. A separate analysis limited to European isolates, thereby minimizing geographical effects, showed significant differences in the proportions of isolates from blood, commensal carriage, and superficial infections among the five most populous clades. The proportion of isolates with low antifungal susceptibility was highest for MAT homozygous a/a types and then alpha/alpha types and was lowest for heterozygous a/alpha types. The tree of clades defined by MLST was not congruent with trees generated from the individual gene fragments sequenced, implying a separate evolutionary history for each fragment. Analysis of nucleic acid variation among loci and within loci supported recombination. Computational haplotype analysis showed a high frequency of recombination events, suggesting that isolates had mixed evolutionary histories resembling those of a sexually reproducing species.     

Two Genotypes of Pedicularis chamissonis (Scrophulariaceae) Distributed at Mt. Gassan, Japan: Additional Genetic and Morphological Studies

Pedicularis chamissonis Steven, which is broadly distributed from northern Honshu of Japan to southwest Alaska, is known to comprise a few different cpDNA haplotypes grouped into two major (“Northern” and “Southern”) clades. This paper presents on the basis of chloroplast and nuclear DNA analyses that no gene flow exists between the two clades at Mt Gassan where two haplotypes D and J of the Northern and Southern clades occur. The two haplotypes are further distinguished from each other in corolla and leaf shape and the number of flower whorls on an inflorescence axis, a fact indicating that individuals of the two clades should be treated as separate species. At Mt. Gassan and in other regions in general the haplotypes of the Southern clade currently grow at higher altitudes than those of the Northern clade. Based on the difference in habitat as well as the lack of gene flow between the two clades, evolutionary processes resulting in the current distribution of the two genotypes are also briefly discussed. Received 3 October 2000/ Accepted in revised form 4 March 2001     

Using Haplotype Trees for Phylogeographic and Species Inference in Fish Populations

Genetic variation is now routinely screened at the DNA sequence level in many studies. If the DNA region being screened has not experienced excessive amounts of recombination, it is often possible to reconstruct the evolutionary history of the genetic variation in the form of a haplotype tree. This tree estimates the evolutionary pathway that interconnects all the different haplotypes (sequence variants) observed in the sample. This haplotype tree can be used to define a series of nested branches (clades) that reflects the relative temporal history of the haplotypes and groups of haplotypes. Geographical information can then be overlaid upon this temporal series to test for significant associations between geography and temporal position in the haplotype tree. This allows a reconstruction of how the genetic variation arose and spread in both space and time. Such reconstructions can yield many insights into the joint roles of recurrent events such as gene flow and of historical events such as fragmentation or range expansion. These points are illustrated with studies on the chub, Leuciscus cephalus. There is also a need to extend such nested phylogeographic analyses to a phylo/reticulate geographic analysis that incorporates both assortment and recombination between and within DNA regions. A preliminary phylo/reticulate geographic analysis is presented of the transferrin locus in the brown trout, Salmo trutta, species complex that reveals the importance of hybridization in the recent evolutionary history of this group. This example shows the inadequacy of a strictly phylogenetic approach and illustrates the need to incorporate reticulate evolution. The results of nested clade phylogeographic analysis and the new phylo/reticulate geographic analysis are then used for inferring species status of the marbled trout. The results indicate that an old hybridization event may have played a role in the origin of the marbled trout. Currently the marbled trout is primarily endangered by hybridization with introduced brown trout. These results show both the positive and negative impacts of hybridization upon biodiversity. Such phylo/reticulate geographic studies will challenge both our concepts of species and our conservation management strategies.     

Evolutionary signatures of common human cis-regulatory haplotypes

Variation in gene expression may give rise to a significant fraction of inter-individual phenotypic variation. Studies searching for the underlying genetic controls for such variation have been conducted in model organisms and humans in recent years. In our previous effort of assessing conserved underlying haplotype patterns across ethnic populations, we constructed common haplotypes using SNPs having conserved linkage disequilibrium (LD) across ethnic populations. These common haplotypes cluster into a simple evolutionary structure based on their frequencies, defining only up to three conserved clusters termed 'haplotype frameworks'. One intriguing preliminary finding was that a significant portion of reported variants strongly associated with cis-regulation tags these globally conserved haplotype frameworks. Here we expand the investigation by collecting genes showing stringently determined cis-association between genotypes and expression phenotypes from major studies. We conducted phylogenetic analysis of current major haplotypes along with the corresponding haplotypes derived from chimpanzee reference sequences. Our analysis reveals that, for the vast majority of such cis-regulatory genes, the tagging SNPs showing the strongest association also tag the haplotype lineages directly separated from ancestry, inferred from either chimpanzee reference sequences or the allele frequency-derived haplotype frameworks, suggesting that the differentially expressed phenotypes were evolved relatively early in human history. Such evolutionary signatures provide keys for a more effective identification of globally-conserved candidate regulatory haplotypes across human genes in future epidemiologic and pharmacogenetic studies.     

An analysis of genetic variation across the MBL2 locus in Dutch Caucasians indicates that 3′ haplotypes could modify circulating levels of mannose-binding lectin

Mannose-binding protein (MBL) is a critical component of innate immunity and provides first-line protection against pathogens. Both circulating MBL serum levels and functional activity have been correlated with common genetic variants in the MBL2 gene. Associations between MBL deficiency and severe infections have been reported in immuno-incompetent patients and for autoimmune disorders; however, measured MBL serum levels do not fully correlate with the ‘secretor haplotypes’. Previously, the MBL2 locus was resequenced and determined that a recombination hotspot divides MBL2 into two haplotype blocks. It was sought to investigate whether additional variants, in either block structure could associate with MBL serum levels. Therefore, 31 common variants were analysed across the locus in 212 DNA samples of healthy Caucasian individuals with known MBL serum concentrations. The additional 5′ variants were in strong linkage to the elements of the ‘secretor haplotypes’; functional alleles B, C and D also lie on restricted haplotypes. Four variants in the 3′ block (Ex4-1483T>C, Ex4-1067G>A, Ex4-901G>A and Ex4-710G>A) are components of a distinct haplotype block. The results of this study suggest that additional 5′ variants as well as markers of distinct 3′ haplotype blocks in MBL2 may contribute to circulating protein levels, but further studies are required to confirm these observations. Last, there could be a selective advantage for diversification of the 3′ region of the gene.Electronic Supplementary Material Supplementary material is available for this article at     

Intraspecific sequence variation of chloroplast DNA inPedicularis chamissonis Steven (Scrophulariaceae) and geographic structuring of the Japanese “Alpine” plants

In order to clarify evolutionary patterns and processes of intraspecific diversification ofPedicularis chamissonis Steven, we analyzed intraspecific variation of the nucleotide sequences of non-coding regions of chloroplast DNA: the intergenic spacers betweentrnT (UGU) andtrnL (UAA) 5′exon,trnL (UAA) 3′exon andtrnF (GAA), andatpB andrbcL. In 24 populations ofP. chamissonis, 33 nucleotide substitutions and 12 insertions/deletions were inferred, and their genetic distances ranged from 0.001 to 0.014. Seventeen distinct cpDNA haplotypes could be recognized and each haplotype was found to be geographically structured. Two major clades (the Northern and Southern clades) were revealed in phylogenetic analyses of cpDNA haplotypes. The haplotypes of the Northern clade had a wider distribution area in the populations of Mts. lide of central Honshu in Japan, northward to Unalaska Island in the Aleutians. Relationships among most haplotypes were unresolved polytomies. On the other hand, the haplotypes of the Southern clade occurred from the populations of Mt. Gassan southwards to Mt. Arakawa of central Honshu. Within this clade, three subclades were clearly recognized. From these results, we concluded that the haplotypes of the Northern and Southern clades inP. chamissonis might have traveled down to Japanese Archipelago from the north in not a single glacial period.     

Polymorphism,haplotype composition,and selection in the <Emphasis Type="Italic">Mhc-DRB</Emphasis> of wild baboons

General patterns of organization in the major histocompatibility complex (MHC) have been successfully explained by the model of birth-and-death evolution, but understanding why certain MHC genes are maintained together into specific haplotypes remains challenging. The haplotype configurations of the functionally important class II DR region have been described in few primates and display important interspecific variability with respect to the extent of allelic variation, the number of loci and/or combinations of loci present. Understanding the evolutionary mechanisms driving such variation is conditional upon characterizing haplotypes in new species and identifying the selective pressures acting on haplotypes. This study explores the variability of haplotype configurations in the Mhc-DRB region (exon 2) for the first time in wild non-human primates, chacma baboons (Papio ursinus). Paur-DRB haplotypes were characterized through segregation studies and linkage disequilibrium. 23 Paur-DRB sequences and 15 haplotype configurations were identified in 199 animals. The Paur-DRB exon 2 is shown to be subjected to intense positive selection and frequent recombination. An approach recently developed for human vaccine studies was used to classify Paur-DRB sequences into supertypes, based on the physico-chemical properties of amino acids that are positively selected, thus most probably involved in antigen recognition. Sequences grouped into the same supertype (thus presumably sharing antigen-binding affinities) are non-randomly distributed within haplotypes, leading to an increased individual diversity of supertypes. Our results suggest that selection favoring haplotypes with complementary sets of DRB supertypes shapes functionally tuned haplotypes in this natural baboon population.     

Linkage disequilibrium and the mapping of complex human traits.

The potential value of haplotypes defined by several single nucleotide polymorphisms has attracted recent interest. With sufficient linkage disequilibrium (LD), haplotypes could be used in association studies to map common alleles that might influence the susceptibility to common diseases, as well as for reconstructing the evolution of the genome. It has been proposed that a globally useful resource need only be based on high frequency variants, identified from a few modest samples. Rapid progress has been made in quantifying the pattern of human LD and haplotypes defined by such common variants within and among populations. However, the quality and utility of the proposed LD-based resource could be seriously compromised if important sampling and analytical factors are overlooked in its design. The LD map should be based on adequately justified criteria defined by sound population genetic principles.     

Mitochondrial cytochrome <Emphasis Type="Italic">b</Emphasis> sequences resolve the taxonomy of field mice (<Emphasis Type="Italic">Apodemus</Emphasis>) in the western Balkan refugium

Apodemus sylvaticus stankovici, described from the topographically rough landscape of the western Balkan glacial refugium, was recently proposed as being either a junior synonym of Apodemus flavicollis or a species on its own right. To untangle this taxonomic vagueness, we sequenced complete cytochrome b gene in 28 field mice collected at 12 locations in the mountains of Bosnia and Herzegovina, Montenegro, western Macedonia and northern Greece. Samples yielded 27 new haplotypes which clustered into two distinct groups. One of these clades also included the reference haplotype of A. flavicollis, while another cluster emerged as being identical with the reference sample for A. sylvaticus. As is common in Apodemus, both species retrieved in our analysis were characterized by low levels of intraspecific variation (0.4–0.9%) as opposed to a high level of differentiation between them (8.0–10.0%); therefore, the taxonomic classification of our material was without doubt. We found no evidence regarding the presence of an additional cryptic species in the mountains of the western Balkans. The very similar values of genetic variability in the two species imply their common evolutionary history of a long-term coexistence in the western Balkan refugium.     

Evolutionary history of the European whitefish Coregonus lavaretus (L.) species complex as inferred from mtDNA phylogeography and gill-raker numbers

We compared mitochondrial DNA and gill-raker number variation in populations of the European whitefish Coregonus lavaretus (L.) species complex to illuminate their evolutionary history, and discuss mechanisms behind diversification. Using single-strand conformation polymorphism (SSCP) and sequencing 528 bp of combined parts of the cytochrome oxidase b (cyt b) and NADH dehydrogenase subunit 3 (ND3) mithochondrial DNA (mtDNA) regions, we documented phylogeographic relationships among populations and phylogeny of mtDNA haplotypes. Demographic events behind geographical distribution of haplotypes were inferred using nested clade analysis (NCA) and mismatch distribution. Concordance between operational taxonomical groups, based on gill-raker numbers, and mtDNA patterns was tested. Three major mtDNA clades were resolved in Europe: a North European clade from northwest Russia to Denmark, a Siberian clade from the Arctic Sea to southwest Norway, and a South European clade from Denmark to the European Alps, reflecting occupation in different glacial refugia. Demographic events inferred from NCA were isolation by distance, range expansion, and fragmentation. Mismatch analysis suggested that clades which colonized Fennoscandia and the Alps expanded in population size 24 500-5800 years before present, with minute female effective population sizes, implying small founder populations during colonization. Gill-raker counts did not commensurate with hierarchical mtDNA clades, and poorly with haplotypes, suggesting recent origin of gill-raker variation. Whitefish designations based on gill-raker numbers were not associated with ancient clades. Lack of congruence in morphology and evolutionary lineages implies that the taxonomy of this species complex should be reconsidered.     

Evidence of recombination in putative ancient asexuals

Ancient asexuals have been considered to be a contradiction of the basic tenets of evolutionary theory. Barred from rearranging genetic variation by recombination, their reduced number of gene arrangements is thought to hamper their response to changing environments. For the same reason, it should be difficult for them to avoid the build-up of deleterious mutations. Several groups of taxonomically diverse organisms are thought to be ancient asexuals, although clear evidence for or against the existence of recombination events is scarce. Several methods have recently been developed for predicting recombination events by analyzing aligned sequences of a given region of DNA that all originate from one species. The methods are based on phylogenetic, substitution, and compatibility analyses. Here we present the results of analyses of sequence data from different loci studied in several groups of evolutionarily distant species that are considered to be ancient asexuals, using seven different types of analysis. The groups of organisms were the arbuscular mycorrhizal fungi (Glomales), Darwinula stevensoni (Darwinuloidea crustacean ostracods) and the bdelloid rotifers (Bdelloidea), which are thought to have been asexual for the last 400, 25-100, and 35-40 Myr, respectively. The seven different analytical methods evaluated the evolutionary relationships among haplotypes, and these methods had previously been shown to be reliable for predicting the occurrence of recombination events. Despite the different degree of genetic variation among the different groups of organisms, at least some evidence for recombination was found in all species groups. In particular, predictions of recombination events in the arbuscular mycorrhizal fungi were frequent. Predictions of recombination were also found for sequence data that have previously been used to infer the absence of recombination in bdelloid rotifers. Although our results have to be taken with some caution because they could signal very ancient recombination events or possibly other genetic variation of nonrecombinant origin, they suggest that some cryptic recombination events may exist in these organisms.     

The R1 resistance gene cluster contains three groups of independently evolving, type I R1 homologues and shows substantial structural variation among haplotypes of Solanum demissum

Cultivated and wild potatoes contain a major disease-resistance cluster on the short arm of chromosome V, including the R1 resistance (R) gene against potato late blight. To explore the functional and evolutionary significance of clustering in the generation of novel disease-resistance genes, we constructed three approximately 1 Mb physical maps in the R1 gene region, one for each of the three genomes (haplotypes) of allohexaploid Solanum demissum, the wild potato progenitor of the R1 locus. Totals of 691, 919 and 559 kb were sequenced for each haplotype, and three distinct resistance-gene families were identified, one homologous to the potato R1 gene and two others homologous to either the Prf or the Bs4 R-gene of tomato. The regions with R1 homologues are highly divergent among the three haplotypes, in contrast to the conserved flanking non-resistance gene regions. The R1 locus shows dramatic variation in overall length and R1 homologue number among the three haplotypes. Sequence comparisons of the R1 homologues show that they form three distinct clades in a distance tree. Frequent sequence exchanges were detected among R1 homologues within each clade, but not among those in different clades. These frequent sequence exchanges homogenized the intron sequences of homologues within each clade, but did not homogenize the coding sequences. Our results suggest that the R1 homologues represent three independent groups of fast-evolving type I resistance genes, characterized by chimeric structures resulting from frequent sequence exchanges among group members. Such genes were first identified among clustered RGC2 genes in lettuce, where they were distinguished from slow-evolving type II R-genes. Our findings at the R1 locus in S. demissum may indicate that a common or similar mechanism underlies the previously reported differentiation of type I and type II R-genes and the differentiation of type I R-genes into distinct groups, identified here.