The Complete Mitochondrial Genome of a Neglected Breed, the Peruvian Creole Cattle (Bos taurus), and Its Phylogenetic Analysis

Characterization of Complete Mitochondrial Genome of Badri Breed of Bos indicus (Bovidae: Bovinae): Selection Pressure and Comparative Analysis

High-altitude mammals are often subject to specific environmental obstacles, which exert selective pressure on their physiological and morphological traits, hence driving their evolutionary processes. It is anticipated that these circumstances will lead to the adaptive evolution of protein-coding genes (PCGs) in the mitochondrial genome, which play a crucial role in the oxidative phosphorylation system. In this study, we have generated the complete mitochondrial genome of the Badri breed of Bos indicus inhabiting a high-altitude environment to test the signatures of adaptive evolution on PCGs and their phylogenetic relationships. The complete mitogenome of the Badri breed is 16,339 bp and most tRNAs showed typical clover-leaf secondary structure with a few exceptions, like trnS1 and trnS2 without DHU arm and trnK without DHU loop. Comparative analysis of PCGs indicated that cox1 is the most conserved, while atp6 is the most variable gene. Moreover, the ratios of non-synonymous to synonymous substitution rates indicated the purifying selection (Ka/Ks < 1) in the protein-coding genes that shape the diversity in mitogenome of Bos indicus. Furthermore, Branch-site model (BSM) suggested that cox1, cox2, nad3, nad4L, and nad6 underwent stronger purifying selection (ω < 1) than other PCGs in 15 breeds of 4 species, including Badri. BSM also detected 10 positive sites in PCGs and one in 13 PCGs concatenated dataset. Additional analyses in Datamonkey indicated 11 positive sites and 23 purifying sites in the concatenated dataset, a relaxation of selection strength in nad3, and no evidence of episodic diversifying selection in any PCGs. Phylogeny revealed the sister relationship of the Badri with other breeds of Bos indicus as well as Bos frontalis (Gayal-2). The mitogenome of the Badri breed is an important genomic resource for conservation genetics of this species and also contributes to the understanding of the adaptive evolution of mitochondrial protein coding genes.

Genetic diversity and population structure of a Peruvian cattle herd using SNP data

New-generation sequencing technologies, among them SNP chips for massive genotyping, are useful for the effective management of genetic resources. To date, molecular studies in Peruvian cattle are still scarce. For the first time, the genetic diversity and population structure of a reproductive nucleus cattle herd of four commercial breeds from a Peruvian institution were determined. This nucleus comprises Brahman (N = 9), Braunvieh (N = 9), Gyr (N = 5), and Simmental (N = 15) breeds. Additionally, samples from a locally adapted creole cattle, the Arequipa Fighting Bull (AFB, N = 9), were incorporated. Female individuals were genotyped with the GGPBovine100K and males with the BovineHD. Quality control, and the proportion of polymorphic SNPs, minor allele frequency, expected heterozygosity, observed heterozygosity, and inbreeding coefficient were estimated for the five breeds. Admixture, principal component analysis (PCA), and discriminant analysis of principal components (DAPC) were performed. Also, a dendrogram was constructed using the Neighbor-Joining clustering algorithm. The genetic diversity indices in all breeds showed a high proportion of polymorphic SNPs, varying from 51.42% in Gyr to 97.58% in AFB. Also, AFB showed the highest expected heterozygosity estimate (0.41 ± 0.01), while Brahman the lowest (0.33 ± 0.01). Besides, Braunvieh possessed the highest observed heterozygosity (0.43 ± 0.01), while Brahman the lowest (0.37 ± 0.02), indicating that Brahman was less diverse. According to the molecular variance analysis, 75.71% of the variance occurs within individuals, whereas 24.29% occurs among populations. The pairwise genetic differentiation estimates (FST) between breeds showed values that ranged from 0.08 (Braunvieh vs. AFB) to 0.37 (Brahman vs. Braunvieh). Similarly, pairwise Reynold's distance ranged from 0.09 (Braunvieh vs. AFB) to 0.46 (Brahman vs. Braunvieh). The dendrogram, similar to the PCA, identified two groups, showing a clear separation between Bos indicus (Brahman and Gyr) and B. taurus breeds (Braunvieh, Simmental, and AFB). Simmental and Braunvieh grouped closely with the AFB cattle. Similar results were obtained for the population structure analysis with K = 2. The results from this study would contribute to the appropriate management, avoiding loss of genetic variability in these breeds and for future improvements in this nucleus. Additional work is needed to speed up the breeding process in the Peruvian cattle system.