Depletion of slow-cycling PDGFRα+ADAM12+ mesenchymal cells promotes antitumor immunity by restricting macrophage efferocytosis.

The capacity to survive and thrive in conditions of limited resources and high inflammation is a major driver of tumor malignancy. Here we identified slow-cycling ADAM12+PDGFRα+ mesenchymal stromal cells (MSCs) induced at the tumor margins in mouse models of melanoma, pancreatic cancer and prostate cancer. Using inducible lineage tracing and transcriptomics, we demonstrated that metabolically altered ADAM12+ MSCs induced pathological angiogenesis and immunosuppression by promoting macrophage efferocytosis and polarization through overexpression of genes such as Gas6, Lgals3 and Csf1. Genetic depletion of ADAM12+ cells restored a functional tumor vasculature, reduced hypoxia and acidosis and normalized CAFs, inducing infiltration of effector T cells and growth inhibition of melanomas and pancreatic neuroendocrine cancer, in a process dependent on TGF-ß. In human cancer, ADAM12 stratifies patients with high levels of hypoxia and innate resistance mechanisms, as well as factors associated with a poor prognosis and drug resistance such as AXL. Altogether, our data show that depletion of tumor-induced slow-cycling PDGFRα+ MSCs through ADAM12 restores antitumor immunity.

sgDI-tector: defective interfering viral genome bioinformatics for detection of coronavirus subgenomic RNAs.

Coronavirus RNA-dependent RNA polymerases produce subgenomic RNAs (sgRNAs) that encode viral structural and accessory proteins. User-friendly bioinformatic tools to detect and quantify sgRNA production are urgently needed to study the growing number of next-generation sequencing (NGS) data of SARS-CoV-2. We introduced sgDI-tector to identify and quantify sgRNA in SARS-CoV-2 NGS data. sgDI-tector allowed detection of sgRNA without initial knowledge of the transcription-regulatory sequences. We produced NGS data and successfully detected the nested set of sgRNAs with the ranking M > ORF3a > N>ORF6 > ORF7a > ORF8 > S > E>ORF7b. We also compared the level of sgRNA production with other types of viral RNA products such as defective interfering viral genomes.

Experimental evolution links post-transcriptional regulation to Leishmania fitness gain.

The protozoan parasite Leishmania donovani causes fatal human visceral leishmaniasis in absence of treatment. Genome instability has been recognized as a driver in Leishmania fitness gain in response to environmental change or chemotherapy. How genome instability generates beneficial phenotypes despite potential deleterious gene dosage effects is unknown. Here we address this important open question applying experimental evolution and integrative systems approaches on parasites adapting to in vitro culture. Phenotypic analyses of parasites from early and late stages of culture adaptation revealed an important fitness tradeoff, with selection for accelerated growth in promastigote culture (fitness gain) impairing infectivity (fitness costs). Comparative genomics, transcriptomics and proteomics analyses revealed a complex regulatory network associated with parasite fitness gain, with genome instability causing highly reproducible, gene dosage-independent and -dependent changes. Reduction of flagellar transcripts and increase in coding and non-coding RNAs implicated in ribosomal biogenesis and protein translation were not correlated to dosage changes of the corresponding genes, revealing a gene dosage-independent, post-transcriptional mechanism of regulation. In contrast, abundance of gene products implicated in post-transcriptional regulation itself correlated to corresponding gene dosage changes. Thus, RNA abundance during parasite adaptation is controled by direct and indirect gene dosage changes. We correlated differential expression of small nucleolar RNAs (snoRNAs) with changes in rRNA modification, providing first evidence that Leishmania fitness gain in culture may be controlled by post-transcriptional and epitranscriptomic regulation. Our findings propose a novel model for Leishmania fitness gain in culture, where differential regulation of mRNA stability and the generation of modified ribosomes may potentially filter deleterious from beneficial gene dosage effects and provide proteomic robustness to genetically heterogenous, adapting parasite populations. This model challenges the current, genome-centric approach to Leishmania epidemiology and identifies the Leishmania transcriptome and non-coding small RNome as potential novel sources for the discovery of biomarkers that may be associated with parasite phenotypic adaptation in clinical settings.

The CovR regulatory network drives the evolution of Group B Streptococcus virulence.

Virulence of the neonatal pathogen Group B Streptococcus is under the control of the master regulator CovR. Inactivation of CovR is associated with large-scale transcriptome remodeling and impairs almost every step of the interaction between the pathogen and the host. However, transcriptome analyses suggested a plasticity of the CovR signaling pathway in clinical isolates leading to phenotypic heterogeneity in the bacterial population. In this study, we characterized the CovR regulatory network in a strain representative of the CC-17 hypervirulent lineage responsible of the majority of neonatal meningitis. Transcriptome and genome-wide binding analysis reveal the architecture of the CovR network characterized by the direct repression of a large array of virulence-associated genes and the extent of co-regulation at specific loci. Comparative functional analysis of the signaling network links strain-specificities to the regulation of the pan-genome, including the two specific hypervirulent adhesins and horizontally acquired genes, to mutations in CovR-regulated promoters, and to variability in CovR activation by phosphorylation. This regulatory adaptation occurs at the level of genes, promoters, and of CovR itself, and allows to globally reshape the expression of virulence genes. Overall, our results reveal the direct, coordinated, and strain-specific regulation of virulence genes by the master regulator CovR and suggest that the intra-species evolution of the signaling network is as important as the expression of specific virulence factors in the emergence of clone associated with specific diseases.

IRF8 regulates efficacy of therapeutic anti-CD20 monoclonal antibodies.

Anti-CD20 monoclonal antibodies such as Rituximab, Ofatumumab, and Obinutuzumab are widely used to treat lymphomas and autoimmune diseases. They act by depleting B cells, mainly through Fc-dependent effectors functions. Some patients develop resistance to treatment but the underlying mechanisms are poorly understood. Here, we performed a genome-wide CRISPR/Cas9 screen to identify genes regulating the efficacy of anti-CD20 antibodies. We used as a model the killing of RAJI B cells by Rituximab through complement-dependent-cytotoxicity (CDC). As expected, the screen identified MS4A1, encoding CD20, the target of Rituximab. Among other identified genes, the role of Interferon Regulatory Factor 8 (IRF8) was validated in two B-cell lines. IRF8 KO also decreased the efficacy of antibody-dependent cellular cytotoxicity and phagocytosis (ADCC and ADCP) induced by anti-CD20 antibodies. We further show that IRF8 is necessary for efficient CD20 transcription. Levels of IRF8 and CD20 RNA or proteins correlated in normal B cells and in hundreds of malignant B cells. Therefore, IRF8 regulates CD20 expression and controls the depleting capacity of anti-CD20 antibodies. Our results bring novel insights into the pathways underlying resistance to CD20-targeting immunotherapies.

The oxidative stress response of pathogenic Leptospira is controlled by two peroxide stress regulators which putatively cooperate in controlling virulence.

Pathogenic Leptospira are the causative agents of leptospirosis, the most widespread zoonotic infectious disease. Leptospirosis is a potentially severe and life-threatening emerging disease with highest burden in sub-tropical areas and impoverished populations. Mechanisms allowing pathogenic Leptospira to survive inside a host and induce acute leptospirosis are not fully understood. The ability to resist deadly oxidants produced by the host during infection is pivotal for Leptospira virulence. We have previously shown that genes encoding defenses against oxidants in L. interrogans are repressed by PerRA (encoded by LIMLP_10155), a peroxide stress regulator of the Fur family. In this study, we describe the identification and characterization of another putative PerR-like regulator (LIMLP_05620) in L. interrogans. Protein sequence and phylogenetic analyses indicated that LIMLP_05620 displayed all the canonical PerR amino acid residues and is restricted to pathogenic Leptospira clades. We therefore named this PerR-like regulator PerRB. In L. interrogans, the PerRB regulon is distinct from that of PerRA. While a perRA mutant had a greater tolerance to peroxide, inactivating perRB led to a higher tolerance to superoxide, suggesting that these two regulators have a distinct function in the adaptation of L. interrogans to oxidative stress. The concomitant inactivation of perRA and perRB resulted in a higher tolerance to both peroxide and superoxide and, unlike the single mutants, a double perRAperRB mutant was avirulent. Interestingly, this correlated with major changes in gene and non-coding RNA expression. Notably, several virulence-associated genes (clpB, ligA/B, and lvrAB) were repressed. By obtaining a double mutant in a pathogenic Leptospira strain, our study has uncovered an interplay of two PerRs in the adaptation of Leptospira to oxidative stress with a putative role in virulence and pathogenicity, most likely through the transcriptional control of a complex regulatory network.

HP1γ binding pre-mRNA intronic repeats modulates RNA splicing decisions.

HP1 proteins are best known as markers of heterochromatin and gene silencing. Yet, they are also RNA-binding proteins and the HP1γ/CBX3 family member is present on transcribed genes together with RNA polymerase II, where it regulates co-transcriptional processes such as alternative splicing. To gain insight in the role of the RNA-binding activity of HP1γ in transcriptionally active chromatin, we have captured and analysed RNAs associated with this protein. We find that HP1γ is specifically targeted to hexameric RNA motifs and coincidentally transposable elements of the SINE family. As these elements are abundant in introns, while essentially absent from exons, the HP1γ RNA association tethers unspliced pre-mRNA to chromatin via the intronic regions and limits the usage of intronic cryptic splice sites. Thus, our data unveil novel determinants in the relationship between chromatin and co-transcriptional splicing.

Maternal age affects equine day 8 embryo gene expression both in trophoblast and inner cell mass.

BACKGROUND: Breeding a mare until she is not fertile or even until her death is common in equine industry but the fertility decreases as the mare age increases. Embryo loss due to reduced embryo quality is partly accountable for this observation. Here, the effect of mare's age on blastocysts' gene expression was explored. Day 8 post-ovulation embryos were collected from multiparous young (YM, 6-year-old, N = 5) and older (OM, > 10-year-old, N = 6) non-nursing Saddlebred mares, inseminated with the semen of one stallion. Pure or inner cell mass (ICM) enriched trophoblast, obtained by embryo bisection, were RNA sequenced. Deconvolution algorithm was used to discriminate gene expression in the ICM from that in the trophoblast. Differential expression was analyzed with embryo sex and diameter as cofactors. Functional annotation and classification of differentially expressed genes and gene set enrichment analysis were also performed. RESULTS: Maternal aging did not affect embryo recovery rate, embryo diameter nor total RNA quantity. In both compartments, the expression of genes involved in mitochondria and protein metabolism were disturbed by maternal age, although more genes were affected in the ICM. Mitosis, signaling and adhesion pathways and embryo development were decreased in the ICM of embryos from old mares. In trophoblast, ion movement pathways were affected. CONCLUSIONS: This is the first study showing that maternal age affects gene expression in the equine blastocyst, demonstrating significant effects as early as 10 years of age. These perturbations may affect further embryo development and contribute to decreased fertility due to aging.

Role of the major determinant of polar flagellation FlhG in the endoflagella-containing spirochete Leptospira.

Spirochetes can be distinguished from other bacteria by their spiral-shaped morphology and subpolar periplasmic flagella. This study focused on FlhF and FlhG, which control the spatial and numerical regulation of flagella in many exoflagellated bacteria, in the spirochete Leptospira. In contrast to flhF which seems to be essential in Leptospira, we demonstrated that flhG- mutants in both the saprophyte L. biflexa and the pathogen L. interrogans were less motile than the wild-type strains in gel-like environments but not hyperflagellated as reported previously in other bacteria. Cryo-electron tomography revealed that the distance between the flagellar basal body and the tip of the cell decreased significantly in the flhG- mutant in comparison to wild-type and complemented strains. Additionally, comparative transcriptome analyses of L. biflexa flhG- and wild-type strains showed that FlhG acts as a negative regulator of transcription of some flagellar genes. We found that the L. interrogans flhG- mutant was attenuated for virulence in the hamster model. Cross-species complementation also showed that flhG is not interchangeable between species. Our results indicate that FlhF and FlhG in Leptospira contribute to governing cell motility but our data support the hypothesis that FlhF and FlhG function differently in each bacterial species, including among spirochetes.

Transcytosis subversion by M cell-to-enterocyte spread promotes Shigella flexneri and Listeria monocytogenes intracellular bacterial dissemination.

Microfold (M) cell host-pathogen interaction studies would benefit from the visual analysis of dynamic cellular and microbial interplays. We adapted a human in vitro M cell model to physiological bacterial infections, expression of fluorescent localization reporters and long-term three-dimensional time-lapse microscopy. This approach allows following key steps of M cell infection dynamics at subcellular resolution, from the apical onset to basolateral epithelial dissemination. We focused on the intracellular pathogen Shigella flexneri, classically reported to transcytose through M cells to initiate bacillary dysentery in humans, while eliciting poorly protective immune responses. Our workflow was critical to reveal that S. flexneri develops a bimodal lifestyle within M cells leading to rapid transcytosis or delayed vacuolar rupture, followed by direct actin motility-based propagation to neighboring enterocytes. Moreover, we show that Listeria monocytogenes, another intracellular pathogen sharing a tropism for M cells, disseminates in a similar manner and evades M cell transcytosis completely. We established that actin-based M cell-to-enterocyte spread is the major dissemination pathway for both pathogens and avoids their exposure to basolateral compartments in our system. Our results challenge the notion that intracellular pathogens are readily transcytosed by M cells to inductive immune compartments in vivo, providing a potential mechanism for their ability to evade adaptive immunity.

The transcriptional response of pathogenic Leptospira to peroxide reveals new defenses against infection-related oxidative stress.

Pathogenic Leptospira spp. are the causative agents of the waterborne zoonotic disease leptospirosis. Leptospira are challenged by numerous adverse conditions, including deadly reactive oxygen species (ROS), when infecting their hosts. Withstanding ROS produced by the host innate immunity is an important strategy evolved by pathogenic Leptospira for persisting in and colonizing hosts. In L. interrogans, genes encoding defenses against ROS are repressed by the peroxide stress regulator, PerR. In this study, RNA sequencing was performed to characterize both the L. interrogans response to low and high concentrations of hydrogen peroxide and the PerR regulon. We showed that Leptospira solicit three main peroxidase machineries (catalase, cytochrome C peroxidase and peroxiredoxin) and heme to detoxify oxidants produced during peroxide stress. In addition, canonical molecular chaperones of the heat shock response and DNA repair proteins from the SOS response were required for Leptospira recovering from oxidative damage. Identification of the PerR regulon upon exposure to H2O2 allowed to define the contribution of this regulator in the oxidative stress response. This study has revealed a PerR-independent regulatory network involving other transcriptional regulators, two-component systems and sigma factors as well as non-coding RNAs that putatively orchestrate, in concert with PerR, the oxidative stress response. We have shown that PerR-regulated genes encoding a TonB-dependent transporter and a two-component system (VicKR) are involved in Leptospira tolerance to superoxide. This could represent the first defense mechanism against superoxide in L. interrogans, a bacterium lacking canonical superoxide dismutase. Our findings provide an insight into the mechanisms required by pathogenic Leptospira to overcome oxidative damage during infection-related conditions. This will participate in framing future hypothesis-driven studies to identify and decipher novel virulence mechanisms in this life-threatening pathogen.

Mosquito-bacteria interactions during larval development trigger metabolic changes with carry-over effects on adult fitness.

In animals with distinct life stages such as holometabolous insects, adult phenotypic variation is often shaped by the environment of immature stages, including their interactions with microbes colonizing larval habitats. Such carry-over effects were previously observed for several adult traits of the mosquito Aedes aegypti after larval exposure to different bacteria, but the mechanistic underpinnings are unknown. Here, we investigated the molecular changes triggered by gnotobiotic larval exposure to different bacteria in Ae. aegypti. We initially screened a panel of 16 bacterial isolates from natural mosquito breeding sites to determine their ability to influence adult life-history traits. We subsequently focused on four bacterial isolates (belonging to Flavobacterium, Lysobacter, Paenibacillus, and Enterobacteriaceae) with significant carry-over effects on adult survival and found that they were associated with distinct transcriptomic profiles throughout mosquito development. Moreover, we detected carry-over effects at the level of gene expression for the Flavobacterium and Paenibacillus isolates. The most prominent transcriptomic changes in gnotobiotic larvae reflected a profound remodelling of lipid metabolism, which translated into phenotypic differences in lipid storage and starvation resistance at the adult stage. Together, our findings indicate that larval exposure to environmental bacteria trigger substantial physiological changes that impact adult fitness, uncovering a possible mechanism underlying carry-over effects of mosquito-bacteria interactions during larval development.

4-Methylcytosine DNA modification is critical for global epigenetic regulation and virulence in the human pathogen Leptospira interrogans.

In bacteria, DNA methylation can be facilitated by 'orphan' DNA methyltransferases lacking cognate restriction endonucleases, but whether and how these enzymes control key cellular processes are poorly understood. The effects of a specific modification, 4-methylcytosine (4mC), are even less clear, as this epigenetic marker is unique to bacteria and archaea, whereas the bulk of epigenetic research is currently performed on eukaryotes. Here, we characterize a 4mC methyltransferase from the understudied pathogen Leptospira spp. Inactivating this enzyme resulted in complete abrogation of CTAG motif methylation, leading to genome-wide dysregulation of gene expression. Mutants exhibited growth defects, decreased adhesion to host cells, higher susceptibility to LPS-targeting antibiotics, and, importantly, were no longer virulent in an acute infection model. Further investigation resulted in the discovery of at least one gene, that of an ECF sigma factor, whose transcription was altered in the methylase mutant and, subsequently, by mutation of the CTAG motifs in the promoter of the gene. The genes that comprise the regulon of this sigma factor were, accordingly, dysregulated in the methylase mutant and in a strain overexpressing the sigma factor. Our results highlight the importance of 4mC in Leptospira physiology, and suggest the same of other understudied species.

Mycobacterium abscessus virulence traits unraveled by transcriptomic profiling in amoeba and macrophages.

Free-living amoebae are thought to represent an environmental niche in which amoeba-resistant bacteria may evolve towards pathogenicity. To get more insights into factors playing a role for adaptation to intracellular life, we characterized the transcriptomic activities of the emerging pathogen Mycobacterium abscessus in amoeba and murine macrophages (Mϕ) and compared them with the intra-amoebal transcriptome of the closely related, but less pathogenic Mycobacterium chelonae. Data on up-regulated genes in amoeba point to proteins that allow M. abscessus to resist environmental stress and induce defense mechanisms, as well as showing a switch from carbohydrate carbon sources to fatty acid metabolism. For eleven of the most upregulated genes in amoeba and/or Mϕ, we generated individual gene knock-out M. abscessus mutant strains, from which ten were found to be attenuated in amoeba and/or Mϕ in subsequence virulence analyses. Moreover, transfer of two of these genes into the genome of M. chelonae increased the intra-Mϕ survival of the recombinant strain. One knock-out mutant that had the gene encoding Eis N-acetyl transferase protein (MAB_4532c) deleted, was particularly strongly attenuated in Mϕ. Taken together, M. abscessus intra-amoeba and intra-Mϕ transcriptomes revealed the capacity of M. abscessus to adapt to an intracellular lifestyle, with amoeba largely contributing to the enhancement of M. abscessus intra-Mϕ survival.

Insights into amebiasis using a human 3D-intestinal model.

Entamoeba histolytica is the causative agent of amebiasis, an infectious disease targeting the intestine and the liver in humans. Two types of intestinal infection are caused by this parasite: silent infection, which occurs in the majority of cases, and invasive disease, which affects 10% of infected persons. To understand the intestinal pathogenic process, several in vitro models, such as cell cultures, human tissue explants or human intestine xenografts in mice, have been employed. Nevertheless, our knowledge on the early steps of amebic intestinal infection and the molecules involved during human-parasite interaction is scarce, in part due to limitations in the experimental settings. In the present work, we took advantage of tissue engineering approaches to build a three-dimensional (3D)-intestinal model that is able to replicate the general characteristics of the human colon. This system consists of an epithelial layer that develops tight and adherens junctions, a mucus layer and a lamina propria-like compartment made up of collagen containing macrophages and fibroblast. By means of microscopy imaging, omics assays and the evaluation of immune responses, we show a very dynamic interaction between E. histolytica and the 3D-intestinal model. Our data highlight the importance of several virulence markers occurring in patients or in experimental models, but they also demonstrate the involvement of under described molecules and regulatory factors in the amoebic invasive process.

Enzymatic synthesis of random sequences of RNA and RNA analogues by DNA polymerase theta mutants for the generation of aptamer libraries.

Nucleic acid aptamers, especially RNA, exhibit valuable advantages compared to protein therapeutics in terms of size, affinity and specificity. However, the synthesis of libraries of large random RNAs is still difficult and expensive. The engineering of polymerases able to directly generate these libraries has the potential to replace the chemical synthesis approach. Here, we start with a DNA polymerase that already displays a significant template-free nucleotidyltransferase activity, human DNA polymerase theta, and we mutate it based on the knowledge of its three-dimensional structure as well as previous mutational studies on members of the same polA family. One mutant exhibited a high tolerance towards ribonucleotides (NTPs) and displayed an efficient ribonucleotidyltransferase activity that resulted in the assembly of long RNA polymers. HPLC analysis and RNA sequencing of the products were used to quantify the incorporation of the four NTPs as a function of initial NTP concentrations and established the randomness of each generated nucleic acid sequence. The same mutant revealed a propensity to accept other modified nucleotides and to extend them in long fragments. Hence, this mutant can deliver random natural and modified RNA polymers libraries ready to use for SELEX, with custom lengths and balanced or unbalanced ratios.

Sequanix: a dynamic graphical interface for Snakemake workflows.

Summary: We designed a PyQt graphical user interface-Sequanix-aimed at democratizing the use of Snakemake pipelines in the NGS space and beyond. By default, Sequanix includes Sequana NGS pipelines (Snakemake format) (http://sequana.readthedocs.io), and is also capable of loading any external Snakemake pipeline. New users can easily, visually, edit configuration files of expert-validated pipelines and can interactively execute these production-ready workflows. Sequanix will be useful to both Snakemake developers in exposing their pipelines and to a wide audience of users. Availability and implementation: Source on http://github.com/sequana/sequana, bio-containers on http://bioconda.github.io and Singularity hub (http://singularity-hub.org). Contact: dimitri.desvillechabrol@pasteur.fr or thomas.cokelaer@pasteur.fr. Supplementary information: Supplementary data are available at Bioinformatics online.

GH16 and GH81 family ß-(1,3)-glucanases in Aspergillus fumigatus are essential for conidial cell wall morphogenesis.

The fungal cell wall is a rigid structure because of fibrillar and branched ß-(1,3)-glucan linked to chitin. Softening of the cell wall is an essential phenomenon during fungal morphogenesis, wherein rigid cell wall structures are cleaved by glycosylhydrolases. During the search for glycosylhydrolases acting on ß-(1,3)-glucan, we identified seven genes in the Aspergillus fumigatus genome coding for potential endo-ß-(1,3)-glucanase. ENG1 (previously characterized and named ENGL1, Mouyna et al., ), belongs to the Glycoside-Hydrolase 81 (GH81) family, while ENG2 to ENG7, to GH16 family. ENG1 and four GH16 genes (ENG2-5) were expressed in the resting conidia as well as during germination, suggesting an essential role during A. fumigatus morphogenesis. Here, we report the effect of sequential deletion of AfENG2-5 (GH16) followed by AfENG1 (GH81) deletion in the Δeng2,3,4,5 mutant. The Δeng1,2,3,4,5 mutant showed conidial defects, with linear chains of conidia unable to separate while the germination rate was not affected. These results show, for the first time in a filamentous fungus, that endo ß-(1,3)-glucanases are essential for proper conidial cell wall assembly and thus segregation of conidia during conidiation.

A missense variant in the coil1A domain of the keratin 25 gene is associated with the dominant curly hair coat trait (Crd) in horse.

BACKGROUND: Curly horses present a variety of curl phenotypes that are associated with various degrees of curliness of coat, mane, tail and ear hairs. Their origin is still a matter of debate and several genetic hypotheses have been formulated to explain the diversity in phenotype, including the combination of autosomal dominant and recessive alleles. Our purpose was to map the autosomal dominant curly hair locus and identify the causal variant using genome-wide association study (GWAS) and whole-genome sequencing approaches. RESULTS: A GWAS was performed using a Bayesian sparse linear mixed model, based on 51 curly and 19 straight-haired French and North American horses from 13 paternal families genotyped on the Illumina EquineSNP50 BeadChip. A single strong signal was observed on equine chromosome 11, in a region that encompasses the type I keratin gene cluster. This region was refined by haplotype analysis to a segment including 36 genes, among which are 10 keratin genes (KRT-10, -12, -20, -23, -24, -25, -26, -27, -28, -222). To comprehensively identify candidate causal variants within all these genes, whole-genome sequences were obtained for one heterozygous curly stallion and its straight-haired son. Among the four non-synonymous candidate variants identified and validated in the curly region, only variant g.21891160G>A in the KRT25 gene (KRT25:p.R89H) was in perfect agreement with haplotype status in the whole pedigree. Genetic association was then confirmed by genotyping a larger population consisting of 353 horses. However, five discordant curly horses were observed, which carried neither the variant nor the main haplotype associated with curliness. Sequencing of KRT25 for two discordant horses did not identify any other deleterious variant, which suggests locus rather than allelic heterogeneity for the curly phenotype. CONCLUSIONS: We identified the KRT25:p.R89H variant as responsible for the dominant curly trait, but a second dominant locus may also be involved in the shape of hairs within North American Curly horses.

RiboTools: a Galaxy toolbox for qualitative ribosome profiling analysis.

MOTIVATION: Ribosome profiling provides genome-wide information about translational regulation. However, there is currently no standard tool for the qualitative analysis of Ribo-seq data. We present here RiboTools, a Galaxy toolbox for the analysis of ribosome profiling (Ribo-seq) data. It can be used to detect translational ambiguities, stop codon readthrough events and codon occupancy. It provides a large number of plots for the visualisation of these events.