Genomics Confirm an Alarming Status of the Genetic Diversity of Belgian Red and Belgian White Red Cattle

Genetic diversity is increasingly important for researchers and society. Small and local populations deserve more attention especially, as they may harbor important characteristics. Moreover, small populations are at greater risk and their genetic management is often more challenging. Likewise, European red cattle populations are threatened, as they are outcompeted by more specialized cattle breeds. In this study, we investigate the genetic diversity of two local Belgian red cattle breeds: Belgian Red and Belgian White Red cattle. A total of 270 animals were genotyped via medium density SNP arrays. Genetic diversity was assessed using runs of homozygosity screening, effective population size estimation and Fst analyses. Genomic inbreeding coefficients based on runs of homozygosity were estimated at 7.0% for Belgian Red and 6.1% for Belgian White Red cattle, and both populations had a low effective population size (68 and 86, respectively). PCA, Fst and admixture analyses revealed the relationship to 52 other international breeds, where they were closest related to some Belgian, French, Scandinavian and Dutch breeds. Moreover, Fst analyses revealed for Belgian Red cattle a signature of selection on BTA6, adjacent to the KIT gene. This study gains important knowledge on the genetic diversity of these two small local red cattle breeds, and will aid in their (genetic) management.

Runs of homozygosity in the Italian goat breeds: impact of management practices in low-input systems

Background Climate and farming systems, several of which are considered as low-input agricultural systems, vary between goat populations from Northern and Southern Italy and have led to different management practices. These processes have impacted genome shaping in terms of inbreeding and regions under selection and resulted in differences between the northern and southern populations. Both inbreeding and signatures of selection can be pinpointed by the analysis of runs of homozygosity (ROH), which provides useful information to assist the management of this species in different rural areas. Results We analyzed the ROH distribution and inbreeding (FROH) in 902 goats from the Italian Goat Consortium2 dataset. We evaluated the differences in individual ROH number and length between goat breeds from Northern (NRD) and Central-southern (CSD) Italy. Then, we identified the signatures of selection that differentiate these two groups using three methods: ROH, ΔROH, and averaged FST. ROH analyses showed that some Italian goat breeds have a lower inbreeding coefficient, which is attributable to their management and history. ROH are longer in breeds that are undergoing non-optimal management or with small population size. In several small breeds, the ROH length classes are balanced, reflecting more accurate mating planning. The differences in climate and management between the NRD and CSD groups have resulted in different ROH lengths and numbers: the NRD populations bred in isolated valleys present more and shorter ROH segments, while the CSD populations have fewer and longer ROH, which is likely due to the fact that they have undergone more admixture events during the horizontal transhumance practice followed by a more recent standardization. We identified four genes within signatures of selection on chromosome 11 related to fertility in the NRD group, and 23 genes on chromosomes 5 and 6 related to growth in the CSD group. Finally, we identified 17 genes on chromosome 12 related to environmental adaptation and body size with high homozygosity in both groups. Conclusions These results show how different management practices have impacted the level of genomic inbreeding in two Italian goat groups and could be useful to assist management in a low-input system while safeguarding the diversity of small populations.

Correlation of Genomic and Pedigree Inbreeding Coefficients in Small Cattle Populations

This study aimed to identify inbreeding coefficient (F) estimators useful for improvement programs in a small Holstein population through the evaluation of different methodologies in the Mexican Holstein population. F was estimated as follows: (a) from pedigree information (Fped); (b) through runs of homozygosity (Froh); (c) from the number of observed and expected homozygotic SNP in the individuals (Fgeno); (d) through the genomic relationship matrix (Fmg). The study included information from 4277 animals with pedigree records and 100,806 SNP. The average and standard deviation values of F were 3.11 ± 2.30 for Fped, −0.02 ± 3.55 for Fgeno, 2.77 ± 0.71 for Froh and 3.03 ± 3.05 for Fmg. The correlations between coefficients varied from 0.30 between Fped and Froh, to 0.96 between Fgeno and Fmg. Differences in the level of inbreeding among the parent’s country of origin were found regardless of the method used. The correlations among genomic inbreeding coefficients were high; however, they were low with Fped, so further research on this topic is required.

Inbreeding Calculated with Runs of Homozygosity Suggests Chromosome-Specific Inbreeding Depression Regions in Line 1 Hereford

The goal of this study was to evaluate inbreeding in a closed beef cattle population and assess phenotype prediction accuracy using inbreeding information. Effects of inbreeding on average daily gain phenotype in the Line 1 Hereford cattle population were assessed in this study. Genomic data were used to calculate inbreeding based on runs of homozygosity (ROH), and pedigree information was used to calculate the probability of an allele being identical by descent. Prediction ability of phenotypes using inbreeding coefficients calculated based on pedigree information and runs of homozygosity over the whole genome was close to 0, even in the case of significant inbreeding coefficient effects. On the other hand, inbreeding calculated per individual chromosomes’ ROH yielded higher accuracies of prediction. Additionally, including only ROH from chromosomes with higher predicting ability further increased prediction accuracy. Phenotype prediction accuracy, inbreeding depression, and the effects of chromosome-specific ROHs varied widely across the genome. The results of this study suggest that inbreeding should be evaluated per individual regions of the genome. Moreover, mating schemes to avoid inbreeding depression should focus more on specific ROH with negative effects. Finally, using ROH as added information may increase prediction of the genetic merit of animals in a genomic selection program.

COVID-19: Impact on linguistic and genetic isolates of India

The rapid expansion of coronavirus SARS-CoV-2 has impacted various ethnic groups all over the world. The burden of infectious diseases including COVID-19 are generally reported to be higher for the Indigenous people. The historical knowledge have also suggested that the indigenous populations suffer more than the general populations in the pandemic. Recently, it has been reported that the indigenous groups of Brazil have been massively affected by COVID-19. Series of studies have shown that many of the indigenous communities reached at the verge of extinction due to this pandemic. Importantly, South Asia also has several indigenous and smaller communities, that are living in isolation. Till date, despite the two consecutive waves in India, there is no report on the impact of COVID-19 for indigenous tribes. Since smaller populations experiencing drift may have greater risk of such pandemic, we have analysed Runs of Homozygosity (ROH) among South Asian populations and identified several populations with longer homozygous segments. The longer runs of homozygosity at certain genomic regions may increases the susceptibility for COVID-19. Thus, we suggest extreme careful management of this pandemic among isolated populations of South Asia.

Runs of homozygosity derived from pool-seq data reveal fine-scale population structures in Western honey bees (Apis mellifera)

Runs of homozygosity (ROH) are continuous homozygous segments that arise through the transmission of haplotypes that are identical by descent (IBD). The length and distribution of ROH segments provide insight into the genetic diversity of populations and are useful to detect selection signatures. Here, we analysed pooled whole-genome sequencing data from 265 Western honey bee colonies from the two subspecies Apis mellifera mellifera and Apis mellifera carnica. Integrating individual ROH patterns and admixture levels in a high-resolution population network visualization allowed us to ascertain major differences between the two subspecies. Within A. m. mellifera, we identified well-defined substructures according to the genetic origin of the colonies and a fair amount of admixed colonies, despite the current applied conservation efforts. In contrast, A. m. carnica colonies were more inbred and could not be differentiated according to the geographical origin. We identified 29 coding genes in overlapping ROH segments within the two subspecies. Genes embedded in A. m. carnica specific homozygosity islands suggested a strong selection for production and behavioural traits, whilst the identified cuticula protein-coding genes (CPR3 and CPR4) were associated with their breed-specific stripe pattern. Local adaption of the two subspecies could be confirmed by the identification of two genes involved in the response to ultraviolet (UV) light. We demonstrated that colony genotypes derived from pooled honey bee workers are reliable to unravel the population dynamics in A. mellifera and provide fundamental information to conserve native honey bees.