Multiple hot-deck imputation for network inference from RNA sequencing data.

Motivation: Network inference provides a global view of the relations existing between gene expression in a given transcriptomic experiment (often only for a restricted list of chosen genes). However, it is still a challenging problem: even if the cost of sequencing techniques has decreased over the last years, the number of samples in a given experiment is still (very) small compared to the number of genes. Results: We propose a method to increase the reliability of the inference when RNA-seq expression data have been measured together with an auxiliary dataset that can provide external information on gene expression similarity between samples. Our statistical approach, hd-MI, is based on imputation for samples without available RNA-seq data that are considered as missing data but are observed on the secondary dataset. hd-MI can improve the reliability of the inference for missing rates up to 30% and provides more stable networks with a smaller number of false positive edges. On a biological point of view, hd-MI was also found relevant to infer networks from RNA-seq data acquired in adipose tissue during a nutritional intervention in obese individuals. In these networks, novel links between genes were highlighted, as well as an improved comparability between the two steps of the nutritional intervention. Availability and implementation: Software and sample data are available as an R package, RNAseqNet, that can be downloaded from the Comprehensive R Archive Network (CRAN). Contact: alyssa.imbert@inra.fr or nathalie.villa-vialaneix@inra.fr. Supplementary information: Supplementary data are available at Bioinformatics online.

Protein quantitative trait locus study in obesity during weight-loss identifies a leptin regulator.

Thousands of genetic variants have been associated with complex traits through genome-wide association studies. However, the functional variants or mechanistic consequences remain elusive. Intermediate traits such as gene expression or protein levels are good proxies of the metabolic state of an organism. Proteome analysis especially can provide new insights into the molecular mechanisms of complex traits like obesity. The role of genetic variation in determining protein level variation has not been assessed in obesity. To address this, we design a large-scale protein quantitative trait locus (pQTL) analysis based on a set of 1129 proteins from 494 obese subjects before and after a weight loss intervention. This reveals 55 BMI-associated cis-pQTLs and trans-pQTLs at baseline and 3 trans-pQTLs after the intervention. We provide evidence for distinct genetic mechanisms regulating BMI-associated proteins before and after weight loss. Finally, by functional analysis, we identify and validate FAM46A as a trans regulator for leptin.

Transcriptome profiling from adipose tissue during a low-calorie diet reveals predictors of weight and glycemic outcomes in obese, nondiabetic subjects.

Background: A low-calorie diet (LCD) reduces fat mass excess, improves insulin sensitivity, and alters adipose tissue (AT) gene expression, yet the relation with clinical outcomes remains unclear.Objective: We evaluated AT transcriptome alterations during an LCD and the association with weight and glycemic outcomes both at LCD termination and 6 mo after the LCD.Design: Using RNA sequencing (RNAseq), we analyzed transcriptome changes in AT from 191 obese, nondiabetic patients within a multicenter, controlled dietary intervention. Expression changes were associated with outcomes after an 8-wk LCD (800-1000 kcal/d) and 6 mo after the LCD. Results were validated by using quantitative reverse transcriptase-polymerase chain reaction in 350 subjects from the same cohort. Statistical models were constructed to classify weight maintainers or glycemic improvers.Results: With RNAseq analyses, we identified 1173 genes that were differentially expressed after the LCD, of which 350 and 33 were associated with changes in body mass index (BMI; in kg/m2) and Matsuda index values, respectively, whereas 29 genes were associated with both endpoints. Pathway analyses highlighted enrichment in lipid and glucose metabolism. Classification models were constructed to identify weight maintainers. A model based on clinical baseline variables could not achieve any classification (validation AUC: 0.50; 95% CI: 0.36, 0.64). However, clinical changes during the LCD yielded better performance of the model (AUC: 0.73; 95% CI: 0.60, 0.87]). Adding baseline expression to this model improved the performance significantly (AUC: 0.87; 95% CI: 0.77, 0.96; Delong's P = 0.012). Similar analyses were performed to classify subjects with good glycemic improvements. Baseline- and LCD-based clinical models yielded similar performance (best AUC: 0.73; 95% CI: 0.60, 0.86). The addition of expression changes during the LCD improved the performance substantially (AUC: 0.80; 95% CI: 0.69, 0.92; P = 0.058).Conclusions: This study investigated AT transcriptome alterations after an LCD in a large cohort of obese, nondiabetic patients. Gene expression combined with clinical variables enabled us to distinguish weight and glycemic responders from nonresponders. These potential biomarkers may help clinicians understand intersubject variability and better predict the success of dietary interventions. This trial was registered at clinicaltrials.gov as NCT00390637.

Aubergine iCLIP Reveals piRNA-Dependent Decay of mRNAs Involved in Germ Cell Development in the Early Embryo.

The Piwi-interacting RNA (piRNA) pathway plays an essential role in the repression of transposons in the germline. Other functions of piRNAs such as post-transcriptional regulation of mRNAs are now emerging. Here, we perform iCLIP with the PIWI protein Aubergine (Aub) and identify hundreds of maternal mRNAs interacting with Aub in the early Drosophila embryo. Gene expression profiling reveals that a proportion of these mRNAs undergo Aub-dependent destabilization during the maternal-to-zygotic transition. Strikingly, Aub-dependent unstable mRNAs encode germ cell determinants. iCLIP with an Aub mutant that is unable to bind piRNAs confirms piRNA-dependent binding of Aub to mRNAs. Base pairing between piRNAs and mRNAs can induce mRNA cleavage and decay that are essential for embryonic development. These results suggest general regulation of maternal mRNAs by Aub and piRNAs, which plays a key developmental role in the embryo through decay and localization of mRNAs encoding germ cell determinants.

MicroRNA-Dependent Transcriptional Silencing of Transposable Elements in Drosophila Follicle Cells.

RNA interference-related silencing mechanisms concern very diverse and distinct biological processes, from gene regulation (via the microRNA pathway) to defense against molecular parasites (through the small interfering RNA and the Piwi-interacting RNA pathways). Small non-coding RNAs serve as specificity factors that guide effector proteins to ribonucleic acid targets via base-pairing interactions, to achieve transcriptional or post-transcriptional regulation. Because of the small sequence complementarity required for microRNA-dependent post-transcriptional regulation, thousands of microRNA (miRNA) putative targets have been annotated in Drosophila. In Drosophila somatic ovarian cells, genomic parasites, such as transposable elements (TEs), are transcriptionally repressed by chromatin changes induced by Piwi-interacting RNAs (piRNAs) that prevent them from invading the germinal genome. Here we show, for the first time, that a functional miRNA pathway is required for the piRNA-mediated transcriptional silencing of TEs in this tissue. Global miRNA depletion, caused by tissue- and stage-specific knock down of drosha (involved in miRNA biogenesis), AGO1 or gawky (both responsible for miRNA activity), resulted in loss of TE-derived piRNAs and chromatin-mediated transcriptional de-silencing of TEs. This specific TE de-repression was also observed upon individual titration (by expression of the complementary miRNA sponge) of two miRNAs (miR-14 and miR-34) as well as in a miR-14 loss-of-function mutant background. Interestingly, the miRNA defects differentially affected TE- and 3' UTR-derived piRNAs. To our knowledge, this is the first indication of possible differences in the biogenesis or stability of TE- and 3' UTR-derived piRNAs. This work is one of the examples of detectable phenotypes caused by loss of individual miRNAs in Drosophila and the first genetic evidence that miRNAs have a role in the maintenance of genome stability via piRNA-mediated TE repression.

Nuclear-receptor-mediated telomere insertion leads to genome instability in ALT cancers.

The breakage-fusion-bridge cycle is a classical mechanism of telomere-driven genome instability in which dysfunctional telomeres are fused to other chromosomal extremities, creating dicentric chromosomes that eventually break at mitosis. Here, we uncover a distinct pathway of telomere-driven genome instability, specifically occurring in cells that maintain telomeres with the alternative lengthening of telomeres mechanism. We show that, in these cells, telomeric DNA is added to multiple discrete sites throughout the genome, corresponding to regions regulated by NR2C/F transcription factors. These proteins drive local telomere DNA addition by recruiting telomeric chromatin. This mechanism, which we name targeted telomere insertion (TTI), generates potential common fragile sites that destabilize the genome. We propose that TTI driven by NR2C/F proteins contributes to the formation of complex karyotypes in ALT tumors.

A user-friendly chromatographic method to purify small regulatory RNAs.

The discovery of the small regulatory RNAs has changed our vision of cellular regulations. Indeed, when loaded on Argonaute proteins they form ribonucleoprotein complexes (RNPs) that target complementary sequences to achieve widespread silencing mechanisms conserved in most eukaryotes. The recent development of deep sequencing approaches highly contributed to their detection. Small RNA isolation from cells and/or tissues remains a crucial stage to generate robust and relevant sequencing data. In 2006, a novel strategy based on anion-exchange chromatography has been proposed as an alternative to the standard size-isolation purification procedure. Using bioinformatic comparative analysis, we here demonstrate that anion-exchange chromatographic RNP purification prior to small RNA extraction unbiasedly enriches datasets in bona fide reads (small regulatory RNA sequences) and depletes endogenous contaminants (ribosomal RNAs and degradation RNA products). The resulting increase in sequencing depth provides a major benefit to study rare populations. We then developed a fast and basic manual procedure to purify such small non-coding RNAs using anion-exchange chromatography at the bench. We validated the efficiency of this new method and used this strategy to purify small RNAs from various tissues and organisms. We moreover determined that our manual purification increases the output of the previously described anion-exchange chromatography procedure.

Maternally deposited germline piRNAs silence the tirant retrotransposon in somatic cells.

Transposable elements (TEs), whose propagation can result in severe damage to the host genome, are silenced in the animal gonad by Piwi-interacting RNAs (piRNAs). piRNAs produced in the ovaries are deposited in the embryonic germline and initiate TE repression in the germline progeny. Whether the maternally transmitted piRNAs play a role in the silencing of somatic TEs is however unknown. Here we show that maternally transmitted piRNAs from the tirant retrotransposon in Drosophila are required for the somatic silencing of the TE and correlate with an increase in histone H3K9 trimethylation an active tirant copy.

piRNA-mediated transgenerational inheritance of an acquired trait.

The maintenance of genome integrity is an essential trait to the successful transmission of genetic information. In animal germ cells, piRNAs guide PIWI proteins to silence transposable elements (TEs) in order to maintain genome integrity. In insects, most TE silencing in the germline is achieved by secondary piRNAs that are produced by a feed-forward loop (the ping-pong cycle), which requires the piRNA-directed cleavage of two types of RNAs: mRNAs of functional euchromatic TEs and heterochromatic transcripts that contain defective TE sequences. The first cleavage that initiates such an amplification loop remains poorly understood. Taking advantage of the existence of strains that are devoid of functional copies of the LINE-like I-element, we report here that in such Drosophila ovaries, the initiation of a ping-pong cycle is exclusively achieved by secondary I-element piRNAs that are produced in the ovary and deposited in the embryonic germline. This unusual secondary piRNA biogenesis, detected in the absence of functional I-element copies, results from the processing of sense and antisense transcripts of several different defective I-element. Once acquired, for instance after ancestor aging, this capacity to produce heterochromatic-only secondary piRNAs is partially transmitted through generations via maternal piRNAs. Furthermore, such piRNAs acting as ping-pong initiators in a chromatin-independent manner confer to the progeny a high capacity to repress the I-element mobility. Our study explains, at the molecular level, the basis for epigenetic memory of maternal immunity that protects females from hybrid dysgenesis caused by transposition of paternally inherited functional I-element.

Phylogenetic study of plant Q-type C2H2 zinc finger proteins and expression analysis of poplar genes in response to osmotic, cold and mechanical stresses.

Plant Q-type C2H2 zinc finger transcription factors play an important role in plant tolerance to various environmental stresses such as drought, cold, osmotic stress, wounding and mechanical loading. To carry out an improved analysis of the specific role of each member of this subfamily in response to mechanical loading in poplar, we identified 16 two-fingered Q-type C2H2-predicted proteins from the poplar Phytozome database and compared their phylogenetic relationships with 152 two-fingered Q-type C2H2 protein sequences belonging to more than 50 species isolated from the NR protein database of NCBI. Phylogenetic analyses of these Q-type C2H2 proteins sequences classified them into two groups G1 and G2, and conserved motif distributions of interest were established. These two groups differed essentially in their signatures at the C-terminus of their two QALGGH DNA-binding domains. Two additional conserved motifs, MALEAL and LVDCHY, were found only in sequences from Group G1 or from Group G2, respectively. Functional significance of these phylogenetic divergences was assessed by studying transcript accumulation of six poplar C2H2 Q-type genes in responses to abiotic stresses; but no group specificity was found in any organ. Further expression analyses focused on PtaZFP1 and PtaZFP2, the two genes strongly induced by mechanical loading in poplars. The results revealed that these two genes were regulated by several signalling molecules including hydrogen peroxide and the phytohormone jasmonate.