Disease Association Studies Research Articles

Simulation study of misclassification bias in association studies employing partial-mouth protocols.

To simulate the exposure misclassification bias potential in studies of perio-systemic disease associations due to the use of partial-mouth recording (PMR) protocols. Using data from 640 participants in the Dental Longitudinal Study, we evaluated distributions of clinical periodontitis parameters to simulate hypothetical outcome probabilities using bootstrap sampling. Logistic regression models were fit using the hypothetical outcome as the dependent variable. Models were run for exposure classifications based on full-mouth recording (FMR) and PMR protocols over 10,000 repetitions. The impact of periodontitis exposure misclassification was dependent on periodontitis severity. Per cent relative bias for simulated ORs of size 1.5, 2 and 4 ranged from 0% to 30% for the effect of severe periodontitis. The magnitude and direction of the bias was dependent on the underlying distribution of the clinical parameters used in the simulation and the size of the association being estimated. Simulated effects of moderate periodontitis were consistently biased towards the null. Exposure misclassification bias occurring through the use of PMR protocols may be dependent on the sensitivity of the classification system applied. Using the CDC-AAP case definition, bias in the estimated effects of severe disease was small, on average. Whereas effects of moderate disease were underestimated to a larger degree.

Genome-wide association studies of inflammatory bowel disease in German shepherd dogs.

Canine Inflammatory Bowel Disease (IBD) is considered a multifactorial disease caused by complex interactions between the intestinal immune system, intestinal microbiota and environmental factors in genetically susceptible individuals. Although IBD can affect any breed, German shepherd dogs (GSD) in the UK are at increased risk of developing the disease. Based on previous evidence, the aim of the present study was to identify single nucleotide polymorphisms (SNPs), which may confer genetic susceptibility or resistance to IBD using a genome-wide association study (GWAS). Genomic DNA was extracted from EDTA blood or saliva samples of 96 cases and 98 controls. Genotyping of cases and controls was performed on the Canine Illumina HD SNP array and data generated was analyzed using PLINK. Several SNPs and regions on chromosomes 7,9,11 and 13 were detected to be associated with IBD using different SNP-by-SNP association methods and FST windows approach. Searching one Mb up-and down-stream of the most significant SNPs, as identified by single SNP analysis as well as 200Kb before and after the start and the end position of the associated regions identified by FST windows approach, we identified 63 genes. Using a combination of pathways analysis and a list of genes that have been reported to be involved in human IBD, we identified 16 candidate genes potentially associated with IBD in GSD.

Allelic and haplotype diversity of HLA-A, HLA-B and HLA-DRB1 gene at high resolution in the Nanning Han population.

The distribution of human leucocyte antigen (HLA) allele and haplotype varied among different ethnic populations. In this study, we investigated the allele and haplotype frequencies of HLA-A, HLA-B and HLA-DRB1 loci in the Nanning Han population who live in Guangxi province of China. We identified 26 HLA-A, 56 HLA-B and 31 HLA-DRB1 alleles in 562 Nanning individuals of Han ethnic group by sequence-based typing method. Of these, the three most common alleles in HLA-A, HLA-B and HLA-DRB1 loci, respectively, were A*11:01 (32.12%), A*02:07 (12.54%), A*24:02 (12.01%); B*46:01 (14.41%), B*15:02 (13.61%), B*40:01 (11.48%); DRB1*15:01 (14.15%), DRB1*16:02 (11.57%) and DRB1*12:02 (10.14%). With the exception of HLA-DRB1, the p values of the HLA-A and HLA-B loci showed that the HLA allelic distribution in this population was in accordance with Hardy-Weinberg expectation (p>0.05). A total of 173 HLA~A-B~DRB1 haplotype with a frequency of >0.1% were presented and the three most common haplotype were HLA-A*33:03~B*58:01~DRB1*03:01 (6.12%), HLA-A*11:01~B*15:02~DRB1*12:02 (3.39%) and HLA-A*11:01~B*15:02~DRB1*15:01 (3.22%). The phylogenetic tree and the principal component analysis suggested that Nanning Han population had a relative close genetic relationship with Chinese Zhuang population and a relative distant genetic relationship with Northern Han Chinese. The information will be useful for anthropological studies, for HLA matching in transplantation and disease association studies in the Chinese population.

Exploring morphological variation in tomato (Solanum lycopersicum): A combined study of disease resistance, genetic divergence and association of characters

Abstract Genetic variation of quantitative traits is a prevalent characteristic among cultivated tomato varieties. Twenty tomato genotypes comprising indigenous varieties and commercial cultivars, cultured in the Western Region of Cameroon were evaluated using fourteen quantitative traits for disease resistance, phenotypic divergence and heritability estimates. The experiment was carried out using a randomized completed blocks design with three replications. Data collections were disease characteristics, plant development features and yield attributes. The analysis of variance revealed significant variation among genotypes for all the experimental quantitative traits. Hybrid varieties had significantly more fruit yield (1066.00 g/plant), single fruit weight (57.28 g), fruit diameter (4.47 cm) and pericarp thickness (0.54 cm) compared to standard and indigenous varieties. These indigenous varieties were significantly more resistant to late blight, alternaria leaf spot and viral diseases. They also had significantly higher collar diameter (16.30 mm), number of primary branches per plant (8.45), number of fruit per plant (31.58) and plant height (88.33 cm) compared to standard and hybrid plants. The genotype local 2 was the third most productive (1576.39 g / plant) after Rio Semagri (1984.80 g/plant) and Sakato F1 (1691.69 g/plant). Heritability and genetic advance estimates were high for twelve of the fourteen studied quantitative traits. Fruit yield showed significant positive correlations with single fruit weight and number of fruit per plant. However, significant negative correlation was found between fruit yield and time to 50 % flowering, 50 % fruiting, 50 % maturity and viral disease. The first three and the first four components in the principal component analysis explained, respectively, 77.85 % and 88.38 % of the total variation observed among genotypes. The first component determined 41.42 % of the total variation, dominated by the collar diameter, the number of primary branches per plant and plant height. This study clearly indicated that indigenous varieties are the most disease resistant genotypes and are having substantial fruit yield (945.30 g/plant) similar to standard varieties and at a touching distance to hybrid cultivars.

Exploring morphological variation in tomato (Solanum lycopersicum): A combined study of disease resistance, genetic divergence and association of characters

Abstract Genetic variation of quantitative traits is a prevalent characteristic among cultivated tomato varieties. Twenty tomato genotypes comprising indigenous varieties and commercial cultivars, cultured in the Western Region of Cameroon were evaluated using fourteen quantitative traits for disease resistance, phenotypic divergence and heritability estimates. The experiment was carried out using a randomized completed blocks design with three replications. Data collections were disease characteristics, plant development features and yield attributes. The analysis of variance revealed significant variation among genotypes for all the experimental quantitative traits. Hybrid varieties had significantly more fruit yield (1066.00 g/plant), single fruit weight (57.28 g), fruit diameter (4.47 cm) and pericarp thickness (0.54 cm) compared to standard and indigenous varieties. These indigenous varieties were significantly more resistant to late blight, alternaria leaf spot and viral diseases. They also had significantly higher collar diameter (16.30 mm), number of primary branches per plant (8.45), number of fruit per plant (31.58) and plant height (88.33 cm) compared to standard and hybrid plants. The genotype local 2 was the third most productive (1576.39 g / plant) after Rio Semagri (1984.80 g/plant) and Sakato F1 (1691.69 g/plant). Heritability and genetic advance estimates were high for twelve of the fourteen studied quantitative traits. Fruit yield showed significant positive correlations with single fruit weight and number of fruit per plant. However, significant negative correlation was found between fruit yield and time to 50 % flowering, 50 % fruiting, 50 % maturity and viral disease. The first three and the first four components in the principal component analysis explained, respectively, 77.85 % and 88.38 % of the total variation observed among genotypes. The first component determined 41.42 % of the total variation, dominated by the collar diameter, the number of primary branches per plant and plant height. This study clearly indicated that indigenous varieties are the most disease resistant genotypes and are having substantial fruit yield (945.30 g/plant) similar to standard varieties and at a touching distance to hybrid cultivars.

Two complement receptor one alleles have opposing associations with cerebral malaria and interact with α+thalassaemia.

Malaria has been a major driving force in the evolution of the human genome. In sub-Saharan African populations, two neighbouring polymorphisms in the Complement Receptor One (CR1) gene, named Sl2 and McCb, occur at high frequencies, consistent with selection by malaria. Previous studies have been inconclusive. Using a large case-control study of severe malaria in Kenyan children and statistical models adjusted for confounders, we estimate the relationship between Sl2 and McCb and malaria phenotypes, and find they have opposing associations. The Sl2 polymorphism is associated with markedly reduced odds of cerebral malaria and death, while the McCb polymorphism is associated with increased odds of cerebral malaria. We also identify an apparent interaction between Sl2 and α+thalassaemia, with the protective association of Sl2 greatest in children with normal α-globin. The complex relationship between these three mutations may explain previous conflicting findings, highlighting the importance of considering genetic interactions in disease-association studies.

Distinguishable brain networks relate disease susceptibility to symptom expression in schizophrenia.

Disease association studies have characterized altered resting-state functional connectivities describing schizophrenia, but failed to model symptom expression well. We developed a model that could account for symptom severity and meanwhile relate this to disease-related functional pathology. We correlated BOLD signal across brain regions and tested separately for associations with disease (disease edges) and with symptom severity (symptom edges) in a prediction-based scheme. We then integrated them in an "edge bi-color" graph, and adopted mediation analysis to test for causality between the disease and symptom networks and symptom scores. For first-episode schizophrenics (FES, 161 drug-naïve patients and 150 controls), the disease network (with inferior frontal gyrus being the hub) and the symptom-network (posterior occipital-parietal cortex being the hub) were found to overlap in the temporal lobe. For chronic schizophrenis (CS, 69 medicated patients and 62 controls), disease network was dominated by thalamocortical connectivities, and overlapped with symptom network in the middle frontal gyrus. We found that symptom network mediates the relationship between disease network and symptom scores in FEP, but was unable to define a relationship between them for the smaller CS population. Our results suggest that the disease network distinguishing core functional pathology in resting-state brain may be responsible for symptom expression in FES through a wider brain network associated with core symptoms. We hypothesize that top-down control from heteromodal prefrontal cortex to posterior transmodal cortex contributes to positive symptoms of schizophrenia. Our work also suggests differences in mechanisms of symptom expression between FES and CS, highlighting a need to distinguish between these groups.

Epigenome-wide association studies for systemic autoimmune diseases: The road behind and the road ahead

Epigenetics is known to be an important mechanism in the pathogenesis of autoimmune diseases. Epigenetic variations can act as integrators of environmental and genetic exposures and propagate activated states in immune cells. Studying epigenetic alterations by means of genome-wide approaches promises to unravel novel molecular mechanisms related to disease etiology, disease progression, clinical manifestations and treatment responses. This paper reviews what we have learned in the last five years from epigenome-wide studies for three systemic autoimmune diseases, namely systemic lupus erythematosus, primary Sjögren's syndrome, and rheumatoid arthritis. We examine the degree of epigenetic sharing between different diseases and the possible mediating role of epigenetic associations in genetic and environmental risks. Finally, we also shed light into the use of epigenetic markers towards a better precision medicine regarding disease prediction, prevention and personalized treatment in systemic autoimmunity.

Molecular Diagnosis of Felis catus Gammaherpesvirus 1 (FcaGHV1) Infection in Cats of Known Retrovirus Status with and without Lymphoma.

The pathogenicity of Felis catus gammaherpesvirus 1 (FcaGHV1), a common infection of domestic cats, is unknown. To explore an association between FcaGHV1 detection and feline lymphoma, a retrospective, cross-sectional, disease-association study was conducted. The infection status of all cats for feline immunodeficiency virus and feline leukaemia virus was determined. Neither a molecular diagnosis of FcaGHV1 nor whole-blood FcaGHV1 load was related to outcome in 122 lymphoma cases compared with 71 controls matched for age and sex. Molecular analysis of lymphoma-derived DNA paired with autologous uninvolved tissue did not suggest restriction of FcaGHV1 DNA to tumour tissue. FcaGHV1 DNA detection was associated with significantly shorter survival in lymphoma cases, an observation that could not be adequately explained by treatment differences. In addition, regressive feline leukaemia virus infection was identified as a risk factor for lymphoma. A history of fighting or roaming was identified as a novel epidemiological risk factor for FcaGHV1 detection, lending support to intercat aggression as a potential route of transmission. Studies investigating the cellular location and expression of FcaGHV1 are indicated to assist in ruling out a lymphomagenic role for this virus. Prospective investigation of FcaGHV1 DNA detection as a prognostic marker in feline lymphoma is warranted.

Method for identification of condition-associated public antigen receptor sequences.

Diverse repertoires of hypervariable immunoglobulin receptors (TCR and BCR) recognize antigens in the adaptive immune system. The development of immunoglobulin receptor repertoire sequencing methods makes it possible to perform repertoire-wide disease association studies of antigen receptor sequences. We developed a statistical framework for associating receptors to disease from only a small cohort of patients, with no need for a control cohort. Our method successfully identifies previously validated Cytomegalovirus and type one diabetes responsive TCR[Formula: see text] sequences .

In-silico study of toxicokinetics and disease association of chemicals present in smokeless tobacco products

Smokeless tobacco (SLT) products are consumed by millions of people in over 130 countries around the world. Consumption of SLT has been estimated to cause a number of diseases accounting to more than 0.65 million deaths per year. There is sufficient epidemiological evidence on the association of SLT products with nicotine addiction, cancers of oral cavity and digestive systems but there is a lack of understanding of the role of toxic chemicals in these diseases. We provide the first comprehensive in-silico analysis of chemical compounds present in different SLT products used worldwide. Many of these compounds are found to have good absorption, solubility and permeability along with mutagenic and toxic properties. They are also found to target more than 350 human proteins involved in a plethora of human biological processes and pathways. Along with all the previously known diseases, the present study has identified the association of compounds of SLT products with a number of unknown diseases like neurodegenerative, immune and cardiac diseases (Left ventricular non compaction, dilated cardiomyopathy etc). These findings indicate far-reaching impact of SLT products on human health than already known which needs further validations using epidemiological, in-vitro and in-vivo methodologies. Thus, this study will provide one stop information for the policy makers in development of regulatory policies on toxic contents of SLT products.

Parkinson's disease genetic risk in a midbrain neuronal cell line

In genome-wide association studies of complex diseases, many risk polymorphisms are found to lie in non-coding DNA and likely confer risk through allele-dependent differences in gene regulatory elements. However, because distal regulatory elements can alter gene expression at various distances on linear DNA, the identity of relevant genes is unknown for most risk loci. In Parkinson's disease, at least some genetic risk is likely intrinsic to a neuronal subpopulation of cells in the brain regions affected. In order to compare neuron-relevant methods of pairing risk polymorphisms to target genes as well as to further characterize a single-cell model of a neurodegenerative disease, we used the portionally-dopaminergic, neuronal, mesencephalic-derived cell line LUHMES to dissect differentiation-specific mechanisms of gene expression. We compared genome-wide gene expression in undifferentiated and differentiated cells with genome-wide histone H3K27ac and CTCF-bound regions. Whereas promoters and CTCF binding were largely consistent between differentiated and undifferentiated cells, enhancers were mostly unique. We matched the differentiation-specific appearance or disappearance of enhancers with changes in gene expression and identified 22,057 enhancers paired with 6388 differentially expressed genes by proximity. These enhancers are enriched with at least 13 transcription factor response elements, driving a cluster of genes involved in neurogenesis. We show that differentiated LUHMES cells, but not undifferentiated cells, show enrichment for PD-risk SNPs. Candidate genes for these loci are largely unrelated, though a subset is linked to synaptic vesicle cycling and transport, implying that PD-related disruption of these pathways is intrinsic to dopaminergic neurons.

The Five Dimensions of Parkinson's Disease Genetic Risk.

Genome-wide association studies of Parkinson's disease have revealed polymorphic variants associated with closely mapped genes of interest. We propose here that those genes may only represent the tip of an iceberg of regulatory effects and do not necessary reflect disease relevance. To usefully interpret a risk locus, one needs to consider 5 dimensions of information, which represent the three-dimensional structure of chromatin (dimensions #1- 3), which is locally variable across time (dimension #4), and, most importantly, dependent on cell type and context (dimension #5).

Genotyping in patients affected by HLA-related diseases. App development for diagnostic support.

HLA typing requests for association studies of immune-mediated diseases are often redundant and inadequate. We designed a series of meta-analyses to evaluate the accuracy of typing and distribution of HLA alleles predisposing to diseases, aiming at developing an app that can help doctors in choosing the most suitable molecular analysis. The first study was on celiac disease (CD) and HLA-DQ in children. We searched all english articles published in the main bibliographic databases up to May 2016. The search strategy has been developed using controlled terms (e.g. MeSH) and free terms. We identified 1885 articles. 1334 abstracts were examined. 46 manuscripts were evaluated, and 13 studies were included in the meta-analysis (740 CD and 943 controls). The risk of developing CD in children with allelic variants encoding the HLA-DQ2.5 and/or HLA-DQ8 molecules has been confirmed. The greatest CD risk resides in carriers of two DQ2.5 molecules, i.e. subjects homozygous for the DQB1*02:01 and DQA1*05 alleles (OR=5.4, 95 % CI=4.1-6.8) compared to any other DQ genotype. Carriers of two DQB1*02:01 (chain β2) alleles and one DQA1*05 (chain α5) allele have the same risk (p=0.8089) of DQ2.5 homozygotes (OR=5.3%, 95 CI=4,1 to 6.5). We found no differences between DQ8/β2 and DQ2.5/DQ8, nor between β2/DQX and DQ2.5/X. We suggest a two-step process: first typing the DQB1*02:01 allele and, in case of a negative result, full typing of HLA-DQ.

Quantitative CrAssphage real-time PCR assay derived from data of multiple geographically distant populations.

After its computational inference from human stool metagenomes, the CrAssphage has proven to be the most prevalent phage in the human gut, with presumably very wide geographic distribution. The currently available molecular assays do not sufficiently reflect the CrAssphage sequence variability. Here, we report a novel real-time PCR assay whose primers and probes are derived from data of multiple CrAssphage strains obtained from gut viral metagenomes of European, Asian, and African subjects. This assay can be useful in analyses of putative bacterial host co-occurence, and in association studies of non-infectious diseases where the phage may modify the content of gut bacteriomes.

Statistical analysis of human microarray data shows that dietary intervention with n-3 fatty acids, flavonoids and resveratrol enriches for immune response and disease pathways.

n-3 Fatty acids, flavonoids and resveratrol are well publicised for their beneficial effects on human health and wellbeing. Identifying common, underlying biological mechanisms targeted by these functional foods would therefore be informative for the public health sector for advising on nutritional health and disease, food and drug product development and consumer interest. The aim of this study was to explore the potential effects of gene expression changes associated with n-3 fatty acids EPA and DHA, flavonoids and resveratrol on modifying biological systems and disease pathways. To test this, publicly available human microarray data for significant gene expression changes associated with dietary intervention with EPA/DHA, flavonoids and resveratrol was subjected to pathway analysis and significance testing for overlap with signals from genome-wide association studies (GWAS) for common non-communicable diseases and biological functions. There was an enrichment of genes implicated in immune responses and disease pathways which was common to all of the treatment conditions tested. Analysis of biological functions and disease pathways indicated anti-tumorigenic properties for EPA/DHA. In line with this, significance testing of the intersection of genes associated with these functional foods and GWAS hits for common biological functions (ageing and cognition) and non-communicable diseases (breast cancer, CVD, diabesity, neurodegeneration and psychiatric disorders) identified significant overlap between the EPA/DHA and breast cancer gene sets. Dietary intervention with EPA/DHA, flavonoids and resveratrol can target important biological and disease pathways suggesting a potentially important role for these bioactive compounds in the prevention and treatment of dietary-related diseases.

High-Throughput Contiguous Full-Length Next-Generation Sequencing of HLA Class I andII Genes from 96 Donors in aSingle MiSeq Run.

The human leukocyte antigen (HLA) genes regulate and drive the immune system, and are among the most polymorphic loci in the human genome. HLA diversity is known to play an important role in transplantation and disease association studies. There are multiple approaches to DNA-based HLA genotyping and recent advances in next-generation sequencing (NGS) technologies have facilitated the development of whole gene sequencing methods. We describe an accurate, efficient, scalable, and cost-effective approach to contiguously amplify and sequence full-length genes of six HLA class I and II loci from 96 individuals on a single Illumina MiSeq run.

Human Leukocyte Antigen-A, B, C, DRB1, and DQB1 Allele and Haplotype Frequencies in a Subset of 237 Donors in the South African Bone Marrow Registry.

Human leukocyte antigen- (HLA-) A, HLA-B, HLA-C, HLA-DRB1, and HLA-DQB1 allele and haplotype frequencies were studied in a subset of 237 volunteer bone marrow donors registered at the South African Bone Marrow Registry (SABMR). Hapl-o-Mat software was used to compute allele and haplotype frequencies from individuals typed at various resolutions, with some alleles in multiple allele code (MAC) format. Four hundred and thirty-eight HLA-A, 235 HLA-B, 234 HLA-DRB1, 41 HLA-DQB1, and 29 HLA-C alleles are reported. The most frequent alleles were A∗02:02g (0.096), B∗07:02g (0.082), C∗07:02g (0.180), DQB1∗06:02 (0.157), and DRB1∗15:01 (0.072). The most common haplotype was A∗03:01g~B∗07:02g~C∗07:02g~DQB1∗06:02~DRB1∗15:01 (0.067), which has also been reported in other populations. Deviations from Hardy-Weinberg equilibrium were observed in A, B, and DRB1 loci, with C~DQB1 being the only locus pair in linkage disequilibrium. This study describes allele and haplotype frequencies from a subset of donors registered at SABMR, the only active bone marrow donor registry in Africa. Although the sample size was small, our results form a key resource for future population studies, disease association studies, and donor recruitment strategies.

DNA markers from genome-wide association studies of cardiovascular diseases

DNA markers from genome-wide association studies of cardiovascular diseases

Genome-Wide Association Studies of Autoimmune Thyroid Diseases, Thyroid Function, and Thyroid Cancer.

Thyroid diseases, including autoimmune thyroid diseases and thyroid cancer, are known to have high heritability. Family and twin studies have indicated that genetics plays a major role in the development of thyroid diseases. Thyroid function, represented by thyroid stimulating hormone (TSH) and free thyroxine (T4), is also known to be partly genetically determined. Before the era of genome-wide association studies (GWAS), the ability to identify genes responsible for susceptibility to thyroid disease was limited. Over the past decade, GWAS have been used to identify genes involved in many complex diseases, including various phenotypes of the thyroid gland. In GWAS of autoimmune thyroid diseases, many susceptibility loci associated with autoimmunity (human leukocyte antigen [HLA], protein tyrosine phosphatase, non-receptor type 22 [PTPN22], cytotoxic T-lymphocyte associated protein 4 [CTLA4], and interleukin 2 receptor subunit alpha [IL2RA]) or thyroid-specific genes (thyroid stimulating hormone receptor [TSHR] and forkhead box E1 [FOXE1]) have been identified. Regarding thyroid function, many susceptibility loci for levels of TSH and free T4 have been identified through genome-wide analyses. In GWAS of differentiated thyroid cancer, associations at FOXE1, MAP3K12 binding inhibitory protein 1 (MBIP)-NK2 homeobox 1 (NKX2-1), disrupted in renal carcinoma 3 (DIRC3), neuregulin 1 (NRG1), and pecanex-like 2 (PCNXL2) have been commonly identified in people of European and Korean ancestry, and many other susceptibility loci have been found in specific populations. Through GWAS of various thyroid-related phenotypes, many susceptibility loci have been found, providing insights into the pathogenesis of thyroid diseases and disease co-clustering within families and individuals.