Search

Context: Active acromegaly is characterized by lipolysis-induced insulin resistance, which suggests adipose tissue (AT) as a primary driver of metabolic aberrations.Objective: To study the gene expression landscape in AT in patients with acromegaly before and after disease control in order to understand the changes and to identify disease-specific biomarkers.Methods: RNA sequencing was performed on paired subcutaneous adipose tissue (SAT) biopsies from six patients with acromegaly at time of diagnosis and after curative surgery. Clustering and pathway analyses were performed in order to identify disease activity-dependent genes. In a larger patient cohort (n = 23), the corresponding proteins were measured in serum by immunoassay. Correlations between growth hormone (GH), insulin-like growth factor I (IGF-I), visceral AT (VAT), SAT, total AT, and serum proteins were analyzed.Results: 743 genes were significantly differentially expressed (P-adjusted < .05) in SAT before and after disease control. The patients clustered according to disease activity. Pathways related to inflammation, cell adhesion and extracellular matrix, GH and insulin signaling, and fatty acid oxidation were differentially expressed.Serum levels of HTRA1, METRNL, S100A8/A9, and PDGFD significantly increased after disease control (P < .05). VAT correlated with HTRA1 (R = 0.73) and S100A8/A9 (R = 0.55) (P < .05 for both).Conclusion: AT in active acromegaly is associated with a gene expression profile of fibrosis and inflammation, which may corroborate the hyper-metabolic state and provide a means for identifying novel biomarkers. Context: Active acromegaly is characterized by lipolysis-induced insulin resistance, which suggests adipose tissue (AT) as a primary driver of metabolic aberrations. Objective: To study the gene expression landscape in AT in patients with acromegaly before and after disease control in order to understand the changes and to identify disease-specific biomarkers. Methods: RNA sequencing was performed on paired subcutaneous adipose tissue (SAT) biopsies from six patients with acromegaly at time of diagnosis and after curative surgery. Clustering and pathway analyses were performed in order to identify disease activity-dependent genes. In a larger patient cohort (n = 23), the corresponding proteins were measured in serum by immunoassay. Correlations between growth hormone (GH), insulin-like growth factor I (IGF-I), visceral AT (VAT), SAT, total AT, and serum proteins were analyzed. Results: 743 genes were significantly differentially expressed (P-adjusted < .05) in SAT before and after disease control. The patients clustered according to disease activity. Pathways related to inflammation, cell adhesion and extracellular matrix, GH and insulin signaling, and fatty acid oxidation were differentially expressed. Serum levels of HTRA1, METRNL, S100A8/A9, and PDGFD significantly increased after disease control (P < .05). VAT correlated with HTRA1 (R = 0.73) and S100A8/A9 (R = 0.55) (P < .05 for both). Conclusion: AT in active acromegaly is associated with a gene expression profile of fibrosis and inflammation, which may corroborate the hyper-metabolic state and provide a means for identifying novel biomarkers.

Background Novel commercial kits for whole genome library preparation for next-generation sequencing on Illumina platforms promise shorter workflows, lower inputs and cost savings. Time savings are achieved by employing enzymatic DNA fragmentation and by combining end-repair and tailing reactions. Fewer cleanup steps also allow greater DNA input flexibility (1 ng-1 μg), PCR-free options from 100 ng DNA, and lower price as compared to the well-established sonication and tagmentation-based DNA library preparation kits. Results We compared the performance of four enzymatic fragmentation-based DNA library preparation kits (from New England Biolabs, Roche, Swift Biosciences and Quantabio) to a tagmentation-based kit (Illumina) using low input DNA amounts (10 ng) and PCR-free reactions with 100 ng DNA. With four technical replicates of each input amount and kit, we compared the kits’ fragmentation sequence-bias as well as performance parameters such as sequence coverage and the clinically relevant detection of single nucleotide and indel variants. While all kits produced high quality sequence data and demonstrated similar performance, several enzymatic fragmentation methods produced library insert sizes which deviated from those intended. Libraries with longer insert lengths performed better in terms of coverage, SNV and indel detection. Lower performance of shorter-insert libraries could be explained by loss of sequence coverage to overlapping paired-end reads, exacerbated by the preferential sequencing of shorter fragments on Illumina sequencers. We also observed that libraries prepared with minimal or no PCR performed best with regard to indel detection. Conclusions The enzymatic fragmentation-based DNA library preparation kits from NEB, Roche, Swift and Quantabio are good alternatives to the tagmentation based Nextera DNA flex kit from Illumina, offering reproducible results using flexible DNA inputs, quick workflows and lower prices. Libraries with insert DNA fragments longer than the cumulative sum of both read lengths avoid read overlap, thus produce more informative data that leads to strongly improved genome coverage and consequently also increased sensitivity and precision of SNP and indel detection. In order to best utilize such enzymatic fragmentation reagents, researchers should be prepared to invest time to optimize fragmentation conditions for their particular samples.

High-throughput sequencing is increasingly favoured to assay the presence and abundance of micro RNAs (miRNAs) in biological samples, even from low RNA amounts, and a number of commercial vendors now offer kits that allow miRNA sequencing from sub-nanogram (ng) inputs. However, although biases introduced during library preparation have been documented, the relative performance of current reagent kits has not been investigated in detail. Here, six commercial kits capable of handling

Objective To investigate the heterogeneity of in vitro expanded chondrocytes used for autologous chondrocyte implantation. Methods Human articular chondrocytes were expanded in vitro for 14 days, sorted into 86 single cells using fluorescence-activated cell sorting and subjected to single-cell RNA sequencing. Principal component, Cross R2 hierarchical clustering, and differential gene expression analyses were used for data evaluation. Flow cytometry and single-cell RT-qPCR (reverse transcriptase quantitative polymerase chain reaction) was used to validate the results of the RNA sequencing data Polyclonal chondrocyte populations from the same donor were differentiated in vitro toward the osteogenic and adipogenic lineages. Results There was considerable variation in gene expression between individual cells, but we found no evidence for separate cell subpopulations based on principal component, hierarchical clustering, and differential gene expression analysis. Most of the cells expressed all the markers defining mesenchymal stem cells, and as polyclonal chondrocyte populations from the same donor were shown to differentiate into osteocytes and adipocytes in vitro, these cells formally qualify as mesenchymal stem cells. Conclusions In vitro expanded chondrocytes consist of one single population of cells with heterogeneity in gene expression between the cells. Dedifferentiated chondrocytes qualify as mesenchymal stem cells as they fulfill all the criteria suggested by the International Society for Cellular Therapy.

Objectives The emerging SARS-CoV-2 variants of concern (VoC), B.1.1.7, B.1.351 and P.1, with increased transmission and/or immune evasion, emphasise the need for broad and rapid variant monitoring. Our high-volume laboratory introduced a PCR variant assay (Variant PCR) in January 2021 based on the protocol by Vogels et al. Study design To assess whether Variant PCR could be used for rapid B.1.1.7, B.1.351 and P.1 screening, all positive SARS-CoV-2 airway samples were prospectively tested in parallel using both the Variant PCR and whole genome sequencing (WGS). Results In total 1,642 SARS-CoV-2 positive samples from individual patients were tested within a time span of 4 weeks. For all samples with valid results from both Variant PCR and WGS, no VoC was missed by Variant PCR (totalling 399 VoC detected). Conversely, all of the samples identified as “other lineages” (i.e., “non-VoC lineages”) by the Variant PCR, were confirmed by WGS. Conclusions The Variant PCR based on the protocol by Vogels et al., is an effective method for rapid screening for VoC, applicable for most diagnostic laboratories within a pandemic setting. WGS is still required to confirm the identity of certain variants and for continuous surveillance of emerging VoC.

Listeria monocytogenes is a food-borne bacterium that give rise to the potentially life-threatening disease listeriosis. Listeriosis has been mandatorily notifiable in Norway since 1991. All clinical L. monocytogenes isolates are sent to the Norwegian Institute of Public Health (NIPH) for typing. Since 2005 Multi-Locus Variable number tandem repeats Analysis (MLVA) has been used for typing but was recently replaced by whole genome sequencing using core genome Multi-Locus Sequence Typing (cgMLST). In the present study, L. monocytogenes isolates collected at salmon processing plants in Norway in 2007 (n = 12) and 2015 (n = 14) were first subject to MLVA. Twelve clinical L. monocytogenes isolates with matching MLVA profile and sampling time were selected from the strain collection at NIPH. Twenty-one isolates from the salmon processing plants and all clinical isolates (n = 12) were whole genome sequenced and compared using cgMLST and in silico detection of virulence genes. cgMLST revealed four pairs of environmental–human isolates with ≤10 allelic differences over 1708 genes, indicating that they may be assigned as clonal, with the implication that they are descended from the same recent ancestor. No relevant difference in carriage of virulence genes was found between environmental or human isolates. The present study shows that L. monocytogenes strains that genetically resemble contemporary isolates from human listeriosis circulate in Norwegian salmon slaughterhouses, and carry the same virulence genes.

This study investigated the effect of feeding Cyberlindnera jadinii yeast on the development of gastrointestinal function and health in piglets during the first two challenging weeks after weaning. Changes in gastrointestinal function were mainly attributed to weaning, and not to dietary treatment. The post-weaning (PW) transcriptome profiles differed between dietary treatments showing an overall higher number of differentially expressed genes (DEGs) in control piglets than in yeast-fed piglets. DEGs in jejunum and ileum were compared between sampling timepoints within each feeding group and divided into clusters with similar expression trends. Pathway enrichment analysis was run on each cluster to reveal PW physiological changes. Weaning induced downregulation of several immune functions in the control piglets, which was not as evident in the yeast fed piglets. The results indicate that feeding C. jadinii yeast can improve PW gut homeostasis and give more robust piglets.

Simple Summary Corticotroph pituitary adenomas present a spectrum of functionality regarding hormonal production, ranging from functioning to silent tumors. Moreover, they show different invasiveness and recurrent behavior profiles, the silent being considered an aggressive type of adenomas. Through analyses of global transcriptome and proteome, we show that both groups expressed genes and protein related to protein synthesis and vesicular transport, and present a distinct pattern of collagen/ extracellular matrix proteins. Endoplasmic reticulum protein processing is a key factor for hormone production in functioning corticotroph adenomas. Furthermore, a distinct cell adhesion profile in silent corticotroph adenomas may explain the aggressive behavior. Together, our findings shed light on the different repertoires of activated signaling pathways in corticotroph pituitary adenomas and may reveal new potential medical targets. Abstract Functioning (FCA) and silent corticotroph (SCA) pituitary adenomas act differently from a clinical perspective, despite both subtypes showing positive TBX19 (TPIT) and/or adrenocorticotropic hormone (ACTH) staining by immunohistochemistry. They are challenging to treat, the former due to functional ACTH production and consequently hypercortisolemia, and the latter due to invasive and recurrent behavior. Moreover, the molecular mechanisms behind their distinct behavior are not clear. We investigated global transcriptome and proteome changes in order to identify signaling pathways that can explain FCA and SCA differences (e.g., hormone production vs. aggressive growth). In the transcriptomic study, cluster analyses of differentially expressed genes revealed two distinct groups in accordance with clinical and histological classification. However, in the proteomic study, a greater degree of heterogeneity within the SCA group was found. Genes and proteins related to protein synthesis and vesicular transport were expressed by both adenoma groups, although different types and a distinct pattern of collagen/extracellular matrix proteins were presented by each group. Moreover, several genes related to endoplasmic reticulum protein processing were overexpressed in the FCA group. Together, our findings shed light on the different repertoires of activated signaling pathways in corticotroph adenomas, namely, the increased protein processing capacity of FCA and a specific pattern of adhesion molecules that may play a role in the aggressiveness of SCA.

Background Hydrogen peroxide (H2O2) is one of the delousing agents used to control sea lice infestations in salmonid aquaculture. However, some Lepeophtheirus salmonis populations have developed resistance towards H2O2. An increased gene expression and activity of catalase, an enzyme that breaks down H2O2, have been detected in resistant lice, being therefore introduced as a resistance marker in the salmon industry. In the present study the aim was to validate the use of catalase expression as a marker and to identify new candidate genes as additional markers to catalase, related to H2O2 resistance in L. salmonis. Methods A sensitive and an H2O2 resistant laboratory strain (P0 generation, not exposed to H2O2 for several years) were batch crossed to generate a cohort with a wide range of H2O2 sensitivities (F2 generation). F2 adult females were then exposed to H2O2 to separate sensitive and resistant individuals. Those F2 lice, the P0 lice and field-collected resistant lice (exposed to H2O2 in the field) were used in an RNA sequencing study. Results Catalase was upregulated in resistant lice exposed to H2O2 compared to sensitive lice. This was, however, not the case for unexposed resistant P0 lice. Several other genes were found differentially expressed between sensitive and resistant lice, but most of them seemed to be related to H2O2 exposure. However, five genes were consistently up- or downregulated in the resistant lice independent of exposure history. The upregulated genes were: one gene in the DNA polymerase family, one gene encoding a Nesprin-like protein and an unannotated gene encoding a small protein. The downregulated genes encoded endoplasmic reticulum resident protein 29 and an aquaporin (Glp1_v2). Conclusions Catalase expression seems to be induced by H2O2 exposure, since it was not upregulated in unexposed resistant lice. This may pose a challenge for its use as a resistance marker. The five new genes associated with resistance are put forward as complementary candidate genes. The most promising was Glp1_v2, an aquaglyceroporin that may serve as a passing channel for H2O2. Lower channel number can reduce the influx or distribution of H2O2 in the salmon louse, being directly involved in the resistance mechanism.

Background The thymus is a highly specialized organ of the immune system where T cell precursors develop and differentiate into self-tolerant CD4+ or CD8+ T cells. No studies to date have investigated how the human transcriptome profiles differ, between T cells still residing in the thymus and T cells in the periphery. Results We have performed high-throughput RNA sequencing to characterize the transcriptomes of primary single positive (SP) CD4+ and CD8+ T cells from infant thymic tissue, as well as primary CD4+ and CD8+ T cells from infant and adult peripheral blood, to enable the comparisons across tissues and ages. In addition, we have assessed the expression of candidate genes related to autoimmune diseases in thymic CD4+ and CD8+ T cells. The thymic T cells showed the largest number of uniquely expressed genes, suggesting a more diverse transcription in thymic T cells. Comparing T cells of thymic and blood origin, revealed more differentially expressed genes, than between infant and adult blood. Functional enrichment analysis revealed an over-representation of genes involved in cell cycle and replication in thymic T cells, whereas infant blood T cells were dominated by immune related terms. Comparing adult and infant blood T cells, the former was enriched for inflammatory response, cytokine production and biological adhesion, while upregulated genes in infant blood T cells were associated with cell cycle, cell death and gene expression. Conclusion This study provides valuable insight into the transcriptomes of the human primary SP T cells still residing within the thymus, and offers a unique comparison to primary blood derived T cells. Interestingly, the majority of autoimmune disease associated genes were expressed in one or more T cell subset, however ~ 11% of these were not expressed in frequently studied adult peripheral blood.

Listeriosis has been mandatory notifiable in Norway since 1991. Clinical isolatesof L. monocytogenes are routinely sent to the Norwegian Institute of Public Health (NIPH) for typing. Since 2005 Multi-Locus Variable number tandem repeats Analysis (MLVA) has been used for typing but was recently replaced by whole genome sequencing (WGS) using core genome Multi-Locus Sequence Typing (cgMLST). L. monocytogenes isolates collected at salmon processing plants in Norway in 2007 (n=12) and 2015 (n=9) where first subject to MLVA. Twelve clinical L. monocytogenes isolates with matching MLVA profile and sampling time were selected from the strain collection at NIPH. Twenty-one isolates from the salmon processing plants and all clinical isolates (n=12) were whole genome sequenced and compared using cgMLST, revealimg four pairs of environmental – human isolates with ≤ 10 allelic differences, indicating that they may be assigned as clonal, with the implication that they are descended from the same recent ancestor.

High-throughput sequencing has emerged as the favoured method to study microRNA (miRNA) expression, but biases introduced during library preparation have been reported. To assist researchers choose the most appropriate library preparation kit, we recently compared the performance of six commercially-available kits on synthetic miRNAs and human RNA, where library preparation was performed by the vendors. We hereby supplement this study with data from two further commonly used kits (NEBNext, NEXTflex) whose manufacturers initially declined to participate. As before, performance was assessed with respect to sensitivity, reliability, titration response and differential expression. Despite NEXTflex employing partially-randomised adapter sequences to minimise bias, we reaffirm that biases in miRNA abundance are kit-specific, complicating the comparison of miRNA datasets generated using different kits.

The IPD-MHC Database represents the official repository for non-human major histocompatibility complex (MHC) sequences, overseen and supported by the Comparative MHC Nomenclature Committee, providing access to curated MHC data and associated analysis tools. IPD-MHC gathers allelic MHC class I and class II sequences from classical and non-classical MHC loci from various non-human animals including pets, farmed and experimental model animals. So far, Atlantic salmon and rainbow trout are the only teleost fish species with MHC class I and class II sequences present. For the remaining teleost or ray-finned species, data on alleles originating from given classical locus is scarce hampering their inclusion in the database. However, a fast expansion of sequenced genomes opens for identification of classical loci where high-throughput sequencing (HTS) will enable typing of allelic variants in a variety of new teleost or ray-finned species. HTS also opens for large-scale studies of salmonid MHC diversity challenging the current database nomenclature and analysis tools. Here we establish an Illumina approach to identify allelic MHC diversity in Atlantic salmon, using animals from an endangered wild population, and alter the salmonid MHC nomenclature to accommodate the expected sequence expansions. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s00251-019-01143-8) contains supplementary material, which is available to authorized users.

ObjectiveReliable biomarkers associated with aggressiveness of non-functioning gonadotroph adenomas (GAs) are lacking. As the growth of tumor remnants is highly variable, molecular markers for growth potential prediction are necessary. We hypothesized that fast- and slow-growing GAs present different gene expression profiles and reliable biomarkers for tumor growth potential could be identified, focusing on the specific role of epithelial-mesenchymal transition (EMT).Design and methodsEight GAs selected for RNA sequencing were equally divided into fast- and slow-growing group by the tumor volume doubling time (TVDT) median (27.75 months). Data were analyzed by tophat2, cufflinks and cummeRbund pipeline. 40 genes were selected for RT-qPCR validation in 20 GAs based on significance, fold-change and pathway analyses. The effect of silencingMTDH(metadherin) andEMCN(endomucin) onin vitromigration of human adenoma cells was evaluated.Results350 genes were significantly differentially expressed (282 genes upregulated and 68 downregulated in the fast group,P-adjusted RP PCDH18, UNC5D, EMCN, MYO1B, GPM6Aand six EMT-related genes (SPAG9, SKIL, MTDH, HOOK1, CNOT6LandPRKACB).MTDH, but notEMCN, demonstrated involvement in cell migration and association with EMT markers.ConclusionsFast- and slow-growing GAs present different gene expression profiles, and genes related to EMT have higher expression in fast-growing tumors. In addition toMTDH, identified as an important contributor to aggressiveness, the other genes might represent markers for tumor growth potential and possible targets for drug therapy.

High throughput sequencing (HTS) is useful for many purposes as exemplified by the other topics included in this special issue. The purpose of this paper is to look into the unique challenges of using this technology in non-model organisms where resources such as genomes, functional genome annotations or genome complexity provide obstacles not met in model organisms. To describe these challenges, we narrow our scope to RNA sequencing used to study differential gene expression in response to pathogen challenge. As a demonstration species we chose Atlantic salmon, which has a sequenced genome with poor annotation and an added complexity due to many duplicated genes. We find that our RNA-seq analysis pipeline deciphers between duplicates despite high sequence identity. However, annotation issues provide problems in linking differentially expressed genes to pathways. Also, comparing results between approaches and species are complicated due to lack of standardized annotation.

The transcription factor c-Myb is involved in early differentiation and proliferation of haematopoietic cells, where it operates as a regulator of self-renewal and multi-lineage differentiation. Deregulated c-Myb plays critical roles in leukaemias and other human cancers. Due to its role as a master regulator, we hypothesized it might function as a pioneer transcription factor. Our approach to test this was to analyse a mutant of c-Myb, D152V, previously reported to cause haematopoietic defects in mice by an unknown mechanism. Our transcriptome data from K562 cells indicates that this mutation specifically affects c-Myb's ability to regulate genes involved in differentiation, causing failure in c-Myb's ability to block differentiation. Furthermore, we see a major effect of this mutation in assays where chromatin opening is involved. We show that each repeat in the minimal DNA-binding domain of c-Myb binds to histones and that D152V disrupts histone binding of the third repeat. ATAC-seq data indicates this mutation impairs the ability of c-Myb to cause chromatin opening at specific sites. Taken together, our findings support that c-Myb acts as a pioneer factor and show that D152V impairs this function. The D152V mutant is the first mutant of a transcription factor specifically destroying pioneer factor function.

Listeria monocytogenes, the causative agent of the serious food-borne disease listeriosis, can rapidly adapt to a wide-range of environmental stresses, including visible light. This study shows that exposure of the L. monocytogenes EGDe strain to low-intensity broad spectrum visible light inhibited bacterial growth and caused altered multicellular behavior during growth on semi-solid agar as compared to when the bacteria were grown in complete darkness. These light-dependent changes were observed regardless of the presence of the blue light receptor (Lmo0799) and the stressosome regulator Sigma B (SigB), which have been suggested to be important for the ability of L. monocytogenes to respond to blue light. A genome wide transcriptional analysis revealed that exposure of L. monocytogenes EGDe to broad-spectrum visible light caused altered expression of 2409 genes belonging to 18 metabolic pathways, as compared to bacteria grown in darkness. The light-dependent differentially expressed genes are involved in functions such as glycan metabolism, cell wall synthesis, chemotaxis, flagellar synthesis and resistance to oxidative stress. Exposure to light conferred reduced bacterial motility in semi-solid agar, which correlates well with the light-dependent reduction in transcript levels of flagellar and chemotaxis genes. Similar light-induced reduction in growth and motility was also observed in two different L. monocytogenes food isolates, suggesting that these responses are typical for L. monocytogenes. Together, the results show that even relatively small doses of broad-spectrum visible light causes genome-wide transcriptional changes, reduced growth and motility in L. monocytogenes. IMPORTANCE Despite major efforts to control L. monocytogenes, this pathogen remains a major problem to the food industry where it poses a continuous risk of food contamination. The ability of L. monocytogenes to sense and adapt to different stressors in the environments enables it to persist in many different niches, including food production facilities and in food products. The present study shows that exposure of L. monocytogenes to low-intensity broad-spectrum visible light, reduces its growth and motility as well as alters its multicellular behavior. Light exposure also caused genome-wide changes in transcript levels, affecting multiple metabolic pathways, which are likely to influence the bacterial physiology and lifestyle. In practical terms, the data presented in this study suggest that broad spectrum visible light is an important environmental variable to consider as a strategy to improve food safety by reducing L. monocytogenes contamination in food production environments.

The IPD-MHC Database represents the official repository for non-human major histocompatibility complex (MHC) sequences, overseen and supported by the Comparative MHC Nomenclature Committee, providing access to curated MHC data and associated analysis tools. IPD-MHC gathers allelic MHC class I and class II sequences from classical and non-classical MHC loci from various non-human animals including pets, farmed and experimental model animals. So far, Atlantic salmon and rainbow trout are the only teleost fish species with MHC class I and class II sequences present. For the remaining teleost or ray-finned species, data on alleles originating from given classical locus is scarce hampering their inclusion in the database. However, a fast expansion of sequenced genomes opens for identification of classical loci where high-throughput sequencing (HTS) will enable typing of allelic variants in a variety of new teleost or ray-finned species. HTS also opens for large-scale studies of salmonid MHC diversity challenging the current database nomenclature and analysis tools. Here we establish an Illumina approach to identify allelic MHC diversity in Atlantic salmon, using animals from an endangered wild population, and alter the salmonid MHC nomenclature to accommodate the expected sequence expansions.

OBJECTIVE: To investigate the heterogeneity of in vitro expanded chondrocytes used for autologous chondrocyte implantation. METHODS: Human articular chondrocytes were expanded in vitro for 14 days, sorted into 86 single cells using fluorescence-activated cell sorting and subjected to single-cell RNA sequencing. Principal component, Cross R(2) hierarchical clustering, and differential gene expression analyses were used for data evaluation. Flow cytometry and single-cell RT-qPCR (reverse transcriptase quantitative polymerase chain reaction) was used to validate the results of the RNA sequencing data Polyclonal chondrocyte populations from the same donor were differentiated in vitro toward the osteogenic and adipogenic lineages. RESULTS: There was considerable variation in gene expression between individual cells, but we found no evidence for separate cell subpopulations based on principal component, hierarchical clustering, and differential gene expression analysis. Most of the cells expressed all the markers defining mesenchymal stem cells, and as polyclonal chondrocyte populations from the same donor were shown to differentiate into osteocytes and adipocytes in vitro, these cells formally qualify as mesenchymal stem cells. CONCLUSIONS: In vitro expanded chondrocytes consist of one single population of cells with heterogeneity in gene expression between the cells. Dedifferentiated chondrocytes qualify as mesenchymal stem cells as they fulfill all the criteria suggested by the International Society for Cellular Therapy.

This study was performed to investigate the effects of dietary inclusion of 20% rapeseed meal (RSM) as an alternative to soybean meal (SBM) in a three-month feeding experiment with growing finishing pigs. Dietary alteration affected growth performance, several carcass traits and transcriptional responses in the skeletal muscle, but did not affect measured meat quality traits. In general, pigs fed the RSM test diet exhibited reduced growth performance compared to pigs on SBM control diet. Significant transcriptional changes in the skeletal muscle of growing pigs fed RSM diet were likely the consequence of an increased amount of fiber and higher polyunsaturated fatty acids, and presence of bioactive phytochemicals, such as glucosinolates. RNAseq pipeline using Tophat2-Cuffdiff identified 57 upregulated and 63 downregulated genes in RSM compared to SBM pigs. Significantly enriched among downregulated pathways was p53–mediated signalling involved in cellular proliferation, while activation of negative growth regulators (IER5, KLF10, BTG2, KLF11, RETREG1, PRUNE2) in RSM fed pigs provided further evidence for reduced proliferation and increased cellular death, in accordance with the observed reduction in performance traits. Upregulation of well-known metabolic controllers (PDK4, UCP3, ESRRG and ESRRB), involved in energy homeostasis (glucose and lipid metabolism, and mitochondrial function), suggested less available energy and nutrients in RSM pigs. Furthermore, several genes supported more pronounced proteolysis (ABTB1, OTUD1, PADI2, SPP1) and reduced protein synthesis (THBS1, HSF4, AP1S2) in RSM muscle tissue. In parallel, higher levels of NR4A3, PDK4 and FGF21, and a drop in adropin, ELOVL6 and CIDEC/FSP27 indicated increased lipolysis and fatty acid oxidation, reflective of lower dressing percentage. Finally, pigs exposed to RSM showed greater expression level of genes responsive to oxidative stress, indicated by upregulation of GPX1, GPX2, and TXNIP.

Advances in next-generation sequencing technologies and the development of genome-reduced representation protocols have opened the way to genome-wide population studies in non-model species. However, species with large genomes remain challenging, hampering the development of genomic resources for a number of taxa including marine arthropods. Here, we developed a genome-reduced representation method for the ecologically important marine copepod Calanus finmarchicus (haploid genome size of 6.34 Gbp). We optimized a capture enrichment-based protocol based on 2656 single-copy genes, yielding a total of 154 087 high-quality SNPs in C. finmarchicus including 62 372 in common among the three locations tested. The set of capture probes was also successfully applied to the congeneric C. glacialis . Preliminary analyses of these markers revealed similar levels of genetic diversity between the two Calanus species, while populations of C. glacialis showed stronger genetic structure compared to C. finmarchicus . Using this powerful set of markers, we did not detect any evidence of hybridization between C. finmarchicus and C. glacialis . Finally, we propose a shortened version of our protocol, offering a promising solution for population genomics studies in non-model species with large genomes.

Resistance towards deltamethrin (DMT) in the crustacean ectoparasite Lepeophtheirus salmonis (Caligidae) is a problem on fish farms lining the North Atlantic Ocean. Two Norwegian strains with different susceptibility towards DMT were crossed in the parental generation (P0), females from a sensitive strain were crossed with males from a resistant strain and vice versa. Individual susceptibility towards DMT was assessed in the second filial generation (F2). DMT resistance was only found in F2 descendants when the P0 females were from the resistant strain, pointing to maternal inheritance. Since maternal inheritance might be linked to the mitochondrial (mt) genome, the nucleotide sequences and the gene expressions of mt-genes were analysed. Twenty non-synonymous single nucleotide polymorphisms (SNPs) were identified in mt-transcripts from resistant F2 parasites, including SNPs in two cytochrome C oxidase subunits (COX1 and COX3) and two subunits of the NADH dehydrogenase complex (ND1 and ND5) previously linked to DMT resistance in the salmon louse. Differential expression analysis between the sensitive and resistant strain revealed strain effect in seven out of twelve mt-genes. The current study also show that DNA fragmentation (indicating apoptosis) was affected by DMT exposure in skeletal muscle tissue and that resistant parasites undergo less apoptosis than sensitive parasites.

Background Teleosts of the genus Seriola, commonly known as amberjacks, are of high commercial value in international markets due to their flesh quality and worldwide distribution. The Seriola species of interest to Mediterranean aquaculture is the greater amberjack (Seriola dumerili). This species holds great potential for the aquaculture industry, but in captivity, reproduction has proved to be challenging, and observed growth dysfunction hinders their domestication. Insights into molecular mechanisms may contribute to a better understanding of traits like growth and sex, but investigations to unravel the molecular background of amberjacks have begun only recently. Findings Illumina HiSeq sequencing generated a high-coverage greater amberjack genome sequence comprising 45 909 scaffolds. Comparative mapping to the Japanese yellowtail (Seriola quinqueriadiata) and to the model species medaka (Oryzias latipes) allowed the generation of in silico groups. Additional gonad transcriptome sequencing identified sex-biased transcripts, including known sex-determining and differentiation genes. Investigation of the muscle transcriptome of slow-growing individuals showed that transcripts involved in oxygen and gas transport were differentially expressed compared with fast/normal-growing individuals. On the other hand, transcripts involved in muscle functions were found to be enriched in fast/normal-growing individuals. Conclusion The present study provides the first insights into the molecular background of male and female amberjacks and of fast- and slow-growing fish. Therefore, valuable molecular resources have been generated in the form of a first draft genome and a reference transcriptome. Sex-biased genes, which may also have roles in sex determination or differentiation, and genes that may be responsible for slow growth are suggested.

Background: Activated T helper type 2 (Th2) cells are believed to play a pivotal role in allergic airway inflammation, but which cells attract and activate Th2 cells locally have not been fully determined. Recently, it was shown in an experimental human model of allergic rhinitis (AR) that activated monocytes rapidly accumulate in the nasal mucosa after local allergen challenge, where they promote recruitment of Th2 cells and eosinophils. Objective: To investigate whether monocytes are recruited to the lungs in paediatric asthma. Methods: Tissue samples obtained from children and adolescents with fatal asthma attack (n = 12), age-matched non-atopic controls (n = 9) and allergen-challenged AR patients (n = 8) were subjected to in situ immunostaining. Results: Monocytes, identified as CD68+S100A8/A9+ cells, were significantly increased in the lower airway mucosa and in the alveoli of fatal asthma patients compared with control individuals. Interestingly, cellular aggregates containing CD68+S100A8/A9+ monocytes obstructing the lumen of bronchioles were found in asthmatics (8 out of 12) but not in controls. Analysing tissue specimens from challenged AR patients, we confirmed that co-staining with CD68 and S100A8/A9 was a valid method to identify recently recruited monocytes. We also showed that the vast majority of accumulating monocytes both in the lungs and in the nasal mucosa expressed matrix metalloproteinase 10, suggesting that this protein may be involved in their migration within the tissue. Conclusions and clinical relevance: Monocytes accumulated in the lungs of children and adolescents with fatal asthma attack. This finding strongly suggests that monocytes are directly involved in the immunopathology of asthma and that these pro-inflammatory cells are potential targets for therapy.

BackgroundAdvances in sequencing technologies and bioinformatics have made the analysis of microbial communities almost routine. Nonetheless, the need remains to improve on the techniques used for gathering such data, including increasing throughput while lowering cost, and benchmarking the techniques so that potential sources of bias can be better characterized.ResultsWe present a triple-index amplicon sequencing strategy that uses a two-stage PCR protocol. The strategy was extensively benchmarked through analysis of a mock community in order to assess biases introduced by sample indexing, number of PCR cycles, and template concentration. We further evaluated the method through re-sequencing of a standardized environmental sample. Finally, we evaluated our protocol on a set of fecal samples from a small cohort of healthy adults, demonstrating good performance in a realistic experimental setting. Between-sample variation was mainly related to batch effects, such as DNA extraction, while sample indexing was also a significant source of bias. PCR cycle number strongly influenced chimera formation and affected relative abundance estimates of species with high GC content. Libraries were sequenced using the Illumina HiSeq and MiSeq platforms to demonstrate that this protocol is highly scalable to sequence thousands of samples at a very low cost.ConclusionsHere, we provide the most comprehensive study of performance and bias inherent to a 16S rRNA gene amplicon sequencing method to date. Triple-indexing greatly reduces the number of long custom DNA oligos required for library preparation, while the inclusion of variable length heterogeneity spacers minimizes the need for PhiX spike-in. This design results in a significant cost reduction of highly multiplexed amplicon sequencing. The biases we characterize highlight the need for highly standardized protocols. Reassuringly, we find that the biological signal is a far stronger structuring factor than the various sources of bias.

Background Advances in sequencing technologies and bioinformatics have made the analysis of microbial communities almost routine. Nonetheless, the need remains to improve on the techniques used for gathering such data, including increasing throughput while lowering cost and benchmarking the techniques so that potential sources of bias can be better characterized. Methods We present a triple-index amplicon sequencing strategy to sequence large numbers of samples at significantly lower c ost and in a shorter timeframe compared to existing methods. The design employs a two-stage PCR protocol, incorpo rating three barcodes to each sample, with the possibility to add a fourth-index. It also includes heterogeneity spacers to overcome low complexity issues faced when sequencing amplicons on Illumina platforms. Results The library preparation method was extensively benchmarked through analysis of a mock community in order to assess biases introduced by sample indexing, number of PCR cycles, and template concentration. We further evaluated the method through re-sequencing of a standardized environmental sample. Finally, we evaluated our protocol on a set of fecal samples from a small cohort of healthy adults, demonstrating good performance in a realistic experimental setting. Between-sample variation was mainly related to batch effects, such as DNA extraction, while sample indexing was also a significant source of bias. PCR cycle number strongly influenced chimera formation and affected relative abundance estimates of species with high GC content. Libraries were sequenced using the Illumina HiSeq and MiSeq platforms to demonstrate that this protocol is highly scalable to sequence thousands of samples at a very low cost. Conclusions Here, we provide the most comprehensive study of performance and bias inherent to a 16S rRNA gene amplicon sequencing method to date. Triple-indexing greatly reduces the number of long custom DNA oligos required for library preparation, while the inclusion of variable length heterogeneity spacers minimizes the need for PhiX spike-in. This design results in a significant cost reduction of highly multiplexed amplicon sequencing. The biases we characterize highlight the need for highly standardized protocols. Reassuringly, we find that the biological signal is a far stronger structuring factor than the various sources of bias.

Antigen presenting cells (APCs) in the thymus play an essential role in the establishment of central tolerance, i.e. the generation of a repertoire of functional and self-tolerant T cells to prevent autoimmunity. In this study, we have compared the transcriptomes of four primary APCs from human thymus (mTECs, CD19+ B cells, CD141+ and CD123+ DCs). We investigated a set of genes including the HLA genes, genes encoding transcriptional regulators and finally, tissue-enriched genes, i.e, genes with a five-fold higher expression in a particular human tissue. We show that thymic CD141+ DCs express the highest levels of all classical HLA genes and 67% (14/21) of the HLA class I and II pathway genes investigated in this study. CD141+ DCs also expressed the highest levels of the transcriptional regulator DEAF1, whereas AIRE and FEZF2 expression were mainly found in primary human mTECs. We found expression of “tissue enriched genes” from the Human Protein Atlas (HPA) in all four APC types, but the mTECs were clearly dominating in the number of uniquely expressed tissue enriched genes (20% in mTECs, 7% in CD19+ B cells, 4% in CD123+ DCs and 2% in CD141+ DCs). The tissue enriched genes also overlapped with reported human autoantigens. This is, to our knowledge, the first study that performs RNA sequencing of mTECs, CD19+ B cells, CD141+ and CD123+ DCs isolated from the same individuals and provides insight into the transcriptomes of these human thymic APCs.

Background Teleosts are exposed to a broad range of external stimuli, which may be either of acute or chronic nature. The larval phase of certain fish species offer a unique opportunity to study the interactions between genes and environmental factors during early life. The present study investigates the effects of early-life events, applied at different time points of early ontogeny (first feeding, flexion and development of all fins; Phase 1) as well as on the subsequent juvenile stage after the application of an additional acute stressor (Phase 2) in the gilthead sea bream (Sparus aurata), a commercially important European aquaculture species. Animal performance, the cortisol response and gene expression patterns during early development as well as on the subsequent phases (juveniles) after the application of additional acute stressors were investigated. Results Significant differences on fish performance were found only for juveniles exposed to early-life events at the phase of the formation of all fins. On the transcriptome level distinct expression patterns were obtained for larvae as well as for juveniles with the most divergent expression pattern found to be again at the phase of the development of all fins, which showed to have also an impact later on in the acute stress response of juveniles. Conclusions The present study showed that applying an early-life protocol, characterized by the unpredictable, variable and moderate intensity of the applied stimuli provides a relative realistic model to evaluate the impact of daily aquaculture practices on fish performance. In addition, the power of investigating global gene expression patterns is shown, providing significant insights regarding the response of early-life events during development and as juveniles after the application of extra acute stressors. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2874-0) contains supplementary material, which is available to authorized users.

Background Several miRNAs are known to control myogenesis in vertebrates. Some of them are specifically expressed in muscle while others have a broader tissue expression but are still involved in establishing the muscle phenotype. In teleosts, water temperature markedly affects embryonic development and larval growth. It has been previously shown that higher embryonic temperatures promoted faster development and increased size of Senegalese sole (Solea senegalensis) larvae relatively to a lower temperature. The role of miRNAs in thermal-plasticity of growth is hitherto unknown. Hence, we have used high-throughput SOLiD sequencing to determine potential changes in the miRNA transcriptome in Senegalese sole embryos that were incubated at 15°C or 21°C until hatching and then reared at a common temperature of 21°C. Results We have identified 320 conserved miRNAs in Senegalese sole, of which 48 had not been previously described in teleosts. mir-17a-5p, mir-26a, mir-130c, mir-206-3p, mir-181a-5p, mir-181a-3p and mir-199a-5p expression levels were further validated by RT- qPCR. The majority of miRNAs were dynamically expressed during early development, with peaks of expression at pre-metamorphosis or metamorphosis. Also, a higher incubation temperature (21°C) was associated with expression of some miRNAs positively related with growth (e.g., miR-17a, miR-181-5p and miR-206) during segmentation and at hatching. Target prediction revealed that these miRNAs may regulate myogenesis through MAPK and mTOR pathways. Expression of miRNAs involved in lipid metabolism and energy production (e.g., miR-122) also differed between temperatures. A miRNA that can potentially target calpain (miR-181-3p), and therefore negatively regulate myogenesis, was preferentially expressed during segmentation at 15°C compared to 21°C. Conclusions Temperature has a strong influence on expression of miRNAs during embryonic and larval development in fish. Higher expression levels of miR-17a, miR-181-5p and miR-206-3p and down-regulation of miR-181a-3p at 21°C may promote myogenesis and are in agreement with previous studies in Senegalese sole, which reported enhanced growth at higher embryonic temperatures compared to 15°C. Moreover, miRNAs involved in lipid metabolism and energy production may also contribute to increased larval growth at 21°C compared to 15°C. Taken together, our data indicate that miRNAs may play a role in temperature-induced phenotypic plasticity of growth in teleosts. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-525) contains supplementary material, which is available to authorized users.

Background Toll-like receptors (Tlrs) are major molecular pattern recognition receptors of the innate immune system. Atlantic cod (Gadus morhua) is the first vertebrate known to have lost most of the mammalian Tlr orthologues, particularly all bacterial recognising and other cell surface Tlrs. On the other hand, its genome encodes a unique repertoire of teleost-specific Tlrs. The aim of this study was to investigate if these duplicate Tlrs have been retained through adaptive evolution to compensate for the lack of other cell surface Tlrs in the cod genome. Results In this study, one tlr21, 12 tlr22 and two tlr23 genes representing the teleost-specific Tlr family have been cloned and characterised in cod. Phylogenetic analysis grouped all tlr22 genes under a single clade, indicating that the multiple cod paralogues have arisen through lineage-specific duplications. All tlrs examined were transcribed in immune-related tissues as well as in stomach, gut and gonads of adult cod and were differentially expressed during early development. These tlrs were also differentially regulated following immune challenge by immersion with Vibrio anguillarum, indicating their role in the immune response. An increase in water temperature from 4 to 12°C was associated with a 5.5-fold down-regulation of tlr22d transcript levels in spleen. Maximum likelihood analysis with different evolution models revealed that tlr22 genes are under positive selection. A total of 24 codons were found to be positively selected, of which 19 are in the ligand binding region of ectodomain. Conclusion Positive selection pressure coupled with experimental evidence of differential expression strongly support the hypothesis that teleost-specific tlr paralogues in cod are undergoing neofunctionalisation and can recognise bacterial pathogen-associated molecular patterns to compensate for the lack of other cell surface Tlrs.

Toll-like receptors (TLRs) play a crucial role in host defence, since they trigger immune response following recognition of pathogen-associated molecular patterns (PAMPs) in potential infectious agents. TLRs have been found in numerous organisms, including mammals, birds and teleosts. Some TLR members are commonly retained across all species, whilst others were lost, gained or diverged independently during evolution. Our knowledge about the evolution and specific functions of tlr21, tlr22 and tlr23 in teleosts are still scarce. Phylogenetic analysis of 18 tlr13, tlr21, tlr22 and tlr23 genes from 9 different fish species divided them in two groups. All tlr21 genes were under the first clade, while the second comprised tlr22, tlr23 and tlr13 from Atlantic salmon. Evidence of positive selection was detected at three sites within the leucine-rich repeat regions of Tlr22, which may influence PAMP recognition. Immunostimulation experiments revealed that expression of zebrafish tlr22 is modulated by several unrelated PAMPs. Up to a 3-fold increase in tlr21 and tlr22 expression was detected in larvae exposed to immunostimulants such as lipopolysaccharide, peptidoglycan or poly I:C. We found that zebrafish tlrs are expressed mainly in immune-related organs, such as spleen and kidney as well as in testis and temperature stress did not have an effect on the expression of tlr21 and tlr22 in the early stages of development in zebrafish larvae. Our data indicates that these teleost tlrs may play a role in innate host defence. In particular, tlr22 is evolving under positive selection, which indicates functional diversification and adaptation of the response to different PAMPs.

LSH, a member of the SNF2 family of chromatin remodeling ATPases encoded by the Hells gene, is essential for normal levels of DNA methylation in the mammalian genome. While the role of LSH in the methylation of repetitive DNA sequences is well characterized, its contribution to the regulation of DNA methylation and the expression of protein-coding genes has not been studied in detail. In this report we investigate genome-wide patterns of DNA methylation at gene promoters in Hells−/− mouse embryonic fibroblasts (MEFs). We find that in the absence of LSH, DNA methylation is lost or significantly reduced at ∼20% of all normally methylated promoter sequences. As a consequence, a large number of genes are misexpressed in Hells−/− MEFs. Comparison of Hells−/− MEFs with wild-type MEFs and embryonic stem (ES) cells suggests that LSH is important for de novo DNA methylation events that accompany the establishment and differentiation of embryonic lineage cells. We further show that the generation of normal DNA methylation patterns and stable gene silencing at specific promoters require cooperation between LSH and the G9a/GLP complex of histone methylases. At such loci, G9a recruitment is compromised when LSH is absent or greatly reduced. Taken together, our data suggest a mechanism whereby LSH promotes binding of DNA methyltransferases and the G9a/GLP complex to specific loci and facilitates developmentally programmed DNA methylation and stable gene silencing during lineage commitment and differentiation.

SummaryPrenatal paracetamol exposure has been associated with neurodevelopmental outcomes in childhood. Pharmacoepigenetic studies show differences in cord blood DNA methylation between paracetamol exposed and unexposed neonates. However, causal implications and impact of long-term prenatal long-term paracetamol exposure on brain development remain unclear. Using a multi-omics approach, we investigated the effects of paracetamol on a model of early human brain development. We exposed human embryonic stem cells undergoing in vitro neuronal differentiation to daily media changes with paracetamol concentrations corresponding to maternal therapeutic doses. Single-cell RNA-seq and ATAC-seq integration identified paracetamol-induced chromatin-opening changes linked to gene expression. Differentially methylated and/or expressed genes were involved in signalling, neurotransmission, and cell fate-determination trajectories. Some genes involved in neuronal injury and development-specific pathways, such asKCNE3, overlapped with differentially methylated genes previously identified in cord blood associated with prenatal paracetamol exposure. Our data suggest that paracetamol may play a causal role in impaired neurodevelopment.Graphical Abstract

Background ChIP-seq is the primary technique used to investigate genome-wide protein-DNA interactions. As part of this procedure, immunoprecipitated DNA must undergo “library preparation” to enable subsequent high-throughput sequencing. To facilitate the analysis of biopsy samples and rare cell populations, there has been a recent proliferation of methods allowing sequencing library preparation from low-input DNA amounts. However, little information exists on the relative merits, performance, comparability and biases inherent to these procedures. Notably, recently developed single-cell ChIP procedures employing microfluidics must also employ library preparation reagents to allow downstream sequencing. Results In this study, seven methods designed for low-input DNA/ChIP-seq sample preparation (Accel-NGS® 2S, Bowman-method, HTML-PCR, SeqPlex™, DNA SMART™, TELP and ThruPLEX®) were performed on five replicates of 1 ng and 0.1 ng input H3K4me3 ChIP material, and compared to a “gold standard” reference PCR-free dataset. The performance of each method was examined for the prevalence of unmappable reads, amplification-derived duplicate reads, reproducibility, and for the sensitivity and specificity of peak calling. Conclusions We identified consistent high performance in a subset of the tested reagents, which should aid researchers in choosing the most appropriate reagents for their studies. Furthermore, we expect this work to drive future advances by identifying and encouraging use of the most promising methods and reagents. The results may also aid judgements on how comparable are existing datasets that have been prepared with different sample library preparation reagents. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-3135-y) contains supplementary material, which is available to authorized users.