Diversity analysis, runs of homozygosity and genomic inbreeding reveal recent selection in Blanco Orejinegro cattle

Selective signatures in composite MONTANA TROPICAL beef cattle reveal potential genomic regions for tropical adaptation

Genomic regions related to tropical adaptability are of paramount importance for animal breeding nowadays, especially in the context of global climate change. Moreover, understanding the genomic architecture of these regions may be very relevant for aiding breeding programs in choosing the best selection scheme for tropical adaptation and/or implementing a crossbreeding scheme. The composite MONTANA TROPICAL® population was developed by crossing cattle of four different biological types to improve production in harsh environments. Pedigree and genotype data (51962 SNPs) from 3215 MONTANA TROPICAL® cattle were used to i) characterize the population structure; ii) identify signatures of selection with complementary approaches, i.e. Integrated Haplotype Score (iHS) and Runs of Homozygosity (ROH); and iii) understand genes and traits related to each selected region. The population structure based on principal components had a weak relationship with the genetic contribution of the different biological types. Clustering analyses (ADMIXTURE) showed different clusters according to the number of generations within the composite population. Considering results of both selection signatures approaches, we identified only one consensus region on chromosome 20 (35399405-40329703 bp). Genes in this region are related to immune function, regulation of epithelial cell differentiation, and cell response to ionizing radiation. This region harbors the slick locus which is related to slick hair and epidermis anatomy, both of which are related to heat stress adaptation. Also, QTLs in this region were related to feed intake, milk yield, mastitis, reproduction, and slick hair coat. The signatures of selection detected here arose in a few generations after crossbreeding between contrasting breeds. Therefore, it shows how important this genomic region may be for these animals to thrive in tropical conditions. Further investigations on sequencing this region can identify candidate genes for animal breeding and/or gene editing to tackle the challenges of climate change.

Differential patterns in runs of homozygosity in two mice lines under divergent selection for environmental variability for birth weight

Runs of homozygosity (ROH) are defined as long continuous homozygous stretches in the genome which are assumed to originate from a common ancestor. It has been demonstrated that divergent selection for variability in mice is possible and that low variability in birth weight is associated with robustness. To analyse ROH patterns and ROH-based genomic inbreeding, two mouse lines that were divergently selected for birth weight variability for 26 generations were used, with: 752 individuals for the high variability line (H-Line), 766 individuals for the low variability line (L-Line) and 74 individuals as a reference population. Individuals were genotyped using the high density Affymetrix Mouse Diversity Genotyping Array. ROH were identified using both the sliding windows (SW) and the consecutive runs (CR) methods. Inbreeding coefficients were calculated based on pedigree (FPED ) information, on ROH identified using the SW method (FROHSW ) and on ROH identified using the CR method (FROHCR ). Differences in genomic inbreeding were not consistent across generations and these parameters did not show clear differences between lines. Correlations between FPED and FROH were high, particularly for FROHSW . Moreover, correlations between FROHSW and FPED were even higher when ROH were identified with no restrictions in the number of heterozygotes per ROH. The comparison of FROH estimates between either of the selected lines were based on significant differences at the chromosome level, mainly in chromosomes 3, 4, 6, 8, 11, 15 and 19. ROH-based inbreeding estimates that were computed using longer homozygous segments had a higher relationship with FPED . Differences in robustness between lines were not attributable to a higher homozygosis in the L-Line, but maybe to the different distribution of ROH at the chromosome level between lines. The analysis identified a set of genomic regions for future research to establish the genomic basis of robustness.

Assessment of runs of homozygosity, heterozygosity-rich regions and genomic inbreeding estimates in a subpopulation of Guzerá (Bos indicus) dual-purpose cattle

For decades, inbreeding in cattle has been evaluated using pedigree information. Nowadays, inbreeding coefficients can be obtained using genomic information such as runs of homozygosity (ROH). The aims of this study were to quantify ROH and heterozygosity-rich regions (HRR) in a subpopulation of Guzerá dual-purpose cattle, to examine ROH and HRR islands, and to compare inbreeding coefficients obtained by ROH with alternative genomic inbreeding coefficients. A subpopulation of 1733 Guzerá animals genotyped for 50k SNPs was used to obtain the ROH and HRR segments. Inbreeding coefficients by ROH (FROH ), by genomic relationship matrix based on VanRaden's method 1 using reference allele frequency in the population (FGRM ), by genomic relationship matrix based on VanRaden's method 1 using allele frequency fixed in 0.5 (FGRM_0.5 ), and by the proportion of homozygous loci (FHOM ) were calculated. A total of 15,660 ROH were identified, and the chromosome with the highest number of ROH was BTA6. A total of 4843 HRRs were identified, and the chromosome with the highest number of HRRs was BTA23. No ROH and HRR islands were identified according to established criteria, but the regions closest to the definition of an island were examined from 64 to 67 Mb of BTA6, from 36 to 37 Mb of BTA2 and from 0.50 to 1.25 Mb of BTA23. The genes identified in ROH islands have previously been associated with dairy and beef traits, while genes identified on HRR islands have previously been associated with reproductive traits and disease resistance. FROH was equal to 0.095 ± 0.084, and its Spearman correlation with FGRM was low (0.44) and moderate-high with FHOM (0.79) and with FGRM_0.5 (0.80). The inbreeding coefficients determined by ROH were higher than other cattle breeds' and higher than pedigree-based inbreeding in the Guzerá breed obtained in previous studies. It is recommended that future studies investigate the effects of inbreeding determined by ROH on the traits under selection in the subpopulation studied.

Genomic insights into the population history and adaptive traits of Latin American Criollo cattle

Criollo cattle, the descendants of animals brought by Iberian colonists to the Americas, have been the subject of natural and human-mediated selection in novel tropical agroecological zones for centuries. Consequently, these breeds have evolved distinct characteristics such as resistance to diseases and exceptional heat tolerance. In addition to European taurine (Bos taurus) ancestry, it has been proposed that gene flow from African taurine and Asian indicine (Bos indicus) cattle has shaped the ancestry of Criollo cattle. In this study, we analysed Criollo breeds from Colombia and Venezuela using whole-genome sequencing (WGS) and single-nucleotide polymorphism (SNP) array data to examine population structure and admixture at high resolution. Analysis of genetic structure and ancestry components provided evidence for African taurine and Asian indicine admixture in Criollo cattle. In addition, using WGS data, we detected selection signatures associated with a myriad of adaptive traits, revealing genes linked to thermotolerance, reproduction, fertility, immunity and distinct coat and skin coloration traits. This study underscores the remarkable adaptability of Criollo cattle and highlights the genetic richness and potential of these breeds in the face of climate change, habitat flux and disease challenges. Further research is warranted to leverage these findings for more effective and sustainable cattle breeding programmes.

Characterization of runs of Homozygosity revealed genomic inbreeding and patterns of selection in indigenous sahiwal cattle

Runs of homozygosity (ROH) are contiguous genomic regions, homozygous across all sites which arise in an individual due to the parents transmitting identical haplotypes to their offspring. The genetic improvement program of Sahiwal cattle after decades of selection needs re-assessment of breeding strategy and population phenomena. Hence, the present study was carried out to optimize input parameters in PLINK for ROH estimates, to explore ROH islands and assessment of pedigree and genome-based inbreeding in Sahiwal cattle. The sliding window approach with parameters standardized to define ROH for the specific population under study was used for the identification of runs. The optimum maximum gap, density, window-snp and window-threshold were 250 Kb, 120 Kb/SNP, 10, 0.05 respectively and ROH patterns were also characterized. ROH islands were defined as the short homozygous genomic regions shared by a large proportion of individuals in a population, containing significantly higher occurrences of ROH than the population specific threshold level. These were identified using the -homozyg-group function of the PLINK v1.9 program. Our results indicated that the Islands of ROH harbor a few candidate genes, ACAD11, RFX4, BANP, UBA5 that are associated with major economic traits. The average FPED (Pedigree based inbreeding coefficient), FROH (Genomic inbreeding coefficient), FHOM (Inbreeding estimated as the ratio of observed and expected homozygous genotypes), FGRM (Inbreeding estimated on genomic relationship method) and FGRM0.5 (Inbreeding estimated from the diagonal of a GRM with allele frequencies near to 0.5) were 0.009, 0.091, 0.035, -0.104 and -0.009, respectively. Our study revealed the optimum parameter setting in PLINK viz. maximal gaps between two SNPs, minimal density of SNPs in a segment (in kb/SNP) and scanning window size to identify ROH segments, which will enable ROH estimation more efficient and comparable across various SNP genotyping-based studies. The result further emphasized the significant role of genomics in unraveling population diversity, selection signatures and inbreeding in the ongoing Sahiwal breed improvement programs.

Genetic Architecture and Signatures of Selection in the Caqueteño Creole (Colombian Native Cattle)

Evolutionary mechanisms have shaped the genomic architecture of Colombian Creole cattle breeds. The mating and selection processes have impacted several traits, promoting differences within and between populations. Studies of population structure and selection signatures in Colombian Creole breeds are scarce, and need more attention to better understand genetic differentiation, gene flow, and genetic distance. This study aimed to analyze the population structure and identify selection imprints in the Criollo Caqueteño (CAQ) population. It used 127 CAQ animals genotyped with Chip HD 777,000 SNPs. The population structure analyses used discriminant principal component analysis (DAPC), integrated haplotype scoring (iHS), and index-fixing (Fst) methodologies to detect selection signals. We can highlight SNP regions on the genes TMPRSS15, PGAM2, and EGFR, identified by the Fst method. Additionally, the iHS regions for cluster 1 identified candidate genes on BTA 3 (CMPK1 and FOXD2), BTA 11 (RCAN1), and BTA 22 (ARPP21). In group 2, we can highlight the genes on BTA 4 (SLC13A4, BRAF), BTA 9 (ULBP), BTA 14 (CSMD3) and BTA 19 (KRTAP9-2). These candidate genes have been associated with fertility traits, precocity, growth, and environmental and disease resistance, indicating a genetic potential in CAQ animals. All this promotes a better understanding of the diversity and genetic structure in the CAQ population. Based on that, our study can significantly assist the sustainable development and conservation of the breed in the Colombian Amazon.

Identification of Runs of Homozygosity Islands and Genomic Estimated Inbreeding Values in Caqueteño Creole Cattle (Colombia)

The Caqueteño Creole (CAQ) is a native breed of cattle from the Caquetá department (Colombia), adapted to tropical conditions, which is extremely important to production systems in those regions. However, CAQ is poorly studied. In this sense, population structure studies associated with runs of homozygosity (ROH) analysis would allow for a better understanding of CAQ. Through ROH analysis, it is possible to reveal genetic relationships between individuals, measure genome inbreeding levels, and identify regions associated with traits of economic interest. Samples from a CAQ population (n = 127) were genotyped with the Bovine HD BeadChip (777,000 SNPs) and analyzed with the PLINK 1.9 program to estimate FROH and ROH islands. We highlighted a decrease in inbreeding frequency for FROH 4−8 Mb, 8−16 Mb, and >16 Mb classes, indicating inbreeding control in recent matings. We also found genomic hotspot regions on chromosomes 3, 5, 6, 8, 16, 20, and 22, where chromosome 20 harbored four hotspots. Genes in those regions were associated with fertility and immunity traits, muscle development, and environmental resistance, which may be present in the CAQ breed due to natural selection. This indicates potential for production systems in tropical regions. However, further studies are necessary to elucidate the CAQ production objective.

Frequency of genotypic markers for genetic disorders, colour, polledness, and major genes in Blanco Orejinegro cattle

The purpose of this work was to establish DNA marker frequencies for genetic disorders, colour, horned/polled trait, and major genes of importance for productive and reproductive traits in Blanco Orejinegro (BON) cattle. The Blanco Orejinegro breed is a Colombian creole breed characterized by a white hair coat on black skin with black ears, black hair on the middle part of the legs, and absence of horns. We genotyped 420 animals of Colombia, 70 with the GGP-Bovine 150 K chip and 350 with the GGP-Bovine 50 K chip. The markers were associated with 50 genetic diseases, 52 major gene variants related to productive traits, and 12 variants related to coat coloration, presence of horns, and adaptation, selected from the information contained in the chips. Genotype frequencies were estimated using the R statistical program. Genetic disorder annotations were derived using the Online Mendelian Inheritance in Animals tool (OMIA) and the average inbreeding coefficient (F) (n = 7799) using the MTDFNRM program. Carriers were found for 16 of the genetic disorders evaluated but with low frequencies (0.24 to 2.46%); no homozygous animals were found for the disorders. Carriers were associated with disorders such as bovine leukocyte adhesion deficiency (BLAD), deficiency of uridine monophosphate synthase (DUMPS), syndactyly, and epidermolysis bullosa (EB). The F was 4.41%. Concerning the genes associated with colour (TYR, MC1R, and PMEL), alleles related to black pigmentation, the absence of horns (polledness), and slick coat (an adaptive trait) were highly frequent (> 81.90%). Markers associated with milk production and quality, yellow fat, and fertility showed variable frequencies, indicating selection potential. Allele frequency of genetic disorders in BON cattle was low, suggesting few genetic disorder problems, with syndactyly being the most frequent condition. The markers associated with colour and polledness were almost fixed, with a frequency at or near 100%. Production and reproduction markers showed variability for selection.

Genomic characterization of the nucleus for conservation of the Chino Santandereano breed using SNP markers

The Chino Santandereano (CHS) creole cattle breed has undergone an adaptation process of more than 500 years to the Colombian mountainous tropics. Despite its advantages, the breed has been replaced by specialized cattle that put it at risk of extinction. The aim of this work was to estimate the genomic diversity and population structure of the nucleus of conservation of the CHS breed. Thirty-seven CHS and 20 Brahman animals were genotyped due to the possible introgression of genes with the GGP Bovine LD v3 chip. Quality control was performed, and linkage disequilibrium (LD), effective population size (Ne), ROH segments, homozygosity, and genomic inbreeding in the breed were estimated. Subsequently, 50 K genomic information of the Holstein (n = 30) and Romosinuano (n = 8) breeds were included to estimate the minor allele frequency (MAF) with common markers and constructing the graphs of the principal component analysis (PCA). Pairwise FSTs were estimated and a neighbor-joining tree was constructed using the IBS matrix. Admixture was used with k = 2 to 10 for the racial composition. LD (r²) was found up to a distance of 0.13 Mb, r² > 0.3 at a distance of 340.3 kb, and Ne of 32 ± 1. ROH inbreeding was 5.36 ± 0.86%, with a higher contribution from recent inbreeding (4.55%). The PCA showed that the creole breeds were closer together, and the Brahman was more distant. The admixture analysis suggested k = 5 possible ancestral groups and shows that within the CHS breed, there seem to be two different groups with little Holstein and Brahman introgression. The genetic diversity parameters obtained in this work show minimal diversity in this breed and reinforce the need to protect this resource and the conservation banks.