Major transcriptomic, epigenetic and metabolic changes underlie the pluripotency continuum in rabbit preimplantation embryos.

Despite the growing interest in the rabbit model for developmental and stem cell biology, the characterization of embryos at the molecular level is still poorly documented. We conducted a transcriptome analysis of rabbit preimplantation embryos from E2.7 (morula stage) to E6.6 (early primitive streak stage) using bulk and single-cell RNA-sequencing. In parallel, we studied oxidative phosphorylation and glycolysis, and analysed active and repressive epigenetic modifications during blastocyst formation and expansion. We generated a transcriptomic, epigenetic and metabolic map of the pluripotency continuum in rabbit preimplantation embryos, and identified novel markers of naive pluripotency that might be instrumental for deriving naive pluripotent stem cell lines. Although the rabbit is evolutionarily closer to mice than to primates, we found that the transcriptome of rabbit epiblast cells shares common features with those of humans and non-human primates.

HEXIM1 and NEAT1 Long Non-coding RNA Form a Multi-subunit Complex that Regulates DNA-Mediated Innate Immune Response.

The DNA-mediated innate immune response underpins anti-microbial defenses and certain autoimmune diseases. Here we used immunoprecipitation, mass spectrometry, and RNA sequencing to identify a ribonuclear complex built around HEXIM1 and the long non-coding RNA NEAT1 that we dubbed the HEXIM1-DNA-PK-paraspeckle components-ribonucleoprotein complex (HDP-RNP). The HDP-RNP contains DNA-PK subunits (DNAPKc, Ku70, and Ku80) and paraspeckle proteins (SFPQ, NONO, PSPC1, RBM14, and MATRIN3). We show that binding of HEXIM1 to NEAT1 is required for its assembly. We further demonstrate that the HDP-RNP is required for the innate immune response to foreign DNA, through the cGAS-STING-IRF3 pathway. The HDP-RNP interacts with cGAS and its partner PQBP1, and their interaction is remodeled by foreign DNA. Remodeling leads to the release of paraspeckle proteins, recruitment of STING, and activation of DNAPKc and IRF3. Our study establishes the HDP-RNP as a key nuclear regulator of DNA-mediated activation of innate immune response through the cGAS-STING pathway.

Single cell transcriptomics of cerebrospinal fluid cells from patients with recent-onset narcolepsy.

Narcolepsy is a rare cause of hypersomnolence and may be associated or not with cataplexy, i.e. sudden muscle weakness. These forms are designated narcolepsy-type 1 (NT1) and -type 2 (NT2), respectively. Notable characteristics of narcolepsy are that most patients carry the HLA-DQB1*06:02 allele and NT1-patients have strongly decreased levels of hypocretin-1 (synonym orexin-A) in the cerebrospinal fluid (CSF). The pathogenesis of narcolepsy is still not completely understood but the strong HLA-bias and increased frequencies of CD4+ T cells reactive to hypocretin in the peripheral blood suggest autoimmune processes in the hypothalamus. Here we analyzed the transcriptomes of CSF-cells from twelve NT1 and two NT2 patients by single cell RNAseq (scRNAseq). As controls, we used CSF cells from patients with multiple sclerosis, radiologically isolated syndrome, and idiopathic intracranial hypertension. From 27,255 CSF cells, we identified 20 clusters of different cell types and found significant differences in three CD4+ T cell and one monocyte clusters between narcolepsy and multiple sclerosis patients. Over 1000 genes were differentially regulated between patients with NT1 and other diseases. Surprisingly, the most strongly upregulated genes in narcolepsy patients as compared to controls were coding for the genome-encoded MTRNR2L12 and MTRNR2L8 peptides, which are homologous to the mitochondria-encoded HUMANIN peptide that is known playing a role in other neurological diseases including Alzheimer's disease.

H3K9me2 genome-wide distribution in the holocentric insect Spodoptera frugiperda (Lepidoptera: Noctuidae).

BACKGROUND: Eukaryotic genomes are packaged by Histone proteins in a structure called chromatin. There are different chromatin types. Euchromatin is typically associated with decondensed, transcriptionally active regions and heterochromatin to more condensed regions of the chromosomes. Methylation of Lysine 9 of Histone H3 (H3K9me) is a conserved biochemical marker of heterochromatin. In many organisms, heterochromatin is usually localized at telomeric as well as pericentromeric regions but can also be found at interstitial chromosomal loci. This distribution may vary in different species depending on their general chromosomal organization. Holocentric species such as Spodoptera frugiperda (Lepidoptera: Noctuidae) possess dispersed centromeres instead of a monocentric one and thus no observable pericentromeric compartment. To identify the localization of heterochromatin in such species we performed ChIP-Seq experiments and analyzed the distribution of the heterochromatin marker H3K9me2 in the Sf9 cell line and whole 4th instar larvae (L4) in relation to RNA-Seq data. RESULTS: In both samples we measured an enrichment of H3K9me2 at the (sub) telomeres, rDNA loci, and satellite DNA sequences, which could represent dispersed centromeric regions. We also observed that density of H3K9me2 is positively correlated with transposable elements and protein-coding genes. But contrary to most model organisms, H3K9me2 density is not correlated with transcriptional repression. CONCLUSION: This is the first genome-wide ChIP-Seq analysis conducted in S. frugiperda for H3K9me2. Compared to model organisms, this mark is found in expected chromosomal compartments such as rDNA and telomeres. However, it is also localized at numerous dispersed regions, instead of the well described large pericentromeric domains, indicating that H3K9me2 might not represent a classical heterochromatin marker in Lepidoptera. (242 words).

TALC: Transcript-level Aware Long-read Correction.

MOTIVATION: Long-read sequencing technologies are invaluable for determining complex RNA transcript architectures but are error-prone. Numerous 'hybrid correction' algorithms have been developed for genomic data that correct long reads by exploiting the accuracy and depth of short reads sequenced from the same sample. These algorithms are not suited for correcting more complex transcriptome sequencing data. RESULTS: We have created a novel reference-free algorithm called Transcript-level Aware Long-Read Correction (TALC) which models changes in RNA expression and isoform representation in a weighted De Bruijn graph to correct long reads from transcriptome studies. We show that transcript-level aware correction by TALC improves the accuracy of the whole spectrum of downstream RNA-seq applications and is thus necessary for transcriptome analyses that use long read technology. AVAILABILITY AND IMPLEMENTATION: TALC is implemented in C++ and available at https://github.com/lbroseus/TALC. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Systemic control of nodule formation by plant nitrogen demand requires autoregulation-dependent and independent mechanisms.

In legumes interacting with rhizobia, the formation of symbiotic organs involved in the acquisition of atmospheric nitrogen gas (N2) is dependent on the plant nitrogen (N) demand. We used Medicago truncatula plants cultivated in split-root systems to discriminate between responses to local and systemic N signaling. We evidenced a strong control of nodule formation by systemic N signaling but obtained no clear evidence of a local control by mineral nitrogen. Systemic signaling of the plant N demand controls numerous transcripts involved in root transcriptome reprogramming associated with early rhizobia interaction and nodule formation. SUPER NUMERIC NODULES (SUNN) has an important role in this control, but we found that major systemic N signaling responses remained active in the sunn mutant. Genes involved in the activation of nitrogen fixation are regulated by systemic N signaling in the mutant, explaining why its hypernodulation phenotype is not associated with higher nitrogen fixation of the whole plant. We show that the control of transcriptome reprogramming of nodule formation by systemic N signaling requires other pathway(s) that parallel the SUNN/CLE (CLAVATA3/EMBRYO SURROUNDING REGION-LIKE PEPTIDES) pathway.

Transcriptome of tambaqui Colossoma macropomum during gonad differentiation: Different molecular signals leading to sex identity.

Tambaqui (Colossoma macropomum) is the major native species in Brazilian aquaculture, and we have shown that females exhibit a higher growth compared to males, opening up the possibility for the production of all-female population. To date, there is no information on the sex determination and differentiation molecular mechanisms of tambaqui. In the present study, transcriptome sequencing of juvenile trunks was performed to understand the molecular network involved in the gonadal sex differentiation. The results showed that before differentiation, components of the Wnt/ß-catenin pathway, fox and fst genes imprint female sex development, whereas antagonistic pathways (gsk3b, wt1 and fgfr2), sox9 and genes for androgen synthesis indicate male differentiation. Hence, in undifferentiated tambaqui, the Wnt/ß-catenin exerts a role on sex differentiation, either upregulated in female-like individuals, or antagonized in male-like individuals.

Seed comparative genomics in three coffee species identify desiccation tolerance mechanisms in intermediate seeds.

In contrast to desiccation-tolerant 'orthodox' seeds, so-called 'intermediate' seeds cannot survive complete drying and are short-lived. All species of the genus Coffea produce intermediate seeds, but they show a considerable variability in seed desiccation tolerance (DT), which may help to decipher the molecular basis of seed DT in plants. We performed a comparative transcriptome analysis of developing seeds in three coffee species with contrasting desiccation tolerance. Seeds of all species shared a major transcriptional switch during late maturation that governs a general slow-down of metabolism. However, numerous key stress-related genes, including those coding for the late embryogenesis abundant protein EM6 and the osmosensitive calcium channel ERD4, were up-regulated during DT acquisition in the two species with high seed DT, C. arabica and C. eugenioides. By contrast, we detected up-regulation of numerous genes involved in the metabolism, transport, and perception of auxin in C. canephora seeds with low DT. Moreover, species with high DT showed a stronger down-regulation of the mitochondrial machinery dedicated to the tricarboxylic acid cycle and oxidative phosphorylation. Accordingly, respiration measurements during seed dehydration demonstrated that intermediate seeds with the highest DT are better prepared to cease respiration and avoid oxidative stresses.

Responses of mature symbiotic nodules to the whole-plant systemic nitrogen signaling.

In symbiotic root nodules of legumes, terminally differentiated rhizobia fix atmospheric N2 producing an NH4+ influx that is assimilated by the plant. The plant, in return, provides photosynthates that fuel the symbiotic nitrogen acquisition. Mechanisms responsible for the adjustment of the symbiotic capacity to the plant N demand remain poorly understood. We have investigated the role of systemic signaling of whole-plant N demand on the mature N2-fixing nodules of the model symbiotic association Medicago truncatula/Sinorhizobium using split-root systems. The whole-plant N-satiety signaling rapidly triggers reductions of both N2 fixation and allocation of sugars to the nodule. These responses are associated with the induction of nodule senescence and the activation of plant defenses against microbes, as well as variations in sugars transport and nodule metabolism. The whole-plant N-deficit responses mirror these changes: a rapid increase of sucrose allocation in response to N-deficit is associated with a stimulation of nodule functioning and development resulting in nodule expansion in the long term. Physiological, transcriptomic, and metabolomic data together provide evidence for strong integration of symbiotic nodules into whole-plant nitrogen demand by systemic signaling and suggest roles for sugar allocation and hormones in the signaling mechanisms.

Transcriptomic Analysis With the Progress of Symbiosis in 'Crack-Entry' Legume Arachis hypogaea Highlights Its Contrast With 'Infection Thread' Adapted Legumes.

In root-nodule symbiosis, rhizobial invasion and nodule organogenesis is host controlled. In most legumes, rhizobia enter through infection threads and nodule primordium in the cortex is induced from a distance. But in dalbergoid legumes like Arachis hypogaea, rhizobia directly invade cortical cells through epidermal cracks to generate the primordia. Herein, we report the transcriptional dynamics with the progress of symbiosis in A. hypogaea at 1 day postinfection (dpi) (invasion), 4 dpi (nodule primordia), 8 dpi (spread of infection in nodule-like structure), 12 dpi (immature nodules containing rod-shaped rhizobia), and 21 dpi (mature nodules with spherical symbiosomes). Expression of putative ortholog of symbiotic genes in 'crack entry' legume A. hypogaea was compared with infection thread-adapted model legumes. The contrasting features were i) higher expression of receptors like LYR3 and EPR3 as compared with canonical Nod factor receptors, ii) late induction of transcription factors like NIN and NSP2 and constitutive high expression of ERF1, EIN2, bHLH476, and iii) induction of divergent pathogenesis-responsive PR-1 genes. Additionally, symbiotic orthologs of SymCRK, ROP6, RR9, SEN1, and DNF2 were not detectable and microsynteny analysis indicated the absence of a RPG homolog in diploid parental genomes of A. hypogaea. The implications are discussed and a molecular framework that guides crack-entry symbiosis in A. hypogaea is proposed.

Plant population dynamics on oceanic islands during the Late Quaternary climate changes: genetic evidence from a tree species (Coffea mauritiana) in Reunion Island.

Past climatic fluctuations have played a major role in shaping the current plant biodiversity. Although harbouring an exceptional biota, oceanic islands have received little attention in studies on species demographic history and past vegetation patterns. We investigated the impact of past climatic changes on the effective population size of a tree (Coffea mauritiana) that is endemic to Reunion Island, located in the south-western Indian Ocean (SWIO). Demographic changes were inferred using summary statistics calculated from genomic data. Using ecological niche modelling and the current distribution of genetic diversity, the paleodistribution of the species was also assessed. A reduction in the effective population size of C. mauritiana during the last glaciation maximum was inferred. The distribution of the species was reduced on the western side of the island, due to low rainfall. It appeared that a major reduction in rainfall and a slight temperature decrease prevailed in the SWIO. Our findings indicated that analyses on the current patterns of intraspecific genetic variations can efficiently contribute to past climatic changes characterisation in remote islands. Identifying area with higher resilience in oceanic islands could provide guidance in forest management and conservation faced to the global climate change.

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.

Control of nerve cord formation by Engrailed and Gooseberry-Neuro: A multi-step, coordinated process.

One way to better understand the molecular mechanisms involved in the construction of a nervous system is to identify the downstream effectors of major regulatory proteins. We previously showed that Engrailed (EN) and Gooseberry-Neuro (GsbN) transcription factors act in partnership to drive the formation of posterior commissures in the central nervous system of Drosophila. In this report, we identified genes regulated by both EN and GsbN through chromatin immunoprecipitation ("ChIP on chip") and transcriptome experiments, combined to a genetic screen relied to the gene dose titration method. The genomic-scale approaches allowed us to define 175 potential targets of EN-GsbN regulation. We chose a subset of these genes to examine ventral nerve cord (VNC) defects and found that half of the mutated targets show clear VNC phenotypes when doubly heterozygous with en or gsbn mutations, or when homozygous. This strategy revealed new groups of genes never described for their implication in the construction of the nerve cord. Their identification suggests that, to construct the nerve cord, EN-GsbN may act at three levels, in: (i) sequential control of the attractive-repulsive signaling that ensures contralateral projection of the commissural axons, (ii) temporal control of the translation of some mRNAs, (iii) regulation of the capability of glial cells to act as commissural guideposts for developing axons. These results illustrate how an early, coordinated transcriptional control may orchestrate the various mechanisms involved in the formation of stereotyped neuronal networks. They also validate the overall strategy to identify genes that play crucial role in axonal pathfinding.

Transcriptome profiling of the honeybee parasite Varroa destructor provides new biological insights into the mite adult life cycle.

BACKGROUND: The parasite Varroa destructor represents a significant threat to honeybee colonies. Indeed, development of Varroa infestation within colonies, if left untreated, often leads to the death of the colony. Although its impact on bees has been extensively studied, less is known about its biology and the functional processes governing its adult life cycle and adaptation to its host. We therefore developed a full life cycle transcriptomic catalogue in adult Varroa females and included pairwise comparisons with males, artificially-reared and non-reproducing females (10 life cycle stages and conditions in total). RESULTS: Extensive remodeling of the Varroa transcriptome was observed, with an upregulation of energetic and chitin metabolic processes during the initial and final phases of the life cycle (e.g. phoretic and post-oviposition stages), whereas during reproductive stages in brood cells genes showing functions related to transcriptional regulation were overexpressed. Several neurotransmitter and neuropeptide receptors involved in behavioural regulation, as well as active compounds of salivary glands, were also expressed at a higher level outside the reproductive stages. No difference was detected between artificially-reared phoretic females and their counterparts in colonies, or between females who failed to reproduce and females who successfully reproduced, indicating that phoretic individuals can be reared outside host colonies without impacting their physiology and that mechanisms underlying reproductive failure occur before oogenesis. CONCLUSIONS: We discuss how these new findings reveal the remarkable adaptation of Varroa to its host biology and notably to the switch from living on adults to reproducing in sealed brood cells. By spanning the entire adult life cycle, our work captures the dynamic changes in the parasite gene expression and serves as a unique resource for deciphering Varroa biology and identifying new targets for mite control.

Insulators recruit histone methyltransferase dMes4 to regulate chromatin of flanking genes.

Chromosomal domains in Drosophila are marked by the insulator-binding proteins (IBPs) dCTCF/Beaf32 and cofactors that participate in regulating long-range interactions. Chromosomal borders are further enriched in specific histone modifications, yet the role of histone modifiers and nucleosome dynamics in this context remains largely unknown. Here, we show that IBP depletion impairs nucleosome dynamics specifically at the promoters and coding sequence of genes flanked by IBP binding sites. Biochemical purification identifies the H3K36 histone methyltransferase NSD/dMes-4 as a novel IBP cofactor, which specifically co-regulates the chromatin accessibility of hundreds of genes flanked by dCTCF/Beaf32. NSD/dMes-4 presets chromatin before the recruitment of transcriptional activators including DREF that triggers Set2/Hypb-dependent H3K36 trimethylation, nucleosome positioning, and RNA splicing. Our results unveil a model for how IBPs regulate nucleosome dynamics and gene expression through NSD/dMes-4, which may regulate H3K27me3 spreading. Our data uncover how IBPs dynamically regulate chromatin organization depending on distinct cofactors.

Gene coexpression network analysis of oil biosynthesis in an interspecific backcross of oil palm.

Global demand for vegetable oils is increasing at a dramatic rate, while our understanding of the regulation of oil biosynthesis in plants remains limited. To gain insights into the mechanisms that govern oil synthesis and fatty acid (FA) composition in the oil palm fruit, we used a multilevel approach combining gene coexpression analysis, quantification of allele-specific expression and joint multivariate analysis of transcriptomic and lipid data, in an interspecific backcross population between the African oil palm, Elaeis guineensis, and the American oil palm, Elaeis oleifera, which display contrasting oil contents and FA compositions. The gene coexpression network produced revealed tight transcriptional coordination of fatty acid synthesis (FAS) in the plastid with sugar sensing, plastidial glycolysis, transient starch storage and carbon recapture pathways. It also revealed a concerted regulation, along with FAS, of both the transfer of nascent FA to the endoplasmic reticulum, where triacylglycerol assembly occurs, and of the production of glycerol-3-phosphate, which provides the backbone of triacylglycerols. Plastid biogenesis and auxin transport were the two other biological processes most tightly connected to FAS in the network. In addition to WRINKLED1, a transcription factor (TF) known to activate FAS genes, two novel TFs, termed NF-YB-1 and ZFP-1, were found at the core of the FAS module. The saturated FA content of palm oil appeared to vary above all in relation to the level of transcripts of the gene coding for ß-ketoacyl-acyl carrier protein synthase II. Our findings should facilitate the development of breeding and engineering strategies in this and other oil crops.

Genomic characterization of Ensifer aridi, a proposed new species of nitrogen-fixing rhizobium recovered from Asian, African and American deserts.

BACKGROUND: Nitrogen fixing bacteria isolated from hot arid areas in Asia, Africa and America but from diverse leguminous plants have been recently identified as belonging to a possible new species of Ensifer (Sinorhizobium). In this study, 6 strains belonging to this new clade were compared with Ensifer species at the genome-wide level. Their capacities to utilize various carbon sources and to establish a symbiotic interaction with several leguminous plants were examined. RESULTS: Draft genomes of selected strains isolated from Morocco (Merzouga desert), Mexico (Baja California) as well as from India (Thar desert) were produced. Genome based species delineation tools demonstrated that they belong to a new species of Ensifer. Comparison of its core genome with those of E. meliloti, E. medicae and E. fredii enabled the identification of a species conserved gene set. Predicted functions of associated proteins and pathway reconstruction revealed notably the presence of transport systems for octopine/nopaline and inositol phosphates. Phenotypic characterization of this new desert rhizobium species showed that it was capable to utilize malonate, to grow at 48 °C or under high pH while NaCl tolerance levels were comparable to other Ensifer species. Analysis of accessory genomes and plasmid profiling demonstrated the presence of large plasmids that varied in size from strain to strain. As symbiotic functions were found in the accessory genomes, the differences in symbiotic interactions between strains may be well related to the difference in plasmid content that could explain the different legumes with which they can develop the symbiosis. CONCLUSIONS: The genomic analysis performed here confirms that the selected rhizobial strains isolated from desert regions in three continents belong to a new species. As until now only recovered from such harsh environment, we propose to name it Ensifer aridi. The presented genomic data offers a good basis to explore adaptations and functionalities that enable them to adapt to alkalinity, low water potential, salt and high temperature stresses. Finally, given the original phylogeographic distribution and the different hosts with which it can develop a beneficial symbiotic interaction, Ensifer aridi may provide new biotechnological opportunities for degraded land restoration initiatives in the future.

Do changes in gene expression contribute to sexual isolation and reinforcement in the house mouse?

Expression divergence, rather than sequence divergence, has been shown to be important in speciation, particularly in the early stages of divergence of traits involved in reproductive isolation. In the two European subspecies of house mice, Mus musculus musculus and Mus musculus domesticus, earlier studies have demonstrated olfactory-based assortative mate preference in populations close to their hybrid zone. It has been suggested that this behaviour evolved following the recent secondary contact between the two taxa (~3,000 years ago) in response to selection against hybridization. To test for a role of changes in gene expression in the observed behavioural shift, we conducted a RNA sequencing experiment on mouse vomeronasal organs. Key candidate genes for pheromone-based subspecies recognition, the vomeronasal receptors, are expressed in these organs. Overall patterns of gene expression varied significantly between samples from the two subspecies, with a large number of differentially expressed genes between the two taxa. In contrast, only ~200 genes were found repeatedly differentially expressed between populations within M. m. musculus that did or did not display assortative mate preferences (close to or more distant from the hybrid zone, respectively), with an overrepresentation of genes belonging to vomeronasal receptor family 2. These receptors are known to play a key role in recognition of chemical cues that handle information about genetic identity. Interestingly, four of five of these differentially expressed receptors belong to the same phylogenetic cluster, suggesting specialization of a group of closely related receptors in the recognition of odorant signals that may allow subspecies recognition and assortative mating.

Paraburkholderia piptadeniae sp. nov. and Paraburkholderia ribeironis sp. nov., two root-nodulating symbiotic species of Piptadenia gonoacantha in Brazil.

During a survey of root-nodulating symbionts of Mimosoid species in the south-east region of Brazil, eight Paraburkholderia isolates were obtained from nodules of the legume species Piptadenia gonoacantha, either from the field or following a soil trapping method with the same plant host. 16S rRNA gene as well as recA and gyrB phylogenetic markers placed these strains in two new clades within the genus Burkholderia sensu lato. DNA-DNA hybridization values and analyses of average nucleotide identities of the whole genome sequence of selected strains in each clade (STM 7183 and STM 7296) showed that the two clades represented novel species of the genus Paraburkholderia. All eight isolates were further characterized using DNA base content determination, chemotaxonomic and biochemical profiling and symbiotic properties, which allowed to distinguish the novel species from known diazotrophic species of the genus Paraburkholderia. Based on genomic and phenotypic data, the names Paraburkholderia piptadeniae sp. nov. with type strain STM 7183T (=DSM 101189T=LMG 29163T) and Paraburkholderia ribeironis sp. nov. with type strain STM 7296T (=DSM 101188T=LMG 29351T) are proposed.

Gene evolution and gene expression after whole genome duplication in fish: the PhyloFish database.

With more than 30,000 species, ray-finned fish represent approximately half of vertebrates. The evolution of ray-finned fish was impacted by several whole genome duplication (WGD) events including a teleost-specific WGD event (TGD) that occurred at the root of the teleost lineage about 350 million years ago (Mya) and more recent WGD events in salmonids, carps, suckers and others. In plants and animals, WGD events are associated with adaptive radiations and evolutionary innovations. WGD-spurred innovation may be especially relevant in the case of teleost fish, which colonized a wide diversity of habitats on earth, including many extreme environments. Fish biodiversity, the use of fish models for human medicine and ecological studies, and the importance of fish in human nutrition, fuel an important need for the characterization of gene expression repertoires and corresponding evolutionary histories of ray-finned fish genes. To this aim, we performed transcriptome analyses and developed the PhyloFish database to provide (i) de novo assembled gene repertoires in 23 different ray-finned fish species including two holosteans (i.e. a group that diverged from teleosts before TGD) and 21 teleosts (including six salmonids), and (ii) gene expression levels in ten different tissues and organs (and embryos for many) in the same species. This resource was generated using a common deep RNA sequencing protocol to obtain the most exhaustive gene repertoire possible in each species that allows between-species comparisons to study the evolution of gene expression in different lineages. The PhyloFish database described here can be accessed and searched using RNAbrowse, a simple and efficient solution to give access to RNA-seq de novo assembled transcripts.