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:Environmental heat stress in dairy cattle leads to poor health, reduced milk production and decreased reproductive efficiency. Multiple genes interact and coordinate the response to overcome the impact of heat stress. The present study identified heat shock regulated genes in the peripheral blood mononuclear cells (PBMC). Genome-wide expression patterns for cellular stress response were compared between two genetically distinct groups of cattle viz., Hariana (B. indicus) and Vrindavani (B. indicus X B. taurus). In addition to major heat shock response genes, oxidative stress and immune response genes were also found to be affected by heat stress. Heat shock proteins such as HSPH1, HSPB8, FKB4, DNAJ4 and SERPINH1 were up-regulated at higher fold change in Vrindavani compared to Hariana cattle. The oxidative stress response genes (HMOX1, BNIP3, RHOB and VEGFA) and immune response genes (FSOB, GADD45B and JUN) were up-regulated in Vrindavani whereas the same were down-regulated in Hariana cattle. The enrichment analysis of dysregulated genes revealed the biological functions and signaling pathways that were affected by heat stress. Overall, these results show distinct cellular responses to heat stress in two different genetic groups of cattle. This also highlight the long-term adaptation of B. indicus (Hariana) to tropical climate as compared to the crossbred (Vrindavani) with mixed genetic makeup (B. indicus X B. taurus).
Project description:Epigenetic variations result from long-term adaptation to environmental factors. The Bos indicus (zebu) adapted to tropical conditions, whereas Bos taurus adapted to temperate conditions; hence native zebu cattle and its crossbred (B indicus × B taurus) show differences in responses to heat stress. The present study evaluated genome-wide DNA methylation profiles of these two breeds of cattle that may explain distinct heat stress responses. Physiological responses to heat stress and estimated values of Iberia heat tolerance coefficient and Benezra's coefficient of adaptability revealed better relative thermotolerance of Hariana compared to the Vrindavani cattle. Genome-wide DNA methylation patterns were different for Hariana and Vrindavani cattle. The comparison between breeds indicated the presence of 4599 significant differentially methylated CpGs with 756 hypermethylated and 3845 hypomethylated in Hariana compared to the Vrindavani cattle. Further, we found 79 genes that showed both differential methylation and differential expression that are involved in cellular stress response functions. Differential methylations in the microRNA coding sequences also revealed their functions in heat stress responses. Taken together, epigenetic differences represent the potential regulation of long-term adaptation of Hariana (B indicus) cattle to the tropical environment and relative thermotolerance.
Project description:BackgroundReference genomes are essential in the analysis of genomic data. As the cost of sequencing decreases, multiple reference genomes are being produced within species to alleviate problems such as low mapping accuracy and reference allele bias in variant calling that can be associated with the alignment of divergent samples to a single reference individual. The latest reference sequence adopted by the scientific community for the analysis of cattle data is ARS_UCD1.2, built from the DNA of a Hereford cow (Bos taurus taurus-B. taurus). A complementary genome assembly, UOA_Brahman_1, was recently built to represent the other cattle subspecies (Bos taurus indicus-B. indicus) from a Brahman cow haplotype to further support analysis of B. indicus data. In this study, we aligned the sequence data of 15 B. taurus and B. indicus breeds to each of these references.ResultsThe alignment of B. taurus individuals against UOA_Brahman_1 detected up to five million more single-nucleotide variants (SNVs) compared to that against ARS_UCD1.2. Similarly, the alignment of B. indicus individuals against ARS_UCD1.2 resulted in one and a half million more SNVs than that against UOA_Brahman_1. The number of SNVs with nearly fixed alternative alleles also increased in the alignments with cross-subspecies. Interestingly, the alignment of B. taurus cattle against UOA_Brahman_1 revealed regions with a smaller than expected number of counts of SNVs with nearly fixed alternative alleles. Since B. taurus introgression represents on average 10% of the genome of Brahman cattle, we suggest that these regions comprise taurine DNA as opposed to indicine DNA in the UOA_Brahman_1 reference genome. Principal component and admixture analyses using genotypes inferred from this region support these taurine-introgressed loci. Overall, the flagged taurine segments represent 13.7% of the UOA_Brahman_1 assembly. The genes located within these segments were previously reported to be under positive selection in Brahman cattle, and include functional candidate genes implicated in feed efficiency, development and immunity.ConclusionsWe report a list of taurine segments that are in the UOA_Brahman_1 assembly, which will be useful for the interpretation of interesting genomic features (e.g., signatures of selection, runs of homozygosity, increased mutation rate, etc.) that could appear in future re-sequencing analysis of indicine cattle.