Project description:Background: The Malnad Gidda are unique dwarf Bos indicus cattle native to heavy rainfall Malnad and coastal areas of Karnataka in India. These cattle are highly adapted to harsh climatic conditions and are more resistant to Foot and Mouth disease as compared to other breeds of B.indicus. Since the first genome reference became available from B.taurus Hereford breed, only a few other breeds have been genotyped using high-throughput platforms. Also despite the known reports on high diversity within indicine breeds as compared to taurine breeds, only one draft genome of Nellore and horn transcriptome of Kankrej breed were sequenced at base level resolution. Because of the special characteristics Malnad Gidda possess, it becomes the choice of breed among many indicine cows to study at molecular level and genotyping. Results: Sequencing mRNA from the PBMCs isolated from blood of one selected Malnad Gidda bull resulted in generation of 55 million paired-end reads of 100bp length. Raw sequencing data is processed to trim the adaptor and low quality bases, and are aligned against the whole genome and transcript assemblies of Bos taurus UMD 3.1 and Bos indicus (Nellore breed) respectively. About 72% of the sequenced reads from our study could be mapped against the B.taurus genome where as only 41% of reads could be mapped against the Bos indicus transcript assembly. Transcript assembly from the alignment carried out against the annotated B.taurus UMD 3.1 genome resulted in identification of ~10,000 genes with significant expression (FPKM>1). In a similar analysis against the B.indicus Kankrej assembled transcripts we could identify only ~6,000 transcripts. From the variant analysis of the sequencing data we found ~10,000 SNPs in coding regions among which ~9,000 are novel and ~6,400 are amino acid changing. Conclusions: For the first time we have genotyped and explored the transcriptome of B.indicus Malnad Gidda breed. A comparative analysis of mapping the RNA-Seq data against the available reference genome and transcript sequences is demonstrated. An enhanced utility of transcript sequencing could be achieved by improving or completing the sequence assembly of any B.indicus breed to better characterize the indicine breeds for productivity features and selective breeding.
Project description:To understand the immune response of cows to the apicomplexan parasite Theileria annulata, we used ex vivo isolate cells derived from two infected calve : Holstein 12886 (Bos taurus), which is known to be susceptible to the disease, and Sahiwal 82H (Bos indicus), which is known to be resistant. The infected bovine macrophages of the two species with Theileria were collected and performed multiome 10X Chromium genomics scRNA-Seq.
Project description:Copy number variations (CNVs) have been demonstrated as crucial substrates for evolution, adaptation and breed formation. Chinese indigenous cattle breeds exhibit a broad geographical distribution and diverse environmental adaptability. Here, we analyzed the population structure and adaptation to high altitude of Chinese indigenous cattle based on genome-wide CNVs derived from the high-density BovineHD SNP array. We successfully detected the genome-wide CNVs of 318 individuals from 24 Chinese indigenous cattle breeds and 37 yaks as outgroups. A total of 5,818 autosomal CNV regions (683 bp - 4,477,860 bp in size), covering ~14.34% of the bovine genome (UMD3.1), were identified, showing abundant CNV resources. Neighbor-joining clustering, principal component analysis (PCA), and population admixture analysis based on these CNVs support that most Chinese cattle breeds are hybrids of Bos taurus taurus (hereinafter to be referred as Bos taurus) and Bos taurus indicus (Bos indicus). The distribution patterns of the CNVs could to some extent be related to the geographical backgrounds of the habitat of the breeds, and admixture among cattle breeds from different districts. We analyzed the selective signatures of CNVs positively involved in high-altitude adaptation using pairwise Fst analysis within breeds with a strong Bos taurus background (taurine-type breeds) and within Bos taurus×Bos indicus hybrids, respectively. CNV-overlapping genes with strong selection signatures (at top 0.5% of Fst value), including LETM1 (Fst = 0.490), TXNRD2 (Fst=0.440) and STUB1 (Fst=0.420) within taurine-type breeds, and NOXA1 (Fst = 0.233), RUVBL1 (Fst=0.222) and SLC4A3 (Fst=0.154) within hybrids, were potentially involved in the adaptation to hypoxia. Thus, we provide a new profile of population structure from the CNV aspects of Chinese indigenous cattle and new insights into high-altitude adaptation in cattle.
Project description:Bovine tropical theileriosis is a major haemoprotozoan disease associated with high rates of morbidity and mortality particularly in exotic and crossbred cattle. It is one of the major constraints for of the livestock development programmes in India and southern Asia. Indigenous cattle (Bos indicus) are less affected by this disease than exotic and crossbred cattle. Genetic basis of resistance to tropical theileriosis in indigenous cattle is not well studied. Recent studies gives an idea that differentially genes expressed in exotic and indigenous breeds play an important role in breed specific resistance to tropical theileriosis. The present study was designed to visualize the global gene expression profiling in PBMCs derived from indigenous (Tharparkar) and crossbred cattle with in vitro infection of T. annulata. T. annulata Parbhani strain, originally isolated from Maharashtra (India) and maintained as cryopreserved stabilates of ground-up tick tissue sporozoite (GUTS) of infected H. anatolicum anatolicum was used as infective material. Two separate microarray experiments were carried out using separately each for crossbred and Tharparkar cattle. The crossbred cattle showed 1082 differentially expressed genes (DEGs). Out of total DEGs, 597 genes were downregulated and 485 were upregulated. Their fold change varies from 2283.93 to -4816.02. Tharparkar cattle showed 875 differentially expressed genes. Out of total DEGs in Tharparkar cattle, 451 genes were downregulated and 424 genes were upregulated. Their fold change varies from 94.93 to -19.20. A subset of genes was validated by quantitative RT-PCR and results correlated well with data obtained from the microarrays indicating that the microarray results gave an accurate report of transcript level. Functional annotation study of differentially expressed genes has confirmed their involvement in various pathways including response to oxidative stress, immune system regulation, cell proliferation, cytoskeletal changes, kinases activity and apoptosis. Gene network analysis of these differentially expressed genes provided an effective way to understand the interaction among them. It is therefore, hypothesised that the dissimilar susceptibility to tropical theileriosis exhibited by indigenous and crossbred cattle is due to breed-specific differences in the interaction of infected cells with other immune cells, which ultimately influences the immune response generated against T. annulata infection. Global gene expression profiling in PBMCs derived from indigenous (Tharparkar) and crossbred cattle were studied after in vitro infection of T. annulata Parbhani strain at 2h time period. Two separate microarray experiments were carried out using Bovine (V2) Gene Expression Microarray, 4x44K (Agilent). Two biological replicate samples were profiled per condition (i.e. replicates samples each in crossbred and Tharparkar cattle).
Project description:Here we describe a genome-wide analysis of copy number variations (CNVs) in Chinese domestic cattle by using array comparative genomic hybridization (array CGH) and quantitative PCR (qPCR). We conducted array CGH analysis on 30 male cattle individuals, animals from consisting of 12 breeds of Bos taurus/Bos indicus, 1 Bos grunniens and and two ones of Bubalus bubalis breeds for with beef, and/or dairy or dual purpose. We identified over 470 candidate CNV regions (CNVRs) in Bos B. taurus/B. indicus; 118 candidate CNV regions (CNVRs) in B. grunniens, 139 CNVRs in B. bubalis. Furthermore, based on the Y haplotypes of B. taurus/ B. indicus, Wwe also identified 69, 337, and 251 candidate CNV regions (CNVRs) in the sub-groups of Y1, Y2 and Y3 haplotypes.
Project description:Bovine tropical theileriosis is a major haemoprotozoan disease associated with high rates of morbidity and mortality particularly in exotic and crossbred cattle. It is one of the major constraints for of the livestock development programmes in India and southern Asia. Indigenous cattle (Bos indicus) are less affected by this disease than exotic and crossbred cattle. Genetic basis of resistance to tropical theileriosis in indigenous cattle is not well studied. Recent studies gives an idea that differentially genes expressed in exotic and indigenous breeds play an important role in breed specific resistance to tropical theileriosis. The present study was designed to visualize the global gene expression profiling in PBMCs derived from indigenous (Tharparkar) and crossbred cattle with in vitro infection of T. annulata. T. annulata Parbhani strain, originally isolated from Maharashtra (India) and maintained as cryopreserved stabilates of ground-up tick tissue sporozoite (GUTS) of infected H. anatolicum anatolicum was used as infective material. Two separate microarray experiments were carried out using separately each for crossbred and Tharparkar cattle. The crossbred cattle showed 1082 differentially expressed genes (DEGs). Out of total DEGs, 597 genes were downregulated and 485 were upregulated. Their fold change varies from 2283.93 to -4816.02. Tharparkar cattle showed 875 differentially expressed genes. Out of total DEGs in Tharparkar cattle, 451 genes were downregulated and 424 genes were upregulated. Their fold change varies from 94.93 to -19.20. A subset of genes was validated by quantitative RT-PCR and results correlated well with data obtained from the microarrays indicating that the microarray results gave an accurate report of transcript level. Functional annotation study of differentially expressed genes has confirmed their involvement in various pathways including response to oxidative stress, immune system regulation, cell proliferation, cytoskeletal changes, kinases activity and apoptosis. Gene network analysis of these differentially expressed genes provided an effective way to understand the interaction among them. It is therefore, hypothesised that the dissimilar susceptibility to tropical theileriosis exhibited by indigenous and crossbred cattle is due to breed-specific differences in the interaction of infected cells with other immune cells, which ultimately influences the immune response generated against T. annulata infection.