Extensive Linkage Disequilibrium Research Articles

Low-Frequency and Rare-Coding Variation Contributes to Multiple Sclerosis Risk

(Cell 175, 1679–1687.e1–e7; November 29, 2018) We the editors of Cell have been informed by the Lead Contact and corresponding author of this paper, Dr. Chris Cotsapas, that further analyses following the publication of the paper have indicated that two variants (rs61999302 encoding PRKRA.D33G and rs62176112 encoding PRKRA.P11L) reported in this paper to be associated with MS risk were in fact spurious associations due to the presence of a dispersed duplication event, in which this region of chromosome 2 is duplicated in the MHC region of chromosome 6. The potential for this dispersed duplication to lead to spurious associations has been previously documented (Wellcome Trust Case Control Consortium et al., 2010Craddock N. Hurles M.E. Cardin N. Pearson R.D. Plagnol V. Robson S. Vukcevic D. Barnes C. Conrad D.F. Giannoulatou E. et al.Wellcome Trust Case Control ConsortiumGenome-wide association study of CNVs in 16,000 cases of eight common diseases and 3,000 shared controls.Nature. 2010; 464: 713-720Google Scholar). Upon learning of the potential problems associated with signals in this region, the authors used the data reported in this paper to impute classical HLA alleles and amino acid variants using SNP2HLA (Jia et al., 2013Jia X. Han B. Onengut-Gumuscu S. Chen W.M. Concannon P.J. Rich S.S. Raychaudhuri S. de Bakker P.I.W. Imputing amino acid polymorphisms in human leukocyte antigens.PLoS ONE. 2013; 8: e64683Google Scholar) and found that the two reported variants in PRKRA showed extensive linkage disequilibrium with genotyped and imputed HLA variants, including known MS risk alleles. Conditioning on any of these known variants was sufficient to explain the reported association, leading to the conclusion that these two signals are driven by the duplication and are not true associations of the reported variants with MS risk. With this note, we and the authors would like to alert readers to this issue. The remaining findings in the paper are unchanged, as are the main conclusions of the paper. The authors would also like to thank Dr. Matthew Hurles (Wellcome Sanger Institute) for bringing this issue to their attention. Furthermore, in the printed version of the experimental model and subject details section of the STAR methods, Köln is listed as one of the patient recruitment centers, but Düsseldorf should have been listed instead. One of the authors, Dr. Clemens Warnke, although currently based in Köln, in fact contributed patient samples collected in Düsseldorf to the study. This has been corrected in the article online, and the authors apologize for any confusion these issues may have caused. Low-Frequency and Rare-Coding Variation Contributes to Multiple Sclerosis RiskMitrovič et al.CellOctober 18, 2018In BriefIn a large multi-cohort study, unexplained heritability for multiple sclerosis is detected in low-frequency coding variants that are missed by GWAS analyses, further underscoring the role of immune genes in MS pathology. Full-Text PDF Open AccessLow-Frequency and Rare-Coding Variation Contributes to Multiple Sclerosis RiskMitrovič et al.CellJanuary 23, 2020In Brief(Cell 175, 1679–1687.e1–e7; November 29, 2018) Full-Text PDF Open Access

Eleven Amino Acids of HLA-DRB1 and Fifteen Amino Acids of HLA-DRB3, 4, and 5 Include Potentially Causal Residues Responsible for the Risk of Childhood Type 1 Diabetes.

Next-generation targeted sequencing of HLA-DRB1 and HLA-DRB3, -DRB4, and -DRB5 (abbreviated as DRB345) provides high resolution of functional variant positions to investigate their associations with type 1 diabetes risk and with autoantibodies against insulin (IAA), GAD65 (GADA), IA-2 (IA-2A), and ZnT8 (ZnT8A). To overcome exceptional DR sequence complexity as a result of high polymorphisms and extended linkage disequilibrium among the DR loci, we applied a novel recursive organizer (ROR) to discover disease-associated amino acid residues. ROR distills disease-associated DR sequences and identifies 11 residues of DRB1, sequences of which retain all significant associations observed by DR genes. Furthermore, all 11 residues locate under/adjoining the peptide-binding groove of DRB1, suggesting a plausible functional mechanism through peptide binding. The 15 residues of DRB345, located respectively in the β49–55 homodimerization patch and on the face of the molecule shown to interact with and bind to the accessory molecule CD4, retain their significant disease associations. Further ROR analysis of DR associations with autoantibodies finds that DRB1 residues significantly associated with ZnT8A and DRB345 residues with GADA. The strongest association is between four residues (β14, β25, β71, and β73) and IA-2A, in which the sequence ERKA confers a risk association (odds ratio 2.15, P = 10−18), and another sequence, ERKG, confers a protective association (odds ratio 0.59, P = 10−11), despite a difference of only one amino acid. Because motifs of identified residues capture potentially causal DR associations with type 1 diabetes, this list of residuals is expected to include corresponding causal residues in this study population.

Heterogeneity in the extent of linkage disequilibrium among exonic, intronic, non-coding RNA and intergenic chromosome regions.

Whole-genome sequence data enable construction of high-resolution linkage disequilibrium (LD) maps revealing the LD structure of functional elements within genic and subgenic sequences. The Malecot-Morton model defines LD map distances in linkage disequilibrium units (LDUs), analogous to the centimorgan scale of linkage maps. For whole-genome sequence-derived LD maps, we introduce the ratio of corresponding map lengths kilobases/LDU to describe the extent of LD within genome components. The extent of LD is highly variable across the genome ranging from ~38 kb for intergenic sequences to ~858 kb for centromeric regions. LD is ~16% more extensive in genic, compared with intergenic sequences, reflecting relatively increased selection and/or reduced recombination in genes. The LD profile across 18,268 autosomal genes reveals reduced extent of LD, consistent with elevated recombination, in exonic regions near the 5' end of genes but more extensive LD, compared with intronic sequences, across more centrally located exons. Genes classified as essential and genes linked to Mendelian phenotypes show more extensive LD compared with genes associated with complex traits, perhaps reflecting differences in selective pressure. Significant differences between exonic, intronic and intergenic components demonstrate that fine-scale LD structure provides important insights into genome function, which cannot be revealed by LD analysis of much lower resolution array-based genotyping and conventional linkage maps.

Signatures of selection in Charolais beef cattle identified by genome-wide analysis.

Charolais cattle are one of the most important breeds for meat production worldwide; in México, its selection is mainly made by live weight traits. One strategy for mapping important genomic regions that might influence productive traits is the identification of signatures of selection. This type of genomic features contains loci with extended linkage disequilibrium (LD) and homozygosity patterns that are commonly associated with sites of quantitative trait locus (QTL). Therefore, the objective of this study was to identify the signatures of selection in Charolais cattle genotyped with the GeneSeek Genomic Profiler Bovine HD panel consisting of 77K single nucleotide polymorphisms (SNPs). A total 61,311 SNPs and 819 samples were used for the analysis. Identification of signatures of selection was carried out using the integrated haplotype score (iHS) methodology implemented in the rehh R package. The top ten SNPs with the highest piHS values were located on BTA 4, 5, 6 and 14. By identifying markers in LD with top ten SNPs, the candidate regions defined were mapped to 52.8-59.3Mb on BTA 4; 67.5-69.3 on BTA 5; 39.5-41.0Mb on BTA 6; and 26.4-29.6Mb on BTA 14. The comparison of these candidate regions with the bovine QTLdb effectively confirmed the association (p<0.05) with QTL related to growth traits and other important productive traits. The genomic regions identified in this study indicated selection for growth traits on the Charolais population via the conservation of haplotypes on various chromosomes. These genomic regions and their associated genes could serve as the basis for haplotype association studies and for the identification of causal genes related to growth traits.

Extensive Ethnic Variation and Linkage Disequilibrium at the FCGR2/3 Locus: Different Genetic Associations Revealed in Kawasaki Disease.

The human Fc-gamma receptors (FcγRs) link adaptive and innate immunity by binding immunoglobulin G (IgG). All human low-affinity FcγRs are encoded by the FCGR2/3 locus containing functional single nucleotide polymorphisms (SNPs) and gene copy number variants. This locus is notoriously difficult to genotype and high-throughput methods commonly used focus on only a few SNPs. We performed multiplex ligation-dependent probe amplification for all relevant genetic variations at the FCGR2/3 locus in >4,000 individuals to define linkage disequilibrium (LD) and allele frequencies in different populations. Strong LD and extensive ethnic variation in allele frequencies was found across the locus. LD was strongest for the FCGR2C-ORF haplotype (rs759550223+rs76277413), which leads to expression of FcγRIIc. In Europeans, the FCGR2C-ORF haplotype showed strong LD with, among others, rs201218628 (FCGR2A-Q27W, r2 = 0.63). LD between these two variants was weaker (r2 = 0.17) in Africans, whereas the FCGR2C-ORF haplotype was nearly absent in Asians (minor allele frequency <0.005%). The FCGR2C-ORF haplotype and rs1801274 (FCGR2A-H131R) were in weak LD (r2 = 0.08) in Europeans. We evaluated the importance of ethnic variation and LD in Kawasaki Disease (KD), an acute vasculitis in children with increased incidence in Asians. An association of rs1801274 with KD was previously shown in ethnically diverse genome-wide association studies. Now, we show in 1,028 European KD patients that the FCGR2C-ORF haplotype, although nearly absent in Asians, was more strongly associated with susceptibility to KD than rs1801274 in Europeans. Our data illustrate the importance of interpreting findings of association studies concerning the FCGR2/3 locus with knowledge of LD and ethnic variation.

Genomic Prediction of Growth and Stem Quality Traits in Eucalyptus globulus Labill. at Its Southernmost Distribution Limit in Chile

The present study was undertaken to examine the ability of different genomic selection (GS) models to predict growth traits (diameter at breast height, tree height and wood volume), stem straightness and branching quality of Eucalyptus globulus Labill. trees using a genome-wide Single Nucleotide Polymorphism (SNP) chip (60 K), in one of the southernmost progeny trials of the species, close to its southern distribution limit in Chile. The GS methods examined were Ridge Regression-BLUP (RRBLUP), Bayes-A, Bayes-B, Bayesian least absolute shrinkage and selection operator (BLASSO), principal component regression (PCR), supervised PCR and a variant of the RRBLUP method that involves the previous selection of predictor variables (RRBLUP-B). RRBLUP-B and supervised PCR models presented the greatest predictive ability (PA), followed by the PCR method, for most of the traits studied. The highest PA was obtained for the branching quality (~0.7). For the growth traits, the maximum values of PA varied from 0.43 to 0.54, while for stem straightness, the maximum value of PA reached 0.62 (supervised PCR). The study population presented a more extended linkage disequilibrium (LD) than other populations of E. globulus previously studied. The genome-wide LD decayed rapidly within 0.76 Mbp (threshold value of r2 = 0.1). The average LD on all chromosomes was r2 = 0.09. In addition, the 0.15% of total pairs of linked SNPs were in a complete LD (r2 = 1), and the 3% had an r2 value &gt;0.5. Genomic prediction, which is based on the reduction in dimensionality and variable selection may be a promising method, considering the early growth of the trees and the low-to-moderate values of heritability found in the traits evaluated. These findings provide new understanding of how develop novel breeding strategies for tree improvement of E. globulus at its southernmost range limit in Chile, which could represent new opportunities for forest planting that can benefit the local economy.

Growth factor gene IGF1 is associated with bill size in the black-bellied seedcracker Pyrenestes ostrinus

Pyrenestes finches are unique among birds in showing a non-sex-determined polymorphism in bill size and are considered a textbook example of disruptive selection. Morphs breed randomly with respect to bill size, and differ in diet and feeding performance relative to seed hardness. Previous breeding experiments are consistent with the polymorphism being controlled by a single genetic factor. Here, we use genome-wide pooled sequencing to explore the underlying genetic basis of bill morphology and identify a single candidate region. Targeted resequencing reveals extensive linkage disequilibrium across a 300 Kb region containing the insulin-like growth factor 1 (IGF1) gene, with a single 5-million-year-old haplotype associating with phenotypic dominance of the large-billed morph. We find no genetic similarities controlling bill size in the well-studied Darwin’s finches (Geospiza). Our results show how a single genetic factor may control bill size and provide a foundation for future studies to examine this phenomenon within and among avian species.

Variants in the Chromosomal Region of the Myostatin Gene and Their Association With Lines, Performance, and Body Measurements of Quarter Horses

The objectives of this study were to evaluate differences between the racing and cutting lines of Quarter Horses in the allele frequencies of polymorphisms identified in the chromosomal region of the myostatin gene (MSTN, ECA18) and to analyze the association of polymorphisms in this chromosomal region with the estimated breeding values (EBVs) for racing performance, given by the maximum speed index, height at withers (HW), thoracic perimeter (TP), and body length (BL) in the racing line, and with EBVs for HW, TP, and BL in the cutting line. Genomic, phenotypic, and pedigree data from 420 horses of both sexes, registered in the Brazilian Association of Quarter Horse breeders, were used (352 of racing line and 68 of cutting line). In the genome region studied (MSTN ± 2 Mb), 46 SNPs and one SINE (ERE1) shared among the cutting and racing lines were identified, and the allele frequencies of 32 SNPs and the one SINE differed significantly between the two lines (P < .05). Furthermore, the part of this region closest to the MSTN (±1 Mb) was found to be less polymorphic in the racing line. Regarding phenotype-associated polymorphisms (unadjusted P < .05), those related to HW in the racing line are highlighted because of their large number and extensive linkage disequilibrium. In addition, these markers exhibited significantly higher frequencies of the favorable allele compared to the cutting line.

On the Relationship Between High-Order Linkage Disequilibrium and Epistasis.

A plausible explanation for statistical epistasis revealed in genome wide association analyses is the presence of high order linkage disequilibrium (LD) between the genotyped markers tested for interactions and unobserved functional polymorphisms. Based on findings in experimental data, it has been suggested that high order LD might be a common explanation for statistical epistasis inferred between local polymorphisms in the same genomic region. Here, we empirically evaluate how prevalent high order LD is between local, as well as distal, polymorphisms in the genome. This could provide insights into whether we should account for this when interpreting results from genome wide scans for statistical epistasis. An extensive and strong genome wide high order LD was revealed between pairs of markers on the high density 250k SNP-chip and individual markers revealed by whole genome sequencing in the Arabidopsis thaliana 1001-genomes collection. The high order LD was found to be more prevalent in smaller populations, but present also in samples including several hundred individuals. An empirical example illustrates that high order LD might be an even greater challenge in cases when the genetic architecture is more complex than the common assumption of bi-allelic loci. The example shows how significant statistical epistasis is detected for a pair of markers in high order LD with a complex multi allelic locus. Overall, our study illustrates the importance of considering also other explanations than functional genetic interactions when genome wide statistical epistasis is detected, in particular when the results are obtained in small populations of inbred individuals.

Genetic fine mapping of systemic lupus erythematosus MHC associations in Europeans and African Americans

Genetic variation within the major histocompatibility complex (MHC) contributes substantial risk for systemic lupus erythematosus, but high gene density, extreme polymorphism and extensive linkage disequilibrium (LD) have made fine mapping challenging. To address the problem, we compared two association techniques in two ancestrally diverse populations, African Americans (AAs) and Europeans (EURs). We observed a greater number of Human Leucocyte Antigen (HLA) alleles in AA consistent with the elevated level of recombination in this population. In EUR we observed 50 different A—C—B—DRB1—DQA—DQB multilocus haplotype sequences per hundred individuals; in the AA sample, these multilocus haplotypes were twice as common compared to Europeans. We also observed a strong narrow class II signal in AA as opposed to the long-range LD observed in EUR that includes class I alleles. We performed a Bayesian model choice of the classical HLA alleles and a frequentist analysis that combined both single nucleotide polymorphisms (SNPs) and classical HLA alleles. Both analyses converged on a similar subset of risk HLA alleles: in EUR HLA– B*08:01 + B*18:01 + (DRB1*15:01 frequentist only) + DQA*01:02 + DQB*02:01 + DRB3*02 and in AA HLA–C*17:01 + B*08:01 + DRB1*15:03 + (DQA*01:02 frequentist only) + DQA*02:01 + DQA*05:01+ DQA*05:05 + DQB*03:19 + DQB*02:02. We observed two additional independent SNP associations in both populations: EUR rs146903072 and rs501480; AA rs389883 and rs114118665. The DR2 serotype was best explained by DRB1*15:03 + DQA*01:02 in AA and by DRB1*15:01 + DQA*01:02 in EUR. The DR3 serotype was best explained by DQA*05:01 in AA and by DQB*02:01 in EUR. Despite some differences in underlying HLA allele risk models in EUR and AA, SNP signals across the extended MHC showed remarkable similarity and significant concordance in direction of effect for risk-associated variants.

Effective Genomic Selection in a Narrow-Genepool Crop with Low-Density Markers: Asian Rapeseed as an Example.

Genomic selection (GS) has revolutionized breeding for quantitative traits in plants, offering potential to optimize resource allocation in breeding programs and increase genetic gain per unit of time. Modern high-density single nucleotide polymorphism (SNP) arrays comprising up to several hundred thousand markers provide a user-friendly technology to characterize the genetic constitution of whole populations and for implementing GS in breeding programs. However, GS does not build upon detailed genotype profiling facilitated by maximum marker density. With extensive genome-wide linkage disequilibrium (LD) being a common characteristic of breeding pools, fewer representative markers from available high-density genotyping platforms could be sufficient to capture the association between a genomic region and a phenotypic trait. To examine the effects of reduced marker density on genomic prediction accuracy, we collected data on three traits across 2 yr in a panel of 203 homozygous Chinese semiwinter rapeseed ( L.) inbred lines, broadly encompassing allelic variability in the Asian genepool. We investigated two approaches to selecting subsets of markers: a trait-dependent strategy based on genome-wide association study (GWAS) significance thresholds and a trait-independent method to detect representative tag SNPs. Prediction accuracies were evaluated using cross-validation with ridge-regression best linear unbiased predictions (rrBLUP). With semiwinter rapeseed as a model species, we demonstrate that low-density marker sets comprising a few hundred to a few thousand markers enable high prediction accuracies in breeding populations with strong LD comparable to those achieved with high-density arrays. Our results are valuable for facilitating routine application of cost-efficient GS in breeding programs.

Insights into population genetics and evolution of polyploids and their ancestors.

We have developed the first comprehensive simulator for polyploid genomes (PolySim) and demonstrated its value by performing large-scale simulations to examine the effect of different population parameters on the evolution of polyploids. PolySim is unlimited in terms of ploidy, population size or number of simulated loci. Our process considered the evolution of polyploids from diploid ancestors, polysomic inheritance, inbreeding, recombination rate change in polyploids and gene flow from lower to higher ploidies. We compared the number of segregating single nucleotide polymorphisms, minor allele frequency, heterozygosity, R2 and average kinship relatedness between different simulated scenarios, and to real data from polyploid species. As expected, allotetraploid populations showed no difference from their ancestral diploids when population size remained constant and there was no gene flow or multivalent (MV) pairing between subgenomes. Autotetraploid populations showed significant differences from their ancestors for most parameters and diverged from their ancestral populations faster than allotetraploids. Autotetraploids can have significantly higher heterozygosity, relatedness and extended linkage disequilibrium compared with allotetraploids. Interestingly, autotetraploids were more sensitive to increasing selfing rate and decreasing population size. MV formation can homogenize allotetraploid subgenomes, but this homogenization requires a higher MV rate than previously proposed. Our results can be considered as the first building block to understand polyploid population evolutionary dynamics. PolySim can be used to simulate a wide variety of polyploid organisms that mimic empirical populations, which, in combination with quantitative genetics tools, can be used to investigate the power of genomewide association, genomic selection or breeding programme designs in these species.

Genetic Association in the HLA Region.

The MHC/HLA region has been consistently associated with a large number of complex traits, including but not limited to, most immune-mediated ones. Efforts to pinpoint drivers of this commonly encountered association peak at the short arm of chromosome 6, however, have been challenging, owing to the high density of genes and the long and extended linkage disequilibrium that are characteristic of this region.The development of methods to impute classical HLA alleles and amino acids from SNP genotyping data has offered an important additional layer of information to the investigators seeking to fine map the signal in the region. As a result, imputation-aided association analyses are now typically employed to shed light on the relationship of this locus with disease susceptibility and response to drugs.In this chapter we discuss how the signal in the HLA region can be interrogated in practice, from performing the imputation to understanding its output and to incorporating it into downstream analysis. In addition, we recount some of the analytical approaches that are commonly used and suggest ways in which the findings from such imputation-aided analyses can be interpreted.

Application of Genome Wide Association and Genomic Prediction for Improvement of Cacao Productivity and Resistance to Black and Frosty Pod Diseases.

Chocolate is a highly valued and palatable confectionery product. Chocolate is primarily made from the processed seeds of the tree species Theobroma cacao. Cacao cultivation is highly relevant for small-holder farmers throughout the tropics, yet its productivity remains limited by low yields and widespread pathogens. A panel of 148 improved cacao clones was assembled based on productivity and disease resistance, and phenotypic single-tree replicated clonal evaluation was performed for 8 years. Using high-density markers, the diversity of clones was expressed relative to 10 known ancestral cacao populations, and significant effects of ancestry were observed in productivity and disease resistance. Genome-wide association (GWA) was performed, and six markers were significantly associated with frosty pod disease resistance. In addition, genomic selection was performed, and consistent with the observed extensive linkage disequilibrium, high predictive ability was observed at low marker densities for all traits. Finally, quantitative trait locus mapping and differential expression analysis of two cultivars with contrasting disease phenotypes were performed to identify genes underlying frosty pod disease resistance, identifying a significant quantitative trait locus and 35 differentially expressed genes using two independent differential expression analyses. These results indicate that in breeding populations of heterozygous and recently admixed individuals, mapping approaches can be used for low complexity traits like pod color cacao, or in other species single gene disease resistance, however genomic selection for quantitative traits remains highly effective relative to mapping. Our results can help guide the breeding process for sustainable improved cacao productivity.

Japanese genome-wide association study identifies a significant colorectal cancer susceptibility locus at chromosome 10p14.

Genome‐wide association studies are a powerful tool for searching for disease susceptibility loci. Several studies identifying single nucleotide polymorphisms (SNP) connected intimately to the onset of colorectal cancer (CRC) have been published, but there are few reports of genome‐wide association studies in Japan. To identify genetic variants that modify the risk of CRC oncogenesis, especially in the Japanese population, we performed a multi‐stage genome‐wide association study using a large number of samples: 1846 CRC cases and 2675 controls. We identified 4 SNP (rs7912831, rs4749812, rs7898455 and rs10905453) in chromosome region 10p14 associated with CRC; however, there are no coding or non‐coding genes within this region of fairly extensive linkage disequilibrium (a 500‐kb block) on 10p14. Our study revealed that the 10p14 locus is significantly correlated with susceptibility to CRC in the Japanese population, in accordance with the results of multiple studies in other races.

Population structure, genetic relatedness and linkage disequilibrium blocks in cultivars of tropical soybean (Glycine max)

Soybean (Glycine max L.) is an annual, self-pollinated species, whose genetic base in Brazil is the result of several cycles of selection and effective recombination among a relatively small number of genotypes selected from the USA cultivars. This frequent selection, admixed population, and the crossing of a small number of cultivars can lead to increase the genetic relationship and affect the patterns of population structure. These factors affect the patterns of linkage disequilibrium (LD) blocks, which can be an effective approach for the screening of target loci for agricultural traits in cultivars of tropical soybean. The objectives of this research were to analyze LD blocks, estimate population structure and relatedness through of genotyping of 169 cultivars of tropical soybean by using a BARCSoySNP6K of Illumina iScan platform. The genotyping revealed a high genetic relatedness and population structure among the cultivars of soybean in Brazil, suggesting the existence of a shared genetic base. Our results provide a help to understand the distribution of genetic variation contained within the Brazilian soybean cultivar collection. We showed that the extensive use of a small number of elite genotypes in Brazilian breeding program further reduced the genetic variability, generate extensive LD and probably increase the haplotype block size. These results are in agreement with results of USDA soybean collection, mainly with accessions of north American, when compared with wild and landraces accessions. We constructed a small haplotype block maps (941 blocks), which may be useful for association studies aimed at the identification of genes controlling traits of economic importance in soybean.

HLA haplotypes in primary sclerosing cholangitis patients of admixed and non-European ancestry.

Primary sclerosing cholangitis (PSC) is strongly associated with several human leukocyte antigen (HLA) haplotypes. Due to extensive linkage disequilibrium and multiple polymorphic candidate genes in the HLA complex, identifying the alleles responsible for these associations has proven difficult. We aimed to evaluate whether studying populations of admixed or non-European descent could help in defining the causative HLA alleles. When assessing haplotypes carrying HLA-DRB1*13:01 (hypothesized to specifically increase the susceptibility to chronic cholangitis), we observed that every haplotype in the Scandinavian PSC population carried HLA-DQB1*06:03. In contrast, only 65% of HLA-DRB1*13:01 haplotypes in an admixed/non-European PSC population carried this allele, suggesting that further assessments of the PSC-associated haplotype HLA-DRB1*13:01-DQA1*01:03-DQB1*06:03 in admixed or multi-ethnic populations could aid in identifying the causative allele.

Divergence and evolution of assortative mating in a polygenic trait model of speciation with gene flow.

Assortative mating is an important driver of speciation in populations with gene flow and is predicted to evolve under certain conditions in few-locus models. However, the evolution of assortment is less understood for mating based on quantitative traits, which are often characterized by high genetic variability and extensive linkage disequilibrium between trait loci. We explore this scenario for a two-deme model with migration, by considering a single polygenic trait subject to divergent viability selection across demes, as well as assortative mating and sexual selection within demes, and investigate how trait divergence is shaped by various evolutionary forces. Our analysis reveals the existence of sharp thresholds of assortment strength, at which divergence increases dramatically. We also study the evolution of assortment via invasion of modifiers of mate discrimination and show that the ES assortment strength has an intermediate value under a range of migration-selection parameters, even in diverged populations, due to subtle effects which depend sensitively on the extent of phenotypic variation within these populations. The evolutionary dynamics of the polygenic trait is studied using the hypergeometric and infinitesimal models. We further investigate the sensitivity of our results to the assumptions of the hypergeometric model, using individual-based simulations.

Notch4 and mhc class II polymorphisms are associated with hcv-related benign and malignant lymphoproliferative diseases.

Mixed cryoglobulinemia (MC), is a HCV-related, clinically benign, lymphoproliferative disorder (LPD) that may evolve into a non Hodgkin's lymphoma (NHL). Significant associations were found between two single nucleotide polymorphisms near NOTCH4 (rs2071286) and the HLA class II (rs9461776) genes and HCV-related MC syndrome (MCS). We analyzed NOTCH4 rs2071286 and HLA-II rs9461776 in 3 HCV-related LPD groups (asymptomatic MC, MCS, NHL) with HCV infection without lymphoproliferative disorders.We found a positive relationship between NOTCH4 rs207186 T minor allele frequency (MAF) and patients with HCV-related LPDs at risk of NHL (Chi-square test for trend = 14.84 p = 0.0001), in accordance with an over-dominant model in the NHL group (CT vs CC + TT, OR=1.88, 95% CI 1.24–2.83, p = 0.0026).Regarding HLA II rs9461776, G MAF increased in patients with HCV-related LPDs at risk of NHL (Chi-square test for trend = 8.40 p = 0.0038), in accordance with a recessive genotypic model in the NHL group (G/G vs A/A + A/G, OR = 11.07, 95% CI 2.37–51.64, p = 0.0022).Both NOTCH4 rs2071286 and HLA-II rs9461776 were present on chromosome 6 and showed D’ and r values of linkage disequilibrium (LD) of about 0.5 values, thereby suggesting there is no extensive LD in the HCV+ population.This data shows that the previously demonstrated association between NOTCH4 rs2071286 and HLA-II rs9461776 polymorphisms and HCV-related MCS could be extended to overall patients with HCV-related LPDs. The significant relationship between rs2071286 and rs9461776 MAF and the increased risk for NHL, suggests their use as non-invasive markers to categorize patients at risk of lymphoma.

Intrapopulation genomics in a model mutualist: Population structure and candidate symbiosis genes under selection in Medicago truncatula.

Bottom-up evolutionary approaches, including geographically explicit population genomic analyses, have the power to reveal the mechanistic basis of adaptation. Here, we conduct a population genomic analysis in the model legume, Medicago truncatula, to characterize population genetic structure and identify symbiosis-related genes showing evidence of spatially variable selection. Using RAD-seq, we generated over 26,000 SNPs from 191 accessions from within three regions of the native range in Europe. Results from STRUCTURE analysis identify five distinct genetic clusters with divisions that separate east and west regions in the Mediterranean basin. Much of the genetic variation is maintained within sampling sites, and there is evidence for isolation by distance. Extensive linkage disequilibrium was identified, particularly within populations. We conducted genetic outlier analysis with FST -based genome scans and a Bayesian modeling approach (PCAdapt). There were 70 core outlier loci shared between these distinct methods with one clear candidate symbiosis related gene, DMI1. This work sets that stage for functional experiments to determine the important phenotypes that selection has acted upon and complementary efforts in rhizobium populations.