Disruption of c-MYC Binding and Chromosomal Looping Involving Genetic Variants Associated With Ankylosing Spondylitis Upstream of the RUNX3 Promoter

Background: Ankylosing Spondylitis (AS) is a common form of inflammatory spinal arthritis with a complex aetiology and high heritability, involving more than 100 genetic associations. These include several AS-associated single nucleotide polymorphisms (SNPs) upstream of RUNX3, which encodes the multifunctional RUNT-related transcription factor (TF) 3. The lead associated SNP rs6600247 (p = 2.6 × 10−15) lies ∼13kb upstream of the RUNX3 promoter adjacent to a c-MYC TF binding-site. The effect of rs6600247 genotype on DNA binding and chromosome looping were investigated by electrophoretic mobility gel shift assays (EMSA), Western blotting-EMSA (WEMSA) and Chromosome Conformation Capture (3C).Results: Interrogation of ENCODE published data showed open chromatin in the region overlapping rs6600247 in primary human CD14+ monocytes, in contrast to the Jurkat T cell line or primary human T-cells. The rs6600247 AS-risk allele is predicted to specifically disrupt a c-MYC binding-site. Using a 50bp DNA probe spanning rs6600247 we consistently observed reduced binding to the AS-risk “C” allele of both purified c-MYC protein and nuclear extracts (NE) from monocyte-like U937 cells. WEMSA on U937 NE and purified c-MYC protein confirmed these differences (n = 3; p < 0.05). 3C experiments demonstrated negligible interaction between the region encompassing rs6600247 and the RUNX3 promoter. A stronger interaction frequency was demonstrated between the RUNX3 promoter and the previously characterised AS-associated SNP rs4648889.Conclusion: The lead SNP rs6600247, located in an enhancer-like region upstream of the RUNX3 promoter, modulates c-MYC binding. However, the region encompassing rs6600247 has rather limited physical interaction with the promoter of RUNX3. In contrast a clear chromatin looping event between the region encompassing rs4648889 and the RUNX3 promoter was observed. These data provide further evidence for complexity in the regulatory elements upstream of the RUNX3 promoter and the involvement of RUNX3 transcriptional regulation in AS.

Multivariate GWAS of Structural Dental Anomalies and Dental Caries in a Multi-Ethnic Cohort

Odontogenesis is a complex process, where disruption can result in dental anomalies and/or increase the risk of developing dental caries. Based on previous studies, certain dental anomalies tend to co-occur in patients, suggesting that these traits may share common genetic and etiological components. The main goal of this study was to implement a multivariate genome wide association study approach to identify genetic variants shared between correlated structural dental anomalies and dental caries. Our cohort (N = 3,579) was derived from the Pittsburgh Orofacial Clefts Study, where multiple dental traits were assessed in both the unaffected relatives of orofacial cleft (OFC) cases (n = 2,187) and unaffected controls (n = 1,392). We identified four multivariate patterns of correlated traits in this data: tooth agenesis, impaction, and rotation (AIR); enamel hypoplasia, displacement, and rotation (HDR); displacement, rotation, and mamelon (DRM); and dental caries, tooth agenesis and enamel hypoplasia (CAH). We analyzed each of these four models using genome-wide multivariate tests of association. No genome-wide statistically significant results were found, but we identified multiple suggestive association signals (P ≤ 10−5) near genes with known biological roles during tooth development, including ADAMTS9 and PRICKLE2 associated with AIR; GLIS3, WDR72, and ROR2 associated with HDR and DRM; ROBO2 associated with DRM; BMP7 associated with HDR; and ROBO1, SMAD2, and MSX2 associated with CAH. This is the first study to investigative genetic associations for multivariate patterns of correlated dental anomalies and dental caries. Further studies are needed to replicate these results in independent cohorts.

Genomics and predictive medicine

Progress in understanding of structural and functional human genome organization and deciphering primary DNA sequence in human cells allowed for hitherto unreachable new capabilities of medical genetics in identifying the causes and mechanisms of inherited and inborn pathology. Implementation of genetics into medicine is progressively advancing along with improvement of molecular analysis of genome. Knowledge of genome and its functions allows to provide more accurate diagnosis, predict, to a considerable extent, the presence of genetic predisposition of a person to pathology, and to assess the chances for developing one or another disease. This approach became the basis for a new area of medical genetics named predictive medicine. The progress of predictive medicine refl ects success in tremendous upgrowth of molecular genetic methods and new capabilities of studying structure and functions of genome. Within less than 15 years after deciphering genome, medical genetics has travelled a long way from a single gene analysis to whole genome studies, from screening of genetic associations to systems genetics of multifactorial diseases, from translational to high-precision genetics, and from genetic passport idea to electronic genetic health records. The development of a genetic passport, prognostic genetic testing, and genomic chart of reproductive health is especially relevant for current practical medicine.

Genetic Factors, Brain Atrophy, and Response to Rehabilitation Therapy After Stroke

Objective Patients show substantial differences in response to rehabilitation therapy after stroke. We hypothesized that specific genetic profiles might explain some of this variance and, secondarily, that genetic factors are related to cerebral atrophy post-stroke. Methods The phase 3 ICARE study examined response to motor rehabilitation therapies. In 216 ICARE enrollees, DNA was analyzed for presence of the BDNF val66met and the ApoE ε4 polymorphism. The relationship of polymorphism status to 12-month change in motor status (Wolf Motor Function Test, WMFT) was examined. Neuroimaging data were also evaluated (n=127). Results Subjects were 61±13 years old (mean±SD) and enrolled 43±22 days post-stroke; 19.7% were BDNF val66met carriers and 29.8% ApoE ε4 carriers. Carrier status for each polymorphism was not associated with WMFT, either at baseline or over 12 months of follow-up. Neuroimaging, acquired 5±11 days post-stroke, showed that BDNF val66met polymorphism carriers had a 1.34-greater degree of cerebral atrophy compared to non-carriers (P=.01). Post hoc analysis found that age of stroke onset was 4.6 years younger in subjects with the ApoE ε4 polymorphism (P=.02). Conclusion Neither the val66met BDNF nor ApoE ε4 polymorphism explained inter-subject differences in response to rehabilitation therapy. The BDNF val66met polymorphism was associated with cerebral atrophy at baseline, echoing findings in healthy subjects, and suggesting an endophenotype. The ApoE ε4 polymorphism was associated with younger age at stroke onset, echoing findings in Alzheimer’s disease and suggesting a common biology. Genetic associations provide insights useful to understanding the biology of outcomes after stroke.

A Network-Based Analysis of Disease Complication Associations for Obstetric Disorders in the UK Biobank

Background: Recent studies have found that women with obstetric disorders are at increased risk for a variety of long-term complications. However, the underlying pathophysiology of these connections remains undetermined. A network-based view incorporating knowledge of other diseases and genetic associations will aid our understanding of the role of genetics in pregnancy-related disease complications. Methods: We built a disease–disease network (DDN) using UK Biobank (UKBB) summary data from a phenome-wide association study (PheWAS) to elaborate multiple disease associations. We also constructed egocentric DDNs, where each network focuses on a pregnancy-related disorder and its neighboring diseases. We then applied graph-based semi-supervised learning (GSSL) to translate the connections in the egocentric DDNs to pathologic knowledge. Results: A total of 26 egocentric DDNs were constructed for each pregnancy-related phenotype in the UKBB. Applying GSSL to each DDN, we obtained complication risk scores for additional phenotypes given the pregnancy-related disease of interest. Predictions were validated using co-occurrences derived from UKBB electronic health records. Our proposed method achieved an increase in average area under the receiver operating characteristic curve (AUC) by a factor of 1.35 from 55.0% to 74.4% compared to the use of the full DDN. Conclusion: Egocentric DDNs hold promise as a clinical tool for the network-based identification of potential disease complications for a variety of phenotypes.

A Common R219K Variant of ATP-Binding Cassette Transporter A1 Gene Alters Atherometabolic Traits in Pregnant Women With Gestational Diabetes Mellitus

BackgroundATP-binding cassette transporter A1 (ABCA1) has important roles in high-density lipoprotein (HDL) metabolism and reverse cholesterol transport, and is implicated in lipid-related disorders. Genetic variants are involved in the pathogenesis of gestational diabetes mellitus (GDM). The objective of this study was to investigate the association of rs2230806 (R219K), a single nucleotide polymorphism (SNP) in the lipid-related gene, with the risk of GDM and related traits.MethodsThe SNP, rs2230806, was genotyped, and clinical and metabolic parameters were determined in 660 GDM patients and 1,097 control subjects. Genetic associations with related traits were also analyzed.ResultsThe genotype distributions were similar in GDM patients and normal controls. However, significant differences in the variables examined in the study subjects were noted across the three genotypes. The genotype at the rs2230806 polymorphism was significantly associated with HDL-cholesterol (HDL-C) levels and atherogenic index (AI) values in GDM patients and total cholesterol (TC) and low-density lipoprotein-cholesterol (LDL-C) levels in control subjects. Subgroup analysis showed that the polymorphism was associated with diastolic blood pressure, in addition to HDL-C levels and AI, in overweight/obese GDM patients, while it was associated with TC levels, AI, pre-pregnancy body mass index (BMI), and BMI at delivery in non-obese GDM patients. In addition, this polymorphism was associated with TC, LDL-C, and apoB levels in overweight/obese control subjects.ConclusionsThe rs2230806 polymorphism in the ABCA1 gene was associated with variations in atherometabolic traits in GDM patients, with characteristics of BMI dependency, but not with GDM. Our findings highlight a link between related phenotypes in women with GDM and genetic factors.

Shared Inherited Genetics of Benign Prostatic Hyperplasia and Prostate Cancer

BackgroundThe association between benign prostatic hyperplasia (BPH) and prostate cancer (PCa) remains controversial, largely due to inherent detection bias in traditional observational studies. The objective of this study is to assess their association using inherited SNPs.MethodsSubjects were White men from the large population-based UK Biobank (UKB). Association between BPH and PCa was tested: 1) phenotypical correlation using chi-square test, 2) genetic correlation (rg) based on 1,126,841 polymorphic SNPs across the genome using linkage disequilibrium score regression (LDSR), and 3) cross-disease genetic associations based on known risk-associated SNPs (15 for BPH and 239 for PCa), individually and cumulatively as measured by genetic risk score (GRS).FindingsAmong 214,717 White men in the UKB, 24,623 (11.47%) and 14,311 (6.67%) had a diagnosis of BPH and PCa, respectively. Diagnoses of these two diseases were significantly correlated, χ2=1862.80, P<1E-299. A significant genetic correlation was found, rg (95% confidence interval (CI))=0.27 (0.15-0.39), P=9.17E-06. In addition, significant cross-disease genetic associations for established risk-associated SNPs were also found. Among the 250 established GWAS-significant SNPs of PCa or BPH, 51 were significantly associated with risk of the other disease at P<0.05, significantly more than expected by chance (N=12), P=3.04E-7 (χ2-test). Furthermore, significant cross-disease GRS associations were also found; GRSBPH was significantly associated with PCa risk (odds ratio (OR)=1.26 (1.18-1.36), P=1.62E-10), and GRSPCa was significantly associated with BPH risk (OR=1.03 (1.02-1.04), P=8.57E-06). Moreover, GRSBPH was significantly and inversely associated with lethal PCa risk in a PCa case-case analysis (OR=0.58 (0.41-0.81), P=1.57E-03). In contrast, GRSPCa was not significantly associated with lethal PCa (OR=0.99 (0.94-1.04), P=0.79).InterpretationBPH and PCa share common inherited genetics which suggests the phenotypical association of these two diseases in observational studies is not entirely caused by detection bias. This novel finding may have implications in disease etiology and risk stratification.FundingNone.

A Precise Method to Evaluate 360 Degree Measures of Optic Cup and Disc Morphology in an African American Cohort and Its Genetic Applications

(1) Background: Vertical cup-to-disc ratio (CDR) is an important measure for evaluating damage to the optic nerve head (ONH) in glaucoma patients. However, this measure often does not fully capture the irregular cupping observed in glaucomatous nerves. We developed and evaluated a method to measure cup-to-disc ratio (CDR) at all 360 degrees of the ONH. (2) Methods: Non-physician graders from the Scheie Reading Center outlined the cup and disc on digital stereo color disc images from African American patients enrolled in the Primary Open-Angle African American Glaucoma Genetics (POAAGG) study. After converting the resultant coordinates into polar representation, the CDR at each 360-degree location of the ONH was obtained. We compared grader VCDR values with clinical VCDR values, using Spearman correlation analysis, and validated significant genetic associations with clinical VCDR, using grader VCDR values. (3) Results: Graders delineated outlines of the cup contour and disc boundaries twice in each of 1815 stereo disc images. For both cases and controls, the mean CDR was highest at the horizontal bisector, particularly in the temporal region, as compared to other degree locations. There was a good correlation between grader CDR at the vertical bisector and clinical VCDR (Spearman Correlation OD: r = 0.78 [95% CI: 0.76–0.79]). An SNP in the MPDZ gene, associated with clinical VCDR in a prior genome-wide association study, showed a significant association with grader VCDR (p = 0.01) and grader CDR area ratio (p = 0.02). (4) Conclusions: The CDR of both glaucomatous and non-glaucomatous eyes varies by degree location, with the highest measurements in the temporal region of the eye. This method can be useful for capturing innate eccentric ONH morphology, tracking disease progression, and identifying genetic associations.