Linkage Disequilibrium Calculations Research Articles

Association between sleep-related phenotypes and gut microbiota: a two-sample bidirectional Mendelian randomization study.

An increasing body of evidence suggests a profound interrelation between the microbiome and sleep-related concerns. Nevertheless, current observational studies can merely establish their correlation, leaving causality unexplored. To ascertain whether specific gut microbiota are causally linked to seven sleep-related characteristics and propose potential strategies for insomnia prevention. The study employed an extensive dataset of gut microbiota genetic variations from the MiBioGen alliance, encompassing 18,340 individuals. Taxonomic classification was conducted, identifying 131 genera and 196 bacterial taxa for analysis. Sleep-related phenotype (SRP) data were sourced from the IEU OpenGWAS project, covering traits such as insomnia, chronotype, and snoring. Instrumental variables (IVs) were selected based on specific criteria, including locus-wide significance, linkage disequilibrium calculations, and allele frequency thresholds. Statistical methods were employed to explore causal relationships, including inverse variance weighted (IVW), MR-Egger, weighted median, and weighted Mode. Sensitivity analyses, pleiotropy assessments, and Bonferroni corrections ensured result validity. Reverse causality analysis and adherence to STROBE-MR guidelines were conducted to bolster the study's rigor. Bidirectional Mendelian randomization (MR) analysis reveals a causative interplay between selected gut microbiota and sleep-related phenotypes. Notably, outcomes from the rigorously Bonferroni-corrected examination illuminate profound correlations amid precise compositions of the intestinal microbiome and slumber-associated parameters. Elevated abundance within the taxonomic ranks of class Negativicutes and order Selenomonadales was markedly associated with heightened susceptibility to severe insomnia (OR = 1.03, 95% CI: 1.02-1.05, p = 0.0001). Conversely, the augmented representation of the phylum Lentisphaerae stands in concord with protracted sleep duration (OR = 1.02, 95% CI: 1.01-1.04, p = 0.0005). Furthermore, heightened exposure to the genus Senegalimassilia exhibits the potential to ameliorate the manifestation of snoring symptoms (OR = 0.98, 95% CI: 0.96-0.99, p = 0.0001). This study has unveiled the causal relationship between gut microbiota and SRPs, bestowing significant latent value upon future endeavors in both foundational research and clinical therapy.

GWAS Follow-up Study Discovers a Novel Genetic Signal on 10q21.2 for Atopic Dermatitis in Chinese Han Population.

Atopic dermatitis (AD) is a common inflammatory skin disease with high heritability. Two susceptibility loci have been confirmed in our previous AD genome-wide association study (GWAS). To look for additional genetic factors in Chinese Han ethnicity, we performed a large-scale GWAS follow-up study. Forty-nine top single nucleotide polymorphisms (SNPs) that had never been reported previously were genotyped using Sequenom Massarray system in an independent cohort, which consist of northern Chinese (1634 cases and 1263 controls) and southern Chinese (2985 cases and 9526 controls). Association analyses were performed using PLINK 2 software. Three SNPs in northern and ten SNPs in southern were found exhibiting association evidence with AD (P < 0.05). Finally, SNP rs224108 on 10q21.2 showed high significance for AD in joint analysis of GWAS and replication study (Pmeta = 4.55 × 10−9, OR = 1.21), and was confirmed as an independent genetic marker by Linkage disequilibrium calculation and conditional logistic regression analysis. Bioinformatics analysis strongly suggested that rs224108 may have the potential to alter the target gene expression through non-coding epigenetic regulation effects. Meanwhile, SNP rs11150780 on 17q25.3 was also found suggestive association with AD (Pmeta = 7.64 × 10−7, OR = 1.18). Our findings confirmed a novel susceptibility signal on 10q21.2 for AD in Chinese Han population and advanced the understanding of the genetic contribution to AD.

Genetic relatedness of previously Plant-Variety-Protected commercial maize inbreds.

The emergence of high-throughput, high-density genotyping methods combined with increasingly powerful computing systems has created opportunities to further discover and exploit the genes controlling agronomic performance in elite maize breeding populations. Understanding the genetic basis of population structure in an elite set of materials is an essential step in this genetic discovery process. This paper presents a genotype-based population analysis of all maize inbreds whose Plant Variety Protection certificates had expired as of the end of 2013 (283 inbreds) as well as 66 public founder inbreds. The results provide accurate population structure information and allow for important inferences in context of the historical development of North American elite commercial maize germplasm. Genotypic data was obtained via genotyping-by-sequencing on 349 inbreds. After filtering for missing data, 77,314 high-quality markers remained. The remaining missing data (average per individual was 6.22 percent) was fully imputed at an accuracy of 83 percent. Calculation of linkage disequilibrium revealed that the average r2 of 0.20 occurs at approximately 1.1 Kb. Results of population genetics analyses agree with previously published studies that divide North American maize germplasm into three heterotic groups: Stiff Stalk, Non-Stiff Stalk, and Iodent. Principal component analysis shows that population differentiation is indeed very complex and present at many levels, yet confirms that division into three main sub-groups is optimal for population description. Clustering based on Nei’s genetic distance provides an additional empirical representation of the three main heterotic groups. Overall fixation index (FST), indicating the degree of genetic divergence between the three main heterotic groups, was 0.1361. Understanding the genetic relationships and population differentiation of elite germplasm may help breeders to maintain and potentially increase the rate of genetic gain, resulting in higher overall agronomic performance.

STUDYING THE STRUCTURE OF A GENE POOL POPULATION OF THE RUSSIAN WHITE CHICKEN BREED BY GENOME-WIDE SNP SCAN

A population of the Russian White chickens, bred at the gene pool farm of ARRIFAGB for 25 generations using individual selection, is characterized by resistance to a lowered temperature in the early postnatal period and white colour of the embryonic down. In 2002-2012, breeding was carried out by panmixia, and by now a new population of the Russian White chickens has been formed on the basis of the surviving stock. Comparison of the genetic variability of this population and the archival DNA of representatives of the 2001 population using microarray screening technology will help to assess the population structure and the preservation of the unique characteristics of its genome. The material for the study was DNA extracted from 162 chicken blood samples. Two groups of the Russian White breed were studied, the 2001 population and the current population. Genome-wide analysis using single nucleotide markers (SNP) included screening by means of the Illumina Chicken 60K SNP iSelect BeadChip microarray. Quality control of genotyping, determination of the population genetic structure by multidimensional scaling (MDS), calculation of linkage disequilibrium (LD) and allele frequency in the groups were carried out using PLINK 1.9 software program. The construction of a cluster delimitation model based on SNP genotypes was carried out using the ADMIXTURE program. According to the MDS analysis results, the current population can be divided into four MDS groups, which, when compared to the data of the pedigree, adequately reflect the origin of the studied individuals. The representatives of the ancestral population were genetically similar to the MDS3 group of the current population. Using the F-statistic of the two-way analysis of variance, a significant effect of the group, chromosome, chromosome in the group, and the distance between SNP markers on LD (r2) values was observed. In the 2001 group, the maximum r2 and the high incidence of LD equal to 1 were observed for all chromosomes, with a distance between SNP markers being 500-1000 Kb. There was also the greatest number of monomorphic alleles in this group. Based on the SNP analysis, we may conclude that the current Russian White chicken population is characterized by the disintegration of long LD regions of the ancestral population. Modelling clusters using the ADMIXTURE program revealed differences between the current population groups determined by MDS analysis. The groups composed of individuals included in MDS1 and MDS2 had a homogeneous structure and differed from each other at K = 4 and K = 5. The MDS4 group formed a genetically heterogeneous cluster different from the MDS1 and MDS2 groups at K of 2-5. The MDS3 group was phylogenetically close to the 2001 population (at K of 2-5). In general, the analysis of the current gene pool population of the Russian White chickens showed its heterogeneity while one of its groups (MDS3) was similar to the ancestral population of 2001, which in turn is characterized by a large number of monomorphic alleles and a high frequency of long LD regions. Thus, SNP scanning allowed evaluating the genetic similarity of individuals and the population structure of the Russian White chicken breed. Understanding the genetic structure is an important point in the panmictic breeding and tracking of historical changes in the molecular organization of the genome of a gene pool population with a limited number of animals.

P057 Measuring linkage disequilibrium with next-generation sequencing HLA typing data

Aim To show how HLA typing data from next-generation sequencing (NGS) can dramatically improve our estimates of linkage disequilibrium associations between different loci. Method Given a set of haplotypes in which loci have been placed in phase, for each pair of alleles (Xi, Yi) in two loci X and Y, a 2 × 2 contingency table (a, b, c, d) is built with the following counts: a = CARD(Xi and Yi); b = CARD(X that is not Xi and Yi); c = CARD(Xi and Y that is not Yi), and d = CARD(Xi that is not Xi and Y that is not Yi), where CARD() is a function that returns the cardinality of a set. From these 2 × 2 contingency tables several measures can be calculated, like the classical measure of linkage disequilibrium delta prime, or ‘mutual information’ from information theory, or the ‘ Ochiai semidistance’ from geometry, or chi square from statistical theory, or the likelihood ratio, as well as various useful conditional probabilities. Several testing sets with loci in phase were used. Linkage disequilibrium calculations were made at the classical high-resolution level (4-digit resolution) and at the full resolution level (4-field resolution). NGS typing was done with the Mia Fora technology from Immucor. Results A set of 755 DRB1-DQB1 blocks was used to calculate linkage disequilibrium between Xi = DQA1 ∗ 01:02 and Yi = DRB1 ∗ 13:02 (4-digit res.) and also between Xi = DQA1 ∗ 01:02:01:04 and Yi = DRB1 ∗ 13:02:01 (4-field res.). The results for the first comparison (4-digit res.) were (a, b, c, d) = (70, 0, 196, 489); and for the second comparison (4-field res.) were (a, b, c, d) = (70, 0, 0, 685). While the conditional probability of having DRB1 ∗ 13:02 in the presence of DQA1 ∗ 01:02 (4-digit res.) is 26%, the conditional probability of having DRB1 ∗ 13:02:01 in the presence of DQA1 ∗ 01:02:01:04 (4-field res.) is 100%. The chi square went from 141 for 4-digit res. to 755 for 4-field res. This is just an example of many such findings which will be presented at ASHI 2016. Conclusions Using 4-field resolution data from NGS often brings linkage disequilibrium to a much stronger level, sometime to a total certainty. The most efficient way to present linkage disequilibrium data is, as detailed here, in the form of the cardinality of the logical possibilities as a 2 × 2 contingency table.

A Replication Study Examining Novel Common Single Nucleotide Polymorphisms Identified Through a Prostate Cancer Genome-wide Association Study in a Japanese Population

Five novel prostate cancer risk loci were identified in a recent genome-wide association study (GWAS) of Japanese persons (Takata et al., Nat Genet. 2010;42(9):751-754). Those authors proposed that apart from population-specific linkage disequilibrium patterns, limitations of GWAS single nucleotide polymorphism (SNP) prioritization and/or study design could explain the lack of identification of these loci in GWAS previously conducted among Caucasians. Thus, the authors undertook a replication study in 1,357 prostate cancer patients and 1,403 healthy Australian males of European descent (2004-2008). The rs12653946 SNP at 5p15 was found to be significantly associated with prostate cancer risk (odds ratio = 1.20, 95% confidence interval: 1.07, 1.34; P = 0.002). On the basis of linkage disequilibrium calculations, the rs12653946 SNP represents an independent locus, distinct from the previously identified TERT-CLPTM1L cancer nexus region. Further, analysis from AceView (Thierry-Mieg and Thierry-Mieg, Genome Biol. 2006;7(suppl 1):S12) indicated that rs12653946 falls within the intron of a testis-expressed gene strongly predicted to translate a conceptual 8.1-kilodalton protein named tojy.aApr07. The authors' findings suggest that follow-up of apparently ethnicity-specific risk associations are warranted in order to highlight risk-associated loci for experimental studies and for incorporation into future risk prediction models for prostate cancer.

The distribution of human leukocyte antigen-A, -B, and -DRB1 alleles and haplotypes based on high-resolution genotyping of 167 families from Jiangsu Province, China

We investigated the human leukocyte antigen (HLA)-A, -B, and -DRB1 allele frequencies, the A–B–DRB1, A–B, B–DRB1, and A–DRB1 haplotype frequencies, and the characteristics of linkage disequilibrium between 2 loci in high resolution based on 167 unrelated families from Jiangsu Province, China. A total of 26 alleles at the A locus, 55 alleles at the B locus, and 34 alleles at the DRB1 locus were reported in this study. The top 5 most frequent HLA alleles at the HLA-A, -B, and -DRB1 loci, respectively, were A*11:01, A*24:02, A*02:01, A*33:03, A*30:01; B*13:02, B*40:01 B*46:01, B*58:01, B*54:01; DRB1*09:01, DRB1*07:01, DRB1*12:02, DRB1*15:01, and DRB1*08:03. Several haplotypes with high frequencies were deduced in this study. The top 3 most common A–B–DRB1 haplotypes observed were A*30:01–B*13:02–DRB1*07:01, A*33:03–B*58:01–DRB1*03:01, and A*02:07–B*46:01–DRB1*09:01. The top 3 most common A–B haplotypes were A*30:01–B*13:02, A*33:03–B*58:01, and A*02:07–B*46:01. The top 4 most common A–DRB1 haplotypes were A*30:01–DRB1*07:01, A*33:03–DRB1*13:02, A*24:02–DRB1*09:01, and A*33:03–DRB1*03:01. Finally, the top 3 most common B–DRB1 haplotypes were B*13:02–DRB1*07:01, B*46:01–DRB1*09:01, and B*58:01–DRB1*03:01. From the linkage disequilibrium calculation, the most prominent associations were A*30:01–B*13:02, B*13:02–DRB1*07:01, and A*01:03–DRB1*01:02. These allele and haplotype frequencies could be useful for finding the best matched donors for patients in the China Marrow Donor Program Jiangsu Branch.

The T-1237C Polymorphism of the Toll-like Receptor-9 Gene Is Associated with Chronic Kidney Disease in a Han Chinese Population

Chronic kidney disease (CKD) is increasingly recognized as a global public health problem. As inflammatory processes and genetic factors are involved in the pathogenesis of CKD, we have investigated the potential genetic contribution of Toll-like receptor (TLR) gene polymorphisms in CKD. In a case-control association study, 149 CKD patients and 429 healthy controls were genotyped by real-time polymerase chain reaction. CKD patients were defined as kidney damage (albuminuria, proteinuria or hematuria) or glomerular filtration rate < 60 ml/min/1.73 m(2) for 3 months or more. Single nucleotide polymorphisms (SNPs) at TLR-2 G2408A, TLR-4 A12874G and C13174T, and TLR-9 T-1237C, T-1486C, and G1635A were assessed, and linkage disequilibrium calculations and haplotype association analysis were undertaken. The functions of TLR-9 have been documented to recognize the viral and bacterial CpG DNA sequences, whereas detects microbe-derived peptidoglycan and lipopeptides and TLR-4 binds lipopolysaccharides. SNPs within the TLR genes may influence promoter activity, mRNA conformation and subcellular localization, and/or protein structure and function. Our results show that only the TLR-9 T-1237C and G1635A gene polymorphisms demonstrate an association with CKD (p = 0.002 and p = 0.04, respectively). The TLR-9 TCA haplotype at T-1237C, T-1486C, and G1635A was associated with a lower risk of CKD, whereas the TTA haplotype was associated with a higher risk of CKD. In the Han Chinese population, those who carry the C and A alleles at SNPs T-1237C and G1635A in the TLR-9 gene appear to be more susceptible to the development of CKD.

Three microsatellites from the T1DGC MHC data set show highly significant association with type 1 diabetes, independent of the HLA‐DRB1, ‐DQA1 and ‐DQB1 genes

The aim of this study was to test the microsatellites in the Type 1 Diabetes Genetics Consortium major histocompatibility complex (MHC) data set for association with type 1 diabetes (T1D) independent of the HLA-DRB1, -DQA1 and -DQB1 genes. The data set was edited to contain only one affected child per family, and broad ethnic subgroups were defined. Genotypes for HLA-DRB1, -DQA1 and -DQB1 were replaced by a haplotype code spanning all three loci, with phase inferred based on common haplotypes. The final data set contained 8190 samples in 2301 families, 59 microsatellites and the DRB1-DQA1-DQB1 haplotype code. Statistical analyses consisted of conditional logistic regression and haplotype estimations and linkage disequilibrium calculations. The data set was screened using a main effects test approach adjusted for DRB1-DQA1-DQB1, and significant results tested for validity. After these procedures, four markers remained significant at the Bonferroni-corrected threshold: D6S2773 (p = 0.00014), DG6S185 (p = 0.00015), DG6S398 (p = 0.00043) and D6S2998 (p = 0.00015). These results were supported by allelic tests conditioned on DRB1-DQA1-DQB1 haplotypes, except for DG6S185, which may contain artefacts. We have identified three microsatellites that mark additional risk factors for T1D at highly significant levels in the MHC. Further analyses are needed to establish the relationship with other possible genetic determinants in this region.

Population stratification and patterns of linkage disequilibrium

Although the importance of selecting cases and controls from the same population has been recognized for decades, the recent advent of genome-wide association studies has heightened awareness of this issue. Because these studies typically deal with large samples, small differences in allele frequencies between cases and controls can easily reach statistical significance. When, unbeknownst to a researcher, cases and controls have different substructures, the number of false-positive findings is inflated. There have been three recent developments of purely statistical approaches to assessing the ancestral comparability of case and control samples: genomic control, structured association, and multivariate reduction analyses. The widespread use of high-throughput technology has allowed the quick and accurate genotyping of the large number of markers required by these methods. Group 13 dealt with four population stratification issues: single-nucleotide polymorphism marker selection, association testing, nonstandard methods, and linkage disequilibrium calculations in stratified or mixed ethnicity samples. We demonstrated that there are continuous axes of ethnic variation in both data sets of Genetic Analysis Workshop 16. Furthermore, ignoring this structure created P-value inflation for a variety of phenotypes. Principal-components analysis (or multidimensional scaling) can control inflation as covariates in a logistic regression. One can weigh for local ancestry estimation and allow the use of related individuals. Problems arise in the presence of extremely high association or unusually strong linkage disequilibrium (e.g., in chromosomal inversions). Our group also reported a method for performing an association test controlling for substructure, when genome-wide markers are not available, to explicitly compute stratification.

Polymorphisms near SOCS3 are associated with obesity and glucose homeostasis traits in Hispanic Americans from the Insulin Resistance Atherosclerosis Family Study

The SOCS3 gene product participates in the feedback inhibition of a range of cytokine signals. Most notably, SOCS3 inhibits the functioning of leptin and downstream steps in insulin signaling after being expressed by terminal transcription factors, such as STAT3 and c-fos. The SOCS3 gene is located in the chromosome region 17q24-17q25, previously linked to body mass index (BMI), visceral adipose tissue (VAT), and waist circumference (WAIST) in Hispanic families in the Insulin Resistance Atherosclerosis Family Study (IRASFS). A high density map of 1,536 single nucleotide polymorphisms (SNPs) was constructed to cover a portion of the 17q linkage interval in 1,425 Hispanic subjects from 90 extended families in IRASFS. Analysis of this dense SNP map data revealed evidence of association of rs9914220 (located 10 kb 5' of the SOCS3 gene) with BMI, VAT, and WAIST (P-value ranging from 0.003 to 0.017). Using a tagging SNP approach, rs9914220 and 22 additional SOCS3 SNPs were genotyped for genetic association analysis with measures of adiposity and glucose homeostasis. The adiposity phenotypes utilized in association analyses included BMI, WAIST, waist to hip ratio (WHR), subcutaneous adipose tissue, VAT, and visceral to subcutaneous ratio (VSR). Linkage disequilibrium calculations revealed three haplotype blocks near SOCS3. Haplotype Block 3 (5' of SOCS3) contained SNPs consistently associated with BMI, WAIST, WHR, and VAT (P-values ranging from 2.00 x 10(-4) to 0.036). Haplotype Block 1 contained single-SNPs that were associated with most adiposity traits except for VSR (P-values ranging from 0.002 to 0.047). When trait associated SNPs were included in linkage analyses as covariates, a reduction of VAT LOD score from 1.26 to 0.76 above the SOCS3 locus (110 cM) was observed. Multi-SNP haplotype testing using the quantitative pedigree disequilibrium test was broadly consistent with the single-SNP associations. In conclusion, these results support a role for SOCS3 genetic variants in human obesity.

Cubic exact solutions for the estimation of pairwise haplotype frequencies: implications for linkage disequilibrium analyses and a web tool 'CubeX'.

BackgroundThe frequency of a haplotype comprising one allele at each of two loci can be expressed as a cubic equation (the 'Hill equation'), the solution of which gives that frequency. Most haplotype and linkage disequilibrium analysis programs use iteration-based algorithms which substitute an estimate of haplotype frequency into the equation, producing a new estimate which is repeatedly fed back into the equation until the values converge to a maximum likelihood estimate (expectation-maximisation).ResultsWe present a program, "CubeX", which calculates the biologically possible exact solution(s) and provides estimated haplotype frequencies, D', r2 and χ2 values for each. CubeX provides a "complete" analysis of haplotype frequencies and linkage disequilibrium for a pair of biallelic markers under situations where sampling variation and genotyping errors distort sample Hardy-Weinberg equilibrium, potentially causing more than one biologically possible solution. We also present an analysis of simulations and real data using the algebraically exact solution, which indicates that under perfect sample Hardy-Weinberg equilibrium there is only one biologically possible solution, but that under other conditions there may be more.ConclusionOur analyses demonstrate that lower allele frequencies, lower sample numbers, population stratification and a possible |D'| value of 1 are particularly susceptible to distortion of sample Hardy-Weinberg equilibrium, which has significant implications for calculation of linkage disequilibrium in small sample sizes (eg HapMap) and rarer alleles (eg paucimorphisms, q < 0.05) that may have particular disease relevance and require improved approaches for meaningful evaluation.

Caffeic acid O-methyltransferase allelic polymorphism characterization and analysis in different maize inbred lines

The caffeic acid O-methyltransferase (COMT) gene is a gene candidate forimproving silage maize digestibility. Before assessing possible functionalconsequences of COMT gene polymorphism, the nucleotide sequence variability atthe COMT locus in 6 normal maize inbred lines was characterized and analyzed.Inaddition, the COMT bm3-3 mutation characterization was specified. Even if theCOMT gene seemed to evolve according to the predictions of the neutral model,estimates of nucleotide diversity indicated that polymorphism distributionwithin the sequence was not uniform and that nonsynonymous sites were the leastvariable. Linkage disequilibrium calculation and estimates of recombinationshowed that polymorphism was a dimorphic segregating distribution caused by alot of recombination events.

Alcohol dehydrogenase genes: restriction fragment length polymorphisms for ADH4 (pi-ADH) and ADH5 (chi-ADH) and construction of haplotypes among different ADH classes.

Of the five human alcohol dehydrogenase (ADH) genes located in the region q21-25 of chromosome 4, genetic markers have been reported previously only for class I enzymes, ADH1-3. Here, new restriction fragment length polymorphisms (RFLPs) are described for the genes of two other classes, ADH4 (pi) and ADH5 (chi or formaldehyde dehydrogenase, FDH). The frequencies and modes of inheritance of these RFLPs were determined with DNA both from unrelated individuals and from families. A polymorphic PstI site is assigned to the fourth intron of the ADH4 gene. Pairwise linkage disequilibrium calculations for these new RFLPs and already known RFLPs at the ADH2 and ADH3 loci establish strong linkage disequilibria between polymorphic MspI and BstXI sites in the ADH5 gene as well as between XbaI and MspI sites in the ADH3 gene. Furthermore, linkage disequilibria were detected between RFLPs of the ADH2 and ADH3 genes as well as between those of the ADH4 and ADH5 genes. The latter disequilibrium implies a hitherto unknown physical proximity of two genes belonging to different ADH classes. The RFLPs were used to construct chromosomal haplotypes that include three ADH classes. Of the 16 possible haplotypes for four RFLP markers used here, 10 were experimentally detected. The potential application of the ADH RFLPs and haplotypes in linkage or association studies of inherited diseases such as familial "alcoholism" is discussed.