miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades

microRNAs (miRNAs) are a large class of small non-coding RNAs which post-transcriptionally regulate the expression of a large fraction of all animal genes and are important in a wide range of biological processes. Recent advances in high-throughput sequencing allow miRNA detection at unprecedented sensitivity, but the computational task of accurately identifying the miRNAs in the background of sequenced RNAs remains challenging. For this purpose, we have designed miRDeep2, a substantially improved algorithm which identifies canonical and non-canonical miRNAs such as those derived from transposable elements and informs on high-confidence candidates that are detected in multiple independent samples. Analyzing data from seven animal species representing the major animal clades, miRDeep2 identified miRNAs with an accuracy of 98.6–99.9% and reported hundreds of novel miRNAs. To test the accuracy of miRDeep2, we knocked down the miRNA biogenesis pathway in a human cell line and sequenced small RNAs before and after. The vast majority of the >100 novel miRNAs expressed in this cell line were indeed specifically downregulated, validating most miRDeep2 predictions. Last, a new miRNA expression profiling routine, low time and memory usage and user-friendly interactive graphic output can make miRDeep2 useful to a wide range of researchers

MicroRNAs targeting oncogenes are down-regulated in pancreatic malignant transformation from benign tumors

BACKGROUND MicroRNA (miRNA) expression profiles have been described in pancreatic ductal adenocarcinoma (PDAC), but these have not been compared with pre-malignant pancreatic tumors. We wished to compare the miRNA expression signatures in pancreatic benign cystic tumors (BCT) of low and high malignant potential with PDAC, in order to identify miRNAs deregulated during PDAC development. The mechanistic consequences of miRNA dysregulation were further evaluated. METHODS Tissue samples were obtained at a tertiary pancreatic unit from individuals with BCT and PDAC. MiRNA profiling was performed using a custom microarray and results were validated using RT-qPCR prior to evaluation of miRNA targets. RESULTS Widespread miRNA down-regulation was observed in PDAC compared to low malignant potential BCT. We show that amongst those miRNAs down-regulated, miR-16, miR-126 and let-7d regulate known PDAC oncogenes (targeting BCL2, CRK and KRAS respectively). Notably, miR-126 also directly targets the KRAS transcript at a "seedless" binding site within its 3'UTR. In clinical specimens, miR-126 was strongly down-regulated in PDAC tissues, with an associated elevation in KRAS and CRK proteins. Furthermore, miR-21, a known oncogenic miRNA in pancreatic and other cancers, was not elevated in PDAC compared to serous microcystic adenoma (SMCA), but in both groups it was up-regulated compared to normal pancreas, implicating early up-regulation during malignant change. CONCLUSIONS Expression profiling revealed 21 miRNAs down-regulated in PDAC compared to SMCA, the most benign lesion that rarely progresses to invasive carcinoma. It appears that miR-21 up-regulation is an early event in the transformation from normal pancreatic tissue. MiRNA expression has the potential to distinguish PDAC from normal pancreas and BCT. Mechanistically the down-regulation of miR-16, miR-126 and let-7d promotes PDAC transformation by post-transcriptional up-regulation of crucial PDAC oncogenes. We show that miR-126 is able to directly target KRAS; re-expression has the potential as a therapeutic strategy against PDAC and other KRAS-driven cancers

Sequence variation between 462 human individuals fine-tunes functional sites of RNA processing

A. Palotie on työryhmän GEUVADIS Consortium jäsen.Recent advances in the cost-efficiency of sequencing technologies enabled the combined DNA-and RNA-sequencing of human individuals at the population-scale, making genome-wide investigations of the inter-individual genetic impact on gene expression viable. Employing mRNA-sequencing data from the Geuvadis Project and genome sequencing data from the 1000 Genomes Project we show that the computational analysis of DNA sequences around splice sites and poly-A signals is able to explain several observations in the phenotype data. In contrast to widespread assessments of statistically significant associations between DNA polymorphisms and quantitative traits, we developed a computational tool to pinpoint the molecular mechanisms by which genetic markers drive variation in RNA-processing, cataloguing and classifying alleles that change the affinity of core RNA elements to their recognizing factors. The in silico models we employ further suggest RNA editing can moonlight as a splicing-modulator, albeit less frequently than genomic sequence diversity. Beyond existing annotations, we demonstrate that the ultra-high resolution of RNA-Seq combined from 462 individuals also provides evidence for thousands of bona fide novel elements of RNA processing-alternative splice sites, introns, and cleavage sites-which are often rare and lowly expressed but in other characteristics similar to their annotated counterparts.Peer reviewe

Global human footprint on the linkage between biodiversity and ecosystem functioning in reef fishes

Copyright: © 2011 Mora et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.Difficulties in scaling up theoretical and experimental results have raised controversy over the consequences of biodiversity loss for the functioning of natural ecosystems. Using a global survey of reef fish assemblages, we show that in contrast to previous theoretical and experimental studies, ecosystem functioning (as measured by standing biomass) scales in a non-saturating manner with biodiversity (as measured by species and functional richness) in this ecosystem. Our field study also shows a significant and negative interaction between human population density and biodiversity on ecosystem functioning (i.e., for the same human density there were larger reductions in standing biomass at more diverse reefs). Human effects were found to be related to fishing, coastal development, and land use stressors, and currently affect over 75% of the world's coral reefs. Our results indicate that the consequences of biodiversity loss in coral reefs have been considerably underestimated based on existing knowledge and that reef fish assemblages, particularly the most diverse, are greatly vulnerable to the expansion and intensity of anthropogenic stressors in coastal areas

The Influence of Number and Timing of Pregnancies on Breast Cancer Risk for Women With BRCA1 or BRCA2 Mutations

International audienceBACKGROUND:Full-term pregnancy (FTP) is associated with a reduced breast cancer (BC) risk over time, but women are at increased BC risk in the immediate years following an FTP. No large prospective studies, however, have examined whether the number and timing of pregnancies are associated with BC risk for BRCA1 and BRCA2 mutation carriers.METHODS:Using weighted and time-varying Cox proportional hazards models, we investigated whether reproductive events are associated with BC risk for mutation carriers using a retrospective cohort (5707 BRCA1 and 3525 BRCA2 mutation carriers) and a prospective cohort (2276 BRCA1 and 1610 BRCA2 mutation carriers), separately for each cohort and the combined prospective and retrospective cohort.RESULTS:For BRCA1 mutation carriers, there was no overall association with parity compared with nulliparity (combined hazard ratio [HRc] = 0.99, 95% confidence interval [CI] = 0.83 to 1.18). Relative to being uniparous, an increased number of FTPs was associated with decreased BC risk (HRc = 0.79, 95% CI = 0.69 to 0.91; HRc = 0.70, 95% CI = 0.59 to 0.82; HRc = 0.50, 95% CI = 0.40 to 0.63, for 2, 3, and ≥4 FTPs, respectively, P trend < .0001) and increasing duration of breastfeeding was associated with decreased BC risk (combined cohort P trend = .0003). Relative to being nulliparous, uniparous BRCA1 mutation carriers were at increased BC risk in the prospective analysis (prospective hazard ration [HRp] = 1.69, 95% CI = 1.09 to 2.62). For BRCA2 mutation carriers, being parous was associated with a 30% increase in BC risk (HRc = 1.33, 95% CI = 1.05 to 1.69), and there was no apparent decrease in risk associated with multiparity except for having at least 4 FTPs vs. 1 FTP (HRc = 0.72, 95% CI = 0.54 to 0.98).CONCLUSIONS:These findings suggest differential associations with parity between BRCA1 and BRCA2 mutation carriers with higher risk for uniparous BRCA1 carriers and parous BRCA2 carriers

Homologous recombination DNA repair defects in PALB2- associated breast cancers

Abstract: Mono-allelic germline pathogenic variants in the Partner And Localizer of BRCA2 (PALB2) gene predispose to a high-risk of breast cancer development, consistent with the role of PALB2 in homologous recombination (HR) DNA repair. Here, we sought to define the repertoire of somatic genetic alterations in PALB2-associated breast cancers (BCs), and whether PALB2-associated BCs display bi-allelic inactivation of PALB2 and/or genomic features of HR-deficiency (HRD). Twenty-four breast cancer patients with pathogenic PALB2 germline mutations were analyzed by whole-exome sequencing (WES, n = 16) or targeted capture massively parallel sequencing (410 cancer genes, n = 8). Somatic genetic alterations, loss of heterozygosity (LOH) of the PALB2 wild-type allele, large-scale state transitions (LSTs) and mutational signatures were defined. PALB2-associated BCs were found to be heterogeneous at the genetic level, with PIK3CA (29%), PALB2 (21%), TP53 (21%), and NOTCH3 (17%) being the genes most frequently affected by somatic mutations. Bi-allelic PALB2 inactivation was found in 16 of the 24 cases (67%), either through LOH (n = 11) or second somatic mutations (n = 5) of the wild-type allele. High LST scores were found in all 12 PALB2-associated BCs with bi-allelic PALB2 inactivation sequenced by WES, of which eight displayed the HRD-related mutational signature 3. In addition, bi-allelic inactivation of PALB2 was significantly associated with high LST scores. Our findings suggest that the identification of bi-allelic PALB2 inactivation in PALB2-associated BCs is required for the personalization of HR-directed therapies, such as platinum salts and/or PARP inhibitors, as the vast majority of PALB2-associated BCs without PALB2 bi-allelic inactivation lack genomic features of HRD

Polygenic risk scores and breast and epithelial ovarian cancer risks for carriers of BRCA1 and BRCA2 pathogenic variants

Purpose We assessed the associations between population-based polygenic risk scores (PRS) for breast (BC) or epithelial ovarian cancer (EOC) with cancer risks forBRCA1andBRCA2pathogenic variant carriers. Methods Retrospective cohort data on 18,935BRCA1and 12,339BRCA2female pathogenic variant carriers of European ancestry were available. Three versions of a 313 single-nucleotide polymorphism (SNP) BC PRS were evaluated based on whether they predict overall, estrogen receptor (ER)-negative, or ER-positive BC, and two PRS for overall or high-grade serous EOC. Associations were validated in a prospective cohort. Results The ER-negative PRS showed the strongest association with BC risk forBRCA1carriers (hazard ratio [HR] per standard deviation = 1.29 [95% CI 1.25-1.33],P = 3x10(-72)). ForBRCA2, the strongest association was with overall BC PRS (HR = 1.31 [95% CI 1.27-1.36],P = 7x10(-50)). HR estimates decreased significantly with age and there was evidence for differences in associations by predicted variant effects on protein expression. The HR estimates were smaller than general population estimates. The high-grade serous PRS yielded the strongest associations with EOC risk forBRCA1(HR = 1.32 [95% CI 1.25-1.40],P = 3x10(-22)) andBRCA2(HR = 1.44 [95% CI 1.30-1.60],P = 4x10(-12)) carriers. The associations in the prospective cohort were similar. Conclusion Population-based PRS are strongly associated with BC and EOC risks forBRCA1/2carriers and predict substantial absolute risk differences for women at PRS distribution extremes.Peer reviewe