Population diversity, admixture, and demographic trend of the Sumba Ongole cattle based on genomic data

OBJECTIVE

Sumba Ongole (SO) cattle are valuable breed due to their important role in the development of Indonesian cattle. Despite rapid advances in molecular technology, no genomic studies on SO cattle have been conducted to date. The aim of this study is to provide genomic profile related to the population diversity, admixture, and demographic trends of SO cattle.

METHODS

Genomic information was gathered from 79 SO cattle using the Illumina Bovine SNP50 v3 Beadchip, and for comparative purposes, additional genotypes from 209 cattle populations worldwide were included. The expected and observed heterozygosity, inbreeding coefficient, pairwise fixation indices between-population, and Nei's genetic distance were examined. Multidimensional scaling, admixture, and treemix analyses were used to investigate the population structure. Based on linkage disequilibrium and effective population size calculations, the demographic trend was observed.

RESULTS

The findings indicated that the genetic diversity of SO cattle was similar to that of other indicine breeds. SO cattle were genetically related to indicines but not to taurines or Bali cattle. The study further confirmed the close relationship between SO, Ongole, and Nellore cattle. Additionally, a small portion of the Ongole mixture were identified dominant in the SO population at the moment. The study also discovered that SO and Bali cattle (Bos javanicus) could have been ancestors in the development of Ongole Grade cattle, which corresponds to the documented history of Ongolization. Our finding indicate that SO cattle have maintained stability and possess unique traits separate from their ancestors.

CONCLUSION

In conclusion, the genetic diversity of the SO cattle has been conserved as a result of the growing significance of the present demographic trend. Consistent endeavors are necessary to uphold the fitness of the breed.

Morphometric and physical characteristics of Indonesian beef cattle

To assess their potential for beef production, this study describes the morphometrics, physical traits, and body weight of 1034 Indonesian beef cattle from eight breeds, namely Bali, Rambon, Madura, Ongole Grade, Kebumen Ongole Grade, Sasra, Jabres, and Pasundan. An analysis of variance in addition to cluster, Euclidean distance, dendrogram, discriminant function, stepwise linear regression, and morphological index analyses were performed to describe the differences in traits among breeds. The morphometric proximity analysis revealed two distinct clusters with a common ancestor, where the first cluster included Jabres, Pasundan, Rambon, Bali, and Madura cattle and the second included Ongole Grade, Kebumen Ongole Grade, and Sasra cattle, with an average suitability value of 93.20 %. This showed that the classification and validation methods can be used to distinguish breeds. The most important factor in estimating body weight was the heart girth circumference. Ongole Grade cattle had the highest cumulative index, followed by Sasra, Kebumen Ongole Grade, Rambon, and Bali cattle. A cumulative index value > 3 can be used as a threshold for determining the type and function of beef cattle.

Comparison of growth curve models for Ongole Grade cattle

Growth curves are important for understanding the growth phase and specific characteristics of farm animals. This study evaluated growth curve models to predict the age at puberty and growth rate of Ongole Grade cattle. The study used the weights of cattle from birth to 5 years old in 2011 to 2019 obtained from the Beef Cattle Research Station. The data were analyzed using four non-linear growth curve models: Brody, von Bertalanffy, logistic, and Gompertz. The coefficients of determination (R²) of the growth equation for the respective models were in the high range. The lowest Akaike's information criterion (AIC) value for males and females was obtained by the Brody and logistic models, respectively. Males with the von Bertalanffy, logistic, and Gompertz models had an estimated age and body weight at puberty of 12.56 months and 198.93 kg, 21.46 months and 298.43 kg, and 15.71 months and 235.49 kg, respectively, while females had an estimated age and body weight of 11.79 months and 166.44 kg, 13.94 months and 164.08 kg, and 11.76 months and 169.57 kg, respectively. The Brody function best fit the growth rate and exhibited the smallest mean deviation.

Association of pleomorphic adenoma gene 1 with body weight and measurement of Bali cattle (Bos javanicus)

Background and Aim:

Pleomorphic adenoma gene 1 (PLAG1) encodes a multifunctional transcription factor that controls many genes and pathways and is associated with cattle body weight and measurements. This study aimed to evaluate the association between PLAG1 polymorphisms with body weight and measurements in Bali cattle.

Materials and Methods:

A total of 87 Bali cattle, consisting of 48 bulls and 39 heifers at the Breeding Center for Bali Cattle, were used as the population in this study. Cattle were 2 years old and kept semi-intensively in the pasture. Phenotype data consisting of body weight, withers height, body length, chest girth, waist height, and chest depth were measured. Birth weight data were obtained from birth records, and weight gain, adjusted weaning weight, and yearling weight were calculated using formulas. Blood samples were taken from the jugular vein as much as 5 mL, and genomic DNA was isolated using the salting-out method. Polymerase chain reaction (PCR) was performed to amplify three target polymorphisms, namely, g.48308 C>T, g.32212 (19 bp indel), and g.45233 T>C. The presence of a 19 bp indel was determined by direct observation of the PCR product on a 2% agarose gel. Two other polymorphisms were detected by PCR-restriction fragment length polymorphism using the restriction endonuclease enzymes SacII and BclI. PLAG1 genotype and phenotype associations were analyzed using a general linear model.

Results:

The results showed that two of the target polymorphisms in PLAG1 did not vary. The DD genotype indicated by 123 bp of PCR product was the only genotype identified for g.32212 19 bp indel, and TT genotype was the only genotype found for g.45233 T>C single-nucleotide polymorphism (SNP). Conversely, g.48308 C>T SNP was found to be polymorphic. In addition, the g.48308 C>T polymorphism of PLAG1 was significantly associated with body length of Bali cattle. Cattle with the CC genotype had a greater body length than the other two genotypes.

Conclusion:

The g.48308 C>T SNP in PLAG1 was associated with Bali cattle body length characteristics. This finding could be used as a basis for selecting Bali cattle based on body length characteristics.

Genomic structure of Bali cattle based on linkage disequilibrium and effective population size analyses using 50K single nucleotide polymorphisms data

Background and aim:

Bali Cattle (Bos j. javanicus) is a local breed originating in Indonesia, accounting for 32.3% of the total cattle population. To date, no studies of the genetic structure and demographic status of Bali cattle have been conducted, even though the breeding of Bali cattle has a long and unique history that is likely to have impacted its genetic diversity. Therefore, a study that used molecular breeding technologies to characterize the demography of Bali cattle would be timely. This study aimed to examine genome diversity in Bali cattle and estimate the linkage disequilibrium (LD) and effective population size (N_e) values in the cattle population.

Materials and Methods:

In this study, we explored the population structure and genetic diversity of Bali cattle using genomic-level analyses. Our study primarily studied cattle that had been bred in livestock breeding centers since these breeds had subsequently spread throughout Indonesia. We focused on characterizing the genetic structure, determining the level of LD present, and estimating the N_e of the Bali cattle population. The genomic data used for this study were obtained from DNA samples of 48 Bali cattle collected at the Breeding Center of Bali Cattle as well as 54 genomic samples from Bali cattle collected elsewhere in Indonesia that had been used in recent publications. This genomic dataset included exclusively 50K single nucleotide polymorphisms (SNP) array (Illumina Bovine 50SNP bead chip, Illumina, USA) data.

Results:

We found that the LD values of Bali cattle from the breeding center and those raised elsewhere were 0.48±0.43 and 0.39±0.40, respectively. Subsequently, the N_e value of Bali cattle from the breeding center and farmers was 151 and 96, respectively.

Conclusion:

Our results suggest that the selection program of the breeding center is beneficial for maintaining the genetic diversity of Bali cattle.

An insight into the evolutionary history of Indonesian cattle assessed by whole genome data analysis

The domestication of Indonesian cattle was investigated through a study of their genetic diversity, up to the genome level. Little documentation exists regarding the history of domestication of Indonesian cattle and questions remain despite a growing body of molecular evidence. In this study, we genotyped seven Indonesian cattle breeds using an Illumina BovineSNP50 Bead Chip to provide insight into their domestication and demographic history in a worldwide population context. Our analyses indicated the presence of hybrid cattle, with Bos javanicus and Bos indicus ancestries being most prevalent, as well as purebred cattle. We revealed that all the breeds were interconnected through several migration events. However, their demographic status varied widely. Although almost all the Indonesian cattle had an effective population size higher than the minimum level required to ensure breed fitness, efforts are still needed to maintain their genetic variability and purity.

Estimation of linkage disequilibrium and analysis of genetic diversity in Korean chicken lines

The development of genetic markers for animal breeding is an effective strategy to reduce the time and cost required to improve economically important traits. To implement genomic selection in the multibreed chicken population of Korea, an understanding of the linkage disequilibrium (LD) status of the target population is essential. In this study, we performed population genetic analyses to investigate LD decay, the effective population size, and breed diversity using 600K high-density single nucleotide polymorphism genotypes of 189 native chickens in 14 lines (including Korean native chicken, imported and adapted purebred and commercial chickens). The results indicated that commercial native chickens have less calculated LD (average, r2 = 0.13-0.26) and purebred native chickens have more calculated LD (average, r2 = 0.24-0.37) across the entire genome. The effective population sizes of the examined lines showed patterns opposite to those of population LD. The phylogeny and admixture analyses showed that commercial and purebred chickens were well distinguished, except for Rhode Island Red (RIR) purebred lines of NC (NIAS_RIR_C) and ND (NIAS_RIR_D). These lines are difficult to distinguish clearly because they originated from the same respective breeds. The results of this study may provide important information for the development of genetic markers that can be used in breeding to improve the economic traits of native chickens.

Genome-wide linkage disequilibrium and past effective population size in three Korean cattle breeds

The routine collection and use of genomic data are useful for effectively managing breeding programs for endangered populations. Linkage disequilibrium (LD) using high-density DNA markers has been widely used to determine population structures and predict the genomic regions that are associated with economic traits in beef cattle. The extent of LD also provides information about historical events, including past effective population size (N_e ), and it allows inferences on the genetic diversity of breeds. The objective of this study was to estimate the LD and N_e in three Korean cattle breeds that are genetically similar but have different coat colors (Brown, Brindle and Jeju Black Hanwoo). Brindle and Jeju Black are endangered breeds with small populations, whereas Brown Hanwoo is the main breeding population in Korea. DNA samples from these cattle breeds were genotyped using the Illumina BovineSNP50 Bead Chip. We examined 13 cattle breeds, including European taurines, African taurines and indicines, and hybrids to compare their LD values. Brown Hanwoo consistently had the lowest mean LD compared to Jeju Black, Brindle and the other 13 cattle breeds (0.13, 0.19, 0.21 and 0.15-0.22 respectively). The high LD values of Brindle and Jeju Black contributed to small N_e values (53 and 60 respectively), which were distinct from that of Brown Hanwoo (531) for 11 generations ago. The differences in LD and N_e for each breed reflect the breeding strategy applied. The N_e for these endangered cattle breeds remain low; thus, effort is needed to bring them back to a sustainable tract.

Validation of Single Nucleotide Polymorphisms Associated with Carcass Traits in a Commercial Hanwoo Population

Four carcass traits, namely carcass weight (CW), eye muscle area (EMA), back fat thickness (BF), and marbling score (MS), are the main price decision parameters used for purchasing Hanwoo beef. The development of DNA markers for these carcass traits for use in a beef management system could result in substantial profit for beef producers in Korea. The objective of this study was to validate the association of highly significant single nucleotide polymorphisms (SNPs) identified in a previous genome-wide association study (GWAS) with the four carcass traits in a commercial Hanwoo population. We genotyped 83 SNPs distributed across all 29 autosomes in 867 steers from a Korean Hanwoo feedlot. Six SNPs, namely ARS-BFGL-NGS-22774 (Chr4, Pos:4889229), ARS-BFGL-NGS-100046 (Chr6, Pos:61917424), ARS-BFGL-NGS-39006 (Chr27, Pos:38059196), ARS-BFGL-NGS-18790 (Chr10, Pos:26489109), ARS-BFGL-NGS-43879 (Chr9, Pos:39964297), and BTB-00775794 (Chr20, Pos:20476265), were found to be associated with CW, EMA, BF, and MS. The ARS-BFGL-NGS-22774, BTB-00775794, and ARS-BFGL-NGS-39006 markers accounted for 1.80%, 1.72%, and 1.35% (p<0.01), respectively, of the phenotypic variance in the commercial Hanwoo population. Many genes located in close proximity to the significant SNPs identified in this study were previously reported to have roles in carcass traits. The results of this study could be useful for marker-assisted selection programs.

Stories and Challenges of Genome Wide Association Studies in Livestock - A Review

Undoubtedly livestock is one of the major contributors to the economy of any country. The economic value of livestock includes meat, dairy products, fiber, fertilizer etc. Understanding and identifying the associations of quantitative trait loci (QTL) with the economically important traits is believed to substantially benefit the livestock industry. The past two decades have seen a flurry of interest in mapping the QTL associated with traits of economic importance on the genome. With the availability of single nucleotide polymorphism chip of various densities it is possible to identify regions, QTL and genes on the genome that explain the association and its effect on the phenotype under consideration. Remarkable advancement has been seen in genome wide association studies (GWAS) since its inception till the present day. In this review we describe the progress and challenges of GWAS in various livestock species.