Loss of GPR75 protects against non-alcoholic fatty liver disease and body fat accumulation

Approximately 1 in 4 people worldwide have non-alcoholic fatty liver disease (NAFLD); however, there are currently no medications to treat this condition. This study investigated the role of adiposity-associated orphan G protein-coupled receptor 75 (GPR75) in liver lipid accumulation. We profiled Gpr75 expression and report that it is most abundant in the brain. Next, we generated the first single-cell-level analysis of Gpr75 and identified a subpopulation co-expressed with key appetite-regulating hypothalamic neurons. CRISPR-Cas9-deleted Gpr75 mice fed a palatable western diet high in fat adjusted caloric intake to remain in energy balance, thereby preventing NAFLD. Consistent with mouse results, analysis of whole-exome sequencing data from 428,719 individuals (UK Biobank) revealed that variants in GPR75 are associated with a reduced likelihood of hepatic steatosis. Here, we provide a significant advance in understanding of the expression and function of GPR75, demonstrating that it is a promising pharmaceutical target for NAFLD treatment.

Investigation of Genetic Alterations Associated With Interval Breast Cancer

Importance

Breast cancers (BCs) diagnosed between 2 screening examinations are called interval cancers (ICs), and they have worse clinicopathological characteristics and poorer prognosis than screen-detected cancers (SDCs). However, the association of rare germline genetic variants with IC have not been studied.

Objective

To evaluate whether rare germline deleterious protein-truncating variants (PTVs) can be applied to discriminate between IC and SDC while considering mammographic density.

Design, Setting, and Participants

This population-based genetic association study was based on women aged 40 to 76 years who were attending mammographic screening in Sweden. All women with a diagnosis of BC between January 2001 and January 2016 were included, together with age-matched controls. Patients with BC were followed up for survival until 2021. Statistical analysis was performed from September 2021 to December 2022.

Exposure

Germline PTVs in 34 BC susceptibility genes as analyzed by targeted sequencing.

Main Outcomes and Measures

Odds ratios (ORs) were used to compare IC with SDC using logistic regression. Hazard ratios were used to investigate BC-specific survival using Cox regression.

Results

All 4121 patients with BC (IC, n = 1229; SDC, n = 2892) were female, with a mean (SD) age of 55.5 (7.1) years. There were 5631 age-matched controls. The PTVs of the ATM, BRCA1, BRCA2, CHEK2, and PALB2 genes were more common in patients with IC compared with SDC (OR, 1.48; 95% CI, 1.06-2.05). This association was primarily influenced by BRCA1/2 and PALB2 variants. A family history of BC together with PTVs of any of these genes synergistically increased the probability of receiving a diagnosis of IC rather than SDC (OR, 3.95; 95% CI, 1.97-7.92). Furthermore, 10-year BC-specific survival revealed that if a patient received a diagnosis of an IC, carriers of PTVs in any of these 5 genes had significantly worse survival compared with patients not carrying any of them (hazard ratio, 2.04; 95% CI, 1.06-3.92). All of these associations were further pronounced in a subset of patients with IC who had a low mammographic density at prior screening examination.

Conclusions and Relevance

The results of this study may be helpful in future optimizations of screening programs that aim to lower mortality as well as the clinical treatment of patients with BC.

The impact of circulating protein levels identified by affinity proteomics on short-term, overall breast cancer risk

OBJECTIVE

Current breast cancer risk prediction scores and algorithms can potentially be further improved by including molecular markers. To this end, we studied the association of circulating plasma proteins using Proximity Extension Assay (PEA) with incident breast cancer risk.

SUBJECTS

In this study, we included 1577 women participating in the prospective KARMA mammographic screening cohort.

RESULTS

In a targeted panel of 164 proteins, we found 8 candidates nominally significantly associated with short-term breast cancer risk (P < 0.05). Similarly, in an exploratory panel consisting of 2204 proteins, 115 were found nominally significantly associated (P < 0.05). However, none of the identified protein levels remained significant after adjustment for multiple testing. This lack of statistically significant findings was not due to limited power, but attributable to the small effect sizes observed even for nominally significant proteins. Similarly, adding plasma protein levels to established risk factors did not improve breast cancer risk prediction accuracy.

CONCLUSIONS

Our results indicate that the levels of the studied plasma proteins captured by the PEA method are unlikely to offer additional benefits for risk prediction of short-term overall breast cancer risk but could provide interesting insights into the biological basis of breast cancer in the future.

Cancer risks among first-degree relatives of women with a genetic predisposition to breast cancer

BACKGROUND

Associations between germline alterations in women and cancer risks among their relatives are largely unknown.

METHODS

We used women from two Swedish cohorts (KARMA and pKARMA), including 28,362 women with genotyping data and 13,226 with sequencing data. Using Swedish Multi-Generation Register, we linked these women to 133,389 first-degree relatives. Associations between protein-truncating variants (PTVs) in 8 risk genes and breast cancer polygenic risk score (PRS) in index women and cancer risks among their relatives were modeled via Cox regression.

RESULTS

Female relatives of index women who were PTV carriers in any of the 8 risk genes had an increased breast cancer risk compared to those of non-carriers (HR1.85, 95% CI: 1.52-2.27), with the strongest association found for PTVs in BRCA1/2. These relatives had a statistically higher risk of early-onset than late-onset breast cancer (P = .001). Elevated breast cancer risk was also observed in female relatives of index women with higher PRS (HR per SD: 1.28, 95% CI: 1.23-1.32). The estimated lifetime risk was 22.3% for female relatives of PTV carriers and 14.4% for those related to women in the top PRS quartile. Moreover, relatives of index women with PTV presence (HR: 1.30, 95% CI: 1.06-1.59) or higher PRS (HR per SD: 1.04, 95% CI: 1.01-1.07) were also at higher risk of non-breast-HBOC cancers, including prostate, ovarian, pancreatic cancer, and melanoma.

CONCLUSIONS

Both PTVs of risk genes and higher PRS in index women are associated with an increased risk of breast and other HBOC-related cancers among relatives.

The glucocorticoid receptor as a master regulator of the Müller cell response to diabetic conditions in mice

Diabetic retinopathy (DR) is considered a primarily microvascular complication of diabetes. Müller glia cells are at the centre of the retinal neurovascular unit and play a critical role in DR. We therefore investigated Müller cell-specific signalling pathways that are altered in DR to identify novel targets for gene therapy. Using a multi-omics approach on purified Müller cells from diabetic db/db mice, we found the mRNA and protein expression of the glucocorticoid receptor (GR) to be significantly decreased, while its target gene cluster was down-regulated. Further, oPOSSUM TF analysis and ATAC- sequencing identified the GR as a master regulator of Müller cell response to diabetic conditions. Cortisol not only increased GR phosphorylation. It also induced changes in the expression of known GR target genes in retinal explants. Finally, retinal functionality was improved by AAV-mediated overexpression of GR in Müller cells. Our study demonstrates an important role of the glial GR in DR and implies that therapeutic approaches targeting this signalling pathway should be aimed at increasing GR expression rather than the addition of more ligand.

An ancient polymorphic regulatory region within the BDNF gene associated with obesity modulates anxiety-like behaviour in mice and humans

Obesity and anxiety are morbidities notable for their increased impact on society during the recent COVID-19 pandemic. Understanding the mechanisms governing susceptibility to these conditions will increase our quality of life and resilience to future pandemics. In the current study, we explored the function of a highly conserved regulatory region (BE5.1) within the BDNF gene that harbours a polymorphism strongly associated with obesity (rs10767664; p = 4.69 × 10-26). Analysis in primary cells suggested that the major T-allele of BE5.1 was an enhancer, whereas the obesity-associated A-allele was not. However, CRISPR/CAS9 deletion of BE5.1 from the mouse genome (BE5.1KO) produced no significant effect on the expression of BDNF transcripts in the hypothalamus, no change in weight gain after 28 days and only a marginally significant increase in food intake. Nevertheless, transcripts were significantly increased in the amygdala of female mice and elevated zero maze and marble-burying tests demonstrated a significant increase in anxiety-like behaviour that could be reversed by diazepam. Consistent with these observations, human GWAS cohort analysis demonstrated a significant association between rs10767664 and anxiousness in human populations. Intriguingly, interrogation of the human GTEx eQTL database demonstrated no effect on BDNF mRNA levels associated with rs10767664 but a highly significant effect on BDNF-antisense (BDNF-AS) gene expression and splicing. The subsequent observation that deletion of BE5.1 also significantly reduced BDNF-AS expression in mice suggests a novel mechanism in the regulation of BDNF expression common to mice and humans, which contributes to the modulation of mood and anxiety in both species.

Evaluation of circulating plasma proteins in breast cancer using Mendelian randomisation

Biomarkers for early detection of breast cancer may complement population screening approaches to enable earlier and more precise treatment. The blood proteome is an important source for biomarker discovery but so far, few proteins have been identified with breast cancer risk. Here, we measure 2929 unique proteins in plasma from 598 women selected from the Karolinska Mammography Project to explore the association between protein levels, clinical characteristics, and gene variants, and to identify proteins with a causal role in breast cancer. We present 812 cis-acting protein quantitative trait loci for 737 proteins which are used as instruments in Mendelian randomisation analyses of breast cancer risk. Of those, we present five proteins (CD160, DNPH1, LAYN, LRRC37A2 and TLR1) that show a potential causal role in breast cancer risk with confirmatory results in independent cohorts. Our study suggests that these proteins should be further explored as biomarkers and potential drug targets in breast cancer.

Prediction of breast cancer risk for sisters of women attending screening

BACKGROUND

Risk assessment is important for breast cancer prevention and early detection. We aimed to examine whether common risk factors, mammographic features, and breast cancer risk prediction scores of a woman were associated with breast cancer risk for her sisters.

METHODS

We included 53 051 women from the Karolinska Mammography Project for Risk Prediction of Breast Cancer (KARMA) study. Established risk factors were derived using self-reported questionnaires, mammograms, and single nucleotide polymorphism genotyping. Using the Swedish Multi-Generation Register, we identified 32 198 sisters of the KARMA women (including 5352 KARMA participants and 26 846 nonparticipants). Cox models were used to estimate the hazard ratios of breast cancer for both women and their sisters, respectively.

RESULTS

A higher breast cancer polygenic risk score, a history of benign breast disease, and higher breast density in women were associated with an increased risk of breast cancer for both women and their sisters. No statistically significant association was observed between breast microcalcifications and masses in women and breast cancer risk for their sisters. Furthermore, higher breast cancer risk scores in women were associated with an increased risk of breast cancer for their sisters. Specifically, the hazard ratios for breast cancer per 1 standard deviation increase in age-adjusted KARMA, Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm (BOADICEA), and Tyrer-Cuzick risk scores were 1.16 (95% confidence interval [CI] = 1.07 to 1.27), 1.23 (95% CI = 1.12 to 1.35), and 1.21 (95% CI = 1.11 to 1.32), respectively.

CONCLUSION

A woman's breast cancer risk factors are associated with her sister's breast cancer risk. However, the clinical utility of these findings requires further investigation.

Risk of estrogen receptor-specific breast cancer by family history of estrogen receptor subtypes and other cancers

BACKGROUND

The extent to which the risk of estrogen receptor (ER)-specific breast cancer is associated with ER status of breast cancer and other cancers among first-degree relatives is unclear.

METHODS

This population-based cohort included 464 707 cancer-free women in Stockholm, Sweden, during 1978-2019. For ER-negative and ER-positive breast cancers, we estimated hazard ratios (HRs) associated with ER status of female first-degree relatives with breast cancer and of other cancers in all first-degree relatives. Associations between ER-negative and ER-positive status by family cancer history were estimated using logistic regression in a case-only design.

RESULTS

Women with familial ER-positive breast cancer had 1.87 times (95% confidence interval [CI] = 1.77 to 1.97) higher risk of ER-positive subtype, whereas the corresponding hazard ratio for ER-negative was 2.54 (95% CI = 2.08 to 3.10) when having familial ER-negative breast cancer. The risk increased with an increasing number of female first-degree relatives having concordant subtypes and younger age at diagnosis (Ptrend <.001 for both). Nonbreast cancers among first-degree relatives were associated with both ER-positive (HR = 1.14, 95% CI = 1.10 to 1.17) and ER-negative (HR = 1.08, 95% CI = 1.01 to 1.16) breast cancers. Compared with women with ER-positive breast cancer, women with ER-negative breast cancer were more likely to have family history of liver (odds ratio [OR] = 1.33, 95% CI = 1.05 to 1.67), ovary (OR = 1.28, 95% CI = 1.01 to 1.61), and testicle cancer (OR = 1.79, 95% CI = 1.01 to 3.16) but less likely to have family history of endometrial cancer (OR = 0.77, 95% CI = 0.60 to 1.00) and leukemia (OR = 0.72, 95% CI = 0.56 to 0.91).

CONCLUSIONS

Risk of ER-specific breast cancer differs according to ER status of female first-degree relatives with breast cancer and some other cancers of first-degree relatives. This family history information should be considered in the individual risk prediction for ER subtypes.

A genome-wide gene-environment interaction study of breast cancer risk for women of European ancestry

BACKGROUND

Genome-wide studies of gene-environment interactions (G×E) may identify variants associated with disease risk in conjunction with lifestyle/environmental exposures. We conducted a genome-wide G×E analysis of ~ 7.6 million common variants and seven lifestyle/environmental risk factors for breast cancer risk overall and for estrogen receptor positive (ER +) breast cancer.

METHODS

Analyses were conducted using 72,285 breast cancer cases and 80,354 controls of European ancestry from the Breast Cancer Association Consortium. Gene-environment interactions were evaluated using standard unconditional logistic regression models and likelihood ratio tests for breast cancer risk overall and for ER + breast cancer. Bayesian False Discovery Probability was employed to assess the noteworthiness of each SNP-risk factor pairs.

RESULTS

Assuming a 1 × 10-5 prior probability of a true association for each SNP-risk factor pairs and a Bayesian False Discovery Probability < 15%, we identified two independent SNP-risk factor pairs: rs80018847(9p13)-LINGO2 and adult height in association with overall breast cancer risk (ORint = 0.94, 95% CI 0.92-0.96), and rs4770552(13q12)-SPATA13 and age at menarche for ER + breast cancer risk (ORint = 0.91, 95% CI 0.88-0.94).

CONCLUSIONS

Overall, the contribution of G×E interactions to the heritability of breast cancer is very small. At the population level, multiplicative G×E interactions do not make an important contribution to risk prediction in breast cancer.

Association of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival

BACKGROUND

Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC-specific survival (BCSS) compared to non-carriers.

AIM

To assessed the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS.

METHODS

Analyses were based on 82,701 women diagnosed with a first primary invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations with treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death.

RESULTS

There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status. The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR (95% CI): 0.66 (0.55-0.78)]. No association was observed with radiotherapy. Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR (95% CI): 1.30 (1.09-1.56)].

CONCLUSION

Systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk.

Spectrum and Frequency of Germline FANCM Protein-Truncating Variants in 44,803 European Female Breast Cancer Cases

FANCM germline protein truncating variants (PTVs) are moderate-risk factors for ER-negative breast cancer. We previously described the spectrum of FANCM PTVs in 114 European breast cancer cases. In the present, larger cohort, we report the spectrum and frequency of four common and 62 rare FANCM PTVs found in 274 carriers detected among 44,803 breast cancer cases. We confirmed that p.Gln1701* was the most common PTV in Northern Europe with lower frequencies in Southern Europe. In contrast, p.Gly1906Alafs*12 was the most common PTV in Southern Europe with decreasing frequencies in Central and Northern Europe. We verified that p.Arg658* was prevalent in Central Europe and had highest frequencies in Eastern Europe. We also confirmed that the fourth most common PTV, p.Gln498Thrfs*7, might be a founder variant from Lithuania. Based on the frequency distribution of the carriers of rare PTVs, we showed that the FANCM PTVs spectra in Southwestern and Central Europe were much more heterogeneous than those from Northeastern Europe. These findings will inform the development of more efficient FANCM genetic testing strategies for breast cancer cases from specific European populations.

Bipolar haemostatic forceps versus standard therapy by haemoclip + / - epinephrine injection as initial endoscopic treatment in active non-variceal upper GI bleeding: study protocol for a prospective, randomized multicentre trial (BeBop-Trial)

BACKGROUND

Patients with active nonvariceal upper gastrointestinal bleeding (NVUGIB) usually require urgent endoscopic treatment. Standard therapy (ST) using haemoclip + / - epinephrine injection is not always successful. Bipolar haemostatic forceps (HemoStat/Pentax®) are an approved medical device for the treatment of gastrointestinal bleeding. However, their use as a primary endoscopic treatment for active NVUGIB has not yet been proven in a randomized prospective study.

METHODS

This is a prospective, randomized, multicentre superiority trial (n ≥ 5). Patients with active NVUGIB will be randomized (1:1) to ST and to experimental therapy (ET) by application of bipolar haemostatic forceps. In the case of failed initial treatment within 15 min, crossover treatment will be attempted first. Rescue treatment (e.g. via over-the-scope-clip) will then be allowed after 30 min. All patients will also receive standard therapy with proton pump inhibitors. Forty-five patients per treatment arm are required to demonstrate an absolute difference of 25.4% with a power of 80% and a significance level of 0.05.

DISCUSSION

The hypothesis of the study is that bipolar haemostatic forceps are superior to ST in terms of successful primary haemostasis and the absence of recurrent bleeding within 30 days (combined endpoint). The 1:1 randomization is also ethically justifiable for this study, as both procedures are approved for the intervention in question. To further increase the safety of the patients in the study, crossover treatment and rescue treatment are planned. The prospective design seems feasible in a reasonable time frame (recruitment period of 12 months), as nonvariceal upper gastrointestinal bleeding is common. Anticoagulants and/or antiplatelet drugs could be an important confounding factor in the statistical analysis that needs to be taken into account and calculated if necessary. In conclusion, this randomized, prospective, multicentre study could make an important contribution to answering the question of whether bipolar haemostatic forceps could be the first-line therapy in the endoscopic treatment of stage Forrest I a + b NVUGIB.

TRIAL REGISTRATION

ClinicalTrials.gov NCT05353062. Registered on April 30 2022.

Association of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival

Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC-specific survival (BCSS) compared to non-carriers. We aimed to assess the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS. Analyses were based on 82,701 women diagnosed with invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations of treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR(95%CI): 0.66 (0.55-0.78)]. No association was observed with radiotherapy. Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR(95%CI) :1.30 (1.09-1.56)]. In conclusion, systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk. (Main MS: 3201 words).

FANCM missense variants and breast cancer risk: a case-control association study of 75,156 European women

Evidence from literature, including the BRIDGES study, indicates that germline protein truncating variants (PTVs) in FANCM confer moderately increased risk of ER-negative and triple-negative breast cancer (TNBC), especially for women with a family history of the disease. Association between FANCM missense variants (MVs) and breast cancer risk has been postulated. In this study, we further used the BRIDGES study to test 689 FANCM MVs for association with breast cancer risk, overall and in ER-negative and TNBC subtypes, in 39,885 cases (7566 selected for family history) and 35,271 controls of European ancestry. Sixteen common MVs were tested individually; the remaining rare 673 MVs were tested by burden analyses considering their position and pathogenicity score. We also conducted a meta-analysis of our results and those from published studies. We did not find evidence for association for any of the 16 variants individually tested. The rare MVs were significantly associated with increased risk of ER-negative breast cancer by burden analysis comparing familial cases to controls (OR = 1.48; 95% CI 1.07-2.04; P = 0.017). Higher ORs were found for the subgroup of MVs located in functional domains or predicted to be pathogenic. The meta-analysis indicated that FANCM MVs overall are associated with breast cancer risk (OR = 1.22; 95% CI 1.08-1.38; P = 0.002). Our results support the definition from previous analyses of FANCM as a moderate-risk breast cancer gene and provide evidence that FANCM MVs could be low/moderate risk factors for ER-negative and TNBC subtypes. Further genetic and functional analyses are necessary to clarify better the increased risks due to FANCM MVs.

Aggregation tests identify new gene associations with breast cancer in populations with diverse ancestry

BACKGROUND

Low-frequency variants play an important role in breast cancer (BC) susceptibility. Gene-based methods can increase power by combining multiple variants in the same gene and help identify target genes.

METHODS

We evaluated the potential of gene-based aggregation in the Breast Cancer Association Consortium cohorts including 83,471 cases and 59,199 controls. Low-frequency variants were aggregated for individual genes' coding and regulatory regions. Association results in European ancestry samples were compared to single-marker association results in the same cohort. Gene-based associations were also combined in meta-analysis across individuals with European, Asian, African, and Latin American and Hispanic ancestry.

RESULTS

In European ancestry samples, 14 genes were significantly associated (q < 0.05) with BC. Of those, two genes, FMNL3 (P = 6.11 × 10-6) and AC058822.1 (P = 1.47 × 10-4), represent new associations. High FMNL3 expression has previously been linked to poor prognosis in several other cancers. Meta-analysis of samples with diverse ancestry discovered further associations including established candidate genes ESR1 and CBLB. Furthermore, literature review and database query found further support for a biologically plausible link with cancer for genes CBLB, FMNL3, FGFR2, LSP1, MAP3K1, and SRGAP2C.

CONCLUSIONS

Using extended gene-based aggregation tests including coding and regulatory variation, we report identification of plausible target genes for previously identified single-marker associations with BC as well as the discovery of novel genes implicated in BC development. Including multi ancestral cohorts in this study enabled the identification of otherwise missed disease associations as ESR1 (P = 1.31 × 10-5), demonstrating the importance of diversifying study cohorts.

Low leukocyte mitochondrial DNA abundance drives atherosclerotic cardiovascular diseases: cohort and Mendelian randomization study

Aim

Mitochondrial DNA dysfunction has been implicated in the pathogenesis of cardiovascular diseases. We aimed to investigate the associations of leukocyte mitochondrial DNA (mtDNA) abundance, as a proxy of mitochondrial copy numbers and function, and coronary artery disease (CAD) and heart failure (HF) in a cohort study and approximate the causal nature of these relationships using Mendelian randomization (MR) in genetic studies.

Methods and Results

Multivariable-adjusted Cox regression analyses were conducted in 273,619 unrelated participants of European ancestry from UK Biobank (UKB). For genetic studies, we first performed MR analyses with individual-level data from the UKB using a weighted genetic risk score (GRS); two-sample MR analyses were subsequently performed using summary-level data from the publicly available three consortia/biobank for CAD and two for HF. MR analyses were performed per database separately and results were subsequently meta-analyzed using fixed-effects models. During a median follow-up of 11.8 years, cox regression restricted cubic spline analyses showed associations between lower mtDNA abundance and higher risk of CAD and HF. Hazard ratios for participants in the lowest quintile of mtDNA abundance compared with those in the highest quintile were 1.08 (95% confidence interval: 1.03, 1.14) and 1.15 (1.05, 1.24) for CAD and HF. Genetically, no evidence was observed for a possible non-linear causal effect using individual-level weighted genetic risk scores calculated in the UKB on the study outcomes; the pooled odds ratios from two-sample MR of genetically predicted per one-SD decrease in mtDNA abundance were 1.09 (1.03, 1.16) for CAD and 0.99 (0.92, 1.08) for HF, respectively.

Conclusion

Our findings support a possible causal role of lower leukocyte mtDNA abundance in higher CAD risk, but not in HF.

Systems genomics in age-related macular degeneration

Genomic studies in age-related macular degeneration (AMD) have identified genetic variants that account for the majority of AMD risk. An important next step is to understand the functional consequences and downstream effects of the identified AMD-associated genetic variants. Instrumental for this next step are 'omics' technologies, which enable high-throughput characterization and quantification of biological molecules, and subsequent integration of genomics with these omics datasets, a field referred to as systems genomics. Single cell sequencing studies of the retina and choroid demonstrated that the majority of candidate AMD genes identified through genomic studies are expressed in non-neuronal cells, such as the retinal pigment epithelium (RPE), glia, myeloid and choroidal cells, highlighting that many different retinal and choroidal cell types contribute to the pathogenesis of AMD. Expression quantitative trait locus (eQTL) studies in retinal tissue have identified putative causal genes by demonstrating a genetic overlap between gene regulation and AMD risk. Linking genetic data to complement measurements in the systemic circulation has aided in understanding the effect of AMD-associated genetic variants in the complement system, and supports that protein QTL (pQTL) studies in plasma or serum samples may aid in understanding the effect of genetic variants and pinpointing causal genes in AMD. A recent epigenomic study fine-mapped AMD causal variants by determing regulatory regions in RPE cells differentiated from induced pluripotent stem cells (iPSC-RPE). Another approach that is being employed to pinpoint causal AMD genes is to produce synthetic DNA assemblons representing risk and protective haplotypes, which are then delivered to cellular or animal model systems. Pinpointing causal genes and understanding disease mechanisms is crucial for the next step towards clinical translation. Clinical trials targeting proteins encoded by the AMD-associated genomic loci C3, CFB, CFI, CFH, and ARMS2/HTRA1 are currently ongoing, and a phase III clinical trial for C3 inhibition recently showed a modest reduction of lesion growth in geographic atrophy. The EYERISK consortium recently developed a genetic test for AMD that allows genotyping of common and rare variants in AMD-associated genes. Polygenic risk scores (PRS) were applied to quantify AMD genetic risk, and may aid in predicting AMD progression. In conclusion, genomic studies represent a turning point in our exploration of AMD. The results of those studies now serve as a driving force for several clinical trials. Expanding to omics and systems genomics will further decipher function and causality from the associations that have been reported, and will enable the development of therapies that will lessen the burden of AMD.

Percutaneous transhepatic or endoscopic ultrasound-guided biliary drainage in malignant distal bile duct obstruction using a self-expanding metal stent: Study protocol for a prospective European multicenter trial (PUMa trial)

Background

Endoscopic ultrasound-guided biliary drainage (EUS-BD) was associated with better clinical success and a lower rate of adverse events (AEs) than fluoroscopy-guided percutaneous transhepatic biliary drainage (PTBD) in recent single center studies with mainly retrospective design and small case numbers (< 50). The aim of this prospective European multicenter study is to compare both drainage procedures using ultrasound-guidance and primary metal stent implantation in patients with malignant distal bile duct obstruction (PUMa Trial).

Methods

The study is designed as a non-randomized, controlled, parallel group, non-inferiority trial. Each of the 16 study centers performs the procedure with the best local expertise (PTBD or EUS-BD). In PTBD, bile duct access is performed by ultrasound guidance. EUS-BD is performed as an endoscopic ultrasound (EUS)-guided hepaticogastrostomy (EUS-HGS), EUS-guided choledochoduodenostomy (EUS-CDS) or EUS-guided antegrade stenting (EUS-AGS). Insertion of a metal stent is intended in both procedures in the first session. Primary end point is technical success. Secondary end points are clinical success, duration pf procedure, AEs graded by severity, length of hospital stay, re-intervention rate and survival within 6 months. The target case number is 212 patients (12 calculated dropouts included).

Discussion

This study might help to clarify whether PTBD is non-inferior to EUS-BD concerning technical success, and whether one of both interventions is superior in terms of efficacy and safety in one or more secondary endpoints. Randomization is not provided as both procedures are rarely used after failed endoscopic biliary drainage and study centers usually prefer one of both procedures that they can perform best.

Trial registration

ClinicalTrials.gov ID: NCT03546049 (22.05.2018).

Assessment of the bi-directional relationship between blood mitochondrial DNA copy number and type 2 diabetes mellitus: a multivariable-adjusted regression and Mendelian randomisation study

AIMS/HYPOTHESIS

Mitochondrial dysfunction, which can be approximated by blood mitochondrial DNA copy number (mtDNA-CN), has been implicated in the pathogenesis of type 2 diabetes mellitus. Thus far, however, insights from prospective cohort studies and Mendelian randomisation (MR) analyses on this relationship are limited. We assessed the association between blood mtDNA-CN and incident type 2 diabetes using multivariable-adjusted regression analyses, and the associations between blood mtDNA-CN and type 2 diabetes and BMI using bi-directional MR.

METHODS

Multivariable-adjusted Cox proportional hazard models were used to estimate the association between blood mtDNA-CN and incident type 2 diabetes in 285,967 unrelated European individuals from UK Biobank free of type 2 diabetes at baseline. Additionally, a cross-sectional analysis was performed to investigate the association between blood mtDNA-CN and BMI. We also assessed the potentially causal relationship between blood mtDNA-CN and type 2 diabetes (N=898,130 from DIAGRAM, N=215,654 from FinnGen) and BMI (N=681,275 from GIANT) using bi-directional two-sample MR.

RESULTS

During a median follow-up of 11.87 years, 15,111 participants developed type 2 diabetes. Participants with a higher level of blood mtDNA-CN are at lower risk of developing type 2 diabetes (HR 0.90 [95% CI 0.89, 0.92]). After additional adjustment for BMI and other confounders, these results attenuated moderately and remained present. The multivariable-adjusted cross-sectional analyses showed that higher blood mtDNA-CN was associated with lower BMI (-0.12 [95% CI -0.14, -0.10]) kg/m². In the bi-directional MR analyses, we found no evidence for causal associations between blood mtDNA-CN and type 2 diabetes, and blood mtDNA-CN and BMI in either direction.

CONCLUSIONS/INTERPRETATION

The results from the present study indicate that the observed association between low blood mtDNA-CN and higher risk of type 2 diabetes is likely not causal.

The association between mitochondrial DNA abundance and stroke: A combination of multivariable-adjusted survival and Mendelian randomization analyse

BACKGROUND AND AIMS

Mitochondrial dysfunction is associated with increased reactive oxygen species (ROS) that are thought to drive disease risk, including stroke. We investigated the association between mtDNA abundance, as a proxy measure of mitochondrial function, and incident stroke, using multivariable-adjusted survival and Mendelian Randomization (MR) analyses.

METHODS

Cox-proportional hazard model analyses were conducted to assess the association between mtDNA abundance, and incident ischemic and hemorrhagic stroke over a maximum of 14-year follow-up in European-ancestry participants from UK Biobank. MR was conducted using independent (R² < 0.001) lead variants for mtDNA abundance (p < 5 × 10^-8) as instrumental variables. Single-nucleotide polymorphism (SNP)-ischemic stroke associations were derived from three published open source European-ancestry results databases (cases/controls): MEGASTROKE (60,341/454,450), UK Biobank (2404/368,771) and FinnGen (10,551/202,223). MR was performed per study, and results were subsequently meta-analyzed.

RESULTS

In total, 288,572 unrelated participants (46% men) with mean (SD) age of 57 (8) years were included in the Cox-proportional hazard analyses. After correction for considered confounders (BMI, hypertension, cholesterol, T2D), no association was found between low versus high mtDNA abundance and ischemic (HR: 1.06 [95% CI: 0.95, 1.18]) or hemorrhagic (HR: 0.97 [95% CI: 0.82, 1.15]) stroke. However, in the MR analyses after removal of platelet count-associated SNPs, we found evidence for an association between genetically-influenced mtDNA abundance and ischemic stroke (odds ratio, 1.17; confidence interval, 1.03, 1.32).

CONCLUSIONS

Although the results from both multivariable-adjusted prospective and basis MR analyses did not show an association between low mtDNA and increased risk of ischemic stroke, in-depth MR sensitivity analyses may suggest evidence for a causal relationship.

Risk of heart disease following treatment for breast cancer: results from a population-based cohort study

Background: There is a rising concern about treatment-associated cardiotoxicities in breast cancer patients. This study aimed to determine the time- and treatment-specific incidence of arrhythmia, heart failure and ischemic heart disease in women diagnosed with breast cancer.

Methods: A register-based matched cohort study was conducted including 8015 breast cancer patients diagnosed from 2001-2008 in the Stockholm-Gotland region and followed-up until 2017. Time-dependent risks of arrhythmia, heart failure and ischemic heart disease in breast cancer patients were assessed using flexible parametric models as compared to matched controls from general population. Treatment-specific effects were estimated in breast cancer patients using Cox model.

Results: Time-dependent analyses revealed long-term increased risks of arrhythmia and heart failure following breast cancer diagnosis. Hazard ratios (HRs) within the first year of diagnosis were 2.14 (95% CI = 1.63-2.81) for arrhythmia and 2.71 (95% CI = 1.70-4.33) for heart failure. HR more than 10 years following diagnosis was 1.42 (95% CI = 1.21-1.67) for arrhythmia and 1.28 (95% CI = 1.03-1.59) for heart failure. The risk for ischemic heart disease was significantly increased only during the first year after diagnosis (HR=1.45, 95% CI = 1.03-2.04). Trastuzumab and anthracyclines were associated with increased risk of heart failure. Aromatase inhibitors, but not tamoxifen, were associated with risk of ischemic heart disease. No increased risk of heart disease was identified following loco-regional radiotherapy.

Conclusions: Administration of systemic adjuvant therapies appears to be associated with increased risks of heart disease. The risk estimates observed in this study may aid adjuvant therapy decision-making and patient counseling in oncology practices.

Funding: This work was supported by the Swedish Research Council [grant no: 2018-02547]; Swedish Cancer Society [grant no: CAN-19-0266] and FORTE [grant no: 2016-00081].

Interval breast cancer is associated with interferon immune response

BACKGROUND

The aggressive nature of breast cancers detected between planned mammographic screens, so-called interval cancers, remains elusive. Here, we aim to characterise underlying molecular features of interval cancer.

METHODS

From 672 patients with invasive breast cancer, we analysed gene expression differences between 90 'true' interval cancer cases (i.e. women with low-dense breasts defined as per cent mammographic density <25%) and 310 screen-detected tumours while accounting for PAM50 subtypes and thus overall tumour aggressiveness. We computed an interval cancer gene expression profile (IC-Gx) in a total of 2270 breast tumours (regardless of interval cancer status) and tested for association with expression-based immune subtypes in breast cancer. In addition, we investigated the contribution of inherited and somatic genetic variants in distinct features of interval cancer.

RESULTS

We identified 8331 genes nominally associated with interval cancer (P-value < 0.05, fold-change > 1.5). Gene set enrichment analysis showed immune-related pathways as key processes altered in interval cancer. Our IC-Gx, based on 47 genes with the strongest associations (false discovery rate < 0.05), was found to be associated mainly with immune subtypes involving interferon response. We isolated an interaction network of interval cancer and interferon genes for which a significant load of somatic and germline variants in class I interferon genes was observed.

CONCLUSION

We identified novel molecular features of interval breast cancer highlighting interferon pathways as a potential target for prevention or treatment.

Seed sequence polymorphism rs2168518 and allele-specific target gene regulation of hsa-miR-4513

MicroRNAs (miRNAs) are small post-transcriptional regulators that offer promising targets for treating complex diseases. To this end, hsa-miR-4513 is an excellent candidate as this gene harbors within its conserved heptametrical seed sequence a frequent polymorphism (rs2168518), which has previously been associated with several complex phenotypes. So far, little is known about the biological mechanism(s) underlying these associations. In an initial step, we now aimed to identify allele-specific target genes of hsa-miR-4513. We performed RNA sequencing in a miRNA overexpression model in human umbilical vein endothelial cells transfected with separated hsa-miR-4513 alleles at rs2168518, namely hsa-miR-4513-G and hsa-miR-4513-A. Genes specifically regulated by the rs2168518 alleles were independently verified by quantitative reverse transcription PCR (qRT-PCR), western blot analysis and allele-specific miRNA binding via a luciferase reporter assay. By a text-based search publicly available databases such as Online Mendelian Inheritance in Man and Mouse Genome Informatics were utilized to link target genes of hsa-miR-4513 to previously described phenotypes. Overall, we identified 23 allele-specific hsa-miR-4513 target genes and replicated 19 of those independently via qRT-PCR. Western blot analysis and luciferase reporter assays conducted for an exemplary subsample further confirmed the allele-specific regulation of these genes by hsa-miR-4513. Remarkably, multiple allele-specific target genes identified are linked via text retrieval to several phenotypes previously reported to be associated with hsa-miR-4513. These genes offer promising candidates for ongoing research on the functional pathobiological impact of hsa-miR-4513 and its seed polymorphism rs2168518. This could give rise to therapeutic applications targeting this miRNA.

Mammographic features are associated with cardiometabolic disease risk and mortality

Aims

In recent years, microcalcifications identified in routine mammograms were found to be associated with cardiometabolic disease in women. Here, we aimed to systematically evaluate the association of microcalcifications and other mammographic features with cardiometabolic disease risk and mortality in a large screening cohort and to understand a potential genetic contribution.

Methods and results

This study included 57 867 women from a prospective mammographic screening cohort in Sweden (KARMA) and 49 583 sisters. Cardiometabolic disease diagnoses and mortality and medication were extracted by linkage to Swedish population registries with virtually no missing data. In the cardiometabolic phenome-wide association study, we found that a higher number of microcalcifications were associated with increased risk for multiple cardiometabolic diseases, particularly in women with pre-existing cardiometabolic diseases. In contrast, dense breasts were associated with a lower incidence of cardiometabolic diseases. Importantly, we observed similar associations in sisters of KARMA women, indicating a potential genetic overlap between mammographic features and cardiometabolic traits. Finally, we observed that the presence of microcalcifications was associated with increased cardiometabolic mortality in women with pre-existing cardiometabolic diseases (hazard ratio and 95% confidence interval: 1.79 [1.24–2.58], P = 0.002) while we did not find such effects in women without cardiometabolic diseases.

Conclusions

We found that mammographic features are associated with cardiometabolic risk and mortality. Our results strengthen the notion that a combination of mammographic features and other breast cancer risk factors could be a novel and affordable tool to assess cardiometabolic health in women attending mammographic screening.

CYP2D6 genotype predicts tamoxifen discontinuation and drug response: a secondary analysis of the KARISMA trial

BACKGROUND

Guidelines regarding whether tamoxifen should be prescribed based on women's cytochrome P450 2D6 (CYP2D6) genotypes are conflicting and have caused confusion. This study aims to investigate if CYP2D6 metabolizer status is associated with tamoxifen-related endocrine symptoms, tamoxifen discontinuation, and mammographic density change.

PATIENTS AND METHODS

We used data from 1440 healthy women who participated the KARISMA dose determination trial. Endocrine symptoms were measured using a modified Functional Assessment of Cancer Therapy - Endocrine Symptoms (FACT-ES) questionnaire. Change in mammographic density was measured and used as a proxy for tamoxifen response. Participants were genotyped and categorized as poor, intermediate, normal, or ultrarapid CYP2D6 metabolizers.

RESULTS

The median endoxifen level per mg oral tamoxifen among poor, intermediate, normal and ultrarapid CYP2D6 metabolizers were 0.18 ng/ml, 0.38 ng/ml, 0.56 ng/ml and 0.67 ng/ml, respectively. Ultrarapid CYP2D6 metabolizers were more likely than other groups to report a clinically relevant change in cold sweats, hot flash, mood swings, being irritable, as well as the overall modified FACT-ES score, after taking tamoxifen. The 6-month tamoxifen discontinuation rates among poor, intermediate, normal, and ultrarapid CYP2D6 metabolizers were 25.7%, 23.6%, 28.6%, and 44.4%, respectively. Among those who continued and finished the 6-month tamoxifen intervention, the mean change in dense area among poor, intermediate, normal, and ultrarapid CYP2D6 metabolizers were -0.8 cm², -4.5 cm², -4.1 cm², and -8.0 cm² respectively.

CONCLUSIONS

Poor CYP2D6 metabolizers are likely to experience an impaired response to tamoxifen, measured through mammographic density reduction. In contrast, ultrarapid CYP2D6 metabolizers are at risk for exaggerated response with pronounced adverse effects that may lead to treatment discontinuation.

Lymph node metastases in breast cancer: Investigating associations with tumor characteristics, molecular subtypes and polygenic risk score using a continuous growth model

We investigate the association between rate of breast cancer lymph node spread and grade, estrogen receptor (ER) status, progesteron receptor status, decision tree derived PAM50 molecular subtype and a polygenic risk score (PRS), using data on 10 950 women included from two different data sources. Lymph node spread was analyzed using a novel continuous tumor progression model that adjusts for tumor volume in a biologically motivated way and that incorporates covariates of interest. Grades 2 and 3 tumors, respectively, were associated with 1.63 and 2.17 times faster rates of lymph node spread than Grade 1 tumors (P < 10^-16 ). ER/PR negative breast cancer was associated with a 1.25/1.19 times faster spread than ER/PR positive breast cancer, respectively (P = .0011 and .0012). Among the molecular subtypes luminal A, luminal B, Her2-enriched and basal-like, Her2-enriched breast cancer was associated with 1.53 times faster spread than luminal A cancer (P = .00072). PRS was not associated with the rate of lymph node spread. Continuous growth models are useful for quantifying associations between lymph node spread and tumor characteristics. These may be useful for building realistic progression models for microsimulation studies used to design individualized screening programs.

Pleiotropic Locus 15q24.1 Reveals a Gender-Specific Association with Neovascular but Not Atrophic Age-Related Macular Degeneration (AMD)

Genome-wide association studies (GWAS) have identified an abundance of genetic loci associated with complex traits and diseases. In contrast, in-depth characterization of an individual genetic signal is rarely available. Here, we focus on the genetic variant rs2168518 in 15q24.1 previously associated with age-related macular degeneration (AMD), but only with suggestive evidence. In a two-step procedure, we initially conducted a series of association analyses to further delineate the association of rs2168518 with AMD but also with other complex phenotypes by using large independent datasets from the International AMD Genomics Consortium (IAMDGC) and the UK Biobank. We then performed a functional annotation with reference to gene expression regulation based on data from the Genotype-Tissue Expression (GTEx) project and RegulomeDB. Association analysis revealed a gender-specific association with male AMD patients and an association predominantly with choroidal neovascularization. Further, the AMD association colocalizes with an association signal of several blood pressure-related phenotypes and with the gene expression regulation of CYP1A1, a member of the cytochrome P450 superfamily of monooxygenases. Functional annotation revealed altered transcription factor (TF) binding sites for gender-specific TFs, including SOX9 and SRY. In conclusion, the pleiotropic 15q24.1 association signal suggests a shared mechanism between blood pressure regulation and choroidal neovascularization with a potential involvement of CYP1A1.

The association of single nucleotide polymorphisms (SNPs) with breast density and breast cancer survival: the Malmö Diet and Cancer Study

Background

Genetic factors are important in determining breast density, and heritable factors account for 60% of the variation. Certain single nucleotide polymorphisms (SNPs) are associated with density and risk of breast cancer but the association with prognosis is not clear.

Purpose

To investigate associations between selected SNPs and breast cancer survival in the Malmö Diet and Cancer Study (MDCS).

Material and Methods

A total of 724 unrelated women with breast cancer and registered radiological and pathological data were identified in MDCS 1991–2007, with genotyping available for 672 women. Associations among 15 SNPs, density, and breast cancer-specific survival were analyzed using logistic/Cox regression, adjusted for factors affecting density and survival. Variants significantly associated with either density or survival were validated in a large independent breast cancer cohort (LIBRO-1).

Results

Minor homozygotes of SNPs rs9383589, CCDC170 and rs6557161, ESR1 were associated with high breast density (adjusted odds ratio [AOR] 8.97, 95% confidence interval [CI] 1.35–59.57; AOR 2.08, 95% CI 1.19–3.65, respectively) and poorer breast cancer survival (adjusted hazard ratio [HR_adj] 6.46, 95% CI 1.95–21.39; HR_adj 2.30, 95% CI 1.33–3.96, respectively) compared to major homozygotes. For SNP rs3757318, ESR1, minor homozygotes (HR_adj 7.46, 95% CI 2.28–24.45) were associated with poorer survival. We confirmed that rs6557161, ESR1 was significantly associated with both density and survival in the LIBRO-1 study.

Conclusion

These findings support a shared genetic basis for density and breast cancer survival. The SNP significantly associated with both density and survival in both cohorts may be of interest in future research investigating polygenic risk scores for breast cancer risk and screening stratification purposes.

Association between breast cancer risk and disease aggressiveness: Characterizing underlying gene expression patterns

The association between breast cancer risk defined by the Tyrer-Cuzick score (TC) and disease prognosis is not well established. Here, we investigated the relationship between 5-year TC and disease aggressiveness and then characterized underlying molecular processes. In a case-only study (n = 2474), we studied the association of TC with molecular subtypes and tumor characteristics. In a subset of patients (n = 672), we correlated gene expression to TC and computed a low-risk TC gene expression (TC-Gx) profile, that is, a profile expected to be negatively associated with risk, which we used to test for association with disease aggressiveness. We performed enrichment analysis to pinpoint molecular processes likely to be altered in low-risk tumors. A higher TC was found to be inversely associated with more aggressive surrogate molecular subtypes and tumor characteristics (P < .05) including Ki-67 proliferation status (P < 5 × 10^-07 ). Our low-risk TC-Gx, based on the weighted sum of 37 expression values of genes strongly correlated with TC, was associated with basal-like (P < 5 × 10^-13 ), HER2-enriched subtype (P < 5 × 10^-07 ) and worse 10-year breast cancer-specific survival (log-rank P < 5 × 10^-04 ). Associations between low-risk TC-Gx and more aggressive molecular subtypes were replicated in an independent cohort from The Cancer Genome Atlas database (n = 975). Gene expression that correlated with low TC was enriched in proliferation and oncogenic signaling pathways (FDR < 0.05). Moreover, higher proliferation was a key factor explaining the association with worse survival. Women who developed breast cancer despite having a low risk were diagnosed with more aggressive tumors and had a worse prognosis, most likely driven by increased proliferation. Our findings imply the need to establish risk factors associated with more aggressive breast cancer subtypes.

A mega-analysis of expression quantitative trait loci in retinal tissue

Significant association signals from genome-wide association studies (GWAS) point to genomic regions of interest. However, for most loci the causative genetic variant remains undefined. Determining expression quantitative trait loci (eQTL) in a disease relevant tissue is an excellent approach to zoom in on disease- or trait-associated association signals and hitherto on relevant disease mechanisms. To this end, we explored regulation of gene expression in healthy retina (n = 311) and generated the largest cis-eQTL data set available to date. Genotype- and RNA-Seq data underwent rigorous quality control protocols before FastQTL was applied to assess the influence of genetic markers on local (cis) gene expression. Our analysis identified 403,151 significant eQTL variants (eVariants) that regulate 3,007 genes (eGenes) (Q-Value < 0.05). A conditional analysis revealed 744 independent secondary eQTL signals for 598 of the 3,007 eGenes. Interestingly, 99,165 (24.71%) of all unique eVariants regulate the expression of more than one eGene. Filtering the dataset for eVariants regulating three or more eGenes revealed 96 potential regulatory clusters. Of these, 31 harbour 130 genes which are partially regulated by the same genetic signal. To correlate eQTL and association signals, GWAS data from twelve complex eye diseases or traits were included and resulted in identification of 80 eGenes with potential association. Remarkably, expression of 10 genes is regulated by eVariants associated with multiple eye diseases or traits. In conclusion, we generated a unique catalogue of gene expression regulation in healthy retinal tissue and applied this resource to identify potentially pleiotropic effects in highly prevalent human eye diseases. Our study provides an excellent basis to further explore mechanisms of various retinal disease etiologies.

The retina is a multilayered and highly specified neural tissue crucial for high-resolution visual perception and spatial orientation. Environmental and genetic insults to the retina result in many blinding diseases, such as age-related macular degeneration or glaucoma. Commonly, many of these diseases are age-related suggesting that minor changes are accumulating over a life-time, with little or no contribution of strong individual effects. Specifically, this is true for genetic factors known to underlie the etiology of complex diseases including the prevalent eye diseases. In our study, we searched for effects on gene expression due to genetic variation using 311 healthy post-mortem retinal tissue samples. We show that 3,007 of the 16,766 genes investigated are regulated in the retina by genetic variations. Of these, 80 genes are potentially associated to one or more of twelve complex eye diseases or retinal traits tested. Interestingly, 10 genes appear to be involved in the development of several eye traits suggesting that cellular mechanisms may act at a common point in the disease process. Consequently, our study provides the basis to further explore retinal disease pathways and is likely to highlight target molecules for future therapeutic applications.

Genome-wide association meta-analysis for early age-related macular degeneration highlights novel loci and insights for advanced disease

BACKGROUND

Advanced age-related macular degeneration (AMD) is a leading cause of blindness. While around half of the genetic contribution to advanced AMD has been uncovered, little is known about the genetic architecture of early AMD.

METHODS

To identify genetic factors for early AMD, we conducted a genome-wide association study (GWAS) meta-analysis (14,034 cases, 91,214 controls, 11 sources of data including the International AMD Genomics Consortium, IAMDGC, and UK Biobank, UKBB). We ascertained early AMD via color fundus photographs by manual grading for 10 sources and via an automated machine learning approach for > 170,000 photographs from UKBB. We searched for early AMD loci via GWAS and via a candidate approach based on 14 previously suggested early AMD variants.

RESULTS

Altogether, we identified 10 independent loci with statistical significance for early AMD: (i) 8 from our GWAS with genome-wide significance (P < 5 × 10- 8), (ii) one previously suggested locus with experiment-wise significance (P < 0.05/14) in our non-overlapping data and with genome-wide significance when combining the reported and our non-overlapping data (together 17,539 cases, 105,395 controls), and (iii) one further previously suggested locus with experiment-wise significance in our non-overlapping data. Of these 10 identified loci, 8 were novel and 2 known for early AMD. Most of the 10 loci overlapped with known advanced AMD loci (near ARMS2/HTRA1, CFH, C2, C3, CETP, TNFRSF10A, VEGFA, APOE), except two that have not yet been identified with statistical significance for any AMD. Among the 17 genes within these two loci, in-silico functional annotation suggested CD46 and TYR as the most likely responsible genes. Presence or absence of an early AMD effect distinguished the known pathways of advanced AMD genetics (complement/lipid pathways versus extracellular matrix metabolism).

CONCLUSIONS

Our GWAS on early AMD identified novel loci, highlighted shared and distinct genetics between early and advanced AMD and provides insights into AMD etiology. Our data provide a resource comparable in size to the existing IAMDGC data on advanced AMD genetics enabling a joint view. The biological relevance of this joint view is underscored by the ability of early AMD effects to differentiate the major pathways for advanced AMD.

Hyperthyroidism is associated with breast cancer risk and mammographic and genetic risk predictors

BACKGROUND

Despite the biological link between thyroid hormones and breast cancer cell proliferation shown in experimental studies, little is known about the association between hyperthyroidism and breast cancer, as well as its association with the most common mammographic and genetic risk predictors for breast cancer.

METHODS

This study estimates the incidence rate ratios (IRRs) of breast cancer among women diagnosed with hyperthyroidism, compared to those who are not, using two cohorts: a Swedish national cohort of the general female population (n = 3,793,492, 2002-2011) and the Karolinska Mammography Project for Risk Prediction of Breast Cancer (KARMA, n = 69,598, 2002-2017). We used logistic regression to estimate the odds ratios (ORs) of hyperthyroidism according to the mammographic and genetic risk predictors for breast cancer.

RESULTS

An increased risk of breast cancer was observed in patients in the national cohort with hyperthyroidism (IRR = 1.23, 95% CI = 1.12-1.36), particularly for toxic nodular goiter (IRR = 1.38, 95% CI = 1.16-1.63). Hyperthyroidism was associated with higher body mass index, early age at first birth, and lower breastfeeding duration. Higher mammographic density was observed in women with toxic nodular goiter, compared to women without hyperthyroidism. Additionally, among genotyped women without breast cancer in the KARMA cohort (N = 11,991), hyperthyroidism was associated with a high polygenic risk score (PRS) for breast cancer overall (OR = 1.98, 95% CI = 1.09-3.60) and for estrogen receptor-positive specific PRS (OR = 1.90, 95% CI = 1.04-3.43).

CONCLUSION

Hyperthyroidism is associated with an increased risk of breast cancer, particularly for patients with toxic nodular goiter. The association could be explained by higher mammographic density among these women, as well as pleiotropic genetic variants determining shared hormonal/endocrine factors leading to the pathology of both diseases.

A Circulating MicroRNA Profile in a Laser-Induced Mouse Model of Choroidal Neovascularization

Choroidal neovascularization (CNV) is a pathological process in which aberrant blood vessels invade the subretinal space of the mammalian eye. It is a characteristic feature of the prevalent neovascular age-related macular degeneration (nAMD). Circulating microRNAs (cmiRNAs) are regarded as potentially valuable biomarkers for various age-related diseases, including nAMD. Here, we investigated cmiRNA expression in an established laser-induced CNV mouse model. Upon CNV induction in C57Bl/6 mice, blood-derived cmiRNAs were initially determined globally by RNA next generation sequencing, and the most strongly dysregulated cmiRNAs were independently replicated by quantitative reverse transcription PCR (RT-qPCR) in blood, retinal, and retinal pigment epithelium (RPE)/choroidal tissue. Our findings suggest that two miRNAs, mmu-mir-486a-5p and mmur-mir-92a-3p, are consistently dysregulated during CNV formation. Furthermore, in functional in vitro assays, a significant impact of mmu-mir-486a-5p and mmu-mir-92a-3p on murine microglial cell viability was observed, while mmu-mir-92a-3p also showed an impact on microglial mobility. Taken together, we report a robust dysregulation of two miRNAs in blood and RPE/choroid after laser-induced initiation of CNV lesions in mice, highlighting their potential role in pathology and eventual therapy of CNV-associated complications.

CYP2D6 Genotype Predicts Tamoxifen Discontinuation and Prognosis in Patients With Breast Cancer

PURPOSE

To examine the association between CYP2D6 genotype, discontinuation of tamoxifen therapy, and prognosis for breast cancer.

PATIENTS AND METHODS

We conducted a prospective-retrospective study linking data from a clinical breast cancer register, the Swedish Prescribed Drug Register, and self-reported questionnaires. We genotyped CYP2D6 in 1,309 patients with breast cancer who were treated with tamoxifen and were diagnosed from 2005 to 2012; they were categorized as poor, intermediate, normal, or ultrarapid CYP2D6 metabolizers. We investigated whether metabolizer status was associated with tamoxifen discontinuation and prognosis for breast cancer using Cox regression analysis.

RESULTS

The 6-month discontinuation rates of tamoxifen among poor, intermediate, normal, and ultrarapid CYP2D6 metabolizers were 7.1%, 7.6%, 6.7%, and 18.8%, respectively. A U-shaped association was found between CYP2D6 metabolizer status and breast cancer–specific mortality, with adjusted hazard ratios of 2.59 (95% CI, 1.01 to 6.67) for poor, 1.48 (95% CI, 0.72 to 3.05) for intermediate, 1 (reference) for normal, and 4.52 (95% CI, 1.42 to 14.37) for ultrarapid CYP2D6 metabolizers.

CONCLUSION

Both poor and ultrarapid CYP2D6 metabolizers of tamoxifen have a worse prognosis for breast cancer compared with normal metabolizers after receiving a standard dose of tamoxifen. This U-shaped association might call for individualized tamoxifen dosage.

Assessment of Novel Genome-Wide Significant Gene Loci and Lesion Growth in Geographic Atrophy Secondary to Age-Related Macular Degeneration

Importance

Age-related macular degeneration (AMD) is a common threat to vision loss in individuals older than 50 years. While neovascular complications in AMD are treatable, there is currently no therapy for geographic atrophy secondary to AMD. Geographic atrophy lesion progression over time shows considerable interindividual variability, but little is known about prognostic factors.

Objective

To elucidate the contribution of common genetic variants to geographic atrophy lesion growth.

Design, Setting, and Participants

This pooled analysis combined 4 independent studies: the Fundus Autofluorescence Imaging in Age-Related Macular Degeneration (FAM) study, the Directional Spread in Geographic Atrophy (DSGA) study, the Age-Related Eye Disease Study (AREDS), and the Geographic Atrophy Treatment Evaluation (GATE) study. Each provided data for geographic atrophy lesion growth in specific designs. Patients with geographic atrophy secondary to AMD were recruited to these studies. Genotypes were retrieved through the database of Genotypes and Phenotypes (for AREDS) or generated at the Cologne Center for Genomics (for FAM, DSGA, and GATE).

Main Outcomes

The correlation between square root-transformed geographic atrophy growth rate and 7 596 219 genetic variants passing quality control was estimated using linear regression. The calculations were adjusted for known factors influencing geographic atrophy growth, such as the presence of bilateral geographic atrophy as well as the number of lesion spots and follow-up times.

Main Outcomes and Measures

Slopes per allele, 95% CIs, and P values of genetic variants correlated with geographic atrophy lesion growth.

Results

A total of 935 patients (mean [SD] age, 74.7 [7.8] years; 547 female participants [59.0%]) were included. Two gene loci with conservative genome-wide significance were identified. Each minor allele of the genome-wide associated variants increased the geographic atrophy growth rate by a mean of about 15% or 0.05 mm per year. Gene prioritization within each locus suggests the protein arginine methyltransferase 6 gene (PRMT6; chromosome 1; slope, 0.046 [95% CI, 0.026-0.066]; P = 4.09 × 10-8) and the lanosterol synthase gene (LSS; chromosome 21; slope, 0.105 [95% CI, 0.068-0.143]; P = 4.07 × 10-7) as the most likely progression-associated genes.

Conclusions and Relevance

These data provide further insight into the genetic architecture of geographic atrophy lesion growth. Geographic atrophy is a clinical outcome with a high medical need for effective therapy. The genes PRMT6 and LSS are promising candidates for future studies aimed at understanding functional aspects of geographic atrophy progression and also for designing novel and targeted treatment options.

Cost-effective molecular inversion probe-based ABCA4 sequencing reveals deep-intronic variants in Stargardt disease

PURPOSE

Stargardt disease (STGD1) is caused by biallelic mutations in ABCA4, but many patients are genetically unsolved due to insensitive mutation-scanning methods. We aimed to develop a cost-effective sequencing method for ABCA4 exons and regions carrying known causal deep-intronic variants.

METHODS

Fifty exons and 12 regions containing 14 deep-intronic variants of ABCA4 were sequenced using double-tiled single molecule Molecular Inversion Probe (smMIP)-based next-generation sequencing. DNAs of 16 STGD1 cases carrying 29 ABCA4 alleles and of four healthy persons were sequenced using 483 smMIPs. Thereafter, DNAs of 411 STGD1 cases with one or no ABCA4 variant were sequenced. The effect of novel noncoding variants on splicing was analyzed using in vitro splice assays.

RESULTS

Thirty-four ABCA4 variants previously identified in 16 STGD1 cases were reliably identified. In 155/411 probands (38%), two causal variants were identified. We identified 11 deep-intronic variants present in 62 alleles. Two known and two new noncanonical splice site variants showed splice defects, and one novel deep-intronic variant (c.4539+2065C>G) resulted in a 170-nt mRNA pseudoexon insertion (p.[Arg1514Lysfs*35,=]).

CONCLUSIONS

smMIPs-based sequence analysis of coding and selected noncoding regions of ABCA4 enabled cost-effective mutation detection in STGD1 cases in previously unsolved cases.

The agonistic TSPO ligand XBD173 attenuates the glial response thereby protecting inner retinal neurons in a murine model of retinal ischemia

BACKGROUND

Ligand-driven modulation of the mitochondrial translocator protein 18 kDa (TSPO) was recently described to dampen the neuroinflammatory response of microglia in a retinal light damage model resulting in protective effects on photoreceptors. We characterized the effects of the TSPO ligand XBD173 in the postischemic retina focusing on changes in the response pattern of the major glial cell types of the retina-microglia and Müller cells.

METHODS

Retinal ischemia was induced by increasing the intraocular pressure for 60 min followed by reperfusion of the tissue in mice. On retinal cell types enriched via immunomagnetic separation expression analysis of TSPO, its ligand diazepam-binding inhibitor (DBI) and markers of glial activation were performed at transcript and protein level using RNA sequencing, qRT-PCR, lipid chromatography-mass spectrometry, and immunofluorescent labeling. Data on cell morphology and numbers were assessed in retinal slice and flatmount preparations. The retinal functional integrity was determined by electroretinogram recordings.

RESULTS

We demonstrate that TSPO is expressed by Müller cells, microglia, vascular cells, retinal pigment epithelium (RPE) of the healthy and postischemic retina, but only at low levels in retinal neurons. While an alleviated neurodegeneration upon XBD173 treatment was found in postischemic retinae as compared to vehicle controls, this neuroprotective effect of XBD173 is mediated putatively by its action on retinal glia. After transient ischemia, TSPO as a marker of activation was upregulated to similar levels in microglia as compared to their counterparts in healthy retinae irrespective of the treatment regimen. However, less microglia were found in XBD173-treated postischemic retinae at 3 days post-surgery (dps) which displayed a more ramified morphology than in retinae of vehicle-treated mice indicating a dampened microglia activation. Müller cells, the major retinal macroglia, show upregulation of the typical gliosis marker GFAP. Importantly, glutamine synthetase was more stably expressed in Müller glia of XBD173-treated postischemic retinae and homeostatic functions such as cellular volume regulation typically diminished in gliotic Müller cells remained functional.

CONCLUSIONS

In sum, our data imply that beneficial effects of XBD173 treatment on the postischemic survival of inner retinal neurons were primarily mediated by stabilizing neurosupportive functions of glial cells.

Conduct and Quality Control of Differential Gene Expression Analysis Using High-Throughput Transcriptome Sequencing (RNASeq)

High-throughput transcriptome sequencing (RNASeq) represents one of the most comprehensive and scalable methods to analyze global gene expression. It allows for absolute quantification of gene expression and also enables the discovery of novel transcripts and alternatively spliced isoforms. This chapter provides hand-on tools and a step-by-step procedure to analyze RNASeq data from punctures of two different retinal tissues (retina and RPE-choroid-sclera) at two different locations (periphery and macular region) from eight individuals. The procedure described in this chapter will use various programs from the free, open-source Tuxedo Suite software package to analyze sequencing data and to ascertain genes that are differentially expressed between retina and RPE-choroid-sclera.

Y chromosome mosaicism is associated with age-related macular degeneration

Age-related macular degeneration (AMD) is the leading cause of blindness in industrialised countries, and thereby a major individual but also a socio-economic burden. Y chromosome loss in nucleated blood cells has been implicated in age-related diseases such as Alzheimer disease and was shown to be caused by increasing age, smoking and genetic factors. Mosaic loss of Y chromosome (mLOY) in peripheral blood was estimated from normalised dosages of genotyping chip data covering the male-specific region of the Y chromosome. After quality control, we assessed the association of mLOY on AMD risk in 5772 male cases and 6732 male controls. In controls the prevalence of mLOY increased significantly with age, which is consistent with previous reports. Importantly, mLOY was associated with late-stage AMD with genome-wide significance (OR: 1.332 [95% CI: 1.206; 1.472], P = 1.60e-08), independent of age, the AMD genetic risk score and the first two principle components of ancestry. Additionally conditioning on smoking behaviour had no influence on the observed association strength. mLOY was strongest associated in individuals aged between 65 and 75 years. Taken together, mLOY is significantly associated with risk for AMD, independent of known and potential confounding factors.

A Deep Learning Algorithm for Prediction of Age-Related Eye Disease Study Severity Scale for Age-Related Macular Degeneration from Color Fundus Photography

PURPOSE

Age-related macular degeneration (AMD) is a common threat to vision. While classification of disease stages is critical to understanding disease risk and progression, several systems based on color fundus photographs are known. Most of these require in-depth and time-consuming analysis of fundus images. Herein, we present an automated computer-based classification algorithm.

DESIGN

Algorithm development for AMD classification based on a large collection of color fundus images. Validation is performed on a cross-sectional, population-based study.

PARTICIPANTS

We included 120 656 manually graded color fundus images from 3654 Age-Related Eye Disease Study (AREDS) participants. AREDS participants were >55 years of age, and non-AMD sight-threatening diseases were excluded at recruitment. In addition, performance of our algorithm was evaluated in 5555 fundus images from the population-based Kooperative Gesundheitsforschung in der Region Augsburg (KORA; Cooperative Health Research in the Region of Augsburg) study.

METHODS

We defined 13 classes (9 AREDS steps, 3 late AMD stages, and 1 for ungradable images) and trained several convolution deep learning architectures. An ensemble of network architectures improved prediction accuracy. An independent dataset was used to evaluate the performance of our algorithm in a population-based study.

MAIN OUTCOME MEASURES

κ Statistics and accuracy to evaluate the concordance between predicted and expert human grader classification.

RESULTS

A network ensemble of 6 different neural net architectures predicted the 13 classes in the AREDS test set with a quadratic weighted κ of 92% (95% confidence interval, 89%-92%) and an overall accuracy of 63.3%. In the independent KORA dataset, images wrongly classified as AMD were mainly the result of a macular reflex observed in young individuals. By restricting the KORA analysis to individuals >55 years of age and prior exclusion of other retinopathies, the weighted and unweighted κ increased to 50% and 63%, respectively. Importantly, the algorithm detected 84.2% of all fundus images with definite signs of early or late AMD. Overall, 94.3% of healthy fundus images were classified correctly.

CONCLUSIONS

Our deep learning algoritm revealed a weighted κ outperforming human graders in the AREDS study and is suitable to classify AMD fundus images in other datasets using individuals >55 years of age.

Investigating the modulation of genetic effects on late AMD by age and sex: Lessons learned and two additional loci

Late-stage age-related macular degeneration (AMD) is the leading cause of visual impairment in the elderly with a complex etiology. The most important non-modifiable risk factors for onset and progression of late AMD are age and genetic risk factors, however, little is known about the interplay between genetics and age or sex. Here, we conducted a large-scale age- and sex-stratified genome-wide association study (GWAS) using 1000 Genomes imputed genome-wide and ExomeChip data (>12 million variants). The data were established by the International Age-related Macular Degeneration Genomics Consortium (IAMDGC) from 16,144 late AMD cases and 17,832 controls. Our systematic search for interaction effects yielded significantly stronger effects among younger individuals at two known AMD loci (near CFH and ARMS2/HTRA1). Accounting for age and gene-age interaction using a joint test identified two additional AMD loci compared to the previous main effect scan. One of these two is a novel AMD GWAS locus, near the retinal clusterin-like protein (CLUL1) gene, and the other, near the retinaldehyde binding protein 1 (RLBP1), was recently identified in a joint analysis of nuclear and mitochondrial variants. Despite considerable power in our data, neither sex-dependent effects nor effects with opposite directions between younger and older individuals were observed. This is the first genome-wide interaction study to incorporate age, sex and their interaction with genetic effects for late AMD. Results diminish the potential for a role of sex in the etiology of late AMD yet highlight the importance and existence of age-dependent genetic effects.

Mutation Spectrum of the ABCA4 Gene in 335 Stargardt Disease Patients From a Multicenter German Cohort-Impact of Selected Deep Intronic Variants and Common SNPs

Purpose

Stargardt disease (STGD1) is an autosomal recessive retinopathy, caused by mutations in the retina-specific ATP-binding cassette transporter (ABCA4) gene. To establish the mutational spectrum and to assess effects of selected deep intronic and common genetic variants on disease, we performed a comprehensive sequence analysis in a large cohort of German STGD1 patients.

Methods

DNA samples of 335 STGD1 patients were analyzed for ABCA4 mutations in its 50 coding exons and adjacent intronic sequences by resequencing array technology or next generation sequencing (NGS). Parts of intron 30 and 36 were screened by Sanger chain-terminating dideoxynucleotide sequencing. An in vitro splicing assay was used to test selected variants for their splicing behavior. By logistic regression analysis we assessed the association of common ABCA4 alleles while a multivariate logistic regression model calculated a genetic risk score (GRS).

Results

Our analysis identified 148 pathogenic or likely pathogenic mutations, of which 48 constitute so far unpublished ABCA4-associated disease alleles. Four rare deep intronic variants were found once in 472 alleles analyzed. In addition, we identified six risk-modulating common variants. Genetic risk score estimates suggest that defined common ABCA4 variants influence disease risk in carriers of a single pathogenic ABCA4 allele.

Conclusions

Our study adds to the mutational spectrum of the ABCA4 gene. Moreover, in our cohort, deep intronic variants in intron 30 and 36 likely play no or only a minor role in disease pathology. Of note, our findings demonstrate a possible modifying effect of common sequence variants on ABCA4-associated disease.

Distinct Genetic Risk Profile of the Rapidly Progressing Diffuse-Trickling Subtype of Geographic Atrophy in Age-Related Macular Degeneration (AMD)

PURPOSE

To genetically characterize a subphenotype of geographic atrophy (GA) in AMD associated with rapid progression and a diffuse-trickling appearance on fundus autofluorescence imaging.

METHODS

Patients from the Fundus Autofluorescence in Age-Related Macular Degeneration Study were phenotyped for diffuse-trickling GA (dt-GA; n = 44). DNA was analyzed for 10 known AMD-associated genetic variants. A genetic risk score (GRS) was calculated and compared with patients with nondiffuse-trickling GA (ndt-GA; n = 311) and individuals from the 1000 genomes project (1000G; n = 267). Given the phenotypic overlap between diffuse-trickling and late-onset retinal degeneration (LORD), all C1QTNF5 exons and their exon/intron boundaries were sequenced.

RESULTS

A statistically significant difference in allele frequencies between dt-GA and ndt-GA were found for CFH:rs1061170 and CFH:rs800292 (Pcorrected = 0.03). The ARMS2 variant rs10490924 was significantly more frequent in dt-GA than in 1000G individuals (Pcorrected < 0.01). The GRS of dt-GA patients was in-between the score of the 1000G individuals and that of patients with ndt-GA, significantly differing from both (Pcorrected <0.01). Sequencing of C1QTNF5 revealed 28 unique variants although none showed a statistically significant association with dt-GA when compared with 1000G individuals.

CONCLUSIONS

The dt-GA phenotype shows a remarkably different genetic risk profile from other GA phenotypes secondary to AMD. Disease-associated C1QTNF5 mutations were not identified. Together, these results suggest that the dt-GA phenotype is associated with a genetic background substantially different from other GA phenotypes and underlines the necessity to refine the clinical phenotyping, specifically when aiming for individualized therapies in AMD.

Genetic pleiotropy between age-related macular degeneration and 16 complex diseases and traits

BACKGROUND

Age-related macular degeneration (AMD) is a common condition of vision loss with disease development strongly influenced by environmental and genetic factors. Recently, 34 loci were associated with AMD at genome-wide significance. So far, little is known about a genetic overlap between AMD and other complex diseases or disease-relevant traits.

METHODS

For each of 60 complex diseases/traits with publicly available genome-wide significant association data, the lead genetic variant per independent locus was extracted and a genetic score was calculated for each disease/trait as the weighted sum of risk alleles. The association with AMD was estimated based on 16,144 AMD cases and 17,832 controls using logistic regression.

RESULTS

Of the respective disease/trait variance, the 60 genetic scores explained on average 4.8% (0.27-20.69%) and 16 of them were found to be significantly associated with AMD (Q-values < 0.01, p values from < 1.0 × 10^-16 to 1.9 × 10^-3). Notably, an increased risk for AMD was associated with reduced risk for cardiovascular diseases, increased risk for autoimmune diseases, higher HDL and lower LDL levels in serum, lower bone-mineral density as well as an increased risk for skin cancer. By restricting the analysis to 1824 variants initially used to compute the 60 genetic scores, we identified 28 novel AMD risk variants (Q-values < 0.01, p values from 1.1 × 10^-7 to 3.0 × 10^-4), known to be involved in cardiovascular disorders, lipid metabolism, autoimmune diseases, anthropomorphic traits, ocular disorders, and neurological diseases. The latter variants represent 20 novel AMD-associated, pleiotropic loci. Genes in the novel loci reinforce previous findings strongly implicating the complement system in AMD pathogenesis.

CONCLUSIONS

We demonstrate a substantial overlap of the genetics of several complex diseases/traits with AMD and provide statistically significant evidence for an additional 20 loci associated with AMD. This highlights the possibility that so far unrelated pathologies may have disease pathways in common.

An Eye on Age-Related Macular Degeneration: The Role of MicroRNAs in Disease Pathology

Age-related macular degeneration (AMD) is the primary cause of blindness in developed countries, and is the third leading cause worldwide. Emerging evidence suggests that beside environmental and genetic factors, epigenetic mechanisms, such as microRNA (miRNA) regulation of gene expression, are relevant to AMD providing an exciting new avenue for research and therapy. MiRNAs are short, non-coding RNAs thought to be imperative for coping with cellular stress. Numerous studies have analyzed miRNA dysregulation in AMD patients, although with varying outcomes. Four studies which profiled dysregulated circulating miRNAs in AMD yielded unique sets, and there is only minimal overlap in ocular miRNA profiling of AMD. Mouse models of AMD, including oxygen-induced retinopathy and laser-induced choroidal neovascularization, showed similarities to some extent with miRNA patterns in AMD. For example, miR-146a is an extensively researched miRNA thought to modulate inflammation, and was found to be upregulated in AMD mice and cellular systems, but also in human AMD retinae and vitreous humor. Similarly, mir-17, miR-125b and miR-155 were dysregulated in multiple AMD mouse models as well as in human AMD plasma or retinae. These miRNAs are thought to regulate angiogenesis, apoptosis, phagocytosis, and inflammation. A promising avenue of research is the modulation of such miRNAs, as the phenotype of AMD mice could be ameliorated with antagomirs or miRNA-mimic treatment. However, before meaningful strides can be made to develop miRNAs as a diagnostic or therapeutic tool, reproducible miRNA profiles need to be established for the various clinical outcomes of AMD.