Pairwise Linkage Disequilibrium Research Articles

Extent of third-order linkage disequilibrium in a composite line of Iberian pigs

BackgroundPrevious studies on linkage disequilibrium have investigated second order linkage disequilibrium in animal and plant populations. The objective of this paper was to investigate the genome-wide levels of third order linkage disequilibrium in a composite line founded by admixture of four Iberian pig strains. A model for the generation of third order linkage disequilibrium by population admixture is proposed. A computer Expectation-Maximization algorithm is developed and applied to the estimation of third order linkage disequilibrium at inter- and intra-chromosomal level using 26,347 SNPs typed in 306 sows. The relationship of third order linkage disequilibrium with physical distance was investigated over 35 million triplets in SSC12. Basic and normalized estimates of inter and intra-chromosomal third order linkage disequilibrium are reported.ResultsGenome-wide analyses revealed that third order linkage disequilibrium is rather common among linked loci in this Iberian pig line. It is shown that population admixture of multiple populations may explain the observed levels of third order linkage disequilibrium although it could be generated by genetic drift. Third order linkage disequilibrium decreases rapidly up to 4 Mb and then declines slowly. The short distances between consecutive markers explain the maintenance of the observed third order linkage disequilibria levels when using a model incorporating the break-up of disequilibrium by recombination. Genome-wide testing also revealed that only 3.6% of the normalized estimates were different from 1, − 1, 0, or from a not well-defined situation in which there is only one possible value for the third order linkage disequilibrium parameter, given allele frequencies and pairwise linkage disequilibria parameters.ConclusionsThird order linkage disequilibrium is common among linked markers in the analyzed pig line and may have been generated by population admixture of multiple populations or by genetic drift. As with second order linkage disequilibrium, the absolute value of the third order linkage disequilibrium decreases with physical distance. Normalization of third order linkage disequilibrium should be avoided for closely linked bi-allelic loci.

Single nucleotide polymorphisms within HLA region are associated with disease relapse for patients with unrelated cord blood transplantation.

Disease relapse occurs in unrelated cord blood transplantation (CBT) even when the alleles of human leukocyte antigen (HLA) are fully matched between donor and recipient. This is similar to that observed in other types of hematopoietic stem cell transplantation. Fourteen single nucleotide polymorphisms (SNPs) within the HLA region have been reported previously by Petersdorf et al. and Piras et al. as transplantation determinants in unrelated hematopoietic cell transplantation. In this study, the genomic sequences within 500 base pairs upstream and downstream of the fourteen transplantation-related SNPs from 53 patients and their HLA-matched unrelated donors were analyzed for determining whether or not genetic variants, conferred by either recipient or donor SNP genotype or by recipient-donor SNP mismatching, were associated with the risk of relapse. Seven SNPs were associated with the risk of relapse in unrelated CBT. These included the donor genotype with the SNPs of rs2523675 and rs2518028 at the telomeric end of HCP5 gene, rs2071479 in the intron of the HLA-DOB gene, and rs2523958 in the MICD gene; and the recipient genotype with SNPs of rs9276982 in the HLA-DOA gene, and rs435766 and rs380924 in the MICD gene. As measured by pair-wise linkage disequilibrium (LD) with D′ as the parameter for normalized standard measurement of LD which compares the observed and expected frequencies of one haplotype comprised by alleles at different loci, rs2523675 had high LD with rs4713466 (D′ = 0.86) and rs2523676 (D′ = 0.91) in the HCP5 gene. The rs2518028 had no LD with all other SNPs except rs2523675 (D′ = 0.76). This study provides the basis for developing a method or algorithm for selecting better unrelated CBT candidate donors.

SNPrune: an efficient algorithm to prune large SNP array and sequence datasets based on high linkage disequilibrium

BackgroundHigh levels of pairwise linkage disequilibrium (LD) in single nucleotide polymorphism (SNP) array or whole-genome sequence data may affect both performance and efficiency of genomic prediction models. Thus, this warrants pruning of genotyping data for high LD. We developed an algorithm, named SNPrune, which enables the rapid detection of any pair of SNPs in complete or high LD throughout the genome.MethodsLD, measured as the squared correlation between phased alleles (r2), can only reach a value of 1 when both loci have the same count of the minor allele. Sorting loci based on the minor allele count, followed by comparison of their alleles, enables rapid detection of loci in complete LD. Detection of loci in high LD can be optimized by computing the range of the minor allele count at another locus for each possible value of the minor allele count that can yield LD values higher than a predefined threshold. This efficiently reduces the number of pairs of loci for which LD needs to be computed, instead of considering all pairwise combinations of loci. The implemented algorithm SNPrune considered bi-allelic loci either using phased alleles or allele counts as input. SNPrune was validated against PLINK on two datasets, using an r2 threshold of 0.99. The first dataset contained 52k SNP genotypes on 3534 pigs and the second dataset contained simulated whole-genome sequence data with 10.8 million SNPs and 2500 animals.ResultsSNPrune removed a similar number of SNPs as PLINK from the pig data but SNPrune was almost 12 times faster than PLINK. From the simulated sequence data with 10.8 million SNPs, SNPrune removed 6.4 and 1.4 million SNPs due to complete and high LD. Results were very similar regardless of whether phased alleles or allele counts were used. Using allele counts and multi-threading with 10 threads, SNPrune completed the analysis in 21 min. Using a sliding window of up to 500,000 SNPs, PLINK removed ~ 43,000 less SNPs (0.6%) in the sequence data and SNPrune was 24 to 170 times faster, using one or ten threads, respectively.ConclusionsThe SNPrune algorithm developed here is able to remove SNPs in high LD throughout the genome very efficiently in large datasets.

Cytokine Gene Variants as Predictors of Type 2 Diabetes Mellitus.

Diabetes is the third widespread after heart disease and cancer. We have investigated genetic polymorphisms in cytokine genes viz. IL-4, IL-1Ra, IL-1β, IL-18, IL-6, TNF-α, IL-10 and ADIPOQ. The aim of study was to investigate the haplotypes, gene-gene interactions and their role in determining individual susceptibility to T2DM of family members with diabetic history. Haplotype analysis of 2 SNPs each in IL-6 and adiponectin genes showing Pairwise Linkage disequilibrium (LD) was done by SHEsis software. Logistic regression was used to study various combinations of gene-gene interactions. The TCGT* set of allele combination appeared to increase the disease risk upto 2 times while TATG* upto 51.4 times when four SNPs are taken together viz. IL-1β-511 C/T, IL-18-607 A/C, ADIPOQ1 +45 G/T and ADIPOQ2 +10211 T/G. Interaction of SNPs in eight genes showed one highly significant combination of alleles, TCGAGCTT* which increased the risk of T2DM upto 7.4 times while CAGAGCGT* allele combination increased the risk upto 4 times. During pedigree analysis in six families with four SNPs, it was interesting to note that susceptible 'AC' genotype of IL-18-607 A/C was frequent in diabetic individuals in almost all families. Moreover, when checked for the presence of risk haplotypes it was observed that TCGT* and TATG* sets of allele combinations were present in most of the diabetic individuals. Individuals with certain abnormal biochemical parameters but not yet diagnosed for T2DM carried the risk genotype or haplotype. This suggested that individuals carrying risk genotypes/haplotypes might be susceptible to T2DM and develop the disease in the future.

Linkage Disequilibrium Estimation in Low Coverage High-Throughput Sequencing Data

High-throughput sequencing methods that multiplex a large number of individuals have provided a cost-effective approach for discovering genome-wide genetic variation in large populations. These sequencing methods are increasingly being utilized in population genetic studies across a diverse range of species. Two side-effects of these methods, however, are (1) sequencing errors and (2) heterozygous genotypes called as homozygous due to only one allele at a particular locus being sequenced, which occurs when the sequencing depth is insufficient. Both of these errors have a profound effect on the estimation of linkage disequilibrium (LD) and, if not taken into account, lead to inaccurate estimates. We developed a new likelihood method, GUS-LD, to estimate pairwise linkage disequilibrium using low coverage sequencing data that accounts for undercalled heterozygous genotypes and sequencing errors. Our findings show that accurate estimates were obtained using GUS-LD, whereas underestimation of LD results if no adjustment is made for the errors.

Relationship between matrilin-1 gene polymorphisms and mandibular retrognathism.

Mandibular retrognathism is a type of malocclusion that refers to an abnormal posterior position of the mandible as a result of a developmental abnormality. From the literature, it is evident that the mandibular growth pattern is determined by the intramembranous ossification of the mandibular body and endochondral ossification of the condyle. Matrilin-1 is a cartilage extracellular matrix protein, and matrilin-1 gene (MATN1) polymorphisms have been found to be involved in dental malocclusions of humans. In this study, we aimed to examine the association between MATN1 polymorphisms and the risk of mandibular retrognathism, in a case-control study with a South Indian population. Eighty-one patients with mandibular retrognathism (SNB,<78°) and 71 controls having an orthognathic mandible (SNB, 80°±2°) were recruited. In both the patient and control groups, subjects with an orthognathic maxilla (SNA, 82°±2°) were included. Three single nucleotide polymorphisms of the MATN1 gene (rs1149048, rs1149042, and rs1065755) were genotyped using polymerase chain reaction-restriction fragment length polymorphism. The statistical association analysis was performed using the chi-square test. Pair-wise linkage disequilibrium was computed, and haplotypes were compared between subjects and controls. Nonparametric tests were used to compare cephalometric measurements between groups. No polymorphic site deviated from Hardy-Weinberg equilibrium in the controls. The rs1149042 genotypes and alleles were found to be associated with reduced risk of mandibular retrognathism. Furthermore, rs1149042 genotypes were associated with mandibular measurements (SNB and ANB). There was no strong and consistent linkage disequilibrium linkage disequilibrium across two different single nucleotide polymorphisms and haplotypes were not associated with mandibular retrognathism. The results of our study suggest an association between the MATN1 gene polymorphisms and mandibular retrognathism.

Population genetics of Enterocytozoon bieneusi in captive giant pandas of China

BackgroundMost studies on Enterocytozoon bieneusi are conducted based on the internal transcribed spacer (ITS) region of the rRNA gene, whereas some have examined E. bieneusi population structures. Currently, the population genetics of this pathogen in giant panda remains unknown. The objective of this study was to determine the E. bieneusi population in captive giant pandas in China.ResultsWe examined 69 E. bieneusi-positive specimens from captive giant pandas in China using five loci (ITS, MS1, MS3, MS4 and MS7) to infer E. bieneusi population genetics. For multilocus genotype (MLG) analysis of E. bieneusi-positive isolates, the MS1, MS3, MS4, and MS7 microsatellite and minisatellite loci were amplified and sequenced in 48, 45, 50 and 47 specimens, respectively, generating ten, eight, nine and five types. We successfully amplified 36 specimens and sequenced all five loci, forming 24 MLGs. Multilocus sequence analysis revealed a strong and significant linkage disequilibrium (LD), indicating a clonal population. This result was further supported by measurements of pairwise intergenic LD and a standardized index of association (ISA) from allelic profile data. The analysis in STRUCTURE suggested three subpopulations in E. bieneusi, further confirmed using right’s fixation index (FST). Subpopulations 1 and 2 exhibited an epidemic structure, whereas subpopulation 3 had a clonal structure.ConclusionsOur results describe E. bieneusi population genetics in giant pandas for the first time, improving the current understanding E. bieneusi epidemiology in the studied region. These data also benefit future studies exploring potential transmission risks from pandas to other animals, including humans.

Assessing the degree of stratification between closely related Holstein-Friesian populations

Genomic information is an important part of the routine evaluation of dairy cattle and provides the wide availability of animals genotyped using single nucleotide polymorphism (SNP) microarrays. We analyzed 2243 Polish and 2294 German Holstein-Friesian bulls genotyped using the Illumina BovineSNP50 BeadChip. For each bull, estimated breeding values (EBVs) calculated from national routine genetic evaluation were available for production traits and for somatic cell score (SCS). Separately for each population, we estimated SNP haplotypes, pairwise linkage disequilibrium (LD), and SNP effects. The SNP genetic covariance between both populations was estimated using a bivariate mixed model. The average LD was lower in the Polish than in the German population and, with increasing genomic distance, LD decays 1.7 times more rapidly in German than in Polish cattle. The comparison of SNP allele frequencies for base populations estimated separately using Polish and German data revealed a very good agreement. The comparison of genetic effects corresponding to various window lengths defined in bp emerged a systematic pattern: regardless of the length of the compared region, few significant differences were found for production traits, while many were observed for SCS. For each trait, the German population had much higher SNP variances than the Polish population and the genetic covariance estimates were all positive. Depending on traits’ inheritance mode, the additive genetic variation can be stored in many genes following the infinitesimal model (like for SCS) or distributed between genes with high effects and the polygenic “background” (like for production traits). Accounting for those differences has implications on the prospective international genomic evaluation.

Vitamin D receptor gene polymorphisms and risk of polycystic ovary syndrome in South Indian women

Polycystic ovary syndrome (PCOS) is the most common endocrine disorder of reproductive age women. Emerging evidence suggests that Vitamin D Receptor (VDR) might be a causal factor for characteristics associated with PCOS such as obesity and type 2 diabetes. Present study investigated association between VDR gene BsmI A/G (rs1544410), ApaI A/C (rs7975232) and TaqI T/C (rs731236) single nucleotide polymorphisms and PCOS risk in South Indian women. Genotyping of VDR gene SNPs was carried out in PCOS patients (n = 95) and controls (n = 130) by PCR-RFLP method and confirmed by sequencing analysis. Haplotype frequencies for multiple loci and the standardized disequilibrium coefficient (D′) for pairwise linkage disequilibrium (LD) were assessed by Haploview software. Results showed significantly increased frequencies of BsmI G/G (p = .0197), ApaI C/C (p = .048), TaqI C/C (p = .044) genotypes and BsmI G (p = .0181), ApaI C (p = .0092), TaqI C (p = .0066) alleles in patients compared to controls. In addition, the frequency of the ‘BsmI G, ApaI C, TaqI C’ haplotype was also significantly elevated in patients (p = .0087). In conclusion, the VDR gene BsmI A/G ApaI A/C TaqI T/C and haplotype may constitute an inheritable risk factor for PCOS in South Indian women.

Genetic analysis of 12 X-STR loci in the Serbian population from Vojvodina Province.

The analysis of 12 X-STR loci included in the Investigator® Argus X-12 kit was performed on a sample of 325 unrelated persons from Vojvodina Province, Republic of Serbia. No significant differences were observed in the allele frequencies in males and females. Heterozygosity values among the studied loci ranged from 67.62 to 94.28%. All loci in female individuals were consistent with the Hardy-Weinberg equilibrium test. The combined power of discrimination values in male and female individuals was 0.9999999994 and 0.999999999999999, respectively. The combined mean exclusion chance was 0.999998 in deficiency cases, 0.9999999977 in normal trio cases, and 0.9999994 in duo cases. Loci DXS10135 and DXS10101 were found to be most polymorphic. The haplotype diversity was found to be greater than 0.993 for all linkage groups. The exact test for pairwise linkage disequilibrium for the 12 loci in the male samples showed significant linkage disequilibrium for the DXS10103-DXS10101 and DXS10134-DXS10146 pairs of loci. The results from the current study confirmed that the panel of 12 X-STR loci is highly polymorphic and informative and can be implemented as a powerful tool in deficient paternity testing and kinship analysis, as well as a useful complement tool of autosomal short tandem repeats (STRs) in forensic investigation. Population differentiation analyses indicated significant differences in genetic structure between the Serbian population and the geographically and ethno-linguistically distant populations, while genetic homogeneity was present in populations with similar geographic origin.

Variance in estimated pairwise genetic distance under high versus low coverage sequencing: The contribution of linkage disequilibrium

The mean pairwise genetic distance among haplotypes is an estimator of the population mutation rate θ and a standard measure of variation in a population. With the advent of next-generation sequencing (NGS) methods, this and other population parameters can be estimated under different modes of sampling. One approach is to sequence individual genomes with high coverage, and to calculate genetic distance over all sample pairs. The second approach, typically used for microbial samples or for tumor cells, is sequencing a large number of pooled genomes with very low individual coverage. With low coverage, pairwise genetic distances are calculated across independently sampled sites rather than across individual genomes. In this study, we show that the variance in genetic distance estimates is reduced with low coverage sampling if the mean pairwise linkage disequilibrium weighted by allele frequencies is positive. Practically, this means that if on average the most frequent alleles over pairs of loci are in positive linkage disequilibrium, low coverage sequencing results in improved estimates of θ, assuming similar per-site read depths. We show that this result holds under the expected distribution of allele frequencies and linkage disequilibria for an infinite sites model at mutation–drift equilibrium. From simulations, we find that the conditions for reduced variance only fail to hold in cases where variant alleles are few and at very low frequency. These results are applied to haplotype frequencies from a lung cancer tumor to compute the weighted linkage disequilibria and the expected error in estimated genetic distance using high versus low coverage.

Joint Analysis of Nuclear and Mitochondrial Variants in Age-Related Macular Degeneration Identifies Novel Loci TRPM1 and ABHD2/RLBP1.

PurposePresently, 52 independent nuclear single nucleotide polymorphisms (nSNPs) have been associated with age-related macular degeneration (AMD) but their effects do not explain all its variance. Genetic interactions between the nuclear and mitochondrial (mt) genome may unearth additional genetic loci previously unassociated with AMD risk.MethodsJoint effects of nSNPs and selected mtSNPs were analyzed by two degree of freedom (2df) joint tests of association in the International AMD Genomics Consortium (IAMDGC) dataset (17,832 controls and 16,144 advanced AMD cases of European ancestry). Subjects were genotyped on the Illumina HumanCoreExome array. After imputation using MINIMAC and the 1000 Genomes Project Phase I reference panel, pairwise linkage disequilibrium pruning, and quality control, 3.9 million nSNPs were analyzed for interaction with mtSNPs chosen based on association in this dataset or publications: A4917G, T5004C, G12771A, and C16069T.ResultsNovel locus TRPM1 was identified with genome-wide significant joint effects (P < 5.0 × 10−8) of two intronic TRPM1 nSNPs and AMD-associated nonsynonymous MT-ND2 mtSNP A4917G. Stratified analysis by mt allele identified an association only in 4917A (major allele) carriers (P = 4.4 × 10−9, odds ratio [OR] = 0.90, 95% confidence interval [CI] = 0.87–0.93). Intronic and intergenic ABHD2/RLBP1 nSNPs demonstrated genome-wide significant joint effects (2df joint test P values from 1.8 × 10−8 to 4.9 × 10−8) and nominally statistically significant interaction effects with MT-ND5 synonymous mtSNP G12771A. Although a positive association was detected in both strata, the association was stronger in 12771A subjects (P = 0.0020, OR = 2.17, 95% CI = 1.34–3.60).ConclusionsThese results show that joint tests of main effects and gene–gene interaction reveal associations at some novel loci that were missed when considering main effects alone.

Association of microRNA-125a and microRNA-499a polymorphisms in chronic periodontitis in a sample south Indian population: A hospital-based genetic association study

Periodontitis is a chronic inflammatory disease, caused by interaction between periodontopathic bacteria and the host immune response. MicroRNAs are small, single-stranded molecules, which play a key role in the regulation of diverse biological processes. Dysregulation of microRNAs function can lead to several diseases such as autoimmune and chronic inflammatory diseases. The objective of the study was to determine the association between selected single nucleotide polymorphisms in miR-125a, miR-499 and LIN28 homology A with chronic periodontitis susceptibility in a sample population from south India. Genotyping of the single nucleotide polymorphisms in miR-125a (rs41275794, rs12976445, rs10404453 and rs12975333), miR-499 (rs3746444) and LIN28 homolog A (rs3811463) was performed in DNA from288 controls (individuals with healthy gingiva) and 262 cases (chronic periodontitis patients) by direct dye-terminator sequencing. Disease association analysis revealed a significant association of the variant alleles of the miR-499a polymorphism (rs3746444) in chronic periodontitis [OR=2.07; 95%CI (1.35–3.17)]. The risk associated C-allele frequency was found to be higher in chronic periodontitis subjects as compared to that of healthy individuals. Similar results were also observed in the dominant model [OR=2.42; 95% CI (1.67–3.51)]. The recessive model for miR-125a polymorphism (rs12976445) was also found to be statistically significant with OR=1.54 and 95% CI (1.03–2.30). The haplotype “GCGGCA” was found to be higher in chronic periodontitis subjects than in healthy individuals. Pairwise linkage disequilibrium analysis exhibited that the polymorphisms, rs41275794 and rs12976445 in miR-125a, were in strong linkage equilibrium (D′=0.97). Epistatic interaction by multifactorial dimensionality reduction analysis revealed that the genotypes of the polymorphisms of miR-125a (rs41275794, rs12976445, rs10404453), miR-499a (rs3746444) and LIN28 (rs3811463) were interacting significantly [OR=2.54 (1.65–3.92)], thereby contributing to the risk of chronic periodontitis.

PPARGC1A and ADIPOQ polymorphisms are associated with aggressive prostate cancer in Mexican-Mestizo men with overweight or obesity.

Obesity constitutes a risk factor for the development of aggressive forms of prostate cancer. It has been proposed, that prostate cancer has a genetic predisposition and that PPARGC1A and ADIPOQ polymorphisms play a role in the development of this condition. To analyse the association of two PPARGC1A and ADIPOQ polymorphisms as well as their haplotypes, with the development of aggressive prostate cancer in Mexican-Mestizo men with overweight or obesity. Two hundred fifty seven men with prostate cancer of Mexican-Mestizo origin were included. Body mass index (BMI) was determined and the degree of prostate cancer aggressiveness by the D'Amico classification. DNA was obtained. Rs7665116 and rs2970870 of PPARGC1A, and rs266729 and rs1501299 of ADIPOQ were studied by real-time PCR allelic discrimination. Pairwise linkage disequilibrium, between single nucleotide polymorphisms was calculated and haplotype analysis was performed. A higher-risk (D'Amico classification) was observed in 21.8% of patients. An association of cancer aggressiveness with rs2970870 of PPARGC1A, and rs501299 of ADIPOQ, as well as with one haplotype of ADIPOQ was documented. This is the first study regarding the relationship of PPARGC1A and ADIPOQ polymorphisms, and the aggressiveness of prostate cancer in men with overweight or obesity.

Abstract 1310: Identification of pleiotropic cancer susceptibility variants from genome-wide association studies reveals functional characteristics

Abstract Background: There exists compelling evidence that some genetic variants are associated with the risk of multiple cancers (i.e., pleiotropy). However, the biological mechanisms of the pleiotropic effects are unclear. Thus, we investigated the functional effects for genetic variants associated with the risk of multiple cancers. Methods: The National Human Genome Research Institute-European Bioinformatics Institute (NHGRI-EBI) GWAS Catalog contains 28,643 variant-trait associations with p &amp;lt; 10-5. We utilized the Experimental Factor Ontology (EFO) to classify cancer traits and obtained all associations between variants and cancer risk. Based on pairwise linkage disequilibrium (LD) determined from the European (EUR) population of Phase 3 of the 1000 Genomes Project in LDlink, correlated variants (R-squared ≥ 0.8) were clustered into groups. Variant groups associated with the risk of multiple cancers were annotated using the Ensembl Variant Effect Predictor (VEP), and tested for functional enrichment using the DAVID Functional Annotation Tool. Results: We identified 1,456 variant-cancer risk associations. The majority (57.1%) of the associations were discovered in European ancestry populations, 19.0% in East Asians, 7.9% in an African Americans or Afro-Caribbeans, and 7.1% in Hispanics or Latin Americans. Removing duplicates, we found 1,034 unique variant-cancer risk associations for 1,005 unique variants and 27 unique cancer sites. After clustering correlated variants, we identified 29 pleiotropic variant groups, of which 2 were associated with risk of five different cancer sites. Variant group rs10936599 and rs12696304 within the MYNN gene (7.8 kb from the TERC gene) was associated with leukemia, multiple myeloma, and colorectal, skin, and bladder cancer. Variant group rs31489, rs31490, rs401681, and rs4975616 within the CLPTM1L gene (22.8 kb from the TERT gene) was associated with leukemia, and lung, pancreatic, skin, and bladder cancer. Among the 42 unique variants that composed the 29 pleiotropic variant groups, 41 variants were within a known gene, and only 1 variant was in an intergenic region. The 41 variants mapped to 26 genes, which were enriched in cellular response to hypoxia (p= 0.0071), establishment of protein localization to telomere (p= 0.0071), and ubl conjugation (p= 0.0072). Conclusion: This study identified and functionally characterized genetic variants showing pleiotropic effects on cancer risk. Our findings improve the understanding of shared biological mechanisms common to different cancers. Clinical implications could include the classification of cancers based on etiology, genetic testing for multiple cancers, and repurposing cancer treatments. Citation Format: Yi-Hsuan Wu, Rebecca E. Graff, Michael N. Passarelli, Thomas J. Hoffmann, Elad Ziv, John S. Witte. Identification of pleiotropic cancer susceptibility variants from genome-wide association studies reveals functional characteristics [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1310. doi:10.1158/1538-7445.AM2017-1310

SNPs, linkage disequilibrium, and chronic mountain sickness in Tibetan Chinese.

Chronic mountain sickness (CMS) is estimated at 1.2% in Tibetans living at the Qinghai–Tibetan Plateau. Eighteen single-nucleotide polymorphisms (SNPs) from nine nuclear genes that have an association with CMS in Tibetans have been analyzed by using pairwise linkage disequilibrium (LD). The SNPs included are the angiotensin-converting enzyme (rs4340), the angiotensinogen (rs699), and the angiotensin II type 1 receptor (AGTR1) (rs5186) from the renin–angiotensin system. A low-density lipoprotein apolipoprotein B (rs693) SNP was also included. From the hypoxia-inducible factor oxygen signaling pathway, the endothetal Per-Arnt-Sim domain protein 1 (EPAS1) and the egl nine homolog 1 (ENGL1) (rs480902) SNPs were included in the study. SNPs from the vascular endothelial growth factor (VEGF) signaling pathway included are the v-akt murine thymoma viral oncogene homolog 3 (rs4590656 and rs2291409), the endothelial cell nitric oxide synthase 3 (rs1007311 and rs1799983), and the (VEGFA) (rs699947, rs34357231, rs79469752, rs13207351, rs28357093, rs1570360, rs2010963, and rs3025039). An increase in LD occurred in 40 pairwise comparisons, whereas a decrease in LD was found in 55 pairwise comparisons between the controls and CMS patients. These changes were found to occur within and between signaling pathways, which suggests that there is an interaction between SNP alleles from different areas of the genome that affect CMS.

Association of ATG16L1 gene haplotype with inflammatory bowel disease in Indians.

Inflammatory bowel disease (IBD) is characterized by multigenic inheritance. Defects in autophagy related genes are considered to show genetic heterogeneity between populations. We evaluated the association of several single nucleotide polymorphisms (SNPs) in the autophagy related 16 like 1 (ATG16L1) gene with IBD in Indians. The ATG16L1 gene was genotyped for ten different SNPs using DNA extracted from peripheral blood of 234 patients with Crohn’s disease (CD), 249 patients with ulcerative colitis (UC) and 393 healthy controls The SNPs rs2241880, rs4663396, rs3792106, rs10210302, rs3792109, rs2241877, rs6737398, rs11682898, rs4663402 and rs4663421 were genotyped using the Sequenom MassArray platform. PLINK was used for the association analysis and pairwise linkage disequilibrium (LD) values. Haplotype analysis was done using Haploview. All SNPs were in Hardy Weinberg equilibrium in cases and controls. The G allele at rs6737398 exhibited a protective association with both CD and UC. The T allele at rs4663402 and C allele at rs4663421 were positively associated with CD and UC. The T allele at rs2241877 exhibited protective association with UC only. The AA genotype at rs4663402 and the GG genotype at rs4663421 were protectively associated with both CD and UC. Haplotype analysis revealed that all the SNPs in tight LD (D’ = 0.76–1.0) and organized in a single haplotype block. Haplotype D was positively associated with IBD (P = 5.8 x 10−6 for CD and 0.002 for UC). SNPs in ATG16L1 were associated with IBD in Indian patients. The relevance to management of individual patients requires further study.

Correction: Genetic Dissection of Drought and Heat Tolerance in Chickpea through Genome-Wide and Candidate Gene-Based Association Mapping Approaches.

[This corrects the article DOI: 10.1371/journal.pone.0096758.].

Relationship of polymorphisms and haplotype in interleukin-16 and adiponectin gene with late-onset Alzheimer's disease risk.

AimsTo investigate the impact of Interleukin-16 (IL- 16) and Adiponectin (ANP) gene single nucleotide polymorphisms (SNPs), gene- gene interactions and haplotype on late-onset Alzheimer’s disease (LOAD) risk.MethodsHardy-Weinberg equilibrium (HWE), haplotype and pairwise linkage disequilibrium (LD) analysis were investigated by using SNPstats (available online at http://bioinfo.iconcologia.net/SNPstats). Generalized multifactor dimensionality reduction (GMDR) was used to examine interaction among 4 SNPs, odds ratio (OR) and 95% confident interval (95%CI) were calculated by logistic regression model.ResultsLOAD risk was significantly higher in carriers of rs266729- G allele than those with CC genotype (CG+ GG versus CC), OR (95%CI) =1.61 (1.26-1.96), and higher in carriers of rs1501299- T allele, OR (95%CI) = 1.62 (1.32-2.12), lower in carriers of rs4072111- T allele, adjusted OR (95%CI) =0.65 (0.44-0.93). We also found a significant gene- gene interaction between rs266729 and rs4072111. Participants with CG or GG of rs266729 and CC of rs4072111 genotype have the highest LOAD risk, OR (95%CI) = 2.62 (1.64 -3.58). Haplotype containing the rs266729- G and rs1501299- T alleles were associated with increased LOAD risk, OR (95%CI)= 1.83 (1.32- 2.43), and haplotype containing the rs1131445- C and rs4072111- T alleles were associated with decreased LOAD risk, OR (95%CI)= 0.53 (0.18- 0.95).ConclusionsWe concluded that rs266729 and rs1501299 minor alleles were associated with increased LOAD risk, but rs4072111 minor allele was associated with decreased LOAD risk. We also found that interaction involving rs266729 and rs4072111, and haplotype combinations were associated with LOAD risk.

Genetic diversity of a New Zealand multi-breed sheep population and composite breeds\u2019 history revealed by a high-density SNP chip

BackgroundKnowledge about the genetic diversity of a population is a crucial parameter for the implementation of successful genomic selection and conservation of genetic resources. The aim of this research was to establish the scientific basis for the implementation of genomic selection in a composite Terminal sheep breeding scheme by providing consolidated linkage disequilibrium (LD) measures across SNP markers, estimating consistency of gametic phase between breed-groups, and assessing genetic diversity measures, such as effective population size (Ne), and population structure parameters, using a large number of animals (n = 14,845) genotyped with a high density SNP chip (606,006 markers). Information generated in this research will be useful for optimizing molecular breeding values predictions and managing the available genetic resources.ResultsOverall, as expected, levels of pairwise LD decreased with increasing distance between SNP pairs. The mean LD r2 between adjacent SNP was 0.26 ± 0.10. The most recent effective population size for all animals (687) and separately per breed-groups: Primera (974), Lamb Supreme (380), Texel (227) and Dual-Purpose (125) was quite variable. The genotyped animals were outbred or had an average low level of inbreeding. Consistency of gametic phase was higher than 0.94 for all breed pairs at the average distance between SNP on the chip (~4.74 kb). Moreover, there was not a clear separation between the breed-groups based on principal component analysis, suggesting that a mixed-breed training population for calculation of molecular breeding values would be beneficial.ConclusionsThis study reports, for the first time, estimates of linkage disequilibrium, genetic diversity and population structure parameters from a genome-wide perspective in New Zealand Terminal Sire composite sheep breeds. The levels of linkage disequilibrium indicate that genomic selection could be implemented with the high density SNP panel. The moderate to high consistency of gametic phase between breed-groups and overlapping population structure support the pooling of the animals in a mixed training population for genomic predictions. In addition, the moderate to high Ne highlights the need to genotype and phenotype a large training population in order to capture most of the haplotype diversity and increase accuracies of genomic predictions. The results reported herein are a first step toward understanding the genomic architecture of a Terminal Sire composite sheep population and for the optimal implementation of genomic selection and genome-wide association studies in this sheep population.