OROSOMUCOID PROTEIN 1 regulation of sphingolipid synthesis is required for nodulation in Aeschynomene evenia.

Legumes establish symbiotic interactions with nitrogen-fixing rhizobia that are accommodated in root-derived organs known as nodules. Rhizobial recognition triggers a plant symbiotic signaling pathway that activates 2 coordinated processes: infection and nodule organogenesis. How these processes are orchestrated in legume species utilizing intercellular infection and lateral root base nodulation remains elusive. Here, we show that Aeschynomene evenia OROSOMUCOID PROTEIN 1 (AeORM1), a key regulator of sphingolipid biosynthesis, is required for nodule formation. Using A. evenia orm1 mutants, we demonstrate that alterations in AeORM1 function trigger numerous early aborted nodules, defense-like reactions, and shorter lateral roots. Accordingly, AeORM1 is expressed during lateral root initiation and elongation, including at lateral root bases where nodule primordium form in the presence of symbiotic bradyrhizobia. Sphingolipidomics revealed that mutations in AeORM1 lead to sphingolipid overaccumulation in roots relative to the wild type, particularly for very long-chain fatty acid-containing ceramides. Taken together, our findings reveal that AeORM1-regulated sphingolipid homeostasis is essential for rhizobial infection and nodule organogenesis, as well as for lateral root development in A. evenia.

Multi-genome comparisons reveal gain-and-loss evolution of anti-Mullerian hormone receptor type 2 as a candidate master sex-determining gene in Percidae.

BACKGROUND: The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. RESULTS: We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary. CONCLUSIONS: Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

The genome of the ant Tetramorium bicarinatum reveals a tandem organization of venom peptides genes allowing the prediction of their regulatory and evolutionary profiles.

BACKGROUND: Venoms have evolved independently over a hundred times in the animal kingdom to deter predators and/or subdue prey. Venoms are cocktails of various secreted toxins, whose origin and diversification provide an appealing system for evolutionary researchers. Previous studies of the ant venom of Tetramorium bicarinatum revealed several Myrmicitoxin (MYRTX) peptides that gathered into seven precursor families suggesting different evolutionary origins. Analysis of the T. bicarinatum genome enabling further genomic approaches was necessary to understand the processes underlying the evolution of these myrmicitoxins. RESULTS: Here, we sequenced the genome of Tetramorium bicarinatum and reported the organisation of 44 venom peptide genes (vpg). Of the eleven chromosomes that make up the genome of T. bicarinatum, four carry the vpg which are organized in tandem repeats. This organisation together with the ML evolutionary analysis of vpg sequences, is consistent with evolution by local duplication of ancestral genes for each precursor family. The structure of the vpg into two or three exons is conserved after duplication events while the promoter regions are the least conserved parts of the vpg even for genes with highly identical sequences. This suggests that enhancer sequences were not involved in duplication events, but were recruited from surrounding regions. Expression level analysis revealed that most vpg are highly expressed in venom glands, although one gene or group of genes is much more highly expressed in each family. Finally, the examination of the genomic data revealed that several genes encoding transcription factors (TFs) are highly expressed in the venom glands. The search for binding sites (BS) of these TFs in the vpg promoters revealed hot spots of GATA sites in several vpg families. CONCLUSION: In this pioneering investigation on ant venom genes, we provide a high-quality assembly genome and the annotation of venom peptide genes that we think can fosters further genomic research to understand the evolutionary history of ant venom biochemistry.

Soluble CD14 produced by bovine mammary epithelial cells modulates their response to full length LPS.

Bovine mastitis remains a major disease in cattle world-wide. In the mammary gland, mammary epithelial cells (MEC) are sentinels equipped with receptors allowing them to detect and respond to the invasion by bacterial pathogens, in particular Escherichia coli. Lipopolysaccharide (LPS) is the major E. coli motif recognized by MEC through its interaction with the TLR4 receptor and the CD14 co-receptor. Previous studies have highlighted the role of soluble CD14 (sCD14) in the efficient recognition of LPS molecules possessing a full-length O-antigen (LPSS). We demonstrate here that MEC are able to secrete CD14 and are likely to contribute to the presence of sCD14 in milk. We then investigated how sCD14 modulates and is required for the response of MEC to LPSS. This study highlights the key role of sCD14 for the full activation of the Myd88-independent pathway by LPSS. We also identified several lncRNA that are activated in MEC in response to LPS, including one lncRNA showing homologies with the mir-99a-let-7c gene (MIR99AHG). Altogether, our results show that a full response to LPS by mammary epithelial cells requires sCD14 and provide detailed information on how milk sCD14 can contribute to an efficient recognition of LPS from coliform pathogens.

The black honey bee genome: insights on specific structural elements and a first step towards pangenomes.

BACKGROUND: The honey bee reference genome, HAv3.1, was produced from a commercial line sample that was thought to have a largely dominant Apis mellifera ligustica genetic background. Apis mellifera mellifera, often referred to as the black bee, has a separate evolutionary history and is the original type in western and northern Europe. Growing interest in this subspecies for conservation and non-professional apicultural practices, together with the necessity of deciphering genome backgrounds in hybrids, triggered the necessity for a specific genome assembly. Moreover, having several high-quality genomes is becoming key for taking structural variations into account in pangenome analyses. RESULTS: Pacific Bioscience technology long reads were produced from a single haploid black bee drone. Scaffolding contigs into chromosomes was done using a high-density genetic map. This allowed for re-estimation of the recombination rate, which was over-estimated in some previous studies due to mis-assemblies, which resulted in spurious inversions in the older reference genomes. The sequence continuity obtained was very high and the only limit towards continuous chromosome-wide sequences seemed to be due to tandem repeat arrays that were usually longer than 10 kb and that belonged to two main families, the 371 and 91 bp repeats, causing problems in the assembly process due to high internal sequence similarity. Our assembly was used together with the reference genome to genotype two structural variants by a pangenome graph approach with Graphtyper2. Genotypes obtained were either correct or missing, when compared to an approach based on sequencing depth analysis, and genotyping rates were 89 and 76% for the two variants. CONCLUSIONS: Our new assembly for the Apis mellifera mellifera honey bee subspecies demonstrates the utility of multiple high-quality genomes for the genotyping of structural variants, with a test case on two insertions and deletions. It will therefore be an invaluable resource for future studies, for instance by including structural variants in GWAS. Having used a single haploid drone for sequencing allowed a refined analysis of very large tandem repeat arrays, raising the question of their function in the genome. High quality genome assemblies for multiple subspecies such as presented here, are crucial for emerging projects using pangenomes.

A pangenome graph reference of 30 chicken genomes allows genotyping of large and complex structural variants.

BACKGROUND: The red junglefowl, the wild outgroup of domestic chickens, has historically served as a reference for genomic studies of domestic chickens. These studies have provided insight into the etiology of traits of commercial importance. However, the use of a single reference genome does not capture diversity present among modern breeds, many of which have accumulated molecular changes due to drift and selection. While reference-based resequencing is well-suited to cataloging simple variants such as single-nucleotide changes and short insertions and deletions, it is mostly inadequate to discover more complex structural variation in the genome. METHODS: We present a pangenome for the domestic chicken consisting of thirty assemblies of chickens from different breeds and research lines. RESULTS: We demonstrate how this pangenome can be used to catalog structural variants present in modern breeds and untangle complex nested variation. We show that alignment of short reads from 100 diverse wild and domestic chickens to this pangenome reduces reference bias by 38%, which affects downstream genotyping results. This approach also allows for the accurate genotyping of a large and complex pair of structural variants at the K feathering locus using short reads, which would not be possible using a linear reference. CONCLUSIONS: We expect that this new paradigm of genomic reference will allow better pinpointing of exact mutations responsible for specific phenotypes, which will in turn be necessary for breeding chickens that meet new sustainability criteria and are resilient to quickly evolving pathogen threats.

A mutant-based analysis of the establishment of Nod-independent symbiosis in the legume Aeschynomene evenia.

Intensive research on nitrogen-fixing symbiosis in two model legumes has uncovered the molecular mechanisms, whereby rhizobial Nod factors activate a plant symbiotic signaling pathway that controls infection and nodule organogenesis. In contrast, the so-called Nod-independent symbiosis found between Aeschynomene evenia and photosynthetic bradyrhizobia, which does not involve Nod factor recognition nor infection thread formation, is less well known. To gain knowledge on how Nod-independent symbiosis is established, we conducted a phenotypic and molecular characterization of A. evenia lines carrying mutations in different nodulation genes. Besides investigating the effect of the mutations on rhizobial symbiosis, we examined their consequences on mycorrhizal symbiosis and in nonsymbiotic conditions. Analyzing allelic mutant series for AePOLLUX, Ca2+/calmodulin dependent kinase, AeCYCLOPS, nodulation signaling pathway 2 (AeNSP2), and nodule inception demonstrated that these genes intervene at several stages of intercellular infection and during bacterial accommodation. We provide evidence that AeNSP2 has an additional nitrogen-dependent regulatory function in the formation of axillary root hairs at lateral root bases, which are rhizobia-colonized infection sites. Our investigation of the recently discovered symbiotic actor cysteine-rich receptor-like kinase specified that it is not involved in mycorrhization; however, it is essential for both symbiotic signaling and early infection during nodulation. These findings provide important insights on the modus operandi of Nod-independent symbiosis and contribute to the general understanding of how rhizobial-legume symbioses are established by complementing the information acquired in model legumes.

Contrasting Gene Decay in Subterranean Vertebrates: Insights from Cavefishes and Fossorial Mammals.

Evolution sometimes proceeds by loss, especially when structures and genes become dispensable after an environmental shift relaxes functional constraints. Subterranean vertebrates are outstanding models to analyze this process, and gene decay can serve as a readout. We sought to understand some general principles on the extent and tempo of the decay of genes involved in vision, circadian clock, and pigmentation in cavefishes. The analysis of the genomes of two Cuban species belonging to the genus Lucifuga provided evidence for the largest loss of eye-specific genes and nonvisual opsin genes reported so far in cavefishes. Comparisons with a recently evolved cave population of Astyanax mexicanus and three species belonging to the Chinese tetraploid genus Sinocyclocheilus revealed the combined effects of the level of eye regression, time, and genome ploidy on eye-specific gene pseudogenization. The limited extent of gene decay in all these cavefishes and the very small number of loss-of-function mutations per pseudogene suggest that their eye degeneration may not be very ancient, ranging from early to late Pleistocene. This is in sharp contrast with the identification of several vision genes carrying many loss-of-function mutations in ancient fossorial mammals, further suggesting that blind fishes cannot thrive more than a few million years in cave ecosystems.

DNAModAnnot: a R toolbox for DNA modification filtering and annotation.

MOTIVATION: Long-read sequencing technologies can be employed to detect and map DNA modifications at the nucleotide resolution on a genome-wide scale. However, published software packages neglect the integration of genomic annotation and comprehensive filtering when analyzing patterns of modified bases detected using Pacific Biosciences (PacBio) or Oxford Nanopore Technologies (ONT) data. Here, we present DNA Modification Annotation (DNAModAnnot), a R package designed for the global analysis of DNA modification patterns using adapted filtering and visualization tools. RESULTS: We tested our package using PacBio sequencing data to analyze patterns of the 6-methyladenine (6mA) in the ciliate Paramecium tetraurelia, in which high 6mA amounts were previously reported. We found P. tetraurelia 6mA genome-wide distribution to be similar to other ciliates. We also performed 5-methylcytosine (5mC) analysis in human lymphoblastoid cells using ONT data and confirmed previously known patterns of 5mC. DNAModAnnot provides a toolbox for the genome-wide analysis of different DNA modifications using PacBio and ONT long-read sequencing data. AVAILABILITY AND IMPLEMENTATION: DNAModAnnot is distributed as a R package available via GitHub (https://github.com/AlexisHardy/DNAModAnnot). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of the master sex determining gene in Northern pike (Esox lucius) reveals restricted sex chromosome differentiation.

Teleost fishes, thanks to their rapid evolution of sex determination mechanisms, provide remarkable opportunities to study the formation of sex chromosomes and the mechanisms driving the birth of new master sex determining (MSD) genes. However, the evolutionary interplay between the sex chromosomes and the MSD genes they harbor is rather unexplored. We characterized a male-specific duplicate of the anti-Müllerian hormone (amh) as the MSD gene in Northern Pike (Esox lucius), using genomic and expression evidence as well as by loss-of-function and gain-of-function experiments. Using RAD-Sequencing from a family panel, we identified Linkage Group (LG) 24 as the sex chromosome and positioned the sex locus in its sub-telomeric region. Furthermore, we demonstrated that this MSD originated from an ancient duplication of the autosomal amh gene, which was subsequently translocated to LG24. Using sex-specific pooled genome sequencing and a new male genome sequence assembled using Nanopore long reads, we also characterized the differentiation of the X and Y chromosomes, revealing a small male-specific insertion containing the MSD gene and a limited region with reduced recombination. Our study reveals an unexpectedly low level of differentiation between a pair of sex chromosomes harboring an old MSD gene in a wild teleost fish population, and highlights both the pivotal role of genes from the amh pathway in sex determination, as well as the importance of gene duplication as a mechanism driving the turnover of sex chromosomes in this clade.

Neurodevelopment vs. the immune system: Complementary contributions of maternally-inherited gene transcripts and proteins to successful embryonic development in fish.

The aim of this study was to investigate the respective contribution of maternally-inherited mRNAs and proteins to egg molecular cargo and to its developmental competence in fish using pikeperch as a model. Our study provides novel insights into the understanding of type-specific roles of maternally-inherited molecules in fish. Here we show, for the first time, that transcripts and proteins have distinct, yet complementary, functions in the egg of teleost fish. Maternally-inherited mRNAs would shape embryo neurodevelopment, while maternally-inherited proteins would rather be responsible for protecting the embryo against pathogens. Additionally, we observed that processes directly preceding ovulation may considerably affect the reproductive success by modifying expression level of genes crucial for proper embryonic development, being novel fish egg quality markers (e.g., smarca4 or h3f3a). These results are of major importance for understanding the influence of external factors on reproductive fitness in both captive and wild-type fish species.

Transcriptome Profiling of the Pacific Oyster Crassostrea gigas Visceral Ganglia over a Reproduction Cycle Identifies Novel Regulatory Peptides.

The neuropeptides involved in the regulation of reproduction in the Pacific oyster (Crassostrea gigas) are quite diverse. To investigate this diversity, a transcriptomic survey of the visceral ganglia (VG) was carried out over an annual reproductive cycle. RNA-seq data from 26 samples corresponding to VG at different stages of reproduction were de novo assembled to generate a specific reference transcriptome of the oyster nervous system and used to identify differentially expressed transcripts. Transcriptome mining led to the identification of novel neuropeptide precursors (NPPs) related to the bilaterian Eclosion Hormone (EH), crustacean female sex hormone/Interleukin 17, Nesfatin, neuroparsin/IGFBP, prokineticins, and urotensin I; to the protostome GNQQN, pleurin, prohormones 3 and 4, prothoracotropic hormones (PTTH), and QSamide/PXXXamide; to the lophotrochozoan CCWamide, CLCCY, HFAamide, and LXRX; and to the mollusk-specific NPPs CCCGS, clionin, FYFY, GNamide, GRWRN, GSWN, GWE, IWMPxxGYxx, LXRYamide, RTLFamide, SLRFamide, and WGAGamide. Among the complete repertoire of NPPs, no sex-biased expression was observed. However, 25 NPPs displayed reproduction stage-specific expression, supporting their involvement in the control of gametogenesis or associated metabolisms.

Venom Peptide Repertoire of the European Myrmicine Ant Manica rubida: Identification of Insecticidal Toxins.

Using an integrated transcriptomic and proteomic approach, we characterized the venom peptidome of the European red ant, Manica rubida. We identified 13 "myrmicitoxins" that share sequence similarities with previously identified ant venom peptides, one of them being identified as an EGF-like toxin likely resulting from a threonine residue modified by O-fucosylation. Furthermore, we conducted insecticidal assays of reversed-phase HPLC venom fractions on the blowfly Lucilia caesar, permitting us to identify six myrmicitoxins (i.e., U3-, U10-, U13-, U20-MYRTX-Mri1a, U10-MYRTX-Mri1b, and U10-MYRTX-Mri1c) with an insecticidal activity. Chemically synthesized U10-MYRTX-Mri1a, -Mri1b, -Mri1c, and U20-MYRTX-Mri1a irreversibly paralyzed blowflies at the highest doses tested (30-125 nmol·g-1). U13-MYRTX-Mri1a, the most potent neurotoxic peptide at 1 h, had reversible effects after 24 h (150 nmol·g-1). Finally, U3-MYRTX-Mri1a has no insecticidal activity, even at up to 55 nmol·g-1. Thus, M. rubida employs a paralytic venom rich in linear insecticidal peptides, which likely act by disrupting cell membranes.

Molecular signatures of muscle growth and composition deciphered by the meta-analysis of age-related public transcriptomics data.

The lean-to-fat ratio is a major issue in the beef meat industry from both carcass and meat production perspectives. This industrial perspective has motivated meat physiologists to use transcriptomics technologies to decipher mechanisms behind fat deposition within muscle during the time course of muscle growth. However, synthetic biological information from this volume of data remains to be produced to identify mechanisms found in various breeds and rearing practices. We conducted a meta-analysis on 10 transcriptomic data sets stored in public databases, from the longissimus thoracis of five different bovine breeds divergent by age. We updated gene identifiers on the last version of the bovine genome (UCD1.2), and the 715 genes common to the 10 studies were subjected to the meta-analysis. Of the 238 genes differentially expressed (DEG), we identified a transcriptional signature of the dynamic regulation of glycolytic and oxidative metabolisms that agrees with a known shift between those two pathways from the animal puberty. We proposed some master genes of the myogenesis, namely MYOG and MAPK14, as probable regulators of the glycolytic and oxidative metabolisms. We also identified overexpressed genes related to lipid metabolism (APOE, LDLR, MXRA8, and HSP90AA1) that may contribute to the expected enhanced marbling as age increases. Lastly, we proposed a transcriptional signature related to the induction (YBX1) or repression (MAPK14, YWAH, ERBB2) of the commitment of myogenic progenitors into the adipogenic lineage. The relationships between the abundance of the identified mRNA and marbling values remain to be analyzed in a marbling biomarkers discovery perspectives.

Sex chromosome and sex locus characterization in goldfish, Carassius auratus (Linnaeus, 1758).

BACKGROUND: Goldfish is an important model for various areas of research, including neural development and behavior and a species of significant importance in aquaculture, especially as an ornamental species. It has a male heterogametic (XX/XY) sex determination system that relies on both genetic and environmental factors, with high temperatures being able to produce female-to-male sex reversal. Little, however, is currently known on the molecular basis of genetic sex determination in this important cyprinid model. Here we used sequencing approaches to better characterize sex determination and sex-chromosomes in an experimental strain of goldfish. RESULTS: Our results confirmed that sex determination in goldfish is a mix of environmental and genetic factors and that its sex determination system is male heterogametic (XX/XY). Using reduced representation (RAD-seq) and whole genome (pool-seq) approaches, we characterized sex-linked polymorphisms and developed male specific genetic markers. These male specific markers were used to distinguish sex-reversed XX neomales from XY males and to demonstrate that XX female-to-male sex reversal could even occur at a relatively low rearing temperature (18 °C), for which sex reversal has been previously shown to be close to zero. We also characterized a relatively large non-recombining region (~ 11.7 Mb) on goldfish linkage group 22 (LG22) that contained a high-density of male-biased genetic polymorphisms. This large LG22 region harbors 373 genes, including a single candidate as a potential master sex gene, i.e., the anti-Mullerian hormone gene (amh). However, no sex-linked polymorphisms were detected in the coding DNA sequence of the goldfish amh gene. CONCLUSIONS: These results show that our goldfish strain has a relatively large sex locus on LG22, which is likely the Y chromosome of this experimental population. The presence of a few XX males even at low temperature also suggests that other environmental factors in addition to temperature could trigger female-to-male sex reversal. Finally, we also developed sex-linked genetic markers, which will be important tools for future research on sex determination in our experimental goldfish population. However, additional work would be needed to explore whether this sex locus is conserved in other populations of goldfish.

Domestication modulates the expression of genes involved in neurogenesis in high-quality eggs of Sander lucioperca.

Pikeperch, Sander lucioperca, is a species of high interest to the aquaculture. The expansion of its production can only be achieved by furthering domestication level. However, the mechanisms driving the domestication process in finfishes are poorly understood. Transcriptome profiling of eggs was found to be a useful tool allowing understanding of the domestication process in teleosts. In this study, using next-generation sequencing, the first pikeperch transcriptome has been generated as well as pikeperch-specific microarray comprising 35,343 unique probes. Next, we performed transcriptome profiling of eggs obtained from wild and domesticated populations. We found 710 differentially expressed genes that were linked mostly to nervous system development. These results provide new insights into processes that are directly involved in the domestication of finfishes. It can be suggested that all the identified processes were predetermined by the maternally derived set of genes contained in the unfertilized eggs. This allows us to suggest that fish behavior, along with many other processes, can be predetermined at the cellular level and may have significant implications on the adaptation of cultured fish to the natural environment. This also allows to suggest that fish behavior should be considered as a very important pikeperch aquaculture selection trait.

Meta-omics Provides Insights into the Impact of Hydrocarbon Contamination on Microbial Mat Functioning.

Photosynthetic microbial mats are stable, self-supported communities. Due to their coastal localization, these mats are frequently exposed to hydrocarbon contamination and are able to grow on it. To decipher how this contamination disturbs the functioning of microbial mats, we compared two mats: a contaminated mat exposed to chronic petroleum contamination and a reference mat. The taxonomic and metabolic structures of the mats in spring and fall were determined using metagenomic and metatranscriptomic approaches. Extremely high contamination disturbed the seasonal variations of the mat. ABC transporters, two-component systems, and type IV secretion system-related genes were overabundant in the contaminated mats. Xenobiotic degradation metabolism was minor in the metagenomes of both mats, and only the expression of genes involved in polycyclic aromatic hydrocarbon degradation was higher in the contaminated mat. Interestingly, the expression rates of genes involved in hydrocarbon activation decreased during the 1-year study period, concomitant with the decrease in easily degradable hydrocarbons, suggesting a transient effect of hydrocarbon contamination. Alteromonadales and Oceanospirillales hydrocarbonoclastic bacteria appeared to be key in hydrocarbon remediation in the contaminated mat. Overall, the contaminated microbial mat was able to cope with hydrocarbon contamination and displayed an adaptive functioning that modified seasonal behaviour.

Attenuated sensing of SHH by Ptch1 underlies evolution of bovine limbs.

The large spectrum of limb morphologies reflects the wide evolutionary diversification of the basic pentadactyl pattern in tetrapods. In even-toed ungulates (artiodactyls, including cattle), limbs are adapted for running as a consequence of progressive reduction of their distal skeleton to symmetrical and elongated middle digits with hoofed phalanges. Here we analyse bovine embryos to establish that polarized gene expression is progressively lost during limb development in comparison to the mouse. Notably, the transcriptional upregulation of the Ptch1 gene, which encodes a Sonic hedgehog (SHH) receptor, is disrupted specifically in the bovine limb bud mesenchyme. This is due to evolutionary alteration of a Ptch1 cis-regulatory module, which no longer responds to graded SHH signalling during bovine handplate development. Our study provides a molecular explanation for the loss of digit asymmetry in bovine limb buds and suggests that modifications affecting the Ptch1 cis-regulatory landscape have contributed to evolutionary diversification of artiodactyl limbs.

Multi-species annotation of transcriptome and chromatin structure in domesticated animals.

BACKGROUND: Comparative genomics studies are central in identifying the coding and non-coding elements associated with complex traits, and the functional annotation of genomes is a critical step to decipher the genotype-to-phenotype relationships in livestock animals. As part of the Functional Annotation of Animal Genomes (FAANG) action, the FR-AgENCODE project aimed to create reference functional maps of domesticated animals by profiling the landscape of transcription (RNA-seq), chromatin accessibility (ATAC-seq) and conformation (Hi-C) in species representing ruminants (cattle, goat), monogastrics (pig) and birds (chicken), using three target samples related to metabolism (liver) and immunity (CD4+ and CD8+ T cells). RESULTS: RNA-seq assays considerably extended the available catalog of annotated transcripts and identified differentially expressed genes with unknown function, including new syntenic lncRNAs. ATAC-seq highlighted an enrichment for transcription factor binding sites in differentially accessible regions of the chromatin. Comparative analyses revealed a core set of conserved regulatory regions across species. Topologically associating domains (TADs) and epigenetic A/B compartments annotated from Hi-C data were consistent with RNA-seq and ATAC-seq data. Multi-species comparisons showed that conserved TAD boundaries had stronger insulation properties than species-specific ones and that the genomic distribution of orthologous genes in A/B compartments was significantly conserved across species. CONCLUSIONS: We report the first multi-species and multi-assay genome annotation results obtained by a FAANG project. Beyond the generation of reference annotations and the confirmation of previous findings on model animals, the integrative analysis of data from multiple assays and species sheds a new light on the multi-scale selective pressure shaping genome organization from birds to mammals. Overall, these results emphasize the value of FAANG for research on domesticated animals and reinforces the importance of future meta-analyses of the reference datasets being generated by this community on different species.

Arctic charr brain transcriptome strongly affected by summer seasonal growth but only subtly by feed deprivation.

BACKGROUND: The Arctic charr (Salvelinus alpinus) has a highly seasonal feeding cycle that comprises long periods of voluntary fasting and a short but intense feeding period during summer. Therefore, the charr represents an interesting species for studying appetite-regulating mechanisms in fish. RESULTS: In this study, we compared the brain transcriptomes of fed and feed deprived charr over a 4 weeks trial during their summer feeding season. Despite prominent differences in body condition between fed and feed deprived charr at the end of the trial, feed deprivation affected the brain transcriptome only slightly. In contrast, the transcriptome differed markedly over time in both fed and feed deprived charr, indicating strong shifts in basic cell metabolic processes possibly due to season, growth, temperature, or combinations thereof. The GO enrichment analysis revealed that many biological processes appeared to change in the same direction in both fed and feed deprived fish. In the feed deprived charr processes linked to oxygen transport and apoptosis were down- and up-regulated, respectively. Known genes encoding for appetite regulators did not respond to feed deprivation. Gene expression of Deiodinase 2 (DIO2), an enzyme implicated in the regulation of seasonal processes in mammals, was lower in response to season and feed deprivation. We further found a higher expression of VGF (non-acronymic) in the feed deprived than in the fed fish. This gene encodes for a neuropeptide associated with the control of energy metabolism in mammals, and has not been studied in relation to regulation of appetite and energy homeostasis in fish. CONCLUSIONS: In the Arctic charr, external and endogenous seasonal factors for example the increase in temperature and their circannual growth cycle, respectively, evoke much stronger responses in the brain than 4 weeks feed deprivation. The absence of a central hunger response in feed deprived charr give support for a strong resilience to the lack of food in this high Arctic species. DIO2 and VGF may play a role in the regulation of energy homeostasis and need to be further studied in seasonal fish.