Identification of novel splice variants in horn cancer by RNA-Seq analysis in Zebu cattle

Cross-population selection signatures in Canchim composite beef cattle

Analyses of livestock genomes have been used to detect selection signatures, which are genomic regions associated with traits under selection leading to a change in allele frequency. The objective of the present study was to characterize selection signatures in Canchim composite beef cattle using cross-population analyses with the founder Nelore and Charolais breeds. High-density single nucleotide polymorphism genotypes were available on 395 Canchim representing the target population, along with genotypes from 809 Nelore and 897 Charolais animals representing the reference populations. Most of the selection signatures were co-located with genes whose functions agree with the expectations of the breeding programs; these genes have previously been reported to associate with meat quality, as well as reproductive traits. Identified genes were related to immunity, adaptation, morphology, as well as behavior, could give new perspectives for understanding the genetic architecture of Canchim. Some selection signatures identified genes that were recently introduced in Canchim, such as the loci related to the polled trait.

Molecular portrait of squamous cell carcinoma of the bovine horn evaluated by high-throughput targeted exome sequencing: a preliminary report

Background

Squamous Cell Carcinoma of horn, also known as horn cancer, is a prevailing type of cancer in cattles especially Bos indicus. It is one of the most prevalent disease in Indian bullocks often resulting in death and huge economic losses to farmers. Here, we have reported the use of targeted exome sequencing to identify variants present in horn cancer affected horn mucosa tissue and blood of the same animal to identify some of the prevalent markers of horn cancer.

Results

We have observed higher number of variants present in tissue as compared to blood as well as among cancer samples compared to samples from normal animals. Eighty six and 1437 cancer-specific variants were identified among the predicted variants in blood and tissue samples, respectively. Total 25 missense variants were observed distributed over 18 genes. KRT8 gene coding for Keratin8, one of the key constituents of horn, displayed 5 missense variants. Additionally, three other genes involved in apoptosis pathway and two genes involved in antigen presentation and processing also contained variants.

Conclusions

Several genes involved in various apoptotic pathways were found to contain non-synonymous mutations. Keratin8 coding for Keratin, a chief constituent of horn was observed to have the highest number of mutations. In all, we present a preliminary report of mutations observed in horn cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-020-02683-y.

RNA-Seq Bayesian Network Exploration of Immune System in Bovine

Background:

The stress is one of main factors effects on production system. Several factors (both genetic and environmental elements) regulate immune response to stress.

Objectives:

In order to determine the major immune system regulatory genes underlying stress responses, a learning Bayesian network approach for those regulatory genes was applied to RNA-Seq data from a bovine leukocyte model system.

Material and Methods:

The transcriptome dataset GSE37447 was used from GEO and a Bayesian network on differentially expressed genes was learned to investigate the gene regulatory network.

Results:

Applying the method produced a strongly interconnected network with four genes (TERF2IP, PDCD10, DDX10 and CENPE) acting as nodes, suggesting these genes may be important in the transcriptome regulation program of stress response. Of these genes TERF2IP has been shown previously to regulate gene expression, act as a regulator of the nuclear factor-kappa B (NF-κB) signalling, and to activate expression of NF-κB target genes; PDCD10 encodes a conserved protein associated with cell apoptosis; DDX10 encodes a DEAD box protein and is believed to be associated with cellular growth and division; and CENPE involves unstable spindle microtubule capture at kinetochores. Together these genes are involved in DNA damage of apoptosis, RNA splicing, DNA repairing, and regulating cell division in the bovine genome. The topology of the learned Bayesian gene network indicated that the genes had a minimal interrelationship with each other. This type of structure, using the publically available computational tool, was also observed on human orthologous genes of the differentially expressed genes.

Conclusions:

Overall, the results might be used in transcriptomic-assisted selection and design of new drug targets to treat stress-related problems in bovines.

Deciphering GRINA/Lifeguard1: Nuclear Location, Ca2+ Homeostasis and Vesicle Transport

The Glutamate Receptor Ionotropic NMDA-Associated Protein 1 (GRINA) belongs to the Lifeguard family and is involved in calcium homeostasis, which governs key processes, such as cell survival or the release of neurotransmitters. GRINA is mainly associated with membranes of the endoplasmic reticulum, Golgi, endosome, and the cell surface, but its presence in the nucleus has not been explained yet. Here we dissect, with the help of different software tools, the potential roles of GRINA in the cell and how they may be altered in diseases, such as schizophrenia or celiac disease. We describe for the first time that the cytoplasmic N-terminal half of GRINA (which spans a Proline-rich domain) contains a potential DNA-binding sequence, in addition to cleavage target sites and probable PY-nuclear localization sequences, that may enable it to be released from the rest of the protein and enter the nucleus under suitable conditions, where it could participate in the transcription, alternative splicing, and mRNA export of a subset of genes likely involved in lipid and sterol synthesis, ribosome biogenesis, or cell cycle progression. To support these findings, we include additional evidence based on an exhaustive review of the literature and our preliminary data of the protein–protein interaction network of GRINA.

Identification and quantification of novel RNA isoforms in horn cancer of Bos indicus by comprehensive RNA-Seq

Horn cancer (HC) is a squamous cell carcinoma of horn, commonly observed in Bos indicus of the Asian countries. To elucidate the complexity of alternative splicing present in the HC, high-throughput sequencing and analysis of HC and matching horn normal (HN) tissue were carried out. A total of 535,067 and 849,077 reads were analysed after stringent quality filtering for HN and HC, respectively. Cufflinks pipeline for transcriptome analysis revealed 4786 novel splice isoforms comprising 2432 exclusively in HC, 2055 exclusively in HN and 298 in both the conditions. Based on pathway clustering and in silico verification, 102 novel splice isoforms were selected and further analysed with respect to change in protein sequence using Blastp. Finally, fourteen novel splicing events supported both by Cufflinks and UCSC genome browser were selected and confirmed expression by RT-qPCR. Future studies targeted at in-depth characterization of these potential candidate splice isoforms might be helpful in the development of relevant biomarkers for early diagnosis of HC. The results reported in this study refine the available information on transcriptome repertoire of bovine species and boost the research in the line of development of relevant biomarkers for early diagnosis of HC.

Assessment of myometrial transcriptome changes associated with spontaneous human labour by high-throughput RNA-seq

The transition of the human uterus from a quiescent to a contractile state takes place over a number of weeks. On such biological time scales, cellular phenotype is modified by changes in the transcriptome, which in turn is under the control of the underlying endocrine, paracrine, and biophysical processes resulting from the ongoing pregnancy. In this study, we characterize the transition of the human myometrial transcriptome at term from not in labour (NIL) to in labour (LAB) using high throughput RNA sequencing (RNA-seq). RNA was isolated from the myometrium of uterine biopsies from patients at term who were not in labour (n = 5) and at term in spontaneous labour (n = 5) without augmentation. A total of 143.6 million separate reads were sequenced, achieving, on average, ∼13 times coverage of the expressed human transcriptome per sample. Principal component analysis indicated that the NIL and LAB transcriptomes could be distinguished as two distinct clusters. A comparison of the NIL and LAB groups, using three different statistical approaches (baySeq, edgeR, and DESeq), demonstrated an overlap of 764 differentially expressed genes. A comparison with currently available microarray data revealed only a partial overlap in differentially expressed genes. We conclude that the described RNA-seq data sets represent the first fully annotated catalogue of expressed mRNAs in human myometrium. When considered together, the full expression repertoire and the differentially expressed gene sets should provide an excellent resource for formulating new hypotheses of physiological function, as well as the discovery of novel therapeutic targets.