Regions Of High Linkage Disequilibrium Research Articles

MAPT Locus in Parkinson's Disease Patients of Ashkenazi Origin: A Stratified Analysis.

Introduction: MAPT locus is associated with Parkinson's disease (PD), which is located within a large inversion region of high linkage disequilibrium (LD). We aimed to determine whether the H2-haplotype protective effect and its effect size depends on the GBA1 or LRRK2 risk allele carrier status, and to further characterize genetic alterations that might contribute to its effect. Methods: LD analysis was performed using whole-genome sequencing data of 202 unrelated Ashkenazi Jewish (AJ) PDs. A haplotype-divergent variant was genotyped in a cohort of 1200 consecutively recruited AJ-PDs. The odd ratios were calculated using AJ-non-neuro cases from the gnomAD database as the controls in an un-stratified and a stratified manner according to the mutation carrier status, and the effect on the Age at Motor Symptom Onset (AMSO) was examined. Expression and splicing quantitative trait locus (eQTL and sQTL) analyses were carried out using brain tissues from a database. Results: The H2 haplotype exhibited significant association with PD protection, with a similar effect size in GBA1 carriers, LRRK2-G2019S carriers, and non-carriers (OR = 0.77, 0.69, and 0.82, respectively), and there was no effect on AMSO. The LD interval was narrowed to approximately 1.2 Mb. The H2 haplotype carried potential variants in candidate genes (MAPT and SPPL2C); structural deletions and segmental duplication (KANSL1); and variants affecting gene expression and intron excision ratio in brain tissues (LRRC37A/2). Conclusions: Our results demonstrate that H2 is associated with PD and its protective effect is not influenced by the GBA1/LRRK2 risk allele carrier status. This effect may be genetically complex, resulting from different levels of variations such as missense mutations in relevant genes, structural variations, epigenetic modifications, and RNA expression changes, which may operate independently or in synergy.

A Tagging SNP Set Method Based on Network Community Partition of Linkage Disequilibrium and Node Centrality

Aims:Solving the tagSNP selection problem by network method and reconstructing unknown individual from tagSNPs by a prediction method.Background:As a genetic marker, SNP has been used for linkage analysis of genetic diseases in genome-wide association studies. The genetic information carried by SNPs is redundant in regions of high linkage disequilibrium in the human genome. Therefore, a subset of informative SNPs (tagSNP set) is sufficient to represent the rest of the SNPs, reducing the genotyping cost and computational complexity greatlyMethod:A novel tagSNP set selection method named NCCRT is proposed, which combines the ideas of the network community partition of the SNP network and node centrality ranking to select tagSNPs of genotype data.Methods:A novel tagSNP set selection method named NCCRT is proposed, which combines the ideas of the network community partition of the SNP network and node centrality ranking to select tagSNPs of genotype data.Results:The method is tested on three data sets, including 176 SNPs, 169 SNPs, and 56 SNPs of gene ASAH1, HTR2A, and OLFM4. The experimental results show that our method achieves the best effect in terms of prediction accuracy and stability for ASAH1 and HTR2A.Conclusion:Compared with random sampling, greedy algorithm, and TSMI algorithm, our method does not rely on causal SNP selection, but it can also quickly identify the tagSNP nodes and improve the prediction accuracy.

Characterisation of the Function of a SINE-VNTR-Alu Retrotransposon to Modulate Isoform Expression at the MAPT Locus.

SINE-VNTR-Alu retrotransposons represent one class of transposable elements which contribute to the regulation and evolution of the primate genome and have the potential to be involved in genetic instability and disease progression. However, these polymorphic elements have not been extensively analysed when addressing the missing heritability of neurodegenerative diseases, including Parkinson’s disease (PD) and amyotrophic lateral sclerosis (ALS). SVA_67, a retrotransposon insertion polymorphism, is located in a 1.8 Mb region of high linkage disequilibrium, called the MAPT locus, which is known to contribute to increased risk of developing PD, frontotemporal dementia and other tauopathies. To investigate the role of SVA_67 in directing differential gene expression at this locus, we characterised the impact of SVA_67 allele dosage on isoform expression of several genes in the MAPT locus using the datasets from both the Parkinson’s Progression Markers Initiative and New York Genome Center Consortium Target ALS cohort. The Parkinson’s data was from gene expression in the blood and the ALS data from a variety of CNS regions and allowed us to demonstrate that SVA_67 presence or absence correlated with both isoform- and tissue-specific expression of multiple genes at this locus. This study highlights the importance of addressing SVA polymorphism in disease genetics to gain insight into a better understanding of the role of these regulatory domains to a variety of neurodegenerative diseases.

Genome-wide Association Study Identifies 2 New Loci Associated With Anti-NMDAR Encephalitis.

Background and ObjectivesTo investigate the genetic determinants of the most common type of antibody-mediated autoimmune encephalitis, anti-NMDA receptor (anti-NMDAR) encephalitis.MethodsWe performed a genome-wide association study in 178 patients with anti-NMDAR encephalitis and 590 healthy controls, followed by a colocalization analysis to identify putatively causal genes.ResultsWe identified 2 independent risk loci harboring genome-wide significant variants (p < 5 × 10−8, OR ≥ 2.2), 1 on chromosome 15, harboring only the LRRK1 gene, and 1 on chromosome 11 centered on the ACP2 and NR1H3 genes in a larger region of high linkage disequilibrium. Colocalization signals with expression quantitative trait loci for different brain regions and immune cell types suggested ACP2, NR1H3, MADD, DDB2, and C11orf49 as putatively causal genes. The best candidate genes in each region are LRRK1, encoding leucine-rich repeat kinase 1, a protein involved in B-cell development, and NR1H3 liver X receptor alpha, a transcription factor whose activation inhibits inflammatory processes.DiscussionThis study provides evidence for relevant genetic determinants of antibody-mediated autoimmune encephalitides outside the human leukocyte antigen (HLA) region. The results suggest that future studies with larger sample sizes will successfully identify additional genetic determinants and contribute to the elucidation of the pathomechanism.

Association genetics of the parameters related to nitrogen use efficiency in Brassica juncea L.

Genome wide association studies allowed prediction of 17 candidate genes for association with nitrogen use efficiency. Novel information obtained may provide better understanding of genomic controls underlying germplasm variations for this trait in Indian mustard. Nitrogen use efficiency (NUE) of Indian mustard (Brassica juncea (L.) Czern & Coss.) is low and most breeding efforts to combine NUE with crop performance have not succeeded. Underlying genetics also remain unexplored. We tested 92 SNP-genotyped inbred lines for yield component traits, N uptake efficiency (NUPEFF), nitrogen utilization efficiency (NUTEFF), nitrogen harvest index (NHI) and NUE for two years at two nitrogen doses (No without added N and N100 added @100kg/ha). Genotypes IC-2489-88, M-633, MCP-632, HUJM 1080, GR-325 and DJ-65 recorded high NUE at low N. These also showed improved crop performance under high N. One determinate mustard genotype DJ-113 DT-3 revealed maximum NUTEFF. Genome wide association studies (GWAS) facilitated recognition of 17 quantitative trait loci (QTLs). Environment specificity was high. B-genome chromosomes (B02, B03, B05, B07 and B08) harbored many useful loci. We also used regional association mapping (RAM) to supplement results from GWAS. Annotation of the genomic regions around peak SNPs helped to predict several gene candidates for root architecture, N uptake, assimilation and remobilization. CAT9 (At1g05940) was consistently envisaged for both NUE and NUPEFF. Major N transporter genes, NRT1.8 and NRT3.1 were predicted for explaining variation for NUTEFFand NUPEFF, respectively. Most significant amino acid transporter gene, AAP1 appeared associated with NUE under limited N conditions. All these candidates were predicted in the regions of high linkage disequilibrium. Sequence information of the predicted candidate genes will permit development of molecular markers to aid breeding for high NUE.

Screening non-MAPT genes of the Chr17q21 H1 haplotype in Parkinson's disease

IntroductionThe microtubule-associated protein tau (MAPT) gene is considered a strong genetic risk factor for Parkinson's disease (PD) in Caucasians. MAPT is located within an inversion region of high linkage disequilibrium designated as H1 and H2 haplotype, and contains eight other genes which have been implicated in neurodegeneration. The aim of the current study was to identify common coding variants in strong linkage disequilibrium (LD) within the associated loci on chr17q21 harboring MAPT. MethodsSanger sequencing of coding exons in 90 Caucasian late-onset PD (LOPD) patients was performed. Specific gene sequencing for LRRC37A, LRRC37A2, ARL17A and ARL17B was not possible given the high homology, presence of pseudogenes and copy number variants that are in the region, and therefore four genes (NSF, KANSL1, SPPL2C, and CRHR1) were included in the analysis. Coding variants from these four genes that did not perfectly tag (r2 = 1) the MAPT H1/H2 haplotype were genotyped in an independent replication series of Caucasian PD cases (N = 851) and controls (N = 730). ResultsIn the 90 LOPD cases we identified 30 coding variants. Eleven non-synonymous variants tagged the MAPT H1/H2 haplotype, including two SPPL2C variants (rs12185233 and rs12373123) that had high pathogenic combined annotation dependent depletion (CADD) scores of >20. In the replication series, the non-synonymous KANSL1 rs17585974 variant was in very strong LD with MAPT H1/H2 and had a high CADD score of 24.7. ConclusionWe have identified several non-synonymous variants across neighboring genes of MAPT that may warrant further genetic and functional investigation within the biological etiology of PD.

Increased burden of ultra-rare structural variants localizing to boundaries of topologically associated domains in schizophrenia

Despite considerable progress in schizophrenia genetics, most findings have been for large rare structural variants and common variants in well-imputed regions with few genes implicated from exome sequencing. Whole genome sequencing (WGS) can potentially provide a more complete enumeration of etiological genetic variation apart from the exome and regions of high linkage disequilibrium. We analyze high-coverage WGS data from 1162 Swedish schizophrenia cases and 936 ancestry-matched population controls. Our main objective is to evaluate the contribution to schizophrenia etiology from a variety of genetic variants accessible to WGS but not by previous technologies. Our results suggest that ultra-rare structural variants that affect the boundaries of topologically associated domains (TADs) increase risk for schizophrenia. Alterations in TAD boundaries may lead to dysregulation of gene expression. Future mechanistic studies will be needed to determine the precise functional effects of these variants on biology.

Linkage Disequilibrium and Population Structure in Wild and Cultivated Populations of Rubber Tree (Hevea brasiliensis).

Among rubber tree species, which belong to the Hevea genus of the Euphorbiaceae family, Hevea brasiliensis (Willd. ex Adr.de Juss.) Muell. Arg. is the main commercial source of natural rubber production worldwide. Knowledge of the population structure and linkage disequilibrium (LD) of this species is essential for the efficient organization and exploitation of genetic resources. Here, we obtained single-nucleotide polymorphisms (SNPs) using a genotyping-by-sequencing (GBS) approach and then employed the SNPs for the following objectives: (i) to identify the positions of SNPs on a genetic map of a segregating mapping population, (ii) to evaluate the population structure of a germplasm collection, and (iii) to detect patterns of LD decay among chromosomes for future genetic association studies in rubber tree. A total of 626 genotypes, including both germplasm accessions (368) and individuals from a genetic mapping population (254), were genotyped. A total of 77,660 and 21,283 SNPs were detected by GBS in the germplasm and mapping populations, respectively. The mapping population, which was previously mapped, was constructed with 1,062 markers, among which only 576 SNPs came from GBS, reducing the average interval between two adjacent markers to 4.4 cM. SNPs from GBS genotyping were used for the analysis of genetic structure and LD estimation in the germplasm accessions. Two groups, which largely corresponded to the cultivated and wild populations, were detected using STRUCTURE and via principal coordinate analysis. LD analysis, also using the mapped SNPs, revealed that non-random associations varied along chromosomes, with regions of high LD interspersed with regions of low LD. Considering the length of the genetic map (4,693 cM) and the mean LD (0.49 for cultivated and 0.02 for wild populations), a large number of evenly spaced SNPs would be needed to perform genome-wide association studies in rubber tree, and the wilder the genotypes used, the more difficult the mapping saturation.

The Landscape of Pervasive Horizontal Pleiotropy in Human Genetic Variation is Driven by Extreme Polygenicity of Human Traits and Diseases

Horizontal pleiotropy, where one variant has independent effects on multiple traits, is vitally important for our understanding of the genetic architecture of human phenotypes. We present a method to quantify horizontal pleiotropy using summary statistics from published genomewide association studies (GWASs), and perform validation using simulations. When applied to 1,564 medical phenotypes measured in 337,119 humans from the UK Biobank, our pleiotropy score detects a significant excess of horizontal pleiotropy. This signal of horizontal pleiotropy is pervasive throughout the human genome and across a wide range of phenotypes, but is especially prominent in regions of high linkage disequilibrium and among highly polygenic phenotypes. We identify thousands of loci with extreme horizontal pleiotropy, a majority of which have never been reported in any published GWAS. Our results highlight the central role horizontal pleiotropy plays in the genetic architecture of human phenotypes, and the importance of modeling horizontal pleiotropy in genomic medicine.

Haplotype-specific MAPT exon 3 expression regulated by common intronic polymorphisms associated with Parkinsonian disorders

BackgroundGenome wide association studies have identified microtubule associated protein tau (MAPT) H1 haplotype single nucleotide polymorphisms (SNPs) as leading common risk variants for Parkinson’s disease, progressive supranuclear palsy and corticobasal degeneration. The MAPT risk variants fall within a large 1.8 Mb region of high linkage disequilibrium, making it difficult to discern the functionally important risk variants. Here, we leverage the strong haplotype-specific expression of MAPT exon 3 to investigate the functionality of SNPs that fall within this H1 haplotype region of linkage disequilibrium.MethodsIn this study, we dissect the molecular mechanisms by which haplotype-specific SNPs confer allele-specific effects on the alternative splicing of MAPT exon 3. Firstly, we use haplotype-hybrid whole-locus genomic MAPT vectors studies to identify functional SNPs. Next, we characterise the RNA-protein interactions at two loci by mass spectrometry. Lastly, we knockdown candidate splice factors to determine their effect on MAPT exon 3 using a novel allele-specific qPCR assay.ResultsUsing whole-locus genomic DNA expression vectors to express MAPT haplotype variants, we demonstrate that rs17651213 regulates exon 3 inclusion in a haplotype-specific manner. We further investigated the functionality of this region using RNA-electrophoretic mobility shift assays to show differential RNA-protein complex formation at the H1 and H2 sequence variants of SNP rs17651213 and rs1800547 and subsequently identified candidate trans-acting splicing factors interacting with these functional SNPs sequences by RNA-protein pull-down experiment and mass spectrometry. Finally, gene knockdown of candidate splice factors identified by mass spectrometry demonstrate a role for hnRNP F and hnRNP Q in the haplotype-specific regulation of exon 3 inclusion.ConclusionsWe identified common splice factors hnRNP F and hnRNP Q regulating the haplotype-specific splicing of MAPT exon 3 through intronic variants rs1800547 and rs17651213. This work demonstrates an integrated approach to characterise the functionality of risk variants in large regions of linkage disequilibrium.

Exceptionally long-range haplotypes in Plasmodium falciparum chromosome 6 maintained in an endemic African population.

BackgroundPrevious genome-wide analyses of single nucleotide variation in Plasmodium falciparum identified evidence of an extended haplotype region on chromosome 6 in West Africa, suggesting recent positive selection. Such a pattern is not seen in samples from East Africa or South East Asia, so it could be marking a selective process specific to West Africa. Analyses of the haplotype structure in samples taken at different times could give clues to possible causes of selection.MethodsThis study investigates chromosome 6 extended haplotypes in The Gambia by analysing alleles at multiple microsatellite loci using genome sequence data previously obtained from clinical isolates collected in 2008, followed by genotyping of 13 loci in 439 isolates from 1984, 1991, 2008 and 2014. Temporal changes in haplotype structure and frequencies were determined.ResultsA region of high linkage disequilibrium spanning over 170 kilobases (kb) was identified with both NGS and laboratory determined microsatellite alleles. Multiple long haplotypes were found in all temporal populations from The Gambia. Two of the haplotypes were detected in samples from 1984 and 1991. The frequency of long-range haplotypes increased in 2008 and 2014 populations. There was higher Fst between older and more recent populations at loci in proximity to genes involved in drug metabolism pathways.ConclusionsThe occurrence of several long haplotypes at intermediate frequencies suggests an unusual mode of selection in chromosome 6, possibly combined with recombination suppression on specific haplotypes. Such selection apparently occurred before the emergence of known anti-malarial drug resistance alleles, and could be due to effects of other drugs or unknown processes that have long been operating in this endemic region.Electronic supplementary materialThe online version of this article (doi:10.1186/s12936-016-1560-7) contains supplementary material, which is available to authorized users.

Abstract P6-14-03: Genome wide association study (GWAS) to identify variants conferring ramucirumab-associated hypertension in the ROSE/TRIO-012 breast cancer trial

Abstract Background: In a candidate single nucleotide polymorphism (SNP) association study, we previously identified VEGFR-1 and VEGFR-2 SNPs strongly associated with treatment emergent hypertension (HT) in the ramucirumab (RAM) and docetaxel (Doc) arm in the ROSE/TRIO-012 study, a double-blinded multinational phase III trial that randomized 1,144 patients with advanced breast cancer to receive first-line Doc in combination with RAM or placebo (Mackey et al, JCO Jan 10, 2015:141-148; Mackey et al, JCO, Volume 33, Issue 15_suppl, May 20, 2015: 547). However, candidate SNP studies limit the number of genes for interrogation and a more comprehensive genome wide search may identify critical variants associated with the phenotype of HT. Preliminary analysis indicated that patients experiencing HT with RAM showed better overall survival (Mackey JR, et al (2015). Reply to H. Lee, et al. JCO; in press). These observations provide the potential to identify those patients with genetic variants for predisposition to RAM-associated HT to inform therapeutic decisions. Methods: Genotyping of samples is underway using Affymetrix SNP 6.0 arrays. Genotype data will be filtered for deviations from Hardy Weinberg Equilibrium and minor allele frequency of &amp;gt;0.05. Study subjects (n=792) provided ethics-committee approved prospective consent for this genetic study of whom 478 subjects were allocated RAM + Doc arm. Toxicity grades 0-1 (n= 394 controls; low toxicity) vs. grade &amp;gt;2 (n= 84 cases, high toxicity) is our binary outcome. Dominant genotypic model is assumed. Chi-square test, FDR and/or 10000 permutation tests will be employed (Golden Helix-SVS v8.3) and p&amp;lt;0.05 considered statistically significant. Population stratification will be identified (EIGENSTRAT) and association statistics will be corrected using Eigenvectors along with age as covariates. Fine mapping of loci showing significant associations will be attempted using imputation tools. Results and conclusions: We expect up to 700,000 SNPs to be retained after filtering based on our previous breast cancer GWAS analyses (Damaraju et al. Cancer Research (suppl); Vol 70 (24), page 258s, 2010 and Sehrawat et al Hum Genet. 2011 Oct;130(4):529-37) and 30,000 SNPs to show significance at a nominal p-value (0.05); these will be analysed for regions of high linkage disequilibrium to narrow down potential loci showing association with HT to serve as candidate markers in further independent validation studies. Cumulative dose to adverse events will be considered in the analysis. Identified loci will be interrogated for potential genes in the flanking regions with biological relevance based on pathway analysis. Identified variants from candidate SNP and GWAS may allow developing predictive tools to enable stratification of patients for therapies. The analysis is expected to be completed by mid- November, 2015. Citation Format: Mackey JR, Lipatov O, Martín M, Webster M, Hegg R, Verma S, Ramos-Vázquez M, Fresco R, Thireau F, Houé V, Press MF, Kumaran M, Damaraju S. Genome wide association study (GWAS) to identify variants conferring ramucirumab-associated hypertension in the ROSE/TRIO-012 breast cancer trial. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr P6-14-03.

TNF Block Gene Variants Associate With Pain Intensity in Black Southern Africans With HIV-associated Sensory Neuropathy.

HIV-associated sensory neuropathy (HIV-SN) is a common neurological complication of HIV infection, and it is often painful. Tumor necrosis factor (TNF)-α is implicated in neuropathic pain, but associations between neuropathic pain and polymorphisms in the TNFA gene have not been identified. The "TNF block" is a region of high linkage disequilibrium within the central major histocompatability complex that contains several genes involved in the regulation of inflammation, including TNFA. Polymorphisms in the block have been associated with an altered risk of HIV-SN, but no investigations into whether this region is associated with the painful symptoms of neuropathy have been undertaken. Therefore, we investigated whether polymorphisms in the TNF block are associated with pain intensity in black Southern Africans with HIV-SN. Single-nucleotide polymorphisms (SNPs) defining TNF block haplotypes and African-specific tagSNPs were genotyped in samples from 150 black Southern Africans with HIV-SN. One SNP allele, rs28445017*A, was significantly associated with an increased pain intensity after correction for age, sex, and the CD4 T-cell count. A common 3-SNP haplotype containing rs28445017*G remained associated with a reduced pain intensity after correction for covariates and multiple comparisons. We identified a novel genetic association between polymorphisms in the TNF block and the pain intensity in black Southern Africans with HIV-SN. Our study implicates rs28445017 in painful HIV-SN, although its precise role and whether it may be causative is unclear. rs28445017 was not associated with the risk for HIV-SN as such, highlighting potential differences between the pathophysiology of the neuropathy and the painful features of the neuropathy.

Correlation of the osteoarthritis susceptibility variants that map to chromosome 20q13 with an expression quantitative trait locus operating on NCOA3 and with functional variation at the polymorphism rs116855380.

ObjectiveTo functionally characterize the osteoarthritis (OA) susceptibility variants that map to a region of high linkage disequilibrium (LD) on chromosome 20q13 marked by the single‐nucleotide polymorphism (SNP) rs6094710 and encompassing NCOA3 and SULF2.MethodsNucleic acids were extracted from the cartilage of OA patients. Overall and allelic expression of NCOA3 and SULF2 were measured by quantitative reverse transcription–polymerase chain reaction and pyrosequencing, respectively. The functional effect of SNPs within the 20q13 locus was assessed in vitro using luciferase reporter constructs and electrophoretic mobility shift assays (EMSAs). The in vivo effect of nuclear receptor coactivator 3 (NCOA3) protein depletion on primary human OA articular cartilage chondrocytes was assessed using RNA interference.ResultsExpression of NCOA3 correlated with the genotype at rs6094710 (P = 0.006), and the gene demonstrated allelic expression imbalance (AEI) in individuals heterozygous for the SNP (mean AEI 1.21; P < 0.0001). In both instances, expression of the OA‐associated allele was reduced. In addition, there was reduced enhancer activity of the OA‐associated allele of rs116855380, a SNP in perfect LD with rs6094710 in luciferase assays (P < 0.001). EMSAs demonstrated a protein complex binding with reduced affinity to this allele. Depletion of NCOA3 led to significant changes (all P < 0.05) in the expression of genes involved in cartilage homeostasis.Conclusion NCOA3 is subject to a cis‐acting expression quantitative trait locus in articular cartilage, which correlates with the OA association signal and with the OA‐associated allele of the functional SNP rs116855380, a SNP that is located only 10.3 kb upstream of NCOA3. These findings elucidate the effect of the association of the 20q13 region on OA cartilage and provide compelling evidence of a potentially causal candidate SNP.

A genome-wide association study reveals a quantitative trait locus for days open on chromosome 2 in Japanese Black cattle.

Days open (DO), which is the interval from calving to conception, is an important trait related to reproductive performance in cattle. To identify quantitative trait loci for DO in Japanese Black cattle, we conducted a genome-wide association study with 33,303 single nucleotide polymorphisms (SNPs) using 459 animals with extreme DO values selected from a larger group of 15,488 animals. We identified a SNP on bovine chromosome 2 (BTA2) that was associated with DO. After imputation using phased haplotype data inferred from 586 812 SNPs of 1041 Japanese Black cattle, six SNPs associated with DO were located in an 8.5-kb region of high linkage disequilibrium on BTA2. These SNPs were located on the telomeric side at a distance of 177 kb from the parathyroid hormone 2 receptor (PTH2R) gene. The association was replicated in a sample of 1778 animals. In the replicated population, the frequency of the reduced-DO allele (Q) was 0.63, and it accounted for 1.72% of the total genetic variance. The effect of a Q-to-q allele substitution on DO was a decrease of 3.74 days. The results suggest that the Q allele could serve as a marker in Japanese Black cattle to select animals with superior DO performance.

Leveling the Playing Field in Homozygosity Mapping Using Map Distances.

Studies of linkage disequilibrium (LD) and its variation in the genome are of central importance for understanding evolutionary history, population structure, and selective sweeps. Extreme forms of the latter may result in runs of homozygosity (ROH). In human gene mapping, long ROHs are the basis for homozygosity mapping (HM) with length measured in terms of Mb (106 base pairs physical distance). LD varies greatly over the human genome so that long ROHs tend to occur preferentially in regions of high LD and ROHs of the same length in different regions are not strictly comparable. Thus, in human gene mapping, LD appears as a confounder that needs to be taken into account in the interpretation of ROHs. The effect of varying LD can be mitigated by working on a scale of centimorgans (cM, genetic distance) instead of Mb. We demonstrate this effect for HapMap 3 data on chromosome 19 and show examples with different ROH lengths depending on whether physical or genetic lengths are used. These results suggest that HM should preferably be done on genetic rather than physical distances.

Massively parallel quantification of the regulatory effects of noncoding genetic variation in a human cohort.

We report a novel high-throughput method to empirically quantify individual-specific regulatory element activity at the population scale. The approach combines targeted DNA capture with a high-throughput reporter gene expression assay. As demonstration, we measured the activity of more than 100 putative regulatory elements from 95 individuals in a single experiment. In agreement with previous reports, we found that most genetic variants have weak effects on distal regulatory element activity. Because haplotypes are typically maintained within but not between assayed regulatory elements, the approach can be used to identify causal regulatory haplotypes that likely contribute to human phenotypes. Finally, we demonstrate the utility of the method to functionally fine map causal regulatory variants in regions of high linkage disequilibrium identified by expression quantitative trait loci (eQTL) analyses.

Genome-wide association studies using haplotypes and individual SNPs in Simmental cattle.

Recent advances in high-throughput genotyping technologies have provided the opportunity to map genes using associations between complex traits and markers. Genome-wide association studies (GWAS) based on either a single marker or haplotype have identified genetic variants and underlying genetic mechanisms of quantitative traits. Prompted by the achievements of studies examining economic traits in cattle and to verify the consistency of these two methods using real data, the current study was conducted to construct the haplotype structure in the bovine genome and to detect relevant genes genuinely affecting a carcass trait and a meat quality trait. Using the Illumina BovineHD BeadChip, 942 young bulls with genotyping data were introduced as a reference population to identify the genes in the beef cattle genome significantly associated with foreshank weight and triglyceride levels. In total, 92,553 haplotype blocks were detected in the genome. The regions of high linkage disequilibrium extended up to approximately 200 kb, and the size of haplotype blocks ranged from 22 bp to 199,266 bp. Additionally, the individual SNP analysis and the haplotype-based analysis detected similar regions and common SNPs for these two representative traits. A total of 12 and 7 SNPs in the bovine genome were significantly associated with foreshank weight and triglyceride levels, respectively. By comparison, 4 and 5 haplotype blocks containing the majority of significant SNPs were strongly associated with foreshank weight and triglyceride levels, respectively. In addition, 36 SNPs with high linkage disequilibrium were detected in the GNAQ gene, a potential hotspot that may play a crucial role for regulating carcass trait components.

Inferring mechanisms of copy number change from haplotype structures at the human DEFA1A3 locus.

BackgroundThe determination of structural haplotypes at copy number variable regions can indicate the mechanisms responsible for changes in copy number, as well as explain the relationship between gene copy number and expression. However, obtaining spatial information at regions displaying extensive copy number variation, such as the DEFA1A3 locus, is complex, because of the difficulty in the phasing and assembly of these regions. The DEFA1A3 locus is intriguing in that it falls within a region of high linkage disequilibrium, despite its high variability in copy number (n = 3–16); hence, the mechanisms responsible for changes in copy number at this locus are unclear.ResultsIn this study, a region flanking the DEFA1A3 locus was sequenced across 120 independent haplotypes with European ancestry, identifying five common classes of DEFA1A3 haplotype. Assigning DEFA1A3 class to haplotypes within the 1000 Genomes project highlights a significant difference in DEFA1A3 class frequencies between populations with different ancestry. The features of each DEFA1A3 class, for example, the associated DEFA1A3 copy numbers, were initially assessed in a European cohort (n = 599) and replicated in the 1000 Genomes samples, showing within-class similarity, but between-class and between-population differences in the features of the DEFA1A3 locus. Emulsion haplotype fusion-PCR was used to generate 61 structural haplotypes at the DEFA1A3 locus, showing a high within-class similarity in structure.ConclusionsStructural haplotypes across the DEFA1A3 locus indicate that intra-allelic rearrangement is the predominant mechanism responsible for changes in DEFA1A3 copy number, explaining the conservation of linkage disequilibrium across the locus. The identification of common structural haplotypes at the DEFA1A3 locus could aid studies into how DEFA1A3 copy number influences expression, which is currently unclear.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-614) contains supplementary material, which is available to authorized users.

Allelic expression analysis of the osteoarthritis susceptibility locus that maps to chromosome 3p21 reveals cis-acting eQTLs at GNL3 and SPCS1.

BackgroundAn osteoarthritis (OA) susceptibility locus has been mapped to chromosome 3p21, to a region of high linkage disequilibrium encompassing twelve genes. Six of these genes are expressed in joint tissues and we therefore assessed whether any of the six were subject to cis-acting regulatory polymorphisms active in these tissues and which could therefore account for the association signal.MethodsWe measured allelic expression using pyrosequencing assays that can distinguish mRNA output from each allele of a transcript single nucleotide polymorphism. We assessed RNA extracted from the cartilage and other joint tissues of OA patients who had undergone elective joint replacement surgery. A two-tailed Mann–Whitney exact test was used to test the significance of any allelic differences.ResultsGNL3 and SPCS1 demonstrated significant allelic expression imbalance (AEI) in OA cartilage (GNL3, mean AEI = 1.04, p = 0.0002; SPCS1, mean AEI = 1.07, p < 0.0001). Similar results were observed in other tissues. Expression of the OA-associated allele was lower than that of the non-associated allele for both genes.Conclusionscis-acting regulatory polymorphisms acting on GNL3 and SPCS1 contribute to the OA association signal at chromosome 3p21, and these genes therefore merit further investigation.