Urinary miRNAs: Technical Updates.

Due to its non-invasive nature and easy accessibility, urine serves as a convenient bi-ological fluid for research purposes. Furthermore, urine samples are uncomplicated to preserve and relatively inexpensive. MicroRNA (miRNAs), small molecules that regulate gene expression post-transcriptionally, play vital roles in numerous cellular processes, including apoptosis, cell differentiation, development, and proliferation. Their dysregulated expression in urine has been proposed as a potential bi-omarker for various human diseases, including bladder cancer. To draw reliable conclusions about the roles of urinary miRNAs in human diseases, it is essential to have dependable and reproducible methods for miRNA extraction and profiling. In this review, we address the technical challenges associated with studying urinary miRNAs and provide an update on the current technologies used for urinary miRNA isolation, quality control assessment, and miRNA profiling, highlighting both their advantages and limitations.

Understanding the Genetics of Early-Onset Obesity in a Cohort of Children From Qatar.

CONTEXT: Monogenic obesity is a rare form of obesity due to pathogenic variants in genes implicated in the leptin-melanocortin signaling pathway and accounts for around 5% of severe early-onset obesity. Mutations in the genes encoding the MC4R, leptin, and leptin receptor are commonly reported in various populations to cause monogenic obesity. Determining the genetic cause has important clinical benefits as novel therapeutic interventions are now available for some forms of monogenic obesity. OBJECTIVE: To unravel the genetic causes of early-onset obesity in the population of Qatar. METHODS: In total, 243 patients with early-onset obesity (above the 95% percentile) and age of onset below 10 years were screened for monogenic obesity variants using a targeted gene panel, consisting of 52 obesity-related genes. RESULTS: Thirty rare variants potentially associated with obesity were identified in 36 of 243 (14.8%) probands in 15 candidate genes (LEP, LEPR, POMC, MC3R, MC4R, MRAP2, SH2B1, BDNF, NTRK2, DYRK1B, SIM1, GNAS, ADCY3, RAI1, and BBS2). Twenty-three of the variants identified were novel to this study and the rest, 7 variants, were previously reported in literature. Variants in MC4R were the most common cause of obesity in our cohort (19%) and the c.485C>T p.T162I variant was the most frequent MC4R variant seen in 5 patients. CONCLUSION: We identified likely pathogenic/pathogenic variants that seem to explain the phenotype of around 14.8% of our cases. Variants in the MC4R gene are the commonest cause of early-onset obesity in our population. Our study represents the largest monogenic obesity cohort in the Middle East and revealed novel obesity variants in this understudied population. Functional studies will be required to elucidate the molecular mechanism of their pathogenicity.

Discovery of new therapeutic targets in ovarian cancer through identifying significantly non-mutated genes.

BACKGROUND: Mutated and non-mutated genes interact to drive cancer growth and metastasis. While research has focused on understanding the impact of mutated genes on cancer biology, understanding non-mutated genes that are essential to tumor development could lead to new therapeutic strategies. The recent advent of high-throughput whole genome sequencing being applied to many different samples has made it possible to calculate if genes are significantly non-mutated in a specific cancer patient cohort. METHODS: We carried out random mutagenesis simulations of the human genome approximating the regions sequenced in the publicly available Cancer Growth Atlas Project for ovarian cancer (TCGA-OV). Simulated mutations were compared to the observed mutations in the TCGA-OV cohort and genes with the largest deviations from simulation were identified. Pathway analysis was performed on the non-mutated genes to better understand their biological function. We then compared gene expression, methylation and copy number distributions of non-mutated and mutated genes in cell lines and patient data from the TCGA-OV project. To directly test if non-mutated genes can affect cell proliferation, we carried out proof-of-concept RNAi silencing experiments of a panel of nine selected non-mutated genes in three ovarian cancer cell lines and one primary ovarian epithelial cell line. RESULTS: We identified a set of genes that were mutated less than expected (non-mutated genes) and mutated more than expected (mutated genes). Pathway analysis revealed that non-mutated genes interact in cancer associated pathways. We found that non-mutated genes are expressed significantly more than mutated genes while also having lower methylation and higher copy number states indicating that they could be functionally important. RNAi silencing of the panel of non-mutated genes resulted in a greater significant reduction of cell viability in the cancer cell lines than in the non-cancer cell line. Finally, as a test case, silencing ANKLE2, a significantly non-mutated gene, affected the morphology, reduced migration, and increased the chemotherapeutic response of SKOV3 cells. CONCLUSION: We show that we can identify significantly non-mutated genes in a large ovarian cancer cohort that are well-expressed in patient and cell line data and whose RNAi-induced silencing reduces viability in three ovarian cancer cell lines. Targeting non-mutated genes that are important for tumor growth and metastasis is a promising approach to expand cancer therapeutic options.

Cell Type-Specific TGF-ß Mediated EMT in 3D and 2D Models and Its Reversal by TGF-ß Receptor Kinase Inhibitor in Ovarian Cancer Cell Lines.

Transcriptome profiling of 3D models compared to 2D models in various cancer cell lines shows differential expression of TGF-ß-mediated and cell adhesion pathways. Presence of TGF-ß in these cell lines shows an increased invasion potential which is specific to cell type. In the present study, we identified exogenous addition of TGF-ß can induce Epithelial to Mesenchymal Transition (EMT) in a few cancer cell lines. RNA sequencing and real time PCR were carried out in different ovarian cancer cell lines to identify molecular profiling and metabolic profiling. Since EMT induction by TGF-ß is cell-type specific, we decided to select two promising ovarian cancer cell lines as model systems to study EMT. TGF-ß modulation in EMT and cancer invasion were successfully depicted in both 2D and 3D models of SKOV3 and CAOV3 cell lines. Functional evaluation in 3D and 2D models demonstrates that the addition of the exogenous TGF-ß can induce EMT and invasion in cancer cells by turning them into aggressive phenotypes. TGF-ß receptor kinase I inhibitor (LY364947) can revert the TGF-ß effect in these cells. In a nutshell, TGF-ß can induce EMT and migration, increase aggressiveness, increase cell survival, alter cell characteristics, remodel the Extracellular Matrix (ECM) and increase cell metabolism favorable for tumor invasion and metastasis. We concluded that transcriptomic and phenotypic effect of TGF-ß and its inhibitor is cell-type specific and not cancer specific.

Silencing of ANKRD12 circRNA induces molecular and functional changes associated with invasive phenotypes.

BACKGROUND: Circular RNAs (circRNAs) that form through non-canonical backsplicing events of pre-mRNA transcripts are evolutionarily conserved and abundantly expressed across species. However, the functional relevance of circRNAs remains a topic of debate. METHODS: We identified one of the highly expressed circRNA (circANKRD12) in cancer cell lines and characterized it validated it by Sanger sequencing, Real-Time PCR. siRNA mediated silencing of the circular junction of circANKRD12 was followed by RNA Seq analysis of circANKRD12 silenced cells and control cells to identify the differentially regulated genes. A series of cell biology and molecular biology techniques (MTS assay, Migration analysis, 3D organotypic models, Real-Time PCR, Cell cycle analysis, Western blot analysis, and Seahorse Oxygen Consumption Rate analysis) were performed to elucidate the function, and underlying mechanisms involved in circANKRD12 silenced breast and ovarian cancer cells. RESULTS: In this study, we identified and characterized a circular RNA derived from Exon 2 and Exon 8 of the ANKRD12 gene, termed here as circANKRD12. We show that this circRNA is abundantly expressed in breast and ovarian cancers. The circANKRD12 is RNase R resistant and predominantly localized in the cytoplasm in contrast to its source mRNA. We confirmed the expression of this circRNA across a variety of cancer cell lines and provided evidence for its functional relevance through downstream regulation of several tumor invasion genes. Silencing of circANKRD12 induces a strong phenotypic change by significantly regulating cell cycle, increasing invasion and migration and altering the metabolism in cancer cells. These results reveal the functional significance of circANKRD12 and provide evidence of a regulatory role for this circRNA in cancer progression. CONCLUSIONS: Our study demonstrates the functional relevance of circANKRD12 in various cancer cell types and, based on its expression pattern, has the potential to become a new clinical biomarker.

Author Correction: Genus-wide sequencing support a two-locus model for sex-determination in Phoenix.

In the original version of this Article, the affiliation of the first author, Maria F. Torres, 'Department of Biological Sciences, University of Cincinnati, Cincinnati, 45221, OH, USA' was incorrectly assigned as a present address and should have been listed as a full affiliation. This error has been corrected in both the PDF and HTML versions of the Article.

Genus-wide sequencing supports a two-locus model for sex-determination in Phoenix.

The date palm tree is a commercially important member of the genus Phoenix whose 14 species are dioecious with separate male and female individuals. To identify sex determining genes we sequenced the genomes of 15 female and 13 male Phoenix trees representing all 14 species. We identified male-specific sequences and extended them using phased single-molecule sequencing or BAC clones. We observed that only four genes contained sequences conserved in all analyzed Phoenix males. Most of these sequences showed similarity to a single genomic locus in the closely related monoecious oil palm. CYP703 and GPAT3, two single copy genes present in males and critical for male flower development in other monocots, were absent in females. A LOG-like gene appears translocated into the Y-linked region and is suggested to play a role in suppressing female flowers. Our data are consistent with a two-mutation model for the evolution of dioecy in Phoenix.

Mutation in noncoding RNA RNU12 causes early onset cerebellar ataxia.

OBJECTIVE: Exome sequences account for only 2% of the genome and may overlook mutations causing disease. To obtain a more complete view, whole genome sequencing (WGS) was analyzed in a large consanguineous family in which members displayed autosomal recessively inherited cerebellar ataxia manifesting before 2 years of age. METHODS: WGS from blood-derived genomic DNA was used for homozygosity mapping and a rare variant search. RNA from isolated blood leukocytes was used for quantitative polymerase chain reaction (PCR), RNA sequencing, and comparison of the transcriptomes of affected and unaffected family members. RESULTS: WGS revealed a point mutation in noncoding RNA RNU12 that was associated with early onset cerebellar ataxia. The U12-dependent minor spliceosome edits 879 known transcripts. Reverse transcriptase PCR demonstrated minor intron retention in all of 9 randomly selected RNAs from this group, and RNAseq showed splicing disruption specific to all U12-type introns detected in blood monocytes from affected individuals. Moreover, 144 minor intron-containing RNAs were differentially expressed, including transcripts for 3 genes previously associated with cerebellar neurodegeneration. INTERPRETATION: Interference with particular spliceosome components, including small nuclear RNAs, cause reproducible uniquely distributed phenotypic and transcript-specific effects, making this an important category of disease-associated mutation. Our approach to differential expression analysis of minor intron-containing genes is applicable to other diseases involving altered transcriptome processing. ANN NEUROL 2017;81:68-78.

Genetic Epidemiology of Glucose-6-Phosphate Dehydrogenase Deficiency in the Arab World.

A systematic search was implemented using four literature databases (PubMed, Embase, Science Direct and Web of Science) to capture all the causative mutations of Glucose-6-phosphate dehydrogenase (G6PD) deficiency (G6PDD) in the 22 Arab countries. Our search yielded 43 studies that captured 33 mutations (23 missense, one silent, two deletions, and seven intronic mutations), in 3,430 Arab patients with G6PDD. The 23 missense mutations were then subjected to phenotypic classification using in silico prediction tools, which were compared to the WHO pathogenicity scale as a reference. These in silico tools were tested for their predicting efficiency using rigorous statistical analyses. Of the 23 missense mutations, p.S188F, p.I48T, p.N126D, and p.V68M, were identified as the most common mutations among Arab populations, but were not unique to the Arab world, interestingly, our search strategy found four other mutations (p.N135T, p.S179N, p.R246L, and p.Q307P) that are unique to Arabs. These mutations were exposed to structural analysis and molecular dynamics simulation analysis (MDSA), which predicting these mutant forms as potentially affect the enzyme function. The combination of the MDSA, structural analysis, and in silico predictions and statistical tools we used will provide a platform for future prediction accuracy for the pathogenicity of genetic mutations.

Correction: Type 2 Diabetes Risk Allele Loci in the Qatari Population.

[This corrects the article DOI: 10.1371/journal.pone.0156834.].

Type 2 Diabetes Risk Allele Loci in the Qatari Population.

BACKGROUND: The prevalence of type 2 diabetes (T2D) is increasing in the Middle East. However, the genetic risk factors for T2D in the Middle Eastern populations are not known, as the majority of studies of genetic risk for T2D are in Europeans and Asians. METHODS: All subjects were ≥3 generation Qataris. Cases with T2D (n = 1,124) and controls (n = 590) were randomly recruited and assigned to the 3 known Qatari genetic subpopulations [Bedouin (Q1), Persian/South Asian (Q2) and African (Q3)]. Subjects underwent genotyping for 37 single nucleotide polymorphisms (SNPs) in 29 genes known to be associated with T2D in Europeans and/or Asian populations, and an additional 27 tag SNPs related to these susceptibility loci. Pre-study power analysis suggested that with the known incidence of T2D in adult Qataris (22%), the study population size would be sufficient to detect significant differences if the SNPs were risk factors among Qataris, assuming that the odds ratio (OR) for T2D SNPs in Qatari's is greater than or equal to the SNP with highest known OR in other populations. RESULTS: Haplotype analysis demonstrated that Qatari haplotypes in the region of known T2D risk alleles in Q1 and Q2 genetic subpopulations were similar to European haplotypes. After Benjamini-Hochberg adjustment for multiple testing, only two SNPs (rs7903146 and rs4506565), both associated with transcription factor 7-like 2 (TCF7L2), achieved statistical significance in the whole study population. When T2D subjects and control subjects were assigned to the known 3 Qatari subpopulations, and analyzed individually and with the Q1 and Q2 genetic subpopulations combined, one of these SNPs (rs4506565) was also significant in the admixed group. No other SNPs associated with T2D in all Qataris or individual genetic subpopulations. CONCLUSIONS: With the caveats of the power analysis, the European/Asian T2D SNPs do not contribute significantly to the high prevalence of T2D in the Qatari population, suggesting that the genetic risks for T2D are likely different in Qataris compared to Europeans and Asians.

Altered expression pattern of circular RNAs in primary and metastatic sites of epithelial ovarian carcinoma.

Recently, a class of endogenous species of RNA called circular RNA (circRNA) has been shown to regulate gene expression in mammals and their role in cellular function is just beginning to be understood. To investigate the role of circRNAs in ovarian cancer, we performed paired-end RNA sequencing of primary sites, peritoneal and lymph node metastases from three patients with stage IIIC ovarian cancer. We developed an in-house computational pipeline to identify and characterize the circRNA expression from paired-end RNA-Seq libraries. This pipeline revealed thousands of circular isoforms in Epithelial Ovarian Carcinoma (EOC). These circRNAs are enriched for potentially effective miRNA seed matches. A significantly larger number of circRNAs are differentially expressed between tumor sites than mRNAs. Circular and linear expression exhibits an inverse trend for many cancer related pathways and signaling pathways like NFkB, PI3k/AKT and TGF-ß typically activated for mRNA in metastases are inhibited for circRNA expression. Further, circRNAs show a more robust expression pattern across patients than mRNA forms indicating their suitability as biomarkers in highly heterogeneous cancer transcriptomes. The consistency of circular RNA expression may offer new candidates for cancer treatment and prognosis.

Preferential Allele Expression Analysis Identifies Shared Germline and Somatic Driver Genes in Advanced Ovarian Cancer.

Identifying genes where a variant allele is preferentially expressed in tumors could lead to a better understanding of cancer biology and optimization of targeted therapy. However, tumor sample heterogeneity complicates standard approaches for detecting preferential allele expression. We therefore developed a novel approach combining genome and transcriptome sequencing data from the same sample that corrects for sample heterogeneity and identifies significant preferentially expressed alleles. We applied this analysis to epithelial ovarian cancer samples consisting of matched primary ovary and peritoneum and lymph node metastasis. We find that preferentially expressed variant alleles include germline and somatic variants, are shared at a relatively high frequency between patients, and are in gene networks known to be involved in cancer processes. Analysis at a patient level identifies patient-specific preferentially expressed alleles in genes that are targets for known drugs. Analysis at a site level identifies patterns of site specific preferential allele expression with similar pathways being impacted in the primary and metastasis sites. We conclude that genes with preferentially expressed variant alleles can act as cancer drivers and that targeting those genes could lead to new therapeutic strategies.