Adaptive evolution of a candidate gene for aging in Drosophila

Examination of the phenotypic effects of specific mutations has been extensively used to identify candidate genes affecting traits of interest. However, such analyses do not reveal anything about the evolutionary forces acting at these loci, or whether standing allelic variation contributes to phenotypic variance in natural populations. The Drosophila gene methuselah (mth) has been proposed as having major effects on organismal stress response and longevity phenotype. Here, we examine patterns of polymorphism and divergence at mth in population level samples of Drosophila melanogaster, D. simulans, and D. yakuba. Mth has experienced an unusually high level of adaptive amino acid divergence concentrated in the intra- and extracellular loop domains of the receptor protein, suggesting the historical action of positive selection on those regions of the molecule that modulate signal transduction. Further analysis of single nucleotide polymorphisms (SNPs) in D. melanogaster provided evidence for contemporary and spatially variable selection at the mth locus. In ten surveyed populations, the most common mth haplotype exhibited a 40% cline in frequency that coincided with population level differences in multiple life-history traits including lifespan. This clinal pattern was not associated with any particular SNP in the coding region, indicating that selection is operating at a closely linked site that may be involved in gene expression. Together, these consistently nonneutral patterns of inter- and intraspecific variation suggest adaptive evolution of a signal transduction pathway that may modulate lifespan in nature.

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.

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.

Y chromosome polymorphism is a strong determinant of male fitness in Drosophila melanogaster

In many species, the Y (or W) chromosome carries relatively few functional genes. This observation motivates the null hypothesis that the Y will be a minor contributor to genetic variation for fitness. Previous data and theory supported the null hypothesis, but evidence presented here shows that the Y of Drosophila melanogaster is a major determinant of a male's total fitness, with standing genetic variation estimated to be 68% of that of an entire X/autosome genomic haplotype. Most Y-linked genes are expressed during spermatogenesis, and correspondingly, we found that the Y influences fitness primarily through its effect on a male's reproductive success (sperm competition and/or mating success) rather than his egg-to-adult viability. But the fitness of a Y highly depended on the genetic makeup of its bearer, reverting from high to low in different genetic backgrounds. This pattern leads to large epistatic (inconsistent among backgrounds) but no additive (consistent among backgrounds) Y-linked genetic variance for fitness. On a microevolutionary scale, the observed large epistatic variation on the Y substantially reduces heritable variation for fitness among males, and on a macroevolutionary scale, the Y produces strong selection for genomic rearrangements that move interacting genes onto the nonrecombining region of the Y.

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.

Transloading of tumor cells with foreign major histocompatibility complex class I peptide ligand: a novel general strategy for the generation of potent cancer vaccines.

The major hurdle to be cleared in active immunotherapy of cancer is the poor immunogenicity of cancer cells. In previous attempts to overcome this problem, whole tumor cells have been used as vaccines, either admixed with adjuvant(s) or genetically engineered to express nonself proteins or immunomodulatory factors before application. We have developed a novel approach to generate an immunogeneic, highly effective vaccine: major histocompatibility complex (MHC) class I-positive cancer cells are administered together with MHC class I-matched peptide ligands of foreign, nonself origin, generated by a procedure we term transloading. Murine tumor lines of the H2-Kd or the H2-Db haplotype, melanoma M-3 and B16-F10, respectively, as well as colon carcinoma CT-26 (H2-Kd), were transloaded with MHC-matched influenza virus-derived peptides and applied as irradiated vaccines. Mice bearing a deposit of live M-3 melanoma cells were efficiently cured by this treatment. In the CT-26 colon carcinoma and the B16-F10 melanoma, high efficacies were obtained against tumor challenge, suggesting the universal applicability of this new type of vaccine. With foreign peptide ligands adapted to the requirements of a desired MHC class I haplotype, this concept may be used for the treatment of human cancers.

Patterns of kinship in groups of free-living sperm whales (Physeter macrocephalus) revealed by multiple molecular genetic analyses.

Mature female sperm whales (Physeter macrocephalus) live in socially cohesive groups of 10-30, which include immature animals of both sexes, and within which there is communal care of the young. We examined kinship in such groups using analyses of microsatellite DNA, mitochondrial DNA sequence, and sex-linked markers on samples of sloughed skin collected noninvasively from animals in three groups off the coast of Ecuador. Social groups were defined through photographic identification of individuals. Each group contained about 26 members, mostly female (79%). Relatedness was greater within groups, as compared to between groups. Particular mitochondrial haplotypes were characteristic of groups, but all groups contained more than one haplotype. The data are generally consistent with each group being comprised of several matrillines from which males disperse at about the age of 6 years. There are indications of paternal relatedness among grouped individuals with different mitochondrial haplotypes, suggesting long-term associations between different matrilines.

Isolation of an immunodominant viral peptide that is endogenously bound to the stress protein GP96/GRP94.

Heat shock protein gp96 primes class I restricted cytotoxic T cells against antigens present in the cells from which it was isolated. Moreover, gp96 derived from certain tumors functions as an effective vaccine, causing complete tumor regressions in in vivo tumor challenge protocols. Because tumor-derived gp96 did not differ from gp96 isolated from normal tissues, a role for gp96 as a peptide carrier has been proposed. To test this hypothesis, we analyzed whether such an association of antigenic peptides with gp96 occurs in a well-defined viral model system. Here we present the full characterization of an antigenic peptide that endogenously associates with the stress protein gp96 in cells infected with vesicular stomatitis virus (VSV). This peptide is identical to the immunodominant peptide of VSV, which is also naturally presented by H-2Kb major histocompatibility complex class I molecules. This peptide associates with gp96 in VSV-infected cells regardless of the major histocompatibility com- plex haplotype of the cell. Our observations provide a biochemical basis for the vaccine function of gp96.

A pre-Columbian Y chromosome-specific transition and its implications for human evolutionary history.

A polymorphic C-->T transition located on the human Y chromosome was found by the systematic comparative sequencing of Y-specific sequence-tagged sites by denaturing high-performance liquid chromatography. The results of genotyping representative global indigenous populations indicate that the locus is polymorphic exclusively within the Western Hemisphere. The pre-Columbian T allele occurs at > 90% frequency within the native South and Central American populations examined, while its occurrence in North America is approximately 50%. Concomitant genotyping at the polymorphic tetranucleotide microsatellite DYS19 locus revealed that the C-->T mutation displayed significant linkage disequilibrium with the 186-bp allele. The data suggest a single origin of linguistically diverse native Americans with subsequent haplotype differentiation within radiating indigenous populations as well as post-Columbian European and African gene flow. The mutation may have originated either in North America at a very early time during the expansion or before it, in the ancestral population(s) from which all Americans may have originated. The analysis of linkage of the DYS199 and the DYS19 tetranucleotide loci suggests that the C-->T mutation may have occurred around 30,000 years ago. We estimate the nucleotide diversity over 4.2 kb of the nonrecombining portion of the Y chromosome to be 0.00014. compared to autosomes, the majority of variation is due to the smaller effective population size of the Y chromosome rather than selective sweeps. There begins to emerge a pattern of pronounced geographical localization of Y-specific nucleotide substitution polymorphisms.

Periodic variation in side-chain polarities of T-cell antigenic peptides correlates with their structure and activity.

We present an analysis that synthesizes information on the sequence, structure, and motifs of antigenic peptides, which previously appeared to be in conflict. Fourier analysis of T-cell antigenic peptides indicates a periodic variation in amino acid polarities of 3-3.6 residues per period, suggesting an amphipathic alpha-helical structure. However, the diffraction patterns of major histocompatibility complex (MHC) molecules indicate that their ligands are in an extended non-alpha-helical conformation. We present two mutually consistent structural explanations for the source of the alpha-helical periodicity, based on an observation that the side chains of MHC-bound peptides generally partition with hydrophobic (hydrophilic) side chains pointing into (out of) the cleft. First, an analysis of haplotype-dependent peptide motifs indicates that the locations of their defining residues tend to force a period 3-4 variation in hydrophobicity along the peptide sequence, in a manner consistent with the spacing of pockets in the MHC. Second, recent crystallographic determination of the structure of a peptide bound to a class II MHC molecule reveals an extended but regularly twisted peptide with a rotation angle of about 130 degrees. We show that similar structures with rotation angles of 100-130 degrees are energetically acceptable and also span the length of the MHC cleft. These results provide a sound physical chemical and structural basis for the existence of a haplotype-independent antigenic motif which can be particularly important in limiting the search time for antigenic peptides.

Staphylococcal enterotoxins bind H-2Db molecules on macrophages.

We screened a panel of monoclonal antibodies against selected macrophage cell surface molecules for their ability to inhibit enterotoxin binding to major histocompatibility complex class II-negative C2D (H-2b) macrophages. Two monoclonal antibodies, HB36 and TIB126, that are specific for the alpha 2 domain of major histocompatibility complex class I, blocked staphylococcal enterotoxins A and B (SEA and SEB, respectively) binding to C2D macrophages in a specific and concentration-dependent manner. Inhibitory activities were haplotype-specific in that SEA and SEB binding to H-2k or H-2d macrophages was not inhibited by either monoclonal antibody. HB36, but not TIB126, inhibited enterotoxin-induced secretion of cytokines by H-2b macrophages. Lastly, passive protection of D-galactosamine-sensitized C2D mice by injection with HB36 antibody prevented SEB-induced death. Therefore, SEA and SEB binding to the alpha 2 domain of the H-2Db molecule induces biological activity and has physiological consequences.

Life-span, T-cell responses, and incidence of lymphomas in congenic mice.

Survival, T-cell functions, and postmortem histopathology were studied in H-2 congenic strains of mice bearing H-2b, H-2k, and H-2d haplotypes. Males lived longer than females in all homozygous and heterozygous combinations except for H-2d homozygotes, which showed no differences between males and females. Association of heterozygosity with longer survival was observed only with H-2b/H-2b and H-2b/H-2d mice. Analysis using classification and regression trees (CART) showed that both males and females of H-2b homozygous and H-2k/H-2b mice had the shortest life-span of the strains studied. In histopathological analyses, lymphomas were noted to be more frequent in females, while hemangiosarcomas and hepatomas were more frequent in males. Lymphomas appeared earlier than hepatomas or hemangiosarcomas. The incidence of lymphomas was associated with the H-2 haplotype--e.g., H-2b homozygous mice had more lymphomas than did mice of the H-2d haplotype. More vigorous T-cell function was maintained with age (27 months) in H-2d, H-2b/H-2d, and H-2d/H-2k mice as compared with H-2b, H-2k, and H-2b/H-2k mice, which showed a decline of T-cell responses with age.

Beta diversity at multiple hierarchical levels: explaining the high diversity of scarab beetles in tropical montane forests

Aim: High gamma diversity in tropical montane forests may be ascribed to high geographical turnover of community composition, resulting from population isolation that leads to speciation. We studied the evolutionary processes responsible for diversity and turnover in assemblages of tropical scarab beetles (Scarabaeidae) by assessing DNA sequence variation at multiple hierarchical levels. Location: A 300-km transect across six montane forests (900–1100 m) in Costa Rica. Methods: Assemblages of Scarabaeidae (subfamilies Dynastinae, Rutelinae, Melolonthinae) including 118 morphospecies and > 500 individuals were sequenced for the cox1 gene to establish species limits with a mixed Yule–coalescent method. A species-level phylogenetic tree was constructed from cox1 and rrnL genes. Total diversity and turnover among assemblages were then assessed at three hierarchical levels: haplotypes, species and higher clades. Results: DNA-based analyses showed high turnover among communities at all hierarchical levels. Turnover was highest at the haplotype level (community similarity 0.02–0.12) and decreased with each step of the hierarchy (species: 0.21–0.46; clades: 0.41–0.43). Both compositional and phylogenetic similarities of communities were geographically structured, but turnover was not correlated with distance among forests. When three major clades were investigated separately, communities of Dynastinae showed consistently higher alpha diversity, larger species ranges and lower turnover than Rutelinae and Melolonthinae. Main conclusions: Scarab communities of montane forests show evidence of evolutionary persistence of communities in relative isolation, presumably tracking suitable habitats elevationally to accommodate climatic changes. Patterns of diversity on all hierarchical levels seem to be determined by restricted dispersal, and differences in Dynastinae could be explained by their greater dispersal ability. Community-wide DNA sequencing across multiple lineages and hierarchical levels reveals the evolutionary processes that led to high beta diversity in tropical montane forests through time.

Genetic association of CD247 (CD3ζ) with SLE in a large-scale multiethnic study

A classic T-cell phenotype in systemic lupus erythematosus (SLE) is the downregulation and replacement of the CD3ζ chain that alters T-cell receptor signaling. However, genetic associations with SLE in the human CD247 locus that encodes CD3ζ are not well established and require replication in independent cohorts. Our aim was therefore to examine, localize and validate CD247-SLE association in a large multiethnic population. We typed 44 contiguous CD247 single-nucleotide polymorphisms (SNPs) in 8922 SLE patients and 8077 controls from four ethnically distinct populations. The strongest associations were found in the Asian population (11 SNPs in intron 1, 4.99 × 10(-4) < P < 4.15 × 10(-2)), where we further identified a five-marker haplotype (rs12141731-rs2949655-rs16859085-rs12144621-rs858554; G-G-A-G-A; P(hap) = 2.12 × 10(-5)) that exceeded the most associated single SNP rs858554 (minor allele frequency in controls = 13%; P = 4.99 × 10(-4), odds ratio = 1.32) in significance. Imputation and subsequent association analysis showed evidence of association (P < 0.05) at 27 additional SNPs within intron 1. Cross-ethnic meta-analysis, assuming an additive genetic model adjusted for population proportions, showed five SNPs with significant P-values (1.40 × 10(-3) < P< 3.97 × 10(-2)), with one (rs704848) remaining significant after Bonferroni correction (P(meta) = 2.66 × 10(-2)). Our study independently confirms and extends the association of SLE with CD247, which is shared by various autoimmune disorders and supports a common T-cell-mediated mechanism.

Multiple regulatory variants located in cell type-specific enhancers within the PKP2 locus form major risk and protective haplotypes for canine atopic dermatitis in German shepherd dogs.

BACKGROUND Canine atopic dermatitis (CAD) is a chronic inflammatory skin disease triggered by allergic reactions involving IgE antibodies directed towards environmental allergens. We previously identified a ~1.5 Mb locus on canine chromosome 27 associated with CAD in German shepherd dogs (GSDs). Fine-mapping indicated association closest to the PKP2 gene encoding plakophilin 2. RESULTS Additional genotyping and association analyses in GSDs combined with control dogs from five breeds with low-risk for CAD revealed the top SNP 27:19,086,778 (p = 1.4 × 10(-7)) and a rare ~48 kb risk haplotype overlapping the PKP2 gene and shared only with other high-risk CAD breeds. We selected altogether nine SNPs (four top-associated in GSDs and five within the ~48 kb risk haplotype) that spanned ~280 kb forming one risk haplotype carried by 35 % of the GSD cases and 10 % of the GSD controls (OR = 5.1, p = 5.9 × 10(-5)), and another haplotype present in 85 % of the GSD cases and 98 % of the GSD controls and conferring a protective effect against CAD in GSDs (OR = 0.14, p = 0.0032). Eight of these SNPs were analyzed for transcriptional regulation using reporter assays where all tested regions exerted regulatory effects on transcription in epithelial and/or immune cell lines, and seven SNPs showed allelic differences. The DNA fragment with the top-associated SNP 27:19,086,778 displayed the highest activity in keratinocytes with 11-fold induction of transcription by the risk allele versus 8-fold by the control allele (pdifference = 0.003), and also mapped close (~3 kb) to an ENCODE skin-specific enhancer region. CONCLUSIONS Our experiments indicate that multiple CAD-associated genetic variants located in cell type-specific enhancers are involved in gene regulation in different cells and tissues. No single causative variant alone, but rather multiple variants combined in a risk haplotype likely contribute to an altered expression of the PKP2 gene, and possibly nearby genes, in immune and epithelial cells, and predispose GSDs to CAD.