Haplotype diversity across 100 candidate genes for inflammation, lipid metabolism, and blood pressure regulation in two populations

Recent studies have suggested that a significant fraction of the human genome is contained in blocks of strong linkage disequilibrium, ranging from ~5 to >100 kb in length, and that within these blocks a few common haplotypes may account for >90% of the observed haplotypes. Furthermore, previous studies have suggested that common haplotypes in candidate genes are generally shared across populations and represent the majority of chromosomes in each population. The conclusions drawn from these preliminary studies, however, are based on an incomplete knowledge of the variation in the regions examined. To bridge this gap in knowledge, we have completely resequenced 100 candidate genes in a population of African descent and one of European descent. Although these genes have been well studied because of their medical importance, we demonstrate that a large amount of sequence variation has not yet been described. We also report that the average number of inferred haplotypes per gene, when complete data is used, is higher than in previous reports and that the number and proportion of all haplotypes represented by common haplotypes per gene is variable. Furthermore, we demonstrate that haplotypes shared between the two populations constitute only a fraction of the total number of haplotypes observed and that these shared haplotypes represent fewer of the African-descent chromosomes than was expected from previous studies. Finally, we show that restricting variation discovery to coding regions does not adequately describe all common haplotypes or the true haplotype block structure observed when all common variation is used to infer haplotypes. These data, derived from complete knowledge of genetic variation in these genes, suggest that the haplotype architecture of candidate genes across the human genome is more complex than previously suggested, with important implications for candidate gene and genomewide association studies.     

Haplotype diversity in mtDNA and Y-chromosome in populations of Altai-Sayan region

Polymorphism of mtDNA was examined in five ethnic populations that belong to the Turkic language group and inhabit the territory of the Altai-Sayan upland (N = 1007). Most of the haplogroups identified in the examined populations belonged to East Eurasian lineages. In all five populations, only three haplogroups, C, D, and F, were prevailing. The frequencies of the other six haplogroups (A, B, G, M, Y, and Z) varied in the range from 1.1 to 6.5%. Among West Eurasian haplogrous, the most common were haplogroups H, J, T, and U. An analysis of Y-chromosome haplogroups in 407 individuals showed that only two haplogroups, N* and R1a1, were present in all five populations examined. Moreover, in different ethnic groups, the highest frequencies were observed for C-M130, N-P43, and N-Tat haplogroups. The differences in the distribution patterns of ancient West Eurasian and East Eurasian haplotypes from Gorny Altai in the present-day populations from the northern part of Eurasia revealed can be explained in terms of the multistage expansion of humans across these territories. The ubiquity of haplotypes from haplogroup H and cluster U across the wide territory from the Yenisei River basin to the Atlantic Ocean can indicate directional human expansion, which most likely occurred out of Central Asia as early as in the Paleolithic era, and took place in several waves with the glacier retreat.     

Common miRNAs,candidate genes and their interaction network across four subtypes of epithelial ovarian cancer

Epithelial ovarian cancer (EOC) is categorized into four major histological subtypes such as clear cell carcinoma (CCC), endometrioid carcinoma (EC), mucinous carcinoma (MC), and serous carcinoma (SC). Heterogeneity of the EOC leads to different clinical outcomes of the disease, although all the subtypes are originated from the same layer of tissue. Therefore, it is of interest to identify the common candidate genes, miRNA and their interaction network in four the subtypes of EOC. A comparative gene expression analysis identified 248 common differentially expressed genes (DEGs) in the four subtypes of EOC. Identified common DEGs were found to be enriched in cancer specific pathways. A protein-protein interaction (PPI) network of the common DEGs were constructed, and subsequent module and survival analyses identified seven key candidate genes (CCNB1, CENPM, CEP55, RACGAP1, TPX2, UBE2C, and ZWINT). We also documented 10 key candidate miRNAs (hsa-mir-16-5p, hsa-mir-23b-3p, hsa-mir-34a-5p, hsa-mir-103a-3p, hsa-mir-107, hsa-mir-124-3p, hsa-mir-129-2-3p, hsa-mir-147a, hsa-mir-205-5p, and hsa-mir-195-5p) linked to the candidate genes. These derived data find application in the understanding of EOC.     

2型糖尿病易感候选基因在世界不同人群中的多样性比较分析

近10年来兴起的全基因组关联分析（Genome-wide association study, GWAS）相关研究结果获得了大量与2型糖尿病相关的候选易感基因,了解这些候选基因在正常人群中的遗传多样性程度以及在不同人群间的遗传差异,不但有助于阐明2型糖尿病的遗传机理,而且对于今后在特定人群中进行2型糖尿病发病机制的深入研究具有指导意义。本研究通过对GWAS数据库和相关文献的搜索和整理确定了170个与2型糖尿病相关的基因或基因区域;随后基于千人基因组计划的全基因组测序数据对这些候选基因在世界范围内14个人群间的遗传多样性进行了比较分析;进一步确定了在人群间存在显著差异的易感基因,并分析了这些基因的多样性特征。在所研究的14个世界人群中,2型糖尿病候选易感基因的遗传多样性与基因组范围的平均水平没有显著差异;但其中8个易感基因IL20RA、RNMTL1-NXN、NOTCH2、ADRA2A-BTBD7P2、TBC1D4、RBM38-HMGB1P1、UBE2E2和PPARD在群体间呈现显著差异,其中最明显的是IL20RA基因 (FST=0.152),该易感基因在非洲人群和非非洲人群间存在显著等位基因频率和单倍型频率差异。14个人群中易感基因遗传结构差异的主要原因是由于非洲人群与非非洲人群之间的群体遗传结构的不同所造成的。进一步比较东西方人群间的2型糖尿病候选基因遗传结构差异,发现在东西方人群中同样存在明显的群体遗传结构差别,其中DGKB-AGMO(FST=0.173)和JAZF1(FST=0.182)是差异最显著的易感基因。本研究通过对群体间2型糖尿病易感基因遗传结构进行比较,鉴别出一些差异特别显著的易感基因,对今后2型糖尿病易感基因与不同人群间发病率和易感性差异的相关研究提供重要参考。     

Haplotype distribution and association of candidate genes with salt tolerance in Indian wild rice germplasm

Key message

The association of natural genetic variations of salt-responsive candidate genes belonging to different gene families with salt-tolerance phenotype and their haplotype variation in different geographic regions.

Abstract

Soil salinity covers a large part of the arable land of the world and is a major factor for yield losses in salt-sensitive crops, such as rice. Different gene families that respond to salinity have been identified in rice, but limited success has been achieved in developing salt-tolerant cultivars. Therefore, 21 salt stress-responsive candidate genes belonging to different gene families were re-sequenced to analyse their genetic variation and association with salt tolerance. The average single nucleotide polymorphism (SNP) density was 16 SNPs per kbp amongst these genes. The identified nucleotide and haplotype diversity showed comparatively higher genetic variation in the transporter family genes. Linkage disequilibrium (LD) analysis showed significant associations of SNPs in BADH2, HsfC1B, MIPS1, MIPS2, MYB2, NHX1, NHX2, NHX3, P5CS1, P5CS2, PIP1, SIK1, SOS1, and SOS2 genes with the salt-tolerant phenotype. A combined analysis of SNPs in the 21 candidate genes and eight other HKT transporter genes produced two separate clusters of tolerant genotypes, carrying unique SNPs in the ion transporter and osmoticum-related genes. Haplotype network analysis showed all the major and few minor alleles distributed over distant geographic regions. Minor haplotypes may be recently evolved alleles which migrated to distant geographic regions and may represent recent expansion of Indian wild rice. The analysis of genetic variation in different gene families identified the relationship between adaptive variations and functional significance of the genes. Introgression of the identified alleles from wild relatives may enhance the salt tolerance and consequently rice production in the salinity-affected areas.

     

An extensive candidate gene approach to speciation: diversity,divergence and linkage disequilibrium in candidate pigmentation genes across the European crow hybrid zone

Colouration patterns have an important role in adaptation and speciation. The European crow system, in which all-black carrion crows and grey-coated hooded crows meet in a narrow hybrid zone, is a prominent example. The marked phenotypic difference is maintained by assortative mating in the absence of neutral genetic divergence, suggesting the presence of few pigmentation genes of major effect. We made use of the rich phenotypic and genetic resources in mammals and identified a comprehensive panel of 95 candidate pigmentation genes for birds. Based on functional annotation, we chose a subset of the most promising 37 candidates, for which we developed a marker system that demonstrably works across the avian phylogeny. In total, we sequenced 107 amplicons (∼3 loci per gene, totalling 60 kb) in population samples of crows (n=23 for each taxon). Tajima''s D, Fu''s F_S, DHEW and HKA (Hudson–Kreitman–Aguade) statistics revealed several amplicons that deviated from neutrality; however, none of these showed significantly elevated differentiation between the two taxa. Hence, colour divergence in this system may be mediated by uncharacterized pigmentation genes or regulatory regions outside genes. Alternatively, the observed high population recombination rate (4N_er∼0.03), with overall linkage disequilibrium dropping rapidly within the order of few 100 bp, may compromise the power to detect causal loci with nearby markers. Our results add to the debate as to the utility of candidate gene approaches in relation to genomic features and the genetic architecture of the phenotypic trait in question.     

Sequence diversity in 36 candidate genes for cardiovascular disorders.

Two strategies involving whole-genome association studies have been proposed for the identification of genes involved in complex diseases. The first one seeks to characterize all common variants of human genes and to test their association with disease. The second one seeks to develop dense maps of single-nucleotide polymorphisms (SNPs) and to detect susceptibility genes through linkage disequilibrium. We performed a molecular screening of the coding and/or flanking regions of 36 candidate genes for cardiovascular diseases. All polymorphisms identified by this screening were further genotyped in 750 subjects of European descent. In the whole set of genes, the lengths explored spanned 53.8 kb in the 5' regions, 68.4 kb in exonic regions, and 13 kb in the 3' regions. The strength of linkage disequilibrium within candidate regions suggests that genomewide maps of SNPs might be efficient ways to identify new disease-susceptibility genes, provided that the maps are sufficiently dense. However, the relatively large number of polymorphisms within coding and regulatory regions of candidate genes raises the possibility that several of them might be functional and that the pattern of genotype-phenotype association might be more complex than initially envisaged, as actually has been observed in some well-characterized genes. These results argue in favor of both genomewide association studies and detailed studies of the overall sequence variation of candidate genes, as complementary approaches.     

Ontogenic differences in sexual size dimorphism across four plover populations

Sexual size dimorphism (SSD) among adults is commonly observed in animals and is considered to be adaptive. However, the ontogenic emergence of SSD, i.e. the timing of divergence in body size between males and females, has only recently received attention. It is widely acknowledged that the ontogeny of SSD may differ between species, but it remains unclear how variable the ontogeny of SSD is within species. Kentish Plovers Charadrius alexandrinus and Snowy Plovers C. nivosus are closely related wader species that exhibit similar, moderate (c. 4%), male‐biased adult SSD. To assess when SSD emerges we recorded tarsus length variation among 759 offspring in four populations of these species. Tarsus length of chicks was measured on the day of hatching and up to three times on recapture before fledging. In one population (Mexico, Snowy Plovers), males and females differed in size from the day of hatching, whereas growth rates differed between the sexes in two populations (Turkey and United Arab Emirates, both Kentish Plovers). In contrast, a fourth population (Cape Verde, Kentish Plovers) showed no significant SSD in juveniles. Our results suggest that adult SSD can emerge at different stages of development (prenatal, postnatal and post‐juvenile) in different populations of the same species. We discuss the proximate mechanisms that may underlie these developmental differences.     

Mitochondrial diversity in human head louse populations across the Americas

Anthropological studies suggest that the genetic makeup of human populations in the Americas is the result of diverse processes including the initial colonization of the continent by the first people plus post‐1492 European migrations. Because of the recent nature of some of these events, understanding the geographical origin of American human diversity is challenging. However, human parasites have faster evolutionary rates and larger population sizes allowing them to maintain greater levels of genetic diversity than their hosts. Thus, we can use human parasites to provide insights into some aspects of human evolution that may be unclear from direct evidence. In this study, we analyzed mitochondrial DNA (mtDNA) sequences from 450 head lice in the Americas. Haplotypes clustered into two well‐supported haplogroups, known as A and B. Haplogroup frequencies differ significantly among North, Central and South America. Within each haplogroup, we found evidence of demographic expansions around 16,000 and 20,000 years ago, which correspond broadly with those estimated for Native Americans. The parallel timing of demographic expansions of human lice and Native Americans plus the contrasting pattern between the distribution of haplogroups A and B through the Americas suggests that human lice can provide additional evidence about the human colonization of the New World. Am J Phys Anthropol 152:118–129, 2013. © 2013 Wiley Periodicals, Inc.     

Mapping spot blotch resistance genes in four barley populations

Bipolaris sorokiniana (teleomorph: Cochliobolus sativus) is the fungal pathogen responsible for spot blotch in barley (Hordeum vulgare L.) and occurs worldwide in warmer, humid growing conditions. Current Australian barley varieties are largely susceptible to this disease and attempts are being made to introduce sources of resistance from North America. In this study we have compared chromosomal locations of spot blotch resistance reactions in four North American two-rowed barley lines; the North Dakota lines ND11231-12 and ND11231-11 and the Canadian lines TR251 and WPG8412-9-2-1. Diversity arrays technology-based PCR, expressed sequence tag and SSR markers have been mapped across four populations derived from crosses between susceptible parental lines and these four resistant parents to determine the location of resistance loci. Quantitative trait loci (QTL) conferring resistance to spot blotch in adult plants (APR) were detected on chromosomes 3HS and 7HS. In contrast, seedling resistance (SLR) was controlled solely by a locus on chromosome 7HS. The phenotypic variance explained by the APR QTL on 3HS was between 16 and 25% and the phenotypic variance explained by the 7HS APR QTL was between 8 and 42% across the four populations. The SLR QTL on 7HS explained between 52 and 64% of the phenotypic variance. An examination of the pedigrees of these resistance sources supports the common identity of resistance in these lines and indicates that only a limited number of major resistance loci are available in current two-rowed germplasm.     

Mapping QTLs and candidate genes for rice root traits under different water-supply conditions and comparative analysis across three populations

To investigate the genetic factors underlying constitutive and adaptive morphological traits of roots under different water-supply conditions, a recombinant inbred line (RIL) population derived from a cross between the lowland rice variety IR1552 and the upland rice variety Azucena with 249 molecular markers, was used in cylindrical-pot experiments. Eighteen QTLs were detected for seminal root length (SRL), adventitious root number (ARN), and lateral root length (LRL) and lateral root number (LRN) on the seminal root at a soil depth of from 3 to 6 cm under flooding and upland conditions. One identical QTL was detected under both flooding and upland conditions. The relative parameters under the two water-supply conditions were also used for QTL analysis. Five QTLs for upland induced variations in the traits were detected with the positive alleles from Azucena. A comparative analysis was performed for the QTLs detected in this study and those reported from two other populations with Azucena as a parent. Several identical QTLs for root elongation were found across the three populations with positive alleles from Azucena. Candidate genes were screened from ESTs and cDNA-AFLP clones for comparative mapping with the detected QTLs. Two genes for cell expansion, OsEXP2 and endo-1,4--D-glucanase EGase, and four cDNA-AFLP clones from root tissues of Azucena, were mapped on the intervals carrying the QTLs for SRL and LRL under upland conditions, respectively.Communicated by H.C. Becker     

Haplotype reconstruction for diploid populations

The inference of haplotype pairs directly from unphased genotype data is a key step in the analysis of genetic variation in relation to disease and pharmacogenetically relevant traits. Most popular methods such as Phase and PL do require either the coalescence assumption or the assumption of linkage between the single-nucleotide polymorphisms (SNPs). We have now developed novel approaches that are independent of these assumptions. First, we introduce a new optimization criterion in combination with a block-wise evolutionary Monte Carlo algorithm. Based on this criterion, the 'haplotype likelihood', we develop two kinds of estimators, the maximum haplotype-likelihood (MHL) estimator and its empirical Bayesian (EB) version. Using both real and simulated data sets, we demonstrate that our proposed estimators allow substantial improvements over both the expectation-maximization (EM) algorithm and Clark's procedure in terms of capacity/scalability and error rate. Thus, hundreds and more ambiguous loci and potentially very large sample sizes can be processed. Moreover, applying our proposed EB estimator can result in significant reductions of error rate in the case of unlinked or only weakly linked SNPs.     

Tests of candidate genes in breed cross populations for QTL mapping in livestock

In recent years, several F₂ crosses between outbred lines of livestock have been developed to identify quantitative trait loci (QTL). These populations are valuable for further genetic analysis, including positional candidate gene loci (CGL). Analysis of CGL in F₂ populations is, however, hindered by extensive between-breed linkage disequilibrium (LD). The objectives here were to develop and evaluate three tests for CGL in simulated F₂ breed-cross populations. 1) A standard association test, based on the fixed effect of CGL genotype. This test was significant for CGL at considerable distances from the QTL. 2) A marker-assisted association test, based on a test at the CGL of the fixed effect of CGL genotype in a breed-cross QTL interval mapping model. This removed the impact of between-breed LD, but was not powerful in detecting CGL closely linked to the QTL, unless the CGL was the QTL. 3) An F-drop test, comparing F ratios for a QTL at the CGL with and without the CGL included as fixed effect. It had low power to distinguish close from distant CGL. Power to distinguish two CGL within 10 cM from the QTL was limited and little improved by including QTL effects associated with markers to remove between-breed LD, although the power was greater when one of the CGL was the causative mutation. Therefore, while we conclude that candidate gene tests in QTL mapping populations must be interpreted with caution, we now have a clearer picture of the value of candidate gene tests in these populations.     

Population structure and genetic diversity of four Henan pig populations

To investigate the population structure and genetic diversity of Henan indigenous pig breeds, samples from a total of 78 pigs of 11 breeds were collected, including four pig populations from Henan Province, three Western commercial breeds, three Chinese native pig breeds from other provinces and one Asian wild boar. The genotyping datasets were obtained by genotyping‐by‐sequencing technology. We found a high degree of polymorphism and rapid linkage disequilibrium decay in Henan pigs. A neighbor‐joining tree, principal component analysis and structure analysis revealed that the Huainan and Erhualian pigs were clustered together and that the Queshan black pigs were clearly grouped together but that the Nanyang and Yuxi pigs were extensively admixed with Western pigs. In addition, heterozygosity values might indicate that Henan indigenous pigs, especially the Queshan black and Huainan pigs, were subjected to little selection during domestication. The results presented here indicate that Henan pig breeds were admixed from Western breeds, especially Nanyang and Yuxi pigs. Therefore, establishment of purification and rejuvenation systems to implement conservation strategies is urgent. In addition, it is also necessary to accelerate genetic resources improvement and utilization using modern breeding technologies, such as genomic selection and genome‐wide association studies.     

The contribution of three strong candidate schizophrenia susceptibility genes in demographically distinct populations

Here we characterize and compare the contribution of three recently identified strong candidate schizophrenia susceptibility genes; G72, neuregulin 1 (NRG1) and dystrobrevin-binding protein 1 (DTNBP1) in two independent datasets of patients with distinct genetic backgrounds. On the basis of corrected P-values from single- and multilocus transmission distortion tests our analysis provides no support for a contribution of G72, NRG1 or DTNBP1 in the tested samples. When transmission of individual haplotypes was considered, a picture more consistent with the original studies emerged, where transmission distortions in the same direction as the original samples and involving the same core haplotypes were observed for G72 and NRG1. Interestingly, whereas the NRG1 gene analysis was dominated by the presence of over-transmitted haplotypes, the G72 gene analysis was consistently dominated in both datasets by under-transmissions. Negative transmissions involved a core haplotype complementary to the originally detected over-transmitted haplotype, suggesting the presence of a protective variant within the G72 locus.     

Haplotype inference in crossbred populations without pedigree information

Background

Current methods for haplotype inference without pedigree information assume random mating populations. In animal and plant breeding, however, mating is often not random. A particular form of nonrandom mating occurs when parental individuals of opposite sex originate from distinct populations. In animal breeding this is called crossbreeding and hybridization in plant breeding. In these situations, association between marker and putative gene alleles might differ between the founding populations and origin of alleles should be accounted for in studies which estimate breeding values with marker data. The sequence of alleles from one parent constitutes one haplotype of an individual. Haplotypes thus reveal allele origin in data of crossbred individuals.

Results

We introduce a new method for haplotype inference without pedigree that allows nonrandom mating and that can use genotype data of the parental populations and of a crossbred population. The aim of the method is to estimate line origin of alleles. The method has a Bayesian set up with a Dirichlet Process as prior for the haplotypes in the two parental populations. The basic idea is that only a subset of the complete set of possible haplotypes is present in the population.

Conclusion

Line origin of approximately 95% of the alleles at heterozygous sites was assessed correctly in both simulated and real data. Comparing accuracy of haplotype frequencies inferred with the new algorithm to the accuracy of haplotype frequencies inferred with PHASE, an existing algorithm for haplotype inference, showed that the DP algorithm outperformed PHASE in situations of crossbreeding and that PHASE performed better in situations of random mating.     

Haplotype reconstruction in connected tetraploid F1 populations

In diploid species, many multiparental populations have been developed to increase genetic diversity and quantitative trait loci (QTL) mapping resolution. In these populations, haplotype reconstruction has been used as a standard practice to increase the power of QTL detection in comparison with the marker-based association analysis. However, such software tools for polyploid species are few and limited to a single biparental F1 population. In this study, a statistical framework for haplotype reconstruction has been developed and implemented in the software PolyOrigin for connected tetraploid F1 populations with shared parents, regardless of the number of parents or mating design. Given a genetic or physical map of markers, PolyOrigin first phases parental genotypes, then refines the input marker map, and finally reconstructs offspring haplotypes. PolyOrigin can utilize single nucleotide polymorphism (SNP) data coming from arrays or from sequence-based genotyping; in the latter case, bi-allelic read counts can be used (and are preferred) as input data to minimize the influence of genotype calling errors at low depth. With extensive simulation we show that PolyOrigin is robust to the errors in the input genotypic data and marker map. It works well for various population designs with ≥30 offspring per parent and for sequences with read depth as low as 10x. PolyOrigin was further evaluated using an autotetraploid potato dataset with a 3 × 3 half-diallel mating design. In conclusion, PolyOrigin opens up exciting new possibilities for haplotype analysis in tetraploid breeding populations.     

Haplotype block structure is conserved across mammals

Genetic variation in genomes is organized in haplotype blocks, and species-specific block structure is defined by differential contribution of population history effects in combination with mutation and recombination events. Haplotype maps characterize the common patterns of linkage disequilibrium in populations and have important applications in the design and interpretation of genetic experiments. Although evolutionary processes are known to drive the selection of individual polymorphisms, their effect on haplotype block structure dynamics has not been shown. Here, we present a high-resolution haplotype map for a 5-megabase genomic region in the rat and compare it with the orthologous human and mouse segments. Although the size and fine structure of haplotype blocks are species dependent, there is a significant interspecies overlap in structure and a tendency for blocks to encompass complete genes. Extending these findings to the complete human genome using haplotype map phase I data reveals that linkage disequilibrium values are significantly higher for equally spaced positions in genic regions, including promoters, as compared to intergenic regions, indicating that a selective mechanism exists to maintain combinations of alleles within potentially interacting coding and regulatory regions. Although this characteristic may complicate the identification of causal polymorphisms underlying phenotypic traits, conservation of haplotype structure may be employed for the identification and characterization of functionally important genomic regions.