Genomic patterns of homozygosity in worldwide human populations

Genome-wide patterns of homozygosity runs and their variation across individuals provide a valuable and often untapped resource for studying human genetic diversity and evolutionary history. Using genotype data at 577,489 autosomal SNPs, we employed a likelihood-based approach to identify runs of homozygosity (ROH) in 1,839 individuals representing 64 worldwide populations, classifying them by length into three classes—short, intermediate, and long—with a model-based clustering algorithm. For each class, the number and total length of ROH per individual show considerable variation across individuals and populations. The total lengths of short and intermediate ROH per individual increase with the distance of a population from East Africa, in agreement with similar patterns previously observed for locus-wise homozygosity and linkage disequilibrium. By contrast, total lengths of long ROH show large interindividual variations that probably reflect recent inbreeding patterns, with higher values occurring more often in populations with known high frequencies of consanguineous unions. Across the genome, distributions of ROH are not uniform, and they have distinctive continental patterns. ROH frequencies across the genome are correlated with local genomic variables such as recombination rate, as well as with signals of recent positive selection. In addition, long ROH are more frequent in genomic regions harboring genes associated with autosomal-dominant diseases than in regions not implicated in Mendelian diseases. These results provide insight into the way in which homozygosity patterns are produced, and they generate baseline homozygosity patterns that can be used to aid homozygosity mapping of genes associated with recessive diseases.     

Comparison of SNPs and microsatellites for assessing the genetic structure of chicken populations

Many studies in human genetics compare informativeness of single‐nucleotide polymorphisms (SNPs) and microsatellites (single sequence repeats; SSR) in genome scans, but it is difficult to transfer the results directly to livestock because of different population structures. The aim of this study was to determine the number of SNPs needed to obtain the same differentiation power as with a given standard set of microsatellites. Eight chicken breeds were genotyped for 29 SSRs and 9216 SNPs. After filtering, only 2931 SNPs remained. The differentiation power was evaluated using two methods: partitioning of the Euclidean distance matrix based on a principal component analysis (PCA) and a Bayesian model‐based clustering approach. Generally, with PCA‐based partitioning, 70 SNPs provide a comparable resolution to 29 SSRs. In model‐based clustering, the similarity coefficient showed significantly higher values between repeated runs for SNPs compared to SSRs. For the membership coefficients, reflecting the proportion to which a fraction segment of the genome belongs to the ith cluster, the highest values were obtained for 29 SSRs and 100 SNPs respectively. With a low number of loci (29 SSRs or ≤100 SNPs), neither marker types could detect the admixture in the Gödöllö Nhx population. Using more than 250 SNPs allowed a more detailed insight into the genetic architecture. Thus, the admixed population could be detected. It is concluded that breed differentiation studies will substantially gain power even with moderate numbers of SNPs.     

Biomarker responses in human populations: towards a worldwide map.

The discipline of epidemiology studies the determinants of diseases in human populations, identifies causes, determines outcomes and develops prevention strategies. Traditional epidemiology is most useful for studies of acute, relatively common diseases with short incubation periods but less so for studies of chronic low incidence diseases with long incubation periods. Molecular epidemiology, which employs biological responses or biomarkers as surrogates of exposures or effects, can help with the latter. For this reason, there is a great interest in developing and validating biomarkers. DNA damage underlies an important group of chronic diseases with long incubation periods, i.e., cancer. Biomarkers may measure the exposures that induce the DNA damage, the damage itself, or individual susceptibility to damage. Before they can be used for human population research, however, these measures must be validated. Biomarker validation critically depends on field studies. This is accomplished through transitional epidemiological studies that 'bridge the gap' between laboratory and field. Transitional epidemiological studies are of three varieties: (i) Developmental, (ii) Characterization, and (iii) Applied. Biomarkers are the dependent variables in transitional studies. An international network of laboratories for human population monitoring requires yet another dimension for validation, i.e., the comparability of results among laboratories must be determined. This will be achieved by sample sharing projects, with workshops to compare results. Only then can results in one population be compared with results in another. Interlaboratory standardization of assays for biomarkers validated by transitional studies will have far-reaching benefits. It will allow development of worldwide databases of background values for the various biomarkers-or biomarker maps. This, in turn, will facilitate problem identification and eventually constitute the baselines for area-specific population monitoring. Biomarker databases so developed can be compared with worldwide databases for cancer and heritable diseases, validating the former as statistical surrogates of the latter.     

Possible ancestral structure in human populations

Determining the evolutionary relationships between fossil hominid groups such as Neanderthals and modern humans has been a question of enduring interest in human evolutionary genetics. Here we present a new method for addressing whether archaic human groups contributed to the modern gene pool (called ancient admixture), using the patterns of variation in contemporary human populations. Our method improves on previous work by explicitly accounting for recent population history before performing the analyses. Using sequence data from the Environmental Genome Project, we find strong evidence for ancient admixture in both a European and a West African population (p ≈ 10⁻⁷), with contributions to the modern gene pool of at least 5%. While Neanderthals form an obvious archaic source population candidate in Europe, there is not yet a clear source population candidate in West Africa.     

Genome-wide SNPs detect fine-scale genetic structure in threatened populations of squirrel glider Petaurus norfolcensis

Australian arboreal mammals are experiencing significant population declines, particularly due to land clearing and resulting habitat fragmentation. The squirrel glider, Petaurus norfolcensis, is a threatened species in New South Wales, with a stronghold population in the Lake Macquarie Local Government Area (LGA) where fragmentation due to urbanization is an ongoing problem for the species conservation. Here we report on the use of squirrel glider mitochondrial (385 bp cytochrome b gene, 70 individuals) and nuclear DNA (6,834 SNPs, 87 individuals) markers to assess their population genetic structure and connectivity across 14 locations sampled in the Lake Macquarie LGA. The mitochondrial DNA sequences detected evidence of a historical genetic bottleneck, while the genome-wide SNPs detected significant population structure in the Lake Macquarie squirrel glider populations at scales as fine as one kilometer. There was no evidence of inbreeding within patches, however there were clear effects of habitat fragmentation and biogeographical barriers on gene flow. A least cost path analysis identified thin linear corridors that have high priority for conservation. These areas should be protected to avoid further isolation of squirrel glider populations and the loss of genetic diversity through genetic drift.

     

A population-genetic perspective on the similarities and differences among worldwide human populations

Recent studies have produced a variety of advances in the investigation of genetic similarities and differences among human populations. Here, I pose a series of questions about human population-genetic similarities and differences, and I then answer these questions by numerical computation with a single shared population-genetic data set. The collection of answers obtained provides an introductory perspective for understanding key results on the features of worldwide human genetic variation.     

Ecologically and evolutionarily important SNPs identified in natural populations

Evolution by natural selection acts on natural populations amidst migration, gene-by-environmental interactions, constraints, and tradeoffs, which affect the rate and frequency of adaptive change. We asked how many and how rapidly loci change in populations subject to severe, recent environmental changes. To address these questions, we used genomic approaches to identify randomly selected single nucleotide polymorphisms (SNPs) with evolutionarily significant patterns in three natural populations of Fundulus heteroclitus that inhabit and have adapted to highly polluted Superfund sites. Three statistical tests identified 1.4-2.5% of SNPs that were significantly different from the neutral model in each polluted population. These nonneutral patterns in populations adapted to highly polluted environments suggest that these loci or closely linked loci are evolving by natural selection. One SNP identified in all polluted populations using all tests is in the gene for the xenobiotic metabolizing enzyme, cytochrome P4501A (CYP1A), which has been identified previously as being refractory to induction in the three highly polluted populations. Extrapolating across the genome, these data suggest that rapid evolutionary change in natural populations can involve hundreds of loci, a few of which will be shared in independent events.     

Allele frequencies of 19 autosomal STR loci in Manchu population of China with phylogenetic structure among worldwide populations

The aim of this study was to estimate the allelic frequencies of the 19 STR loci with the Goldeneye™ DNA ID system 20A kit in a sample of 150 Manchu individuals from China to be used for forensic purposes and population studies. The observed heterozygosity(HO)values of these 19 STR loci ranged from 0.600 (D3S1358) to 0.914 (D18S51), the expected (HE) ranged from 0.615 (TPOX) to 0.876 (D16S1043). The power of discrimination (PD) values were found to range from 0.793 (TPOX) to 0.950 (D16S1043) and the probability of exclusion (PE) varies between 0.291 (D3S1358) and 0.825 (D18S51 and Penta E). Among all the 19 loci, D16S1043 had the highest polymorphism (PIC = 0.860), whereas TPOX had the lowest (PIC = 0.550). For the 19 loci, the combined power of discrimination and the combined probability of exclusion are 0.9999999999999999999942 and 0.999999996777, respectively. The phylogenetic tree established among worldwide population shows different populations who say the same language usually have a close genetic relationship with each other across the three language families studied (Sino-Tibetan, Altaic and Arabic).     

Assortative mating: Biodemographical structure of human populations

Concerning marriage, norms are existing in the collective conscience. These norms, to which everyone will try to agree, may vary in function of the societies, but in any way the choice of a partner will be always limited. The number of the potential partners is very limited. The assortative mating of the anthropological characters, studied in samples with limited geographical and socio-economic variability, is not always positive: the variability of the criteria at the level of the choice of a partner and the variability of the social norms relative to the marriage, can explain the variability observed in the coefficients of correlation. This variability relative to the assortative mating of physical characters is chiefly evident in the studies of non-European and non-industrial populations. In the European societies, even when the sample has a reduced variability, the assortative mating is on the contrary highly positive.     

Genomewide identification and annotation of SNPs in Bubalus bubalis

The present study was carried out to identify and annotate the genome wide SNPs in Murrah buffalo genome. A total of 21.2 million raw reads from 4 pooled female Murrah buffalo samples were obtained using restriction enzyme digestion followed by sequencing with Illumina Hiseq 2000. After quality filtration, the reads were aligned to Murrah buffalo genome (ICAR-NBAGR) and Water buffalo genome (UMD_CASPUR_WB_2.0) which resulted in 99.37% and 99.67% of the reads aligning, respectively. A total of 130,688 high quality SNPs along with 35,110 indels were identified versus the Murrah bufffalo genome. Similarly 219,856 high quality SNPs along with 15,201 indels were identified versus the Water buffalo genome. We report 483 SNPs in 66 genes affecting Milk Production, 436 SNPs in 38 genes affecting fertility and 559 SNPs in 72 genes affecting other major traits. The average genome coverage was 13.4% and 14.8% versus the Murrah and Water buffalo genomes, respectively.     

Linkage disequilibria and haplotype structure of four SNPs of the interleukin 1 gene cluster in seven Asian Indian populations

Variation at four single nucleotide polymorphism (SNP) sites of the interleukin 1 (IL1) gene cluster was investigated among 280 unrelated individuals, representing 7 caste groups from the state of Karnataka, India, and one European American community of Boston, Massachusetts. Allele and haplotype frequencies, strength of linkage disequilibrium, and signatures of recombination varied considerably among populations. Variable community sizes and traditions of consanguinity may account for the observed variation.     

Fixation biases affecting human SNPs

Under neutrality all classes of mutation have an equal probability of becoming fixed in a population. In this article, we describe our analysis of the frequency distributions of >5000 human SNPs and provide evident of biases in the process of fixation of certain classes of point mutation that are most likely to be attributable to biased gene conversion. The results indicate an increased fixation probability of mutations that result in the incorporation of a GC base pair. Furthermore, in transcribed regions this process exhibits strand asymmetry, and is biased towards preserving a G base on the coding strand. Biased gene conversion has the potential to explain both existence of isochores and the compositional asymmetry in mammalian transcribed regions.