Centre Inria de Lyon

By pairing to messenger RNAs (mRNAs for short), microRNAs (miRNAs) regulate gene expression in animals and plants. Accurately identifying which mRNAs interact with a given miRNA and the precise location of the interaction sites is crucial to reaching a more complete view of the regulatory network of an organism. Only a few experimental approaches, however, allow the identification of both within a single experiment. Computational predictions of miRNA-mRNA interactions thus remain generally the first step used, despite their drawback of a high rate of false-positive predictions. The major computational approaches available rely on a diversity of features, among which anchoring the miRNA seed and measuring mRNA accessibility are the key ones, with the first being universally used, while the use of the second remains controversial. Revisiting the importance of each is the aim of this paper, which uses Cross-Linking, Ligation, And Sequencing of Hybrids (CLASH) datasets to achieve this goal. Contrary to what might be expected, the results are more ambiguous regarding the use of the seed match as a feature, while accessibility appears to be a feature worth considering, indicating that, at least under some conditions, it may favour anchoring by miRNAs.

MicroRNAs (miRNAs) are small noncoding RNAs that are key players in the regulation of gene expression. In the past decade, with the increasing accessibility of high-throughput sequencing technologies, different methods have been developed to identify miRNAs, most of which rely on preexisting reference genomes. However, when a reference genome is absent or is not of high quality, such identification becomes more difficult. In this context, we developed BrumiR, an algorithm that is able to discover miRNAs directly and exclusively from small RNA (sRNA) sequencing (sRNA-seq) data. We benchmarked BrumiR with datasets encompassing animal and plant species using real and simulated sRNA-seq experiments. The results demonstrate that BrumiR reaches the highest recall for miRNA discovery, while at the same time being much faster and more efficient than the state-of-the-art tools evaluated. The latter allows BrumiR to analyze a large number of sRNA-seq experiments, from plants or animal species. Moreover, BrumiR detects additional information regarding other expressed sequences (sRNAs, isomiRs, etc.), thus maximizing the biological insight gained from sRNA-seq experiments. Additionally, when a reference genome is available, BrumiR provides a new mapping tool (BrumiR2reference) that performs an a posteriori exhaustive search to identify the precursor sequences. Finally, we also provide a machine learning classifier based on a random forest model that evaluates the sequence-derived features to further refine the prediction obtained from the BrumiR-core. The code of BrumiR and all the algorithms that compose the BrumiR toolkit are freely available at https://github.com/camoragaq/BrumiR.

AbstractMutation rates vary widely along genomes and across inheritance systems. This suggests that complex traits-resulting from the contributions of multiple determinants-might be composite in terms of the underlying mutation rates. Here we investigate through mathematical modeling whether such a heterogeneity may drive changes in a trait's architecture, especially in fluctuating environments, where phenotypic instability can be beneficial. We first identify a convexity principle related to the shape of the trait's fitness function, setting conditions under which composite architectures should be adaptive or, conversely and more commonly, should be selected against. Simulations reveal, however, that applying this principle to realistic evolving populations requires taking into account pervasive epistatic interactions that take place in the system. Indeed, the fate of a mutation affecting the architecture depends on the (epi)genetic background, which itself depends on the current architecture in the population. We tackle this problem by borrowing the adaptive dynamics framework from evolutionary ecology-where it is routinely used to deal with such resident/mutant dependencies-and find that the principle excluding composite architectures generally prevails. Yet the predicted evolutionary trajectories will typically depend on the initial architecture, possibly resulting in historical contingencies. Finally, by relaxing the large population size assumption, we unexpectedly find that not only the strength of selection on a trait's architecture but also its direction depend on population size, revealing a new occurrence of the recently identified phenomenon coined "sign inversion."

Cancer cells are known to express the Warburg effect-increased glycolysis and formation of lactic acid even in the presence of oxygen-as well as high glutamine uptake. In tumors, cancer cells are surrounded by collagen, immune cells, and neoangiogenesis. Whether collagen formation, neoangiogenesis, and inflammation in cancer are associated with the Warburg effect needs to be established. Metabolic modelling has proven to be a tool of choice to understand biological reality better and make in silico predictions. Elementary Flux Modes (EFMs) are essential for conducting an unbiased decomposition of a metabolic model into its minimal functional units. EFMs can be investigated using our tool, aspefm, an innovative approach based on logic programming where biological constraints can be incorporated. These constraints allow networks to be characterized regardless of their size. Using a metabolic model of the human cell containing collagen, neoangiogenesis, and inflammation markers, we derived a subset of EFMs of biological relevance to the Warburg effect. Within this model, EFMs analysis provided more adequate results than parsimonious flux balance analysis and flux sampling. Upon further inspection, the EFM with the best linear regression fit to cancer cell lines exometabolomics data was selected. The minimal pathway, presenting the Warburg effect, collagen synthesis, angiogenesis, and release of inflammation markers, showed that collagen production was possible directly de novo from glutamine uptake and without extracellular import of glycine and proline, collagen's main constituents.

The major economic and health consequences of COVID-19 called for various protective measures and mass vaccination campaigns. A previsional model was used to predict the future impacts of various measure combinations on COVID-19 mortality over a 400-day period in France. Calibrated on previous national hospitalization and mortality data, an agent-based epidemiological model was used to predict individual and combined effects of booster doses, vaccination of refractory adults, and vaccination of children, according to infection severity, immunity waning, and graded non-pharmaceutical interventions (NPIs). Assuming a 1.5 hospitalization hazard ratio and rapid immunity waning, booster doses would reduce COVID-19-related deaths by 50-70% with intensive NPIs and 93% with moderate NPIs. Vaccination of initially-refractory adults or children ≥5 years would half the number of deaths whatever the infection severity or degree of immunity waning. Assuming a 1.5 hospitalization hazard ratio, rapid immunity waning, moderate NPIs and booster doses, vaccinating children ≥12 years, ≥5 years, and ≥6 months would result in 6212, 3084, and 3018 deaths, respectively (vs. 87,552, 64,002, and 48,954 deaths without booster, respectively). In the same conditions, deaths would be 2696 if all adults and children ≥12 years were vaccinated and 2606 if all adults and children ≥6 months were vaccinated (vs. 11,404 and 3624 without booster, respectively). The model dealt successfully with single measures or complex combinations. It can help choosing them according to future epidemic features, vaccination extensions, and population immune status.

International audience

Abstract Mutation rates vary widely along genomes and across inheritance systems. This suggests that complex traits – resulting from the contributions of multiple determinants – might be composite in terms of the underlying mutation rates. Here we investigate through mathematical modeling whether such an heterogeneity may drive changes in a trait’s architecture, especially in fluctuating environments where phenotypic instability can be beneficial. We first identify a convexity principle, related to the shape of the trait’s fitness function, setting conditions under which composite architectures should be adaptive or, conversely and more commonly, should be selected against. Simulations reveal, however, that applying this principle to realistic evolving populations requires taking into account pervasive epistatic interactions that take place in the system. Indeed, the fate of a mutation affecting the architecture depends on the (epi)genetic background, itself depending upon the current architecture in the population. We tackle this problem by borrowing the adaptive dynamics framework from evolutionary ecology – where it is routinely used to deal with such resident/mutant dependencies – and find that the principle excluding composite architectures generally prevails. Yet, the predicted evolutionary trajectories will typically depend on the initial architecture, possibly resulting in historical contingencies. Finally, by relaxing the large population size assumption, we unexpectedly find that not only the strength of selection on a trait’s architecture, but also its direction, depend on population size, revealing a new occurrence of the recently coined phenomenon of ‘sign inversion’.