Characterisation of genome-wide association epistasis signals for serum uric acid in human population isolates

Genome-wide association (GWA) studies have identified a number of loci underlying variation in human serum uric acid (SUA) levels with the SLC2A9 gene having the largest effect identified so far. Gene-gene interactions (epistasis) are largely unexplored in these GWA studies. We performed a full pair-wise genome scan in the Italian MICROS population (n = 1201) to characterise epistasis signals in SUA levels. In the resultant epistasis profile, no SNP pairs reached the Bonferroni adjusted threshold for the pair-wise genome-wide significance. However, SLC2A9 was found interacting with multiple loci across the genome, with NFIA - SLC2A9 and SLC2A9 - ESRRAP2 being significant based on a threshold derived for interactions between GWA significant SNPs and the genome and jointly explaining 8.0% of the phenotypic variance in SUA levels (3.4% by interaction components). Epistasis signal replication in a CROATIAN population (n = 1772) was limited at the SNP level but improved dramatically at the gene ontology level. In addition, gene ontology terms enriched by the epistasis signals in each population support links between SUA levels and neurological disorders. We conclude that GWA epistasis analysis is useful despite relatively low power in small isolated populations.     

UniFrag and GenomePrimer: selection of primers for genome-wide production of unique amplicons

The complementary programs UniFrag and GenomePrimer were developed to provide a reliable high-throughput method to select the most unique regions within genomic DNA sequence(s) and design primers therein, involving minimal user intervention and maximum flexibility.     

Delineation of signals required for thymocyte positive selection

Peptide/MHC complexes capable of inducing positive selection in mouse fetal thymic organ cultures fail to do so in suspension culture. Furthermore, this type of culture does not promote initial stages of differentiation, such as coreceptor down-modulation, unless peptides used for stimulation have (at least) weak agonist activity. We show in this study that signals provided in suspension culture by nonagonist peptide/MHC complexes on the surface of macrophages, even though apparently silent, are sufficient to promote complete phenotypic differentiation when CD4+CD8+ thymocytes are subsequently placed in a proper anatomical setting. Furthermore, the synergistic actions of suboptimal concentrations of phorbol esters and nonagonist peptide/MHC complexes can make the initial stages of positive selection visible, without converting maturation into negative selection. Thus, the correlation between efficiency of positive selection and the degree of coreceptor down-modulation on CD4+CD8+ thymocytes is not linear. Furthermore, these results suggest that the unique role of thymic stromal cells in positive selection is related not to presentation of self-peptide/MHC complexes, but most likely to another ligand.     

Sex-specific population structure,natural selection,and linkage disequilibrium in a wild bird population as revealed by genome-wide microsatellite analyses

Background  

Sexual dimorphism in ecologically important traits is widespread, yet the differences in the genomic architecture between the two sexes are largely unexplored. We employed a genome-wide multilocus approach to examine the sexual differences in population subdivision, natural selection and linkage disequilibrium (LD) in a wild Siberian jay (Perisoreus infaustus) population, using genotypes at a total of 107 autosomal and Z-chromosomal microsatellites.     

Exploration of signals of positive selection derived from genotype-based human genome scans using re-sequencing data

We have investigated whether regions of the genome showing signs of positive selection in scans based on haplotype structure also show evidence of positive selection when sequence-based tests are applied, whether the target of selection can be localized more precisely, and whether such extra evidence can lead to increased biological insights. We used two tools: simulations under neutrality or selection, and experimental investigation of two regions identified by the HapMap2 project as putatively selected in human populations. Simulations suggested that neutral and selected regions should be readily distinguished and that it should be possible to localize the selected variant to within 40 kb at least half of the time. Re-sequencing of two ~300 kb regions (chr4:158Mb and chr10:22Mb) lacking known targets of selection in HapMap CHB individuals provided strong evidence for positive selection within each and suggested the micro-RNA gene hsa-miR-548c as the best candidate target in one region, and changes in regulation of the sperm protein gene SPAG6 in the other.     

Haplotype structure and evidence for positive selection at the human IL13 locus

Interleukin-13 (IL13) is believed to play an important role in the pathogenesis of atopy and allergic asthma. To better understand genetic variation at the IL13 locus, we resequenced a 5.1-kb genomic region spanning the entire locus and identified 26 single-nucleotide polymorphisms (SNPs) in 74 individuals from three major populations-Chinese, Caucasian, and African. Our survey suggests exceptionally high and significant geographic structure at the IL13 locus between African and outside Africa populations. This unusual pattern suggests that positive selection that acts in some local populations may have played a role on the IL13 locus. In support of this suggestion, we found a significant excess of high frequency-derived SNPs in the Chinese population and Caucasian population, respectively, as expected after a recent episode of positive selection. Further, the unusual haplotype structure indicates that different scenarios of the action of positive selection on the IL13 locus in different populations may exist. In the Caucasian population, the skewed haplotype distribution dominated by one common haplotype supports the hypothesis of simple directional selection. Whereas, in the Chinese population, the two-round hitchhiking hypothesis may explain the skewed haplotype structure with three dominant ones. These findings may provide insight into the likely relative roles of selection and population history in establishing present-day variation at the IL13 locus, and, motivate further studies of this locus as an important candidate in common diseases association studies.     

Mitochondrial cytochrome b sequence variations and population structure of Siberian chipmunk (Tamias sibiricus) in Northeastern Asia and population substructure in South Korea

Twenty-five chipmunk species occur in the world, of which only the Siberian chipmunk, Tamias sibiricus, inhabits Asia. To investigate mitochondrial cytochrome b sequence variations and population structure of the Siberian chipmunk in northeastern Asia, we examined mitochondrial cytochrome b sequences (1140 bp) from 3 countries. Analyses of 41 individuals from South Korea and 33 individuals from Russia and northeast China resulted in 37 haplotypes and 27 haplotypes, respectively. There were no shared haplotypes between South Korea and Russia--northeast China. Phylogenetic trees and network analysis showed 2 major maternal lineages for haplotypes, referred to as the S and R lineages. Haplotype grouping in each cluster was nearly coincident with its geographic affinity. In particular, 3 distinct groups were found that mostly clustered in the northern, central and southern parts of South Korea. Nucleotide diversity of the S lineage was twice that of lineage R. The divergence between S and R lineages was estimated to be 2.98-0.98 Myr. During the ice age, there may have been at least 2 refuges in South Korea and Russia--northeast China. The sequence variation between the S and R lineages was 11.3% (K2P), which is indicative of specific recognition in rodents. These results suggest that T. sibiricus from South Korea could be considered a separate species. However, additional information, such as details of distribution, nuclear genes data or morphology, is required to strengthen this hypothesis.     

Evidence for recent positive selection at the human AIM1 locus in a European population

Two missense polymorphisms (E272K and L374F) of the AIM1 locus,encoding a melanocyte differentiation antigen, were shown tohave a clear association with human ethnicities. These two nonpathogenicsingle nucleotide polymorphisms (SNPs) may be associated withhuman pigmentation variation. In this study, we investigatedsequence variation in the coding region and exon-flanking sequenceand found low genetic variation only in subjects of Europeandescent. All four statistical tests applied to the 7.55-kb regionsurrounding the L374F polymorphism detected statistically significantdeviations from selective neutrality in Europeans. In addition,haplotype analysis revealed that one haplotype carrying 374Fwas overrepresented in this population, and the low rate ofvariation, with some features of selective sweeps, was shownto be statistically significant. These results suggest thatpositive selection recently has been acting or has acted onat least this region of the melanogenic gene and that an advantageoushaplotype spread rapidly in Europe.     

Shared and unique features of evolutionary diversification

A fundamental question in evolutionary biology asks whether organisms experiencing similar selective pressures will evolve similar solutions or whether historical contingencies dominate the evolutionary process and yield disparate evolutionary outcomes. It is perhaps most likely that both shared selective forces as well as unique histories play key roles in the course of evolution. Consequently, when multiple species face a common environmental gradient, their patterns of divergence might exhibit both shared and unique elements. Here we describe a general framework for investigating and evaluating the relative importance of these contrasting features of diversification. We examined morphological diversification in three species of livebearing fishes across a predation gradient. All species (Gambusia affinis from the United States of America, Brachyrhaphis rhabdophora from Costa Rica, and Poecilia reticulata from Trinidad) exhibited more elongate bodies, a larger caudal peduncle, and a relatively lower position of the eye in predator populations. This shared response suggests that common selective pressures generated parallel outcomes within three different species. However, each species also exhibited unique features of divergence, which might reflect phylogenetic tendencies, chance events, or localized environmental differences. In this system, we found that shared aspects of divergence were of larger magnitude than unique elements, suggesting common natural selective forces have played a greater role than unique histories in producing the observed patterns of morphological diversification. Assessing the nature and relative importance of shared and unique responses should aid in elucidating the relative generality or peculiarity in evolutionary divergence.     

A genome-wide metabolic QTL analysis in Europeans implicates two loci shaped by recent positive selection

We have performed a metabolite quantitative trait locus (mQTL) study of the (1)H nuclear magnetic resonance spectroscopy ((1)H NMR) metabolome in humans, building on recent targeted knowledge of genetic drivers of metabolic regulation. Urine and plasma samples were collected from two cohorts of individuals of European descent, with one cohort comprised of female twins donating samples longitudinally. Sample metabolite concentrations were quantified by (1)H NMR and tested for association with genome-wide single-nucleotide polymorphisms (SNPs). Four metabolites' concentrations exhibited significant, replicable association with SNP variation (8.6×10(-11)相似文献   

Dissecting the genetic structure of Korean population using genome-wide SNP arrays

Genome-wide SNP arrays have generated unprecedented quantities of data allow the detection of human evolutionary history and dense genome-wide data also enable the identification of distance ancestry among individuals or ethnic groups. To explain wider aspects of the genetic structure of Koreans and the East Asian population, we analyzed 79 individuals from the Korean HapMap project at 555,352 common single-nucleotide polymorphism loci, and compared this data with the worldwide population groups with the 53 ethnic groups from Human Genome Diversity Panel (HGDP-CEPH). Population differentiation (F_ST), Principal Component Analyses, STRUCTURE and ADMIXTURE are examined. In general, all the individual samples studies here were classified into subset of ethnic groups according to their geographical origins. Korean HapMap individuals were grouped together with East Asian populations from HGDP panel. Recently, a sub-population structure within Korean population has been reported. Our result, however, revealed the genetic homogeneity of Korean population. The ADMIXTURE analysis showed that, overall the Korean populations derive 79 % of their genomic ancestry from southern Asia and have relatively little northern Asian ancestry (21 %). The present work, therefore, provide the evidence that the male-biased southern-to-northern migration influenced not only for the genetic make up of the Y chromosome in the Korean population but also, its autosomal composition.     

Population structure and population history of Anopheles dirus mosquitoes in Southeast Asia

Separating the confounding effects of long-term population history from gene flow can be difficult. Here, we address the question of what inferences about gene flow can be made from mitochondrial sequence data in three closely related species of mosquitoes, Anopheles dirus species A, C, and D, from southeast Asia. A total of 84 sequences of 923 bp of the mitochondrial cytochrome oxidase I gene were obtained from 14 populations in Thailand, Myanmar, and Bangladesh. The genealogy of sequences obtained from two populations of AN: dirus C indicates no contemporary gene flow between them. The F(ST) value of 0.421 therefore probably represents a recent common history, perhaps involving colonization events. Anopheles dirus A and D are parapatric, yet no differentiation was seen either within or between species. The starlike genealogy of their haplotypes, smooth unimodal mismatch distributions, and excess of low frequency mutations indicate population expansion in An. dirus A and D. This, rather than widespread gene flow, explains their low within-species F(ST) values (0.018 and 0.022). The greater genetic diversity of An. dirus D suggests that expansion occurred first in species D and subsequently in species A. The current geographical separation and low hybrid fitness of these species also argue against ongoing interspecific gene flow. They suggest instead either historical introgression of mtDNA from An. dirus D into species A followed by independent range expansions, or a selective sweep of mtDNA that originated in An. dirus D. While not excluding contemporary gene flow, historical population processes are sufficient to explain the data in An. dirus A and D. The genealogical relationships between haplotypes could not be used to make inferences of gene flow because of extensive homoplasy due to hypervariable sites and possibly also recombination. However, it is concluded that this approach, rather than the use of fixation indices, is required in the future to understand contemporary gene flow in these mosquitoes. The implications of these results for understanding gene flow in another important and comparable group of malaria vector mosquitoes in Africa, the An. gambiae complex, are also discussed.     

Contrasting genetic structure in Plasmodium vivax populations from Asia and South America

Populations of Plasmodium falciparum show striking differences in linkage disequilibrium, population differentiation and diversity, but only fragmentary data exists on the genetic structure of Plasmodium vivax. We genotyped nine tandem repeat loci bearing 2-8 bp motifs from 345 P. vivax infections collected from three Asian countries and from five locations in Colombia. We observed 9-37 alleles per locus and high diversity (He=0.72-0.79, mean=0.75) in all countries. Numbers of multiple clone infections varied considerably: these were rare in Colombia and India, but > 60% of isolates carried multiple alleles in at least one locus in Thailand and Laos. However, only one or two of the nine loci show >1 allele in many samples, suggesting that mutation within infections may result in overestimation of true multiple carriage rates. Identical nine-locus genotypes were frequently found in Colombian populations, contributing to strong linkage disequilibrium. These identical genotypes were strongly clustered in time, consistent with epidemic transmission of clones and subsequent breakdown of allelic associations, suggesting high rates of inbreeding and low effective recombination rates in this country. In contrast, identical genotypes were rare and loci were randomly associated in all three Asian populations, consistent with higher rates of outcrossing and recombination. We observed low but significant differentiation between different Asian countries (standardized FST = 0.13-0.45). In comparison, we see greater differentiation between collection locations within Colombia (standardized FST = 0.4-0.7), and strong differentiation between continents (standardized FST = 0.48-0.79). The observed heterogeneity in multiple clone carriage rates, linkage disequilibrium and population differentiation are similar in some, but not all, respects to those observed in P. falciparum, and have important implications for the design of association mapping studies, and interpretation of P. vivax epidemiology.     

Shared and unique patterns of embryo development in extremophile poeciliids

Background

Closely related lineages of livebearing fishes have independently adapted to two extreme environmental factors: toxic hydrogen sulphide (H₂S) and perpetual darkness. Previous work has demonstrated in adult specimens that fish from these extreme habitats convergently evolved drastically increased head and offspring size, while cave fish are further characterized by reduced pigmentation and eye size. Here, we traced the development of these (and other) divergent traits in embryos of Poecilia mexicana from benign surface habitats (“surface mollies”) and a sulphidic cave (“cave mollies”), as well as in embryos of the sister taxon, Poecilia sulphuraria from a sulphidic surface spring (“sulphur mollies”). We asked at which points during development changes in the timing of the involved processes (i.e., heterochrony) would be detectible.

Methods and Results

Data were extracted from digital photographs taken of representative embryos for each stage of development and each type of molly. Embryo mass decreased in convergent fashion, but we found patterns of embryonic fat content and ovum/embryo diameter to be divergent among all three types of mollies. The intensity of yellow colouration of the yolk (a proxy for carotenoid content) was significantly lower in cave mollies throughout development. Moreover, while relative head size decreased through development in surface mollies, it increased in both types of extremophile mollies, and eye growth was arrested in mid-stage embryos of cave mollies but not in surface or sulphur mollies.

Conclusion

Our results clearly demonstrate that even among sister taxa convergence in phenotypic traits is not always achieved by the same processes during embryo development. Furthermore, teleost development is crucially dependent on sufficient carotenoid stores in the yolk, and so we discuss how the apparent ability of cave mollies to overcome this carotenoid-dependency may represent another potential mechanism explaining the lack of gene flow between surface and cave mollies.     

Estimation of life cycle components of selection in an experimental plant population

Viability and fertility components of selection associated with linkage blocks marked by four electrophoretically detectable loci were estimated in an experimental population of barley [Composite Cross V (CCV)]. The intensity of selection affecting the distribution of pollen types in the outcross pool was also estimated and comparisons were made between the selective values of genes in the pools of uniting ovules and pollen. The estimates show that selection was intense at various stages of the life cycle and that viability and fertility components often opposed one another. Estimates of viability and fertility components of selection were also extended to the three-locus level. The multilocus estimates reveal large differences in viability and fertility among homozygous genotypes. It is concluded that strong selection operates at all life cycle stages in CCV, although often in differing directions.     

Low-pass genome-wide sequencing and variant inference using identity-by-descent in an isolated human population

Whole-genome sequencing in an isolated population with few founders directly ascertains variants from the population bottleneck that may be rare elsewhere. In such populations, shared haplotypes allow imputation of variants in unsequenced samples without resorting to complex statistical methods as in studies of outbred cohorts. We focus on an isolated population cohort from the Pacific Island of Kosrae, Micronesia, where we previously collected SNP array and rich phenotype data for the majority of the population. We report identification of long regions with haplotypes co-inherited between pairs of individuals and methodology to leverage such shared genetic content for imputation. Our estimates show that sequencing as few as 40 personal genomes allows for inference in up to 60% of the 3000-person cohort at the average locus. We ascertained a pilot data set of whole-genome sequences from seven Kosraean individuals, with average 5× coverage. This assay identified 5,735,306 unique sites of which 1,212,831 were previously unknown. Additionally, these variants are unusually enriched for alleles that are rare in other populations when compared to geographic neighbors (published Korean genome SJK). We used the presence of shared haplotypes between the seven Kosraen individuals to estimate expected imputation accuracy of known and novel homozygous variants at 99.6% and 97.3%, respectively. This study presents whole-genome analysis of a homogenous isolate population with emphasis on optimal rare variant inference.