MicroRNA-guided prioritization of genome-wide association signals reveals the importance of microRNA-target gene networks for complex traits in cattle

Signatures of positive selection after the introduction of genomic selection in the Finnish Ayrshire population

The Finnish Ayrshire (FAY) belongs to the Nordic Red breeds and is characterized by high milk yield, high milk components, good fertility, and functional conformation. The FAY breeding program is based on genomic selection. Despite the benefits of selection on breeding values, autozygosity in the genome may increase due to selection, and increased autozygosity may cause inbreeding depression in selected traits. However, there is lack of studies concerning selection signatures in the FAY after genomic selection introduction. The aim of this study was to identify signatures of selection in FAY after the introduction of genomic selection. Genomic data included 45,834 SNPs. The genotyped animals were divided into 2 groups: animals born before genomic selection introduction (6,108 cows) and animals born after genomic selection introduction (47,361 cows). We identified the selection signatures using 3 complementary methods: 2 based on identification of selection signatures from runs of homozygosity (ROH) islands and one based on the decay of site-specific extended haplotype between populations at SNP sites (Rsb). In total, we identified 34 ROH islands on chromosomes 1, 3, 6, 8, 12-15, 17, 19, 22, and 26 in FAY animals born before genomic selection (between 1980 and 2011) and 30 ROH islands on chromosomes 1-3, 13-17, 22, and 25-26 in FAY animals born after genomic selection introduction (between 2015 and 2020). We additionally detected 22 ΔROH islands on chromosomes 2-3, 11, 13, 14, 16, 18, 20, and 25-26. Finally, a total of 31 Rsb regions on chromosomes 2, 3, 14, 18, 20, and 25 were identified. Based on the results, genomic selection has favored certain alleles and haplotypes on genomic regions related to traits relevant in the FAY breeding program: milk production, fertility, growth, beef production traits, and feed efficiency. Several genes related to these traits (e.g., PLA2G4A, MECR, CHUK, COX15, RICTOR, SHISA9, and SEMA4G) overlapped or partially overlapped the observed selection signature regions. The association of genotypes within these regions and their effects on traits relevant in the FAY breeding program should be studied and genetic regions undergoing selection monitored in the FAY population.

Active Fraction of Polyrhachis Vicina Roger (AFPR) Ameliorate Depression Induced Inflammation Response by FTO/miR-221-3p/SOCS1 Axis

Purpose

Neuroinflammation is a significant etiological factor in the development of depression. Traditional Chinese medicine (TCM) has demonstrated notable efficacy in the treatment of inflammation. Our previous study surfaces that the active fraction of Polyrhachis vicina Roger (AFPR) has antidepressant and anti-neuroinflammatory effects, but the specific mechanisms remain to be elucidated. The objective of this study was to examine the impact of AFPR on inflammation in depression via the FTO/miR-221-3p/SOCS1 axis.

Methods

Chronic unpredictable stress (CUMS)-induced rats and LPS-induced BV2 cells were employed to simulate depression models in vivo and in vitro. The levels of inflammatory factors were detected using the ELISA assay. The expression of genes and proteins was detected using qRT-PCR and Western blot. Gene interactions were detected using the dual luciferase reporter gene. Protein-RNA interactions were investigated using RNA methylation immunoprecipitation (MeRIP) and RNA immunoprecipitation (RIP). Neuroinflammation in the brain was examined through H&E staining, while neuronal apoptosis was assessed using TUNEL staining.

Results

The results showed that AFPR ameliorated depression induced inflammation by increasing SOCS1 expression. However, SOCS1 was identified as a target of miR-221-3p. Overexpression of miR-221-3p decreased the expression of SOCS1 and increased the levels of NF-κB, IL-7, and IL-6. In addition, we found that miR-221-3p was regulated by FTO-mediated m6A modification through MeRIP and RIP experiments. Interference with miR-221-3p and overexpression of FTO resulted in increased SOCS1 gene expression and decreased levels of NF-κB, IL-7, and IL-6, which were reversed by AFPR.

Conclusion

AFPR inhibits the maturation of pri-miR-221-3p through FTO-mediated m6A modification, reduces the production of miR-221-3p, increases the expression of SOCS1, and reduces the level of inflammation, thereby improving depressive symptoms.

Evaluation of circulating microRNA profiles in blood as potential candidate biomarkers in a subacute ruminal acidosis cow model - a pilot study

BACKGROUND

Subacute ruminal acidosis (SARA) is a metabolic disorder often observed in high-yielding dairy cows, that are fed diets high in concentrates. We hypothesized that circulating miRNAs in blood of cows could serve as potential candidate biomarkers to detect animals with metabolic dysbalances such as SARA. MicroRNAs (miRNAs) are a class of small non-coding RNAs, serving as regulators of a plethora of molecular processes. To test our hypothesis, we performed a pilot study with non-lactating Holstein-Friesian cows fed a forage diet (FD; 0% concentrate, n = 4) or a high-grain diet (HG; 65% concentrate, n = 4) to induce SARA. Comprehensive profiling of miRNA expression in plasma and leucocytes were performed by next generation sequencing (NGS). The success of our model to induce SARA was evaluated based on ruminal pH and was evidenced by increased time spent with a pH threshold of 5.8 for an average period of 320 min/d.

RESULTS

A total of 520 and 730 miRNAs were found in plasma and leucocytes, respectively. From these, 498 miRNAs were shared by both plasma and leucocytes, with 22 miRNAs expressed exclusively in plasma and 232 miRNAs expressed exclusively in leucocytes. Differential expression analysis revealed 10 miRNAs that were up-regulated and 2 that were down-regulated in plasma of cows when fed the HG diet. A total of 63 circulating miRNAs were detected exclusively in the plasma of cows with SARA, indicating that these animals exhibited a higher number and diversity of circulating miRNAs. Considering the total read counts of miRNAs expressed when fed the HG diet, differentially expressed miRNAs ( log2 fold change) and known function, we have identified bta-miR-11982, bta-miR-1388-5p, bta-miR-12034, bta-miR-2285u, and bta-miR-30b-3p as potential candidates for SARA-biomarker in cows by NGS. These were further subjected to validation using small RNA RT-qPCR, confirming the promising role of bta-miR-30b-3p and bta-miR-2285.

CONCLUSION

Our data demonstrate that dietary change impacts the release and expression of miRNAs in systemic circulation, which may modulate post-transcriptional gene expression in cows undergoing SARA. Particularly, bta-miR-30b-3p and bta-miR-2285 might serve as promising candidate biomarker predictive for SARA and should be further validated in larger cohorts.

In silico prioritisation of microRNA-associated common variants in multiple sclerosis

BACKGROUND

Genome-wide association studies (GWAS) have highlighted over 200 autosomal variants associated with multiple sclerosis (MS). However, variants in non-coding regions such as those encoding microRNAs have not been explored thoroughly, despite strong evidence of microRNA dysregulation in MS patients and model organisms. This study explores the effect of microRNA-associated variants in MS, through the largest publicly available GWAS, which involved 47,429 MS cases and 68,374 controls.

METHODS

We identified SNPs within the coordinates of microRNAs, ± 5-kb microRNA flanking regions and predicted 3'UTR target-binding sites using miRBase v22, TargetScan 7.0 RNA22 v2.0 and dbSNP v151. We established the subset of microRNA-associated SNPs which were tested in the summary statistics of the largest MS GWAS by intersecting these datasets. Next, we prioritised those microRNA-associated SNPs which are among known MS susceptibility SNPs, are in strong linkage disequilibrium with the former or meet a microRNA-specific Bonferroni-corrected threshold. Finally, we predicted the effects of those prioritised SNPs on their microRNAs and 3'UTR target-binding sites using TargetScan v7.0, miRVaS and ADmiRE.

RESULTS

We have identified 30 candidate microRNA-associated variants which meet at least one of our prioritisation criteria. Among these, we highlighted one microRNA variant rs1414273 (MIR548AC) and four 3'UTR microRNA-binding site variants within SLC2A4RG (rs6742), CD27 (rs1059501), MMEL1 (rs881640) and BCL2L13 (rs2587100). We determined changes to the predicted microRNA stability and binding site recognition of these microRNA and target sites.

CONCLUSIONS

We have systematically examined the functional, structural and regulatory effects of candidate MS variants among microRNAs and 3'UTR targets. This analysis allowed us to identify candidate microRNA-associated MS SNPs and highlights the value of prioritising non-coding RNA variation in GWAS. These candidate SNPs could influence microRNA regulation in MS patients. Our study is the first thorough investigation of both microRNA and 3'UTR target-binding site variation in multiple sclerosis using GWAS summary statistics.

MicroRNA sequence variation can impact interactions with target mRNA in cattle

BACKGROUND

Small non-coding microRNAs (miRNAs) are important modulators at post-transcriptional levels. Single-nucleotide polymorphisms (SNPs) located in miRNA genes can alter the secondary structure of pre-miRNA to either impair or promote the miRNA maturation processes. Furthermore, SNPs located in the miRNA seed regions can stabilize or disturb miRNA-target interactions, thereby, quantitatively influence the expression of target genes. Therefore, the main objective of this study was to detect SNPs in bovine miRNAs using the whole-genome re-sequencing datasets of 1632 cattle of five breeds from the 1000 bull genomes project.

RESULTS

In total, our study identified 1109, 334, and 130 SNPs in the miRNA precursor, mature, and seed regions, respectively. The heterozygosity values were generally less than 0.3, and the minor allele frequencies (MAFs) were mainly less than 0.1. Most SNPs were in Hardy-Weinberg equilibrium (HWE) (HWE-P > 0.05). Furthermore, we found that the majority of SNPs (MAF > 0.1 and HWE-P > 0.05) in the miRNA seed regions altered the repertoire of target genes, which in turn were enriched in different KEGG pathways or GO terms. Thus target prediction for bta-miR-2888 revealed loss of 309 target genes and gain of 691 target genes. The 691 gained target genes were significantly enriched in 60 KEGG pathways and 21 GO terms.

CONCLUSION

In summary, our study identified candidate SNPs in miRNA precursor, mature, and seed regions that are likely to affect target RNA interactions, thereby potentially influencing cattle phenotypic traits.

Emerging roles of noncoding micro RNAs and circular RNAs in bovine mastitis: Regulation, breeding, diagnosis, and therapy

Bovine mastitis is one of the most troublesome and costly problems in the modern dairy industry, which is not only difficult to monitor, but can also cause economic losses while having significant implications on public health. However, efficacious preventative methods and therapy are still lacking. Moreover, new drugs and therapeutic targets are in increasing demand due to antibiotic restrictions. In recent years, noncoding RNAs have gained popularity as a topic in pathological and genetic studies. Meanwhile, there is growing evidence that they play a role in regulating various biological processes and developing novel treatment platforms. In light of this, this review focuses on two types of noncoding RNAs, micro RNAs and circular RNAs, and summarizes their characterizations, relationships, potential applications as selection markers, diagnostic or treatment targets and potential applications in RNA-based therapy, in order to shed new light on further research.

Functional annotation of regulatory elements in cattle genome reveals the roles of extracellular interaction and dynamic change of chromatin states in rumen development during weaning

We profiled landscapes of bovine regulatory elements and explored dynamic changes of chromatin states in rumen development during weaning. The regulatory elements (15 chromatin states) and their coordinated activities in cattle were defined through genome-wide profiling of four histone modifications, CTCF-binding, DNA accessibility, DNA methylation, and transcriptome in rumen epithelial tissues. Each chromatin state presented specific enrichment for sequence ontology, methylation, trait-associated variants, transcription, gene expression-associated variants, selection signatures, and evolutionarily conserved elements. During weaning, weak enhancers and flanking active transcriptional start sites (TSS) were the most dynamic chromatin states and occurred in tandem with significant variations in gene expression and DNA methylation, significantly associated with stature, production, and reproduction economic traits. By comparing with in vitro cultured epithelial cells and in vivo rumen tissues, we showed the commonness and uniqueness of these results, especially the roles of cell interactions and mitochondrial activities in tissue development.

The nature of coxofemoral joint pathology across family Canidae

We evaluated coxofemoral joints from museum specimens of: Vulpes lagopus; Vulpes vulpes; Vulpes velox; Nyctereutes procyonoides; Urocyon cinereoargenteus; Aenocyon [Canis] dirus; Canis latrans; Canis lupus lupus; Canis lupus familiaris; C. l. familiaris × latrans; and Canis dingo. Acetabular components included: fossa; articular surface; medial and lateral articular margins; and periarticular surfaces. Acetabular components variably revealed: osteophyte-like features; varying appearance of articular margin rims (especially contour changes); rough bone surfaces (especially fossa and articular surface); and surface wear. Proximal femoral components included: articular surface; articular margin; periarticular surfaces; and joint capsule attachment. Femoral components variably revealed: rough bone surface; bone loss; articular margin osteophyte-like features; caudal post-developmental mineralized prominence; and enthesophytes along the joint capsule attachment. Non-metric multidimensional scaling was used to analyze right-left asymmetric relationships between observed traits, across taxa. Significantly different acetabular trait asymmetry involved only C. latrans-C. l. familiaris; V. vulpes-N. procyonoides, and U. cinereoargenteus-N. procyonoides. There were no significant lateralized differences in proximal femoral traits involving modern canids, ancient and modern C. l. familiaris, or modern vulpines. Thus, the observations were strongly bilateral. We hypothesized high similarity of traits across taxa. The data confirm the hypothesis and strongly suggest broad and deep morphological and mechanistic conservation that almost certainly pre-existed (at least) all modern canids. Further zoological studies are needed to evaluate phylogenic implications in greater detail.

Integrated Small RNA Sequencing, Transcriptome and GWAS Data Reveal microRNA Regulation in Response to Milk Protein Traits in Chinese Holstein Cattle

Milk protein is one of the most important economic traits in the dairy industry. Yet, the regulatory network of miRNAs for the synthesis of milk protein in mammary is poorly understood. Samples from 12 Chinese Holstein cows with three high ( ≥ 3.5%) and three low ( ≤ 3.0%) phenotypic values for milk protein percentage in lactation and non-lactation were examined through deep small RNA sequencing. We characterized 388 known and 212 novel miRNAs in the mammary gland. Differentially expressed analysis detected 28 miRNAs in lactation and 52 miRNAs in the non-lactating period with a highly significant correlation with milk protein concentration. Target prediction and correlation analysis identified some key miRNAs and their targets potentially involved in the synthesis of milk protein. We analyzed for enrichments of GWAS signals in miRNAs and their correlated targets. Our results demonstrated that genomic regions harboring DE miRNA genes in lactation were significantly enriched with GWAS signals for milk protein percentage traits and that enrichments within DE miRNA targets were significantly higher than in random gene sets for the majority of milk production traits. This integrated study on the transcriptome and posttranscriptional regulatory profiles between significantly differential phenotypes of milk protein concentration provides new insights into the mechanism of milk protein synthesis, which should reveal the regulatory mechanisms of milk secretion.

Identification of loci associated with pathological outcomes in Holstein cattle infected with Mycobacterium avium subsp. paratuberculosis using whole-genome sequence data

Bovine paratuberculosis (PTB), caused by Mycobacterium avium subsp. paratuberculosis (MAP), is a chronic granulomatous enteritis that affects cattle worldwide. According to their severity and extension, PTB-associated histological lesions have been classified into the following groups; focal, multifocal, and diffuse. It is unknown whether these lesions represent sequential stages or divergent outcomes. In the current study, the associations between host genetic and pathology were explored by genotyping 813 Spanish Holstein cows with no visible lesions (N = 373) and with focal (N = 371), multifocal (N = 33), and diffuse (N = 33) lesions in gut tissues and regional lymph nodes. DNA from peripheral blood samples of these animals was genotyped with the bovine EuroG MD Bead Chip, and the corresponding genotypes were imputed to whole-genome sequencing (WGS) data using the 1000 Bull genomes reference population. A genome-wide association study (GWAS) was performed using the WGS data and the presence or absence of each type of histological lesion in a case–control approach. A total of 192 and 92 single nucleotide polymorphisms (SNPs) defining 13 and 9 distinct quantitative trait loci (QTLs) were highly-associated (P ≤ 5 × 10⁻⁷) with the multifocal (heritability = 0.075) and the diffuse (heritability = 0.189) lesions, respectively. No overlap was seen in the SNPs controlling these distinct pathological outcomes. The identified QTLs overlapped with some QTLs previously associated with PTB susceptibility, bovine tuberculosis susceptibility, clinical mastitis, somatic cell score, bovine respiratory disease susceptibility, tick resistance, IgG level, and length of productive life. Pathway analysis with candidate genes overlapping the identified QTLs revealed a significant enrichment of the keratinization pathway and cholesterol metabolism in the animals with multifocal and diffuse lesions, respectively. To test whether the enrichment of SNP variants in candidate genes involved in the cholesterol metabolism was associated with the diffuse lesions; the levels of total cholesterol were measured in plasma samples of cattle with focal, multifocal, or diffuse lesions or with no visible lesions. Our results showed reduced levels of plasma cholesterol in cattle with diffuse lesions. Taken together, our findings suggested that the variation in MAP-associated pathological outcomes might be, in part, genetically determined and indicative of distinct host responses.

Characterization of microRNA expression in B cells derived from Japanese black cattle naturally infected with bovine leukemia virus by deep sequencing

Bovine leukemia virus (BLV) is the causative agent of enzootic bovine leukosis (EBL), a malignant B cell lymphoma. However, the mechanisms of BLV-associated lymphomagenesis remain poorly understood. Here, after deep sequencing, we performed comparative analyses of B cell microRNAs (miRNAs) in cattle infected with BLV and those without BLV. In BLV-infected cattle, BLV-derived miRNAs (blv-miRNAs) accounted for 38% of all miRNAs in B cells. Four of these blv-miRNAs (blv-miR-B1-5p, blv-miR-B2-5p, blv-miR-B4-3p, and blv-miR-B5-5p) had highly significant positive correlations with BLV proviral load (PVL). The read counts of 90 host-derived miRNAs (bta-miRNAs) were significantly down-regulated in BLV-infected cattle compared to those in uninfected cattle. Only bta-miR-375 had a positive correlation with PVL in BLV-infected cattle and was highly expressed in the B cell lymphoma tissue of EBL cattle. There were a few bta-miRNAs that correlated with BLV tax/rex gene expression; however, BLV AS1 expression had a significant negative correlation with many of the down-regulated bta-miRNAs that are important for tumor development and/or tumor suppression. These results suggest that BLV promotes lymphomagenesis via AS1 and blv-miRNAs, rather than tax/rex, by down-regulating the expression of bta-miRNAs that have a tumor-suppressing function, and this downregulation is linked to increased PVL.

Integrative genomics of the mammalian alveolar macrophage response to intracellular mycobacteria

Background

Bovine TB (bTB), caused by infection with Mycobacterium bovis, is a major endemic disease affecting global cattle production. The key innate immune cell that first encounters the pathogen is the alveolar macrophage, previously shown to be substantially reprogrammed during intracellular infection by the pathogen. Here we use differential expression, and correlation- and interaction-based network approaches to analyse the host response to infection with M. bovis at the transcriptome level to identify core infection response pathways and gene modules. These outputs were then integrated with genome-wide association study (GWAS) data sets to enhance detection of genomic variants for susceptibility/resistance to M. bovis infection.

Results

The host gene expression data consisted of RNA-seq data from bovine alveolar macrophages (bAM) infected with M. bovis at 24 and 48 h post-infection (hpi) compared to non-infected control bAM. These RNA-seq data were analysed using three distinct computational pipelines to produce six separate gene sets: 1) DE genes filtered using stringent fold-change and P-value thresholds (DEG-24: 378 genes, DEG-48: 390 genes); 2) genes obtained from expression correlation networks (CON-24: 460 genes, CON-48: 416 genes); and 3) genes obtained from differential expression networks (DEN-24: 339 genes, DEN-48: 495 genes). These six gene sets were integrated with three bTB breed GWAS data sets by employing a new genomics data integration tool—gwinteR. Using GWAS summary statistics, this methodology enabled detection of 36, 102 and 921 prioritised SNPs for Charolais, Limousin and Holstein-Friesian, respectively.

Conclusions

The results from the three parallel analyses showed that the three computational approaches could identify genes significantly enriched for SNPs associated with susceptibility/resistance to M. bovis infection. Results indicate distinct and significant overlap in SNP discovery, demonstrating that network-based integration of biologically relevant transcriptomics data can leverage substantial additional information from GWAS data sets. These analyses also demonstrated significant differences among breeds, with the Holstein-Friesian breed GWAS proving most useful for prioritising SNPS through data integration. Because the functional genomics data were generated using bAM from this population, this suggests that the genomic architecture of bTB resilience traits may be more breed-specific than previously assumed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07643-w.

The micro-RNA content of unsorted cryopreserved bovine sperm and its relation to the fertility of sperm after sex-sorting

BACKGROUND

The use of sex-sorted sperm in cattle assisted reproduction is constantly increasing. However, sperm fertility can substantially differ between unsorted (conventional) and sex-sorted semen batches of the same sire. Sperm microRNAs (miRNA) have been suggested as promising biomarkers of bull fertility the last years. In this study, we hypothesized that the miRNA profile of cryopreserved conventional sperm is related to bull fertility after artificial insemination with X-bearing sperm. For this purpose, we analyzed the miRNA profile of 18 conventional sperm samples obtained from nine high- (HF) and nine low-fertility (LF) bulls that were contemporaneously used to produce conventional and sex-sorted semen batches. The annual 56-day non-return rate for each semen type (NRRconv and NRRss, respectively) was recorded for each bull.

RESULTS

In total, 85 miRNAs were detected. MiR-34b-3p and miR-100-5p were the two most highly expressed miRNAs with their relative abundance reaching 30% in total. MiR-10a-5p and miR-9-5p were differentially expressed in LF and HF samples (false discovery rate < 10%). The expression levels of miR-9-5p, miR-34c, miR-423-5p, miR-449a, miR-5193-5p, miR-1246, miR-2483-5p, miR-92a, miR-21-5p were significantly correlated to NRRss but not to NRRconv. Based on robust regression analysis, miR-34c, miR-7859 and miR-342 showed the highest contribution to the prediction of NRRss.

CONCLUSIONS

A set of miRNAs detected in conventionally produced semen batches were linked to the fertilizing potential of bovine sperm after sex-sorting. These miRNAs should be further evaluated as potential biomarkers of a sire's suitability for the production of sex-sorted sperm.

Multi-Omics Approach Reveals miR-SNPs Affecting Muscle Fatty Acids Profile in Nelore Cattle

MicroRNAs (miRNAs) are key regulators of gene expression, potentially affecting several biological processes, whose function can be altered by sequence variation. Hence, the integration of single nucleotide polymorphisms (SNP) and miRNAs can explain individual differences in economic traits. To provide new insights into the effects of SNPs on miRNAs and their related target genes, we carried out a multi-omic analysis to identify SNPs in miRNA mature sequences (miR-SNPs) associated with fatty acid (FA) composition in the Nelore cattle. As a result, we identified 3 miR-SNPs in different miRNAs (bta-miR-2419-3p, bta-miR-193a-2, and bta-miR-1291) significantly associated with FA traits (p-value < 0.02, Bonferroni corrected). Among these, the rs110817643C>T, located in the seed sequence of the bta-miR-1291, was associated with different ω6 FAs, polyunsaturated FA, and polyunsaturated:saturated FA ratios. Concerning the other two miR-SNPs, the rs43400521T>C (located in the bta-miR-2419-3p) was associated with C12:0 and C18:1 cis-11 FA, whereas the rs516857374A>G (located in the bta-miR-193a-2) was associated with C18:3 ω6 and ratio of ω6/ω3 traits. Additionally, to identify potential biomarkers for FA composition, we described target genes affected by these miR-SNPs at the mRNA or protein level. Our multi-omics analysis outlines the effects of genetic polymorphism on miRNA, and it highlights miR-SNPs and target candidate genes that control beef fatty acid composition.

MicroRNA-345-5p regulates depression by targeting suppressor of cytokine signaling 1

BACKGROUND/AIMS

MicroRNA(miR)-345-5p plays a key role in various cellular functions. However, the function of miR-345-5p in resistant depression (TRD) is unclear. The aim of this study was to evaluate the role and mechanism of miR-345-5p in the treatment of resistance depression (TRD).

METHODS

RT-qPCR was used to detect the expression of miR-345-5p in BV-2 microglia. CCK-8 method and flow cytometry were used for cell viability and apoptosis of microglia. Target gene prediction and screening, and luciferase reporter assays were used to verify the downstream target gene of miR-345-5p. Western blot was used to analyze the protein expression of related proteins.

RESULTS

miR-345-5p increased the cell viability of BV-2 microglia and the expression level of pro-inflammatory cytokines. In addition, the conditioned medium of microglia treated with miR-345-5p reduced the cell viability of HT22 hippocampal cells and caused S-phase arrest. The miR-345-5p-treated microglia induced apoptosis by regulating the expression levels of Bax, Bcl-2, pro-caspase-3, and cleaved caspase-3. Furthermore, SOCS1 was a direct target of miR-345-5p, and overexpression of SOCS1 was able to reverse the proapoptotic effect of miR-345-5p on activation of microglia on hippocampal neurons.

CONCLUSION

miR-345-5p induced inflammatory damage in hippocampal neurons by activating microglia. MiR-345-5p may be an effective target for TRD therapy.

Integrating RNA-Seq with GWAS reveals novel insights into the molecular mechanism underpinning ketosis in cattle

Background

Ketosis is a common metabolic disease during the transition period in dairy cattle, resulting in long-term economic loss to the dairy industry worldwide. While genetic selection of resistance to ketosis has been adopted by many countries, the genetic and biological basis underlying ketosis is poorly understood.

Results

We collected a total of 24 blood samples from 12 Holstein cows, including 4 healthy and 8 ketosis-diagnosed ones, before (2 weeks) and after (5 days) calving, respectively. We then generated RNA-Sequencing (RNA-Seq) data and seven blood biochemical indicators (bio-indicators) from leukocytes and plasma in each of these samples, respectively. By employing a weighted gene co-expression network analysis (WGCNA), we detected that 4 out of 16 gene-modules, which were significantly engaged in lipid metabolism and immune responses, were transcriptionally (FDR < 0.05) correlated with postpartum ketosis and several bio-indicators (e.g., high-density lipoprotein and low-density lipoprotein). By conducting genome-wide association signal (GWAS) enrichment analysis among six common health traits (ketosis, mastitis, displaced abomasum, metritis, hypocalcemia and livability), we found that 4 out of 16 modules were genetically (FDR < 0.05) associated with ketosis, among which three were correlated with postpartum ketosis based on WGCNA. We further identified five candidate genes for ketosis, including GRINA, MAF1, MAFA, C14H8orf82 and RECQL4. Our phenome-wide association analysis (Phe-WAS) demonstrated that human orthologues of these candidate genes were also significantly associated with many metabolic, endocrine, and immune traits in humans. For instance, MAFA, which is involved in insulin secretion, glucose response, and transcriptional regulation, showed a significantly higher association with metabolic and endocrine traits compared to other types of traits in humans.

Conclusions

In summary, our study provides novel insights into the molecular mechanism underlying ketosis in cattle, and highlights that an integrative analysis of omics data and cross-species mapping are promising for illustrating the genetic architecture underpinning complex traits.

Epigenomics and genotype-phenotype association analyses reveal conserved genetic architecture of complex traits in cattle and human

Background

Lack of comprehensive functional annotations across a wide range of tissues and cell types severely hinders the biological interpretations of phenotypic variation, adaptive evolution, and domestication in livestock. Here we used a combination of comparative epigenomics, genome-wide association study (GWAS), and selection signature analysis, to shed light on potential adaptive evolution in cattle.

Results

We cross-mapped 8 histone marks of 1300 samples from human to cattle, covering 178 unique tissues/cell types. By uniformly analyzing 723 RNA-seq and 40 whole genome bisulfite sequencing (WGBS) datasets in cattle, we validated that cross-mapped histone marks captured tissue-specific expression and methylation, reflecting tissue-relevant biology. Through integrating cross-mapped tissue-specific histone marks with large-scale GWAS and selection signature results, we for the first time detected relevant tissues and cell types for 45 economically important traits and artificial selection in cattle. For instance, immune tissues are significantly associated with health and reproduction traits, multiple tissues for milk production and body conformation traits (reflecting their highly polygenic architecture), and thyroid for the different selection between beef and dairy cattle. Similarly, we detected relevant tissues for 58 complex traits and diseases in humans and observed that immune and fertility traits in humans significantly correlated with those in cattle in terms of relevant tissues, which facilitated the identification of causal genes for such traits. For instance, PIK3CG, a gene highly specifically expressed in mononuclear cells, was significantly associated with both age-at-menopause in human and daughter-still-birth in cattle. ICAM, a T cell-specific gene, was significantly associated with both allergic diseases in human and metritis in cattle.

Conclusion

Collectively, our results highlighted that comparative epigenomics in conjunction with GWAS and selection signature analyses could provide biological insights into the phenotypic variation and adaptive evolution. Cattle may serve as a model for human complex traits, by providing additional information beyond laboratory model organisms, particularly when more novel phenotypes become available in the near future.

qgg: an R package for large-scale quantitative genetic analyses

Summary

Here, we present the R package qgg, which provides an environment for large-scale genetic analyses of quantitative traits and diseases. The qgg package provides an infrastructure for efficient processing of large-scale genetic data and functions for estimating genetic parameters, and performing single and multiple marker association analyses and genomic-based predictions of phenotypes.

Availability and implementation

The qgg package is freely available. For the latest updates, user guides and example scripts, consult the main page http://psoerensen.github.io/qgg. The current release is available from CRAN (https://CRAN.R-project.org/package=qgg) for all major operating systems.

Supplementary information

Supplementary data are available at Bioinformatics online.

Isolation of microglia-derived extracellular vesicles: towards miRNA signatures and neuroprotection

The functional preservation of the central nervous system (CNS) is based on the neuronal plasticity and survival. In this context, the neuroinflammatory state plays a key role and involves the microglial cells, the CNS-resident macrophages. In order to better understand the microglial contribution to the neuroprotection, microglia-derived extracellular vesicles (EVs) were isolated and molecularly characterized to be then studied in neurite outgrowth assays. The EVs, mainly composed of exosomes and microparticles, are an important cell-to-cell communication process as they exhibit different types of mediators (proteins, lipids, nucleic acids) to recipient cells. The medicinal leech CNS was initially used as an interesting model of microglia/neuron crosstalk due to their easy collection for primary cultures. After the microglia-derived EV isolation following successive methods, we developed their large-scale and non-targeted proteomic analysis to (i) detect as many EV protein markers as possible, (ii) better understand the biologically active proteins in EVs and (iii) evaluate the resulting protein signatures in EV-activated neurons. The EV functional properties were also evaluated in neurite outgrowth assays on rat primary neurons and the RNAseq analysis of the microglia-derived EVs was performed to propose the most representative miRNAs in microglia-derived EVs. This strategy allowed validating the EV isolation, identify major biological pathways in EVs and corroborate the regenerative process in EV-activated neurons. In parallel, six different miRNAs were originally identified in microglia-derived EVs including 3 which were only known in plants until now. The analysis of the neuronal proteins under the microglial EV activation suggested possible miRNA-dependent regulation mechanisms. Taken together, this combination of methodologies showed the leech microglial EVs as neuroprotective cargos across species and contributed to propose original EV-associated miRNAs whose functions will have to be evaluated in the EV-dependent dialog between microglia and neurons.

Functional annotation of the cattle genome through systematic discovery and characterization of chromatin states and butyrate-induced variations

BACKGROUND

The functional annotation of genomes, including chromatin accessibility and modifications, is important for understanding and effectively utilizing the increased amount of genome sequences reported. However, while such annotation has been well explored in a diverse set of tissues and cell types in human and model organisms, relatively little data are available for livestock genomes, hindering our understanding of complex trait variation, domestication, and adaptive evolution. Here, we present the first complete global landscape of regulatory elements in cattle and explore the dynamics of chromatin states in rumen epithelial cells induced by the rumen developmental regulator-butyrate.

RESULTS

We established the first global map of regulatory elements (15 chromatin states) and defined their coordinated activities in cattle, through genome-wide profiling for six histone modifications, RNA polymerase II, CTCF-binding sites, DNA accessibility, DNA methylation, and transcriptome in rumen epithelial primary cells (REPC), rumen tissues, and Madin-Darby bovine kidney epithelial cells (MDBK). We demonstrated that each chromatin state exhibited specific enrichment for sequence ontology, transcription, methylation, trait-associated variants, gene expression-associated variants, selection signatures, and evolutionarily conserved elements, implying distinct biological functions. After butyrate treatments, we observed that the weak enhancers and flanking active transcriptional start sites (TSS) were the most dynamic chromatin states, occurred concomitantly with significant alterations in gene expression and DNA methylation, which was significantly associated with heifer conception rate and stature economic traits.

CONCLUSION

Our results demonstrate the crucial role of functional genome annotation for understanding genome regulation, complex trait variation, and adaptive evolution in livestock. Using butyrate to induce the dynamics of the epigenomic landscape, we were able to establish the correlation among nutritional elements, chromatin states, gene activities, and phenotypic outcomes.

Integrating Signals from Sperm Methylome Analysis and Genome-Wide Association Study for a Better Understanding of Male Fertility in Cattle

Decreased male fertility is a big concern in both human society and the livestock industry. Sperm DNA methylation is commonly believed to be associated with male fertility. However, due to the lack of accurate male fertility records (i.e., limited mating times), few studies have investigated the comprehensive impacts of sperm DNA methylation on male fertility in mammals. In this study, we generated 10 sperm DNA methylomes and performed a preliminary correlation analysis between signals from sperm DNA methylation and signals from large-scale (n = 27,214) genome-wide association studies (GWAS) of 35 complex traits (including 12 male fertility-related traits). We detected genomic regions, which experienced DNA methylation alterations in sperm and were associated with aging and extreme fertility phenotypes (e.g., sire-conception rate or SCR). In dynamic hypomethylated regions (HMRs) and partially methylated domains (PMDs), we found genes (e.g., HOX gene clusters and microRNAs) that were involved in the embryonic development. We demonstrated that genomic regions, which gained rather than lost methylations during aging, and in animals with low SCR were significantly and selectively enriched for GWAS signals of male fertility traits. Our study discovered 16 genes as the potential candidate markers for male fertility, including SAMD5 and PDE5A. Collectively, this initial effort supported a hypothesis that sperm DNA methylation may contribute to male fertility in cattle and revealed the usefulness of functional annotations in enhancing biological interpretation and genomic prediction for complex traits and diseases.

Comparative analyses of sperm DNA methylomes among human, mouse and cattle provide insights into epigenomic evolution and complex traits

Sperm DNA methylation is crucial for fertility and viability of offspring but epigenome evolution in mammals is largely understudied. By comparing sperm DNA methylomes and large-scale genome-wide association study (GWAS) signals between human and cattle, we aimed to examine the DNA methylome evolution and its associations with complex phenotypes in mammals. Our analysis revealed that genes with conserved non-methylated promoters (e.g., ANKS1A and WNT7A) among human and cattle were involved in common system and embryo development, and enriched for GWAS signals of body conformation traits in both species, while genes with conserved hypermethylated promoters (e.g., TCAP and CD80) were engaged in immune responses and highlighted by immune-related traits. On the other hand, genes with human-specific hypomethylated promoters (e.g., FOXP2 and HYDIN) were engaged in neuron system development and enriched for GWAS signals of brain-related traits, while genes with cattle-specific hypomethylated promoters (e.g., LDHB and DGAT2) mainly participated in lipid storage and metabolism. We validated our findings using sperm-retained nucleosome, preimplantation transcriptome, and adult tissue transcriptome data, as well as sequence evolutionary features, including motif binding sites, mutation rates, recombination rates and evolution signatures. In conclusion, our results demonstrate important roles of epigenome evolution in shaping the genetic architecture underlying complex phenotypes, hence enhance signal prioritization in GWAS and provide valuable information for human neurological disorders and livestock genetic improvement.

Genetic and epigenetic architecture of paternal origin contribute to gestation length in cattle

The length of gestation can affect offspring health and performance. Both maternal and fetal effects contribute to gestation length; however, paternal contributions to gestation length remain elusive. Using genome-wide association study (GWAS) in 27,214 Holstein bulls with millions of gestation records, here we identify nine paternal genomic loci associated with cattle gestation length. We demonstrate that these GWAS signals are enriched in pathways relevant to embryonic development, and in differentially methylated regions between sperm samples with long and short gestation length. We reveal that gestation length shares genetic and epigenetic architecture in sperm with calving ability, body depth, and conception rate. While several candidate genes are detected in our fine-mapping analysis, we provide evidence indicating ZNF613 as a promising candidate for cattle gestation length. Collectively, our findings support that the paternal genome and epigenome can impact gestation length potentially through regulation of the embryonic development.

Lingzhao Fang et al. studied the paternal genetic variants that affect gestational length in cattle. They found that paternal genes from pathways involved in embryonic development were associated with gestation length, and that these were often found in differentially methylated regions of the genome.