Haplotype Structure Research Articles

Determination and use of haplotypes: Ethnic comparison and association of the lipoprotein lipase gene and coronary artery disease in Mexican-Americans

Purpose: To illustrate an approach of deriving haplotypes for genetic association studies, using the lipoprotein lipase (LPL) gene and coronary artery disease.Methods: Six polymorphisms sufficient to distinguish the most common haplotypes in the 3′ end of LPL were identified by genotyping 10 polymorphisms in a small pilot population. These were used to haplotype LPL in large family samples of Mexican-Americans and non-Hispanic Caucasians. A case-control association study was performed comparing Mexican-Americans with and without coronary artery disease.Results: The two ethnic groups exhibited significant genetic differences. Among Mexican-Americans, homozygosity for LPL haplotype 1 was protective against coronary artery disease (OR = 0.50, 95% CI 0.27–0.91).Conclusion: This study outlines the haplotype structure of the LPL gene, illustrates the utility of haplotype-based analysis in association studies, and demonstrates the importance of defining haplotype frequencies for different ethnic groups.

Further insight into the global variability of the OCA2-HERC2 locus for human pigmentation from multiallelic markers

The OCA2-HERC2 locus is responsible for the greatest proportion of eye color variation in humans. Numerous studies extensively described both functional SNPs and associated patterns of variation over this region. The goal of our study is to examine how these haplotype structures and allelic associations vary when highly variable markers such as microsatellites are used. Eleven microsatellites spanning 357 Kb of OCA2-HERC2 genes are analyzed in 3029 individuals from worldwide populations. We found that several markers display large differences in allele frequency (10% to 35% difference) among Europeans, East Asians and Africans. In Europe, the alleles showing increased frequency can also discriminate individuals with (IrisPlex) predicted blue and brown eyes. Distinct haplotypes are identified around the variants C and T of the functional SNP rs12913832 (associated to blue eyes), with linkage disequilibrium r2 values significant up to 237 Kb. The haplotype carrying the allele rs12913832 C has high frequency (76%) in blue eye predicted individuals (30% in brown eye predicted individuals), while the haplotype associated to the allele rs12913832 T is restricted to brown eye predicted individuals. Finally, homozygosity values reach levels of 91% near rs12913832. Odds ratios show values of 4.2, 7.4 and 10.4 for four markers around rs12913832 and 7.1 for their core haplotype. Hence, this study provides an example on the informativeness of multiallelic markers that, despite their current limited potential contribution to forensic eye color prediction, supports the use of microsatellites for identifying causing variants showing similar genetic features and history.

JAK2 GGCC (46/1) Haplotype in Unprovoked Venous Thrombotic Events

Background: JAK2 GGCC 46/1 haplotype can be represented by four main SNPs (rs3780367, rs10974944, rs12343867, and rs1159782) which replace one cytosine and three thymidines by two guanosines and two cytosines, generating a “GGCC” combination. These four SNPs located on JAK2 introns 10, 12, 14, and 15, respectively, and are always inherited together, being in complete linkage disequilibrium. The 46/1 component of the name came from Jones et al. study where the haplotype structure of the JAK2 gene was mapped using 14 SNPs genotyped by the Wellcome Trust Case Control Consortium (WTCCC) in 1500 healthy blood donors. Two haplotypes (numbers 46 and 1) were found to be identical except for one SNP, and they have a combined frequency of 0.24 in healthy individuals. Numerous observational studies associate this haplotype with myeloproliferative neoplasms (MPNs), as well as splanchnic vein thrombosis (SVT) and non-splanchnic vein thrombosis (non-SVT). In contrast to 24% frequency noted in healthy population, the frequency goes up to 40-80% in JAK2 V617F mutated MPN, and in 64% of those with JAK2 exon 12 mutations (Anelli et al. IJMS, 2018). We herein report our study of JAK2 GGCC (46/1) Haplotype in unprovoked Venous Thrombotic Events (VTE) in patients with negative thrombophilia workup, including negative JAK2 V617F mutation.Methods: We retrospectively identified patients positive for one of the two SNPs (rs12343867 and rs10974900) and unprovoked venous thrombotic among adult patients with negative thrombophilia workup (including JAK2 mutation) treated at tertiary care center from January 2018 to January 2021.Results: We have identified 8 patients, Table (1), that were positive for JAK2 46/1 haplotype SNPs, of whom 62.5% were homozygous 2/2, 25% heterozygous 1/2, while only 12.5% harbor homozygous 1/1 (a normal variant of JAK2 haplotype). The median age 48.5 years (23-65), and the majority (87.5%) were females. Thrombosis site was noted to be SVT in half of the patients, while non-SVT was noted in the other half (12.5% had cerebral vein thrombosis, 12.5% had deep venous thrombosis, 12.5% had a pulmonary embolism, and 12.5% had jugular vein thrombosis). Half of the patients had more than one site venous thrombosis and the other half had only one site. Around 37.50% of the patients had recurrent venous thrombosis on top of therapeutic anticoagulation. Two patients (25%) had high hemoglobin (17.4/16.7) g/dl, but did not fulfill the criteria for polycythemia vera diagnosis (of whom one is a male smoker and one was a female). None of the patients had leukocytosis or thrombocytosis. By imaging, one patient had mild splenomegaly which could be related to SVT.Conclusion: We report on a potential correlation between unprovoked thrombotic events, mainly venous thrombotic events, with JAK2 46/1 haplotype in patients with a negative thrombophilia workup, a finding that merit further investigation. [Display omitted] DisclosuresNo relevant conflicts of interest to declare.

Genetic Polymorphisms of the TGFB1 Signal Peptide and Promoter Region: Role in Wilms Tumor Susceptibility?

The aim of the present study was to investigate the rs1800468 (G-800A), rs1800469 (C-509T), rs1800470 (C29T), and rs1800471 (G74C) TGFB1 genetic polymorphisms and their haplotype structures in patients with Wilms Tumor (WT) and neoplasia-free controls. The genomic DNA was extracted from 35 WT patients and 160 neoplasia-free children, and the TGFB1 polymorphisms were genotyped by polymerase chain reaction, followed by restriction fragment length polymorphism. The haplotype structures were inferred, and permutation and logistic regression tests were performed to check for differences in haplotype distribution between the control and WT individuals. Positive associations were found in the recessive model for rs1800469 T allele (OR: 8.417; 95% CI: 3.177 to 22.297; P < 0.001) and for the rs1800470 C allele (OR: 3.000; 95% CI: 1.296 to 6.944; P = 0.01). Haplotype analysis revealed a significant negative association between GCTG and WT (OR: 0.236, 95% CI: 0.105 to 0.534; P = 0.0002); by contrast, the GTTG haplotype was associated with increased risk for WT (OR: 12.0; 95% CI: 4.202 to 34.270; P < 0.001). Furthermore, rs1800469 was negatively correlated with tumor size and a trend toward a positive correlation for capsular invasion was observed in the dominant model (Tau-b: −0.43, P = 0.02 and tau-b: 0.5, P = 0.06, respectively). This is the first study with rs1800468, rs1800469, rs1800470, and rs1800471 TGFB1 polymorphisms in WT, and our results suggest that the TGFB1 promoter and signal peptide region polymorphisms may be associated with WT susceptibility and clinical presentation.

Bee Trypanosomatids: First Steps in the Analysis of the Genetic Variation and Population Structure of Lotmaria passim, Crithidia bombi and Crithidia mellificae

Trypanosomatids are among the most prevalent parasites in bees but, despite the fact that their impact on the colonies can be quite important and that their infectivity may potentially depend on their genotypes, little is known about the population diversity of these pathogens. Here we cloned and sequenced three non-repetitive single copy loci (DNA topoisomerase II, glyceraldehyde-3-phosphate dehydrogenase and RNA polymerase II large subunit, RPB1) to produce new genetic data from Crithidia bombi, C. mellificae and Lotmaria passim isolated from honeybees and bumblebees. These were analysed by applying population genetic tools in order to quantify and compare their variability within and between species, and to obtain information on their demography and population structure. The general pattern for the three species was that (1) they were subject to the action of purifying selection on nonsynonymous variants, (2) the levels of within species diversity were similar irrespective of the host, (3) there was evidence of recombination among haplotypes and (4) they showed no haplotype structuring according to the host. C. bombi exhibited the lowest levels of synonymous variation (πS= 0.06 ± 0.04 %) — and a mutation frequency distribution compatible with a population expansion after a bottleneck — that contrasted with the extensive polymorphism displayed by C. mellificae (πS= 2.24 ± 1.00 %), which likely has a more ancient origin. L. passim showed intermediate values (πS= 0.40 ± 0.28 %) and an excess of variants a low frequencies probably linked to the spread of this species to new geographical areas.

Haplotype associated RNA expression (HARE) improves prediction of complex traits in maize.

Genomic prediction typically relies on associations between single-site polymorphisms and traits of interest. This representation of genomic variability has been successful for predicting many complex traits. However, it usually cannot capture the combination of alleles in haplotypes and it has generated little insight about the biological function of polymorphisms. Here we present a novel and cost-effective method for imputing cis haplotype associated RNA expression (HARE), studied their transferability across tissues, and evaluated genomic prediction models within and across populations. HARE focuses on tightly linked cis acting causal variants in the immediate vicinity of the gene, while excluding trans effects from diffusion and metabolism. Therefore, HARE estimates were more transferrable across different tissues and populations compared to measured transcript expression. We also showed that HARE estimates captured one-third of the variation in gene expression. HARE estimates were used in genomic prediction models evaluated within and across two diverse maize panels–a diverse association panel (Goodman Association panel) and a large half-sib panel (Nested Association Mapping panel)–for predicting 26 complex traits. HARE resulted in up to 15% higher prediction accuracy than control approaches that preserved haplotype structure, suggesting that HARE carried functional information in addition to information about haplotype structure. The largest increase was observed when the model was trained in the Nested Association Mapping panel and tested in the Goodman Association panel. Additionally, HARE yielded higher within-population prediction accuracy as compared to measured expression values. The accuracy achieved by measured expression was variable across tissues, whereas accuracy by HARE was more stable across tissues. Therefore, imputing RNA expression of genes by haplotype is stable, cost-effective, and transferable across populations.

Trans-ethnic fine-mapping of the MHC region on Parkinson's disease identified shared genetic features associated with the risk

Although different studies reported the genetic predisposition of human leukocyte antigen (HLA) to the risk of Parkinson's disease (PD), the complex haplotype structure and highly polymorphic feature of the major histocompatibility complex (MHC) region has hampered a unified insight. In addition, the evidence for non-European populations was scarce. We conducted trans-ethnic fine-mapping for large cohorts to elucidate shared genetic features of the MHC region associated with the PD risk.

Unraveling the Genomic Architecture of the CYP3A Locus and ADME Genes for Personalized Tacrolimus Dosing.

Tacrolimus (TAC) is an immunosuppressant widely prescribed following an allogenic organ transplant. Due to wide interindividual pharmacokinetic (PK) variability, optimizing TAC dosing based on genetic factors is required to minimize nephrotoxicity and acute rejections. We enrolled 1133 participants receiving TAC from 4 cohorts, consisting of 3 with kidney transplant recipients and 1 with healthy males from clinical trials. The effects of clinical factors were estimated to appropriately control confounding variables. A genome-wide association study, haplotype analysis, and a gene-based association test were conducted using the Korea Biobank Array or targeted sequencing for 114 pharmacogenes. Genome-wide association study verified that CYP3A5*3 is the only common variant associated with TAC PK variability in Koreans. We detected several CYP3A5 and CYP3A4 rare variants that could potentially affect TAC metabolism. The haplotype structure of CYP3A5 stratified by CYP3A5*3 was a significant factor for CYP3A5 rare variant interpretation. CYP3A4 rare variant carriers among CYP3A5 intermediate metabolizers displayed higher TAC trough levels. Gene-based association tests in the 61 absorption, distribution, metabolism, and excretion genes revealed that CYP1A1 are associated with additional TAC PK variability: CYP1A1 rare variant carriers among CYP3A5 poor metabolizers showed lower TAC trough levels than the noncarrier controls. Our study demonstrates that rare variant profiling of CYP3A5 and CYP3A4, combined with the haplotype structures of CYP3A locus, provide additive value for personalized TAC dosing. We also identified a novel association between CYP1A1 rare variants and TAC PK variability in the CYP3A5 nonexpressers that needs to be further investigated.

Haplotype heterogeneity and low linkage disequilibrium reduce reliable prediction of genotypes for the ‑α3.7I form of α-thalassaemia using genome-wide microarray data

Background: The -α 3.7I-thalassaemia deletion is very common throughout Africa because it protects against malaria. When undertaking studies to investigate human genetic adaptations to malaria or other diseases, it is important to account for any confounding effects of α-thalassaemia to rule out spurious associations. Methods: In this study, we have used direct α-thalassaemia genotyping to understand why GWAS data from a large malaria association study in Kilifi Kenya did not identify the α-thalassaemia signal. We then explored the potential use of a number of new approaches to using GWAS data for imputing α-thalassaemia as an alternative to direct genotyping by PCR. Results: We found very low linkage-disequilibrium of the directly typed data with the GWAS SNP markers around α-thalassaemia and across the haemoglobin-alpha ( HBA) gene region, which along with a complex haplotype structure, could explain the lack of an association signal from the GWAS SNP data. Some indirect typing methods gave results that were in broad agreement with those derived from direct genotyping and could identify an association signal, but none were sufficiently accurate to allow correct interpretation compared with direct typing, leading to confusing or erroneous results. Conclusions: We conclude that going forwards, direct typing methods such as PCR will still be required to account for α-thalassaemia in GWAS studies.

Haplotype heterogeneity and low linkage disequilibrium reduce reliable prediction of genotypes for the ‑α 3.7I form of α-thalassaemia using genome-wide microarray data

Background: The -α 3.7I-thalassaemia deletion is very common throughout Africa because it protects against malaria. When undertaking studies to investigate human genetic adaptations to malaria or other diseases, it is important to account for any confounding effects of α-thalassaemia to rule out spurious associations.Methods: In this study, we have used direct α-thalassaemia genotyping to understand why GWAS data from a large malaria association study in Kilifi Kenya did not identify the α-thalassaemia signal. We then explored the potential use of a number of new approaches to using GWAS data for imputing α-thalassaemia as an alternative to direct genotyping by PCR.Results: We found very low linkage-disequilibrium of the directly typed data with the GWAS SNP markers around α-thalassaemia and across the haemoglobin-alpha ( HBA) gene region, which along with a complex haplotype structure, could explain the lack of an association signal from the GWAS SNP data. Some indirect typing methods gave results that were in broad agreement with those derived from direct genotyping and could identify an association signal, but none were sufficiently accurate to allow correct interpretation compared with direct typing, leading to confusing or erroneous results.Conclusions: We conclude that going forwards, direct typing methods such as PCR will still be required to account for α-thalassaemia in GWAS studies.

Transforming growth factor beta 1 (TGFβ1) plasmatic levels and haplotype structures in obesity: a role for TGFβ1 in steatosis development.

Obesity is considered a chronic inflammatory disease and transforming growth factor beta 1 (TGFβ1) might exert important roles in disease pathogenesis regulating adipocyte differentiation and immune-inflammatory environment. However, the role of this cytokine as a biomarker in obesity is poorly addressed. Therefore, the present study aimed to evaluate the impact of TGFB1 polymorphisms and TGFβ1 plasmatic levels in obesity METHODS AND RESULTS: TGFB1 promoter region polymorphisms (rs1800468, G-800A and rs1800469, C-509T) were evaluated in 75 obese patients and 45 eutrophic patients through PCR-RFLP and plasmatic TGFβ1 was quantified through ELISA from 37 of the obese patients, and correlations with clinical and biochemical parameters were tested. Despite no association was found between TGFB1 polymorphisms and obesity susceptibility, several correlations with clinical data were noted. Among others, AC haplotype negatively correlated with plasmatic TGFβ1, while plasmatic TGFβ1 negatively correlated with C-reactive protein and positively correlated with liver abnormalities on ultrasound and, specifically, with steatosis presence and degree. Conversely, GT haplotype, which associates with higher TGFβ1 production, was also positively correlated with the same parameters of liver abnormalities. Further, plasmatic vitamin D negatively correlated with TGFβ1, while positively correlated with AC haplotype. Overall, the results indicate that TGFβ1 might exert important roles in obesity pathophysiology and correlate with biochemical and clinical parameters both at systemic protein as well as at genetic level. Importantly, the consistent positive correlation at both levels with steatosis might suggest this cytokine as a biomarker for this hepatic abnormality in obese patients.

Risk homozygous haplotype regions for autism identifies population-specific ten genes for numerous pathways

BackgroundRecessive homozygous haplotype (rHH) mapping is a reliable tool for identifying recessive genes by detecting homozygous segments of identical haplotype structures. These are shared at a higher frequency amongst probands compared to parental controls. Finding out such rHH blocks in autism subjects can help in deciphering the disorder etiology.ObjectivesThe study aims to detect rHH segments of identical haplotype structure shared at a higher frequency in autism subjects than controls to identify recessive genes responsible for autism manifestation.MethodsIn the present study, 426 unrelated autism genotyped probands with 232 parents (116 trios) were obtained from Gene Expression Omnibus (GEO) Database. Homozygosity mapping analyses have been performed on the samples using standardized algorithms using the Affymetrix GeneChip® 500K SNP Nsp and Sty mapping arrays datasets.ResultsA total of 38 homozygous haplotype blocks were revealed across sample datasets. Upon downstream analysis, 10 autism genes were identified based on selected autism candidate genes criteria. Further, expressive Quantitative Trait Loci (QTL) analysis of SNPs revealed various binding sites for regulatory proteins BX3, FOS, BACH1, MYC, JUND, MAFK, POU2F2, RBBP5, RUNX3, and SMARCA4 impairing essential autism genes CEP290, KITLG, CHD8, and INS2. Pathways and processes such as adherens junction, dipeptidase activity, and platelet-derived growth factor—vital to autism manifestation were identified with varied protein-protein clustered interactions.ConclusionThese findings bring various population clusters with significant rHH genes. It is suggestive of the existence of common but population-specific risk alleles in related autism subjects.

A 24,482-bp deletion is associated with increased seed weight in Brassica napus L.

A major QTL for seed weight was fine-mapped in rapeseed, and a 24,482-bp deletion likely mediates the effect through multiple pathways. Exploration of the genes controlling seed weight is critical to the improvement of crop yield and elucidation of the mechanisms underlying seed formation in rapeseed (Brassica napus L.). We previously identified the quantitative trait locus (QTL) qSW.C9 for the thousand-seed weight (TSW) in a double haploid population constructed from F1 hybrids between the parental accessions HZ396 and Y106. Here, we confirmed the phenotypic effects associated with qSW.C9 in BC3F2 populations and fine-mapped the candidate causal locus to a 266-kb interval. Sequence and expression analyses revealed that a 24,482-bp deletion in HZ396 containing six predicted genes most likely underlies qSW.C9. Differential gene expression analysis and cytological observations suggested that qSW.C9 affects both cell proliferation and cell expansion through multiple signaling pathways. After genotyping of a rapeseed diversity panel to define the haplotype structure, it could be concluded that the selection of germplasm with two specific markers may be effective in improving the seed weight of rapeseed. This study provides a solid foundation for the identification of the causal gene of qSW.C9 and offers a promising target for the breeding of higher-yielding rapeseed.

Natural variation in fecundity is correlated with species-wide levels of divergence in Caenorhabditis elegans.

Life history traits underlie the fitness of organisms and are under strong natural selection. A new mutation that positively impacts a life history trait will likely increase in frequency and become fixed in a population (e.g., a selective sweep). The identification of the beneficial alleles that underlie selective sweeps provides insights into the mechanisms that occurred during the evolution of a species. In the global population of Caenorhabditis elegans, we previously identified selective sweeps that have drastically reduced chromosomal-scale genetic diversity in the species. Here, we measured the fecundity of 121 wild C. elegans strains, including many recently isolated divergent strains from the Hawaiian islands and found that strains with larger swept genomic regions have significantly higher fecundity than strains without evidence of the recent selective sweeps. We used genome-wide association (GWA) mapping to identify three quantitative trait loci (QTL) underlying the fecundity variation. In addition, we mapped previous fecundity data from wild C. elegans strains and C. elegans recombinant inbred advanced intercross lines that were grown in various conditions and detected eight QTL using GWA and linkage mappings. These QTL show the genetic complexity of fecundity across this species. Moreover, the haplotype structure in each GWA QTL region revealed correlations with recent selective sweeps in the C. elegans population. North American and European strains had significantly higher fecundity than most strains from Hawaii, a hypothesized origin of the C. elegans species, suggesting that beneficial alleles that caused increased fecundity could underlie the selective sweeps during the worldwide expansion of C. elegans.

Trans-Ethnic Fine-Mapping of the Major Histocompatibility Complex Region Linked to Parkinson's Disease.

ABSTRACTBackgroundDespite evidence for the role of human leukocyte antigen (HLA) in the genetic predisposition to Parkinson's disease (PD), the complex haplotype structure and highly polymorphic feature of the major histocompatibility complex (MHC) region has hampered a unified insight on the genetic risk of PD. In addition, a majority of the reports focused on Europeans, lacking evidence on other populations.ObjectivesThe aim of this study is to elucidate the genetic features of the MHC region associated with PD risk in trans‐ethnic cohorts.MethodsWe conducted trans‐ethnic fine‐mapping of the MHC region for European populations (14,650 cases and 1,288,625 controls) and East Asian populations (7712 cases and 27,372 controls). We adopted a hybrid fine‐mapping approach including both HLA genotype imputation of genome‐wide association study (GWAS) data and direct imputation of HLA variant risk from the GWAS summary statistics.ResultsThrough trans‐ethnic MHC fine‐mapping, we identified the strongest associations at amino acid position 13 of HLA‐DRβ1 (P = 6.0 × 10−15), which explains the majority of the risk in HLA‐DRB1. In silico prediction revealed that HLA‐DRB1 alleles with histidine at amino acid position 13 (His13) had significantly weaker binding affinity to an α‐synuclein epitope than other alleles (P = 9.6 × 10−4). Stepwise conditional analysis suggested additional independent associations at Ala69 in HLA‐B (P = 1.0 × 10−7). A subanalysis in Europeans suggested additional independent associations at non‐HLA genes in the class III MHC region (EHMT2; P = 2.5 × 10−7).ConclusionsOur study highlights the shared and distinct genetic features of the MHC region in patients with PD across ethnicities. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society

GAMIBHEAR: whole-genome haplotype reconstruction from Genome Architecture Mapping data.

MotivationGenome Architecture Mapping (GAM) was recently introduced as a digestion- and ligation-free method to detect chromatin conformation. Orthogonal to existing approaches based on chromatin conformation capture (3C), GAM’s ability to capture both inter- and intra-chromosomal contacts from low amounts of input data makes it particularly well suited for allele-specific analyses in a clinical setting. Allele-specific analyses are powerful tools to investigate the effects of genetic variants on many cellular phenotypes including chromatin conformation, but require the haplotypes of the individuals under study to be known a priori. So far, however, no algorithm exists for haplotype reconstruction and phasing of genetic variants from GAM data, hindering the allele-specific analysis of chromatin contact points in non-model organisms or individuals with unknown haplotypes.ResultsWe present GAMIBHEAR, a tool for accurate haplotype reconstruction from GAM data. GAMIBHEAR aggregates allelic co-observation frequencies from GAM data and employs a GAM-specific probabilistic model of haplotype capture to optimize phasing accuracy. Using a hybrid mouse embryonic stem cell line with known haplotype structure as a benchmark dataset, we assess correctness and completeness of the reconstructed haplotypes, and demonstrate the power of GAMIBHEAR to infer accurate genome-wide haplotypes from GAM data.Availability and implementationGAMIBHEAR is available as an R package under the open-source GPL-2 license at https://bitbucket.org/schwarzlab/gamibhear.Supplementary information Supplementary data are available at Bioinformatics online.

CaBagE: A Cas9-based Background Elimination strategy for targeted, long-read DNA sequencing.

A substantial fraction of the human genome is difficult to interrogate with short-read DNA sequencing technologies due to paralogy, complex haplotype structures, or tandem repeats. Long-read sequencing technologies, such as Oxford Nanopore's MinION, enable direct measurement of complex loci without introducing many of the biases inherent to short-read methods, though they suffer from relatively lower throughput. This limitation has motivated recent efforts to develop amplification-free strategies to target and enrich loci of interest for subsequent sequencing with long reads. Here, we present CaBagE, a method for target enrichment that is efficient and useful for sequencing large, structurally complex targets. The CaBagE method leverages the stable binding of Cas9 to its DNA target to protect desired fragments from digestion with exonuclease. Enriched DNA fragments are then sequenced with Oxford Nanopore's MinION long-read sequencing technology. Enrichment with CaBagE resulted in a median of 116X coverage (range 39-416) of target loci when tested on five genomic targets ranging from 4-20kb in length using healthy donor DNA. Four cancer gene targets were enriched in a single reaction and multiplexed on a single MinION flow cell. We further demonstrate the utility of CaBagE in two ALS patients with C9orf72 short tandem repeat expansions to produce genotype estimates commensurate with genotypes derived from repeat-primed PCR for each individual. With CaBagE there is a physical enrichment of on-target DNA in a given sample prior to sequencing. This feature allows adaptability across sequencing platforms and potential use as an enrichment strategy for applications beyond sequencing. CaBagE is a rapid enrichment method that can illuminate regions of the 'hidden genome' underlying human disease.

Interleukin-8 genetic diversity, haplotype structure and production differ in two ethnically distinct South African populations

A single nucleotide polymorphism (SNP), 251 bases upstream from the IL-8 transcription start (−251A>T, rs4073), has been extensively investigated in cancers and inflammatory and infectious diseases in predominantly European and Asian populations. We sequenced the IL-8 gene of 109 black and 32 white South African (SA) individuals and conducted detailed characterization of gene variation and haplotype structure. IL-8 production in phytohaemagglutinin (PHA)-stimulated peripheral blood mononuclear cells (PBMCs) of a subset (black: N = 22; white: N = 32) of these individuals was measured using ELISA. Select variants were genotyped for additional black individuals (N = 141), and data from the 1000 Genomes Project were used for haplotype analysis and comparative purposes. In white individuals, the −251A>T SNP formed part of a prevalent six-variant haplotype [haplotype frequency (HF): 61%], Hap-1C, involving the following variants: −251A>T; +394T>G (rs2227307); +780C>T (rs2227306); +1240->A (rs2227541); +1635C>T (rs2227543) and +2770A>T (rs2227543). Hap-1C (−251T+394T+780C+1240+A+1635C+2770A) was composed of two three-variant sub-haplotypes [Hap-1Ca: −251T+394T+1240+A; Hap-1Cb: +780C+1635C+2770A) sharing similarities with haplotypes identified in the black population. Hap-1C was found to be present in European, East and South Asian populations. Four haplotypes were identified in the black population with the two prevalent haplotypes each comprised of two variants: Hap-1B [−251A>T and +1240->A; −251T+1240+A; HF: 14%] and Hap-2B [−743T>C (rs2227532) and +2452A>C (rs2227545); −743C+2452C; HF: 13%]. Populations did not differ in unstimulated PBMC IL-8 production. Upon PHA stimulation, PBMCs from white individuals produced more IL-8 (P = 0.04), suggesting the −251T allele is responsible for higher production, however further analysis revealed that Hap-1C (and constituent sub-haplotypes), did not associate with IL-8 production. Populations did however differ in monocyte number with the white population having significantly more monocytes compared to the black population (P = 0.025), and furthermore monocyte number strongly correlated with IL-8 production in both population groups (black: p = 0.0002, r = 0.71; white: P = 0.0005, r = 0.59). Hap-1B, Hap-2B, and a SNP located one base pair upstream of the IL-8 ATG start codon, +100C>T SNP (rs2227538), all associated with higher IL-8 production in the black population - individuals harbouring at least one of these haplotypes/variant associated with higher IL-8 production (P = 0.003) compared to individuals without. The black population was enriched for individuals harbouring Hap-1B and/or Hap-2B compared to the 1000 Genomes project sub-Saharan African population (P = 0.006), suggesting that SA black individuals may be high IL-8 producers. Given the paucity of IL-8-related studies that have been conducted in populations from sub-Saharan Africa, this study has significantly increased our understanding of this important chemokine in the South African population.

Haplotype structure defines effects of common DPYD variants c.85T > C (rs1801265) and c.496A > G (rs2297595) on dihydropyrimidine dehydrogenase activity: Implication for 5-fluorouracil toxicity.

AimsThe aim of this study was to identify risk variants and haplotypes that impair dihydropyrimidine dehydrogenase (DPD) activity and are, therefore, candidate risk variants for severe toxicity to 5‐fluorouracil (5‐FU) chemotherapy.MethodsPlasma dihydrouracil/uracil (UH2/U) ratios were measured as a population marker for DPD activity in a total of 1382 subjects from 4 independent studies. Genotype and haplotype correlations with UH2/U ratios were assessed.ResultsSignificantly lower UH2/U ratios (panova < 2 × 10−16) were observed in carriers of the 4 well‐studied 5‐FU toxicity risk variants with mean differences (MD) of −43.7% for DPYD c.1905 + 1G > A (rs3918290), −46.0% for DPYD c.1679T > G (rs55886062), −37.1%, for DPYD c.2846A > T (rs67376798), and −13.2% for DPYD c.1129‐5923C > G (rs75017182). An additional variant, DPYD c.496A > G (rs2297595), was also associated with lower UH2/U ratios (P < .0001, MD: −12.6%). A haplotype analysis was performed for variants in linkage disequilibrium with c.496A > G, which consisted of the common variant c.85T > C (rs1801265) and the risk variant c.1129‐5923C > G. Both haplotypes carrying c.496A > G were associated with decreased UH2/U ratios (H3, P = .003, MD: −9.6%; H5, P = .002, MD: −16.9%). A haplotype carrying only the variant c.85T > C (H2) was associated with elevated ratios (P = .004, MD: +8.6%).ConclusionsBased on our data, DPYD‐c.496A > G is a strong candidate risk allele for 5‐FU toxicity. Our data suggest that DPYD‐c.85T > C might be protective; however, the deleterious impacts of the linked alleles c.496A > G and c.1129‐5923C > G likely limit this effect in patients. The possible protective effect of c.85T > C and linkage disequilibrium with c.496A > G and c.1129‐5923C > G may have hampered prior association studies and should be considered in future clinical studies.

Polymorphisms at Myostatin Gene (MSTN) and the Associations with Sport Performances in Anglo-Arabian Racehorses.

Simple SummaryMyostatin is a protein which plays a key role in the regulation and differentiation of cells. The encoding gene, MSTN, and its mutations have been studied in cattle because of the association with increased muscle masses. In horses, some polymorphisms at MSTN have favorable effects on sport performance, as evidenced for the Thoroughbred breed. The Anglo-Arabian horse is a very common breed in many Mediterranean countries because of its extreme adaptability, but genetic investigations are rare in the literature. The present study evidenced the relatively high variability of this breed at MSTN. Many phenotypic effects, such as sex and genealogy, affected the results but, overall, the same polymorphic sites already reported for Thoroughbred horses affected race performances also in Anglo-Arabians.One hundred and eighty Anglo-Arabian horses running 1239 races were sampled for the present study. DNA was extracted from the blood and myostatin gene, MSTN, was genotyped. Moreover, prizes won and places were achieved for the 1239 races to perform association analyses between the different genotypes and sport traits. Two SNPs already reported in previous studies regarding the Thoroughbred breed, rs69472472 and rs397152648, were revealed as polymorphic. The linkage disequilibrium analysis investigating the haplotype structure of MSTN did not evidence any association block. Polymorphism at SNP rs397152648, previously known as g.66493737 T>C, significantly influenced sport traits, with heterozygous horses TC showing better results than homozygotes TT. The portion of variance due to the random effect of the individual animal, and the other phenotypic effects of sex, percentage of Arabian blood and race distance, computed together with the genotype at MSTN in the statistical models, exerted a significant influence. Hence, this information is useful to improve knowledge of the genetic profile of Anglo-Arabian horses and a possible selection for better sport performance.