MeSH ORA framework: R/Bioconductor packages to support MeSH over-representation analysis

Comparative analysis of machine learning approaches for predicting respiratory virus infection and symptom severity

Respiratory diseases are among the major health problems causing a burden on hospitals. Diagnosis of infection and rapid prediction of severity without time-consuming clinical tests could be beneficial in preventing the spread and progression of the disease, especially in countries where health systems remain incapable. Personalized medicine studies involving statistics and computer technologies could help to address this need. In addition to individual studies, competitions are also held such as Dialogue for Reverse Engineering Assessment and Methods (DREAM) challenge which is a community-driven organization with a mission to research biology, bioinformatics, and biomedicine. One of these competitions was the Respiratory Viral DREAM Challenge, which aimed to develop early predictive biomarkers for respiratory virus infections. These efforts are promising, however, the prediction performance of the computational methods developed for detecting respiratory diseases still has room for improvement. In this study, we focused on improving the performance of predicting the infection and symptom severity of individuals infected with various respiratory viruses using gene expression data collected before and after exposure. The publicly available gene expression dataset in the Gene Expression Omnibus, named GSE73072, containing samples exposed to four respiratory viruses (H1N1, H3N2, human rhinovirus (HRV), and respiratory syncytial virus (RSV)) was used as input data. Various preprocessing methods and machine learning algorithms were implemented and compared to achieve the best prediction performance. The experimental results showed that the proposed approaches obtained a prediction performance of 0.9746 area under the precision-recall curve (AUPRC) for infection (i.e., shedding) prediction (SC-1), 0.9182 AUPRC for symptom class prediction (SC-2), and 0.6733 Pearson correlation for symptom score prediction (SC-3) by outperforming the best leaderboard scores of Respiratory Viral DREAM Challenge (a 4.48% improvement for SC-1, a 13.68% improvement for SC-2, and a 13.98% improvement for SC-3). Additionally, over-representation analysis (ORA), which is a statistical method for objectively determining whether certain genes are more prevalent in pre-defined sets such as pathways, was applied using the most significant genes selected by feature selection methods. The results show that pathways associated with the 'adaptive immune system' and 'immune disease' are strongly linked to pre-infection and symptom development. These findings contribute to our knowledge about predicting respiratory infections and are expected to facilitate the development of future studies that concentrate on predicting not only infections but also the associated symptoms.

Impact of aging on meningeal gene expression

BACKGROUND

The three-layered meninges cover and protect the central nervous system and form the interface between cerebrospinal fluid and the brain. They are host to a lymphatic system essential for maintaining fluid dynamics inside the cerebrospinal fluid-filled subarachnoid space and across the brain parenchyma via their connection to glymphatic structures. Meningeal fibroblasts lining and traversing the subarachnoid space have direct impact on the composition of the cerebrospinal fluid through endocytotic uptake as well as extensive protein secretion. In addition, the meninges are an active site for immunological processes and act as gatekeeper for immune cells entering the brain. During aging in mice, lymphatic drainage from the brain is less efficient contributing to neurodegenerative processes. Aging also affects the immunological status of the meninges, with increasing numbers of T cells, changing B cell make-up, and altered macrophage complement.

METHODS

We employed RNASeq to measure gene expression and to identify differentially expressed genes in meninges isolated from young and aged mice. Using Ingenuity pathway, GO term, and MeSH analyses, we identified regulatory pathways and cellular functions in meninges affected by aging.

RESULTS

Aging had profound impact on meningeal gene expression. Pathways related to innate as well as adaptive immunity were affected. We found evidence for increasing numbers of T and B lymphocytes and altered activity profiles for macrophages and other myeloid cells. Furthermore, expression of pro-inflammatory cytokine and chemokine genes increased with aging. Similarly, the complement system seemed to be more active in meninges of aged mice. Altered expression of solute carrier genes pointed to age-dependent changes in cerebrospinal fluid composition. In addition, gene expression for secreted proteins showed age-dependent changes, in particular, genes related to extracellular matrix composition and organization were affected.

CONCLUSIONS

Aging has profound effects on meningeal gene expression; thereby affecting the multifaceted functions meninges perform to maintain the homeostasis of the central nervous system. Thus, age-dependent neurodegenerative processes and cognitive decline are potentially in part driven by altered meningeal function.

DSEATM: drug set enrichment analysis uncovering disease mechanisms by biomedical text mining

Disease pathogenesis is always a major topic in biomedical research. With the exponential growth of biomedical information, drug effect analysis for specific phenotypes has shown great promise in uncovering disease-associated pathways. However, this method has only been applied to a limited number of drugs. Here, we extracted the data of 4634 diseases, 3671 drugs, 112 809 disease-drug associations and 81 527 drug-gene associations by text mining of 29 168 919 publications. On this basis, we proposed a 'Drug Set Enrichment Analysis by Text Mining (DSEATM)' pipeline and applied it to 3250 diseases, which outperformed the state-of-the-art method. Furthermore, diseases pathways enriched by DSEATM were similar to those obtained using the TCGA cancer RNA-seq differentially expressed genes. In addition, the drug number, which showed a remarkable positive correlation of 0.73 with the AUC, plays a determining role in the performance of DSEATM. Taken together, DSEATM is an auspicious and accurate disease research tool that offers fresh insights.

Transcriptomic Analysis of Fumarate Compounds Identifies Unique Effects of Isosorbide Di-(Methyl Fumarate) on NRF2, NF-kappaB and IRF1 Pathway Genes

Dimethyl fumarate (DMF) has emerged as a first-line therapy for relapsing-remitting multiple sclerosis (RRMS). This treatment, however, has been limited by adverse effects, which has prompted development of novel derivatives with improved tolerability. We compared the effects of fumarates on gene expression in astrocytes. Our analysis included diroximel fumarate (DRF) and its metabolite monomethyl fumarate (MMF), along with a novel compound isosorbide di-(methyl fumarate) (IDMF). Treatment with IDMF resulted in the largest number of differentially expressed genes. The effects of DRF and MMF were consistent with NRF2 activation and NF-κB inhibition, respectively. IDMF responses, however, were concordant with both NRF2 activation and NF-κB inhibition, and we confirmed IDMF-mediated NF-κB inhibition using a reporter assay. IDMF also down-regulated IRF1 expression and IDMF-decreased gene promoters were enriched with IRF1 recognition sequences. Genes altered by each fumarate overlapped significantly with those near loci from MS genetic association studies, but IDMF had the strongest overall effect on MS-associated genes. These results show that next-generation fumarates, such as DRF and IDMF, have effects differing from those of the MMF metabolite. Our findings support a model in which IDMF attenuates oxidative stress via NRF2 activation, with suppression of NF-κB and IRF1 contributing to mitigation of inflammation and pyroptosis.

The evolution of knowledge on genes associated with human diseases

Thousands of biomedical scientific articles, including those describing genes associated with human diseases, are published every week. Computational methods such as text mining and machine learning algorithms are now able to automatically detect these associations. In this study, we used a cognitive computing text-mining application to construct a knowledge network comprising 3,723 genes and 99 diseases. We then tracked the yearly changes on these networks to analyze how our knowledge has evolved in the past 30 years. Our systems approach helped to unravel the molecular bases of diseases and detect shared mechanisms between clinically distinct diseases. It also revealed that multi-purpose therapeutic drugs target genes that are commonly associated with several psychiatric, inflammatory, or infectious disorders. By navigating this knowledge tsunami, we were able to extract relevant biological information and insights about human diseases.

Genome-wide association study of resistance/susceptibility to infectious bovine keratoconjunctivitis in Brazilian Hereford cattle

Genome-wide association studies were conducted to identify the more informative genomic regions and SNPs, as well as to identify candidate genes associated with infectious bovine keratoconjunctivitis (IBK) resistance/susceptibility in Hereford cattle. A Bayes B statistical approach was initially applied in genome-wide association studies by using deregressed estimated breeding values for IBK resistance/susceptibility. To estimate the combined effect of a genomic region that is potentially associated with QTL, 2504 non-overlapping 1-Mb windows that varied in SNP number were defined, with the most informative 24 windows including 427 SNPs and explaining more than 20% of the estimated genetic variance for IBK resistance/susceptibility. These regions were explored with respect to their biological functions through functional analysis to map potential candidate genes. The significant SNPs were mapped on chromosomes 1, 3, 5, 6, 7, 8, 10, 11, 12, 13, 14, 15, 18, 20, 23, and 28, and candidate genes were detected as related to the IBK. Most informative SNPs in term of genetic variance were located in proximity of genes related to phenotypic expression of lesions and biological processes associated to the IBK. Knowledge about phenotypic and genomic variation generated in the present study can be used to on design selection strategies to improve the resistance to IBK of Hereford cattle herds.

Genome-wide association and genomic prediction for scrotal circumference in Hereford and Braford bulls

Scrotal circumference (SC) is widely used as a selection criterion for bulls in breeding programs, since it is easily assessed and correlated with several desirable reproductive traits. The objectives of this study were: to perform a genome-wide association study (GWAS) to identify genomic regions associated with SC adjusted for age (SCa) and for both age and weight (SCaw); to select Tag SNPs from GWAS to construct low-density panel for genomic prediction; and to compare the prediction accuracy of the SC through different methods for Braford and Hereford bulls from the same genetic breeding program. Data of SC from 18,172 bulls (30.4 ± 3.7 cm) and of genotypes from 131 sires and 3,545 animals were used. From GWAS, the top 1% of 1-Mb windows were observed on chromosome (BTA) 2, 20, 7, 8, 15, 3, 16, 27, 6 and 8 for SCa and on BTA 8, 15, 16, 21, 19, 2, 6, 5 and 10 for SCaw, representing 17.4% and 18.8% of the additive genetic variance of SCa and SCaw, respectively. The MeSH analysis was able to translate genomic information providing biological meanings of more specific gene functions related to the SCa and SCaw. The genomic enhancement methods, especially single step GBLUP, that combined phenotype and pedigree data with direct genomic values generated gains in accuracy in relation to pedigree BLUP, suggesting that genomic predictions should be applied to improve genetic gain and to narrow the generation interval compared to traditional methods. The proposed Tag-SNP panels may be useful for lower-cost commercial genomic prediction applications in the future, when the number of bulls in the reference population increases for SC in Hereford and Braford breeds.

Chamber-enriched gene expression profiles in failing human hearts with reduced ejection fraction

Heart failure with reduced ejection fraction (HFrEF) constitutes 50% of HF hospitalizations and is characterized by high rates of mortality. To explore the underlying mechanisms of HFrEF etiology and progression, we studied the molecular and cellular differences in four chambers of non-failing (NF, n = 10) and HFrEF (n = 12) human hearts. We identified 333 genes enriched within NF heart subregions and often associated with cardiovascular disease GWAS variants. Expression analysis of HFrEF tissues revealed extensive disease-associated transcriptional and signaling alterations in left atrium (LA) and left ventricle (LV). Common left heart HFrEF pathologies included mitochondrial dysfunction, cardiac hypertrophy and fibrosis. Oxidative stress and cardiac necrosis pathways were prominent within LV, whereas TGF-beta signaling was evident within LA. Cell type composition was estimated by deconvolution and revealed that HFrEF samples had smaller percentage of cardiomyocytes within the left heart, higher representation of fibroblasts within LA and perivascular cells within the left heart relative to NF samples. We identified essential modules associated with HFrEF pathology and linked transcriptome discoveries with human genetics findings. This study contributes to a growing body of knowledge describing chamber-specific transcriptomics and revealed genes and pathways that are associated with heart failure pathophysiology, which may aid in therapeutic target discovery.

Discovering lethal alleles across the turkey genome using a transmission ratio distortion approach

Deviation from Mendelian inheritance expectations (transmission ratio distortion, TRD) has been observed in several species, including the mouse and humans. In this study, TRD was characterized in the turkey genome using both allelic (specific- and unspecific-parent TRD) and genotypic (additive- and dominance-TRD) parameterizations within a Bayesian framework. In this study, we evaluated TRD for 23 243 genotyped Turkeys across 56 393 autosomal SNPs. The analyses included 500 sires, 2013 dams and 11 047 offspring (trios). Three different haplotype sliding windows of 4, 10 and 20 SNPs were used across the autosomal chromosomes. Based on the genotypic parameterizations, 14 haplotypes showed additive and dominance TRD effects highlighting regions with a recessive TRD pattern. In contrast, the allelic model uncovered 12 haplotype alleles with the allelic TRD pattern which showed an underrepresentation of heterozygous offspring in addition to the absence of homozygous animals. For regions with the allelic pattern, only one particular region showed a parent-specific TRD where the penetrance was high via the dam, but low via the sire. The gene set analysis uncovered several gene ontology functional terms, Reactome pathways and several Medical Subject Headings that showed significant enrichment of genes associated with TRD. Many of these gene ontology functional terms (e.g. mitotic spindle assembly checkpoint, DRM complex and Aneuploidy), Reactome pathways (e.g. Mismatch repair) and Medical Subject Headings (e.g. Adenosine monophosphate) are known to be related to fertility, embryo development and lethality. The results of this study revealed potential novel candidate lethal haplotypes, functional terms and pathways that may enhance breeding programs in Turkeys through reducing mortality and improving reproduction rate.

pyMeSHSim: an integrative python package for biomedical named entity recognition, normalization, and comparison of MeSH terms

Background

Many disease causing genes have been identified through different methods, but there have been no uniform annotations of biomedical named entity (bio-NE) of the disease phenotypes of these genes yet. Furthermore, semantic similarity comparison between two bio-NE annotations has become important for data integration or system genetics analysis.

Results

The package pyMeSHSim recognizes bio-NEs by using MetaMap which produces Unified Medical Language System (UMLS) concepts in natural language process. To map the UMLS concepts to Medical Subject Headings (MeSH), pyMeSHSim is embedded with a house-made dataset containing the main headings (MHs), supplementary concept records (SCRs), and their relations in MeSH. Based on the dataset, pyMeSHSim implemented four information content (IC)-based algorithms and one graph-based algorithm to measure the semantic similarity between two MeSH terms. To evaluate its performance, we used pyMeSHSim to parse OMIM and GWAS phenotypes. The pyMeSHSim introduced SCRs and the curation strategy of non-MeSH-synonymous UMLS concepts, which improved the performance of pyMeSHSim in the recognition of OMIM phenotypes. In the curation of 461 GWAS phenotypes, pyMeSHSim showed recall > 0.94, precision > 0.56, and F1 > 0.70, demonstrating better performance than the state-of-the-art tools DNorm and TaggerOne in recognizing MeSH terms from short biomedical phrases. The semantic similarity in MeSH terms recognized by pyMeSHSim and the previous manual work was calculated by pyMeSHSim and another semantic analysis tool meshes, respectively. The result indicated that the correlation of semantic similarity analysed by two tools reached as high as 0.89–0.99.

Conclusions

The integrative MeSH tool pyMeSHSim embedded with the MeSH MHs and SCRs realized the bio-NE recognition, normalization, and comparison in biomedical text-mining.

Fine mapping of genomic regions associated with female fertility in Nellore beef cattle based on sequence variants from segregating sires

Background

Impaired fertility in cattle limits the efficiency of livestock production systems. Unraveling the genetic architecture of fertility traits would facilitate their improvement by selection. In this study, we characterized SNP chip haplotypes at QTL blocks then used whole-genome sequencing to fine map genomic regions associated with reproduction in a population of Nellore (Bos indicus) heifers.

Methods

The dataset comprised of 1337 heifers genotyped using a GeneSeek® Genomic Profiler panel (74677 SNPs), representing the daughters from 78 sires. After performing marker quality control, 64800 SNPs were retained. Haplotypes carried by each sire at six previously identified QTL on BTAs 5, 14 and 18 for heifer pregnancy and BTAs 8, 11 and 22 for antral follicle count were constructed using findhap software. The significance of the contrasts between the effects of every two paternally-inherited haplotype alleles were used to identify sires that were heterozygous at each QTL. Whole-genome sequencing data localized to the haplotypes from six sires and 20 other ancestors were used to identify sequence variants that were concordant with the haplotype contrasts. Enrichment analyses were applied to these variants using KEGG and MeSH libraries.

Results

A total of six (BTA 5), six (BTA 14) and five (BTA 18) sires were heterozygous for heifer pregnancy QTL whereas six (BTA 8), fourteen (BTA 11), and five (BTA 22) sires were heterozygous for number of antral follicles’ QTL. Due to inadequate representation of many haplotype alleles in the sequenced animals, fine mapping analysis could only be reliably performed for the QTL on BTA 5 and 14, which had 641 and 3733 concordant candidate sequence variants, respectively. The KEGG “Circadian rhythm” and “Neurotrophin signaling pathway” were significantly associated with the genes in the QTL on BTA 5 whereas 32 MeSH terms were associated with the QTL on BTA 14. Among the concordant sequence variants, 0.2% and 0.3% were classified as missense variants for BTAs 5 and 14, respectively, highlighting the genes MTERF2, RTMB, ENSBTAG00000037306 (miRNA), ENSBTAG00000040351, PRKDC, and RGS20. The potential causal mutations found in the present study were associated with biological processes such as oocyte maturation, embryo development, placenta development and response to reproductive hormones.

Conclusions

The identification of heterozygous sires by positionally phasing SNP chip data and contrasting haplotype effects for previously detected QTL can be used for fine mapping to identify potential causal mutations and candidate genes. Genomic variants on genes MTERF2, RTBC, miRNA ENSBTAG00000037306, ENSBTAG00000040351, PRKDC, and RGS20, which are known to have influence on reproductive biological processes, were detected.

Single-step genome-wide association for longitudinal traits of Canadian Ayrshire, Holstein, and Jersey dairy cattle

Estimating single nucleotide polymorphism (SNP) effects over time is essential to identify and validate candidate genes (or quantitative trait loci) associated with time-dependent variation of economically important traits and to better understand the underlying mechanisms of lactation biology. Therefore, in this study, we aimed to estimate time-dependent effects of SNP and identifying candidate genes associated with milk (MY), fat (FY), and protein (PY) yields, and somatic cell score (SCS) in the first 3 lactations of Canadian Ayrshire, Holstein, and Jersey breeds, as well as suggest their potential pattern of phenotypic effect over time. Random regression coefficients for the additive direct genetic effect were estimated for each animal using single-step genomic BLUP, based on 2 random regression models: one considering MY, FY, and PY in the first 3 lactations and the other considering SCS in the first 3 lactations. Thereafter, SNP solutions were obtained for random regression coefficients, which were used to estimate the SNP effects over time (from 5 to 305 d in lactation). The top 1% of SNP that showed a high magnitude of SNP effect in at least 1 d in lactation were selected as relevant SNP for further analyses of candidate genes, and clustered according to the trajectory of their SNP effects over time. The majority of SNP selected for MY, FY, and PY increased the magnitude of their effects over time, for all breeds. In contrast, for SCS, most selected SNP decreased the magnitude of their effects over time, especially for the Holstein and Jersey breeds. In general, we identified a different set of candidate genes for each breed, and similar genes were found across different lactations for the same trait in the same breed. For some of the candidate genes, the suggested pattern of phenotypic effect changed among lactations. Among the lactations, candidate genes (and their suggested phenotypic effect over time) identified for the second and third lactations were more similar to each other than for the first lactation. Well-known candidate genes with major effects on milk production traits presented different suggested patterns of phenotypic effect across breeds, traits, and lactations in which they were identified. The candidate genes identified in this study can be used as target genes in studies of gene expression.

Genome-wide association for milk production traits and somatic cell score in different lactation stages of Ayrshire, Holstein, and Jersey dairy cattle

We performed genome-wide association analyses for milk, fat, and protein yields and somatic cell score based on lactation stages in the first 3 parities of Canadian Ayrshire, Holstein, and Jersey cattle. The genome-wide association analyses were performed considering 3 different lactation stages for each trait and parity: from 5 to 95, from 96 to 215, and from 216 to 305 d in milk. Effects of single nucleotide polymorphisms (SNP) for each lactation stage, trait, parity, and breed were estimated by back-solving the direct breeding values estimated using the genomic best linear unbiased predictor and single-trait random regression test-day models containing only the fixed population average curve and the random genomic curves. To identify important genomic regions related to the analyzed lactation stages, traits, parities and breeds, moving windows (SNP-by-SNP) of 20 adjacent SNP explaining more than 0.30% of total genetic variance were selected for further analyses of candidate genes. A lower number of genomic windows with a relatively higher proportion of the explained genetic variance was found in the Holstein breed compared with the Ayrshire and Jersey breeds. Genomic regions associated with the analyzed traits were located on 12, 8, and 15 chromosomes for the Ayrshire, Holstein, and Jersey breeds, respectively. Especially for the Holstein breed, many of the identified candidate genes supported previous reports in the literature. However, well-known genes with major effects on milk production traits (e.g., diacylglycerol O-acyltransferase 1) showed contrasting results among lactation stages, traits, and parities of different breeds. Therefore, our results suggest evidence of differential sets of candidate genes underlying the phenotypic expression of the analyzed traits across breeds, parities, and lactation stages. Further functional studies are needed to validate our findings in independent populations.

GM604 regulates developmental neurogenesis pathways and the expression of genes associated with amyotrophic lateral sclerosis

Background

Amyotrophic lateral sclerosis (ALS) is currently an incurable disease without highly effective pharmacological treatments. The peptide drug GM604 (GM6 or Alirinetide) was developed as a candidate ALS therapy, which has demonstrated safety and good drug-like properties with a favorable pharmacokinetic profile. GM6 is hypothesized to bolster neuron survival through the multi-target regulation of developmental pathways, but mechanisms of action are not fully understood.

Methods

This study used RNA-seq to evaluate transcriptome responses in SH-SY5Y neuroblastoma cells following GM6 treatment (6, 24 and 48 h).

Results

We identified 2867 protein-coding genes with expression significantly altered by GM6 (FDR < 0.10). Early (6 h) responses included up-regulation of Notch and hedgehog signaling components, with increased expression of developmental genes mediating neurogenesis and axon growth. Prolonged GM6 treatment (24 and 48 h) altered the expression of genes contributing to cell adhesion and the extracellular matrix. GM6 further down-regulated the expression of genes associated with mitochondria, inflammatory responses, mRNA processing and chromatin organization. GM6-increased genes were located near GC-rich motifs interacting with C2H2 zinc finger transcription factors, whereas GM6-decreased genes were located near AT-rich motifs associated with helix-turn-helix homeodomain factors. Such motifs interacted with a diverse network of transcription factors encoded by GM6-regulated genes (STAT3, HOXD11, HES7, GLI1). We identified 77 ALS-associated genes with expression significantly altered by GM6 treatment (FDR < 0.10), which were known to function in neurogenesis, axon guidance and the intrinsic apoptosis pathway.

Conclusions

Our findings support the hypothesis that GM6 acts through developmental-stage pathways to influence neuron survival. Gene expression responses were consistent with neurotrophic effects, ECM modulation, and activation of the Notch and hedgehog neurodevelopmental pathways. This multifaceted mechanism of action is unique among existing ALS drug candidates and may be applicable to multiple neurodegenerative diseases.

Electronic supplementary material

The online version of this article (10.1186/s40035-018-0135-7) contains supplementary material, which is available to authorized users.

Integrating bioinformatics approaches for a comprehensive interpretation of metabolomics datasets

Access to high quality metabolomics data has become a routine component for biological studies. However, interpreting those datasets in biological contexts remains a challenge, especially because many identified metabolites are not found in biochemical pathway databases. Starting from statistical analyses, a range of new tools are available, including metabolite set enrichment analysis, pathway and network visualization, pathway prediction, biochemical databases and text mining. Integrating these approaches into comprehensive and unbiased interpretations must carefully consider both caveats of the metabolomics dataset itself as well as the structure and properties of the biological study design. Special considerations need to be taken when adopting approaches from genomics for use in metabolomics. R and Python programming language are enabling an easier exchange of diverse tools to deploy integrated workflows. This review summarizes the key ideas and latest developments in regards to these approaches.

Genetic dissection of bull fertility in US Jersey dairy cattle

The service sire has been recognized as an important factor affecting herd fertility in dairy cattle. Recent studies suggest that genetic factors explain part of the difference in fertility among Holstein sires. The main objective of this study was to dissect the genetic architecture of sire fertility in US Jersey cattle. The dataset included 1.5 K Jersey bulls with sire conception rate (SCR) records and 96 K single nucleotide polymorphism (SNP) markers spanning the whole genome. The analysis included whole‐genome scans for both additive and non‐additive effects and subsequent functional enrichment analyses using KEGG Pathway, Gene Ontology (GO) and Medical Subject Headings (MeSH) databases. Ten genomic regions located on eight different chromosomes explained more than 0.5% of the additive genetic variance for SCR. These regions harbor genes, such as PKDREJ,EPB41L2,PDGFD,STX2,SLC25A20 and IP6K1, that are directly implicated in testis development and spermatogenesis, sperm motility and the acrosome reaction. In addition, the genomic scan for non‐additive effects identified two regions on BTA11 and BTA25 with marked recessive effects. These regions harbor three genes—FER1L5,CNNM4 and DNAH3—with known roles in sperm biology. Moreover, the gene‐set analysis revealed terms associated with calcium regulation and signaling, membrane fusion, sperm cell energy metabolism, GTPase activity and MAPK signaling. These gene sets are directly implicated in sperm physiology and male fertility. Overall, this integrative genomic study unravels genetic variants and pathways affecting Jersey bull fertility. These findings may contribute to the development of novel genomic strategies for improving sire fertility in Jersey cattle.

RNA-Seq reveals novel genes and pathways involved in bovine mammary involution during the dry period and under environmental heat stress

The bovine dry period is a dynamic non-lactating phase where the mammary gland undergoes extensive cellular turnover. Utilizing RNA sequencing, we characterized novel genes and pathways involved in this process and determined the impact of dry period heat stress. Mammary tissue was collected before and during the dry period (−3, 3, 7, 14, and 25 days relative to dry-off [day 0]) from heat-stressed (HT, n = 6) or cooled (CL, n = 6) late-gestation Holstein cows. We identified 3,315 differentially expressed genes (DEGs) between late lactation and early involution, and 880 DEGs later in the involution process. DEGs, pathways, and upstream regulators during early involution support the downregulation of functions such as anabolism and milk component synthesis, and upregulation of cell death, cytoskeleton degradation, and immune response. The impact of environmental heat stress was less significant, yet genes, pathways, and upstream regulators involved in processes such as ductal branching morphogenesis, cell death, immune function, and protection against tissue stress were identified. Our research advances understanding of the mammary gland transcriptome during the dry period, and under heat stress insult. Individual genes, pathways, and upstream regulators highlighted in this study point towards potential targets for dry period manipulation and mitigation of the negative consequences of heat stress on mammary function.

RNA-Seq Analysis of IL-1B and IL-36 Responses in Epidermal Keratinocytes Identifies a Shared MyD88-Dependent Gene Signature

IL-36 cytokines have recently emerged as mediators of inflammation in autoimmune conditions including psoriasis vulgaris (PsV) and generalized pustular psoriasis (GPP). This study used RNA-seq to profile the transcriptome of primary epidermal keratinocytes (KCs) treated with IL-1B, IL-36A, IL-36B, or IL-36G. We identified some early IL-1B-specific responses (8 h posttreatment), but nearly all late IL-1B responses were replicated by IL-36 cytokines (24 h posttreatment). Type I and II interferon genes exhibited time-dependent response patterns, with early induction (8 h) followed by no response or repression (24 h). Altogether, we identified 225 differentially expressed genes (DEGs) with shared responses to all 4 cytokines at both time points (8 and 24 h). These involved upregulation of ligands (IL1A, IL1B, and IL36G) and activating proteases (CTSS) but also upregulation of inhibitors such as IL1RN and IL36RN. Shared IL-1B/IL-36 DEGs overlapped significantly with genes altered in PsV and GPP skin lesions, as well as genes near GWAS loci linked to autoimmune and autoinflammatory diseases (e.g., PsV, psoriatic arthritis, inflammatory bowel disease, and primary biliary cholangitis). Inactivation of MyD88 adapter protein using CRISPR/Cas9 completely abolished expression responses of such DEGs to IL-1B and IL-36G stimulation. These results provide a global view of IL-1B and IL-36 expression responses in epidermal KCs with fine-scale characterization of time-dependent and cytokine-specific response patterns. Our findings support an important role for IL-1B and IL-36 in autoimmune or autoinflammatory conditions and show that MyD88 adaptor protein mediates shared IL-1B/IL-36 responses.

Genomic study and Medical Subject Headings enrichment analysis of early pregnancy rate and antral follicle numbers in Nelore heifers

Zebu animals () are known to take longer to reach puberty compared with taurine animals (), limiting the supply of animals for harvest or breeding and impacting profitability. Genomic information can be a helpful tool to better understand complex traits and improve genetic gains. In this study, we performed a genomewide association study (GWAS) to identify genetic variants associated with reproductive traits in Nelore beef cattle. Heifer pregnancy (HP) was recorded for 1,267 genotyped animals distributed in 12 contemporary groups (CG) with an average pregnancy rate of 0.35 (±0.01). Disregarding one of these CG, the number of antral follicles (NF) was also collected for 937 of these animals, with an average of 11.53 (±4.43). The animals were organized in CG: 12 and 11 for HP and NF, respectively. Genes in linkage disequilibrium (LD) with the associated variants can be considered in a functional enrichment analysis to identify biological mechanisms involved in fertility. Medical Subject Headings (MeSH) were detected using the MESHR package, allowing the extraction of broad meanings from the gene lists provided by the GWAS. The estimated heritability for HP was 0.28 ± 0.07 and for NF was 0.49 ± 0.09, with the genomic correlation being -0.21 ± 0.29. The average LD between adjacent markers was 0.23 ± 0.01, and GWAS identified genomic windows that accounted for >1% of total genetic variance on chromosomes 5, 14, and 18 for HP and on chromosomes 2, 8, 11, 14, 15, 16, and 22 for NF. The MeSH enrichment analyses revealed significant ( < 0.05) terms associated with HP-"Munc18 Proteins," "Fucose," and "Hemoglobins"-and with NF-"Cathepsin B," "Receptors, Neuropeptide," and "Palmitic Acid." This is the first study in Nelore cattle introducing the concept of MeSH analysis. The genomic analyses contributed to a better understanding of the genetic control of the reproductive traits HP and NF and provide new selection strategies to improve beef production.

Using scale and feather traits for module construction provides a functional approach to chicken epidermal development

Gene co-expression network analysis has been a research method widely used in systematically exploring gene function and interaction. Using the Weighted Gene Co-expression Network Analysis (WGCNA) approach to construct a gene co-expression network using data from a customized 44K microarray transcriptome of chicken epidermal embryogenesis, we have identified two distinct modules that are highly correlated with scale or feather development traits. Signaling pathways related to feather development were enriched in the traditional KEGG pathway analysis and functional terms relating specifically to embryonic epidermal development were also enriched in the Gene Ontology analysis. Significant enrichment annotations were discovered from customized enrichment tools such as Modular Single-Set Enrichment Test (MSET) and Medical Subject Headings (MeSH). Hub genes in both trait-correlated modules showed strong specific functional enrichment toward epidermal development. Also, regulatory elements, such as transcription factors and miRNAs, were targeted in the significant enrichment result. This work highlights the advantage of this methodology for functional prediction of genes not previously associated with scale- and feather trait-related modules.

Medical Subject Heading (MeSH) annotations illuminate maize genetics and evolution

BACKGROUND

High-density marker panels and/or whole-genome sequencing, coupled with advanced phenotyping pipelines and sophisticated statistical methods, have dramatically increased our ability to generate lists of candidate genes or regions that are putatively associated with phenotypes or processes of interest. However, the speed with which we can validate genes, or even make reasonable biological interpretations about the principles underlying them, has not kept pace. A promising approach that runs parallel to explicitly validating individual genes is analyzing a set of genes together and assessing the biological similarities among them. This is often achieved via gene ontology analysis, a powerful tool that involves evaluating publicly available gene annotations. However, additional resources such as Medical Subject Headings (MeSH) can also be used to evaluate sets of genes to make biological interpretations.

RESULTS

In this manuscript, we describe utilizing MeSH terms to make biological interpretations in maize. MeSH terms are assigned to PubMed-indexed manuscripts by the National Library of Medicine, and can be directly mapped to genes to develop gene annotations. Once mapped, these terms can be evaluated for enrichment in sets of genes or similarity between gene sets to provide biological insights. Here, we implement MeSH analyses in five maize datasets to demonstrate how MeSH can be leveraged by the maize and broader crop-genomics community.

CONCLUSIONS

We demonstrate that MeSH terms can be effectively leveraged to generate hypotheses and make biological interpretations in maize, and we provide a pipeline that enables the use of MeSH terms in other plant species.

Unravelling the genomic architecture of bull fertility in Holstein cattle

BACKGROUND

Fertility is considered an important economic trait in dairy cattle. Most studies have investigated cow fertility while bull fertility has received much less consideration. The main objective of this study was to perform a comprehensive genomic analysis in order to unravel the genomic architecture underlying sire fertility in Holstein dairy cattle. The analysis included the application of alternative genome-wide association mapping approaches and the subsequent use of diverse gene set enrichment tools.

RESULTS

The association analyses identified at least eight genomic regions strongly associated with bull fertility. Most of these regions harbor genes, such as KAT8, CKB, TDRD9 and IGF1R, with functions related to sperm biology, including sperm development, motility and sperm-egg interaction. Moreover, the gene set analyses revealed many significant functional terms, including fertilization, sperm motility, calcium channel regulation, and SNARE proteins. Most of these terms are directly implicated in sperm physiology and male fertility.

CONCLUSIONS

This study contributes to the identification of genetic variants and biological processes underlying sire fertility. These findings can provide opportunities for improving bull fertility via marker-assisted selection.

MeSH-Informed Enrichment Analysis and MeSH-Guided Semantic Similarity Among Functional Terms and Gene Products in Chicken

Biomedical vocabularies and ontologies aid in recapitulating biological knowledge. The annotation of gene products is mainly accelerated by Gene Ontology (GO), and more recently by Medical Subject Headings (MeSH). Here, we report a suite of MeSH packages for chicken in Bioconductor, and illustrate some features of different MeSH-based analyses, including MeSH-informed enrichment analysis and MeSH-guided semantic similarity among terms and gene products, using two lists of chicken genes available in public repositories. The two published datasets that were employed represent (i) differentially expressed genes, and (ii) candidate genes under selective sweep or epistatic selection. The comparison of MeSH with GO overrepresentation analyses suggested not only that MeSH supports the findings obtained from GO analysis, but also that MeSH is able to further enrich the representation of biological knowledge and often provide more interpretable results. Based on the hierarchical structures of MeSH and GO, we computed semantic similarities among vocabularies, as well as semantic similarities among selected genes. These yielded the similarity levels between significant functional terms, and the annotation of each gene yielded the measures of gene similarity. Our findings show the benefits of using MeSH as an alternative choice of annotation in order to draw biological inferences from a list of genes of interest. We argue that the use of MeSH in conjunction with GO will be instrumental in facilitating the understanding of the genetic basis of complex traits.

Generalized enrichment analysis improves the detection of adverse drug events from the biomedical literature

Background

Identification of associations between marketed drugs and adverse events from the biomedical literature assists drug safety monitoring efforts. Assessing the significance of such literature-derived associations and determining the granularity at which they should be captured remains a challenge. Here, we assess how defining a selection of adverse event terms from MeSH, based on information content, can improve the detection of adverse events for drugs and drug classes.

Results

We analyze a set of 105,354 candidate drug adverse event pairs extracted from article indexes in MEDLINE. First, we harmonize extracted adverse event terms by aggregating them into higher-level MeSH terms based on the terms’ information content. Then, we determine statistical enrichment of adverse events associated with drug and drug classes using a conditional hypergeometric test that adjusts for dependencies among associated terms. We compare our results with methods based on disproportionality analysis (proportional reporting ratio, PRR) and quantify the improvement in signal detection with our generalized enrichment analysis (GEA) approach using a gold standard of drug-adverse event associations spanning 174 drugs and four events. For single drugs, the best GEA method (Precision: .92/Recall: .71/F1-measure: .80) outperforms the best PRR based method (.69/.69/.69) on all four adverse event outcomes in our gold standard. For drug classes, our GEA performs similarly (.85/.69/.74) when increasing the level of abstraction for adverse event terms. Finally, on examining the 1609 individual drugs in our MEDLINE set, which map to chemical substances in ATC, we find signals for 1379 drugs (10,122 unique adverse event associations) on applying GEA with p < 0.005.

Conclusions

We present an approach based on generalized enrichment analysis that can be used to detect associations between drugs, drug classes and adverse events at a given level of granularity, at the same time correcting for known dependencies among events. Our study demonstrates the use of GEA, and the importance of choosing appropriate abstraction levels to complement current drug safety methods. We provide an R package for exploration of alternative abstraction levels of adverse event terms based on information content.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1080-z) contains supplementary material, which is available to authorized users.

Genome-wide association mapping and pathway analysis of leukosis incidence in a US Holstein cattle population

Bovine leukosis virus is an oncogenic virus that infects B cells, causing bovine leukosis disease. This disease is known to have a negative impact on dairy cattle production and, because no treatment or vaccine is available, finding a possible genetic solution is important. Our objective was to perform a comprehensive genetic analysis of leukosis incidence in dairy cattle. Data on leukosis occurrence, pedigree and molecular information were combined into multitrait GBLUP models with milk yield (MY) and somatic cell score (SCS) to estimate genetic parameters and to perform whole-genome scans and pathway analysis. Leukosis data were available for 11 554 Holsteins daughters of 3002 sires from 112 herds in 16 US states. Genotypes from a 60K SNP panel were available for 961 of those bulls as well as for 2039 additional bulls. Heritability for leukosis incidence was estimated at about 8%, and the genetic correlations of leukosis disease incidence with MY and SCS were moderate at 0.18 and 0.20 respectively. The genome-wide scan indicated that leukosis is a complex trait, possibly modulated by many genes. The gene set analysis identified many functional terms that showed significant enrichment of genes associated with leukosis. Many of these terms, such as G-Protein Coupled Receptor Signaling Pathway, Regulation of Nucleotide Metabolic Process and different calcium-related processes, are known to be related to retrovirus infection. Overall, our findings contribute to a better understanding of the genetic architecture of this complex disease. The functional categories associated with leukosis may be useful in future studies on fine mapping of genes and development of dairy cattle breeding strategies.

Transcriptomic Analysis of the Mouse Mammary Gland Reveals New Insights for the Role of Serotonin in Lactation

Serotonin regulates numerous processes in the mammary gland. Our objective was to discover novel genes, pathways and functions which serotonin modulates during lactation. The rate limiting enzyme in the synthesis of non-neuronal serotonin is tryptophan-hydroxylase (TPH1). Therefore, we used TPH1 deficient dams (KO; serotonin deficient, n = 4) and compared them to wild-type (WT; n = 4) and rescue (RC; KO + 100 mg/kg 5-hydroxytryptophan injected daily, n = 4) dams. Mammary tissues were collected on day 10 of lactation. Total RNA extraction, amplification, library preparation and sequencing were performed following the Illumina mRNA-Seq. Overall, 97 and 204 genes (false discovery rate, FDR ≤ 0.01) exhibited a minimum of a 2-fold expression difference between WT vs. KO and WT vs. RC dams, respectively. Most differentially expressed genes were related to calcium homeostasis, apoptosis regulation, cell cycle, cell differentiation and proliferation, and the immune response. Additionally, gene set enrichment analysis using Gene Ontology and Medical Subject Headings databases revealed the alteration of several biological processes (FDR ≤ 0.01) including fat cell differentiation and lipid metabolism, regulation of extracellular signal-related kinase and mitogen-activated kinase cascades, insulin resistance, nuclear transport, membrane potential regulation, and calcium release from the endoplasmic reticulum into the cytosol. The majority of the biological processes and pathways altered in the KO dams are central for mammary gland homeostasis. Increasing peripheral serotonin in the RC dams affects specific pathways that favor lactation. Our data confirms the importance of serotonin during lactation in the mammary gland.

An application of MeSH enrichment analysis in livestock

An integral part of functional genomics studies is to assess the enrichment of specific biological terms in lists of genes found to be playing an important role in biological phenomena. Contrasting the observed frequency of annotated terms with those of the background is at the core of overrepresentation analysis (ORA). Gene Ontology (GO) is a means to consistently classify and annotate gene products and has become a mainstay in ORA. Alternatively, Medical Subject Headings (MeSH) offers a comprehensive life science vocabulary including additional categories that are not covered by GO. Although MeSH is applied predominantly in human and model organism research, its full potential in livestock genetics is yet to be explored. In this study, MeSH ORA was evaluated to discern biological properties of identified genes and contrast them with the results obtained from GO enrichment analysis. Three published datasets were employed for this purpose, representing a gene expression study in dairy cattle, the use of SNPs for genome‐wide prediction in swine and the identification of genomic regions targeted by selection in horses. We found that several overrepresented MeSH annotations linked to these gene sets share similar concepts with those of GO terms. Moreover, MeSH yielded unique annotations, which are not directly provided by GO terms, suggesting that MeSH has the potential to refine and enrich the representation of biological knowledge. We demonstrated that MeSH can be regarded as another choice of annotation to draw biological inferences from genes identified via experimental analyses. When used in combination with GO terms, our results indicate that MeSH can enhance our functional interpretations for specific biological conditions or the genetic basis of complex traits in livestock species.