Genetic Variations in MicroRNA-Binding Sites Affect MicroRNA-Mediated Regulation of Several Genes Associated With Cardio-metabolic Phenotypes

Fibroblast activation protein: Pivoting cancer/chemotherapeutic insight towards heart failure

An important mechanism for cancer progression is degradation of the extracellular matrix (ECM) which is accompanied by the emergence and proliferation of an activated fibroblast, termed the cancer associated fibroblast (CAF). More specifically, an enzyme pathway identified to be amplified with local cancer progression and proliferation of the CAF, is fibroblast activation protein (FAP). The development and progression of heart failure (HF) irrespective of the etiology is associated with left ventricular (LV) remodeling and changes in ECM structure and function. As with cancer, HF progression is associated with a change in LV myocardial fibroblast growth and function, and expresses a protein signature not dissimilar to the CAF. The overall goal of this review is to put forward the postulate that scientific discoveries regarding FAP in cancer as well as the development of specific chemotherapeutics could be pivoted to target the emergence of FAP in the activated fibroblast subtype and thus hold translationally relevant diagnostic and therapeutic targets in HF.

Phenome-wide association study on miRNA-related sequence variants: the UK Biobank

BACKGROUND

Genetic variants in the coding region could directly affect the structure and expression levels of genes and proteins. However, the importance of variants in the non-coding region, such as microRNAs (miRNAs), remain to be elucidated. Genetic variants in miRNA-related sequences could affect their biogenesis or functionality and ultimately affect disease risk. Yet, their implications and pleiotropic effects on many clinical conditions remain unknown.

METHODS

Here, we utilised genotyping and hospital records data in the UK Biobank (N = 423,419) to investigate associations between 346 genetic variants in miRNA-related sequences and a wide range of clinical diagnoses through phenome-wide association studies. Further, we tested whether changes in blood miRNA expression levels could affect disease risk through colocalisation and Mendelian randomisation analysis.

RESULTS

We identified 122 associations for six variants in the seed region of miRNAs, nine variants in the mature region of miRNAs, and 27 variants in the precursor miRNAs. These included associations with hypertension, dyslipidaemia, immune-related disorders, and others. Nineteen miRNAs were associated with multiple diagnoses, with six of them associated with multiple disease categories. The strongest association was reported between rs4285314 in the precursor of miR-3135b and celiac disease risk (odds ratio (OR) per effect allele increase = 0.37, P = 1.8 × 10-162). Colocalisation and Mendelian randomisation analysis highlighted potential causal role of miR-6891-3p in dyslipidaemia.

CONCLUSIONS

Our study demonstrates the pleiotropic effect of miRNAs and offers insights to their possible clinical importance.

Association of Circulating Long Noncoding 7S RNA with Deep Vein Thrombosis

Mitochondrial dysfunction is a recognized factor in the pathogenesis of deep vein thrombosis (DVT). The role of 7S RNA, a long noncoding RNA that plays an important role in mitochondrial function, in DVT remains unclear. In this study, we aimed to investigate the potential use of 7S RNA as a biomarker in DVT. Plasma samples were obtained from 237 patients (aged 16-95 years) with suspected DVT recruited in a prospective multicenter management study (SCORE) where 53 patients were objectively confirmed with a diagnosis of DVT and the rest were diagnosed as non-DVT. 7S RNA was measured using quantitative real-time polymerase chain reaction in plasma samples. The plasma expression of 7S RNA was significantly lower in DVT compared with non-DVT (0.50 vs. 0.95, p = 0.043). With the linear regression analysis, we showed that the association between the plasma expression of 7S RNA and DVT (β = -0.72, p = 0.007) was independent of potential confounders. Receiver-operating characteristic curve analysis showed the area under the curve values of 0.60 for 7S RNA. The findings of the present study showed a notable association between 7S RNA and DVT. However, further investigations are needed to fully elucidate the exact role of 7S RNA in the pathophysiology of DVT and its diagnostic value.

PTPA variants and impaired PP2A activity in early-onset parkinsonism with intellectual disability

The protein phosphatase 2A complex (PP2A), the major Ser/Thr phosphatase in the brain, is involved in a number of signalling pathways and functions, including the regulation of crucial proteins for neurodegeneration, such as alpha-synuclein, tau and LRRK2. Here, we report the identification of variants in the PTPA/PPP2R4 gene, encoding a major PP2A activator, in two families with early-onset parkinsonism and intellectual disability. We carried out clinical studies and genetic analyses, including genome-wide linkage analysis, whole-exome sequencing, and Sanger sequencing of candidate variants. We next performed functional studies on the disease-associated variants in cultured cells and knock-down of ptpa in Drosophila melanogaster. We first identified a homozygous PTPA variant, c.893T>G (p.Met298Arg), in patients from a South African family with early-onset parkinsonism and intellectual disability. Screening of a large series of additional families yielded a second homozygous variant, c.512C>A (p.Ala171Asp), in a Libyan family with a similar phenotype. Both variants co-segregate with disease in the respective families. The affected subjects display juvenile-onset parkinsonism and intellectual disability. The motor symptoms were responsive to treatment with levodopa and deep brain stimulation of the subthalamic nucleus. In overexpression studies, both the PTPA p.Ala171Asp and p.Met298Arg variants were associated with decreased PTPA RNA stability and decreased PTPA protein levels; the p.Ala171Asp variant additionally displayed decreased PTPA protein stability. Crucially, expression of both variants was associated with decreased PP2A complex levels and impaired PP2A phosphatase activation. PTPA orthologue knock-down in Drosophila neurons induced a significant impairment of locomotion in the climbing test. This defect was age-dependent and fully reversed by L-DOPA treatment. We conclude that bi-allelic missense PTPA variants associated with impaired activation of the PP2A phosphatase cause autosomal recessive early-onset parkinsonism with intellectual disability. Our findings might also provide new insights for understanding the role of the PP2A complex in the pathogenesis of more common forms of neurodegeneration.

Differential Expression of Non-Coding RNAs in Stem Cell Development and Therapeutics of Bone Disorders

Stem cells' self-renewal and multi-lineage differentiation are regulated by a complex network consisting of signaling factors, chromatin regulators, transcription factors, and non-coding RNAs (ncRNAs). Diverse role of ncRNAs in stem cell development and maintenance of bone homeostasis have been discovered recently. The ncRNAs, such as long non-coding RNAs, micro RNAs, circular RNAs, small interfering RNA, Piwi-interacting RNAs, etc., are not translated into proteins but act as essential epigenetic regulators in stem cells' self-renewal and differentiation. Different signaling pathways are monitored efficiently by the differential expression of ncRNAs, which function as regulatory elements in determining the fate of stem cells. In addition, several species of ncRNAs could serve as potential molecular biomarkers in early diagnosis of bone diseases, including osteoporosis, osteoarthritis, and bone cancers, ultimately leading to the development of new therapeutic strategies. This review aims to explore the specific roles of ncRNAs and their effective molecular mechanisms in the growth and development of stem cells, and in the regulation of osteoblast and osteoclast activities. Furthermore, we focus on and explore the association of altered ncRNA expression with stem cells and bone turnover.

In silico prioritisation of microRNA-associated common variants in multiple sclerosis

BACKGROUND

Genome-wide association studies (GWAS) have highlighted over 200 autosomal variants associated with multiple sclerosis (MS). However, variants in non-coding regions such as those encoding microRNAs have not been explored thoroughly, despite strong evidence of microRNA dysregulation in MS patients and model organisms. This study explores the effect of microRNA-associated variants in MS, through the largest publicly available GWAS, which involved 47,429 MS cases and 68,374 controls.

METHODS

We identified SNPs within the coordinates of microRNAs, ± 5-kb microRNA flanking regions and predicted 3'UTR target-binding sites using miRBase v22, TargetScan 7.0 RNA22 v2.0 and dbSNP v151. We established the subset of microRNA-associated SNPs which were tested in the summary statistics of the largest MS GWAS by intersecting these datasets. Next, we prioritised those microRNA-associated SNPs which are among known MS susceptibility SNPs, are in strong linkage disequilibrium with the former or meet a microRNA-specific Bonferroni-corrected threshold. Finally, we predicted the effects of those prioritised SNPs on their microRNAs and 3'UTR target-binding sites using TargetScan v7.0, miRVaS and ADmiRE.

RESULTS

We have identified 30 candidate microRNA-associated variants which meet at least one of our prioritisation criteria. Among these, we highlighted one microRNA variant rs1414273 (MIR548AC) and four 3'UTR microRNA-binding site variants within SLC2A4RG (rs6742), CD27 (rs1059501), MMEL1 (rs881640) and BCL2L13 (rs2587100). We determined changes to the predicted microRNA stability and binding site recognition of these microRNA and target sites.

CONCLUSIONS

We have systematically examined the functional, structural and regulatory effects of candidate MS variants among microRNAs and 3'UTR targets. This analysis allowed us to identify candidate microRNA-associated MS SNPs and highlights the value of prioritising non-coding RNA variation in GWAS. These candidate SNPs could influence microRNA regulation in MS patients. Our study is the first thorough investigation of both microRNA and 3'UTR target-binding site variation in multiple sclerosis using GWAS summary statistics.

Genetic Polymorphisms in the 3'-Untranslated Regions of SMAD5, FN3KRP, and RUNX-1 Are Associated with Recurrent Pregnancy Loss

Recurrent pregnancy loss (RPL) is typically defined as two or more consecutive pregnancy losses prior to 20 weeks of gestation. Although the causes of idiopathic RPL are not completely understood, vascular development and glucose concentration were reported to correlate with the pregnancy loss. The TGF-β signaling pathway which plays a significant role in pregnancy is activated by the interaction between high glucose and SMAD signaling and affects the vascular cells. SMAD5 and RUNX-1 are involved in the TGF-β signaling pathway and contribute to advanced glycation end products (AGEs) production and vascular development. FN3KRP, a newly described gene, is also associated with vascular diseases and suggested to relate to AGEs. Therefore, in the present study, we investigated associations between RPL risk and genetic polymorphisms of SMAD5, FN3KRP, and RUNX-1 in 388 women with RPL and 280 healthy control women of Korean ethnicity. Participants were genotyped using real-time polymerase chain reaction and restriction fragment length polymorphism assay to determine the frequency of SMAD5 rs10515478 C>G, FN3KRP rs1046875 G>A, and RUNX-1 rs15285 G>A polymorphisms. We found that women with RPL had lower likelihoods of the FN3KRP rs1046875 AA genotype (adjusted odds ratio (AOR), 0.553; p = 0.010) and recessive model (AOR, 0.631; p = 0.017). Furthermore, combination analysis showed that SMAD5 rs10515478 C>G and FN3KRP rs1046875 G>A mutant alleles were together associated with reduced RPL risk. These findings suggest that the FN3KRP rs1046875 G>A polymorphism has a significant role on the prevalence of RPL in Korean women. Considering that it is the first study indicating a significant association between FN3KRP and pregnancy disease, RPL, our results suggest the need for further investigation of the role of FN3KRP in pregnancy loss.

MRKNs: Gene, Functions, and Role in Disease and Infection

The makorin RING finger protein (MKRN) gene family encodes proteins (makorins) with a characteristic array of zinc-finger motifs present in a wide array from invertebrates to vertebrates. MKRNs (MKRN1, MKRN2, MKRN3, MKRN4) as RING finger E3 ligases that mediate substrate degradation are related with conserved RING finger domains that control multiple cellular components via the ubiquitin-proteasome system (UPS), including p53, p21, FADD, PTEN, p65, Nptx1, GLK, and some viral or bacterial proteins. MKRNs also served as diverse roles in disease, like MKRN1 in transcription regulation, metabolic disorders, and tumors; MKRN2 in testis physiology, neurogenesis, apoptosis, and mutation of MKRN2 regulation signals transduction, inflammatory responses, melanoma, and neuroblastoma; MKRN3 in central precocious puberty (CPP) therapy; and MKRN4 firstly reported as a novel E3 ligase instead of a pseudogene to contribute to systemic lupus erythematosus (SLE). Here, we systematically review advances in the gene’s expression, function, and role of MKRNs orthologs in disease and pathogens infection. Further, MKRNs can be considered targets for the host’s innate intracellular antiviral defenses and disease therapy.

MicroRNA-136-5p protects cardiomyocytes from coronary microembolization through the inhibition of pyroptosis

This study investigated how miR-136-5p partially affected cardiomyocyte pyroptosis in rats with coronary microembolization (CME). The cardiac function and structure of rats with CME were evaluated using echocardiography, hematoxylin and eosin staining, Masson staining, and troponin I level. Pyroptosis was induced by lipopolysaccharide (LPS) in isolated rat cardiomyocytes and evaluated by the expression of caspase-1, NOD-like receptor family pyrin domain-containing 3, interleukin-1β, and gasdermin D-N. After cell transfection, the expression of Ataxin-1 like (ATXN1L), pyrin domain-containing 1 (PYDC1), and pyroptosis-related proteins was assessed. Dual-luciferase reporter and immunoprecipitation assays were used to verify the relationships among miR-136-5p, ATXN1L, and capicua (CIC). MiR-136-5p was under-expressed, whereas ATXN1L was overexpressed in rats with CME and in LPS-treated primary cardiomyocytes. MiR-136-5p targeted ATXN1L, and ATXN1L bound to CIC to suppress PYDC1 expression. MiR-136-5p overexpression suppressed pyroptosis by inhibiting the binding of ATXN1L with CIC and promoting PYDC1 expression, which was reversed by simultaneous elevation of ATXN1L. In conclusion, miR-136-5p suppressed pyroptosis by upregulating PYDC1 via ATXN1L/CIC axis, thereby attenuating cardiac damage caused by CME.

Physiological Roles of Mammalian Transmembrane Adenylyl Cyclase Isoforms

Adenylyl cyclases (ACs) catalyze the conversion of ATP to the ubiquitous second messenger cAMP. Mammals possess nine isoforms of transmembrane ACs, dubbed AC1-9, that serve as major effector enzymes of G protein-coupled receptors. The transmembrane ACs display varying expression patterns across tissues, giving potential for them having a wide array of physiologic roles. Cells express multiple AC isoforms, implying that ACs have redundant functions. Furthermore, all transmembrane ACs are activated by Gα_s so it was long assumed that all ACs are activated by Gα_s-coupled GPCRs. AC isoforms partition to different microdomains of the plasma membrane and form prearranged signaling complexes with specific GPCRs that contribute to cAMP signaling compartments. This compartmentation allows for a diversity of cellular and physiological responses by enabling unique signaling events to be triggered by different pools of cAMP. Isoform specific pharmacological activators or inhibitors are lacking for most ACs, making knockdown and overexpression the primary tools for examining the physiological roles of a given isoform. Much progress has been made in understanding the physiological effects mediated through individual transmembrane ACs. GPCR-AC-cAMP signaling pathways play significant roles in regulating functions of every cell and tissue, so understanding each AC isoform's role holds potential for uncovering new approaches for treating a vast array of pathophysiological conditions.

Association of Genetic Variants in miR-217 Gene with Risk of Coronary Artery Disease: A Case-Control Study

Objective

To evaluate the associations of genetic variants of the miR-217 gene with coronary artery disease (CAD) risk, as well as plasma level of vascular endothelial growth factor (VEGF).

Methods

A case-control study with 498 CAD patients and 499 frequency-matched healthy controls was conducted to evaluate the associations of four tagSNPs of the miR-217 gene, including rs6724872, rs4999828, rs10206823, and rs41291177, with CAD risk and plasma level of VEGF.

Results

SNP rs6724872 and rs4999828 were significantly associated with increased risk of CAD (P value was smaller than 0.05 even after Bonferroni multiple adjustment). Compared with the G allele, C allele of rs6724872 was significantly associated with 1.73-fold increased risk of CAD (95% CI: 1.25-2.39; P = 0.001). While C allele of rs4999828 was significantly associated with 1.75-fold increased risk of CAD, compared with T allele (95% CI: 1.34-2.29; P = 4 × 10^-5). Meanwhile, rs6724872 and rs4999828 were also significantly associated with higher level of VEGF (P < 0.001).

Conclusion

These findings highlighted the important role of genetic variants of the miR-217 gene in the pathogenesis of CAD and potential targets for intervention.

Oncogenic Landscape of Somatic Mutations Perturbing Pan-Cancer lncRNA-ceRNA Regulation

Competing endogenous RNAs (ceRNA) are transcripts that communicate with and co-regulate each other by competing for the binding of shared microRNAs (miRNAs). Long non-coding RNAs (lncRNAs) as a type of ceRNA constitute a competitive regulatory network determined by miRNA response elements (MREs). Mutations in lncRNA MREs destabilize their original regulatory pathways. Study of the effects of lncRNA somatic mutations on ceRNA mechanisms can clarify tumor mechanisms and contribute to the development of precision medicine. Here, we used somatic mutation profiles collected from TCGA to characterize the role of lncRNA somatic mutations in the ceRNA regulatory network in 33 cancers. The 31,560 mutation sites identified by TargetScan and miRanda affected the balance of 70,811 ceRNA regulatory pathways. Putative mutations were categorized as high or low based on mutation frequencies. Multivariate multiple regression revealed a significant effect of 162 high-frequency mutations in six cancer types on the expression levels of target mRNAs (ceMs) through the ceRNA mechanism. Low-frequency mutations in multiple cancers perturbing 1624 ceM have been verified by Student’s t-test, indicating a significant mechanism of changes in the expression level of oncogenic genes. Oncogenic signaling pathway studies involving ceMs indicated functional heterogeneity of multiple cancers. Furthermore, we identified that lncRNA, perturbing ceMs associated with patient survival, have potential as biomarkers. Our collective findings revealed individual differences in somatic mutations perturbing ceM expression and impacting tumor heterogeneity.

Hematological findings associated with tubulin-folding cofactors D-related encephalopathy: Expanding the phenotype

The dysfunction of microtubules (α/β-tubulin polymers) underlies a wide range of nervous system genetic abnormalities. Defects in TBCD, a tubulin-folding cofactor, cause diseases highlighted with early-onset encephalopathy with or without neurodegeneration, intellectual disability, seizures, microcephaly and tetraparaperesis. Utilizing various molecular methods, we describe nine patients from four unrelated families with two novel exon 18 variants in TBCD exhibiting the typical neurological phenotype of the disease. Interestingly, all the investigated patients had previously unreported hematological findings in the form of neutropenia and mild degree of anemia and thrombocytopenia. In addition to delineating the neurological phenotype in several patients with TBCD variants, our study stresses on the new association of neutropenia, in particular, with the disease.

Diagnostic utility of circulating plasma microRNA-101a in severity of coronary heart disease

BACKGROUND

For evaluating the severity of coronary heart disease (CHD), coronary arteriography may not be available everywhere due to technical limitations. MicroRNA-101a (miR-101a) associated with inflammation and cholesterol homeostasis. However, whether it related to presence and stratification of CHD is still unknown.

AIM

We aim to evaluate the value of miR-101a in stratifying CHD patients.

METHODS

We enrolled 200 CHD patients and 100 controls, and 200 CHD patients were divided into two groups of low and high SYNTAX score (SYNTAX score ≤ 22 versus SYNTAX score ≥ 33). Intergroup comparisons of miR-101a level were compared among the controls and two groups of low and high SYNTAX score. Correlation between miR-101a and blood lipid profiles was analyzed. The logistic regression analysis were conducted to evaluate the risk factors of CHD.

RESULTS

Relative level of miR-101a in the controls, SYNTAX score ≤ 22 and SYNTAX score ≥ 33 group were 4.61 (1.24-8.91), 3.28 (0.58-6.75) and 2.29 (1.04-3.62), respectively (p < 0.001). All lipid profiles significantly associated with miR-101a expression (all p < 0.001). The odds ratio (OR) of miR-101a in univariate analysis was 0.41 (95% CI, 0.33-0.52). After adjusting for the traditional risk factors, such as blood profiles and history of smoking, the odds ratio of miR-101a was 0.63 (95% CI, 0.47-0.43), which closely associated with CHD (p = 0.002).

CONCLUSIONS

Circulating miR-101a may be considered as a novel biomarker for evaluating the presence and severity of CHD.

Genetic association between CDKN2B/CDKN2B-AS1 gene polymorphisms with primary glaucoma in a North Indian cohort: an original study and an updated meta-analysis

Background

Variants in CDKN2B/CDKN2B-AS1 have been reported to modulate glaucoma risk in several GWAS across different populations. CDKN2B/CDKN2A encodes tumor suppressor proteins p16^INK4A/p15^INK4B which influences cell proliferation/senescence in RGCs, the degeneration of which is a risk factor for glaucoma. CDKN2B-AS1 codes a long non-coding RNA in antisense direction and is involved in influencing nearby CDKN2A/CDKN2B via regulatory mechanisms.

Methods

Current study investigated four SNPs (rs2157719, rs3217992, rs4977756, rs1063192) of aforementioned genes in a case–control study in a North Indian cohort. Genotyping was done with Taqman chemistry. In addition, an updated meta-analysis was performed.

Results

Two SNPs, rs3217992 and rs2157719 were found to be significantly associated with the disease. The frequency of ‘T’ allele of rs3217992 was significantly lower in cases (POAG/PACG) [p = 0.045; OR = 0.80(CI = 0.65–0.99) and p = 0.024; OR = 0.73(CI = 0.55–0.96)], respectively than in controls. Genetic model analysis revealed that TT + CT genotype confers 0.73-fold protection against POAG [p = 0.047; OR = 0.73(CI = 0.54–0.99)] and trend assumed additive model gives 0.53 times higher protection against PACG progression. However the association of rs3217992 with POAG and PACG did not remain significant after Bonferroni correction. For rs2157719, the ‘C’ allele was found to be less prevalent among cases (POAG/PACG) with respect to controls. Cochran Armitage trend test assuming additive model revealed 0.77 and 0.64-fold protection against POAG and PACG respectively. Bonferroni correction (p_corr= 0.003) was applied and the association of rs2157719 remained significant in PACG cases but not among POAG cases (p = 0.024). The ‘CC’ genotype also confers protection against primary glaucoma (POAG/PACG) among males and female subjects. The frequency rs1063192 and rs4977756 did not vary significantly among subjects, however the haplotype ‘CATA’ was found to be associated with increased glaucoma risk. An updated meta-analysis conducted on pooled studies on POAG cases and controls revealed significant association between rs1063192, rs2157719, rs4977756 and POAG except rs3217992.

Conclusion

The study concludes significant association between INK4 variants and primary glaucoma in the targeted North Indian Punjabi cohort. We believe that deep-sequencing of INK4 locus may help in identifying novel variants modifying susceptibility to glaucoma. Functional studies can further delineate the role of CDKN2B and CDKN2B-AS1 in primary glaucoma for therapeutic intervention.

FN3K expression in COPD: a potential comorbidity factor for cardiovascular disease

INTRODUCTION

Cigarette smoking and oxidative stress are common risk factors for the multi-morbidities associated with chronic obstructive pulmonary disease (COPD). Elevated levels of advanced glycation endproducts (AGE) increase the risk of cardiovascular disease (CVD) comorbidity and mortality. The enzyme fructosamine-3-kinase (FN3K) reduces this risk by lowering AGE levels.

METHODS

The distribution and expression of FN3K protein in lung tissues from stable COPD and control subjects, as well as an animal model of COPD, was assessed by immunohistochemistry. Serum FN3K protein and AGE levels were assessed by ELISA in patients with COPD exacerbations receiving metformin. Genetic variants within the FN3K and FN3K-RP genes were evaluated for associations with cardiorespiratory function in the Subpopulations and Intermediate Outcome Measures in COPD Study cohort.

RESULTS

This pilot study demonstrates that FN3K expression in the blood and human lung epithelium is distributed at either high or low levels irrespective of disease status. The percentage of lung epithelial cells expressing FN3K was higher in control smokers with normal lung function, but this induction was not observed in COPD patients nor in a smoking model of COPD. The top five nominal FN3K polymorphisms with possible association to decreased cardiorespiratory function (p<0.008-0.02), all failed to reach the threshold (p<0.0028) to be considered highly significant following multi-comparison analysis. Metformin enhanced systemic levels of FN3K in COPD subjects independent of their high-expression or low-expression status.

DISCUSSION

The data highlight that low and high FN3K expressors exist within our study cohort and metformin induces FN3K levels, highlighting a potential mechanism to reduce the risk of CVD comorbidity and mortality.

SNPs in Sites for DNA Methylation, Transcription Factor Binding, and miRNA Targets Leading to Allele-Specific Gene Expression and Contributing to Complex Disease Risk: A Systematic Review

INTRODUCTION

The complex genetic diversity among human populations results from an assortment of factors acting at various sequential levels, including mutations, population migrations, genetic drift, and selection. Although there are a plethora of DNA sequence variations identified through genome-wide association studies (GWAS), the challenge remains to explain the mechanisms underlying interindividual phenotypic disparity accounting for disease susceptibility. Single nucleotide polymorphisms (SNPs) present in the sites for DNA methylation, transcription factor (TF) binding, or miRNA targets can alter the gene expression. The systematic review aimed to evaluate the complex crosstalk among SNPs, miRNAs, DNA methylation, and TFs for complex multifactorial disease risk.

METHODS

PubMed and Scopus databases were used from inception until May 15, 2019. Initially, screening of articles involved studies assessing the interaction of SNPs with TFs, DNA methylation, or miRNAs resulting in allele-specific gene expression in complex multifactorial diseases. We also included the studies which provided experimental validation of the interaction of SNPs with each of these factors. The results from various studies on multifactorial diseases were assessed.

RESULTS

A total of 11 articles for SNPs interacting with DNA methylation, 30 articles for SNPs interacting with TFs, and 11 articles for SNPs in miRNA binding sites were selected. The interactions of SNPs with epigenetic factors were found to be implicated in different types of cancers, autoimmune diseases, cardiovascular diseases, diabetes, and asthma.

CONCLUSION

The systematic review provides evidence for the interplay between genetic and epigenetic risk factors through allele-specific gene expression in various complex multifactorial diseases.

Polymorphism analysis of miR182 and CDKN2B genes in Greek patients with primary open angle glaucoma

Glaucoma is a progressive optic neuropathy resulting from retinal ganglion cells death; it represents one of the leading causes of irreversible blindness worldwide. Although, primary open angle glaucoma (POAG) is the most common type of the disease, the pathogenesis of POAG and the genetic factors contributing to disease development remain poorly understood. The aim of this study was to investigate whether the polymorphisms rs76481776 in miR182 gene and rs3217992 in cyclin-dependent kinase inhibitor-2B (CDKN2B) gene are risk factors for POAG in a series of patients of Greek origin. A case-control study was conducted including 120 patients with POAG and 113 unaffected healthy controls of Greek origin, surveyed for polymorphisms with potential correlation to POAG. DNA from each individual was tested for the miR182 rs76481776 and CDKN2B rs3217992 polymorphisms. Regarding the miR182 rs76481776 polymorphism, the T allele occurred with significantly higher frequency in POAG patients compared to controls (OR: 2.62, 95% CI: 1.56-4.39; p = 0.0002). The CDKN2B rs3217992 A allele frequency was found significantly increased in POAG patients compared to healthy individuals (OR: 1.72, 95% CI: 1.18-2.49; p = 0.005). Therefore, both rs76481776 polymorphism in miR182 gene and rs3217992 polymorphism in CDKN2B gene seem to be associated with the development of POAG in a Greek population. The carriers of the T allele of rs76481776 in miR182 and the carriers of the A allele of rs3217992 in CDKN2B have an increased risk of developing POAG.

Overexpression of miR-375 Protects Cardiomyocyte Injury following Hypoxic-Reoxygenation Injury

The aim of the study was to evaluate the clinical significance of microRNA-375 in acute myocardial infarction patients and its mimic action in hypoxia/reoxygenation- (H/R-) induced ventricular cardiomyocyte H9c2 injury. In the current study, 90 ST-elevated acute MI patients (STEMI), 75 non-ST-elevated acute MI patients (NSTEMI), 90 healthy subjects, 14 weeks old mice, and ventricular cardiomyocyte H9c2 were included. The expressions of plasma microRNA-375 in patients with STEMI and NSTEMI and AMI mouse models were remarkably decreased than in controls (P < 0.001). The areas under the curve (AUC) of plasma microRNA-375 were revealed 0.939 in STEMI and 0.935 in NSTEMI subjects. Moreover, microRNA-375 levels in H/R-exposed cardiac H9c2 cells were evidently downregulated and significantly increased apoptosis rate and caspase-3 activity levels, while overexpression of miR-375 remarkably reduced apoptosis percentage and caspase-3 levels as compared with normal cells. Furthermore, this study also demonstrated that Nemo-like kinase (NLK), NLK mRNA, and protein expression levels were significantly downregulated in H/R-injured H9c2 cells, on the contrary, H9c2 cells transfected with mimic-miR-375 greatly upregulated NLK mRNA and protein expression. Plasma microRNA-375 may serve as an essential clinical biomarker for diagnosis of early-stage AMI. Mimic expression of miR-375 significantly prevented H/R-induced cardiomyocyte injury by decreasing caspase-3 activity through upregulation of the NLK gene, recommended as a new therapeutic option for AMI patient.

Circulating microRNA expression profiling and bioinformatics analysis of patients with coronary artery disease by RNA sequencing

BACKGROUND

MicroRNAs play a vital role in coronary artery disease. Abnormal expression of microRNAs has been found to be associated with the occurrence of CAD.

METHODS

We identified significantly differentially expressed microRNAs in plasma between 40 patients with CAD and 10 controls with NCA using RNA sequencing. The differentially expressed microRNAs were analyzed for Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment.

RESULTS

Fifty cDNA libraries were constructed and sequenced, and a total of 1871.82 M raw reads were obtained, and 2135 microRNAs were found. Compared to the expressed microRNAs of NCA controls, 159 microRNAs were differentially expressed in CAD patients, including 119 upregulated microRNAs and 40 downregulated microRNAs. The top 10 upregulated miRNAs were miR-144-3p, miR-34a-5p, miR-15b-3p, miR-22-3p, miR-29b-3p, miR-1270, miR-6891-5p, miR-106a-5p, miR-15b-5p, and hsa-miR-499b-3p. The top ten downregulated miRNAs were miR-4437, miR-6842-3p, miR-4664-3p, miR-671-3p, miR-219a-1-3p, miR-7848-3p, miR-664a-3p, miR-1284, miR-361-3p, and miR-6780a-5p. The target genes of differentially expressed microRNAs were related to many basic biological terms, such as biological process, cellular component, and molecular function. According to the KEGG pathway analysis, the most enriched pathways of the differentially expressed microRNAs were endocytosis, focal adhesion, axon guidance, and so on. Furthermore, six upregulated and two downregulated microRNAs were detected by qRT-PCR (Quantitative Real-time PCR) and ROC analysis for diagnosing CAD.

CONCLUSION

The results suggest that the expression levels of some microRNAs may play a vital role in the physiological and pathological course of CAD. Our study may provide useful information for the diagnosis and treatment of CAD.

MicroRNA-8073: Tumor suppressor and potential therapeutic treatment

The comprehensive screening of intracellular and extracellular microRNAs was performed to identify novel tumor suppressors. We found that miR-8073 was present in exosome and predominantly exported from colorectal cancer cells. Treatment with a synthetic miR-8073 mimic resulted in a dramatic decrease in the proliferation of various types of cancer cells, which was not observed in similarly treated normal cells. As little is known about the biological functions of miR-8073, its target mRNAs were analyzed by both mRNA expression and in silico sequence analyses, leading to five probable target candidates (FOXM1, MBD3, CCND1, KLK10, and CASP2) that enhance survival during the regulation of the cell cycle, cell proliferation, and apoptosis. We experimentally confirmed that miR-8073 binds the 3’-UTR of each of these mRNA target candidates and that the introduction of a synthetic miR-8073 mimic into cancer cells reduced levels of protein expression. Finally, the antiproliferative effects of miR-8073 were validated in vivo: the subcutaneous injection of a synthetic miR-8073 mimic suppressed colorectal tumor volume to 43% in tumor-bearing xenografted mice. These results suggest that because miR-8073 binds, and thus reduces the levels of, these oncogenic targets, cancer cells must actively downregulate miR-8073 as a survival mechanism. The introduction of miR-8073 into tumors could thus inhibit tumor growth, indicating its great potential for cancer therapeutics.

Cell Cycle Regulation of Stem Cells by MicroRNAs

MicroRNAs (miRNAs) are a class of small non-coding RNA molecules involved in the regulation of gene expression. They are involved in the fine-tuning of fundamental biological processes such as proliferation, differentiation, survival and apoptosis in many cell types. Emerging evidence suggests that miRNAs regulate critical pathways involved in stem cell function. Several miRNAs have been suggested to target transcripts that directly or indirectly coordinate the cell cycle progression of stem cells. Moreover, previous studies have shown that altered expression levels of miRNAs can contribute to pathological conditions, such as cancer, due to the loss of cell cycle regulation. However, the precise mechanism underlying miRNA-mediated regulation of cell cycle in stem cells is still incompletely understood. In this review, we discuss current knowledge of miRNAs regulatory role in cell cycle progression of stem cells. We describe how specific miRNAs may control cell cycle associated molecules and checkpoints in embryonic, somatic and cancer stem cells. We further outline how these miRNAs could be regulated to influence cell cycle progression in stem cells as a potential clinical application.

Evaluation of the role of CDKN2B gene in type 2 diabetes mellitus and hypertension in ethnic Saudi Arabs

Background

Coronary heart disease (CAD) is a multiple with several contributory risk traits, including type 2 diabetes and hypertension, which may share common genetic risk variants with the disease. Genome-wide association studies (GWASs) have yielded a wealth of information suggesting that CAD, the extent of contributory variants may differ according to genetic locus. The present study aimed at verifying whether the cyclin-dependent kinase 4 inhibitor B (CDKN2B) genomic region strongly associated with coronary artery disease (CAD)/myocardial infarction (MI) may also constitute risk for its risk factors type 2 diabetes mellitus (T2DM) and hypertension (HTN) in ethnic Saudi Arabs.

Methodology

We genotyped eight CDKN2B SNPs for cardiovascular risk in a total of 4650 Saudi Arabs, (3049 male and 1601 female) by Taqman assay. Of these individuals, 3732 had primary hypertension and 2576 had type 2 diabetes mellitus.

Results

Out of the eight studied SNPs, two, rs10757274_A [0.915 (0.840-1.00); p = 0.042], rs1333045_T [0.92(0.84-1.00); p = 0.048] were initially associated with type 2 diabetes but lost the association after multivariate adjustments for CAD, hypertension and MI, while rs10757274_A showed borderline association with hypertension.

Conclusions

Our finding does not support the notion of a critical role for the CDKN2B gene locus as a HTN or T2DM cardiovascular risk in ethnic Arabs. The study also demonstrates the importance of replication studies in ascertaining the role of a genomic sequence in disease.

The role of CDKN2B in cardiovascular risk in ethnic Saudi Arabs: A validation study

BACKGROUND

Genome-wide association studies (GWASs) have yielded a wealth of information furnishing support for the variability in genetic predisposition to disease. However, the actual impact of such findings on any particular ethnic population needs to be validated through replication studies. In the present study, we verified recent findings of a GWAS demonstrating a strong association for the cyclin-dependent kinase 4 inhibitor B (CDKN2B) genomic region with coronary artery disease (CAD)/myocardial infarction (MI) in ethnic Saudi Arabs.

METHODOLOGY

We genotyped 8 CDKN2B SNPs for cardiovascular risk in 4650 Saudi Arabs, comprising 2429 CAD cases (1860 males; 569 female) and 2221 controls (1189 male; 1032 female) by Taqman assay.

RESULTS

Four SNPs, rs4977574_A [0.56(0.50-0.63); p < 0.0001], rs10757274_A [0.87(0.77-0.97); p = 0.014], rs10738607_A [0.89(0.80-1.00); p = 0.043] and rs1333045_T [0.54(0.48-0.61); p < 0.0001] residing on the CDKN2B gene were significantly associated with CAD following multivariate adjustments for MI, HTN and DM, while four others were weakly associated with the disease. Likewise, three SNPs, rs1412829_G [0.84(0.72-0.97); p = 0.019], rs564398_C [0.81(0.70-0.94); p = 0.006], rs4977756_G [0.87(0.76-0.99); p = 0.036] were significantly associated with MI after multivariate adjustments for CAD, HTN and DM, while the other five displayed borderline associations.

CONCLUSIONS

Our findings strongly support the notion of a critical role for the CDKN2B gene locus as a cardiovascular risk in ethnic Arabs. The study also demonstrates the importance of replication studies in ascertaining the role of a genomic sequence in disease.

Integrative Analysis Identifies Genetic Variants Associated With Autoimmune Diseases Affecting Putative MicroRNA Binding Sites

Genome-wide and fine mapping studies have shown that more than 90% of genetic variants associated with autoimmune diseases (AID) are located in non-coding regions of the human genome and especially in regulatory sequences, including microRNAs (miRNA) target sites. MiRNAs are small endogenous noncoding RNAs that modulate gene expression at the post-transcriptional level. Single nucleotide polymorphisms (SNPs) located within the 3' untranslated region of their target mRNAs (miRSNP) can alter miRNA binding sites. Yet, little is known about their effect on regulation by miRNA and the consequences for AID. Conversely, it is well known that two or more AID may share part of their genetic background. Here, we hypothesized that miRSNPs could be associated with more than one AID. To identify miRSNPs associated with AID, we integrated results from three different prediction tools (Polymirts, miRSNP, and miRSNPscore) using a naïve Bayes classifier approach to identify miRSNPs predicted to affect binding sites of miRNAs. Further, to detect miRSNPs associated with two or more AID, we integrated predictions with summary statistics from 12 AID studies. In addition, to prioritize miRSNPs, miRNAs and AID-associated target genes, we used public expression quantitative trait locus (eQTL) data and mRNA-seq and small RNA-seq data. We identified 34 miRNSPs associated with at least two AID. Furthermore, we found 86 miRNAs predicted to target 18 of the associated gene's mRNAs. Our integrative approach revealed new insights into miRNAs and AID associated target genes. Thus, it helped to prioritize AID noncoding risk SNPs that might be involved in the causal mechanisms, providing valuable information for further functional studies.

Alteration in microRNA-155 level correspond to severity of coronary heart disease

Cardiovascular diseases are a consequence of genetic and epigenetic interactions. Inflammation contributes toward the initiation and progression of atherosclerotic lesions. Previous studies have shown that microRNA (miR) 155 plays a role in cardiovascular disease, including the prevention of inflammatory infiltration, regulation of autophagy, and participation of immunoreactions. However, the change of miR-155 level in the development of atherosclerosis remains to be determined. The initial objective of this study was that CHD patients would have altered serum miR-155 level. We also aim to identify whether circulating miR-155 content could be used as a predictor for severity of coronary atherosclerosis. Sample was collected from 300 CHD patients and 100 controls. Quantitative real-time PCR analysis was utilized on RNA isolated from plasma. Expression of miR-155 was identified on the basis of the quartiles of the Gensini score, and association between the microRNA and CHD was analyzed. CHD patients had higher miR-155 level in comparison to controls (p < .001), and the miRNA content significantly increased following an increasing Gensini score (p < .001). Gensini score was significantly associated with miR-155 expression (r = 0.6124, p < .001). Our findings suggest that interaction between circulating miR-155 expressions with classical risk factors of atherosclerotic lesions, and serum miR-155 content may serve as a novel biomarker for evaluating severity of CHD.

Life-Course Genome-wide Association Study Meta-analysis of Total Body BMD and Assessment of Age-Specific Effects

Bone mineral density (BMD) assessed by DXA is used to evaluate bone health. In children, total body (TB) measurements are commonly used; in older individuals, BMD at the lumbar spine (LS) and femoral neck (FN) is used to diagnose osteoporosis. To date, genetic variants in more than 60 loci have been identified as associated with BMD. To investigate the genetic determinants of TB-BMD variation along the life course and test for age-specific effects, we performed a meta-analysis of 30 genome-wide association studies (GWASs) of TB-BMD including 66,628 individuals overall and divided across five age strata, each spanning 15 years. We identified variants associated with TB-BMD at 80 loci, of which 36 have not been previously identified; overall, they explain approximately 10% of the TB-BMD variance when combining all age groups and influence the risk of fracture. Pathway and enrichment analysis of the association signals showed clustering within gene sets implicated in the regulation of cell growth and SMAD proteins, overexpressed in the musculoskeletal system, and enriched in enhancer and promoter regions. These findings reveal TB-BMD as a relevant trait for genetic studies of osteoporosis, enabling the identification of variants and pathways influencing different bone compartments. Only variants in ESR1 and close proximity to RANKL showed a clear effect dependency on age. This most likely indicates that the majority of genetic variants identified influence BMD early in life and that their effect can be captured throughout the life course.

New variant identified in major susceptibility locus to tuberculosis on chromosomal region 8q12-q13 in Moroccan population: a case control study

Background

Tuberculosis (TB) remains a global health problem. Several studies have implicated genetic host factors in predisposing populations to TB disease. In this study, we have selected NSMAF (Neutral Sphingomyelinase Activation Associated Factor) as a candidate gene to evaluate its level of association with TB disease in a Moroccan population for two reasons: first, this gene is located in a major susceptibility locus on chromosomal region 8q12-q13 in the Moroccan population, closely linked to the CYP7A1 gene, which was previously shown to be associated with TB disease; second, NSMAF has an important role in immune system function.

Methods

We conducted a case-control study including 269 genomic DNA samples extracted from pulmonary TB (PTB) patients and healthy controls (HC). We genotyped three selected SNPs (rs2228505, rs36067275 and rs10505004) using TaqMan® allelic discrimination assays.

Results

Only the rs1050504 C > T genotype was observed to be significantly associated with an increased risk for developing pulmonary TB (41.8% vs 27%, OR 1.95, 95% CI 1.16–3.27; p = 0.01). In contrast, the TT genotype was significantly associated with resistance to PTB (4.1% vs 15.6%, OR 0.23, 95% CI 0.08–0.63; p = 0.002).

Conclusion

Our findings suggest that genetic variations in the NSMAF gene could modulate the risk of PTB development in a Moroccan population. Further functional studies are needed to confirm these findings.

A Genome-Wide Scan for MicroRNA-Related Genetic Variants Associated With Primary Open-Angle Glaucoma

Purpose

To identify microRNAs (miRNAs) involved in primary open-angle glaucoma (POAG), using genetic data. MiRNAs are small noncoding RNAs that posttranscriptionally regulate gene expression. Genetic variants in miRNAs or miRNA-binding sites within gene 3'-untranslated regions (3'UTRs) are expected to affect miRNA function and contribute to disease risk.

Methods

Data from the recent genome-wide association studies on intraocular pressure, vertical cup-to-disc ratio (VCDR), cupa area and disc area were used to investigate the association of miRNAs with POAG endophenotypes. Putative targets of the associated miRNAs were studied according to their association with POAG and tested in cell line by transfection experiments for regulation by the miRNAs.

Results

Of 411 miRNA variants, rs12803915:A/G in the terminal loop of pre-miR-612 and rs2273626:A/C in the seed sequence of miR-4707 were significantly associated with VCDR and cup area (P values < 1.2 × 10-4). The first variant is demonstrated to increase the miR-612 expression. We showed that the second variant does not affect the miR-4707 biogenesis, but reduces the binding of miR-4707-3p to CARD10, a gene known to be involved in glaucoma. Moreover, of 72,052 miRNA-binding-site variants, 47 were significantly associated with four POAG endophenotypes (P value < 6.9 × 10-6). Of these, we highlighted 10 variants that are more likely to affect miRNA-mediated gene regulation in POAG. These include rs3217992 and rs1063192, which have been shown experimentally to affect miR-138-3p- and miR-323b-5p-mediated regulation of CDKN2B.

Conclusions

We identified a number of miRNAs that are associated with POAG endophenotypes. The identified miRNAs and their target genes are candidates for future studies on miRNA-related therapies for POAG.

The Rotterdam Study: 2018 update on objectives, design and main results

The Rotterdam Study is a prospective cohort study ongoing since 1990 in the city of Rotterdam in The Netherlands. The study targets cardiovascular, endocrine, hepatic, neurological, ophthalmic, psychiatric, dermatological, otolaryngological, locomotor, and respiratory diseases. As of 2008, 14,926 subjects aged 45 years or over comprise the Rotterdam Study cohort. Since 2016, the cohort is being expanded by persons aged 40 years and over. The findings of the Rotterdam Study have been presented in over 1500 research articles and reports (see www.erasmus-epidemiology.nl/rotterdamstudy ). This article gives the rationale of the study and its design. It also presents a summary of the major findings and an update of the objectives and methods.

Multi-omic network-based interrogation of rat liver metabolism following gastric bypass surgery featuring SWATH proteomics

Morbidly obese patients often elect for Roux-en-Y gastric bypass (RYGB), a form of bariatric surgery that triggers a remarkable 30% reduction in excess body weight and reversal of insulin resistance for those who are type II diabetic. A more complete understanding of the underlying molecular mechanisms that drive the complex metabolic reprogramming post-RYGB could lead to innovative non-invasive therapeutics that mimic the beneficial effects of the surgery, namely weight loss, achievement of glycemic control, or reversal of non-alcoholic steatohepatitis (NASH). To facilitate these discoveries, we hereby demonstrate the first multi-omic interrogation of a rodent RYGB model to reveal tissue-specific pathway modules implicated in the control of body weight regulation and energy homeostasis. In this study, we focus on and evaluate liver metabolism three months following RYGB in rats using both SWATH proteomics, a burgeoning label free approach using high resolution mass spectrometry to quantify protein levels in biological samples, as well as MRM metabolomics. The SWATH analysis enabled the quantification of 1378 proteins in liver tissue extracts, of which we report the significant down-regulation of Thrsp and Acot13 in RYGB as putative targets of lipid metabolism for weight loss. Furthermore, we develop a computational graph-based metabolic network module detection algorithm for the discovery of non-canonical pathways, or sub-networks, enriched with significantly elevated or depleted metabolites and proteins in RYGB-treated rat livers. The analysis revealed a network connection between the depleted protein Baat and the depleted metabolite taurine, corroborating the clinical observation that taurine-conjugated bile acid levels are perturbed post-RYGB.

HDL-cholesterol concentration in pregnant Chinese Han women of late second trimester associated with genetic variants in CETP, ABCA1, APOC3, and GALNT2

Objective

To investigate whether HDL-C level in pregnant Chinese Han women of late second trimester correlated with loci in high-density lipoprotein-cholesterol (HDL-C)-related genes found in genome-wide association studies (GWAS).

Methods

Seven single-nucleotide polymorphisms (rs3764261 in CETP, rs1532085 in LIPC, rs7241918 in LIPG, rs1883025 in ABCA1, rs4225 in APOC3, rs1059611 in LPL, and rs16851339 in GALNT2) were genotyped using the Sequenom MassArray system for 1,884 pregnant women.

Results

The following polymorphisms were statistically associated with HDL-C level after adjusting for age, gestational week, pre-pregnancy BMI and state of GDM or HOMAIR: (i) rs3764261 (b = -0.055 mmol/L, 95% CI -0.101 to -0.008, p = 0.021), (ii) rs1883025 (b = -0.054 mmol/L, 95% CI -0.097 to -0.012, p = 0.013), (iii) rs4225 (b = -0.071 mmol/L, 95% CI -0.116 to -0.027, p = 1.79E-3) and (iv) rs16851339 (b = -0.064 mmol/L, 95% CI -0.120 to -0.008, p = 0.025). The more risk alleles the pregnant women have, the lower the plasma HDL-C levels of the subjects are.

Conclusions

Several risk alleles found to be related to HDL-C in GWAS are also associated with HDL-C levels in pregnant Chinese Han women and these risk loci contribute additively to low HDL-C levels.

Functional annotation of Alzheimer's disease associated loci revealed by GWASs

Genome-wide association studies (GWASs) discovered a number of SNPs and genes associated with Alzheimer's disease (AD). However, how these SNPs and genes influence the liability to AD is not fully understood. We deployed computational approaches to explore the function and action mechanisms of AD -related SNPs and genes identified by GWASs, including the effects of 195 GWAS lead SNPs and 338 proxy SNPs on miRNAs binding and protein phosphorylation, their RegulomeDB and 3DSNP scores, and gene ontology, pathway enrichment and protein-protein interaction network of 126 AD-associated genes. Our computational analysis identified 6 lead SNPs (rs10119, rs1048699, rs148763909, rs610932, rs6857 and rs714948) and 2 proxy SNPs (rs12539172 and rs2847655) that potentially impacted the miRNA binding. Lead SNP rs2296160 and proxy SNPs rs679620 and rs2228145 were identified as PhosSNPs potentially influencing protein phosphorylation. AD-associated genes showed enrichment of “regulation of beta-amyloid formation”, “regulation of neurofibrillary tangle assembly”, “leukocyte mediated immunity” and “protein-lipid complex assembly” signaling pathway. Protein-protein interaction network and functional module analyses identified highly-interconnected “hub” genes (APOE, PICALM, BIN1, ABCA7, CD2AP, CLU, CR1, MS4A4E and MS4A6A) and bottleneck genes (APOE, TOMM40, NME8, PICALM, CD2AP, ZCWPW1, FAM180B, GAB2 and PTK2B) that created three tight subnetworks. Our results provided the targets for further experimental assessment and further insight on AD pathophysiology.

The roles of RNA processing in translating genotype to phenotype

A goal of human genetics studies is to determine the mechanisms by which genetic variation produces phenotypic differences that affect human health. Efforts in this respect have previously focused on genetic variants that affect mRNA levels by altering epigenetic and transcriptional regulation. Recent studies show that genetic variants that affect RNA processing are at least equally as common as, and are largely independent from, those variants that affect transcription. We highlight the impact of genetic variation on pre-mRNA splicing and polyadenylation, and on the stability, translation and structure of mRNAs as mechanisms that produce phenotypic traits. These results emphasize the importance of including RNA processing signals in analyses to identify functional variants.

The roles of RNA processing in translating genotype to phenotype

A goal of human genetics studies is to determine the mechanisms by which genetic variation produces phenotypic differences that affect human health. Efforts in this respect have previously focused on genetic variants that affect mRNA levels by altering epigenetic and transcriptional regulation. Recent studies show that genetic variants that affect RNA processing are at least equally as common as, and are largely independent from, those variants that affect transcription. We highlight the impact of genetic variation on pre-mRNA splicing and polyadenylation, and on the stability, translation and structure of mRNAs as mechanisms that produce phenotypic traits. These results emphasize the importance of including RNA processing signals in analyses to identify functional variants.

Functional single nucleotide polymorphisms within the cyclin-dependent kinase inhibitor 2A/2B region affect pancreatic cancer risk

The CDKN2A (p16) gene plays a key role in pancreatic cancer etiology. It is one of the most commonly somatically mutated genes in pancreatic cancer, rare germline mutations have been found to be associated with increased risk of developing familiar pancreatic cancer and CDKN2A promoter hyper-methylation has been suggested to play a critical role both in pancreatic cancer onset and prognosis. In addition several unrelated SNPs in the 9p21.3 region, that includes the CDNK2A, CDNK2B and the CDNK2B-AS1 genes, are associated with the development of cancer in various organs. However, association between the common genetic variability in this region and pancreatic cancer risk is not clearly understood. We sought to fill this gap in a case-control study genotyping 13 single nucleotide polymorphisms (SNPs) in 2,857 pancreatic ductal adenocarcinoma (PDAC) patients and 6,111 controls in the context of the Pancreatic Disease Research (PANDoRA) consortium. We found that the A allele of the rs3217992 SNP was associated with an increased pancreatic cancer risk (ORhet=1.14, 95% CI 1.01-1.27, p=0.026, ORhom=1.30, 95% CI 1.12-1.51, p=0.00049). This pleiotropic variant is reported to be a mir-SNP that, by changing the binding site of one or more miRNAs, could influence the normal cell cycle progression and in turn increase PDAC risk. In conclusion, we observed a novel association in a pleiotropic region that has been found to be of key relevance in the susceptibility to various types of cancer and diabetes suggesting that the CDKN2A/B locus could represent a genetic link between diabetes and pancreatic cancer risk.

Genome-wide identification of microRNA-related variants associated with risk of Alzheimer's disease

MicroRNAs (miRNAs) serve as key post-transcriptional regulators of gene expression. Genetic variation in miRNAs and miRNA-binding sites may affect miRNA function and contribute to disease risk. Here, we investigated the extent to which variants within miRNA-related sequences could constitute a part of the functional variants involved in developing Alzheimer’s disease (AD), using the largest available genome-wide association study of AD. First, among 237 variants in miRNAs, we found rs2291418 in the miR-1229 precursor to be significantly associated with AD (p-value = 6.8 × 10⁻⁵, OR = 1.2). Our in-silico analysis and in-vitro miRNA expression experiments demonstrated that the variant’s mutant allele enhances the production of miR-1229-3p. Next, we found miR-1229-3p target genes that are associated with AD and might mediate the miRNA function. We demonstrated that miR-1229-3p directly controls the expression of its top AD-associated target gene (SORL1) using luciferase reporter assays. Additionally, we showed that miR-1229-3p and SORL1 are both expressed in the human brain. Second, among 42,855 variants in miRNA-binding sites, we identified 10 variants (in the 3′ UTR of 9 genes) that are significantly associated with AD, including rs6857 that increases the miR-320e-mediated regulation of PVRL2. Collectively, this study shows that miRNA-related variants are associated with AD and suggests miRNA-dependent regulation of several AD genes.

Analysis of the Relationship Between ADIPOR1 Variants and the Susceptibility of Chronic Metabolic Diseases in a Northeast Han Chinese Population

Objective: Shared genetic variants in ADIPOR1 have been identified as closely related to coronary artery disease (CAD), type 2 diabetes (T2D), and T2D with CAD susceptibility, suggesting that these variants are strong candidates for the common soil hypothesis. Therefore, it is essential to analyze the relationship between ADIPOR1 variants and the susceptibility to CAD, T2D, and T2D with CAD in other populations.

Materials and Methods: A case–control study was conducted which included three case cohorts [CAD (n = 316), T2D (n = 295), T2D with CAD (n = 302)], and a control cohort (n = 268) from a population in northeast China. Six ADIPOR1 single-nucleotide polymorphisms were genotyped by high-resolution melting and polymerase chain reaction–restriction fragment length polymorphism.

Results: We confirmed that the shared variant, rs3737884*G, in ADIPOR1 is associated with CAD, T2D, and T2D with CAD (p-value range: 6.54E-6–1.82E-5, odds ratio [OR] range: 1.770–1.844) and that rs16850797*C is associated with T2D and T2D with CAD (p-value range: 0.001–0.001, OR range: 1.529–1.571). We also found that a novel shared variant, rs7514221*C, is associated with an increased susceptibility to CAD, T2D, and T2D with CAD (p-value range: 0.002–0.004, OR range: 1.194–2.382) in this population.

Conclusions:ADPOR1 variants, rs3737884*G and rs7514221*C, may be shared risk factors associated with CAD, T2D, and T2D with CAD in a population of northeast China.

A meta-analysis of 120 246 individuals identifies 18 new loci for fibrinogen concentration

Genome-wide association studies have previously identified 23 genetic loci associated with circulating fibrinogen concentration. These studies used HapMap imputation and did not examine the X-chromosome. 1000 Genomes imputation provides better coverage of uncommon variants, and includes indels. We conducted a genome-wide association analysis of 34 studies imputed to the 1000 Genomes Project reference panel and including ∼120 000 participants of European ancestry (95 806 participants with data on the X-chromosome). Approximately 10.7 million single-nucleotide polymorphisms and 1.2 million indels were examined. We identified 41 genome-wide significant fibrinogen loci; of which, 18 were newly identified. There were no genome-wide significant signals on the X-chromosome. The lead variants of five significant loci were indels. We further identified six additional independent signals, including three rare variants, at two previously characterized loci: FGB and IRF1. Together the 41 loci explain 3% of the variance in plasma fibrinogen concentration.

The Rotterdam Study: 2016 objectives and design update

The Rotterdam Study is a prospective cohort study ongoing since 1990 in the city of Rotterdam in The Netherlands. The study targets cardiovascular, endocrine, hepatic, neurological, ophthalmic, psychiatric, dermatological, otolaryngological, locomotor, and respiratory diseases. As of 2008, 14,926 subjects aged 45 years or over comprise the Rotterdam Study cohort. The findings of the Rotterdam Study have been presented in over 1200 research articles and reports (see www.erasmus-epidemiology.nl/rotterdamstudy ). This article gives the rationale of the study and its design. It also presents a summary of the major findings and an update of the objectives and methods.