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: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.
Project description:BackgroundThe incidence of poor semen quality and sub-fertility/infertility is higher in crossbred as compared to Zebu males. Several attempts have been made to understand the possible reasons for higher incidence of fertility problems in crossbred males, at sperm phenotype, proteome and genome level but with variable results. Since the quality of the ejaculated spermatozoa is determined by the testicular environment, assessing the testicular transcriptome between these breeds would help in identifying the possible mechanisms associated with infertility in crossbred bulls. However, such information is not available. We performed global transcriptomic profiling of testicular tissue from crossbred and Zebu bulls using Agilent Bos taurus GXP 8X60k AMADID: 29411 array. To the best of our knowledge, this is the first study comparing the testicular mRNAs between crossbred and Zebu bulls.ResultsOut of the 14,419 transcripts detected in bovine testis, 1466 were differentially expressed between crossbred and Zebu bulls, in which 1038 were upregulated and 428 were downregulated in crossbred bulls. PI4KB and DPY19L2 genes, reported to be involved in sperm capacitation and acrosome formation respectively, were among the top 10 downregulated transcripts in crossbred testis. Genes involved in ubiquitination and proteolysis were upregulated, while genes involved in cell proliferation, stem cell differentiation, stem cell population maintenance, steroidogenesis, WNT signalling, protein localization to plasma membrane, endocannabinoid signalling, heparin binding, cAMP metabolism and GABA receptor activity were downregulated in crossbred testis. Among the 10 genes validated using qPCR, expression of CCNYL, SOX2, MSMB, SPATA7, TNP1, TNP2 and CRISP2 followed the same trend as observed in microarray analysis with SPATA7 being significantly downregulated and transition proteins (TNP1, TNP2) being significantly upregulated in crossbred bulls.ConclusionsAbundant proteolysis by ubiquitination and downregulation of WNT signaling, cell proliferation, differentiation and steroidogenesis might be associated with higher incidence of poor semen quality and/or sub-fertility/infertility in crossbred bulls as compared to Zebu bulls. Downregulation of SPATA7 (Spermatogenesis Associated 7) and upregulation of transition proteins (TNP1 and TNP2) in crossbred bull testis might be associated with impaired spermatogenesis processes including improper chromatin compaction in crossbred bulls.
Project description:BackgroundBos taurus and Bos indicus are two main sub-species of cattle. However, the differential copy number variations (CNVs) between them are not yet well studied.ResultsBased on the new high-quality cattle reference genome ARS-UCD1.2, we identified 13,234 non-redundant CNV regions (CNVRs) from 73 animals of 10 cattle breeds (4 Bos taurus and 6 Bos indicus), by integrating three detection strategies. While 6990 CNVRs (52.82%) were shared by Bos taurus and Bos indicus, large CNV differences were discovered between them and these differences could be used to successfully separate animals into two subspecies. We found that 2212 and 538 genes uniquely overlapped with either indicine-specific CNVRs and or taurine-specific CNVRs, respectively. Based on FST, we detected 16 candidate lineage-differential CNV segments (top 0.1%) under selection, which overlapped with eight genes (CTNNA1, ENSBTAG00000004415, PKN2, BMPER, PDE1C, DNAJC18, MUSK, and PLCXD3). Moreover, we obtained 1.74 Mbp indicine-specific sequences, which could only be mapped on the Bos indicus reference genome UOA_Brahman_1. We found these sequences and their associated genes were related to heat resistance, lipid and ATP metabolic process, and muscle development under selection. We further analyzed and validated the top significant lineage-differential CNV. This CNV overlapped genes related to muscle cell differentiation, which might be generated from a retropseudogene of CTH but was deleted along Bos indicus lineage.ConclusionsThis study presents a genome wide CNV comparison between Bos taurus and Bos indicus. It supplied essential genome diversity information for understanding of adaptation and phenotype differences between the Bos taurus and Bos indicus populations.