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.

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.

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.

Temporal dysregulation of genes in Lamb testis cell during sheeppox virus infection

The present study was aimed to elucidate the host-virus interactions using RNA-Seq analysis at 1h and 8h of post-infection of SPPV in LTC. The differentially expressed genes (DEGs) and the underlying mechanisms linked to the host immune responses were obtained. The protein-protein interaction (PPI) network analysis and Ingenuity pathway analysis (IPA) illustrated the interaction between the DEGs and their involvement in cell signalling responses. Highly connected hubs viz. AURKA, CHEK1, CCNB2, CDC6, and MAPK14 were identified through PPI network analysis. IPA analysis showed that IL-6 and ERK5 mediated signalling pathways were highly enriched at both time points. The TP53 gene was identified to be the leading upstream regulator that directly responded to SPPV infection, resulting in downregulation at both time points. The study provides an overview of how the lamb testis genes and their underlying mechanisms link to growth and immune response during SPPV infection.

Selection and validation of suitable reference gene for qPCR gene expression analysis in lamb testis cells under Sheep pox virus infection

Virus infection alters host gene expression, therefore ideal and stable reference housekeeping genes are required to normalise the expression of other expressed host genes in quantitative real-time PCR (qRT-PCR). The suitable reference gene may vary in response to different viral infections in different hosts or cells. In the present study, we cultured primary lamb testis cells (LTC) and assessed the expression stability of seven widely used housekeeping genes (B2M, HMBS, HPRT1, HSP-90, POLR2A, 18s_RNA, GAPDH) as reference genes in Sheeppox virus (SPPV) infected and control (uninfected-0h) LTC at 0.5h, 4.0h, 8.0h, and 12.0h post-infection) using NormFinder, Bestkeeper, geNorm, and the comparative ΔCT method in RefFinder based on their expression levels. Analysis revealed that HSP90, 18s_RNA, HPRT, POLR2A, and B2M were the most stable genes from the panel in the individual analysis group in 0h, 0.5h, 4.0h, 8.0h, and 12.0h, respectively. Furthermore, B2M was shown to be the most stable reference gene in the combined control with the respective and overall infected groups, except the control group of 4.0hpi of SPPV infection. In this study, we selected the most suitable reference genes in LTC for particular time points of SPPV infection. The identified most suitable housekeeping gene can be used during normalization of expression of other targeted genes at aspecific time point of SPPV infection.

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.

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 N₆ of Adenosine called m⁶A 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 m⁶A in various types of cancers are being widely recognized. Here, we give a brief about the enzymes that install the m⁶A modification (= m⁶A writers), that remove it (= m⁶A erasers) and certain RNA binding proteins (= m⁶A readers) which affect the fate of the m⁶A-containing RNA by recruiting various proteins. We also discuss the relevance of m⁶A in ncRNAs in various cancer types, followed by a discussion on the role of m⁶A 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.