Sudeep Sahadevan

The Expression Pattern of microRNAs in Granulosa Cells of Subordinate and Dominant Follicles during the Early Luteal Phase of the Bovine Estrous Cycle

This study aimed to investigate the miRNA expression patterns in granulosa cells of subordinate (SF) and dominant follicle (DF) during the early luteal phase of the bovine estrous cycle. For this, miRNA enriched total RNA isolated from granulosa cells of SF and DF obtained from heifers slaughtered at day 3 and day 7 of the estrous cycle was used for miRNAs deep sequencing. The results revealed that including 17 candidate novel miRNAs, several known miRNAs (n = 291–318) were detected in SF and DF at days 3 and 7 of the estrous cycle of which 244 miRNAs were common to all follicle groups. The let-7 families, bta-miR-10b, bta-miR-26a, bta-miR-99b and bta-miR-27b were among abundantly expressed miRNAs in both SF and DF at both days of the estrous cycle. Further analysis revealed that the expression patterns of 16 miRNAs including bta-miR-449a, bta-miR-449c and bta-miR-222 were differentially expressed between the granulosa cells of SF and DF at day 3 of the estrous cycle. However, at day 7 of the estrous cycle, 108 miRNAs including bta-miR-409a, bta-miR-383 and bta-miR-184 were differentially expressed between the two groups of granulosa cell revealing the presence of distinct miRNA expression profile changes between the two follicular stages at day 7 than day 3 of the estrous cycle. In addition, unlike the SF, marked temporal miRNA expression dynamics was observed in DF groups between day 3 and 7 of the estrous cycle. Target gene prediction and pathway analysis revealed that major signaling associated with follicular development including Wnt signaling, TGF-beta signaling, oocyte meiosis and GnRH signaling were affected by differentially expressed miRNAs. Thus, this study highlights the miRNA expression patterns of granulosa cells in subordinate and dominant follicles that could be associated with follicular recruitment, selection and dominance during the early luteal phase of the bovine estrous cycle.

MicroRNA Expression Profile in Bovine Granulosa Cells of Preovulatory Dominant and Subordinate Follicles during the Late Follicular Phase of the Estrous Cycle

In bovine, ovarian follicles grow in a wave-like fashion with commonly 2 or 3 follicular waves emerging per estrous cycle. The dominant follicle of the follicular wave which coincides with the LH-surge becomes ovulatory, leaving the subordinate follicles to undergo atresia. These physiological processes are controlled by timely and spatially expressed genes and gene products, which in turn are regulated by post-transcriptional regulators. MicroRNAs, a class of short non-coding RNA molecules, are one of the important posttranscriptional regulators of genes associated with various cellular processes. Here we investigated the expression pattern of miRNAs in granulosa cells of bovine preovulatory dominant and subordinate follicles during the late follicular phase of bovine estrous cycle using Illumina miRNA deep sequencing. In addition to 11 putative novel miRNAs, a total of 315 and 323 known miRNAs were detected in preovulatory dominant and subordinate follicles, respectively. Moreover, in comparison with the subordinate follicles, a total of 64 miRNAs were found to be differentially expressed in preovulatory dominant follicles, of which 34 miRNAs including the miR-132 and miR-183 clusters were significantly enriched, and 30 miRNAs including the miR-17-92 cluster, bta-miR-409a and bta-miR-378 were significantly down regulated in preovulatory dominant follicles. In-silico pathway analysis revealed that canonical pathways related to oncogenesis, cell adhesion, cell proliferation, apoptosis and metabolism were significantly enriched by the predicted target genes of differentially expressed miRNAs. Furthermore, Luciferase reporter assay analysis showed that one of the differentially regulated miRNAs, the miR-183 cluster miRNAs, were validated to target the 3´-UTR of FOXO1 gene. Moreover FOXO1 was highly enriched in granulosa cells of subordinate follicles in comparison with the preovulatory dominant follicles demonstrating reciprocal expression pattern with miR-183 cluster miRNAs. In conclusion, the presence of distinct sets of miRNAs in granulosa cells of preovulatory dominant and subordinate follicles supports the potential role of miRNAs in post-transcriptional regulation of genes involved in bovine follicular development during the late follicular phase of the estrous cycle.

RNA Deep Sequencing Reveals Novel Candidate Genes and Polymorphisms in Boar Testis and Liver Tissues with Divergent Androstenone Levels

Boar taint is an unpleasant smell and taste of pork meat derived from some entire male pigs. The main causes of boar taint are the two compounds androstenone (5α-androst-16-en-3-one) and skatole (3-methylindole). It is crucial to understand the genetic mechanism of boar taint to select pigs for lower androstenone levels and thus reduce boar taint. The aim of the present study was to investigate transcriptome differences in boar testis and liver tissues with divergent androstenone levels using RNA deep sequencing (RNA-Seq). The total number of reads produced for each testis and liver sample ranged from 13,221,550 to 33,206,723 and 12,755,487 to 46,050,468, respectively. In testis samples 46 genes were differentially regulated whereas 25 genes showed differential expression in the liver. The fold change values ranged from −4.68 to 2.90 in testis samples and −2.86 to 3.89 in liver samples. Differentially regulated genes in high androstenone testis and liver samples were enriched in metabolic processes such as lipid metabolism, small molecule biochemistry and molecular transport. This study provides evidence for transcriptome profile and gene polymorphisms of boars with divergent androstenone level using RNA-Seq technology. Digital gene expression analysis identified candidate genes in flavin monooxygenease family, cytochrome P450 family and hydroxysteroid dehydrogenase family. Moreover, polymorphism and association analysis revealed mutation in IRG6, MX1, IFIT2, CYP7A1, FMO5 and KRT18 genes could be potential candidate markers for androstenone levels in boars. Further studies are required for proving the role of candidate genes to be used in genomic selection against boar taint in pig breeding programs.

Identification of the Novel Candidate Genes and Variants in Boar Liver Tissues with Divergent Skatole Levels Using RNA Deep Sequencing

Boar taint is the unpleasant odour of meat derived from non-castrated male pigs, caused by the accumulation of androstenone and skatole in fat. Skatole is a tryptophan metabolite produced by intestinal bacteria in gut and catabolised in liver. Since boar taint affects consumer’s preference, the aim of this study was to perform transcriptome profiling in liver of boars with divergent skatole levels in backfat by using RNA-Seq. The total number of reads produced for each liver sample ranged from 11.8 to 39.0 million. Approximately 448 genes were differentially regulated (p-adjusted <0.05). Among them, 383 genes were up-regulated in higher skatole group and 65 were down-regulated (p<0.01, FC>1.5). Differentially regulated genes in the high skatole liver samples were enriched in metabolic processes such as small molecule biochemistry, protein synthesis, lipid and amino acid metabolism. Pathway analysis identified the remodeling of epithelial adherens junction and TCA cycle as the most dominant pathways which may play important roles in skatole metabolism. Differential gene expression analysis identified candidate genes in ATP synthesis, cytochrome P450, keratin, phosphoglucomutase, isocitrate dehydrogenase and solute carrier family. Additionally, polymorphism and association analysis revealed that mutations in ATP5B, KRT8, PGM1, SLC22A7 and IDH1 genes could be potential markers for skatole levels in boars. Furthermore, expression analysis of exon usage of three genes (ATP5B, KRT8 and PGM1) revealed significant differential expression of exons of these genes in different skatole levels. These polymorphisms and exon expression differences may have impacts on the gene activity ultimately leading to skatole variation and could be used as genetic marker for boar taint related traits. However, further validation is required to confirm the effect of these genetic markers in other pig populations in order to be used in genomic selection against boar taint in pig breeding programs.

Pathway based analysis of genes and interactions influencing porcine testis samples from boars with divergent androstenone content in back fat

One of the primary factors contributing to boar taint is the level of androstenone in porcine adipose tissues. A majority of the studies performed to identify candidate biomarkers for the synthesis of androstenone in testis tissues follow a reductionist approach, identifying and studying the effect of biomarkers individually. Although these studies provide detailed information on individual biomarkers, a global picture of changes in metabolic pathways that lead to the difference in androstenone synthesis is still missing. The aim of this work was to identify major pathways and interactions influencing steroid hormone synthesis and androstenone biosynthesis using an integrative approach to provide a bird’s eye view of the factors causing difference in steroidogenesis and androstenone biosynthesis. For this purpose, we followed an analysis procedure merging together gene expression data from boars with divergent levels of androstenone and pathway mapping and interaction network retrieved from KEGG database. The interaction networks were weighted with Pearson correlation coefficients calculated from gene expression data and significant interactions and enriched pathways were identified based on these networks. The results show that 1,023 interactions were significant for high and low androstenone animals and that a total of 92 pathways were enriched for significant interactions. Although published articles show that a number of these enriched pathways were activated as a result of downstream signaling of steroid hormones, we speculate that the significant interactions in pathways such as glutathione metabolism, sphingolipid metabolism, fatty acid metabolism and significant interactions in cAMP-PKA/PKC signaling might be the key factors determining the difference in steroidogenesis and androstenone biosynthesis between boars with divergent androstenone levels in our study. The results and assumptions presented in this study are from an in-silico analysis done at the gene expression level and further laboratory experiments at genomic, proteomic or metabolomic level are necessary to validate these findings.

Transcriptome profile of lung dendritic cells after in vitro porcine reproductive and respiratory syndrome virus (PRRSV) infection

The porcine reproductive and respiratory syndrome (PRRS) is an infectious disease that leads to high financial and production losses in the global swine industry. The pathogenesis of this disease is dependent on a multitude of factors, and its control remains problematic. The immune system generally defends against infectious diseases, especially dendritic cells (DCs), which play a crucial role in the activation of the immune response after viral infections. However, the understanding of the immune response and the genetic impact on the immune response to PRRS virus (PRRSV) remains incomplete. In light of this, we investigated the regulation of the host immune response to PRRSV in porcine lung DCs using RNA-sequencing (RNA-Seq). Lung DCs from two different pig breeds (Pietrain and Duroc) were collected before (0 hours) and during various periods of infection (3, 6, 9, 12, and 24 hours post infection (hpi)). RNA-Seq analysis revealed a total of 20,396 predicted porcine genes, which included breed-specific differentially expressed immune genes. Pietrain and Duroc infected lung DCs showed opposite gene expression courses during the first time points post infection. Duroc lung DCs reacted more strongly and distinctly than Pietrain lung DCs during these periods (3, 6, 9, 12 hpi). Additionally, cluster analysis revealed time-dependent co-expressed groups of genes that were involved in immune-relevant pathways. Key clusters and pathways were identified, which help to explain the biological and functional background of lung DCs post PRRSV infection and suggest IL-1β1 as an important candidate gene. RNA-Seq was also used to characterize the viral replication of PRRSV for each breed. PRRSV was able to infect and to replicate differently in lung DCs between the two mentioned breeds. These results could be useful in investigations on immunity traits in pig breeding and enhancing the health of pigs.

Exploring novel mechanistic insights in Alzheimer’s disease by assessing reliability of protein interactions

Protein interaction networks are widely used in computational biology as a graphical means of representing higher-level systemic functions in a computable form. Although, many algorithms exist that seamlessly collect and measure protein interaction information in network models, they often do not provide novel mechanistic insights using quantitative criteria. Measuring information content and knowledge representation in network models about disease mechanisms becomes crucial particularly when exploring new target candidates in a well-defined functional context of a potential disease mechanism. To this end, we have developed a knowledge-based scoring approach that uses literature-derived protein interaction features to quantify protein interaction confidence. Thereby, we introduce the novel concept of knowledge cliffs, regions of the interaction network where a significant gap between high scoring and low scoring interactions is observed, representing a divide between established and emerging knowledge on disease mechanism. To show the application of this approach, we constructed and assessed reliability of a protein-protein interaction model specific to Alzheimer’s disease, which led to screening, and prioritization of four novel protein candidates. Evaluation of the identified candidates showed that two of them are already followed in clinical trials for testing potential AD drugs.

Preliminary study of FMO1, FMO5, CYP21, ESR1, PLIN2 and SULT2A1 as candidate gene for compounds related to boar taint

An association study between polymorphisms of six genes and boar taint related compounds androstenone, skatole and indole was performed in a boar population (n=370). Significant association (P<0.05) was detected for SNP of FMO5 (g.494A>G) with all boar taint compounds, SNP of CYP21 (g.3911T>C) with skatole and indole, and SNP of ESR1 (g.672C>T) with androstenone and indole. mRNA expression of CYP21 and ESR1 was higher in CAB (castrated boar) compared to non-castrated boars; whereas, the expression of FMO5 and ESR1 was higher in LBT (low boar taint) compared to HBT (high boar taint) in liver tissue. FMO5, CYP21 and ESR1 proteins were less detectable in HBT compared with LBT and CAB in liver tissues. These findings suggest that FMO5, CYP21 and ESR1 gene variants might have effects on the boar taint compounds.

Identification of gene co-expression clusters in liver tissues from multiple porcine populations with high and low backfat androstenone phenotype

BackgroundBoar taint is principally caused by accumulation of androstenone and skatole in adipose tissues. Studies have shown high heritability estimates for androstenone whereas skatole production is mainly dependent on nutritional factors. Androstenone is a lipophilic steroid mainly metabolized in liver. Majority of the studies on hepatic androstenone metabolism focus only on a single breed and very few studies account for population similarities/differences in gene expression patterns. In this work, we concentrated on population similarities in gene expression to identify the common genes involved in hepatic androstenone metabolism of multiple pig populations. Based on androstenone measurements, publicly available gene expression datasets from three porcine populations were compiled into either low or high androstenone dataset. Gene expression correlation coefficients from these datasets were converted to rank ratios and joint probabilities of these rank ratios were used to generate dataset specific co-expression clusters. Finally, these networks were clustered using a graph clustering technique.ResultsCluster analysis identified a number of statistically significant co-expression clusters in the dataset. Further enrichment analysis of these clusters showed that one of the clusters from low androstenone dataset was highly enriched for xenobiotic, drug, cholesterol and lipid metabolism and cytochrome P450 associated metabolism of drugs and xenobiotics. Literature references revealed that a number of genes in this cluster were involved in phase I and phase II metabolism. Physical and functional similarity assessment showed that the members of this cluster were dispersed across multiple clusters in high androstenone dataset, possibly indicating a weak co-expression of these genes in high androstenone dataset.ConclusionsBased on these results we hypothesize that majority of the genes in this cluster forms a signature co-expression cluster in low androstenone dataset in our experiment and that majority of the members of this cluster might be responsible for hepatic androstenone metabolism across all the three populations used in our study. We propose these results as a background work towards understanding breed similarities in hepatic androstenone metabolism. Additional large scale experiments using data from multiple porcine breeds are necessary to validate these findings.