Read-depth based approach on whole genome resequencing data reveals important insights into the copy number variation (CNV) map of major global buffalo breeds.

BACKGROUND: Elucidating genome-wide structural variants including copy number variations (CNVs) have gained increased significance in recent times owing to their contribution to genetic diversity and association with important pathophysiological states. The present study aimed to elucidate the high-resolution CNV map of six different global buffalo breeds using whole genome resequencing data at two coverages (10X and 30X). Post-quality control, the sequence reads were aligned to the latest draft release of the Bubaline genome. The genome-wide CNVs were elucidated using a read-depth approach in CNVnator with different bin sizes. Adjacent CNVs were concatenated into copy number variation regions (CNVRs) in different breeds and their genomic coverage was elucidated. RESULTS: Overall, the average size of CNVR was lower at 30X coverage, providing finer details. Most of the CNVRs were either deletion or duplication type while the occurrence of mixed events was lesser in number on a comparative basis in all breeds. The average CNVR size was lower at 30X coverage (0.201 Mb) as compared to 10X (0.013 Mb) with the finest variants in Banni buffaloes. The maximum number of CNVs was observed in Murrah (2627) and Pandharpuri (25,688) at 10X and 30X coverages, respectively. Whereas the minimum number of CNVs were scored in Surti at both coverages (2092 and 17,373). On the other hand, the highest and lowest number of CNVRs were scored in Jaffarabadi (833 and 10,179 events) and Surti (783 and 7553 events) at both coverages. Deletion events overnumbered duplications in all breeds at both coverages. Gene profiling of common overlapped genes and longest CNVRs provided important insights into the evolutionary history of these breeds and indicate the genomic regions under selection in respective breeds. CONCLUSION: The present study is the first of its kind to elucidate the high-resolution CNV map in major buffalo populations using a read-depth approach on whole genome resequencing data. The results revealed important insights into the divergence of major global buffalo breeds along the evolutionary timescale.

Changes in m6A RNA methylation of goat lung following PPRV infection.

Peste des petits ruminants (PPR) is an acute, highly contagious viral disease of goats and sheep, caused by the Peste des petits ruminants virus (PPRV). Earlier studies suggest the involvement of diverse regulatory mechanisms in PPRV infection. Methylation at N6 of Adenosine called m6A is a type RNA modification that influences various physiological and pathological phenomena. As the lung tissue represents the primary target organ of PPRV, the present study explored the m6A changes and their functional significance in PPRV disease pathogenesis. m6A-seq analysis revealed 1289 m6A peaks to be significantly altered in PPRV infected lung in comparison to normal lung, out of which 975 m6A peaks were hypomethylated and 314 peaks were hypermethylated. Importantly, hypomethylated genes were enriched in Interleukin-4 and Interleukin-13 signaling and various processes associated with extracellular matrix organization. Further, of the 843 differentially m6A-containing cellular transcripts, 282 transcripts were also found to be differentially expressed. Functional analysis revealed that these 282 transcripts are significantly enriched in signaling by Interleukins, extracellular matrix organization, cytokine signaling in the immune system, signaling by receptor tyrosine kinases, and Toll-like Receptor Cascades. We also found m6A reader HNRNPC and the core component of methyltransferase complex METTL14 to be highly upregulated than the m6A readers - HNRNPA2B1 and YTHDF1 at the transcriptome level. These findings suggest that alteration in the m6A landscape following PPRV is implicated in diverse processes including Interleukin signaling.

Review: genetic background of milk fatty acid synthesis in bovines.

Milk fat composition is an important trait for the dairy industry as it directly influences the nutritional and technological properties of milk and other dairy products. The synthesis of milk fat is a complex process regulated by a network of genes. Thus, understanding the genetic variation and molecular mechanisms regulating milk fat synthesis will help to improve the nutritional quality of dairy products. In this review, we provide an overview of milk fat synthesis in bovines along with the candidate genes involved in the pathway. We also discuss de novo synthesis of fatty acids (ACSS, ACACA, FASN), uptake of FAs (FATP, FAT, LPL), intracellular activation and channelling of FAs (ACSL, FABP), elongation (EVOLV6), desaturation (SCD, FADS), formation of triglycerides (GPAM, AGPAT, LIPIN, DGAT), and milk lipid secretion (BTN1A1, XDH, PLIN2). The genetic variability of individual fatty acids will help to develop selection strategies for obtaining a healthier milk fat profile in bovines. Thus, this review will offer a potential understanding of the molecular mechanisms that regulate milk fat synthesis in bovines.

An insight in proteome profiling of Tuta absoluta larvae after entomopathogenic fungal infection.

Tuta absoluta (L.) (Lepidoptera: Gelechiidae), a major pest of solanaceous plant species, causes serious losses in the agriculture sector around the globe. For better pest management, entomopathogenic fungi such as Beauveria bassiana and Purpureocillium lilacinum, play an efficient role in suppressing the pest population. The present study was carried out to analyse the effects post fungal infections through proteome profiling using an Orbitrap Fusion Tribrid mass spectrometer. A total of 2,201 proteins were identified from the fourth instar larvae of T. absoluta, of which 442 and 423 proteins were significantly dysregulated upon infection with P. lilacinum and B. bassiana respectively. The potential proteins related to immune systems as well as detoxification processes showed significant alterations after post fungal infection. Studies on T. absoluta proteomics and genomics as well as the consequences of entomopathogenic fungal infection on the immune response of this insect could provide an initial framework for exploring more fungus-host interactions for the development of better strategies for integrated pest management.

Integrated analysis of long-noncoding RNA and circular RNA expression in Peste-des-Petits-Ruminants Virus (PPRV) infected marmoset B lymphocyte (B95a) cells.

Peste-des-Petits-Ruminants (PPR) or goat plague is an important viral disease of sheep and goats caused by the small ruminant morbilli virus or PPR virus (PPRV). Long non coding RNAs (lncRNA) and circular RNAs (circRNA) play a pivotal role in several biological processes including regulation of virus-host interactions. The present study explored the expression of lncRNA, circRNA and their functions in PPRV infected B-lymphocyte (B95a) cells. The results revealed a total of 4531 lncRNA and 2348 circRNA expression in both mock and PPRV infected samples. Analysis of differentially expressed (DE) RNA identified 123 DE-lncRNA and 39 DE-circRNA as significantly dysregulated. Functional analysis of cis-target genes of DE-lncRNA indicated activation of TCF dependent WNT signaling and PKN1 stimulated transcription process. Interactions (sponging) of microRNA (miRNA) revealed 344 DE-lncRNA-miRNA and 93 DE-circRNA-miRNA pairs. The competing endogenous RNA (ceRNA) network of lncRNA/circRNA-miRNA-mRNA in PPRV infected B95a cells was represented by 69 ceRNA pairs. We validated the DE-circRNA by targeted amplification and sequencing of back spliced junctions (BSJs). The present study revealed a profile of lncRNA, circRNA and their potential ceRNA network in PPRV infection. The results provide insight for better understanding of PPRV-host interactions.

Genome-wide transcriptome profiling of CSF virus challenged monocyte-derived macrophages provides distinct insights into immune response of Landrace and indigenous Ghurrah pigs.

The present study was undertaken to characterize the distinct immune response in indigenous Ghurrah and exotic Landrace pigs by challenging monocyte-derived macrophages (MDMs) with CSF virus under in-vitro conditions and assessing the variations in the transcriptome profile at 48 h post-infection (hpi). RNA-sequencing was carried out in infected and non-infected MDMs of Ghurrah (n = 3) and Landrace (n = 3) piglets prior- as well as post-stimulation. MDMs of Ghurrah showed greater immune regulation in response to CSF infection with 518 significantly differentially expressed genes (DEG) in infected versus non-infected MDMs, as compared to only 31 DEGs in Landrace MDMs. In Landrace, the principal regulators of inflammation (IL1α, IL1ß and TNF) were upregulated in infected cells while in Ghurrah, these were downregulated. Overall, macrophages from indigenous Ghurrah showed more immunological dysregulation in response to virulent CSF virus infection as compared to the exotic Landrace pigs.

Porcine Circovirus type 2 infected myocardial tissue transcriptome signature.

The goal of this study was to compare the global gene expression profile in cardiac tissues of pig infected with porcine circovirus 2 (PCV2) to that of healthy cells. Since PCV2 infection causes severe cardiovascular lesions, the myocardial tissue model was chosen for this study. In High-throughput transcriptome analysis, DESeq2 and CLC genomics workbench analyses revealed a total of 196 significantly differentially expressed genes (DEGs) (p-value < 0.05). Furthermore, 194 transcripts were upregulated, while only two were downregulated (HSPA6 and DNAJA1), with fold changes ranging from 16.293 to -10.002. Among the KEGG canonical pathways targeted by the DEGs in the functional analysis, adrenergic signalling in cardiomyocytes, Cardiac Muscle Contraction, Hypertrophic Cardiomyopathy (HCM), and Dilated Cardiomyopathy (DCM) tends to be enriched. The differentially expressed highly connected (DEHC) biomarker genes in pathogenicity of PCV2 infection, such as LDB3, MYOZ2, CASQ2, TNNT2, MLC2V, MYBPC3, ACTC1, TCAP, TNNI3, TRDN, CSRP3, MYL3, RYR2, LMOD2, MYH7, etc., were identified using protein-protein interaction (PPI) network analysis. The study might provide detailed information on the dysregulated genes and biological pathways in infected myocardial tissues that may be essential for PCV2-related heart pathology.

Genome-wide elucidation of CNV regions and their association with production and reproduction traits in composite Vrindavani cattle.

The present study was aimed to analyze the genome-wide copy number variations (CNVs) in Vrindavani composite cattle and concatenate them into CNV regions (CNVRs), and finally test the association of CNVRs with different production and reproduction traits. Genotypic data, generated on BovineSNP50 Beadchip (v3) array for 96 Vrindavani animals, was used to elucidate the CNVs at the genome level. Intensity data covering over 53,218 SNP genotypes on bovine genome was used. Algorithm based on Hidden Markov Model was employed in PennCNV program to detect, normalize and filter CNVs across the genome. 252 putative CNVs, detected via PennCNV program, in different individuals were concatenated into 71 CNV regions (CNVRs) using CNVRuler program. Association of CNVRs with important (re)production traits in Vrindavani animals was assessed using linear regression. Five CNVRs were found to be significantly associated with ten important (re)production traits. The genes harbored in these regions provided useful insights into the association of CNVRs with genes and ultimately the variation at phenotype level. Important genes that overlapped with CNVRs included WASHC4, HS6ST3, MBNL2, TOLLIP, PIDD1 and TSPAN4. Furthermore, the CNVRs were found to overlap with important QTLs available in AnimalQTL database which affect milk yield and composition along with reproduction and immune function traits. The copy number states of three enes were validated using digital droplet PCR technique. The results from the present study significantly enhance the understanding about CNVs in Vrindavani cattle and should help establish its CNV map. The study will also enable further investigation on association of these variants with important traits of economic interest including disease incidence.

Differential expression of long non-coding RNAs under Peste des petits ruminants virus (PPRV) infection in goats.

Peste des petits ruminants (PPR) characterized by fever, sore mouth, conjunctivitis, gastroenteritis, and pneumonia, is an acute, highly contagious viral disease of sheep and goats. The role of long non-coding RNAs (lncRNAs) in PPRV infection has not been explored to date. In this study, the transcriptome profiles of virulent Peste des petits ruminants virus (PPRV) infected goat tissues - lung and spleen were analyzed to identify the role of lncRNAs in PPRV infection. A total of 13,928 lncRNA transcripts were identified, out of which 170 were known lncRNAs. Intergenic lncRNAs (7625) formed the major chunk of the novel lncRNA transcripts. Differential expression analysis revealed that 15 lncRNAs (11 downregulated and 4 upregulated) in the PPRV infected spleen samples and 16 lncRNAs (13 downregulated and 3 upregulated) in PPRV infected lung samples were differentially expressed as compared to control. The differentially expressed lncRNAs (DElncRNAs) possibly regulate various immunological processes related to natural killer cell activation, antigen processing and presentation, and B cell activity, by regulating the expression of mRNAs through the cis- or trans-regulatory mechanism. Functional enrichment analysis of differentially expressed mRNAs (DEmRNAs) revealed enrichment of immune pathways and biological processes in concordance with the pathways in which correlated lncRNA-neighboring genes were enriched. The results suggest that a coordinated immune response is raised in both lung and spleen tissues of the goat through mRNA-lncRNA crosstalk.

A Multi-copy Nucleic Acid-Based Diagnostic Test for Bovine Tropical Theileriosis.

BACKGROUND: Bovine tropical theileriosis (BTT) is a haemoprotozoan tick-borne disease that implicates huge losses to livestock in terms of considerable mortality and morbidity in tropical and subtropical regions of the globe. Currently available diagnostic methods have less specificity and sensitivity towards the detection of Theileria species. Therefore, an attempt was made to diagnose Theileria annulata by targeting a multi-copy gene, viz. mitochondrially encoded cytochrome b (MT-CYB) gene via polymerase chain reaction (PCR) in different agro-zones of India. METHODS AND RESULTS: 129 cattle blood samples were collected from major livestock rearing regions of India and processed for both molecular and microscopic techniques. Screening of Giemsa-stained thin blood smears was able to detect 14 samples (10.85%) as positive for T. annulata. However, the MT-CYB gene-based PCR assay detected 107 samples (82.94%) positive for T. annulata out of 129 samples. Furthermore, the MT-CYB gene-based PCR assay was standardized in terms of its sensitivity and specificity. Specificity of PCR assay was evaluated against other common haemoprotozoan parasites of tropical countries viz. Babesia bigemina, Anaplasma marginale and Trypanosoma evansi. The multi-copy MT-CYB gene-based PCR assay provided an optimum level of sensitivity (up to the level of 10 femtogram) and high specificity. Haematological examination (Hb, PCV and TLC) of 113 samples revealed significantly (p < 0.05) decreased Hb and PCV levels in positive animals in comparison with the control group of healthy animals. However, the control group had significantly higher (p < 0.001) TLC levels than the positive group. CONCLUSION: The MT-CYB gene-based PCR assay was found to be highly sensitive that can accurately detect the occurrence of T. annulata infection in carrier animals which are potential infection sources to healthier populations in naive demographic locations through infected ticks.

Systems biology under heat stress in Indian cattle.

Transcriptome profiling of Vrindavani and Tharparkar cattle (n = 5 each) revealed that more numbers of genes were dysregulated in Vrindavani than in Tharparkar. A contrast in gene expression was observed with 18.9 % of upregulated genes in Vrindavani downregulated in Tharparkar and 17.8% upregulated genes in Tharparkar downregulated in Vrindavani. Functional annotation of genes differentially expressed in Tharparkar and Vrindavani revealed that the systems biology in Tharparkar is moving towards counteracting the effects due to heat stress. Unlike Vrindavani, Tharparkar is not only endowed with higher expression of the scavengers (UBE2G1, UBE2S, and UBE2H) of misfolded proteins but also with protectors (VCP, Serp1, and CALR) of naïve unfolded proteins. Further, higher expression of the antioxidants in Tharparkar enables it to cope up with higher levels of free radicals generated as a result of heat stress. In this study, we found relevant genes dysregulated in Tharparkar in the direction that can counter heat stress.

m6A modification of RNA and its role in cancer, with a special focus on lung cancer.

Epitranscriptomics involves functionally relevant biochemical modifications of RNA taking place at the transcriptome level without a change in the sequence of ribonucleotides. Several types of modifications that affect the processing and function of differentRNA types have been reported. Methylation at N6 of Adenosine called m6A is one such modification, quite widespread in occurrence and reported in snRNAs, lncRNAs, circRNAs, rRNAs, miRNAs, and most abundantly, in mRNAs. The significant implications of m6A in various types of cancers are being widely recognized. Here, we give a brief about the enzymes that install the m6A modification (= m6A writers), that remove it (= m6A erasers) and certain RNA binding proteins (= m6A readers) which affect the fate of the m6A-containing RNA by recruiting various proteins. We also discuss the relevance of m6A in ncRNAs in various cancer types, followed by a discussion on the role of m6A of mRNA and ncRNA in lung cancer.

Canine Parvovirus and Its Non-Structural Gene 1 as Oncolytic Agents: Mechanism of Action and Induction of Anti-Tumor Immune Response.

The exploration into the strategies for the prevention and treatment of cancer is far from complete. Apart from humans, cancer has gained considerable importance in animals because of increased awareness towards animal health and welfare. Current cancer treatment regimens are less specific towards tumor cells and end up harming normal healthy cells. Thus, a highly specific therapeutic strategy with minimal side effects is the need of the hour. Oncolytic viral gene therapy is one such specific approach to target cancer cells without affecting the normal cells of the body. Canine parvovirus (CPV) is an oncolytic virus that specifically targets and kills cancer cells by causing DNA damage, caspase activation, and mitochondrial damage. Non-structural gene 1 (NS1) of CPV, involved in viral DNA replication is a key mediator of cytotoxicity of CPV and can selectively cause tumor cell lysis. In this review, we discuss the oncolytic properties of Canine Parvovirus (CPV or CPV2), the structure of the NS1 protein, the mechanism of oncolytic action as well as role in inducing an antitumor immune response in different tumor models.

Apoptin as a Tumor-Specific Therapeutic Agent: Current Perspective on Mechanism of Action and Delivery Systems.

Cancer remains one of the leading causes of death worldwide in humans and animals. Conventional treatment regimens often fail to produce the desired outcome due to disturbances in cell physiology that arise during the process of transformation. Additionally, development of treatment regimens with no or minimum side-effects is one of the thrust areas of modern cancer research. Oncolytic viral gene therapy employs certain viral genes which on ectopic expression find and selectively destroy malignant cells, thereby achieving tumor cell death without harming the normal cells in the neighborhood. Apoptin, encoded by Chicken Infectious Anemia Virus' VP3 gene, is a proline-rich protein capable of inducing apoptosis in cancer cells in a selective manner. In normal cells, the filamentous Apoptin becomes aggregated toward the cell margins, but is eventually degraded by proteasomes without harming the cells. In malignant cells, after activation by phosphorylation by a cancer cell-specific kinase whose identity is disputed, Apoptin accumulates in the nucleus, undergoes aggregation to form multimers, and prevents the dividing cancer cells from repairing their DNA lesions, thereby forcing them to undergo apoptosis. In this review, we discuss the present knowledge about the structure of Apoptin protein, elaborate on its mechanism of action, and summarize various strategies that have been used to deliver it as an anticancer drug in various cancer models.

Contrasting Gene Expression Profiles of Monocytes and Lymphocytes From Peste-Des-Petits-Ruminants Virus Infected Goats.

In this study, transcriptome analysis of PPRV infected PBMC subsets-T helper cells, T cytotoxic cells, monocytes, and B lymphocytes was done to delineate their role in host response. PPRV was found to infect lymphocytes and not monocytes. The established receptor for PPRV-SLAM was found downregulated in lymphocytes and non-differentially expressed in monocytes. A profound deviation in the global gene expression profile with a large number of unique upregulated genes (851) and downregulated genes (605) was observed in monocytes in comparison to lymphocytes. ISGs-ISG15, Mx1, Mx2, RSAD2, IFIT3, and IFIT5 that play a role in antiviral response and the genes for viral sensors-MDA5, LGP2, and RIG1, were found to be upregulated in lymphocytes and downregulated in monocytes. The transcription factors-IRF-7 and STAT-1 that regulate expression of most of the ISGs were found activated in lymphocytes and not in monocytes. Interferon signaling pathway and RIG1 like receptor signaling pathway were found activated in lymphocytes and not in monocytes. This contrast in gene expression profiles and signaling pathways indicated the predominant role of lymphocytes in generating the antiviral response against PPRV in goats, thus, giving us new insights into host response to PPRV.

Dysregulated miRNAome and Proteome of PPRV Infected Goat PBMCs Reveal a Coordinated Immune Response.

In this study, the miRNAome and proteome of virulent Peste des petits ruminants virus (PPRV) infected goat peripheral blood mononuclear cells (PBMCs) were analyzed. The identified differentially expressed miRNAs (DEmiRNAs) were found to govern genes that modulate immune response based on the proteome data. The top 10 significantly enriched immune response processes were found to be governed by 98 genes. The top 10 DEmiRNAs governing these 98 genes were identified based on the number of genes governed by them. Out of these 10 DEmiRNAs, 7 were upregulated, and 3 were downregulated. These include miR-664, miR-2311, miR-2897, miR-484, miR-2440, miR-3533, miR-574, miR-210, miR-21-5p, and miR-30. miR-664 and miR-484 with proviral and antiviral activities, respectively, were upregulated in PPRV infected PBMCs. miR-210 that inhibits apoptosis was downregulated. miR-21-5p that decreases the sensitivity of cells to the antiviral activity of IFNs and miR-30b that inhibits antigen processing and presentation by primary macrophages were downregulated, indicative of a strong host response to PPRV infection. miR-21-5p was found to be inhibited on IPA upstream regulatory analysis of RNA-sequencing data. This miRNA that was also highly downregulated and was found to govern 16 immune response genes in the proteome data was selected for functional validation vis-a-vis TGFBR2 (TGF-beta receptor type-2). TGFBR2 that regulates cell differentiation and is involved in several immune response pathways was found to be governed by most of the identified immune modulating DEmiRNAs. The decreased luciferase activity in Dual Luciferase Reporter Assay indicated specific binding of miR-21-5p and miR-484 to their target thus establishing specific binding of the miRNAs to their targets.This is the first report on the miRNAome and proteome of virulent PPRV infected goat PBMCs.

Selection and validation of suitable reference genes for qPCR gene expression analysis in goats and sheep under Peste des petits ruminants virus (PPRV), lineage IV infection.

Identification of suitable candidate reference genes is an important prerequisite for validating the gene expression data obtained from downstream analysis of RNA sequencing using quantitative real time PCR (qRT-PCR). Though existence of a universal reference gene is myth, commonly used reference genes can be assessed for expression stability to confer their suitability to be used as candidate reference genes in gene expression studies. In this study, we evaluated the expression stability of ten most commonly used reference genes (GAPDH, ACTB, HSP90, HMBS, 18S rRNA, B2M, POLR2A, HPRT1, ACAC, YWHAZ) in fourteen different Peste des petits ruminants virus (PPRV) infected tissues of goats and sheep. RefFinder and RankAggreg software were used to deduce comprehensive ranking of reference genes. Our results suggested HMBS and B2M in goats and HMBS and HPRT1 in sheep can be used as suitable endogenous controls in gene expression studies of PPRV infection irrespective of tissues and condition as a whole, thus eliminating the use of tissue specific/ condition specific endogenous controls. We report for the first time suitable reference genes for gene expression studies in PPRV infected tissues. The reference genes determined here can be useful for future studies on gene expression in sheep and goat infected with PPRV, thus saving extra efforts and time of repeating the reference gene determination and validation.