Genome-wide scans for footprints of natural selection

Molecular signatures of lineage-specific adaptive evolution in a unique sea basin: the example of an anadromous goby Leucopsarion petersii

Climate changes on various time scales often shape genetic novelty and adaptive variation in many biotas. We explored molecular signatures of directional selection in populations of the ice goby Leucopsarion petersii inhabiting a unique sea basin, the Sea of Japan, where a wide variety of environments existed in the Pleistocene in relation to shifts in sea level by repeated glaciations. This species consisted of two historically allopatric lineages, the Japan Sea (JS) and Pacific Ocean (PO) lineages, and these have lived under contrasting marine environments that are expected to have imposed different selection regimes caused by past climatic and current oceanographic factors. We applied a limited genome-scan approach using seven candidate genes for phenotypic differences between two lineages in combination with 100 anonymous microsatellite loci. Neuropeptide Y (NPY) gene, which is an important regulator of food intake and potent orexigenic agent, and three anonymous microsatellites were identified as robust outliers, that is, candidate loci potentially under directional selection, by multiple divergence- and diversity-based outlier tests in comparisons focused on multiple populations of the JS vs. PO lineages. For these outlier loci, populations of the JS lineage had putative signals of selective sweeps. Additionally, real-time quantitative PCR analysis using fish reared in a common environment showed a higher expression level for NPY gene in the JS lineage. Thus, this study succeeded in identifying candidate genomic regions under selection across populations of the JS lineage and provided evidence for lineage-specific adaptive evolution in this unique sea basin.

Conservation genomics of threatened animal species

The genomics era has opened up exciting possibilities in the field of conservation biology by enabling genomic analyses of threatened species that previously were limited to model organisms. Next-generation sequencing (NGS) and the collection of genome-wide data allow for more robust studies of the demographic history of populations and adaptive variation associated with fitness and local adaptation. Genomic analyses can also advance management efforts for threatened wild and captive populations by identifying loci contributing to inbreeding depression and disease susceptibility, and predicting fitness consequences of introgression. However, the development of genomic tools in wild species still carries multiple challenges, particularly those associated with computational and sampling constraints. This review provides an overview of the most significant applications of NGS and the implications and limitations of genomic studies in conservation.

Evidence for selection at HIV host susceptibility genes in a West Central African human population

Background

HIV-1 derives from multiple independent transfers of simian immunodeficiency virus (SIV) strains from chimpanzees to human populations. We hypothesized that human populations in west central Africa may have been exposed to SIV prior to the pandemic, and that previous outbreaks may have selected for genetic resistance to immunodeficiency viruses. To test this hypothesis, we examined the genomes of Biaka Western Pygmies, who historically resided in communities within the geographic range of the central African chimpanzee subspecies (Pan troglodytes troglodytes) that carries strains of SIV ancestral to HIV-1.

Results

SNP genotypes of the Biaka were compared to those of African human populations who historically resided outside the range of P. t. troglodytes, including the Mbuti Eastern Pygmies. Genomic regions showing signatures of selection were compared to the genomic locations of genes reported to be associated with HIV infection or pathogenesis. In the Biaka, a strong signal of selection was detected at CUL5, which codes for a component of the vif-mediated APOBEC3 degradation pathway. A CUL5 allele protective against AIDS progression was fixed in the Biaka. A signal of selection was detected at TRIM5, which codes for an HIV post-entry restriction factor. A protective mis-sense mutation in TRIM5 had the highest frequency in Biaka compared to other African populations, as did a protective allele for APOBEC3G, which codes for an anti-HIV-1 restriction factor. Alleles protective against HIV-1 for APOBEC3H, CXCR6 and HLA-C were at higher frequencies in the Biaka than in the Mbuti. Biaka genomes showed a strong signal of selection at TSG101, an inhibitor of HIV-1 viral budding.

Conclusions

We found protective alleles or evidence for selection in the Biaka at a number of genes associated with HIV-1 infection or progression. Pygmies have also been reported to carry genotypes protective against HIV-1 for the genes CCR5 and CCL3L1. Our hypothesis that HIV-1 may have shaped the genomes of some human populations in West Central Africa appears to merit further investigation.

Population genomic scan for candidate signatures of balancing selection to guide antigen characterization in malaria parasites

Acquired immunity in vertebrates maintains polymorphisms in endemic pathogens, leading to identifiable signatures of balancing selection. To comprehensively survey for genes under such selection in the human malaria parasite Plasmodium falciparum, we generated paired-end short-read sequences of parasites in clinical isolates from an endemic Gambian population, which were mapped to the 3D7 strain reference genome to yield high-quality genome-wide coding sequence data for 65 isolates. A minority of genes did not map reliably, including the hypervariable var, rifin, and stevor families, but 5,056 genes (90.9% of all in the genome) had >70% sequence coverage with minimum read depth of 5 for at least 50 isolates, of which 2,853 genes contained 3 or more single nucleotide polymorphisms (SNPs) for analysis of polymorphic site frequency spectra. Against an overall background of negatively skewed frequencies, as expected from historical population expansion combined with purifying selection, the outlying minority of genes with signatures indicating exceptionally intermediate frequencies were identified. Comparing genes with different stage-specificity, such signatures were most common in those with peak expression at the merozoite stage that invades erythrocytes. Members of clag, PfMC-2TM, surfin, and msp3-like gene families were highly represented, the strongest signature being in the msp3-like gene PF10_0355. Analysis of msp3-like transcripts in 45 clinical and 11 laboratory adapted isolates grown to merozoite-containing schizont stages revealed surprisingly low expression of PF10_0355. In diverse clonal parasite lines the protein product was expressed in a minority of mature schizonts (<1% in most lines and ∼10% in clone HB3), and eight sub-clones of HB3 cultured separately had an intermediate spectrum of positive frequencies (0.9 to 7.5%), indicating phase variable expression of this polymorphic antigen. This and other identified targets of balancing selection are now prioritized for functional study.

Microsatellites as targets of natural selection

The ability to survey polymorphism on a genomic scale has enabled genome-wide scans for the targets of natural selection. Theory that connects patterns of genetic variation to evidence of natural selection most often assumes a diallelic locus and no recurrent mutation. Although these assumptions are suitable to selection that targets single nucleotide variants, fundamentally different types of mutation generate abundant polymorphism in genomes. Moreover, recent empirical results suggest that mutationally complex, multiallelic loci including microsatellites and copy number variants are sometimes targeted by natural selection. Given their abundance, the lack of inference methods tailored to the mutational peculiarities of these types of loci represents a notable gap in our ability to interrogate genomes for signatures of natural selection. Previous theoretical investigations of mutation-selection balance at multiallelic loci include assumptions that limit their application to inference from empirical data. Focusing on microsatellites, we assess the dynamics and population-level consequences of selection targeting mutationally complex variants. We develop general models of a multiallelic fitness surface, a realistic model of microsatellite mutation, and an efficient simulation algorithm. Using these tools, we explore mutation-selection-drift equilibrium at microsatellites and investigate the mutational history and selective regime of the microsatellite that causes Friedreich's ataxia. We characterize microsatellite selective events by their duration and cost, note similarities to sweeps from standing point variation, and conclude that it is premature to label microsatellites as ubiquitous agents of efficient adaptive change. Together, our models and simulation algorithm provide a powerful framework for statistical inference, which can be used to test the neutrality of microsatellites and other multiallelic variants.

Shared ancestral susceptibility to colorectal cancer and other nutrition related diseases

Background

The majority of non-syndromic colorectal cancers (CRCs) can be described as a complex disease. A two-stage case–control study on CRC susceptibility was conducted to assess the influence of the ancestral alleles in the polymorphisms previously associated with nutrition-related complex diseases.

Methods

In stage I, 28 single nucleotide polymorphisms (SNPs) were genotyped in a hospital-based Czech population (1025 CRC cases, 787 controls) using an allele-specific PCR-based genotyping system (KASPar®). In stage II, replication was carried out for the five SNPs with the lowest p values. The replication set consisted of 1798 CRC cases and 1810 controls from a population-based German study (DACHS). Odds ratios (ORs) and 95% confidence intervals (CIs) for associations between genotypes and CRC risk were estimated using logistic regression. To identify signatures of selection, Fay-Wu’s H and Integrated Haplotype Score (iHS) were estimated.

Results

In the Czech population, carriers of the ancestral alleles of AGT rs699 and CYP3A7 rs10211 showed an increased risk of CRC (OR 1.26 and 1.38, respectively; two-sided p≤0.05), whereas carriers of the ancestral allele of ENPP1 rs1044498 had a decreased risk (OR 0.79; p≤0.05). For rs1044498, the strongest association was detected in the Czech male subpopulation (OR 0.61; p=0.0015). The associations were not replicated in the German population. Signatures of selection were found for all three analyzed genes.

Conclusions

Our study showed evidence of association for the ancestral alleles of polymorphisms in AGT and CYP3A7 and for the derived allele of a polymorphism in ENPP1 with an increased risk of CRC in Czechs, but not in Germans. The ancestral alleles of these SNPs have previously been associated with nutrition-related diseases hypertension (AGT and CYP3A7) and insulin resistance (ENPP1). Future studies may shed light on the complex genetic and environmental interactions between different types of nutrition-related diseases.

Comparative analysis of FLC homologues in Brassicaceae provides insight into their role in the evolution of oilseed rape

We identified nine FLOWERING LOCUS C homologues (BnFLC) in Brassica napus and found that the coding sequences of all BnFLCs were relatively conserved but the intronic and promoter regions were more divergent. The BnFLC homologues were mapped to six of 19 chromosomes. All of the BnFLC homologues were located in the collinear region of FLC in the Arabidopsis genome except BnFLC.A3b and BnFLC.C3b, which were mapped to noncollinear regions of chromosome A3 and C3, respectively. Four of the homologues were associated significantly with quantitative trait loci for flowering time in two mapping populations. The BnFLC homologues showed distinct expression patterns in vegetative and reproductive organs, and at different developmental stages. BnFLC.A3b was differentially expressed between the winter-type and semi-winter-type cultivars. Microsynteny analysis indicated that BnFLC.A3b might have been translocated to the present segment in a cluster with other flowering-time regulators, such as a homologue of FRIGIDA in Arabidopsis. This cluster of flowering-time genes might have conferred a selective advantage to Brassica species in terms of increased adaptability to diverse environments during their evolution and domestication process.

Human pigmentation genes under environmental selection

Genome-wide association studies and comparative genomics have established major loci and specific polymorphisms affecting human skin, hair and eye color. Environmental changes have had an impact on selected pigmentation genes as populations have expanded into different regions of the globe.

Genome patterns of selection and introgression of haplotypes in natural populations of the house mouse (Mus musculus)

General parameters of selection, such as the frequency and strength of positive selection in natural populations or the role of introgression, are still insufficiently understood. The house mouse (Mus musculus) is a particularly well-suited model system to approach such questions, since it has a defined history of splits into subspecies and populations and since extensive genome information is available. We have used high-density single-nucleotide polymorphism (SNP) typing arrays to assess genomic patterns of positive selection and introgression of alleles in two natural populations of each of the subspecies M. m. domesticus and M. m. musculus. Applying different statistical procedures, we find a large number of regions subject to apparent selective sweeps, indicating frequent positive selection on rare alleles or novel mutations. Genes in the regions include well-studied imprinted loci (e.g. Plagl1/Zac1), homologues of human genes involved in adaptations (e.g. alpha-amylase genes) or in genetic diseases (e.g. Huntingtin and Parkin). Haplotype matching between the two subspecies reveals a large number of haplotypes that show patterns of introgression from specific populations of the respective other subspecies, with at least 10% of the genome being affected by partial or full introgression. Using neutral simulations for comparison, we find that the size and the fraction of introgressed haplotypes are not compatible with a pure migration or incomplete lineage sorting model. Hence, it appears that introgressed haplotypes can rise in frequency due to positive selection and thus can contribute to the adaptive genomic landscape of natural populations. Our data support the notion that natural genomes are subject to complex adaptive processes, including the introgression of haplotypes from other differentiated populations or species at a larger scale than previously assumed for animals. This implies that some of the admixture found in inbred strains of mice may also have a natural origin.

A critical assessment of storytelling: gene ontology categories and the importance of validating genomic scans

In the age of whole-genome population genetics, so-called genomic scan studies often conclude with a long list of putatively selected loci. These lists are then further scrutinized to annotate these regions by gene function, corresponding biological processes, expression levels, or gene networks. Such annotations are often used to assess and/or verify the validity of the genome scan and the statistical methods that have been used to perform the analyses. Furthermore, these results are frequently considered to validate "true-positives" if the identified regions make biological sense a posteriori. Here, we show that this approach can be potentially misleading. By simulating neutral evolutionary histories, we demonstrate that it is possible not only to obtain an extremely high false-positive rate but also to make biological sense out of the false-positives and construct a sensible biological narrative. Results are compared with a recent polymorphism data set from Drosophila melanogaster.

Methods to detect selection on noncoding DNA

Vast tracts of noncoding DNA contain elements that regulate gene expression in higher eukaryotes. Describing these regulatory elements and understanding how they evolve represent major challenges for biologists. Advances in the ability to survey genome-scale DNA sequence data are providing unprecedented opportunities to use evolutionary models and computational tools to identify functionally important elements and the mode of selection acting on them in multiple species. This chapter reviews some of the current methods that have been developed and applied on noncoding DNA, what they have shown us, and how they are limited. Results of several recent studies reveal that a significantly larger fraction of noncoding DNA in eukaryotic organisms is likely to be functional than previously believed, implying that the functional annotation of most noncoding DNA in these organisms is largely incomplete. In Drosophila, recent studies have further suggested that a large fraction of noncoding DNA divergence observed between species may be the product of recurrent adaptive substitution. Similar studies in humans have revealed a more complex pattern, with signatures of recurrent positive selection being largely concentrated in conserved noncoding DNA elements. Understanding these patterns and the extent to which they generalize to other organisms awaits the analysis of forthcoming genome-scale polymorphism and divergence data from more species.

Genome-wide detection of natural selection in African Americans pre- and post-admixture

It is particularly meaningful to investigate natural selection in African Americans (AfA) due to the high mortality their African ancestry has experienced in history. In this study, we examined 491,526 autosomal single nucleotide polymorphisms (SNPs) genotyped in 5210 individuals and conducted a genome-wide search for selection signals in 1890 AfA. Several genomic regions showing an excess of African or European ancestry, which were considered the footprints of selection since population admixture, were detected based on a commonly used approach. However, we also developed a new strategy to detect natural selection both pre- and post-admixture by reconstructing an ancestral African population (AAF) from inferred African components of ancestry in AfA and comparing it with indigenous African populations (IAF). Interestingly, many selection-candidate genes identified by the new approach were associated with AfA-specific high-risk diseases such as prostate cancer and hypertension, suggesting an important role these disease-related genes might have played in adapting to a new environment. CD36 and HBB, whose mutations confer a degree of protection against malaria, were also located in the highly differentiated regions between AAF and IAF. Further analysis showed that the frequencies of alleles protecting against malaria in AAF were lower than those in IAF, which is consistent with the relaxed selection pressure of malaria in the New World. There is no overlap between the top candidate genes detected by the two approaches, indicating the different environmental pressures AfA experienced pre- and post-population admixture. We suggest that the new approach is reasonably powerful and can also be applied to other admixed populations such as Latinos and Uyghurs.

Genomic toolboxes for conservation biologists

Conservation genetics is expanding its research horizon with a genomic approach, by incorporating the modern techniques of next-generation sequencing (NGS). Application of NGS overcomes many limitations of conservation genetics. First, NGS allows for genome-wide screening of markers, which may lead to a more representative estimation of genetic variation within and between populations. Second, NGS allows for distinction between neutral and non-neutral markers. By screening populations on thousands of single nucleotide polymorphism markers, signals of selection can be found for some markers. Variation in these markers will give insight into functional rather than neutral genetic variation. Third, NGS facilitates the study of gene expression. Conservation genomics will increase our insight in how the environment and genes interact to affect phenotype and fitness. In addition, the NGS approach opens a way to study processes such as inbreeding depression and local adaptation mechanistically. Conservation genetics programs are directed to a fundamental understanding of the processes involved in conservation genetics and should preferably be started in species for which large databases on ecology, demography and genetics are available. Here, we describe and illustrate the connection between the application of NGS technologies and the research questions in conservation. The perspectives of conservation genomics programs are also discussed.

Recent advances in understanding the role of nutrition in human genome evolution

Dietary transitions in human history have been suggested to play important roles in the evolution of mankind. Genetic variations caused by adaptation to diet during human evolution could have important health consequences in current society. The advance of sequencing technologies and the rapid accumulation of genome information provide an unprecedented opportunity to comprehensively characterize genetic variations in human populations and unravel the genetic basis of human evolution. Series of selection detection methods, based on various theoretical models and exploiting different aspects of selection signatures, have been developed. Their applications at the species and population levels have respectively led to the identification of human specific selection events that distinguish human from nonhuman primates and local adaptation events that contribute to human diversity. Scrutiny of candidate genes has revealed paradigms of adaptations to specific nutritional components and genome-wide selection scans have verified the prevalence of diet-related selection events and provided many more candidates awaiting further investigation. Understanding the role of diet in human evolution is fundamental for the development of evidence-based, genome-informed nutritional practices in the era of personal genomics.

A general mechanistic model for admixture histories of hybrid populations

Admixed populations have been used for inferring migrations, detecting natural selection, and finding disease genes. These applications often use a simple statistical model of admixture rather than a modeling perspective that incorporates a more realistic history of the admixture process. Here, we develop a general model of admixture that mechanistically accounts for complex historical admixture processes. We consider two source populations contributing to the ancestry of a hybrid population, potentially with variable contributions across generations. For a random individual in the hybrid population at a given point in time, we study the fraction of genetic admixture originating from a specific one of the source populations by computing its moments as functions of time and of introgression parameters. We show that very different admixture processes can produce identical mean admixture proportions, but that such processes produce different values for the variance of the admixture proportion. When introgression parameters from each source population are constant over time, the long-term limit of the expectation of the admixture proportion depends only on the ratio of the introgression parameters. The variance of admixture decreases quickly over time after the source populations stop contributing to the hybrid population, but remains substantial when the contributions are ongoing. Our approach will facilitate the understanding of admixture mechanisms, illustrating how the moments of the distribution of admixture proportions can be informative about the historical admixture processes contributing to the genetic diversity of hybrid populations.

Determinants of hierarchical genetic structure in Atlantic salmon populations: environmental factors vs. anthropogenic influences

Disentangling the effects of natural environmental features and anthropogenic factors on the genetic structure of endangered populations is an important challenge for conservation biology. Here, we investigated the combined influences of major environmental features and stocking with non-native fish on the genetic structure and local adaptation of Atlantic salmon (Salmo salar) populations. We used 17 microsatellite loci to genotype 975 individuals originating from 34 French rivers. Bayesian analyses revealed a hierarchical genetic structure into five geographically distinct clusters. Coastal distance, geological substrate and river length were strong predictors of population structure. Gene flow was higher among rivers with similar geologies, suggesting local adaptation to geological substrate. The effect of river length was mainly owing to one highly differentiated population that has the farthest spawning grounds off the river mouth (up to 900km) and the largest fish, suggesting local adaptation to river length. We detected high levels of admixture in stocked populations but also in neighbouring ones, implying large-scale impacts of stocking through dispersal of non-native individuals. However, we found relatively few admixed individuals suggesting a lower fitness of stocked fish and/or some reproductive isolation between wild and stocked individuals. When excluding stocked populations, genetic structure increased as did its correlation with environmental factors. This study overall indicates that geological substrate and river length are major environmental factors influencing gene flow and potential local adaptation among Atlantic salmon populations but that stocking with non-native individuals may ultimately disrupt these natural patterns of gene flow among locally adapted populations.

In silico mining and characterization of 12 EST-SSRs for the invasive slipper limpet Crepidula fornicata

In silico mining of an expressed sequence tags (ESTs) library was found to be efficient at isolating simple sequence repeats (SSRs) loci in the non-indigenous marine mollusc Crepidula fornicata. Twelve SSR loci were developed for routine genotyping. Cross-species amplification to 8 other Crepidula species showed that the 12 loci are highly specific for C. fornicata. Mendelian inheritance was shown for 11 of them (1 being monomorphic in the analyzed offspring array). The genetic diversity for 88 adults was found to be variable across the 12 loci (2-40 alleles, expected heterozygosity between 0.023 and 0.898) with a high overall exclusion probability of 0.99. The degree of genetic polymorphism found here is similar to that shown for 7 anonymous SSRs previously developed and here used on the same samples. This set of 12 specific loci is relevant to perform reliable population and relatedness analyses in Crepidula fornicata.

Gene gain and loss during evolution of obligate parasitism in the white rust pathogen of Arabidopsis thaliana

Biotrophic eukaryotic plant pathogens require a living host for their growth and form an intimate haustorial interface with parasitized cells. Evolution to biotrophy occurred independently in fungal rusts and powdery mildews, and in oomycete white rusts and downy mildews. Biotroph evolution and molecular mechanisms of biotrophy are poorly understood. It has been proposed, but not shown, that obligate biotrophy results from (i) reduced selection for maintenance of biosynthetic pathways and (ii) gain of mechanisms to evade host recognition or suppress host defence. Here we use Illumina sequencing to define the genome, transcriptome, and gene models for the obligate biotroph oomycete and Arabidopsis parasite, Albugo laibachii. A. laibachii is a member of the Chromalveolata, which incorporates Heterokonts (containing the oomycetes), Apicomplexa (which includes human parasites like Plasmodium falciparum and Toxoplasma gondii), and four other taxa. From comparisons with other oomycete plant pathogens and other chromalveolates, we reveal independent loss of molybdenum-cofactor-requiring enzymes in downy mildews, white rusts, and the malaria parasite P. falciparum. Biotrophy also requires “effectors” to suppress host defence; we reveal RXLR and Crinkler effectors shared with other oomycetes, and also discover and verify a novel class of effectors, the “CHXCs”, by showing effector delivery and effector functionality. Our findings suggest that evolution to progressively more intimate association between host and parasite results in reduced selection for retention of certain biosynthetic pathways, and particularly reduced selection for retention of molybdopterin-requiring biosynthetic pathways. These mechanisms are not only relevant to plant pathogenic oomycetes but also to human pathogens within the Chromalveolata.

Plant pathogens that cannot grow except on their hosts are called obligate biotrophs. How such biotrophy evolves is poorly understood. In this study, we sequenced the genome of the obligate biotroph white rust pathogen (Albugo laibachii, Oomycota) of Arabidopsis. From comparisons with other oomycete plant pathogens, diatoms, and the human pathogen Plasmodium falciparum, we reveal a loss of important metabolic enzymes. We also reveal the appearance of defence-suppressing “effectors”, some carrying motifs known from other oomycete effectors, and discover and experimentally verify a novel class of effectors that share a CHXC motif within 50 amino acids of the signal peptide cleavage site. Obligate biotrophy involves an intimate association within host cells at the haustorial interface (where the parasite penetrates the host cell's cell wall), where nutrients are acquired from the host and effectors are delivered to the host. We found that A. laibachii, like Hyaloperonospora arabidopsidis and Plasmodium falciparum, lacks molybdopterin-requiring biosynthetic pathways, suggesting relaxed selection for retention of, or even selection against, this pathway. We propose that when defence suppression becomes sufficiently effective, hosts become such a reliable source of nutrients that a free-living phase can be lost. These mechanisms leading to obligate biotrophy and host specificity are relevant not only to plant pathogenic oomycetes but also to human pathogens.

Analysis of >1000 single nucleotide polymorphisms in geographically matched samples of landrace and wild barley indicates secondary contact and chromosome-level differences in diversity around domestication genes

Barley is a model species for the investigation of the evolution, adaptation and spread of the world's important crops. In this article, we describe the first application of an oligonucleotide pool assay single nucleotide polymorphism (SNP) platform to assess the evolution of barley in a portion of the Fertile Crescent, a key region in the development of farming. A large collection of >1000 genetically mapped, genome-wide SNPs was assayed in geographically matched landrace and wild barley accessions (N=448) from Jordan and Syria. Landrace and wild barley categories were clearly genetically differentiated, but a limited degree of secondary contact was evident. Significant chromosome-level differences in diversity between barley types were observed around genes known to be involved in the evolution of cultivars. The region of Jordan and southern Syria, compared with the north of Syria, was supported by SNP data as a more likely domestication origin. Our data provide evidence for hybridization as a possible mechanism for the continued adaptation of landrace barley under cultivation, indicate regions of the genome that may be subject to selection processes and suggest limited origins for the development of the cultivated crop.

Applications and implications of neutral versus non-neutral markers in molecular ecology

The field of molecular ecology has expanded enormously in the past two decades, largely because of the growing ease with which neutral molecular genetic data can be obtained from virtually any taxonomic group. However, there is also a growing awareness that neutral molecular data can provide only partial insight into parameters such as genetic diversity, local adaptation, evolutionary potential, effective population size, and taxonomic designations. Here we review some of the applications of neutral versus adaptive markers in molecular ecology, discuss some of the advantages that can be obtained by supplementing studies of molecular ecology with data from non-neutral molecular markers, and summarize new methods that are enabling researchers to generate data from genes that are under selection.

Integration of the genetic map and genome assembly of fugu facilitates insights into distinct features of genome evolution in teleosts and mammals

The compact genome of fugu (Takifugu rubripes) has been used widely as a reference genome for understanding the evolution of vertebrate genomes. However, the fragmented nature of the fugu genome assembly has restricted its use for comparisons of genome architecture in vertebrates. To extend the contiguity of the assembly to the chromosomal level, we have generated a comprehensive genetic map of fugu and anchored the scaffolds of the assembly to the 22 chromosomes of fugu. The map consists of 1,220 microsatellite markers that provide anchor points to 697 scaffolds covering 86% of the genome assembly (http://www.fugu-sg.org/). The integrated genome map revealed a higher recombination rate in fugu compared with other vertebrates and a wide variation in the recombination rate between sexes and across chromosomes of fugu. We used the extended assembly to explore recent rearrangement events in the lineages of fugu, Tetraodon, and medaka and compared them with rearrangements in three mammalian (human, mouse, and opossum) lineages. Between the two pufferfishes, fugu has experienced fewer chromosomal rearrangements than Tetraodon. The gene order is more highly conserved in the three teleosts than in mammals largely due to a lower rate of interchromosomal rearrangements in the teleosts. These results provide new insights into the distinct patterns of genome evolution between teleosts and mammals. The consolidated genome map and the genetic map of fugu are valuable resources for comparative genomics of vertebrates and for elucidating the genetic basis of the phenotypic diversity of ~25 species of Takifugu that evolved within the last 5 My.

The population genetics of chronic kidney disease: insights from the MYH9-APOL1 locus

Many rare kidney disorders exhibit a monogenic, Mendelian pattern of inheritance. Population-based genetic studies have identified many genetic variants associated with an increased risk of developing common kidney diseases. Strongly associated variants have potential clinical uses as predictive markers and may advance our understanding of disease pathogenesis. These principles are elegantly illustrated by a region within chromosome 22q12 that has a strong association with common forms of kidney disease. Researchers had identified DNA sequence variants in this locus that were highly associated with an increased prevalence of common chronic kidney diseases in people of African ancestry. Initial research concentrated on MYH9 as the most likely candidate gene; however, population-based whole-genome analysis enabled two independent research teams to discover more strongly associated mutations in the neighboring APOL1 gene. The powerful evolutionary selection pressure of an infectious pathogen in West Africa favored the spread of APOL1 variants that protect against a lethal form of African sleeping sickness but are highly associated with an increased risk of kidney disease. We describe the data sources, process of discovery, and reasons for initial misidentification of the candidate gene, as well as the lessons that can be learned for future population genetics research.

Duffy-null-associated low neutrophil counts influence HIV-1 susceptibility in high-risk South African black women

BACKGROUND

The Duffy-null trait and ethnic netropenia are both highly prevalent in Africa. The influence of pre-seroconversion levels of peripheral blood cell counts (PBCs) on the risk of acquiring human immunodeficiency virus (HIV)-1 infection among Africans is unknown.

METHODS

The triangular relationship among pre-seroconversion PBC counts, host genotypes, and risk of HIV acquisition was determined in a prospective cohort of black South African high-risk female sex workers. Twenty-seven women had seroconversion during follow-up, and 115 remained HIV negative for 2 years, despite engaging in high-risk activity.

RESULTS

Pre-seroconversion neutrophil counts in women who subsequently had seroconversion were significantly lower, whereas platelet counts were higher, compared with those who remained HIV negative. Comprising 27% of the cohort, subjects with pre-seroconversion neutrophil counts of <2500 cells/mm(3) had a ∼3-fold greater risk of acquiring HIV infection. In a genome-wide association analyses, an African-specific polymorphism (rs2814778) in the promoter of Duffy Antigen Receptor for Chemokines (DARC -46T > C) was significantly associated with neutrophil counts (P = 7.9 × 10(-11)). DARC -46C/C results in loss of DARC expression on erthyrocytes (Duffy-null) and resistance to Plasmodium vivax malaria, and in our cohort, only subjects with this genotype had pre-seroconversion neutrophil counts of <2500 cells/mm(3). The risk of acquiring HIV infection was ∼3-fold greater in those with the trait of Duffy-null-associated low neutrophil counts, compared with all other study participants.

CONCLUSIONS

Pre-seroconversion neutrophil and platelet counts influence risk of HIV infection. The trait of Duffy-null-associated low neutrophil counts influences HIV susceptibility. Because of the high prevalence of this trait among persons of African ancestry, it may contribute to the dynamics of the HIV epidemic in Africa.

A regression-based approach to selection mapping

Selection mapping applies the population genetics theory of hitchhiking to the localization of genomic regions containing genes under selection. This approach predicts that neutral loci linked to genes under positive selection will have reduced diversity due to their shared history with a selected locus, and thus, genome scans of diversity levels can be used to identify regions containing selected loci. Most previous approaches to this problem ignore the spatial genomic pattern of diversity expected under selection. The regression-based approach advocated in this paper takes into account the expected pattern of decreasing genetic diversity with increased proximity to a selected locus. Simulated data are used to examine the patterns of diversity under different scenarios, in order to assess the power of a regression-based approach to the identification of regions under selection. Application of this method to both simulated and empirical data demonstrates its potential to detect selection. In contrast to some other methods, the regression approach described in this paper can be applied to any marker type. Results also suggest that this approach may give more precise estimates of the location of the selected locus than alternative methods, although the power is slightly lower in some cases.

Estimating the correlation of non-allele descents along chromosomes

The pattern of the correlation of non-allele descents among linked sites is an important aspect for an insight into the genomic evolution at the population level. Here, we present a new statistical method for estimating two types of non-allele descent correlations. One is the standardized parental descent disequilibrium termed by Cockerham & Weir (1973), the other is the standardized disequilibrium between non-allele descent segments from the same chromosome. Essential to this analysis is the partitioning of the joint identity-by-state probability for a random pair of non-allele gametes into the different components of identity by descents at the two or three sites. We consider the samples of phased haplotypes of single nucleotide polymorphism (SNP) markers and the weighted least square method for fast parameter estimation. Monte Carlo simulations demonstrate that robustly unbiased estimates with appropriate precisions can be obtained with certain sample sizes, ~100 diploids, under the impacts of allele frequency distributions and linkage disequilibrium. This method can be used to construct the maps of non-allele descent correlation blocks for the population whose genetic pedigree is not required on a prior basis.

Comparative genetic approaches to the evolution of human brain and behavior

With advances in genomic technologies, the amount of genetic data available to scientists today is vast. Genomes are now available or planned for 14 different primate species and complete resequencing of numerous human individuals from numerous populations is underway. Moreover, high-throughput deep sequencing is quickly making whole genome efforts within the reach of single laboratories allowing for unprecedented studies. Comparative genetic approaches to the identification of the underlying basis of human brain, behavior, and cognitive ability are moving to the forefront. Two approaches predominate: inter-species divergence comparisons and intra-species polymorphism studies. These methodological differences are useful for different time scales of evolution and necessarily focus on different evolutionary events in the history of primate and hominin evolution. Inter-species divergence is more useful in studying large scale primate, or hominoid, evolution whereas intra-species polymorphism can be more illuminating of recent hominin evolution. These differences in methodological utility also extend to studies of differing genetic substrates; current divergence studies focus primarily on protein evolution whereas polymorphism studies are substrate ambivalent. Some of the issues inherent in these studies can be ameliorated by current sequencing capabilities whereas others remain intractable. New avenues are also being opened that allow for the incorporation of novel substrates and approaches. In the post-genomic era, the study of human evolution, specifically as it relates to the brain, is becoming more complete focusing increasingly on the totality of the system and better conceptualizing the entirety of the genetic changes that have lead to the human phenotype today.

Using human demographic history to infer natural selection reveals contrasting patterns on different families of immune genes

Detecting regions of the human genome that are, or have been, influenced by natural selection remains an important goal for geneticists. Many methods are used to infer selection, but there is a general reliance on an accurate understanding of how mutation and recombination events are distributed, and the well-known link between these processes and their evolutionary transience introduces uncertainty into inferences. Here, we present and apply two new, independent approaches; one based on single nucleotide polymorphisms (SNPs) that exploits geographical patterns in how humans lost variability as we colonized the world, the other based on the relationship between microsatellite repeat number and heterozygosity. We show that the two methods give concordant results. Of these, the SNP-based method is both widely applicable and detects selection over a well-defined time interval, the last 50 000 years. Analysis of all human genes by their Gene Ontology codes reveals how accelerated and decelerated loss of variability are both preferentially associated with immune genes. Applied to 168 immune genes used as the focus of a previous study, we show that members of the same gene family tend to yield similar indices of selection, even when located on different chromosomes. We hope our approach will provide a useful tool with which to infer where selection has acted to shape the human genome.

USING POPULATION GENOMICS TO DETECT SELECTION IN NATURAL POPULATIONS: KEY CONCEPTS AND METHODOLOGICAL CONSIDERATIONS

Natural selection shapes patterns of genetic variation among individuals, populations, and species, and it does so differentially across genomes. The field of population genomics provides a comprehensive genome-scale view of the action of selection, even beyond traditional model organisms. However, even with nearly complete genomic sequence information, our ability to detect the signature of selection on specific genomic regions depends on choosing experimental and analytical tools appropriate to the biological situation. For example, processes that occur at different timescales, such as sorting of standing genetic variation, mutation-selection balance, or fixed interspecific divergence, have different consequences for genomic patterns of variation. Inappropriate experimental or analytical approaches may fail to detect even strong selection or falsely identify a signature of selection. Here we outline the conceptual framework of population genomics, relate genomic patterns of variation to evolutionary processes, and identify major biological factors to be considered in studies of selection. As data-gathering technology continues to advance, our ability to understand selection in natural populations will be limited more by conceptual and analytical weaknesses than by the amount of molecular data. Our aim is to bring critical biological considerations to the fore in population genomics research and to spur the development and application of analytical tools appropriate to diverse biological systems.

Worldwide distribution of the MYH9 kidney disease susceptibility alleles and haplotypes: evidence of historical selection in Africa

MYH9 was recently identified as renal susceptibility gene (OR 3-8, p < 10(-8)) for major forms of kidney disease disproportionately affecting individuals of African descent. The risk haplotype (E-1) occurs at much higher frequencies in African Americans (> or = 60%) than in European Americans (< 4%), revealing a genetic basis for a major health disparity. The population distributions of MYH9 risk alleles and the E-1 risk haplotype and the demographic and selective forces acting on the MYH9 region are not well explored. We reconstructed MYH9 haplotypes from 4 tagging single nucleotide polymorphisms (SNPs) spanning introns 12-23 using available data from HapMap Phase II, and by genotyping 938 DNAs from the Human Genome Diversity Panel (HGDP). The E-1 risk haplotype followed a cline, being most frequent within sub-Saharan African populations (range 50-80%), less frequent in populations from the Middle East (9-27%) and Europe (0-9%), and rare or absent in Asia, the Americas, and Oceania. The fixation indexes (F(ST)) for pairwise comparisons between the risk haplotypes for continental populations were calculated for MYH9 haplotypes; F(ST) ranged from 0.27-0.40 for Africa compared to other continental populations, possibly due to selection. Uniquely in Africa, the Yoruba population showed high frequency extended haplotype length around the core risk allele (C) compared to the alternative allele (T) at the same locus (rs4821481, iHs = 2.67), as well as high population differentiation (F(ST(CEU vs. YRI)) = 0.51) in HapMap Phase II data, also observable only in the Yoruba population from HGDP (F(ST) = 0.49), pointing to an instance of recent selection in the genomic region. The population-specific divergence in MYH9 risk allele frequencies among the world's populations may prove important in risk assessment and public health policies to mitigate the burden of kidney disease in vulnerable populations.

Population structure with localized haplotype clusters

We propose a multilocus version of F(ST) and a measure of haplotype diversity using localized haplotype clusters. Specifically, we use haplotype clusters identified with BEAGLE, which is a program implementing a hidden Markov model for localized haplotype clustering and performing several functions including inference of haplotype phase. We apply this methodology to HapMap phase 3 data. With this haplotype-cluster approach, African populations have highest diversity and lowest divergence from the ancestral population, East Asian populations have lowest diversity and highest divergence, and other populations (European, Indian, and Mexican) have intermediate levels of diversity and divergence. These relationships accord with expectation based on other studies and accepted models of human history. In contrast, the population-specific F(ST) estimates obtained directly from single-nucleotide polymorphisms (SNPs) do not reflect such expected relationships. We show that ascertainment bias of SNPs has less impact on the proposed haplotype-cluster-based F(ST) than on the SNP-based version, which provides a potential explanation for these results. Thus, these new measures of F(ST) and haplotype-cluster diversity provide an important new tool for population genetic analysis of high-density SNP data.