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.

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.

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.

Identification of additional loci associated with antibody response to Mycobacterium avium ssp. Paratuberculosis in cattle by GSEA-SNP analysis.

Mycobacterium avium subsp. paratuberculosis: (MAP) causes a contagious chronic infection results in Johne's disease in a wide range of animal species, including cattle. Several genome-wide association studies (GWAS) have been carried out to identify loci putatively associated with MAP susceptibility by testing each marker separately and identifying SNPs that show a significant association with the phenotype, while SNP with modest effects are usually ignored. The objective of this study was to identify modest-effect genes associated with MAP susceptibility using a pathway-based approach. The Illumina BovineSNP50 BeadChip was used to genotype 966 Holstein cows, 483 positive and 483 negative for antibody response to MAP, data were then analyzed using novel SNP-based Gene Set Enrichment Analysis (GSEA-SNP) and validated with Adaptive Rank Truncated Product methodology. An allele-based test was carried out to estimate the statistical association for each marker with the phenotype, subsequently SNPs were mapped to the closest genes, considering for each gene the single variant with the highest value within a window of 50 kb, then pathway-statistics were tested using the GSEA-SNP method. The GO biological process "embryogenesis and morphogenesis" was most highly associated with antibody response to MAP. Within this pathway, five genes code for proteins which play a role in the immune defense relevant to response to bacterial infection. The immune response genes identified would not have been considered using a standard GWAS, thus demonstrating that the pathway approach can extend the interpretation of genome-wide association analyses and identify additional candidate genes for target traits.

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

Burkitt lymphoma (BL) is characterized by 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 a favorable prognosis and spontaneous regression. The aim of our study was to deeply characterize the immune landscape of 7 cases of EBV + BL with granulomatous reaction compared to 8 cases of EBV + BL and 8 EBV- BL, both with typical "starry sky" pattern, by Gene expression profiling performed on the NanoString nCounter platform. Subsequently, the data were validated by 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, while BLs with "starry sky" were characterized by up-regulation of M2- polarization and pro-tumor response. Moreover, the analysis of additional signatures revealed an up-regulation of Dark zone-signature and epigenetic-signature in BL with typical "starry sky". Tumor associated macrophages (TAM) and epigenetic regulators may be promising targets for additional therapies in BL lymphoma opening novel immunotherapeutic strategies.

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.

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.

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.

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.

Maternal undernutrition during periconceptional 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 periconceptional period (42 d in total: -14/+28 from mating), would impact CpG methylation in muscle tissue (Longissimus dorsi) of adult offspring (11.5 mo old). Reduced representation bisulfite sequencing 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 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.

The formation of the epigenetic pattern of an organism is highly sensitive to environmental factors, especially during early mammalian development, when epigenetic reprogramming of the whole genome takes place. In utero adverse conditions experienced during early pregnancy, such as maternal undernutrition, may induce long-lasting epigenetic changes in the resulting offspring. This study investigated the CpG methylation variations in muscle tissue of adult offspring induced by differences in the diet of their mothers during pregnancy. Our data show that undernutrition during pregnancy leads to epigenetic alterations in the muscle of the offspring, with a potential impact on animal health and productivity.

The identification of TCF1+ progenitor exhausted T cells in THRLBCL may predict a better response to PD-1/PD-L1 blockade.

T-cell/histiocyte-rich large B-cell lymphoma (THRLBCL) is a rare and aggressive variant of diffuse large B-cell lymphoma (DLBCL) that usually affects young to middle-aged patients, with disseminated disease at presentation. The tumor microenvironment (TME) plays a key role in THRLBCL due to its peculiar cellular composition (<10% neoplastic B cells interspersed in a cytotoxic T-cell/histiocyte-rich background). A significant percentage of THRLBCL is refractory to rituximab, cyclophosphamide, doxorubicin, vincristine, and prednisone (RCHOP)-based regimens and to chimeric antigen receptor T-cell therapy; thus, the development of a specific therapeutic approach for these patients represents an unmet clinical need. To better understand the interaction of immune cells in THRLBCL TME and identify more promising therapeutic strategies, we compared the immune gene expression profiles of 12 THRLBCL and 10 DLBCL samples, and further corroborated our findings in an extended in silico set. Gene coexpression network analysis identified the predominant role of the programmed cell death protein 1 (PD-1)/programmed cell death ligand 1 (PD-L1) axis in the modulation of the immune response. Furthermore, the PD-1/PD-L1 activation was flanked by the overexpression of 48 genes related to the functional exhaustion of T cells. Globally, THRLBCL TME was highly interferon-inflamed and severely exhausted. The immune gene profiling findings strongly suggest that THRLBCL may be responsive to anti-PD-1 therapy but also allowed us to take a step forward in understanding THRLBCL TME. Of therapeutic relevance, we validated our results by immunohistochemistry, identifying a subset of TCF1+ (T cell-specific transcription factor 1, encoded by the TCF7 gene) progenitor exhausted T cells enriched in patients with THRLBCL. This subset of TCF1+ exhausted T cells correlates with good clinical response to immune checkpoint therapy and may improve prediction of anti-PD-1 response in patients with THRLBCL.

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.

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.

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.

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.

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

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.

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.

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.

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.

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.