Allelic Imbalance in Regulation of ANRIL through Chromatin Interaction at 9p21 Endometriosis Risk Locus

An Unanticipated Modulation of Cyclin-Dependent Kinase Inhibitors: The Role of Long Non-Coding RNAs

It is now definitively established that a large part of the human genome is transcribed. However, only a scarce percentage of the transcriptome (about 1.2%) consists of RNAs that are translated into proteins, while the large majority of transcripts include a variety of RNA families with different dimensions and functions. Within this heterogeneous RNA world, a significant fraction consists of sequences with a length of more than 200 bases that form the so-called long non-coding RNA family. The functions of long non-coding RNAs range from the regulation of gene transcription to the changes in DNA topology and nucleosome modification and structural organization, to paraspeckle formation and cellular organelles maturation. This review is focused on the role of long non-coding RNAs as regulators of cyclin-dependent kinase inhibitors’ (CDKIs) levels and activities. Cyclin-dependent kinases are enzymes necessary for the tuned progression of the cell division cycle. The control of their activity takes place at various levels. Among these, interaction with CDKIs is a vital mechanism. Through CDKI modulation, long non-coding RNAs implement control over cellular physiology and are associated with numerous pathologies. However, although there are robust data in the literature, the role of long non-coding RNAs in the modulation of CDKIs appears to still be underestimated, as well as their importance in cell proliferation control.

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

APOBEC mediated mutagenesis drives genomic heterogeneity in endometriosis

Endometriosis is a benign gynecologic condition, acting as a precursor of certain histological subtypes of ovarian cancers. The epithelial cells of endometriotic tissues and normal uterine endometrium accumulated somatic mutations in cancer-associated genes such as phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) and Kirsten rat sarcoma (KRAS) proto-oncogene. To determine the genomic characteristic of endometriotic epithelial cells and normal uterine endometrium and to identify the predominant mutational process acting on them, we studied the somatic mutation profiles obtained from whole exome sequencing of 14 endometriotic epithelium and 11 normal uterine endometrium tissues and classified them into mutational signatures. We observed that single base substitutions 2/13 (SBS), attributed to Apolipoprotein B mRNA Editing Enzyme Catalytic Subunit (APOBEC) induced mutagenesis, were significant in endometriotic tissues, but not in the normal uterine endometrium. Additionally, the larger number and wider allele frequency distribution of APOBEC signature mutations, compared to cancer-associated driver mutations in endometriotic epithelium suggested APOBEC mutagenesis as an important source of mutational burden and heterogeneity in endometriosis. Further, the relative risk of enriched APOBEC signature mutations was higher in endometriosis patients who were carriers of APOBEC3A/3B germline deletion, a common polymorphism in East Asians which involves the complete loss of APOBEC3B coding region. Our results illustrate the significance of APOBEC induced mutagenesis in driving the genomic heterogeneity of endometriosis.

Chronic lymphocytic leukemia (CLL) risk is mediated by multiple enhancer variants within CLL risk loci

Chronic lymphocytic leukemia (CLL) is the most common adult leukemia in Western countries. It has a strong genetic basis, showing a ~ 8-fold increased risk of CLL in first-degree relatives. Genome-wide association studies (GWAS) have identified 41 risk variants across 41 loci. However, for a majority of the loci, the functional variants and the mechanisms underlying their causal roles remain undefined. Here, we examined the genetic and epigenetic features associated with 12 index variants, along with any correlated (r2 ≥ 0.5) variants, at the CLL risk loci located outside of gene promoters. Based on publicly available ChIP-seq and chromatin accessibility data as well as our own ChIP-seq data from CLL patients, we identified six candidate functional variants at six loci and at least two candidate functional variants at each of the remaining six loci. The functional variants are predominantly located within enhancers or super-enhancers, including bi-directionally transcribed enhancers, which are often restricted to immune cell types. Furthermore, we found that, at 78% of the functional variants, the alternative alleles altered the transcription factor binding motifs or histone modifications, indicating the involvement of these variants in the change of local chromatin state. Finally, the enhancers carrying functional variants physically interacted with genes enriched in the type I interferon signaling pathway, apoptosis, or TP53 network that are known to play key roles in CLL. These results support the regulatory roles for inherited noncoding variants in the pathogenesis of CLL.

Epigenetic Regulation of the Vascular Endothelium by Angiogenic LncRNAs

The functional properties of the vascular endothelium are diverse and heterogeneous between vascular beds. This is especially evident when new blood vessels develop from a pre-existing closed cardiovascular system, a process termed angiogenesis. Endothelial cells are key drivers of angiogenesis as they undergo a highly choreographed cascade of events that has both exogenous (e.g., hypoxia and VEGF) and endogenous regulatory inputs. Not surprisingly, angiogenesis is critical in health and disease. Diverse therapeutics target proteins involved in coordinating angiogenesis with varying degrees of efficacy. It is of great interest that recent work on non-coding RNAs, especially long non-coding RNAs (lncRNAs), indicates that they are also important regulators of the gene expression paradigms that underpin this cellular cascade. The protean effects of lncRNAs are dependent, in part, on their subcellular localization. For instance, lncRNAs enriched in the nucleus can act as epigenetic modifiers of gene expression in the vascular endothelium. Of great interest to genetic disease, they are undergoing rapid evolution and show extensive inter- and intra-species heterogeneity. In this review, we describe endothelial-enriched lncRNAs that have robust effects in angiogenesis.

Genomics of Endometriosis: From Genome Wide Association Studies to Exome Sequencing

This review aims at better understanding the genetics of endometriosis. Endometriosis is a frequent feminine disease, affecting up to 10% of women, and characterized by pain and infertility. In the most accepted hypothesis, endometriosis is caused by the implantation of uterine tissue at ectopic abdominal places, originating from retrograde menses. Despite the obvious genetic complexity of the disease, analysis of sibs has allowed heritability estimation of endometriosis at ~50%. From 2010, large Genome Wide Association Studies (GWAS), aimed at identifying the genes and loci underlying this genetic determinism. Some of these loci were confirmed in other populations and replication studies, some new loci were also found through meta-analyses using pooled samples. For two loci on chromosomes 1 (near CCD42) and chromosome 9 (near CDKN2A), functional explanations of the SNP (Single Nucleotide Polymorphism) effects have been more thoroughly studied. While a handful of chromosome regions and genes have clearly been identified and statistically demonstrated as at-risk for the disease, only a small part of the heritability is explained (missing heritability). Some attempts of exome sequencing started to identify additional genes from families or populations, but are still scarce. The solution may reside inside a combined effort: increasing the size of the GWAS designs, better categorize the clinical forms of the disease before analyzing genome-wide polymorphisms, and generalizing exome sequencing ventures. We try here to provide a vision of what we have and what we should obtain to completely elucidate the genetics of this complex disease.

The lncRNA ANRIL Gene rs2151280 GG Genotype is Associated with Increased Susceptibility to Recurrent Miscarriage in a Southern Chinese Population

Background

Genetic factors may play an important role in susceptibility to recurrent miscarriage. Some cardiovascular disease-related candidate genes have been shown to be associated with recurrent miscarriage. Long noncoding RNA ANRIL has been confirmed to be associated with susceptibility to various diseases, such as cardiovascular disease. However, it remains unclear whether the ANRIL gene polymorphism is related to recurrent miscarriage susceptibility.

Methods

Three ANRIL gene polymorphisms (rs2151280, rs1063192 and rs564398) were genotyped in 819 controls and 610 recurrent miscarriage patients through TaqMan real-time polymerase chain reaction. The odds ratios and 95% confidence intervals (CIs) were used to assess the strength of each association.

Results

Our results showed that the ANRIL rs2151280 GG genotype was associated with increased susceptibility to recurrent miscarriage (GG vs AA: adjusted OR=1.527, 95% CI=1.051–2.218, p=0.0262; GG vs AG/AA adjusted OR=1.460, 95% CI=1.021–2.089, p=0.0381). By combining the analysis of the risk genotypes in the three SNPs, we found that individuals with 2–3 risk genotypes had a significantly increased risk of recurrent miscarriage compared with those with a 0–1 risk genotype (adjusted OR=1.728, 95% CI=1.112–2.683, p=0.0149). This risk was more significant in subgroups of women less than 35–40 years of age and women with 2–3 miscarriages.

Conclusion

These results suggested that a specific SNP in the ANRIL gene may be associated with increased susceptibility to recurrent miscarriage in a southern Chinese population.

Integrative analysis of liver-specific non-coding regulatory SNPs associated with the risk of coronary artery disease

Genetic factors underlying coronary artery disease (CAD) have been widely studied using genome-wide association studies (GWASs). However, the functional understanding of the CAD loci has been limited by the fact that a majority of GWAS variants are located within non-coding regions with no functional role. High cholesterol and dysregulation of the liver metabolism such as non-alcoholic fatty liver disease confer an increased risk of CAD. Here, we studied the function of non-coding single-nucleotide polymorphisms in CAD GWAS loci located within liver-specific enhancer elements by identifying their potential target genes using liver cis-eQTL analysis and promoter Capture Hi-C in HepG2 cells. Altogether, 734 target genes were identified of which 121 exhibited correlations to liver-related traits. To identify potentially causal regulatory SNPs, the allele-specific enhancer activity was analyzed by (1) sequence-based computational predictions, (2) quantification of allele-specific transcription factor binding, and (3) STARR-seq massively parallel reporter assay. Altogether, our analysis identified 1,277 unique SNPs that display allele-specific regulatory activity. Among these, susceptibility enhancers near important cholesterol homeostasis genes (APOB, APOC1, APOE, and LIPA) were identified, suggesting that altered gene regulatory activity could represent another way by which genetic variation regulates serum lipoprotein levels. Using CRISPR-based perturbation, we demonstrate how the deletion/activation of a single enhancer leads to changes in the expression of many target genes located in a shared chromatin interaction domain. Our integrative genomics approach represents a comprehensive effort in identifying putative causal regulatory regions and target genes that could predispose to clinical manifestation of CAD by affecting liver function.

Effect of indel variants within the sorting nexin 29 (SNX29) gene on growth traits of goats

The sorting nexin 29 gene (SNX29) is a well-known regulator of myocyte differentiation and proliferation. In this work, two indels (17-bp and 21-bp) were identified in the goat SNX29 gene, and their effects on the growth traits of 1,759 Shaanbei white cashmere (SBWC) goats were analyzed. Both indels had three genotypes [homozygote wild type (II), heterozygote (ID), and homozygote mutation (DD)] and displayed medium genetic diversity (0.25 < polymorphism information content (PIC) < 0.50) in the population. The 17-bp indel was significantly associated with chest width (p = 0.009), body weight (p = 0.021), and chest depth (p = 0.032), with the II genotype dominant. The 21-bp indel was significantly associated with chest width (p = 0.001), chest depth (p = 4.8E-5), heart girth (p = 0.007), and hip width (p = 0.002). Because the two indels were in the upstream (17-bp) and intron (21-bp) regions of the SNX29 gene, transcription factor binding sites were predicted. The IRF5 and MYC could bind with the 17-bp indel and 21-bp indel sequences, respectively. This study indicates that SNX29 is a promising candidate gene that can be used to improve meat production in goat breeding.

Should Genetics Now Be Considered the Pre-eminent Etiologic Factor in Endometriosis?

Endometriosis is a common, estrogen-dependent, inflammatory disorder characterized by the growth of endometrial-like tissue at extrauterine locations. Its pathogenesis and mechanisms underlying its pathophysiology are poorly understood, although genetic variation is strongly implicated in these processes. Genetic studies reveal that approximately 50% of risk for endometriosis is due to genetic factors and the other 50% likely owing to environmental factors. As with other complex diseases, genetic variants in the DNA sequence increasing endometriosis risk all have small effects, unlike most single-gene disorders. It is the combinations of these variants adding together that contribute to higher risks for individual women. In addition, recent data on disease lesions demonstrate a high frequency of somatic (likely acquired) mutations, some of which are present in the eutopic endometrium and specifically in the epithelial cell compartment, raising the possibility that abnormal epithelial progenitors in the eutopic endometrium give rise to ectopic disease. Discovery in this field is occurring at a rapid pace, and further definitions of genetic (germline) and environmental (somatic) contributions to the pathogenesis and pathophysiology of this disorder are anticipated soon. These discoveries are expected to increase diagnostic, therapeutic, and preventive strategies to minimize disease and its associated morbidities.

Pig Hsd17b3: Alternative splice variants expression, insertion/deletion (indel) in promoter region and their associations with male reproductive traits

Hydroxysteroid 17-Beta Dehydrogenase 3 (Hsd17b3), primarily expressed in Leydig cells (LCs) of the mammalian testes, is essential for testosterone biosynthesis and male fertility. The aim of our study was to profile the expression, splice variants (SV) and novel insertion/deletion (indel) of Hsd17b3 in boars. Quantitative analysis showed that the expression level of Hsd17b3 in the testis was significantly highest. Among different testicular cell types, the Hsd17b3 mRNA expression level of LCs was significantly higher than that of SSCs (spermatogonial stem cells) and SCs (Sertoli cells). Furthermore, the SV was firstly identified in pigs and it was highly expressed in LCs comparing with SSCs and SCs. In addition, two mutations were identified in pig Hsd17b3 gene promotor and intron, respectively, which were associated with male reproductive traits (P <  0.05). In conclusion, both transcripts of Hsd17b3 gene were highly expressed in pig testes and LCs; the two novel indel variants of Hsd17b3 gene can be used as potential DNA makers for the marker-assisted selection in pigs. All these findings would enrich the study of Hsd17b3 gene in pig genetic breeding.

Clonal Expansion and Diversification of Cancer-Associated Mutations in Endometriosis and Normal Endometrium

Endometriosis is characterized by ectopic endometrial-like epithelium and stroma, of which molecular characteristics remain to be fully elucidated. We sequenced 107 ovarian endometriotic and 82 normal uterine endometrial epithelium samples isolated by laser microdissection. In both endometriotic and normal epithelium samples, numerous somatic mutations were identified within genes frequently mutated in endometriosis-associated ovarian cancers. KRAS is frequently mutated in endometriotic epithelium, with a higher mutant allele frequency (MAF) accompanied by arm-level allelic imbalances. Analyses of MAF, combined with multiregional sequencing, illuminated spatiotemporal evolution of the endometriosis and uterine endometrium genomes. We sequenced 109 single endometrial glands and found that each gland carried distinct cancer-associated mutations, demonstrating the heterogeneity of the genomic architecture of endometrial epithelium. Remarkable increases in MAF of mutations in cancer-associated genes in endometriotic epithelium suggest retrograde flow of endometrial cells already harboring cancer-associated mutations, with selective advantages at ectopic sites, leading to the development of endometriosis.

Unveiling the gene regulatory landscape in diseases through the identification of DNase I-hypersensitive sites

DNase I-hypersensitive sites (DHSs) serve key roles in the regulation of gene transcription as markers of cis-regulatory elements (CREs). Recent advances in next-generation sequencing have enabled the genome-wide location and annotation of DHSs in a variety of cells. Numerous studies have confirmed that DHSs are involved in several processes in cell fate decision and development. DHSs have also been indicated in cancer and inherited diseases as driver distal regulatory elements. Here, the definition of DHSs is reviewed, in addition to high-throughput methods of DHS identification. Furthermore, the function of DHSs in gene expression is probed. The roles of DHSs in disease occurrence are also reviewed and discussed. Concomitant advances in the identification of essential roles of DHSs will assist in disclosing the underlying molecular mechanisms, supplementing gene transcription and enlarging the molecular basis of DHS-related bioprocesses, phenotypes, distinct traits and diseases.

Pig SOX9: Expression profiles of Sertoli cell (SCs) and a functional 18 bp indel affecting testis weight

Sex determining region Y-box 9 (SOX9), an important member of the SRY- type HMGbox (SOX) gene family, plays an important role in the regulation of mammalian reproduction, including sex differentiation during the embryonic development stage and spermatogenesis after birth. To explore the roles of polymorphism and expression of the SOX9 gene in the development of testes, we analyzed the indel of SOX9 in pigs and the corresponding expression level of the SOX9 gene in 7-day and 5-month-old porcine Sertoli cells. Results revealed that the DD haplotype of SOX9 gene as well as the ID genotype were significantly associated with larger testicular weight, while the II haplotype was closely related to the smaller testicular weight. More importantly, the SOX9 gene expression of ID genotyped group was significantly higher than that in II genotyped group. Our results first revealed that the indel polymorphism and expression of SOX9 were significantly associated with pig reproduction traits indicating the critical roles of SOX9 gene in testes development. The study provides a new clue for understanding the regulation of animal reproductive activities.

Comprehensive analysis of circular RNAs in pathological states: biogenesis, cellular regulation, and therapeutic relevance

Circular RNAs (circRNAs) are members of the non-coding transcriptome; however, some of them are translated into proteins. These transcripts have important roles in both physiological and pathological mechanisms due to their ability to directly influence cellular signaling pathways. Specifically, circRNAs are regulators of transcription, translation, protein interaction, and signal transduction. An increased knowledge within their area is observed over the last few years, concomitant with the development of next-generation sequencing techniques. circRNAs are mostly tissue and disease specific with the ability of specifically changing the biological behavior of cells. The altered expression profile is currently investigated as novel minimally invasive diagnosis/prognosis tool and also therapeutic target in human disease. The diagnosis approach is based on their level modification within pathological states, especially cancer, where circRNAs' therapies are intensively explored in anti-aging strategies, diabetes, cardiovascular diseases, and malignant pathologies, and are relying on the restoration of homeostatic profiles.

Defining the genetic profile of endometriosis

Endometriosis is a pathological condition which has been extensively studied, since its pathophysiology stems from a broad spectrum of environmental influences and genetic factors. Familial studies aim at defining inheritance trends, while linkage analysis studies focus on the identification of genetic sites related to endometriosis susceptibility. Genetic association studies take into account candidate genes and single nucleotide polymorphisms, and hence target at unraveling the association between disease severity and genetic variation. The common goal of various types of studies is, through genetic mapping methods, the timely identification of therapeutic strategies for disease symptoms, including pelvic pain and infertility, as well as efficient counselling. While genome-wide association studies (GWAS) play a primary role in depicting genetic contributions to disease development, they entail a certain bias as regards the case-control nature of their design and the reproducibility of the results. Nevertheless, genetic-oriented studies and the implementation of the results through clinical tests, hold a considerable advantage in proper disease management. In this review article, we present information about gene-gene and gene-environment interactions involved in endometriosis and discuss the effectiveness of GWAS in identitying novel potential therapeutic targets in an attempt to develop novel therapeutic strategies for a better management and treatment of patients with endometriosis.

Goat PDGFRB: unique mRNA expression profile in gonad and significant association between genetic variation and litter size

β-Type platelet-derived growth factor receptor (PDGFRB) is a typical tyrosine kinase, as a candidate gene associated with reproduction. Its main roles include regulation of gonocytes (migration and proliferation) and of the cell cycle. The objectives of this study were to identify mRNA expression of the goat PDGFRB gene, as well as insertion/deletion (indel) variants and their association with litter size in 1122 healthy Shaanbei white cashmere goats. The results revealed that PDGFRB was widely expressed in all tested tissues, and the expression levels in testes at different developmental stages indicated a potential association with the mitosis-to-meiosis transition. Furthermore, the expression of PDGFRB was relatively higher in the ovary tissue of mothers of two lambs compared with mothers of single lamb. These results implied that PDGFRB was related to goat fertility. Meanwhile, two intronic indels, 5 bp (n = 501) and 10 bp (n = 1122), were identified. Statistical analysis revealed that only the 10 bp indel was associated with first-born litter size (n = 1122, p = 6.030 × 10^-5), and that individuals of the genotype insertion/deletion had larger litter sizes than those of genotype insertion/insertion. Overall, these results indicated that the 10 bp indel of PDGFRB could be used in marker-assisted selection during goat genetic breeding.

Long Noncoding RNA ANRIL: Lnc-ing Genetic Variation at the Chromosome 9p21 Locus to Molecular Mechanisms of Atherosclerosis

Ever since the first genome-wide association studies (GWAS) on coronary artery disease (CAD), the Chr9p21 risk locus has emerged as a top signal in GWAS of atherosclerotic cardiovascular disease, including stroke and peripheral artery disease. The CAD risk SNPs on Chr9p21 lie within a stretch of 58 kilobases of non-protein-coding DNA, containing the gene body of the long noncoding RNA (lncRNA) antisense non coding RNA in the INK4 locus (ANRIL). How risk is affected by the Chr9p21 locus in molecular detail is a matter of ongoing research. Here we will review recent advances in the understanding that ANRIL serves as a key risk effector molecule of atherogenesis at the locus. One focus of this review is the shift in understanding that genetic variation at Chr9p21 not only affects the abundance of ANRIL, and in some cases expression of the adjacent CDKN2A/B tumor suppressors, but also impacts ANRIL splicing, such that 3′-5′-linked circular noncoding ANRIL RNA species are produced. We describe how the balance of linear and circular ANRIL RNA, determined by the Chr9p21 genotype, regulates molecular pathways and cellular functions involved in atherogenesis. We end with an outlook on how manipulating circular ANRIL abundance may be exploited for therapeutic purposes.

Long Non-Coding RNAs in Obesity-Induced Cancer

Many mechanisms of obesity-induced cancers have been proposed. However, it remains unclear whether or not long non-coding RNAs (lncRNAs) play any role in obesity-induced cancers. In this article, we briefly discuss the generally accepted hypotheses explaining the mechanisms of obesity-induced cancers, summarize the latest evidence for the expression of a number of well-known cancer-associated lncRNAs in obese subjects, and propose the potential contribution of lncRNAs to obesity-induced cancers. We hope this review can serve as an inspiration to scientists to further explore the regulatory roles of lncRNAs in the development of obesity-induced cancers. Those findings will be fundamental in the development of effective therapeutics or interventions to combat this life-threatening adverse effect of obesity.

Identification of LINC01279 as a cell cycle‑associated long non‑coding RNA in endometriosis with GBA analysis

Endometriosis affects 6‑10% of women of reproductive age. Though a significant amount of research has explored the pathogenesis of endometriosis, little is clear. Elucidating the mechanisms is urgently required for improving the therapeutic efficiency of endometriosis treatment. Long non‑coding RNAs (lncRNAs) have recently acquired extensive attention as regulatory components in a variety of biological processes and diseases. However, the functions of many lncRNAs in endometriosis are poorly understood. Therefore, the exploration of the dysregulated genes in endometriosis, particularly lncRNAs, is of importance. In the present study, datasets for endometriosis, including GSE7305, GSE7846, GSE29981 and E‑MTAB‑694, were downloaded from Gene Expression Omnibus and ArrayExpress. Then, the limma and Affy packages were used to analyze the CEL file. The RankProd method was used to conduct meta‑analysis. Long intergenic non‑protein coding RNA 1279 (LINC01279) was significantly upregulated in the three datasets, and was the most upregulated lncRNA as determined by the RankProd method. Gene set enrichment and Gene Ontology analyses were conducted, which revealed that LINC01279 is likely to function as a cell cycle mediator in endometriosis. Finally, it was identified that LINC01279 is strongly associated with certain previously identified key factors in the development of endometriosis, including cyclin‑dependent kinase 14 and C‑X‑C motif chemokine ligand 12. Thus, it was demonstrated that LINC01279 may be associated with the pathogenesis of endometriosis. This may potentially represent a target in the therapy of endometriosis.

Genome-wide haplotype association study identifies risk genes for non-small cell lung cancer

Lung cancer is the leading cause of cancer-related death worldwide. Most lung cancer is non-small cell lung cancer (NSCLC), in which malignant cells form in the lung epithelium. Mutations in multiple genes and environmental factors both contribute to NSCLC, and although some NSCLC susceptibility genes have been characterized, the pathogenesis of this disease remains unclear. To identify genes conferring NSCLC risk and determine their associated pathological mechanism, we combined genome-wide haplotype associated analysis with gene prioritization using 224,677 SNPs in 37 NSCLC cell lines and 116 unrelated European individuals. Five candidate genes were identified: ESR1, TGFBR1, INSR, CDH3, and MAP3K5. All of these have previously been implicated in NSCLC, with the exception of CDH3, which can therefore be considered a novel indicator of NSCLC risk. Functional annotation confirmed the relationship between these five genes and NSCLC. Our findings are indicative of the underlying pathological mechanisms of NSCLC and provide information to support future directions in diagnosing and treating NSCLC.

ANRIL: A lncRNA at the CDKN2A/B Locus With Roles in Cancer and Metabolic Disease

The CDKN2A/B genomic locus is associated with risk of human cancers and metabolic disease. Although the locus contains several important protein-coding genes, studies suggest disease roles for a lesser-known antisense lncRNA encoded at this locus, called ANRIL. ANRIL is a complex gene containing at least 21 exons in simians, with many reported linear and circular isoforms. Like other genes, abundance of ANRIL is regulated by epigenetics, classic transcription regulation, splicing, and post-transcriptional influences such as RNA stability and microRNAs. Known molecular functions of ANRIL include in cis and in trans gene regulation through chromatin modification complexes, and influence over microRNA signaling networks. Polymorphisms at the ANRIL gene are linked to risk for many different cancers, as well as risk of atherosclerotic cardiovascular disease, bone mass, obesity and type 2 diabetes. A broad array of variable reported impacts of polymorphisms on ANRIL abundance, splicing and function suggests that ANRIL has cell-type and context-dependent regulation and actions. In cancer cells, ANRIL gain of function increases proliferation, metastasis, cell survival and epithelial-mesenchymal transformation, whereas ANRIL loss of function decreases tumor size and growth, invasion and metastasis, and increases apoptosis and senescence. In metabolic disease, polymorphisms at the ANRIL gene are linked to risk of type 2 diabetes, coronary artery disease, coronary artery calcium score, myocardial infarction, and stroke. Intriguingly, with the exception of one polymorphism in exon 2 of ANRIL, the single nucleotide polymorphisms (SNPs) associated with atherosclerosis and diabetes are non-overlapping. Evidence suggests that ANRIL gain of function increases atherosclerosis; in diabetes, a risk-SNP reduced the pancreatic beta cell proliferation index. Studies are limited by the uncertain relevance of rodent models to ANRIL studies, since most ANRIL exons do not exist in mouse. Diverse cell-type-dependent results suggest it is necessary to perform studies in the relevant primary human tissue for each disease. Much remains to be learned about the biology of ANRIL in human health and disease; this research area may lead to insight into disease mechanisms and therapeutic approaches.

A novel indel within goat casein alpha S1 gene is significantly associated with litter size

The ruminant casein gene family (CSNs, link as CSN1S1-CSN2-CSN1S2-CSN3) is characterized by diverse variations and has been extensively studied for ruminant milk traits, however, studies on insertion/deletion (indel) mutations within this gene family and their effects on prolificacy are extremely limited. This study aimed to detect possible novel indels within CSNs in an indigenous Chinese goat breed-Shaanbei White Cashmere goat (SBWC, n = 3047) and four other Chinese goat breeds (n = 1136) with varied litter size rates (105%-283%) and different estrus types (seasonal vs. perennial), as well as exploring the association between these potential indels and litter size. Only one novel 11-bp indel within the CSN1S1 gene was found. The association analyses uncovered that this novel indel was related to the first-birth litter size of SBWC population (n = 2690) (P < 1.0 E-8). Individuals with the II genotype (n = 676) had the best litter size when compared with those ID genotype (n = 1098) and DD genotype (n = 916) individuals. Animals with the II genotype were found to have higher relative expression level of CSN1S1 gene in the ovary (P < 0.01). Besides, Chi-square tests for different litter size and estrous cycle breeds showed that perennial-estrus breeds and multi-kids breeds had higher "I" allelic frequencies and "II" genotypic frequencies. These findings suggest the 11-bp indel within the CSN1S1 gene is significantly associated with reproduction traits and can be an effective molecular marker for litter size of goat breeding.

Insertion/Deletion Within the KDM6A Gene Is Significantly Associated With Litter Size in Goat

A previous whole-genome association analysis identified lysine demethylase 6A (KDM6A), which encodes a type of histone demethylase, as a candidate gene associated to goat fecundity. KDM6A gene knockout mouse disrupts gametophyte development, suggesting that it has a critical role in reproduction. In this study, goat KDM6A mRNA expression profiles were determined, insertion/deletion (indel) variants in the gene identified, indel variants effect on KDM6A gene expression assessed, and their association with first-born litter size analyzed in 2326 healthy female Shaanbei white cashmere goats. KDM6A mRNA was expressed in all tissues tested (heart, liver, spleen, lung, kidney, muscle, brain, skin and testis); the expression levels in testes at different developmental stages [1-week-old (wk), 2, 3 wk, 1-month-old (mo), 1.5 and 2 mo] indicated a potential association with the mitosis-to-meiosis transition, implying that KDM6A may have an essential role in goat fertility. Meanwhile, two novel intronic indels of 16 bp and 5 bp were identified. Statistical analysis revealed that only the 16 bp indel was associated with first-born litter size (P < 0.01), and the average first-born litter size of individuals with an insertion/insertion genotype higher than that of those with the deletion/deletion genotype (P < 0.05). There was also a significant difference in genotype distributions of the 16 bp indel between mothers of single-lamb and multi-lamb litters in the studied goat population (P = 0.001). Consistently, the 16 bp indel also had a significant effect on KDM6A gene expression. Additionally, there was no significant linkage disequilibrium (LD) between these two indel loci, consistent with the association analysis results. Together, these findings suggest that the 16 bp indel in KDM6A may be useful for marker-assisted selection (MAS) of goats.

ANRIL regulates production of extracellular matrix proteins and vasoactive factors in diabetic complications

noncoding RNAs (lncRNAs) have gained widespread interest due to their prevailing presence in various diseases. lncRNA ANRIL (a. k. a. CDKN2B-AS1) is located on human chromosome 9 (p21.3) and transcribed in opposite direction to the INK4b-ARF-INK4a gene cluster. It has been identified as a highly susceptible region for diseases such as coronary artery diseases and type 2 diabetes. Here, we explored its regulatory role in diabetic nephropathy (DN) and diabetic cardiomyopathy (DCM) in association with epigenetic modifiers p300 and polycomb repressive complex 2 (PRC2) complex. We used an ANRIL-knockout (ANRILKO) mouse model for this study. The wild-type and ANRILKO animals with or without streptozotocin-induced diabetes were monitored for 2 min. At the end of the time point, urine and tissues were collected. The tissues were measured for fibronectin (FN), type IV collagen (Col1α4), and VEGF mRNA and protein expressions. Renal function was determined by the measurement of 24-h urine volume and albumin/creatinine ratio at euthanasia. Renal and cardiac structures were investigated using periodic acid-Schiff stain and/or immunohistochemical analysis. Elevated expressions of extracellular matrix (ECM) proteins were prevented in ANRILKO diabetic animals. Furthermore, ANRILKO had a protective effect on diabetic mouse kidneys, as evidenced by lowering of urine volume and urine albumin levels in comparison with the wild-type diabetic animals. These alterations regulated by ANRIL may be mediated by p300 and enhancer of zeste 2 (EZH2) of the PRC2 complex. Our study concludes that ANRIL regulates functional and structural alterations in the kidneys and hearts in diabetes through controlling the expressions of ECM proteins and VEGF.

Molecular roles and function of circular RNAs in eukaryotic cells

Protein-coding and noncoding genes in eukaryotes are typically expressed as linear messenger RNAs, with exons arranged colinearly to their genomic order. Recent advances in sequencing and in mapping RNA reads to reference genomes have revealed that thousands of genes express also covalently closed circular RNAs. Many of these circRNAs are stable and contain exons, but are not translated into proteins. Here, we review the emerging understanding that both, circRNAs produced by co- and posttranscriptional head-to-tail "backsplicing" of a downstream splice donor to a more upstream splice acceptor, as well as circRNAs generated from intronic lariats during colinear splicing, may exhibit physiologically relevant regulatory functions in eukaryotes. We describe how circRNAs impact gene expression of their host gene locus by affecting transcriptional initiation and elongation or splicing, and how they partake in controlling the function of other molecules, for example by interacting with microRNAs and proteins. We conclude with an outlook how circRNA dysregulation affects disease, and how the stability of circRNAs might be exploited in biomedical applications.

Intrinsic MYH7 expression regulation contributes to tissue level allelic imbalance in hypertrophic cardiomyopathy

HCM, the most common inherited cardiac disease, is mainly caused by mutations in sarcomeric genes. More than a third of the patients are heterozygous for mutations in the MYH7 gene encoding for the β-myosin heavy chain. In HCM-patients, expression of the mutant and the wildtype allele can be unequal, thus leading to fractions of mutant and wildtype mRNA and protein which deviate from 1:1. This so-called allelic imbalance was detected in whole tissue samples but also in individual cells. There is evidence that the severity of HCM not only depends on the functional effect of the mutation itself, but also on the fraction of mutant protein in the myocardial tissue. Allelic imbalance has been shown to occur in a broad range of genes. Therefore, we aimed to examine whether the MYH7-alleles are intrinsically expressed imbalanced or whether the allelic imbalance is solely associated with the disease. We compared the expression of MYH7-alleles in non-HCM donors and in HCM-patients with different MYH7-missense mutations. In the HCM-patients, we identified imbalanced as well as equal expression of both alleles. Also at the protein level, allelic imbalance was determined. Most interestingly, we also discovered allelic imbalance and balance in non-HCM donors. Our findings therefore strongly indicate that apart from mutation-specific mechanisms, also non-HCM associated allelic-mRNA expression regulation may account for the allelic imbalance of the MYH7 gene in HCM-patients. Since the relative amount of mutant mRNA and protein or the extent of allelic imbalance has been associated with the severity of HCM, individual analysis of the MYH7-allelic expression may provide valuable information for the prognosis of each patient.

Long non-coding RNA LINC00261 inhibits cell growth and migration in endometriosis

AIM

A previous study reported that LINC00261 is significantly downregulated in human ectopic endometrial tissues. The present study aimed to explore whether LINC00261 is functional in endometriosis cell proliferation, apoptosis and migration.

METHODS

By transfecting human endometriosis cell line CRL-7566 with plasmids containing LINC00261, we successfully established the cell CRL-7566/LINC00261 with a high LINC00261 expression level. Cell-counting kit-8 and colony formation assays were conducted to evaluate the effect of LINC00261 on cell proliferation, and flow cytometry analysis and transwell migration assay were conducted to evaluate its effect on cell apoptosis and cell migration, respectively.

RESULTS

Cell-counting kit-8 and colony formation assays both indicated that LINC00261 could inhibit cell proliferation in CRL-7566. Flow cytometry analysis confirmed that LINC00261 mediated inhibition of cell proliferation, which might be a consequence of inducting apoptosis. Furthermore, transwell migration assay indicated that LINC00261 could inhibit cell migration in endometriosis.

CONCLUSION

LncRNA LINC00261 is capable of inhibiting cell growth and migration in endometriosis.

Meta-analysis identifies five novel loci associated with endometriosis highlighting key genes involved in hormone metabolism

Endometriosis is a heritable hormone-dependent gynecological disorder, associated with severe pelvic pain and reduced fertility; however, its molecular mechanisms remain largely unknown. Here we perform a meta-analysis of 11 genome-wide association case-control data sets, totalling 17,045 endometriosis cases and 191,596 controls. In addition to replicating previously reported loci, we identify five novel loci significantly associated with endometriosis risk (P<5 × 10⁻⁸), implicating genes involved in sex steroid hormone pathways (FN1, CCDC170, ESR1, SYNE1 and FSHB). Conditional analysis identified five secondary association signals, including two at the ESR1 locus, resulting in 19 independent single nucleotide polymorphisms (SNPs) robustly associated with endometriosis, which together explain up to 5.19% of variance in endometriosis. These results highlight novel variants in or near specific genes with important roles in sex steroid hormone signalling and function, and offer unique opportunities for more targeted functional research efforts.

Endometriosis is a major cause of infertility. Molecular mechanisms underlying the disease involve genetic and environmental risk factors. In a meta-analysis of eleven GWA studies, Sapkota and colleagues identify five novel risk loci, implicating steroid sex hormone pathways in the pathogenesis.

Role of the LF-SINE-Derived Distal ISL1 Enhancer in Patients with Classic Bladder Exstrophy

A genome-wide association study and meta-analysis identified ISL1 as the first genome-wide significant susceptibility gene for classic bladder exstrophy (CBE). A short interspersed repetitive element (SINE), first detected in lobe-finned fishes (LF-SINE), was shown to drive Isl1 expression in embryonic mouse genital eminence. Hence, we assumed this enhancer a conclusive target for mutations associated with CBE formation and analyzed a cohort of 200 CBE patients. Although we identified two enhancer variants in five CBE patients, their clinical significance seems unlikely, implying that sequence variants in the ISL1 LF-SINE enhancer are not frequently associated with CBE.

Genetic Risk Factors for Endometriosis

Advances in genetics and genomics are driving progress in understanding genetic risk factors for endometriosis. Genome-wide association scans (GWAS) in endometriosis have identified 11 genomic regions associated with increased risk of disease. Many of the regions contain interesting candidate genes, but the risk alleles may not always act through the obvious candidates. Functional evidence to identify the causal gene(s) will require multiple steps including better mapping precision, genetic studies on gene expression and epigenetic marks, chromatin looping and functional studies. Evidence from gene expression studies in endometrium and chromatin looping experiments implicate CDC42 on chromosome 1, CDKN2B-AS1 on chromosome 9 and VEZT on chromosome 12 as likely causal genes in these regions. Confirming the causal gene(s) in these and other regions will identify the important pathways increasing risk for endometriosis and identify novel targets for interventions to improve diagnosis and treatment.