Genetic complexity of diagnostically unresolved Ehlers-Danlos syndrome

BACKGROUND

The Ehlers-Danlos syndromes (EDS) are heritable disorders of connective tissue (HDCT), reclassified in the 2017 nosology into 13 subtypes. The genetic basis for hypermobile Ehlers-Danlos syndrome (hEDS) remains unknown.

METHODS

Whole exome sequencing (WES) was undertaken on 174 EDS patients recruited from a national diagnostic service for complex EDS and a specialist clinic for hEDS. Patients had already undergone expert phenotyping, laboratory investigation and gene sequencing, but were without a genetic diagnosis. Filtered WES data were reviewed for genes underlying Mendelian disorders and loci reported in EDS linkage, transcriptome and genome-wide association studies (GWAS). A genetic burden analysis (Minor Allele Frequency (MAF) <0.05) incorporating 248 Avon Longitudinal Study of Parents and Children (ALSPAC) controls sequenced as part of the UK10K study was undertaken using TASER methodology.

RESULTS

Heterozygous pathogenic (P) or likely pathogenic (LP) variants were identified in known EDS and Loeys-Dietz (LDS) genes. Multiple variants of uncertain significance where segregation and functional analysis may enable reclassification were found in genes associated with EDS, LDS, heritable thoracic aortic disease (HTAD), Mendelian disorders with EDS symptomatology and syndromes with EDS-like features. Genetic burden analysis revealed a number of novel loci, although none reached the threshold for genome-wide significance. Variants with biological plausibility were found in genes and pathways not currently associated with EDS or HTAD.

CONCLUSIONS

We demonstrate the clinical utility of large panel-based sequencing and WES for patients with complex EDS in distinguishing rare EDS subtypes, LDS and related syndromes. Although many of the P and LP variants reported in this cohort would be identified with current panel testing, they were not at the time of this study, highlighting the use of extended panels and WES as a clinical tool for complex EDS. Our results are consistent with the complex genetic architecture of EDS and suggest a number of novel hEDS and HTAD candidate genes and pathways.

All-cause and liver-related mortality risk factors in excessive drinkers: Analysis of data from the UK biobank

BACKGROUND

High alcohol intake is associated with increased mortality. We aimed to identify factors affecting mortality in people drinking extreme amounts of alcohol.

METHODS

We obtained information from the UK Biobank on approximately 500,000 participants aged 40-70 years at baseline assessment in 2006-2010. Habitual alcohol intake, lifestyle and physiological data, laboratory test results, and hospital diagnoses and death certificate data (to June 2020) for 5136 men (2.20% of male participants) and 1504 women (0.60%) who reported consuming ≥80 or ≥50 g/day, respectively, were used in survival analysis.

RESULTS

Mortality hazard ratios for these excessive drinkers, compared to all other participants, were 2.02 (95% CI 1.89-2.17) for all causes, 1.89 (1.69-2.12) for any cancer, 1.87 (1.61-2.17) for any circulatory disease, and 9.40 (7.00-12.64) for any liver disease. Liver disease diagnosis or abnormal liver function tests predicted not only deaths attributed to liver disease but also those from cancers or circulatory diseases. Mortality among excessive drinkers was also associated with quantitative alcohol intake; diagnosed alcohol dependence, harmful use, or withdrawal syndrome; and current smoking at assessment.

CONCLUSIONS

People with chronic excessive alcohol intake experience decreased average survival, but there is substantial variation in their mortality, with liver abnormality and alcohol dependence or other alcohol use disorders associated with a worse prognosis. Clinically, patients with these risk factors and high alcohol intake should be considered for early or intensive management. Research can usefully focus on the factors predisposing to dependence or liver abnormality.

The influence of HLA genotype on the development of metal hypersensitivity following joint replacement

Background

Over five million joint replacements are performed across the world each year. Cobalt chrome (CoCr) components are used in most of these procedures. Some patients develop delayed-type hypersensitivity (DTH) responses to CoCr implants, resulting in tissue damage and revision surgery. DTH is unpredictable and genetic links have yet to be definitively established.

Methods

At a single site, we carried out an initial investigation to identify HLA alleles associated with development of DTH following metal-on-metal hip arthroplasty. We then recruited patients from other centres to train and validate an algorithm incorporating patient age, gender, HLA genotype, and blood metal concentrations to predict the development of DTH. Accuracy of the modelling was assessed using performance metrics including time-dependent receiver operator curves.

Results

Using next-generation sequencing, here we determine the HLA genotypes of 606 patients. 176 of these patients had experienced failure of their prostheses; the remaining 430 remain asymptomatic at a mean follow up of twelve years. We demonstrate that the development of DTH is associated with patient age, gender, the magnitude of metal exposure, and the presence of certain HLA class II alleles. We show that the predictive algorithm developed from this investigation performs to an accuracy suitable for clinical use, with weighted mean survival probability errors of 1.8% and 3.1% for pre-operative and post-operative models respectively.

Conclusions

The development of DTH following joint replacement appears to be determined by the interaction between implant wear and a patient's genotype. The algorithm described in this paper may improve implant selection and help direct patient surveillance following surgery. Further consideration should be given towards understanding patient-specific responses to different biomaterials.

A multiancestry genome-wide association study of unexplained chronic ALT elevation as a proxy for nonalcoholic fatty liver disease with histological and radiological validation

Nonalcoholic fatty liver disease (NAFLD) is a growing cause of chronic liver disease. Using a proxy NAFLD definition of chronic elevation of alanine aminotransferase (cALT) levels without other liver diseases, we performed a multiancestry genome-wide association study (GWAS) in the Million Veteran Program (MVP) including 90,408 cALT cases and 128,187 controls. Seventy-seven loci exceeded genome-wide significance, including 25 without prior NAFLD or alanine aminotransferase associations, with one additional locus identified in European American-only and two in African American-only analyses (P < 5 × 10-8). External replication in histology-defined NAFLD cohorts (7,397 cases and 56,785 controls) or radiologic imaging cohorts (n = 44,289) replicated 17 single-nucleotide polymorphisms (SNPs) (P < 6.5 × 10-4), of which 9 were new (TRIB1, PPARG, MTTP, SERPINA1, FTO, IL1RN, COBLL1, APOH and IFI30). Pleiotropy analysis showed that 61 of 77 multiancestry and all 17 replicated SNPs were jointly associated with metabolic and/or inflammatory traits, revealing a complex model of genetic architecture. Our approach integrating cALT, histology and imaging reveals new insights into genetic liability to NAFLD.

Macrophage scavenger receptor 1 mediates lipid-induced inflammation in non-alcoholic fatty liver disease

BACKGROUND & AIMS

Obesity-associated inflammation is a key player in the pathogenesis of non-alcoholic fatty liver disease (NAFLD). However, the role of macrophage scavenger receptor 1 (MSR1, CD204) remains incompletely understood.

METHODS

A total of 170 NAFLD liver biopsies were processed for transcriptomic analysis and correlated with clinicopathological features. Msr1^-/- and wild-type mice were subjected to a 16-week high-fat and high-cholesterol diet. Mice and ex vivo human liver slices were treated with a monoclonal antibody against MSR1. Genetic susceptibility was assessed using genome-wide association study data from 1,483 patients with NAFLD and 430,101 participants of the UK Biobank.

RESULTS

MSR1 expression was associated with the occurrence of hepatic lipid-laden foamy macrophages and correlated with the degree of steatosis and steatohepatitis in patients with NAFLD. Mice lacking Msr1 were protected against diet-induced metabolic disorder, showing fewer hepatic foamy macrophages, less hepatic inflammation, improved dyslipidaemia and glucose tolerance, and altered hepatic lipid metabolism. Upon induction by saturated fatty acids, MSR1 induced a pro-inflammatory response via the JNK signalling pathway. In vitro blockade of the receptor prevented the accumulation of lipids in primary macrophages which inhibited the switch towards a pro-inflammatory phenotype and the release of cytokines such as TNF-ɑ. Targeting MSR1 using monoclonal antibody therapy in an obesity-associated NAFLD mouse model and human liver slices resulted in the prevention of foamy macrophage formation and inflammation. Moreover, we identified that rs41505344, a polymorphism in the upstream transcriptional region of MSR1, was associated with altered serum triglycerides and aspartate aminotransferase levels in a cohort of over 400,000 patients.

CONCLUSIONS

Taken together, our data suggest that MSR1 plays a critical role in lipid-induced inflammation and could thus be a potential therapeutic target for the treatment of NAFLD.

LAY SUMMARY

Non-alcoholic fatty liver disease (NAFLD) is a chronic disease primarily caused by excessive consumption of fat and sugar combined with a lack of exercise or a sedentary lifestyle. Herein, we show that the macrophage scavenger receptor MSR1, an innate immune receptor, mediates lipid uptake and accumulation in Kupffer cells, resulting in liver inflammation and thereby promoting the progression of NAFLD in humans and mice.

A genetic risk score and diabetes predict development of alcohol-related cirrhosis in drinkers

BACKGROUND & AIMS

Only a minority of excess alcohol drinkers develop cirrhosis. We developed and evaluated risk stratification scores to identify those at highest risk.

METHODS

Three cohorts (GenomALC-1: n = 1,690, GenomALC-2: n = 3,037, UK Biobank: relevant n = 6,898) with a history of heavy alcohol consumption (≥80 g/day (men), ≥50 g/day (women), for ≥10 years) were included. Cases were participants with alcohol-related cirrhosis. Controls had a history of similar alcohol consumption but no evidence of liver disease. Risk scores were computed from up to 8 genetic loci identified previously as associated with alcohol-related cirrhosis and 3 clinical risk factors. Score performance for the stratification of alcohol-related cirrhosis risk was assessed and compared across the alcohol-related liver disease spectrum, including hepatocellular carcinoma (HCC).

RESULTS

A combination of 3 single nucleotide polymorphisms (SNPs) (PNPLA3:rs738409, SUGP1-TM6SF2:rs10401969, HSD17B13:rs6834314) and diabetes status best discriminated cirrhosis risk. The odds ratios (ORs) and (95% CIs) between the lowest (Q1) and highest (Q5) score quintiles of the 3-SNP score, based on independent allelic effect size estimates, were 5.99 (4.18-8.60) (GenomALC-1), 2.81 (2.03-3.89) (GenomALC-2), and 3.10 (2.32-4.14) (UK Biobank). Patients with diabetes and high risk scores had ORs of 14.7 (7.69-28.1) (GenomALC-1) and 17.1 (11.3-25.7) (UK Biobank) compared to those without diabetes and with low risk scores. Patients with cirrhosis and HCC had significantly higher mean risk scores than patients with cirrhosis alone (0.76 ± 0.06 vs. 0.61 ± 0.02, p = 0.007). Score performance was not significantly enhanced by information on additional genetic risk variants, body mass index or coffee consumption.

CONCLUSIONS

A risk score based on 3 genetic risk variants and diabetes status enables the stratification of heavy drinkers based on their risk of cirrhosis, allowing for the provision of earlier preventative interventions.

LAY SUMMARY

Excessive chronic drinking leads to cirrhosis in some people, but so far there is no way to identify those at high risk of developing this debilitating disease. We developed a genetic risk score that can identify patients at high risk. The risk of cirrhosis is increased >10-fold with just two risk factors - diabetes and a high genetic risk score. Risk assessment using this test could enable the early and personalised management of this disease in high-risk patients.

Genome-wide association study identifies susceptibility loci for acute myeloid leukemia

Acute myeloid leukemia (AML) is a hematological malignancy with an undefined heritable risk. Here we perform a meta-analysis of three genome-wide association studies, with replication in a fourth study, incorporating a total of 4018 AML cases and 10488 controls. We identify a genome-wide significant risk locus for AML at 11q13.2 (rs4930561; P = 2.15 × 10-8; KMT5B). We also identify a genome-wide significant risk locus for the cytogenetically normal AML sub-group (N = 1287) at 6p21.32 (rs3916765; P = 1.51 × 10-10; HLA). Our results inform on AML etiology and identify putative functional genes operating in histone methylation (KMT5B) and immune function (HLA).

An international genome-wide meta-analysis of primary biliary cholangitis: Novel risk loci and candidate drugs

BACKGROUNDS & AIMS

Primary biliary cholangitis (PBC) is a chronic liver disease in which autoimmune destruction of the small intrahepatic bile ducts eventually leads to cirrhosis. Many patients have inadequate response to licensed medications, motivating the search for novel therapies. Previous genome-wide association studies (GWAS) and meta-analyses (GWMA) of PBC have identified numerous risk loci for this condition, providing insight into its aetiology. We undertook the largest GWMA of PBC to date, aiming to identify additional risk loci and prioritise candidate genes for in silico drug efficacy screening.

METHODS

We combined new and existing genotype data for 10,516 cases and 20,772 controls from 5 European and 2 East Asian cohorts.

RESULTS

We identified 56 genome-wide significant loci (20 novel) including 46 in European, 13 in Asian, and 41 in combined cohorts; and a 57th genome-wide significant locus (also novel) in conditional analysis of the European cohorts. Candidate genes at newly identified loci include FCRL3, INAVA, PRDM1, IRF7, CCR6, CD226, and IL12RB1, which each play key roles in immunity. Pathway analysis reiterated the likely importance of pattern recognition receptor and TNF signalling, JAK-STAT signalling, and differentiation of T helper (TH)1 and TH17 cells in the pathogenesis of this disease. Drug efficacy screening identified several medications predicted to be therapeutic in PBC, some of which are well-established in the treatment of other autoimmune disorders.

CONCLUSIONS

This study has identified additional risk loci for PBC, provided a hierarchy of agents that could be trialled in this condition, and emphasised the value of genetic and genomic approaches to drug discovery in complex disorders.

LAY SUMMARY

Primary biliary cholangitis (PBC) is a chronic liver disease that eventually leads to cirrhosis. In this study, we analysed genetic information from 10,516 people with PBC and 20,772 healthy individuals recruited in Canada, China, Italy, Japan, the UK, or the USA. We identified several genetic regions associated with PBC. Each of these regions contains several genes. For each region, we used diverse sources of evidence to help us choose the gene most likely to be involved in causing PBC. We used these 'candidate genes' to help us identify medications that are currently used for treatment of other conditions, which might also be useful for treatment of PBC.

The influence of HLA genotype on the severity of COVID-19 infection

The impact of COVID‐19 varies markedly, not only between individual patients but also between different populations. We hypothesised that differences in HLA genes might influence this variation. Using next generation sequencing, we analysed the class I and class II classical HLA genes of 147 individuals of European descent experiencing variable clinical outcomes following COVID‐19 infection. Forty‐nine of these patients were admitted to hospital with severe respiratory disease. They had no significant pre‐existing comorbidities. We compared the results to those obtained from a group of 69 asymptomatic hospital workers who evidence of COVID exposure based on blood antibody testing. Allele frequencies in both the severe and asymptomatic groups were compared to local and national healthy controls with adjustments made for age and sex. With the inclusion of hospital staff who had reported localised symptoms only (limited to loss of smell/taste, n = 13) or systemic symptoms not requiring hospital treatment (n = 16), we carried out ordinal logistic regression modelling to determine the relative influence of age, BMI, sex and the presence of specific HLA genes on symptomatology. We found a significant difference in the allele frequency of HLA‐DRB1*04:01 in the severe patient compared to the asymptomatic staff group (5.1% vs. 16.7%, P = .003 after adjustment for age and sex). There was a significantly lower frequency of the haplotype DQA1*01:01‐DQB1*05:01‐DRB1*01:01 in the asymptomatic group compared to the background population (P = .007). Ordinal logistic regression modelling confirmed the significant influence of DRB1*04:01 on the clinical severity of COVID‐19 observed in the cohorts. These alleles are found in greater frequencies in the North Western European population. This regional study provides evidence that HLA genotype influences clinical outcome in COVID‐19 infection. Validation studies must take account of the complex genetic architecture of the immune system across different geographies and ethnicities.

Genome-wide Association Study and Meta-analysis on Alcohol-Associated Liver Cirrhosis Identifies Genetic Risk Factors

BACKGROUND AND AIMS

Only a minority of heavy drinkers progress to alcohol-associated cirrhosis (ALC). The aim of this study was to identify common genetic variants that underlie risk for ALC.

APPROACH AND RESULTS

We analyzed data from 1,128 subjects of European ancestry with ALC and 614 heavy-drinking subjects without known liver disease from Australia, the United States, the United Kingdom, and three countries in Europe. A genome-wide association study (GWAS) was performed, adjusting for principal components and clinical covariates (alcohol use, age, sex, body mass index, and diabetes). We validated our GWAS findings using UK Biobank. We then performed a meta-analysis combining data from our study, the UK Biobank, and a previously published GWAS. Our GWAS found genome-wide significant risk association of rs738409 in patatin-like phospholipase domain containing 3 (PNPLA3) (odds ratio [OR] = 2.19 [G allele], P = 4.93 × 10^-17 ) and rs4607179 near HSD17B13 (OR = 0.57 [C allele], P = 1.09 × 10^-10 ) with ALC. Conditional analysis accounting for the PNPLA3 and HSD17B13 loci identified a protective association at rs374702773 in Fas-associated factor family member 2 (FAF2) (OR = 0.61 [del(T) allele], P = 2.56 × 10^-8 ) for ALC. This association was replicated in the UK Biobank using conditional analysis (OR = 0.79, P = 0.001). Meta-analysis (without conditioning) confirmed genome-wide significance for the identified FAF2 locus as well as PNPLA3 and HSD17B13. Two other previously known loci (SERPINA1 and SUGP1/TM6SF2) were also genome-wide significant in the meta-analysis. GeneOntology pathway analysis identified lipid droplets as the target for several identified genes. In conclusion, our GWAS identified a locus at FAF2 associated with reduced risk of ALC among heavy drinkers. Like the PNPLA3 and HSD17B13 gene products, the FAF2 product has been localized to fat droplets in hepatocytes.

CONCLUSIONS

Our genetic findings implicate lipid droplets in the biological pathway(s) underlying ALC.

A PDCD1 Role in the Genetic Predisposition to NAFLD-HCC?

Obesity and non-alcoholic fatty liver disease (NAFLD) are contributing to the global rise in deaths from hepatocellular carcinoma (HCC). The pathogenesis of NAFLD-HCC is not well understood. The severity of hepatic steatosis, steatohepatitis and fibrosis are key pathogenic mechanisms, but animal studies suggest altered immune responses are also involved. Genetic studies have so far highlighted a major role of gene variants promoting fat deposition in the liver (PNPLA3 rs738409; TM6SF2 rs58542926). Here, we have considered single-nucleotide polymorphisms (SNPs) in candidate immunoregulatory genes (MICA rs2596542; CD44 rs187115; PDCD1 rs7421861 and rs10204525), in 594 patients with NAFLD and 391 with NAFLD-HCC, from three European centres. Associations between age, body mass index, diabetes, cirrhosis and SNPs with HCC development were explored. PNPLA3 and TM6SF2 SNPs were associated with both progression to cirrhosis and NAFLD-HCC development, while PDCD1 SNPs were specifically associated with NAFLD-HCC risk, regardless of cirrhosis. PDCD1 rs7421861 was independently associated with NAFLD-HCC development, while PDCD1 rs10204525 acquired significance after adjusting for other risks, being most notable in the smaller numbers of women with NAFLD-HCC. The study highlights the potential impact of inter individual variation in immune tolerance induction in patients with NAFLD, both in the presence and absence of cirrhosis.

Transcriptomic profiling across the nonalcoholic fatty liver disease spectrum reveals gene signatures for steatohepatitis and fibrosis

The mechanisms that drive nonalcoholic fatty liver disease (NAFLD) remain incompletely understood. This large multicenter study characterized the transcriptional changes that occur in liver tissue across the NAFLD spectrum as disease progresses to cirrhosis to identify potential circulating markers. We performed high-throughput RNA sequencing on a discovery cohort comprising histologically characterized NAFLD samples from 206 patients. Unsupervised clustering stratified NAFLD on the basis of disease activity and fibrosis stage with differences in age, aspartate aminotransferase (AST), type 2 diabetes mellitus, and carriage of PNPLA3 rs738409, a genetic variant associated with NAFLD. Relative to early disease, we consistently identified 25 differentially expressed genes as fibrosing steatohepatitis progressed through stages F2 to F4. This 25-gene signature was independently validated by logistic modeling in a separate replication cohort (n = 175), and an integrative analysis with publicly available single-cell RNA sequencing data elucidated the likely relative contribution of specific intrahepatic cell populations. Translating these findings to the protein level, SomaScan analysis in more than 300 NAFLD serum samples confirmed that circulating concentrations of proteins AKR1B10 and GDF15 were strongly associated with disease activity and fibrosis stage. Supporting the biological plausibility of these data, in vitro functional studies determined that endoplasmic reticulum stress up-regulated expression of AKR1B10, GDF15, and PDGFA, whereas GDF15 supplementation tempered the inflammatory response in macrophages upon lipid loading and lipopolysaccharide stimulation. This study provides insights into the pathophysiology of progressive fibrosing steatohepatitis, and proof of principle that transcriptomic changes represent potentially tractable and clinically relevant markers of disease progression.

Amino acid residues in five separate HLA genes can explain most of the known associations between the MHC and primary biliary cholangitis

Primary Biliary Cholangitis (PBC) is a chronic autoimmune liver disease characterised by progressive destruction of intrahepatic bile ducts. The strongest genetic association is with HLA-DQA1*04:01, but at least three additional independent HLA haplotypes contribute to susceptibility. We used dense single nucleotide polymorphism (SNP) data in 2861 PBC cases and 8514 controls to impute classical HLA alleles and amino acid polymorphisms using state-of-the-art methodologies. We then demonstrated through stepwise regression that association in the HLA region can be largely explained by variation at five separate amino acid positions. Three-dimensional modelling of protein structures and calculation of electrostatic potentials for the implicated HLA alleles/amino acid substitutions demonstrated a correlation between the electrostatic potential of pocket P6 in HLA-DP molecules and the HLA-DPB1 alleles/amino acid substitutions conferring PBC susceptibility/protection, highlighting potential new avenues for future functional investigation.

Primary Biliary Cholangitis (PBC) is a chronic autoimmune liver disease that exhibits strong genetic associations, especially with variants in the human leukocyte antigen (HLA) gene region. Here we use dense single nucleotide polymorphism (SNP) data from the largest PBC study to date (2861 cases, 8514 controls) to investigate the likely underlying causes of this association, via performing imputation of HLA classical alleles and amino acids. We show that the HLA association can be largely explained by variation at five separate amino acid positions, one of which shows functional relevance to electrostatic potentials of HLA-DP molecules.

Transcriptional regulation of PNPLA3 and its impact on susceptibility to nonalcoholic fatty liver Disease (NAFLD) in humans

The increased expression of PNPLA3^148M leads to hepatosteatosis in mice. This study aims to investigate the genetic control of hepatic PNPLA3 transcription and to explore its impact on NAFLD risk in humans. Through a locus-wide expression quantitative trait loci (eQTL) mapping in two human liver sample sets, a PNPLA3 intronic SNP, rs139051 A>G was identified as a significant eQTL (p = 6.6×10^-8) influencing PNPLA3 transcription, with the A allele significantly associated with increased PNPLA3 mRNA. An electrophoresis mobility shift assay further demonstrated that the A allele has enhanced affinity to nuclear proteins than the G allele. The impact of this eQTL on NAFLD risk was further tested in three independent populations. We found that rs139051 did not independently affect the NAFLD risk, whilst rs738409 did not significantly modulate PNPLA3 transcription but was associated with NAFLD risk. The A-G haplotype associated with higher transcription of the disease-risk rs738409 G allele conferred similar risk for NAFLD compared to the G-G haplotype that possesses a lower transcription level. Our study suggests that the pathogenic role of PNPLA3^148M in NAFLD is independent of the gene transcription in humans, which may be attributed to the high endogenous transcription level of PNPLA3 gene in human livers.

Examination of previously identified associations within the Genetic Analysis Workshop 19 data

We investigate the possible replication of “known” associated single-nucleotide polymorphisms (SNPs) with blood pressure and expression phenotypes. Previous studies have provided a list of 95 SNPs thought to be associated with blood pressure phenotypes, of which 44 were present in the Genetic Analysis Workshop 19 (GAW19) family-imputed genome-wide association studies (GWAS) data and 4 in the GAW19 unrelateds sequence data. Using only the real (not simulated) GAW19 data, we show through the use of statistical tests that account for family relatedness, using FaST-LMM (Factored Spectrally Transformed Linear Mixed Model), that none of our candidate SNPs yields a significant p value. Furthermore, a study of epistasis, aiming to detect statistical interactions between loci with respect to their association with transcription levels has provided a list of 30 associated interacting SNP pairs, of which 13 are present in the GAW19 family GWAS and expression data. We show for this set of results, using the program GEMMA (genome-wide efficient mixed-model analysis) to account for family relatedness, that there is evidence of replication within the real GAW19 data. Two individual SNP pairs reach significance, and the set of remaining results give a combined p value of 0.017 that at least 1 of these remaining SNP pairs interacts to influence an expression phenotype.

Methylation quantitative trait locus analysis of osteoarthritis links epigenetics with genetic risk

Osteoarthritis (OA) is a common, painful and debilitating disease of articulating joints resulting from the age-associated loss of cartilage. Well-powered genetic studies have identified a number of DNA polymorphisms that are associated with OA susceptibility. Like most complex trait loci, these OA loci are thought to influence disease susceptibility through the regulation of gene expression, so-called expression quantitative loci, or eQTLs. One mechanism through which eQTLs act is epigenetic, by modulating DNA methylation. In such cases, there are quantitative differences in DNA methylation between the two alleles of the causal polymorphism, with the association signal referred to as a methylation quantitative trait locus, or meQTL. In this study, we aimed to investigate whether the OA susceptibility loci identified to date are functioning as meQTLs by integrating genotype data with whole genome methylation data of cartilage DNA. We investigated potential genotype-methylation correlations within a 1.0-1.5 Mb region surrounding each of 16 OA-associated single-nucleotide polymorphisms (SNPs) in 99 cartilage samples and identified four that function as meQTLs. Three of these replicated in an additional cohort of up to 62 OA patients. These observations suggest that OA susceptibility loci regulate the level of DNA methylation in cis and provide a mechanistic explanation as to how these loci impact upon OA susceptibility, further increasing our understanding of the role of genetics and epigenetics in this common disease.

Linkage Analysis in Autoimmune Addison's Disease: NFATC1 as a Potential Novel Susceptibility Locus

Background

Autoimmune Addison’s disease (AAD) is a rare, highly heritable autoimmune endocrinopathy. It is possible that there may be some highly penetrant variants which confer disease susceptibility that have yet to be discovered.

Methods

DNA samples from 23 multiplex AAD pedigrees from the UK and Norway (50 cases, 67 controls) were genotyped on the Affymetrix SNP 6.0 array. Linkage analysis was performed using Merlin. EMMAX was used to carry out a genome-wide association analysis comparing the familial AAD cases to 2706 UK WTCCC controls. To explore some of the linkage findings further, a replication study was performed by genotyping 64 SNPs in two of the four linked regions (chromosomes 7 and 18), on the Sequenom iPlex platform in three European AAD case-control cohorts (1097 cases, 1117 controls). The data were analysed using a meta-analysis approach.

Results

In a parametric analysis, applying a rare dominant model, loci on chromosomes 7, 9 and 18 had LOD scores >2.8. In a non-parametric analysis, a locus corresponding to the HLA region on chromosome 6, known to be associated with AAD, had a LOD score >3.0. In the genome-wide association analysis, a SNP cluster on chromosome 2 and a pair of SNPs on chromosome 6 were associated with AAD (P <5x10^-7). A meta-analysis of the replication study data demonstrated that three chromosome 18 SNPs were associated with AAD, including a non-synonymous variant in the NFATC1 gene.

Conclusion

This linkage study has implicated a number of novel chromosomal regions in the pathogenesis of AAD in multiplex AAD families and adds further support to the role of HLA in AAD. The genome-wide association analysis has also identified a region of interest on chromosome 2. A replication study has demonstrated that the NFATC1 gene is worthy of future investigation, however each of the regions identified require further, systematic analysis.

The heritability of abortion in pedigree Charollais flocks

Foetal death, or abortion at term, in sheep is of major significance to the livestock industry, accounting for more than £24million lost per annum. We have investigated whether there is a genetic component to abortion within two flocks of pedigree Charollais sheep, one followed from 1989 to 2006, the other from 1992 to 2006. Abortion occurred at a rate of 5.74-8.78% per annum against a total mortality rate of 14-24%. By model covariate analysis we have shown that 15.5% aborting ewes went on to have one or more abortions and that this risk increased with parity (p=0.006). Heritability estimates were approximately 0.08 as calculated by SOLAR, pedigreemm and ASReml3, with sire and dam components of 0.046 and 0.048, respectively. Where the lamb was aborted, heritability estimates were highly variable according to the method employed, 0.046-0.378, with sex of the lamb being a significant covariate. This variability indicated one or more underlying, significant factors that were not measured in these analyses, potentially including infectious agents that may be involved. Nevertheless, the ASReml3 estimate (0.179) resolved to 0.074 variance attributable to the sire and 0.092 attributable to the dam, which, while not significant, was suggestive that genetic variants passed by the dam to the lamb may be of more weight than that from the sire in determining whether a lamb will abort.

Genotype/phenotype correlations in complement factor H deficiency arising from uniparental isodisomy

We report a male infant who presented at 8 months of age with atypical hemolytic uremic syndrome (aHUS) responsive to plasma therapy. Investigation showed him to have complement factor H (CFH) deficiency associated with a homozygous CFH mutation (c.2880delT [p.Phe960fs]). Mutation screening of the child's parents revealed that the father was heterozygous for this change but that it was not present in his mother. Chromosome 1 uniparental isodisomy of paternal origin was confirmed by genotyping chromosome 1 SNPs. CD46 SNP genotyping was undertaken in this individual and another patient with CFH deficiency associated with chromosome 1 uniparental isodisomy. This showed a homozygous aHUS risk haplotype (CD46GGAAC) in the patient with aHUS and a homozygous glomerulonephritis risk haplotype (CD46AAGGT) in the patient with endocapillary glomerulonephritis. We also showed that FHL-1 (factor H-like protein 1) was present in the patient with aHUS and absent in the patient with glomerulonephritis. This study emphasizes that modifiers such as CD46 and FHL-1 may determine the kidney phenotype of patients who present with homozygous CFH deficiency.

A new genome scan for primary nonsyndromic vesicoureteric reflux emphasizes high genetic heterogeneity and shows linkage and association with various genes already implicated in urinary tract development

Primary vesicoureteric reflux (VUR), the retrograde flow of urine from the bladder toward the kidneys, results from a developmental anomaly of the vesicoureteric valve mechanism, and is often associated with other urinary tract anomalies. It is the most common urological problem in children, with an estimated prevalence of 1–2%, and is a major cause of hypertension in childhood and of renal failure in childhood or adult life. We present the results of a genetic linkage and association scan using 900,000 markers. Our linkage results show a large number of suggestive linkage peaks, with different results in two groups of families, suggesting that VUR is even more genetically heterogeneous than previously imagined. The only marker achieving P < 0.02 for linkage in both groups of families is 270 kb from EMX2. In three sibships, we found recessive linkage to KHDRBS3, previously reported in a Somali family. In another family we discovered sex-reversal associated with VUR, implicating PRKX, for which there was weak support for dominant linkage in the overall data set. Several other candidate genes are suggested by our linkage or association results, and four of our linkage peaks are within copy-number variants recently found to be associated with renal hypodysplasia. Undoubtedly there are many genes related to VUR. Our study gives support to some loci suggested by earlier studies as well as suggesting new ones, and provides numerous indications for further investigations.

Genome-wide association study identifies loci on 12q24 and 13q32 associated with tetralogy of Fallot

We conducted a genome-wide association study to search for risk alleles associated with Tetralogy of Fallot (TOF), using a northern European discovery set of 835 cases and 5159 controls. A region on chromosome 12q24 was associated (P = 1.4 × 10(-7)) and replicated convincingly (P = 3.9 × 10(-5)) in 798 cases and 2931 controls [per allele odds ratio (OR) = 1.27 in replication cohort, P = 7.7 × 10(-11) in combined populations]. Single nucleotide polymorphisms in the glypican 5 gene on chromosome 13q32 were also associated (P = 1.7 × 10(-7)) and replicated convincingly (P = 1.2 × 10(-5)) in 789 cases and 2927 controls (per allele OR = 1.31 in replication cohort, P = 3.03 × 10(-11) in combined populations). Four additional regions on chromosomes 10, 15 and 16 showed suggestive association accompanied by nominal replication. This study, the first genome-wide association study of a congenital heart malformation phenotype, provides evidence that common genetic variation influences the risk of TOF.

Contribution of global rare copy-number variants to the risk of sporadic congenital heart disease

Previous studies have shown that copy-number variants (CNVs) contribute to the risk of complex developmental phenotypes. However, the contribution of global CNV burden to the risk of sporadic congenital heart disease (CHD) remains incompletely defined. We generated genome-wide CNV data by using Illumina 660W-Quad SNP arrays in 2,256 individuals with CHD, 283 trio CHD-affected families, and 1,538 controls. We found association of rare genic deletions with CHD risk (odds ratio [OR] = 1.8, p = 0.0008). Rare deletions in study participants with CHD had higher gene content (p = 0.001) with higher haploinsufficiency scores (p = 0.03) than they did in controls, and they were enriched with Wnt-signaling genes (p = 1 × 10(-5)). Recurrent 15q11.2 deletions were associated with CHD risk (OR = 8.2, p = 0.02). Rare de novo CNVs were observed in ~5% of CHD trios; 10 out of 11 occurred on the paternally transmitted chromosome (p = 0.01). Some of the rare de novo CNVs spanned genes known to be involved in heart development (e.g., HAND2 and GJA5). Rare genic deletions contribute ~4% of the population-attributable risk of sporadic CHD. Second to previously described CNVs at 1q21.1, deletions at 15q11.2 and those implicating Wnt signaling are the most significant contributors to the risk of sporadic CHD. Rare de novo CNVs identified in CHD trios exhibit paternal origin bias.

Phenotype-specific effect of chromosome 1q21.1 rearrangements and GJA5 duplications in 2436 congenital heart disease patients and 6760 controls

Recurrent rearrangements of chromosome 1q21.1 that occur via non-allelic homologous recombination have been associated with variable phenotypes exhibiting incomplete penetrance, including congenital heart disease (CHD). However, the gene or genes within the ~1 Mb critical region responsible for each of the associated phenotypes remains unknown. We examined the 1q21.1 locus in 948 patients with tetralogy of Fallot (TOF), 1488 patients with other forms of CHD and 6760 ethnically matched controls using single nucleotide polymorphism genotyping arrays (Illumina 660W and Affymetrix 6.0) and multiplex ligation-dependent probe amplification. We found that duplication of 1q21.1 was more common in cases of TOF than in controls [odds ratio (OR) 30.9, 95% confidence interval (CI) 8.9-107.6); P = 2.2 × 10(-7)], but deletion was not. In contrast, deletion of 1q21.1 was more common in cases of non-TOF CHD than in controls [OR 5.5 (95% CI 1.4-22.0); P = 0.04] while duplication was not. We also detected rare (n = 3) 100-200 kb duplications within the critical region of 1q21.1 in cases of TOF. These small duplications encompassed a single gene in common, GJA5, and were enriched in cases of TOF in comparison to controls [OR = 10.7 (95% CI 1.8-64.3), P = 0.01]. These findings show that duplication and deletion at chromosome 1q21.1 exhibit a degree of phenotypic specificity in CHD, and implicate GJA5 as the gene responsible for the CHD phenotypes observed with copy number imbalances at this locus.

Whole-genome linkage and association scan in primary, nonsyndromic vesicoureteric reflux

Primary vesicoureteric reflux accounts for approximately 10% of kidney failure requiring dialysis or transplantation, and sibling studies suggest a large genetic component. Here, we report a whole-genome linkage and association scan in primary, nonsyndromic vesicoureteric reflux and reflux nephropathy. We used linkage and family-based association approaches to analyze 320 white families (661 affected individuals, generally from families with two affected siblings) from two populations (United Kingdom and Slovenian). We found modest evidence of linkage but no clear overlap with previous studies. We tested for but did not detect association with six candidate genes (AGTR2, HNF1B, PAX2, RET, ROBO2, and UPK3A). Family-based analysis detected associations with one single-nucleotide polymorphism (SNP) in the UK families, with three SNPs in the Slovenian families, and with three SNPs in the combined families. A case-control analysis detected associations with three additional SNPs. The results of this study, which is the largest to date investigating the genetics of reflux, suggest that major loci may not exist for this common renal tract malformation within European populations.