Convergent Evolution of Slick Coat in Cattle through Truncation Mutations in the Prolactin Receptor

The prolactin receptor: A cross-species comparison of gene structure, transcriptional regulation, tissue-specificity, and genetic variation

The conserved and multifaceted functions of prolactin (PRL) are coordinated through varied distribution and expression of its cell-surface receptor (PRLR) across a range of tissues and physiological states. The resultant heterogeneous expression of PRLR mRNA and protein across different organs and cell types supports a wide range of PRL-regulated processes including reproduction, lactation, development, and homeostasis. Genetic variation within the PRLR gene also accounts for several phenotypes impacting agricultural production and human pathology. The goal of this review is to highlight the many elements that control differential expression of the PRLR across tissues, and the various phenotypes that exist across species due to variation in the PRLR gene.

CRISPR-mediated editing of β-lactoglobulin (BLG) gene in buffalo

Milk is a good source of nutrition but is also a source of allergenic proteins such as α-lactalbumin, β-lactoglobulin (BLG), casein, and immunoglobulins. The Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/Cas technology has the potential to edit any gene, including milk allergens. Previously, CRISPR/Cas has been successfully employed in dairy cows and goats, but buffaloes remain unexplored for any milk trait. In this study, we utilized the CRISPR/Cas9 system to edit the major milk allergen BLG gene in buffaloes. First, the editing efficiency of designed sgRNAs was tested in fibroblast cells using the T7E assay and Sanger sequencing. The most effective sgRNA was selected to generate clonal lines of BLG-edited cells. Analysis of 15 single-cell clones, through TA cloning and Sanger sequencing, revealed that 7 clones exhibited bi-allelic (-/-) heterozygous, bi-allelic (-/-) homozygous, and mono-allelic (-/+) disruptions in BLG. Bioinformatics prediction analysis confirmed that non-multiple-of-3 edited nucleotide cell clones have frame shifts and early truncation of BLG protein, while multiple-of-3 edited nucleotides resulted in slightly disoriented protein structures. Somatic cell nuclear transfer (SCNT) method was used to produce blastocyst-stage embryos that have similar developmental rates and quality with wild-type embryos. This study demonstrated the successful bi-allelic editing (-/-) of BLG in buffalo cells through CRISPR/Cas, followed by the production of BLG-edited blastocyst stage embryos using SCNT. With CRISPR and SCNT methods described herein, our long-term goal is to generate gene-edited buffaloes with BLG-free milk.

Horizon scanning of potential environmental applications of terrestrial animals, fish, algae and microorganisms produced by genetic modification, including the use of new genomic techniques

With scientific progress and the development of new genomic techniques (NGTs), the spectrum of organisms modified for various purposes is rapidly expanding and includes a wide range of taxonomic groups. An improved understanding of which newly developed products may be introduced into the market and released into the environment in the near and more distant future is of particular interest for policymakers, regulatory authorities, and risk assessors. To address this information need, we conducted a horizon scanning (HS) of potential environmental applications in four groups of organisms: terrestrial animals (excluding insects and applications with gene drives), fish, algae and microorganisms. We applied a formal scoping review methodology comprising a structured search of the scientific literature followed by eligibility screening, complemented by a survey of grey literature, and regulatory websites and databases. In all four groups of organisms we identified a broad range of potential applications in stages of basic as well as advanced research, and a limited number of applications which are on, or ready to be placed on, the market. Research on GM animals including fish is focused on farmed animals and primarily targets traits which increase performance, influence reproduction, or convey resistance against diseases. GM algae identified in the HS were all unicellular, with more than half of the articles concerning biofuel production. GM algae applications for use in the environment include biocontrol and bioremediation, which are also the main applications identified for GM microorganisms. From a risk assessor's perspective these potential applications entail a multitude of possible pathways to harm. The current limited level of experience and limited amount of available scientific information could constitute a significant challenge in the near future, for which risk assessors and competent authorities urgently need to prepare.

Thermotolerance capabilities, blood metabolomics, and mammary gland hemodynamics and transcriptomic profiles of slick-haired Holstein cattle during mid lactation in Puerto Rico

Holstein cattle carrying a prolactin receptor gene mutation (SLICK) exhibit short and sleek hair coats (short-haired Holstein [SLK]) enhancing thermotolerance and productivity compared with wild type-haired Holstein (WT) under tropical conditions. The objectives were to unravel the physiological and molecular mechanisms that confer an advantage to this slick genotype in Puerto Rico and determine potential correlations between metabolites and physiological variables. At 160 ± 3 DIM we compared vaginal temperatures (VT) and voluntary solar radiation exposure (VSRE) during 48 h between 9 SLK and 9 WT Holsteins, whereas a subsample of 7 SLK and 7 WT were used to assess udder skin temperature, mammary gland hemodynamics and transcriptomics, and blood plasma untargeted metabolomics at a single time point. The SLK cattle showed lower VT throughout the day and greater VSRE at 1000 h and 1100 h compared with their WT counterparts. Total mammary blood flow (MBF) was greater in SLK Holsteins compared with WT. The metabolite 9-nitrooctadecenoic acid was identified as a potential biomarker for MBF; moreover, SLK cattle had greater amounts of this metabolite in their plasma. Prostaglandin D2 synthase (PTGS) was upregulated in the slick mammary gland, while plasma prostaglandin D2 was positively correlated with milk yield and increased in SLK Holsteins compared with WT. Interestingly, the arachidonic acid metabolism pathway was enriched in the mammary gland transcriptome and perturbed in the blood metabolome in the SLK Holsteins. In conclusion, SLK Holsteins exhibited lower body temperatures, greater VSRE, enhanced blood supply to the mammary gland, and alterations in genes and metabolites involved in arachidonic acid metabolism at the mammary gland and blood plasma. The usage of the SLK Holstein cattle genetics in dairy operations could be a feasible alternative to mitigate the adverse consequences of heat stress.

Candidate signatures of positive selection for environmental adaptation in indigenous African cattle: A review

Environmental adaptation traits of indigenous African cattle are increasingly being investigated to respond to the need for sustainable livestock production in the context of unpredictable climatic changes. Several studies have highlighted genomic regions under positive selection probably associated with adaptation to environmental challenges (e.g. heat stress, trypanosomiasis, tick and tick-borne diseases). However, little attention has focused on pinpointing the candidate causative variant(s) controlling the traits. This review compiled information from 22 studies on signatures of positive selection in indigenous African cattle breeds to identify regions under positive selection. We highlight some key candidate genome regions and genes of relevance to the challenges of living in extreme environments (high temperature, high altitude, high infectious disease prevalence). They include candidate genes involved in biological pathways relating to innate and adaptive immunity (e.g. BoLAs, SPAG11, IL1RL2 and GFI1B), heat stress (e.g. HSPs, SOD1 and PRLH) and hypoxia responses (e.g. BDNF and INPP4A). Notably, the highest numbers of candidate regions are found on BTA3, BTA5 and BTA7. They overlap with genes playing roles in several biological functions and pathways. These include but are not limited to growth and feed intake, cell stability, protein stability and sweat gland development. This review may further guide targeted genome studies aiming to assess the importance of candidate causative mutations, within regulatory and protein-coding genome regions, to further understand the biological mechanisms underlying African cattle's unique adaption.

Global dispersal and adaptive evolution of domestic cattle: a genomic perspective

Domestic cattle have spread across the globe and inhabit variable and unpredictable environments. They have been exposed to a plethora of selective pressures and have adapted to a variety of local ecological and management conditions, including UV exposure, diseases, and stall-feeding systems. These selective pressures have resulted in unique and important phenotypic and genetic differences among modern cattle breeds/populations. Ongoing efforts to sequence the genomes of local and commercial cattle breeds/populations, along with the growing availability of ancient bovid DNA data, have significantly advanced our understanding of the genomic architecture, recent evolution of complex traits, common diseases, and local adaptation in cattle. Here, we review the origin and spread of domestic cattle and illustrate the environmental adaptations of local cattle breeds/populations.

Perspectives of gene editing for cattle farming in tropical and subtropical regions

Cattle productivity in tropical and subtropical regions can be severely affected by the environment. Reproductive performance, milk and meat production are compromised by the heat stress imposed by the elevated temperature and humidity. The resulting low productivity contributes to reduce the farmer's income and to increase the methane emissions per unit of animal protein produced and the pressure on land usage. The introduction of highly productive European cattle breeds as well as crossbreeding with local breeds have been adopted as strategies to increase productivity but the positive effects have been limited by the low adaptation of European animals to hot climates and by the reduction of the heterosis effect in the following generations. Gene editing tools allow precise modifications in the animal genome and can be an ally to the cattle industry in tropical and subtropical regions. Alleles associated with production or heat tolerance can be shifted between breeds without the need of crossbreeding. Alongside assisted reproductive biotechnologies and genome selection, gene editing can accelerate the genetic gain of indigenous breeds such as zebu cattle. This review focuses on some of the potential applications of gene editing for cattle farming in tropical and subtropical regions, bringing aspects related to heat stress, milk yield, bull reproduction and methane emissions.

Physiological responses of Holstein calves and heifers carrying the SLICK1 allele to heat stress in California and Florida dairy farms

Inheritance of the SLICK1 allele of the prolactin receptor gene improves thermotolerance of lactating Holstein cows under humid heat stress conditions. The aim of this study was to investigate whether pre- and postweaning Holstein heifers carrying the SLICK1 allele would show physiological responses indicative of higher tolerance to heat stress in high- and low-humidity climates. A total of 101 heifer calves of two age groups heterozygous for the SLICK1 allele and 103 wild-type half-siblings were evaluated during July 2020 in 3 dairy farms in central California and 2 in south Florida. Dry bulb temperature and relative humidity data were recorded during evaluation and used to calculate the temperature-humidity index (THI). Physiological measurements were obtained between 1600 and 1900 h in California, and 1200 and 1400 h in Florida and included rectal temperature, respiration rate, skin temperature, and sweating rate. Data were analyzed via Generalized Linear Mixed Models including the main effects of genotype, state, group, sire, farm within state, and interactions, with THI included as a covariate. The correlations between THI and dependent variables were analyzed via linear regression. The average 24-h THI was higher in Florida compared with California (90 vs. 72, respectively); the main driver of the higher THI in Florida was the high relative humidity (average 85.6% in Florida vs. 36.7% in California). In Florida, the rectal temperature of slick calves was 0.4°C lower than non-slick calves (39.5 ± 0.1 vs 39.9 ± 0.1°C); no differences were detected between slick and non-slick calves in California. Regardless of genotype, heifer calves in Florida had higher respiration rate, higher rectal and skin temperatures, and lower sweating rate than in California. This study is the first to evaluate physiological responses of calves carrying the SLICK1 allele under heat stress conditions in different climates. Our findings demonstrate that the presence of this allele is associated with lower rectal temperatures in pre- and post-weaning Holstein females. According to the physiological parameters evaluated, calves raised in Florida appeared to be under more severe heat stress; in those conditions, the SLICK1 allele was advantageous to confer thermotolerance as evidenced by lower rectal temperature in slick animals.

Effects of the bovine SLICK1 mutation in PRLR on sweat gland area, FOXA1 abundance, and global gene expression in skin

The SLICK1 mutation in the prolactin receptor (PRLR) results in a short-hair coat and increased ability to regulate body temperature during heat stress. It is unclear whether the mutation affects capacity for sweating. The objective of this observational study was to evaluate whether the SLICK1 mutation in PRLR alters characteristics of skin related to sweat gland abundance or function. Skin biopsies from 31 Holstein heifers, including 14 wild-type (SL^-/-) and 17 heterozygous slick (SL^+/-), were subjected to histological analysis to determine the percent of the surface area of skin sections that are occupied by sweat glands. We detected no effect of genotype on this variable. Immunohistochemical analysis of the forkhead transcription factor A1 (FOXA1), a protein essential for sweating in mice, from 6 SL^-/- and 6 SL^+/- heifers indicated twice as much FOXA1 in sweat glandular epithelia of SL^+/- heifers as in SL^-/- heifers. Results from RNA sequencing of skin biopsies from 5 SL^-/- and 7 SL^+/- heifers revealed few genes that were differentially expressed and none that have been associated with sweat gland development or function. In conclusion, results do not support the idea that the SLICK1 mutation changes the abundance of sweat glands in skin, but do show that functional properties of sweat glands, as indicated by increased abundance of immunoreactive FOXA1, are modified by inheritance of the mutation in PRLR.

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.

Cytoplasmic Injection of Zygotes to Genome Edit Naturally Occurring Sequence Variants Into Bovine Embryos

Genome editing provides opportunities to improve current cattle breeding strategies through targeted introduction of natural sequence variants, accelerating genetic gain. This can be achieved by harnessing homology-directed repair mechanisms following editor-induced cleavage of the genome in the presence of a repair template. Introducing the genome editors into zygotes and editing in embryos has the advantage of uncompromised development into live animals and alignment with contemporary embryo-based improvement practices. In our study, we investigated the potential to introduce sequence variants, known from the pre-melanosomal protein 17 (PMEL) and prolactin receptor (PRLR) genes, and produce non-mosaic, edited embryos, completely converted into the precision genotype. Injection of gRNA/Cas9 editors into bovine zygotes to introduce a 3 bp deletion variant into the PMEL gene produced up to 11% fully converted embryos. The conversion rate was increased to up to 48% with the use of TALEN but only when delivered by plasmid. Testing three gRNA/Cas9 editors in the context of several known PRLR sequence variants, different repair template designs and delivery as DNA, RNA or ribonucleoprotein achieved full conversion rates up to 8%. Furthermore, we developed a biopsy-based screening strategy for non-mosaic embryos which has the potential for exclusively producing non-mosaic animals with intended precision edits.

Effects of the SLICK1 mutation in PRLR on regulation of core body temperature and global gene expression in liver in cattle

The SLICK1 mutation in bovine PRLR (c.1382del; rs517047387) is a deletion mutation resulting in a protein with a truncated intracellular domain. Cattle carrying at least one allele have a phenotype characterized by a short hair coat (slick phenotype) and increased resistance to heat stress. Given the pleiotropic nature of prolactin, the mutation may affect other physiological characteristics. The liver is one organ that could potentially be affected because of the expression of PRLR. The mutation is a dominant allele, and heterozygous animals have a similar hair coat to that of animals homozygous for the mutation. Present objectives were to determine whether inheritance of the SLICK1 mutation affects liver gene expression and if animals homozygous for the SLICK1 allele differ from heterozygotes in liver gene expression and regulation of body temperature during heat stress. In one experiment, rectal and ruminal temperatures were less for Holstein heifers that were heterozygous for the SLICK1 allele compared with wildtype heifers. There were 71 differentially expressed genes in liver, with 13 upregulated and 58 downregulated in SLICK1 heterozygotes. Among the ontologies characteristic of differentially expressed genes were those related to immune function and fatty acid and amino acid metabolism. In a prospective cohort study conducted with adult Senepol cattle, body temperature and hepatic gene expression were compared between animals heterozygous or homozygous for the SLICK1 mutation. There were no differences in ruminal temperatures between genotypes, rectal temperature was higher in animals homozygous for the SLICK1 mutation, and there was only one gene in liver that was differentially expressed. It was concluded that inheritance of the SLICK1 allele can exert functional changes beyond those related to hair growth although changes in liver gene expression were not extensive. Results are also consistent with the SLICK1 allele being dominant because there were few differences in phenotype between animals inheriting one or two copies of the allele.

Coat score. A possible explanation for the zebuine selective sweep located on bovine chromosome 5: 47,670,001-48,100,000 bp

Over 65% of the world's cattle population resides in warm areas where heat stress conditions limit the breed of European taurine cattle. Composite breeds were developed to retain the main traits of both parental breeds. The skin plays a central role in animal response to heat stress. Research on the genetic architecture of skin traits has identified genes and regions related to warm resistance skin features. The aim of this study was to determine whether the indicine proportion accounted for coat type or whether there were genes of large effect segregating in Brangus. Bulls (n = 108) were genotyped using microarrays and their coat score and hair length were evaluated. Indicine-taurine genome-wide composition was estimated and GWAS was performed. Although significant correlations between indicine proportion and traits were not observed, four windows of SNPs on BTA4 and BTA5 explained more than 2% of the trait variance. The GWAS for coat score in summer showed the main peak on BTA5:46,941,446-48,030,219 bp, accounting for 4.65% of the variance. Our results suggest that the variation in coat score and undercoat hair length in Argentinian Brangus bulls is associated with the presence of some particular gene variants, rather than with the whole indicine genetic content.

Towards progressive regulatory approaches for agricultural applications of animal biotechnology

Traditional breeding techniques, applied incrementally over thousands of years, have yielded huge benefits in the characteristics of agricultural animals. This is a result of significant, measurable changes to the genomes of those animal species and breeds. Genome editing techniques may now be applied to achieve targeted DNA sequence alterations, with the potential to affect traits of interest to production of agricultural animals in just one generation. New opportunities arise to improve characteristics difficult to achieve or not amenable to traditional breeding, including disease resistance, and traits that can improve animal welfare, reduce environmental impact, or mitigate impacts of climate change. Countries and supranational institutions are in the process of defining regulatory approaches for genome edited animals and can benefit from sharing approaches and experiences to institute progressive policies in which regulatory oversight is scaled to the particular level of risk involved. To facilitate information sharing and discussion on animal biotechnology, an international community of researchers, developers, breeders, regulators, and communicators recently held a series of seven virtual workshop sessions on applications of biotechnology for animal agriculture, food and environmental safety assessment, regulatory approaches, and market and consumer acceptance. In this report, we summarize the topics presented in the workshop sessions, as well as discussions coming out of the breakout sessions. This is framed within the context of past and recent scientific and regulatory developments. This is a pivotal moment for determination of regulatory approaches and establishment of trust across the innovation through-chain, from researchers, developers, regulators, breeders, farmers through to consumers.

Seasonal variation of the estrous cycle length, corpus luteum area, and size of the pre-ovulatory follicle in Criollo Limonero heifers

In order to determine whether seasonal variations may influence the estrous cycle length (ECL), corpus luteum size (CLS), maximum area of CL (MACL), day of cycle with maximum area of CL (DCMACL), and pre-ovulatory follicles size (PFS), ten Criollo Limonero heifers were subjected to daily ultrasound ovary scanning throughout their estrous cycles during three seasons: hot-dry (HD), hot-humid (HH), and wind-rain (WR). The effect of season on ECL, MACL, DCMACL, and PFS was analyzed with an ANOVA (PROC GLM, SAS), whereas, for the effect of season on CLS, an ANOVA with repeated measures (PROC MIXED, SAS) was used. Results showed no effect (P > 0.05) of season on ECL, MACL, and DCMACL. However, size of PFS was larger (P < 0.02) during the WR season and the CLS tended (P < 0.09) to be lower during the HH. In conclusion, the relative stability of ECL, MACL, DCMACL, PFS, and CLS measures suggests no major seasonal variations which could imply adaptation capability of Criollo Limonero cattle to the tropical environment.

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.

Inheritance of the SLICK1 allele of PRLR in cattle

The slick-hair phenotype in cattle is due to one of a series of mutations in the prolactin receptor (PRLR) that cause truncation of the C-terminal region of the protein involved in JAK2/STAT5 activation during prolactin signaling. Here we evaluated whether the inheritance of the SLICK1 allele, the first slick mutation discovered, is inherited in a fashion consistent with Hardy-Weinberg equilibrium. It was hypothesized that any deleterious effect of inheriting the allele on embryonic or fetal function would result in reduced frequency of the allele in offspring. A total of 525 Holstein and Senepol cattle produced from matings involving one or both parents with the SLICK1 allele were genotyped. The observed frequency of the SLICK1 allele (0.247) was not significantly different than the expected frequency of 0.269. These results support the idea that inheritance of the SLICK1 allele does not act in the embryo or fetus to modify its competence to complete development to term.

The Quest for Genes Involved in Adaptation to Climate Change in Ruminant Livestock

Livestock radiated out from domestication centres to most regions of the world, gradually adapting to diverse environments, from very hot to sub-zero temperatures and from wet and humid conditions to deserts. The climate is changing; generally global temperature is increasing, although there are also more extreme cold periods, storms, and higher solar radiation. These changes impact livestock welfare and productivity. This review describes advances in the methodology for studying livestock genomes and the impact of the environment on animal production, giving examples of discoveries made. Sequencing livestock genomes has facilitated genome-wide association studies to localize genes controlling many traits, and population genetics has identified genomic regions under selection or introgressed from one breed into another to improve production or facilitate adaptation. Landscape genomics, which combines global positioning and genomics, has identified genomic features that enable animals to adapt to local environments. Combining the advances in genomics and methods for predicting changes in climate is generating an explosion of data which calls for innovations in the way big data sets are treated. Artificial intelligence and machine learning are now being used to study the interactions between the genome and the environment to identify historic effects on the genome and to model future scenarios.

Application of Gene Editing for Climate Change in Agriculture

Climate change imposes a severe threat to agricultural systems, food security, and human nutrition. Meanwhile, efforts in crop and livestock gene editing have been undertaken to improve performance across a range of traits. Many of the targeted phenotypes include attributes that could be beneficial for climate change adaptation. Here, we present examples of emerging gene editing applications and research initiatives that are aimed at the improvement of crops and livestock in response to climate change, and discuss technical limitations and opportunities therein. While only few applications of gene editing have been translated to agricultural production thus far, numerous studies in research settings have demonstrated the potential for potent applications to address climate change in the near future.

Evidence for and localization of proposed causative variants in cattle and pig genomes

BACKGROUND

This paper reviews the localization of published potential causative variants in contemporary pig and cattle reference genomes, and the evidence for their causality. In spite of the difficulties inherent to the identification of causative variants from genetic mapping and genome-wide association studies, researchers in animal genetics have proposed putative causative variants for several traits relevant to livestock breeding.

RESULTS

For this review, we read the literature that supports potential causative variants in 13 genes (ABCG2, DGAT1, GHR, IGF2, MC4R, MSTN, NR6A1, PHGK1, PRKAG3, PLRL, RYR1, SYNGR2 and VRTN) in cattle and pigs, and localized them in contemporary reference genomes. We review the evidence for their causality, by aiming to separate the evidence for the locus, the proposed causative gene and the proposed causative variant, and report the bioinformatic searches and tactics needed to localize the sequence variants in the cattle or pig genome.

CONCLUSIONS

Taken together, there is usually good evidence for the association at the locus level, some evidence for a specific causative gene at eight of the loci, and some experimental evidence for a specific causative variant at six of the loci. We recommend that researchers who report new potential causative variants use referenced coordinate systems, show local sequence context, and submit variants to repositories.

Methylome Patterns of Cattle Adaptation to Heat Stress

Heat stress has a detrimental impact on cattle health, welfare and productivity by affecting gene expression, metabolism and immune response, but little is known on the epigenetic mechanisms mediating the effect of temperature at the cellular and organism level. In this study, we investigated genome-wide DNA methylation in blood samples collected from 5 bulls of the heat stress resilient Nellore breed and 5 bulls of the Angus that are more heat stress susceptible, exposed to the sun and high temperature-high humidity during the summer season of the Brazilian South-East region. The methylomes were analyzed during and after the exposure by Reduced Representation Bisulfite Sequencing, which provided genome-wide single-base resolution methylation profiles. Significant methylation changes between stressful and recovery periods were observed in 819 genes. Among these, 351 were only seen in Angus, 366 were specific to Nellore, and 102 showed significant changes in methylation patterns in both breeds. KEGG and Gene Ontology (GO) enrichment analyses showed that responses were breed-specific. Interestingly, in Nellore significant genes and pathways were mainly involved in stress responses and cellular defense and were under methylated during heat stress, whereas in Angus the response was less focused. These preliminary results suggest that heat challenge induces changes in methylation patterns in specific loci, which should be further scrutinized to assess their role in heat tolerance.

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

Blanco Orejinegro (BON) cattle have 500 years of adaptation to the Colombian tropic, but little is known about their genetic history. Our aim was to estimate levels of linkage disequilibrium (LD), effective population size (Ne), genomic inbreeding for runs of homozygosity (F_ROH ), genomic relation matrix (F_GRM ), excess of homozygotes (F_HOM ) and pedigree information (F_PEDCOMP ) and to characterize the runs of homozygosity (ROH), searching for selection signatures. A total of 419 BON animals were genotyped, 70 with a 150K chip and 349 with a 50K chip. Next, an imputation to 50K was performed, and, after editing, databases of 40K were obtained. The PLINK v1.90 and R programs were used to estimate LD, ROH, F_ROH and F_HOM . The SNeP v1.1 program was used to obtain Ne, and PreGSf90 was used to elaborate the scaled G matrix. The MTDFNRM program was used to estimate F_PEDCOMP . The LD mean as r² at 1 Mb was 0.21 (r²  > 0.30 at a distance of 96.72kb), and Ne was 123 ± 1. A total of 7,652 homozygous segments were obtained, with a mean of 18.35 ± 0.55 ROH/animal. Most of the genome was covered by long ROHs (ROH_{>8 Mb}  = 4.86%), indicating significant recent inbreeding. The average inbreeding coefficient for F_PEDCOM , F_GRM , F_HOM and F_ROH was 4.41%, 4.18%, 5.58% and 6.78%, respectively. The highest correlation was observed between F_HOM and F_ROH (0.95). ROH hotspots/islands were defined using the extreme values of a box plot that was generated, and correspond to QTLs related to milk yield (55.11%), external appearance (13.47%), production (13.30%), reproduction (8.15%), health (5.24%) and meat carcass (4.74%).

Genetic Diversity and Signatures of Selection for Thermal Stress in Cattle and Other Two Bos Species Adapted to Divergent Climatic Conditions

Understanding the biological mechanisms of climatic adaptation is of paramount importance for the optimization of breeding programs and conservation of genetic resources. The aim of this study was to investigate genetic diversity and unravel genomic regions potentially under selection for heat and/or cold tolerance in thirty-two worldwide cattle breeds, with a focus on Chinese local cattle breeds adapted to divergent climatic conditions, Datong yak (Bos grunniens; YAK), and Bali (Bos javanicus) based on dense SNP data. In general, moderate genetic diversity levels were observed in most cattle populations. The proportion of polymorphic SNP ranged from 0.197 (YAK) to 0.992 (Mongolian cattle). Observed and expected heterozygosity ranged from 0.023 (YAK) to 0.366 (Sanhe cattle; SH), and from 0.021 (YAK) to 0.358 (SH), respectively. The overall average inbreeding (±SD) was: 0.118 ± 0.028, 0.228 ± 0.059, 0.194 ± 0.041, and 0.021 ± 0.004 based on the observed versus expected number of homozygous genotypes, excess of homozygosity, correlation between uniting gametes, and runs of homozygosity (ROH), respectively. Signatures of selection based on multiple scenarios and methods (F ST, HapFLK, and ROH) revealed important genomic regions and candidate genes. The candidate genes identified are related to various biological processes and pathways such as heat-shock proteins, oxygen transport, anatomical traits, mitochondrial DNA maintenance, metabolic activity, feed intake, carcass conformation, fertility, and reproduction. This highlights the large number of biological processes involved in thermal tolerance and thus, the polygenic nature of climatic resilience. A comprehensive description of genetic diversity measures in Chinese cattle and YAK was carried out and compared to 24 worldwide cattle breeds to avoid potential biases. Numerous genomic regions under positive selection were detected using three signature of selection methods and candidate genes potentially under positive selection were identified. Enriched function analyses pinpointed important biological pathways, molecular function and cellular components, which contribute to a better understanding of the biological mechanisms underlying thermal tolerance in cattle. Based on the large number of genomic regions identified, thermal tolerance has a complex polygenic inheritance nature, which was expected considering the various mechanisms involved in thermal stress response.

The Identification of Runs of Homozygosity Gives a Focus on the Genetic Diversity and Adaptation of the "Charolais de Cuba" Cattle

Simple Summary

The Charolais de Cuba cattle is a tropical adapted breed founded in Cuba around 120 years ago from Charolais French specimens. Nowadays, it is still a closed breed and remains as a small population. In this work, we analyzed the inbreeding and diversity patterns, as well as the population size, of this recent adapted breed via a run of homozygosity (ROH) analysis. We found that the genomic inbreeding levels are higher in the Charolais de Cuba breed compared to French and British Charolais populations. Nevertheless, we detected that the effective population size experienced a very similar decline during the last century in the three Charolais populations studied. Finally, a number of regions with exceptional patterns of long homozygosity were identified in this breed, and these could be related to processes of adaptation to tropical conditions.

Abstract

Inbreeding and effective population size (Ne) are fundamental indicators for the management and conservation of genetic diversity in populations. Genomic inbreeding gives accurate estimates of inbreeding, and the Ne determines the rate of the loss of genetic variation. The objective of this work was to study the distribution of runs of homozygosity (ROHs) in order to estimate genomic inbreeding (F_ROH) and an effective population size using 38,789 Single Nucleotide Polymorphisms (SNPs) from the Illumina Bovine 50K BeadChip in 86 samples from populations of Charolais de Cuba (n = 40) cattle and to compare this information with French (n = 20) and British Charolais (n = 26) populations. In the Cuban, French, and British Charolais populations, the average estimated genomic inbreeding values using the F_ROH statistics were 5.7%, 3.4%, and 4%, respectively. The dispersion measured by variation coefficient was high at 43.9%, 37.0%, and 54.2%, respectively. The effective population size experienced a very similar decline during the last century in Charolais de Cuba (from 139 to 23 individuals), in French Charolais (from 142 to 12), and in British Charolais (from 145 to 14) for the ~20 last generations. However, the high variability found in the ROH indicators and F_ROH reveals an opportunity for maintaining the genetic diversity of this breed with an adequate mating strategy, which can be favored with the use of molecular markers. Moreover, the detected ROH were compared to previous results obtained on the detection of signatures of selection in the same breed. Some of the observed signatures were confirmed by the ROHs, emphasizing the process of adaptation to tropical climate experienced by the Charolais de Cuba population.

Association analysis of loci implied in "buffering" epistasis

The existence of buffering mechanisms is an emerging property of biological networks, and this results in the buildup of robustness through evolution. So far, there are no explicit methods to find loci implied in buffering mechanisms. However, buffering can be seen as interaction with genetic background. Here we develop this idea into a tractable model for quantitative genetics, in which the buffering effect of one locus with many other loci is condensed into a single statistical effect, multiplicative on the total additive genetic effect. This allows easier interpretation of the results and simplifies the problem of detecting epistasis from quadratic to linear in the number of loci. Using this formulation, we construct a linear model for genome-wide association studies that estimates and declares the significance of multiplicative epistatic effects at single loci. The model has the form of a variance components, norm reaction model and likelihood ratio tests are used for significance. This model is a generalization and explanation of previous ones. We test our model using bovine data: Brahman and Tropical Composite animals, phenotyped for body weight at yearling and genotyped at high density. After association analysis, we find a number of loci with buffering action in one, the other, or both breeds; these loci do not have a significant statistical additive effect. Most of these loci have been reported in previous studies, either with an additive effect or as footprints of selection. We identify buffering epistatic SNPs present in or near genes reported in the context of signatures of selection in multi-breed cattle population studies. Prominent among these genes are those associated with fertility (INHBA, TSHR, ESRRG, PRLR, and PPARG), growth (MSTN, GHR), coat characteristics (KIT, MITF, PRLR), and heat resistance (HSPA6 and HSPA1A). In these populations, we found loci that have a nonsignificant statistical additive effect but a significant epistatic effect. We argue that the discovery and study of loci associated with buffering effects allow attacking the difficult problems, among others, of the release of maintenance variance in artificial and natural selection, of quick adaptation to the environment, and of opposite signs of marker effects in different backgrounds. We conclude that our method and our results generate promising new perspectives for research in evolutionary and quantitative genetics based on the study of loci that buffer effect of other loci.

Genome-wide association study identifies variants associated with hair length in Brangus cattle

Thermal stress limits beef cattle production and a shorter hair coat is a key thermoregulative adaptation that allows cattle to lose heat more efficiently. The objective of this study was to identify genetic variants associated with the length of the undercoat and topcoat of cattle utilizing 1456 Brangus heifers genotyped with the Bovine GGP F250 array. Seven SNPs in the PCCA gene were significantly associated with undercoat length. PCCA belongs to the biotin transport and metabolism pathway. Biotin deficiency has been reported to cause hair loss. Four SNPs in an 110 kb including a missense mutation in the PRLR gene were significantly associated with topcoat length. Whereas the association of this polymorphism with hair length is novel, the SLICK mutation in PRLR has previously been demonstrated to significantly impact hair length in cattle. These newly detected genetic variants may contribute to a shorter hair coat and more thermotolerant animals.

Prospects for gene introgression or gene editing as a strategy for reduction of the impact of heat stress on production and reproduction in cattle

In cattle, genetic variation exists in regulation of body temperature and stabilization of cellular function during heat stress. There are opportunities to reduce the impact of heat stress on cattle production by identifying the causative mutations responsible for genetic variation in thermotolerance and transferring specific alleles that confer thermotolerance to breeds not adapted to hot climates. An example of a mutation conferring superior ability to regulate body temperature is the group of frame-sift mutations in the prolactin receptor gene (PRLR) that lead to a truncated receptor and development of cattle with a short, sleek hair coat. Slick mutations in PRLR have been found in several extant breeds derived from criollo cattle. The slick mutation in Senepol cattle has been introgressed into dairy cattle in Puerto Rico, Florida and New Zealand. An example of a mutation that confers cellular protection against elevated body temperature is a deletion mutation in the promoter region of a heat shock protein 70 gene called HSPA1L. Inheritance of the mutation results in amplification of the transcriptional response of HSPA1L to heat shock and increased cell survival. The case of PRLR provides a promising example of the efficacy of the genetic approach outlined in this paper. Identification of other mutations conferring thermotolerance at the whole-animal or cellular level will lead to additional opportunities for using genetic solutions to reduce the impact of heat stress.

Application of genome editing in farm animals: cattle

Milk and meat from cattle and buffaloes contribute 45% of the global animal protein supply, followed by chickens (31%), and pigs (20%). In 2016, the global cattle population of 1.0 billion head produced 6.5 billion tons of cows' milk, and 66 million tons of beef. In the past century, cattle breeding programs have greatly increased the yield per animal with a resultant decrease in the emissions intensity per unit of milk or beef, but this has not been true in all regions. Genome editing research in cattle to date has focused on disease resistance (e.g. tuberculosis), production (e.g. myostatin knockout; production of all-male offspring), elimination of allergens (e.g. beta-lactoglobulin knockout) and welfare (e.g. polled or hornlessness) traits. Modeling has revealed how the use of genome editing to introduce beneficial alleles into cattle breeds could maintain or even accelerate the rate of genetic gain accomplished by conventional breeding programs, and is a superior approach to the lengthy process of introgressing those same alleles from distant breeds. Genome editing could be used to precisely introduce useful alleles (e.g. heat tolerance, disease resistance) and haplotypes into native locally-adapted cattle breeds, thereby helping to improve their productivity. As with earlier genetic engineering approaches, whether breeders will be able to employ genome editing in cattle genetic improvement programs will very much depend upon global decisions around the regulatory framework and governance of genome editing for food animals.

Genetic parameters for hair characteristics and core body temperature in a multibreed Brahman-Angus herd1

Thermal stress in hot humid conditions limits cattle production. The objectives for this study were to estimate genetic parameters for hair characteristics and core body temperature under low and high temperature humidity index (THI) conditions. Hair samples were collected and measured for length and diameter. Core body temperature was measured as vaginal temperature every 15 min over a 5-d period using an iButton temperature measuring device implanted in a blank CIDR in 336 heifers from the University of Florida multibreed herd (ranging from 100% Angus to 100% Brahman). Restricted maximum likelihood procedures were used to estimate heritabilities from multiple bivariate animal models using the WOMBAT program. Estimates of heritability for hair diameter, undercoat length, topcoat length, body temperature under low THI conditions, and body temperature under high THI conditions were 0.50, 0.67, 0.42, 0.32, and 0.26, respectively. The genetic parameters estimated in this study indicate a large, exploitable genetic variance which can be selected upon to improve tolerance in cattle. Breed effects for differing compositions of Brahman and Angus were also estimated. As Brahman breed composition increased by 25% undercoat length, topcoat length, body temperature under low THI conditions, and body temperature under high THI conditions decreased by 1.32 mm, 2.94 mm, 0.11 °C, and 0.14 °C, respectively. Under both low and high THI conditions, cattle with 25% Brahman breed composition or greater maintained a significantly lower body temperature than the 100% Angus breed group. The incorporation of Brahman germplasm is recommended for herds that often experience heat stress conditions in order to increase resilience to heat stress.

Reproductive physiology of the heat-stressed dairy cow: implications for fertility and assisted reproduction

Heat stress causes a large decline in pregnancy success per insemination during warm times of the year. Improvements in fertility are possible by exploiting knowledge about how heat stress affects the reproductive process. The oocyte can be damaged by heat stress at the earliest stages of folliculogenesis and remains sensitive to heat stress in the peri-ovulatory period. Changes in oocyte quality due to heat stress are the result of altered patterns of folliculogenesis and, possibly, direct effects of elevated body temperature on the oocyte. While adverse effects of elevated temperature on the oocyte have been observed in vitro, local cooling of the ovary and protective effects of follicular fluid may limit these actions in vivo. Heat stress can also compromise fertilization rate. The first seven days of embryonic development are very susceptible to disruption by heat stress. During these seven days, the embryo undergoes a rapid change in sensitivity to heat stress from being very sensitive (2- to 4-cell stage) to largely resistant (by the morulae stage). Direct actions of elevated temperature on the embryo are likely to be an important mechanism for reduction in embryonic survival caused by heat stress. An effective way to avoid effects of heat stress on the oocyte, fertilization, and early embryo is to bypass the effects through embryo transfer because embryos are typically transferred into females after acquisition of thermal resistance. There may be some opportunity to mitigate effects of heat stress by feeding antioxidants or regulating the endocrine environment of the cow but neither approach has been reduced to practice. The best long-term solution to the problem of heat stress may be to increase genetic resistance of cows to heat stress. Thermotolerance genes exist within dairy breeds and additional genes can be introgressed from other breeds by traditional means or gene editing.

Genomic clues of the evolutionary history of Bos indicus cattle

Together with their sister subspecies Bos taurus, zebu cattle (Bos indicus) have contributed to important socioeconomic changes that have shaped modern civilizations. Zebu cattle were domesticated in the Indus Valley 8000 years before present (YBP). From the domestication site, they expanded to Africa, East Asia, southwestern Asia and Europe between 4000 and 1300 YBP, intercrossing with B. taurus to form clinal variations of zebu ancestry across the landmass of Afro-Eurasia. In the past 150 years, zebu cattle reached the Americas and Oceania, where they have contributed to the prosperity of emerging economies. The zebu genome is characterized by two mitochondrial haplogroups (I1 and I2), one Y chromosome haplogroup (Y3) and three major autosomal ancestral groups (Indian-Pakistani, African and Chinese). Phenotypically, zebu animals are recognized by their hump, large ears and excess skin. They are rustic, resilient to parasites and capable of bearing the hot and humid climates of the tropics. Many resources are available to study the zebu genome, including commercial arrays of SNP, reference assemblies and publicly available genotypes and whole-genome sequences. Nevertheless, many of these resources were initially developed to support research and subsidize industrial applications in B. taurus, and therefore they can produce bias in data analysis. The combination of genomics with precision agriculture holds great promise for the identification of genetic variants affecting economically important traits such as tick resistance and heat tolerance, which were naturally selected for millennia and played a major role in the evolution of B. indicus cattle.

Impact of slick hair trait on physiological and reproductive performance in beef heifers consuming ergot alkaloids from endophyte-infected tall fescue1

Fescue toxicosis is a multifaceted syndrome common in cattle grazing endophyte-infected tall fescue. The objective of this study was to evaluate the impact of the slick hair trait on physiological and reproductive parameters in heifers experiencing fescue toxicosis. Angus × Senepol heifers (n = 31) were blocked by weight (393.5 ± 17.3 kg) and phenotype relative to hair coat at birth, and randomly fed novel endophyte fescue (EN) or endophyte-infected fescue (EI) haylage in a total mixed ration for 91 d. Weekly measurements were collected to monitor heifer growth and response during ergot alkaloids exposure. Following 28 d of treatment, estrus was synchronized and heifers were inseminated. Ovary mapping and AI pregnancy rate were examined via transrectal ultrasonography. Blood samples were taken for genotyping: slick (S) or wildtype (W). Data were analyzed using repeated measures in PROC MIXED of SAS including fescue treatment (EN vs. EI), genotype (S vs. W), and sample collection time as main effects. Body condition scores were decreased for W heifers compared with S heifers (5.48 vs. 5.66, respectively; P < 0.0001). Surface temperature was greater for EI-W heifers (37.2 °C) compared with other groups (36.4, 36.6, 36.7 °C for EN-S, EN-W, EI-S, respectively; P < 0.05). Serum PRL concentrations were reduced for EI heifers compared with EN heifers (133.5 vs. 163.1 ng/mL, respectively; P < 0.05). The average number of 2 to 4 mm follicles were greater in EI-W heifers (13.8 follicles) compared with other groups (12.2, 10.6, and 11.1 for EN-S, EN-W, and EI-S, respectively; P < 0.0001). However, the average number of preovulatory follicles (≥9 mm) were reduced in EI-W heifers (0.52 follicles) compared with other heifer groups (0.94, 0.88, and 0.85 ± 0.04 for EN-S, EN-W, and EI-S, respectively; P < 0.05). Ovulatory follicle size was smaller in EI-W heifers compared with EN-W heifers (9.14 vs. 11.57 mm, respectively; P = 0.05). Corpus luteum area was reduced in EI-W heifers (235.1 mm2) compared with other heifer groups (297.2, 272.7, and 276.8 mm2 for EN-S, EN-W, and EI-S, respectively; P < 0.05). Concentrations of P4 were greater for EN heifers compared with EI heifers (2.7 vs. 1.8 ng/mL, respectively; P < 0.05). Pregnancy was not established in EI-W heifers (0%) compared with other heifer groups (37.5%, 57.1%, and 62.5% for EN-S, EN-W, and EI-S, respectively; P < 0.05). Overall, the slick hair mutation appears to aid in offsetting the physiological symptoms associated with fescue toxicosis and helps to improve reproductive performance.

Impact of Ergot Alkaloids on Female Reproduction in Domestic Livestock Species

Fescue toxicosis is a multifaceted syndrome that elicits many negative effects on livestock consuming ergot alkaloids produced by endophyte-infected tall fescue. The economic losses associated with fescue toxicosis are primarily due to reproductive failure including altered cyclicity, suppressed hormone secretion, reduced pregnancy rates, agalactia, and reduced offspring birth weights. For decades, a multitude of research has investigated the physiological and cellular mechanisms of these reproductive failures associated with fescue toxicosis. This review will summarize the various effects of ergot alkaloids on female reproduction in grazing livestock species.

Prolactin - a pleiotropic factor in health and disease

The principal role of prolactin in mammals is the regulation of lactation. Prolactin is a hormone that is mainly synthesized and secreted by lactotroph cells in the anterior pituitary gland. Prolactin signalling occurs via a unique transmembrane prolactin receptor (PRL-R). The structure of the PRL-R has now been elucidated and is similar to that of many biologically fundamental receptors of the class 1 haematopoietic cytokine receptor family such as the growth hormone receptor. The PRL-R is expressed in a wide array of tissues, and a growing number of biological processes continue to be attributed to prolactin. In this Review, we focus on the newly discovered roles of prolactin in human health and disease, particularly its involvement in metabolic homeostasis including body weight control, adipose tissue, skin and hair follicles, pancreas, bone, the adrenal response to stress, the control of lactotroph cell homeostasis and maternal behaviour. New data concerning the pathological states of hypoprolactinaemia and hyperprolactinaemia will also be presented and discussed.

Livestock Genomics for Developing Countries - African Examples in Practice

African livestock breeds are numerous and diverse, and typically well adapted to the harsh environment conditions under which they perform. They have been used over centuries to provide livelihoods as well as food and nutritional security. However, African livestock systems are dynamic, with many small- and medium-scale systems transforming, to varying degrees, to become more profitable. In these systems the women and men livestock keepers are often seeking new livestock breeds or genotypes - typically those that increase household income through having enhanced productivity in comparison to traditional breeds while maintaining adaptedness. In recent years genomic approaches have started to be utilized in the identification and development of such breeds, and in this article we describe a number of examples to this end from sub-Saharan Africa. These comprise case studies on: (a) dairy cattle in Kenya and Senegal, as well as sheep in Ethiopia, where genomic approaches aided the identification of the most appropriate breed-type for the local productions systems; (b) a cross-breeding program for dairy cattle in East Africa incorporating genomic selection as well as other applications of genomics; (c) ongoing work toward creating a new cattle breed for East Africa that is both productive and resistant to trypanosomiasis; and (d) the use of African cattle as resource populations to identify genomic variants of economic or ecological significance, including a specific case where the discovery data was from a community based breeding program for small ruminants in Ethiopia. Lessons learnt from the various case studies are highlighted, and the concluding section of the paper gives recommendations for African livestock systems to increasingly capitalize on genomic technologies.