Genome-wide SNP profiling of worldwide goat populations reveals strong partitioning of diversity and highlights post-domestication migration routes

BACKGROUND

Goat populations that are characterized within the AdaptMap project cover a large part of the worldwide distribution of this species and provide the opportunity to assess their diversity at a global scale. We analysed genome-wide 50 K single nucleotide polymorphism (SNP) data from 144 populations to describe the global patterns of molecular variation, compare them to those observed in other livestock species, and identify the drivers that led to the current distribution of goats.

RESULTS

A high degree of genetic variability exists among the goat populations studied. Our results highlight a strong partitioning of molecular diversity between and within continents. Three major gene pools correspond to goats from Europe, Africa and West Asia. Dissection of sub-structures disclosed regional gene pools, which reflect the main post-domestication migration routes. We also identified several exchanges, mainly in African populations, and which often involve admixed and cosmopolitan breeds. Extensive gene flow has taken place within specific areas (e.g., south Europe, Morocco and Mali-Burkina Faso-Nigeria), whereas elsewhere isolation due to geographical barriers (e.g., seas or mountains) or human management has decreased local gene flows.

CONCLUSIONS

After domestication in the Fertile Crescent in the early Neolithic era (ca. 12,000 YBP), domestic goats that already carried differentiated gene pools spread to Europe, Africa and Asia. The spread of these populations determined the major genomic background of the continental populations, which currently have a more marked subdivision than that observed in other ruminant livestock species. Subsequently, further diversification occurred at the regional level due to geographical and reproductive isolation, which was accompanied by additional migrations and/or importations, the traces of which are still detectable today. The effects of breed formation were clearly detected, particularly in Central and North Europe. Overall, our results highlight a remarkable diversity that occurs at the global scale and is locally partitioned and often affected by introgression from cosmopolitan breeds. These findings support the importance of long-term preservation of goat diversity, and provide a useful framework for investigating adaptive introgression, directing genetic improvement and choosing breeding targets.

Strengthening insights into host responses to mastitis infection in ruminants by combining heterogeneous microarray data sources

Background

Gene expression profiling studies of mastitis in ruminants have provided key but fragmented knowledge for the understanding of the disease. A systematic combination of different expression profiling studies via meta-analysis techniques has the potential to test the extensibility of conclusions based on single studies. Using the program Pointillist, we performed meta-analysis of transcription-profiling data from six independent studies of infections with mammary gland pathogens, including samples from cattle challenged in vivo with S. aureus, E. coli, and S. uberis, samples from goats challenged in vivo with S. aureus, as well as cattle macrophages and ovine dendritic cells infected in vitro with S. aureus. We combined different time points from those studies, testing different responses to mastitis infection: overall (common signature), early stage, late stage, and cattle-specific.

Results

Ingenuity Pathway Analysis of affected genes showed that the four meta-analysis combinations share biological functions and pathways (e.g. protein ubiquitination and polyamine regulation) which are intrinsic to the general disease response. In the overall response, pathways related to immune response and inflammation, as well as biological functions related to lipid metabolism were altered. This latter observation is consistent with the milk fat content depression commonly observed during mastitis infection. Complementarities between early and late stage responses were found, with a prominence of metabolic and stress signals in the early stage and of the immune response related to the lipid metabolism in the late stage; both mechanisms apparently modulated by few genes, including XBP1 and SREBF1.

The cattle-specific response was characterized by alteration of the immune response and by modification of lipid metabolism. Comparison of E. coli and S. aureus infections in cattle in vivo revealed that affected genes showing opposite regulation had the same altered biological functions and provided evidence that E. coli caused a stronger host response.

Conclusions

This meta-analysis approach reinforces previous findings but also reveals several novel themes, including the involvement of genes, biological functions, and pathways that were not identified in individual studies. As such, it provides an interesting proof of principle for future studies combining information from diverse heterogeneous sources.

New Insights on Water Buffalo Genomic Diversity and Post-Domestication Migration Routes From Medium Density SNP Chip Data

The domestic water buffalo is native to the Asian continent but through historical migrations and recent importations, nowadays has a worldwide distribution. The two types of water buffalo, i.e., river and swamp, display distinct morphological and behavioral traits, different karyotypes and also have different purposes and geographical distributions. River buffaloes from Pakistan, Iran, Turkey, Egypt, Romania, Bulgaria, Italy, Mozambique, Brazil and Colombia, and swamp buffaloes from China, Thailand, Philippines, Indonesia and Brazil were genotyped with a species-specific medium-density 90K SNP panel. We estimated the levels of molecular diversity and described population structure, which revealed historical relationships between populations and migration events. Three distinct gene pools were identified in pure river as well as in pure swamp buffalo populations. Genomic admixture was seen in the Philippines and in Brazil, resulting from importations of animals for breed improvement. Our results were largely consistent with previous archeological, historical and molecular-based evidence for two independent domestication events for river- and swamp-type buffaloes, which occurred in the Indo-Pakistani region and close to the China/Indochina border, respectively. Based on a geographical analysis of the distribution of diversity, our evidence also indicated that the water buffalo spread out of the domestication centers followed two major divergent migration directions: river buffaloes migrated west from the Indian sub-continent while swamp buffaloes migrated from northern Indochina via an east-south-eastern route. These data suggest that the current distribution of water buffalo diversity has been shaped by the combined effects of multiple migration events occurred at different stages of the post-domestication history of the species.

Response of the goat mammary gland to infection with Staphylococcus aureus revealed by gene expression profiling in milk somatic and white blood cells

Background

S. aureus is one of the main pathogens responsible for the intra-mammary infection in dairy ruminants. Although much work has been carried out to understand the complex physiological and cellular events that occur in the mammary gland in response to S. aureus, the protective mechanisms are still poorly understood. The objectives of the present study were to investigate gene expression during the early response of the goat mammary gland to an experimental challenge with S. aureus, in order to better understand the local and systemic response and to compare them in two divergent lines of goat selected for high and low milk somatic cell scores.

Results

No differences in gene expression were found between high and low SCS (Somatic Cells Score) selection lines. Analysing the two groups together, an expression of 300 genes were found to change from T0 before infection, and T4 at 24 hours and T5 at 30 hours following challenge. In blood derived white blood cells 8 genes showed increased expression between T0 and T5 and 1 gene has reduced expression. The genes showing the greatest increase in expression following challenge (5.65 to 3.16 fold change) play an important role in (i) immune and inflammatory response (NFKB1, TNFAIP6, BASP1, IRF1, PLEK, BATF3); (ii) the regulation of innate resistance to pathogens (PTX3); and (iii) the regulation of cell metabolism (CYTH4, SLC2A6, ARG2). The genes with reduced expression (−1.5 to −2.5 fold) included genes involved in (i) lipid metabolism (ABCG2, FASN), (ii) chemokine, cytokine and intracellular signalling (SPPI), and (iii) cell cytoskeleton and extracellular matrix (KRT19).

Conclusions

Analysis of genes with differential expression following infection showed an inverse relationship between immune response and lipid metabolism in the early response of the mammary gland to the S. aureus challenge. PTX3 showed a large change in expression in both milk and blood, and is therefore a candidate for further studies on immune response associated with mastitis.

OsWRKY22, a monocot WRKY gene, plays a role in the resistance response to blast

With the aim of identifying novel regulators of host and nonhost resistance to fungi in rice, we carried out a systematic mutant screen of mutagenized lines. Two mutant wrky22 knockout lines revealed clear-cut enhanced susceptibility to both virulent and avirulent Magnaporthe oryzae strains and altered cellular responses to nonhost Magnaporthe grisea and Blumeria graminis fungi. In addition, the analysis of the pathogen responses of 24 overexpressor OsWRKY22 lines revealed enhanced resistance phenotypes on infection with virulent M. oryzae strain, confirming that OsWRKY22 is involved in rice resistance to blast. Bioinformatic analyses determined that the OsWRKY22 gene belongs to a well-defined cluster of monocot-specific WRKYs. The co-regulatory analysis revealed no significant co-regulation of OsWRKY22 with a representative panel of OsWRKYs, supporting its unique role in a series of transcriptional responses. In contrast, inquiring a subset of biotic stress-related Affymetrix data, a large number of resistance and defence-related genes were found to be putatively co-expressed with OsWRKY22. Taken together, all gathered experimental evidence places the monocot-specific OsWRKY22 gene at the convergence point of signal transduction circuits in response to both host and nonhost fungi encountering rice plants.

Molecular diversity, population structure, and linkage disequilibrium in a worldwide collection of tobacco (Nicotiana tabacum L.) germplasm

BACKGROUND

The goals of our study were to assess the phylogeny and the population structure of tobacco accessions representing a wide range of genetic diversity; identify a subset of accessions as a core collection capturing most of the existing genetic diversity; and estimate, in the tobacco core collection, the extent of linkage disequilibrium (LD) in seven genomic regions using simple sequence repeat (SSR) markers. To this end, a collection of accessions were genotyped with SSR markers. Molecular diversity was evaluated and LD was analyzed across seven regions of the genome.

RESULTS

A genotyping database for 312 tobacco accessions was profiled with 49 SSR markers. Principal Coordinate Analysis (PCoA) and Bayesian cluster analysis revealed structuring of the tobacco population with regard to commercial classes and six main clades were identified, which correspond to "Oriental", Flue-Cured", "Burley", "Dark", "Primitive", and "Other" classes. Pairwise kinship was calculated between accessions, and an overall low level of co-ancestry was observed. A set of 89 genotypes was identified that captured the whole genetic diversity detected at the 49 loci. LD was evaluated on these genotypes, using 422 SSR markers mapping on seven linkage groups. LD was estimated as squared correlation of allele frequencies (r2). The pattern of intrachromosomal LD revealed that in tobacco LD extended up to distances as great as 75 cM with r2 > 0.05 or up to 1 cM with r2 > 0.2. The pattern of LD was clearly dependent on the population structure.

CONCLUSIONS

A global population of tobacco is highly structured. Clustering highlights the accessions with the same market class. LD in tobacco extends up to 75 cM and is strongly dependent on the population structure.

Adaptive introgression from indicine cattle into white cattle breeds from Central Italy

Cattle domestication occurred at least twice independently and gave rise to the modern taurine and indicine cattle breeds. European cattle diversity is generally dominated by taurine cattle, although elevated levels of indicine ancestry have been recorded in several breeds from southern Europe. Here we use genome-wide high-density SNP genotyping data to investigate the taurine and indicine ancestry in southern European cattle, based on a dataset comprising 508 individuals from 23 cattle breeds of taurine, indicine and mixed ancestry, including three breeds from Central Italy known to exhibit the highest levels of indicine introgression among southern European breeds. Based on local genomic ancestry analyses, we reconstruct taurine and indicine ancestry genome-wide and along chromosomes. We scrutinise local genomic introgression signals and identify genomic regions that have introgressed from indicine into taurine cattle under positive selection, harbouring genes with functions related to body size and feed efficiency. These findings suggest that indicine-derived traits helped enhance Central Italian cattle through adaptive introgression. The identified genes could provide genomic targets for selection for improved cattle performance. Our findings elucidate the key role of adaptive introgression in shaping the phenotypic features of modern cattle, aided by cultural and livestock exchange among historic human societies.

Integrative analysis of methylomic and transcriptomic data in fetal sheep muscle tissues in response to maternal diet during pregnancy

BACKGROUND

Numerous studies have established a link between maternal diet and the physiological and metabolic phenotypes of their offspring. In previous studies in sheep, we demonstrated that different maternal diets altered the transcriptome of fetal tissues. However, the mechanisms underlying transcriptomic changes are poorly understood. DNA methylation is an epigenetic mark regulating transcription and is largely influenced by dietary components of the one-carbon cycle that generate the methyl group donor, SAM. Therefore, in the present study, we tested whether different maternal diets during pregnancy would alter the DNA methylation and gene expression patterns in fetal tissues.

RESULTS

Pregnant ewes were randomly divided into two groups which received either hay or corn diet from mid-gestation (day 67 ± 5) until day 131 ± 1 when fetuses were collected by necropsy. A total of 1516 fetal longissimus dorsi (LD) tissues were used for DNA methylation analysis and gene expression profiling. Whole genome DNA methylation using methyl-binding domain enrichment analysis revealed 60 differentially methylated regions (DMRs) between hay and corn fetuses with 39 DMRs more highly methylated in the hay fetuses vs. 21 DMRs more highly methylated in the corn fetuses. Three DMRs (LPAR3, PLIN5-PLIN4, and the differential methylation of a novel lincRNA) were validated using bisulfite sequencing. These DMRs were associated with differential gene expression. Additionally, significant DNA methylation differences were found at the single CpG level. Integrative methylome and transcriptome analysis revealed an association between gene expression and inter-/intragenic methylated regions. Furthermore, intragenic DMRs were found to be associated with expression of neighboring genes.

CONCLUSIONS

The findings of this study imply that maternal diet from mid- to late-gestation can shape the epigenome and the transcriptome of fetal tissues, and putatively affect phenotypes of the lambs.

Methylome Patterns of Cattle Adaptation to Heat Stress

Heat stress has a detrimental impact on cattle health, welfare and productivity by affecting gene expression, metabolism and immune response, but little is known on the epigenetic mechanisms mediating the effect of temperature at the cellular and organism level. In this study, we investigated genome-wide DNA methylation in blood samples collected from 5 bulls of the heat stress resilient Nellore breed and 5 bulls of the Angus that are more heat stress susceptible, exposed to the sun and high temperature-high humidity during the summer season of the Brazilian South-East region. The methylomes were analyzed during and after the exposure by Reduced Representation Bisulfite Sequencing, which provided genome-wide single-base resolution methylation profiles. Significant methylation changes between stressful and recovery periods were observed in 819 genes. Among these, 351 were only seen in Angus, 366 were specific to Nellore, and 102 showed significant changes in methylation patterns in both breeds. KEGG and Gene Ontology (GO) enrichment analyses showed that responses were breed-specific. Interestingly, in Nellore significant genes and pathways were mainly involved in stress responses and cellular defense and were under methylated during heat stress, whereas in Angus the response was less focused. These preliminary results suggest that heat challenge induces changes in methylation patterns in specific loci, which should be further scrutinized to assess their role in heat tolerance.

Combining Landscape Genomics and Ecological Modelling to Investigate Local Adaptation of Indigenous Ugandan Cattle to East Coast Fever

East Coast fever (ECF) is a fatal sickness affecting cattle populations of eastern, central, and southern Africa. The disease is transmitted by the tick Rhipicephalus appendiculatus, and caused by the protozoan Theileria parva parva, which invades host lymphocytes and promotes their clonal expansion. Importantly, indigenous cattle show tolerance to infection in ECF-endemically stable areas. Here, the putative genetic bases underlying ECF-tolerance were investigated using molecular data and epidemiological information from 823 indigenous cattle from Uganda. Vector distribution and host infection risk were estimated over the study area and subsequently tested as triggers of local adaptation by means of landscape genomics analysis. We identified 41 and seven candidate adaptive loci for tick resistance and infection tolerance, respectively. Among the genes associated with the candidate adaptive loci are PRKG1 and SLA2. PRKG1 was already described as associated with tick resistance in indigenous South African cattle, due to its role into inflammatory response. SLA2 is part of the regulatory pathways involved into lymphocytes' proliferation. Additionally, local ancestry analysis suggested the zebuine origin of the genomic region candidate for tick resistance.

VarGoats project: a dataset of 1159 whole-genome sequences to dissect Capra hircus global diversity

Background

Since their domestication 10,500 years ago, goat populations with distinctive genetic backgrounds have adapted to a broad variety of environments and breeding conditions. The VarGoats project is an international 1000-genome resequencing program designed to understand the consequences of domestication and breeding on the genetic diversity of domestic goats and to elucidate how speciation and hybridization have modeled the genomes of a set of species representative of the genus Capra.

Findings

A dataset comprising 652 sequenced goats and 507 public goat sequences, including 35 animals representing eight wild species, has been collected worldwide. We identified 74,274,427 single nucleotide polymorphisms (SNPs) and 13,607,850 insertion-deletions (InDels) by aligning these sequences to the latest version of the goat reference genome (ARS1). A Neighbor-joining tree based on Reynolds genetic distances showed that goats from Africa, Asia and Europe tend to group into independent clusters. Because goat breeds from Oceania and Caribbean (Creole) all derive from imported animals, they are distributed along the tree according to their ancestral geographic origin.

Conclusions

We report on an unprecedented international effort to characterize the genome-wide diversity of domestic goats. This large range of sequenced individuals represents a unique opportunity to ascertain how the demographic and selection processes associated with post-domestication history have shaped the diversity of this species. Data generated for the project will also be extremely useful to identify deleterious mutations and polymorphisms with causal effects on complex traits, and thus will contribute to new knowledge that could be used in genomic prediction and genome-wide association studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00659-6.

Genome-Wide DNA Methylation and Gene Expression Profiles in Cows Subjected to Different Stress Level as Assessed by Cortisol in Milk

Dairy cattle health, wellbeing and productivity are deeply affected by stress. Its influence on metabolism and immune response is well known, but the underlying epigenetic mechanisms require further investigation. In this study, we compared DNA methylation and gene expression signatures between two dairy cattle populations falling in the high- and low-variant tails of the distribution of milk cortisol concentration (MC), a neuroendocrine marker of stress in dairy cows. Reduced Representation Bisulfite Sequencing was used to obtain a methylation map from blood samples of these animals. The high and low groups exhibited similar amounts of methylated CpGs, while we found differences among non-CpG sites. Significant methylation changes were detected in 248 genes. We also identified significant fold differences in the expression of 324 genes. KEGG and Gene Ontology (GO) analysis showed that genes of both groups act together in several pathways, such as nervous system activity, immune regulatory functions and glucocorticoid metabolism. These preliminary results suggest that, in livestock, cortisol secretion could act as a trigger for epigenetic regulation and that peripheral changes in methylation can provide an insight into central nervous system functions.

Ontology-oriented retrieval of putative microRNAs in Vitis vinifera via GrapeMiRNA: a web database of de novo predicted grape microRNAs

BACKGROUND

Two complete genome sequences are available for Vitis vinifera Pinot noir. Based on the sequence and gene predictions produced by the IASMA, we performed an in silico detection of putative microRNA genes and of their targets, and collected the most reliable microRNA predictions in a web database. The application is available at http://www.itb.cnr.it/ptp/grapemirna/.

DESCRIPTION

The program FindMiRNA was used to detect putative microRNA genes in the grape genome. A very high number of predictions was retrieved, calling for validation. Nine parameters were calculated and, based on the grape microRNAs dataset available at miRBase, thresholds were defined and applied to FindMiRNA predictions having targets in gene exons. In the resulting subset, predictions were ranked according to precursor positions and sequence similarity, and to target identity. To further validate FindMiRNA predictions, comparisons to the Arabidopsis genome, to the grape Genoscope genome, and to the grape EST collection were performed. Results were stored in a MySQL database and a web interface was prepared to query the database and retrieve predictions of interest.

CONCLUSION

The GrapeMiRNA database encompasses 5,778 microRNA predictions spanning the whole grape genome. Predictions are integrated with information that can be of use in selection procedures. Tools added in the web interface also allow to inspect predictions according to gene ontology classes and metabolic pathways of targets. The GrapeMiRNA database can be of help in selecting candidate microRNA genes to be validated.

Newly-Discovered Neural Features Expand the Pathobiological Knowledge of Blastic Plasmacytoid Dendritic Cell Neoplasm

Simple Summary

For the first time, neuronal features are described in blastic plasmacytoid dendritic cell neoplasm (BPDCN) by a complex array of molecular techniques, including microRNA and gene expression profiling, RNA and Chromatin immunoprecipitation sequencing, and immunohistochemistry. The discovery of unexpected neural features in BPDCN may change our vision of this disease, leading to the designing of a new BPDCN cell model and to re-thinking the relations occurring between BPDCN and nervous system. The observed findings contribute to explaining the extreme tumor aggressiveness and also to propose novel therapeutic targets. In view of this, the identification, in this work of new potential neural metastatic inducers might open the way to therapeutic approaches for BPDCN patients based on the use of anti-neurogenic agents.

Abstract

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is a rare and highly aggressive hematologic malignancy originating from plasmacytoid dendritic cells (pDCs). The microRNA expression profile of BPDCN was compared to that of normal pDCs and the impact of miRNA dysregulation on the BPDCN transcriptional program was assessed. MiRNA and gene expression profiling data were integrated to obtain the BPDCN miRNA-regulatory network. The biological process mainly dysregulated by this network was predicted to be neurogenesis, a phenomenon raising growing interest in solid tumors. Neurogenesis was explored in BPDCN by querying different molecular sources (RNA sequencing, Chromatin immunoprecipitation-sequencing, and immunohistochemistry). It was shown that BPDCN cells upregulated neural mitogen genes possibly critical for tumor dissemination, expressed neuronal progenitor markers involved in cell migration, exchanged acetylcholine neurotransmitter, and overexpressed multiple neural receptors that may stimulate tumor proliferation, migration and cross-talk with the nervous system. Most neural genes upregulated in BPDCN are currently investigated as therapeutic targets.

Epigenetic Variability Among Saffron Crocus (Crocus sativus L.) Accessions Characterized by Different Phenotypes

This work represents the first epigenomic study carried out on saffron crocus. Five accessions of saffron, showing differences in tepal pigmentation, yield of saffron and flowering time, were analyzed at the epigenetic level by applying a methylation-sensitive restriction enzyme-sequencing (MRE-seq) approach. Five accession-specific hypomethylomes plus a reference hypomethylome, generated by combining the sequence data from the single accessions, were obtained. Assembled sequences were annotated against existing online databases. In the absence of the Crocus genome, the rice genome was mainly used as the reference as it is the best annotated genome among monocot plants. Comparison of the hypomethylomes revealed many differentially methylated regions, confirming the high epigenetic variability present among saffron accessions, including sequences encoding for proteins that could be good candidates to explain the accessions' alternative phenotypes. In particular, transcription factors involved in flowering process (MADS-box and TFL) and for the production of pigments (MYB) were detected. Finally, by comparing the generated sequences of the different accessions, a high number of SNPs, likely having arisen as a consequence of the prolonged vegetative propagation, were detected, demonstrating surprisingly high genetic variability. Gene ontology (GO) was performed to map and visualize sequence polymorphisms located within the GOs and to compare their distributions among different accessions. As well as suggesting the possible existence of alternative phenotypes with a genetic basis, a clear difference in polymorphic GO is present among accessions based on their geographic origin, supporting a possible signature of selection in the Indian accession with respect to the Spanish ones.

Comparison of the response of mammary gland tissue from two divergent lines of goat with high and low milk somatic cell scores to an experimental Staphylococcus aureus infection

Mastitis represents one of the major economic and health threats to the livestock sector associated with reduction in milk quality, loss of production and is a major reason for culling. Somatic cell score (SCS) is used as a criterion in breeding programmes to select cows genetically less susceptible to mastitis. The relevance of SCS as a predictor of udder health and susceptibility to mastitis is still untested in goats. In this study, two lines of French Alpine goats selected for extreme breeding values for somatic cell scores, one line with high SCS (HSCS) and the other with low SCS (LSCS), were used to test the hypothesis that the mammary response and function differed between the lines. The aim of the present study was to investigate differences in the early immune response in caprine mammary gland tissues challenged with Staphylococcus aureus, one of the main pathogens responsible for the intra-mammary infection in small ruminants, using transcriptomic and histopathology analyses. The comparison between HSCS and LSCS goat lines, showed differences in the response at the histological level for inflammation, presence of neutrophils and micro-abscess formation, and at the molecular level in the expression of CXCL8, IL-6, NFKBIZ and IL-1β. CXCL8 and CXCL2 genes, which showed a higher level of expression in the experimentally infected HSCS line. The molecular data and histopathology both suggested that following S. aureus infection, mobilization, recruitment, infiltration, and chemotaxis of neutrophil, leads to a more severe inflammation in the HSCS compared to LSCS animals. Our results represent an initial basis for further studies to unravel the genetic basis of early mastitis inflammatory responses and the selection of dairy animals more resistant to bacterial mastitis.

Genome-Wide Epigenetic Characterization of Tissues from Three Germ Layers Isolated from Sheep Fetuses

DNA methylation of regulatory and growth-related genes contributes to fetal programming which is important for maintaining the correct development of three germ layers of the embryo that develope into different tissues and organs, and which persists into adult life. In this study, a preliminary epigenetic screen was performed to define genomic regions that are involved in fetal epigenome remodeling. Embryonic ectodermic tissues (origin of nervous tissue), mesenchymal tissues (origin of connective and muscular tissues), and foregut endoderm tissues (origin of epithelial tissue), from day 28 sheep fetuses were collected and the distribution of methylated CpGs was analyzed using whole-genome bisulfite sequencing. Patterns of methylation among the three tissues showed a high level of conservation of hypo-methylated CpG islands CGIs, and a consistent level of methylation in regulatory genetic elements. Analysis of tissue specific differentially methylated regions, revealed that 20% of the total CGIs differed between tissues. A proportion of the methylome was remodeled in gene bodies, 5′ UTRs and 3′ UTRs (7, 11, and 11%, respectively). Genes with overlapping differentially methylated regions in gene bodies and CGIs showed a significant enrichment for tissue morphogenesis and development pathways. The data presented here provides a “reference” for the epigenetic status of genes potentially involved in the maintenance and regulation of fetal developmental during early life, a period expected to be particularly prone to epigenetic alterations induced by environmental and nutritional stressors.

Short Communication: Maternal undernutrition during peri-conceptional period affects whole genome ovine muscle methylation in adult offspring

Experimental and epidemiological studies suggest that maternal nutritional status during early pregnancy, including the period around the time of conception, may induce long-lasting epigenetic changes in the offspring. However, this remains largely unexplored in livestock. Therefore, the objective of this study was to evaluate if modification of the maternal diet of sheep (CTR: control; UND: 50% undernutrition) during the peri-conceptional period (42 days in total: -14/+28 from mating), would impact CpG methylation in muscle tissue (Longissimus dorsi) of adult offspring (11.5 months old). Reduced Representation Bisulfite Sequencing (RRBS), identified 262 (Edge-R, FDR<0.05) and 686 (Logistic Regression, FDR <0.001) differentially methylated regions (DMRs) between the UND and CTR groups. Gene ontology (GO) analysis identified genes related to development, functions of the muscular system and steroid hormone receptor activity within the DMRs. The data reported here show that nutritional stress during early pregnancy leads to epigenetic modifications in the muscle of the resulting offspring, with possible implications for cardiac dysfunction, muscle physiology and meat production.

A Microenvironment-Related Nine-Gene Signature May Predict Survival in Mycosis Fungoides Patients at Diagnosis

Mycosis fungoides (MF) is the most common cutaneous lymphoma characterized by an indolent course. Prognosis is stage-based but this approach does not reflect the different outcomes within stages. Considering that tumor microenvironment is known to be involved in MF pathogenesis and progression, we decided to investigate 99 MF cases by using the PanCancer Immune Profiling Panel. We identified and validated a signature of 9 genes able to predict MF survival and distinguish a high-risk group with a worse outcome from a low-risk group of cases with a better outcome. At the molecular level, low-risk vs. high-risk cases reported a global upregulation of immune genes, enriched in cytokines, and a higher density of dendritic cells and mast cells, possibly associated with a more favorable clinical course.

Immune Signatures Identify Patient Subsets Deriving Long-Term Benefit From First-Line Rituximab in Follicular Lymphoma

Background

The role of first-line single-agent rituximab immunotherapy in follicular lymphoma (FL) remains debated, as most patients eventually undergo chemotherapy.

Methods

In this study, we retrospectively analyzed 81 FL patients treated with first-line single-agent rituximab monotherapy with (n = 53) or without (n = 28) consolidation. Fifty-one patients (63%) were high-tumor burden according to Group d'Etude des Lymphomes Folliculaires (GELF) criteria.

Results

After a median follow-up of 11 years, overall survival (OS) and progression-free survival (PFS) rates were 85% and 32%, respectively. Targeted gene expression profiling (T-GEP) was performed in 40 patients, revealing a 26-gene expression signature distinguishing complete responders and non-responders. This signature included genes involved in T-regulatory (Treg) and natural-killer cell activity, and interleukin-17 signaling. A simplified 14-gene prognostic score (ImSig) enabled accurate outcome stratification in terms of PFS. These data were validated in silico using two independent publicly available cohorts of FL patients treated with chemoimmunotherapy. Deconvolution analyses demonstrated an enrichment in Treg cells in high-risk ImSig patients, which was validated by immunohistochemistry.

Conclusions

These findings demonstrate that the efficacy of front-line anti-CD20 immunotherapy may depend on microenvironment-related factors, and that specific immune signatures could identify patient subsets obtaining long-term benefit from a chemo-free immunotherapeutic approach.

Trial Registration

The authors have confirmed clinical trial registration is not needed for this submission.

The Human EST Ontology Explorer: a tissue-oriented visualization system for ontologies distribution in human EST collections

BACKGROUND

The NCBI dbEST currently contains more than eight million human Expressed Sequenced Tags (ESTs). This wide collection represents an important source of information for gene expression studies, provided it can be inspected according to biologically relevant criteria. EST data can be browsed using different dedicated web resources, which allow to investigate library specific gene expression levels and to make comparisons among libraries, highlighting significant differences in gene expression. Nonetheless, no tool is available to examine distributions of quantitative EST collections in Gene Ontology (GO) categories, nor to retrieve information concerning library-dependent EST involvement in metabolic pathways. In this work we present the Human EST Ontology Explorer (HEOE) http://www.itb.cnr.it/ptp/human_est_explorer, a web facility for comparison of expression levels among libraries from several healthy and diseased tissues.

RESULTS

The HEOE provides library-dependent statistics on the distribution of sequences in the GO Direct Acyclic Graph (DAG) that can be browsed at each GO hierarchical level. The tool is based on large-scale BLAST annotation of EST sequences. Due to the huge number of input sequences, this BLAST analysis was performed with the aid of grid computing technology, which is particularly suitable to address data parallel task. Relying on the achieved annotation, library-specific distributions of ESTs in the GO Graph were inferred. A pathway-based search interface was also implemented, for a quick evaluation of the representation of libraries in metabolic pathways. EST processing steps were integrated in a semi-automatic procedure that relies on Perl scripts and stores results in a MySQL database. A PHP-based web interface offers the possibility to simultaneously visualize, retrieve and compare data from the different libraries. Statistically significant differences in GO categories among user selected libraries can also be computed.

CONCLUSION

The HEOE provides an alternative and complementary way to inspect EST expression levels with respect to approaches currently offered by other resources. Furthermore, BLAST computation on the whole human EST dataset was a suitable test of grid scalability in the context of large-scale bioinformatics analysis. The HEOE currently comprises sequence analysis from 70 non-normalized libraries, representing a comprehensive overview on healthy and unhealthy tissues. As the analysis procedure can be easily applied to other libraries, the number of represented tissues is intended to increase.

TOSCA: an automated Tumor Only Somatic CAlling workflow for somatic mutation detection without matched normal samples

Motivation

Accurate classification of somatic variants in a tumor sample is often accomplished by utilizing a paired normal tissue sample from the same patient to enable the separation of private germline mutations from somatic variants. However, a paired normal sample is not always available, making a reliable somatic variant calling more challenging. In silico screening of variants against public or private databases and other filtering approaches are often used in absence of a paired normal sample. Nevertheless, difficulties in performing a tumor-only calling with sufficient accuracy and lack of open-source software have limited their applications in clinical research.

Results

To address these limitations, we developed TOSCA, the first automated tumor-only somatic calling workflow in whole-exome sequencing and targeted panel sequencing data which performs an end-to-end analysis from raw read files, via quality checks, alignment and variant calling to functional annotation, databases filtering, tumor purity and ploidy estimation and variant classification. Application of our workflow to tumor-only data provides estimates of somatic and germline variants that are consistent with results from paired analyses.

Availability and implementation

TOSCA is a Snakemake-based workflow and freely available at https://github.com/mdelcorvo/TOSCA.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Tumor microenvironment of Burkitt lymphoma: different immune signatures with different clinical behavior

ABSTRACT

Burkitt lymphoma (BL) is characterized by a tumor microenvironment (TME) in which macrophages represent the main component, determining a distinct histological appearance known as "starry sky" pattern. However, in some instances, BL may exhibit a granulomatous reaction that has been previously linked to favorable prognosis and spontaneous regression. The aim of our study was to deeply characterize the immune landscape of 7 cases of Epstein-Barr virus-positive (EBV+) BL with granulomatous reaction compared with 8 cases of EBV+ BL and 8 EBV-negative (EBV-) BL, both with typical starry sky pattern, by Gene expression profiling performed on the NanoString nCounter platform. Subsequently, the data were validated using multiplex and combined immunostaining. Based on unsupervised clustering of differentially expressed genes, BL samples formed 3 distinct clusters differentially enriched in BL with a diffuse granulomatous reaction (cluster 1), EBV+ BL with typical starry sky pattern (cluster 2), EBV- BL with typical "starry sky" (cluster 3). We observed variations in the immune response signature among BL with granulomatous reaction and BL with typical "starry sky," both EBV+ and EBV-. The TME signature in BL with diffuse granulomatous reaction showed a proinflammatory response, whereas BLs with "starry sky" were characterized by upregulation of M2 polarization and protumor response. Moreover, the analysis of additional signatures revealed an upregulation of the dark zone signature and epigenetic signature in BL with a typical starry sky. Tumor-associated macrophages and epigenetic regulators may be promising targets for additional therapies for BL lymphoma, opening novel immunotherapeutic strategies.