The DEAD box RNA helicase DDX42 is an intrinsic inhibitor of positive-strand RNA viruses.

Genome-wide screens are powerful approaches to unravel regulators of viral infections. Here, a CRISPR screen identifies the RNA helicase DDX42 as an intrinsic antiviral inhibitor of HIV-1. Depletion of endogenous DDX42 increases HIV-1 DNA accumulation and infection in cell lines and primary cells. DDX42 overexpression inhibits HIV-1 infection, whereas expression of a dominant-negative mutant increases infection. Importantly, DDX42 also restricts LINE-1 retrotransposition and infection with other retroviruses and positive-strand RNA viruses, including CHIKV and SARS-CoV-2. However, DDX42 does not impact the replication of several negative-strand RNA viruses, arguing against an unspecific effect on target cells, which is confirmed by RNA-seq analysis. Proximity ligation assays show DDX42 in the vicinity of viral elements, and cross-linking RNA immunoprecipitation confirms a specific interaction of DDX42 with RNAs from sensitive viruses. Moreover, recombinant DDX42 inhibits HIV-1 reverse transcription in vitro. Together, our data strongly suggest a direct mode of action of DDX42 on viral ribonucleoprotein complexes. Our results identify DDX42 as an intrinsic viral inhibitor, opening new perspectives to target the life cycle of numerous RNA viruses.

Transcriptomic profiles of consistent risk-taking behaviour across time and contexts in European sea bass.

Bolder individuals have greater access to food sources and reproductive partners but are also at increased risk of predation. Boldness is believed to be consistent across time and contexts, but few studies have investigated the stability of this trait across variable environments, such as varying stress loads or long periods of time. Moreover, the underlying molecular components of boldness are poorly studied. Here, we report that boldness of 1154 European sea bass, evaluated using group risk-taking tests, is consistent over seven months and for individuals subjected to multiple environments, including a chronically stressful environment. Differences in risk-taking behaviour were further supported by differences observed in the responses to a novel environment test: shy individuals displayed more group dispersion, more thigmotaxic behaviour and lower activity levels. Transcriptomic analyses performed on extreme phenotypes revealed that bold individuals display greater expression for genes involved in social and exploration behaviours, and memory in the pituitary, and genes involved in immunity and responses to stimuli in the head kidney. This study demonstrates that personality traits come with an underpinning molecular signature, especially in organs involved in the endocrine and immune systems. As such, our results help to depict state-behaviour feedback mechanisms, previously proposed as key in shaping animal personality.

ISoLDE: a data-driven statistical method for the inference of allelic imbalance in datasets with reciprocal crosses.

MOTIVATION: Allelic imbalance (AI), i.e. the unequal expression of the alleles of the same gene in a single cell, affects a subset of genes in diploid organisms. One prominent example of AI is parental genomic imprinting, which results in parent-of-origin-dependent, mono-allelic expression of a limited number of genes in metatherian and eutherian mammals and in angiosperms. Currently available methods for identifying AI rely on data modeling and come with the associated limitations. RESULTS: We have designed ISoLDE (Integrative Statistics of alleLe Dependent Expression), a novel nonparametric statistical method that takes into account both AI and the characteristics of RNA-seq data to infer allelic expression bias when at least two biological replicates are available for reciprocal crosses. ISoLDE learns the distribution of a specific test statistic from the data and calls genes 'allelically imbalanced', 'bi-allelically expressed' or 'undetermined'. Depending on the number of replicates, predefined thresholds or permutations are used to make calls. We benchmarked ISoLDE against published methods, and showed that ISoLDE compared favorably with respect to sensitivity, specificity and robustness to the number of replicates. Using ISoLDE on different RNA-seq datasets generated from hybrid mouse tissues, we did not discover novel imprinted genes (IGs), confirming the most conservative estimations of IG number. AVAILABILITY AND IMPLEMENTATION: ISoLDE has been implemented as a Bioconductor package available at http://bioconductor.org/packages/ISoLDE/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of Plagl1/Zac1 binding sites and target genes establishes its role in the regulation of extracellular matrix genes and the imprinted gene network.

PLAGL1/ZAC1 undergoes parental genomic imprinting, is paternally expressed, and is a member of the imprinted gene network (IGN). It encodes a zinc finger transcription factor with anti-proliferative activity and is a candidate tumor suppressor gene on 6q24 whose expression is frequently lost in various neoplasms. Conversely, gain of PLAGL1 function is responsible for transient neonatal diabetes mellitus, a rare genetic disease that results from defective pancreas development. In the present work, we showed that Plagl1 up-regulation was not associated with DNA damage-induced cell cycle arrest. It was rather associated with physiological cell cycle exit that occurred with contact inhibition, growth factor withdrawal, or cell differentiation. To gain insights into Plagl1 mechanism of action, we identified Plagl1 target genes by combining chromatin immunoprecipitation and genome-wide transcriptomics in transfected cell lines. Plagl1-elicited gene regulation correlated with multiple binding to the proximal promoter region through a GC-rich motif. Plagl1 target genes included numerous genes involved in signaling, cell adhesion, and extracellular matrix composition, including collagens. Plagl1 targets also included 22% of the 409 genes that make up the IGN. Altogether, this work identified Plagl1 as a transcription factor that coordinated the regulation of a subset of IGN genes and controlled extracellular matrix composition.

Changes in chromatin state reveal ARNT2 at a node of a tumorigenic transcription factor signature driving glioblastoma cell aggressiveness.

Although a growing body of evidence indicates that phenotypic plasticity exhibited by glioblastoma cells plays a central role in tumor development and post-therapy recurrence, the master drivers of their aggressiveness remain elusive. Here we mapped the changes in active (H3K4me3) and repressive (H3K27me3) histone modifications accompanying the repression of glioblastoma stem-like cells tumorigenicity. Genes with changing histone marks delineated a network of transcription factors related to cancerous behavior, stem state, and neural development, highlighting a previously unsuspected association between repression of ARNT2 and loss of cell tumorigenicity. Immunohistochemistry confirmed ARNT2 expression in cell sub-populations within proliferative zones of patients' glioblastoma. Decreased ARNT2 expression was consistently observed in non-tumorigenic glioblastoma cells, compared to tumorigenic cells. Moreover, ARNT2 expression correlated with a tumorigenic molecular signature at both the tissue level within the tumor core and at the single cell level in the patients' tumors. We found that ARNT2 knockdown decreased the expression of SOX9, POU3F2 and OLIG2, transcription factors implicated in glioblastoma cell tumorigenicity, and repressed glioblastoma stem-like cell tumorigenic properties in vivo. Our results reveal ARNT2 as a pivotal component of the glioblastoma cell tumorigenic signature, located at a node of a transcription factor network controlling glioblastoma cell aggressiveness.

In Vitro Corticogenesis from Embryonic Stem Cells Recapitulates the In Vivo Epigenetic Control of Imprinted Gene Expression.

In vitro corticogenesis from embryonic stem cells (ESCs) is an attractive model of cortical development and a promising tool for cortical therapy. It is unknown to which extent epigenetic mechanisms crucial for cortex development and function, such as parental genomic imprinting, are recapitulated by in vitro corticogenesis. Here, using genome-wide transcriptomic and methylation analyses on hybrid mouse tissues and cells, we find a high concordance of imprinting status between in vivo and ESC-derived cortices. Notably, in vitro corticogenesis strictly reproduced the in vivo parent-of-origin-dependent expression of 41 imprinted genes (IGs), including Mest and Cdkn1c known to control corticogenesis. Parent-of-origin-dependent DNA methylation was also conserved at 14 of 18 imprinted differentially methylated regions. The least concordant imprinted locus was Gpr1-Zdbf2, where the aberrant bi-allelic expression of Zdbf2 and Adam23 was concomitant with a gain of methylation on the maternal allele in vitro. Combined, our data argue for a broad conservation of the epigenetic mechanisms at imprinted loci in cortical cells derived from ESCs. We propose that in vitro corticogenesis helps to define the still poorly understood mechanisms that regulate imprinting in the brain and the roles of IGs in cortical development.

Deciphering the Wisent Demographic and Adaptive Histories from Individual Whole-Genome Sequences.

As the largest European herbivore, the wisent (Bison bonasus) is emblematic of the continent wildlife but has unclear origins. Here, we infer its demographic and adaptive histories from two individual whole-genome sequences via a detailed comparative analysis with bovine genomes. We estimate that the wisent and bovine species diverged from 1.7 × 106 to 850,000 years before present (YBP) through a speciation process involving an extended period of limited gene flow. Our data further support the occurrence of more recent secondary contacts, posterior to the Bos taurus and Bos indicus divergence (∼150,000 YBP), between the wisent and (European) taurine cattle lineages. Although the wisent and bovine population sizes experienced a similar sharp decline since the Last Glacial Maximum, we find that the wisent demography remained more fluctuating during the Pleistocene. This is in agreement with a scenario in which wisents responded to successive glaciations by habitat fragmentation rather than southward and eastward migration as for the bovine ancestors. We finally detect 423 genes under positive selection between the wisent and bovine lineages, which shed a new light on the genome response to different living conditions (temperature, available food resource, and pathogen exposure) and on the key gene functions altered by the domestication process.

A duplicated NUCLEOLIN gene with antagonistic activity is required for chromatin organization of silent 45S rDNA in Arabidopsis.

In plants as well as in animals, hundreds to thousands of 45S rRNA gene copies localize in Nucleolus Organizer Regions (NORs), and the activation or repression of specific sets of rDNA depends on epigenetic mechanisms. Previously, we reported that the Arabidopsis thaliana nucleolin protein NUC1, an abundant and evolutionarily conserved nucleolar protein in eukaryotic organisms, is required for maintaining DNA methylation levels and for controlling the expression of specific rDNA variants in Arabidopsis. Interestingly, in contrast with animal or yeast cells, plants contain a second nucleolin gene. Here, we report that Arabidopsis NUC1 and NUC2 nucleolin genes are both required for plant growth and survival and that NUC2 disruption represses flowering. However, these genes seem to be functionally antagonistic. In contrast with NUC1, disruption of NUC2 induces CG hypermethylation of rDNA and NOR association with the nucleolus. Moreover, NUC2 loss of function triggers major changes in rDNA spatial organization, expression, and transgenerational stability. Our analyses indicate that silencing of specific rRNA genes is mostly determined by the active or repressed state of the NORs and that nucleolin proteins play a key role in the developmental control of this process.

FliZ is a global regulatory protein affecting the expression of flagellar and virulence genes in individual Xenorhabdus nematophila bacterial cells.

Heterogeneity in the expression of various bacterial genes has been shown to result in the presence of individuals with different phenotypes within clonal bacterial populations. The genes specifying motility and flagellar functions are coordinately regulated and form a complex regulon, the flagellar regulon. Complex interplay has recently been demonstrated in the regulation of flagellar and virulence gene expression in many bacterial pathogens. We show here that FliZ, a DNA-binding protein, plays a key role in the insect pathogen, Xenorhabdus nematophila, affecting not only hemolysin production and virulence in insects, but efficient swimming motility. RNA-Seq analysis identified FliZ as a global regulatory protein controlling the expression of 278 Xenorhabdus genes either directly or indirectly. FliZ is required for the efficient expression of all flagellar genes, probably through its positive feedback loop, which controls expression of the flhDC operon, the master regulator of the flagellar circuit. FliZ also up- or downregulates the expression of numerous genes encoding non-flagellar proteins potentially involved in key steps of the Xenorhabdus lifecycle. Single-cell analysis revealed the bimodal expression of six identified markers of the FliZ regulon during exponential growth of the bacterial population. In addition, a combination of fluorescence-activated cell sorting and RT-qPCR quantification showed that this bimodality generated a mixed population of cells either expressing ("ON state") or not expressing ("OFF state") FliZ-dependent genes. Moreover, studies of a bacterial population exposed to a graded series of FliZ concentrations showed that FliZ functioned as a rheostat, controlling the rate of transition between the "OFF" and "ON" states in individuals. FliZ thus plays a key role in cell fate decisions, by transiently creating individuals with different potentials for motility and host interactions.

In vivo investigation of Lcr35® anti-candidiasis properties in Caenorhabditis elegans reveals the involvement of highly conserved immune pathways.

Lactic acid bacteria, including the microorganisms formerly designated as Lactobacillus, are the major representatives of Live Biotherapeutic Microorganisms (LBM) when used for therapeutic purposes. However, in most cases, the mechanisms of action remain unknown. The antifungal potential of LBM has already been demonstrated using preclinical models (cell cultures, laboratory animals). Understanding their mechanisms of action is strategic for the development of new therapeutics for humans. Here, Caenorhabditis elegans was used as an in vivo model to analyze pro-longevity, anti-aging and anti-candidiasis effects of the LBM Lacticaseibacillus rhamnosus (formerly Lactobacillus rhamnosus) Lcr35®. A high-throughput transcriptomic analysis revealed a specific response of C. elegans depending on whether it is in the presence of the LBM L. rhamnosus Lcr35® (structural response), the yeast Candida albicans (metabolic response) or both (structural and metabolic responses) in a preventive and a curative conditions. Studies on C. elegans mutants demonstrated that the p38 MAPK (sek-1, skn-1) and the insulin-like (daf-2, daf-16) signaling pathways were involved in the extended lifespan provided by L. rhamnosus Lcr35® strain whereas the JNK pathway was not involved (jnk-1). In addition, the anti C. albicans effect of the bacterium requires the daf-16 and sek-1 genes while it is independent of daf-2 and skn-1. Moreover, the anti-aging effect of Lcr35®, linked to the extension of longevity, is not due to protection against oxidative stress (H2O2). Taken together, these results formally show the involvement of the p38 MAP kinase and insulin-like signaling pathways for the longevity extension and anti-Candida albicans properties of Lcr35® with, however, differences in the genes involved. Overall, these findings provide new insight for understanding the mechanisms of action of a probiotic strain with antimicrobial potential.

Computing biological functions using BioPsi, a formal description of biological processes based on elementary bricks of actions.

MOTIVATION: In the available databases, biological processes are described from molecular and cellular points of view, but these descriptions are represented with text annotations that make it difficult to handle them for computation. Consequently, there is an obvious need for formal descriptions of biological processes. RESULTS: We present a formalism that uses the BioPsi concepts to model biological processes from molecular details to networks. This computational approach, based on elementary bricks of actions, allows us to calculate on biological functions (e.g. process comparison, mapping structure-function relationships, etc.). We illustrate its application with two examples: the functional comparison of proteases and the functional description of the glycolysis network. This computational approach is compatible with detailed biological knowledge and can be applied to different kinds of systems of simulation. AVAILABILITY: www.sysdiag.cnrs.fr/publications/supplementary-materials/BioPsi_Manager/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

BioNetCAD: design, simulation and experimental validation of synthetic biochemical networks.

MOTIVATION: Synthetic biology studies how to design and construct biological systems with functions that do not exist in nature. Biochemical networks, although easier to control, have been used less frequently than genetic networks as a base to build a synthetic system. To date, no clear engineering principles exist to design such cell-free biochemical networks. RESULTS: We describe a methodology for the construction of synthetic biochemical networks based on three main steps: design, simulation and experimental validation. We developed BioNetCAD to help users to go through these steps. BioNetCAD allows designing abstract networks that can be implemented thanks to CompuBioTicDB, a database of parts for synthetic biology. BioNetCAD enables also simulations with the HSim software and the classical Ordinary Differential Equations (ODE). We demonstrate with a case study that BioNetCAD can rationalize and reduce further experimental validation during the construction of a biochemical network. AVAILABILITY AND IMPLEMENTATION: BioNetCAD is freely available at http://www.sysdiag.cnrs.fr/BioNetCAD. It is implemented in Java and supported on MS Windows. CompuBioTicDB is freely accessible at http://compubiotic.sysdiag.cnrs.fr/.

Differential immune gene expression associated with contemporary range expansion in two invasive rodents in Senegal.

Biological invasions are major anthropogenic changes associated with threats to biodiversity and health. However, what determines the successful establishment and spread of introduced populations remains unclear. Here, we explore several hypotheses linking invasion success and immune phenotype traits, including those based on the evolution of increased competitive ability concept. We compared gene expression profiles between anciently and recently established populations of two major invading species, the house mouse Mus musculus domesticus and the black rat Rattus rattus, in Senegal (West Africa). Transcriptome analyses identified differential expression between anciently and recently established populations for 364 mouse genes and 83 rat genes. All immune-related genes displaying differential expression along the mouse invasion route were overexpressed at three of the four recently invaded sites studied. Complement activation pathway genes were overrepresented among these genes. By contrast, no particular immunological process was found to be overrepresented among the differentially expressed genes of black rat. Changes in transcriptome profiles were thus observed along invasion routes, but with different specific patterns between the two invasive species. These changes may be driven by increases in infection risks at sites recently invaded by the house mouse, and by stochastic events associated with colonization history for the black rat. These results constitute a first step toward the identification of immune eco-evolutionary processes potentially involved in the invasion success of these two rodent species.

Functional and molecular heterogeneity of D2R neurons along dorsal ventral axis in the striatum.

Action control is a key brain function determining the survival of animals in their environment. In mammals, neurons expressing dopamine D2 receptors (D2R) in the dorsal striatum (DS) and the nucleus accumbens (Acb) jointly but differentially contribute to the fine regulation of movement. However, their region-specific molecular features are presently unknown. By combining RNAseq of striatal D2R neurons and histological analyses, we identified hundreds of novel region-specific molecular markers, which may serve as tools to target selective subpopulations. As a proof of concept, we characterized the molecular identity of a subcircuit defined by WFS1 neurons and evaluated multiple behavioral tasks after its temporally-controlled deletion of D2R. Consequently, conditional D2R knockout mice displayed a significant reduction in digging behavior and an exacerbated hyperlocomotor response to amphetamine. Thus, targeted molecular analyses reveal an unforeseen heterogeneity in D2R-expressing striatal neuronal populations, underlying specific D2R's functional features in the control of specific motor behaviors.

Transcriptome Profiles of Nod Factor-independent Symbiosis in the Tropical Legume Aeschynomene evenia.

Nod factors (NF) were assumed to be indispensable for the establishment of a rhizobium-legume symbiosis until the discovery that certain Bradyrhizobium strains interacting with certain Aeschynomene species lack the canonical nodABC genes required for their synthesis. So far, the molecular dialogue between Aeschynomene and its symbionts remains an open question. Here we report a time course transcriptional analysis of Aeschynomene evenia in response to inoculation with Bradyrhizobium ORS278. The NF-independent symbiotic process was monitored at five time points between bacterial infection and nodule maturity. The five time points correspond to three specific events, root infection by crack entry, nodule organogenesis, and the establishment of the nitrogen fixing process. During the third stage, about 80 NCR-like genes and eight symbiotic genes known to be involved in signaling, bacterial infection or nodulation regulation were highly expressed. Comparative gene expression analyses at the five time points also enabled the selection of genes with an expression profile that makes them promising markers to monitor early plant responses to bacteria. Such markers could be used in bioassays to identify the nature of the bacterial signal(s). Our data represent valuable resources for investigation of this Nod factor-independent symbiosis.

Comparative Genomics of Glossina palpalis gambiensis and G. morsitans morsitans to Reveal Gene Orthologs Involved in Infection by Trypanosoma brucei gambiense.

Blood-feeding Glossina palpalis gambiense (Gpg) fly transmits the single-celled eukaryotic parasite Trypanosoma brucei gambiense (Tbg), the second Glossina fly African trypanosome pair being Glossina morsitans/T.brucei rhodesiense. Whatever the T. brucei subspecies, whereas the onset of their developmental program in the zoo-anthropophilic blood feeding flies does unfold in the fly midgut, its completion is taking place in the fly salivary gland where does emerge a low size metacyclic trypomastigote population displaying features that account for its establishment in mammals-human individuals included. Considering that the two Glossina-T. brucei pairs introduced above share similarity with respect to the developmental program of this African parasite, we were curious to map on the Glossina morsitans morsitans (Gmm), the Differentially Expressed Genes (DEGs) we listed in a previous study. Briefly, using the gut samples collected at days 3, 10, and 20 from Gpg that were fed or not at day 0 on Tbg-hosting mice, these DGE lists were obtained from RNA seq-based approaches. Here, post the mapping on the quality controlled DEGs on the Gmm genome, the identified ortholog genes were further annotated, the resulting datasets being compared. Around 50% of the Gpg DEGs were shown to have orthologs in the Gmm genome. Under one of the three Glossina midgut sampling conditions, the number of DEGs was even higher when mapping on the Gmm genome than initially recorded. Many Gmm genes annotated as "Hypothetical" were mapped and annotated on many distinct databases allowing some of them to be properly identified. We identify Glossina fly candidate genes encoding (a) a broad panel of proteases as well as (b) chitin-binding proteins, (c) antimicrobial peptide production-Pro3 protein, transferrin, mucin, atttacin, cecropin, etc-to further select in functional studies, the objectives being to probe and validated fly genome manipulation that prevents the onset of the developmental program of one or the other T. brucei spp. stumpy form sampled by one of the other bloodfeeding Glossina subspecies.

Cell Type-Specific mRNA Dysregulation in Hippocampal CA1 Pyramidal Neurons of the Fragile X Syndrome Mouse Model.

Fragile X syndrome (FXS) is a genetic disorder due to the silencing of the Fmr1 gene, causing intellectual disability, seizures, hyperactivity, and social anxiety. All these symptoms result from the loss of expression of the RNA binding protein fragile X mental retardation protein (FMRP), which alters the neurodevelopmental program to abnormal wiring of specific circuits. Aberrant mRNAs translation associated with the loss of Fmr1 product is widely suspected to be in part the cause of FXS. However, precise gene expression changes involved in this disorder have yet to be defined. The objective of this study was to identify the set of mistranslated mRNAs that could contribute to neurological deficits in FXS. We used the RiboTag approach and RNA sequencing to provide an exhaustive listing of genes whose mRNAs are differentially translated in hippocampal CA1 pyramidal neurons as the integrative result of FMRP loss and subsequent neurodevelopmental adaptations. Among genes differentially regulated between adult WT and Fmr1-/y mice, we found enrichment in FMRP-binders but also a majority of non-FMRP-binders. Interestingly, both up- and down-regulation of specific gene expression is relevant to fully understand the molecular deficiencies triggering FXS. More importantly, functional genomic analysis highlighted the importance of genes involved in neuronal connectivity. Among them, we show that Klk8 altered expression participates in the abnormal hippocampal dendritic spine maturation observed in a mouse model of FXS.

Ribosomal Protein S6 Phosphorylation Is Involved in Novelty-Induced Locomotion, Synaptic Plasticity and mRNA Translation.

The phosphorylation of the ribosomal protein S6 (rpS6) is widely used to track neuronal activity. Although it is generally assumed that rpS6 phosphorylation has a stimulatory effect on global protein synthesis in neurons, its exact biological function remains unknown. By using a phospho-deficient rpS6 knockin mouse model, we directly tested the role of phospho-rpS6 in mRNA translation, plasticity and behavior. The analysis of multiple brain areas shows for the first time that, in neurons, phospho-rpS6 is dispensable for overall protein synthesis. Instead, we found that phospho-rpS6 controls the translation of a subset of mRNAs in a specific brain region, the nucleus accumbens (Acb), but not in the dorsal striatum. We further show that rpS6 phospho-mutant mice display altered long-term potentiation (LTP) in the Acb and enhanced novelty-induced locomotion. Collectively, our findings suggest a previously unappreciated role of phospho-rpS6 in the physiology of the Acb, through the translation of a selective subclass of mRNAs, rather than the regulation of general protein synthesis.

Neural stem cell-encoded temporal patterning delineates an early window of malignant susceptibility in Drosophila.

Pediatric neural tumors are often initiated during early development and can undergo very rapid transformation. However, the molecular basis of this early malignant susceptibility remains unknown. During Drosophila development, neural stem cells (NSCs) divide asymmetrically and generate intermediate progenitors that rapidly differentiate in neurons. Upon gene inactivation, these progeny can dedifferentiate and generate malignant tumors. Here, we find that intermediate progenitors are prone to malignancy only when born during an early window of development while expressing the transcription factor Chinmo, and the mRNA-binding proteins Imp/IGF2BP and Lin-28. These genes compose an oncogenic module that is coopted upon dedifferentiation of early-born intermediate progenitors to drive unlimited tumor growth. In late larvae, temporal transcription factor progression in NSCs silences the module, thereby limiting mitotic potential and terminating the window of malignant susceptibility. Thus, this study identifies the gene regulatory network that confers malignant potential to neural tumors with early developmental origins.

Regulatory divergence between parental alleles determines gene expression patterns in hybrids.

Both hybridization and allopolyploidization generate novel phenotypes by conciliating divergent genomes and regulatory networks in the same cellular context. To understand the rewiring of gene expression in hybrids, the total expression of 21,025 genes and the allele-specific expression of over 11,000 genes were quantified in interspecific hybrids and their parental species, Coffea canephora and Coffea eugenioides using RNA-seq technology. Between parental species, cis- and trans-regulatory divergences affected around 32% and 35% of analyzed genes, respectively, with nearly 17% of them showing both. The relative importance of trans-regulatory divergences between both species could be related to their low genetic divergence and perennial habit. In hybrids, among divergently expressed genes between parental species and hybrids, 77% was expressed like one parent (expression level dominance), including 65% like C. eugenioides. Gene expression was shown to result from the expression of both alleles affected by intertwined parental trans-regulatory factors. A strong impact of C. eugenioides trans-regulatory factors on the upregulation of C. canephora alleles was revealed. The gene expression patterns appeared determined by complex combinations of cis- and trans-regulatory divergences. In particular, the observed biased expression level dominance seemed to be derived from the asymmetric effects of trans-regulatory parental factors on regulation of alleles. More generally, this study illustrates the effects of divergent trans-regulatory parental factors on the gene expression pattern in hybrids. The characteristics of the transcriptional response to hybridization appear to be determined by the compatibility of gene regulatory networks and therefore depend on genetic divergences between the parental species and their evolutionary history.