Genetic effects on coat colour in cattle: dilution of eumelanin and phaeomelanin pigments in an F2-Backcross Charolais x Holstein population

Genomic regions associated with coat color in Gir cattle

Indicine cattle breeds are adapted to the tropical climate, and their coat plays an important role in this process. Coat color influences thermoregulation and the adhesion of ectoparasites and may be associated with productive and reproductive traits. Furthermore, coat color is used for breed qualification, with breeders preferring certain colors. The Gir cattle is characterized by a wide variety of coat colors. Therefore, we performed genome-wide association studies to identify candidate genes for coat color in Gir cattle. Different phenotype scenarios were considered in the analyses and regions were identified on eight chromosomes. Some regions and many candidate genes are influencing coat color in the Gir cattle, which was found to be a polygenic trait. The candidate genes identified have been associated with white spotting patterns and base coat color in cattle and other species. In addition, a possible epistatic effect on coat color determination in the Gir cattle was suggested. This is the first published study that identified genomic regions and listed candidate genes associated with coat color in Gir cattle. The findings provided a better understanding of the genetic architecture of the trait in the breed and will allow to guide future fine-mapping studies for the development of genetic markers for selection.

Uncovering the architecture of selection in two Bos taurus cattle breeds

Directional selection alters the genome via hard sweeps, soft sweeps, and polygenic selection. However, mapping polygenic selection is difficult because it does not leave clear signatures on the genome like a selective sweep. In populations with temporally stratified genotypes, the Generation Proxy Selection Mapping (GPSM) method identifies variants associated with generation number (or appropriate proxy) and thus variants undergoing directional allele frequency changes. Here, we use GPSM on two large datasets of beef cattle to detect associations between an animal's generation and 11 million imputed SNPs. Using these datasets with high power and dense mapping resolution, GPSM detected a total of 294 unique loci actively under selection in two cattle breeds. We observed that GPSM has a high power to detect selection in the very recent past (<10 years), even when allele frequency changes are small. Variants identified by GPSM reside in genomic regions associated with known breed-specific selection objectives, such as fertility and maternal ability in Red Angus, and carcass merit and coat color in Simmental. Over 60% of the selected loci reside in or near (<50 kb) annotated genes. Using haplotype-based and composite selective sweep statistics, we identify hundreds of putative selective sweeps that likely occurred earlier in the evolution of these breeds; however, these sweeps have little overlap with recent polygenic selection. This makes GPSM a complementary approach to sweep detection methods when temporal genotype data are available. The selected loci that we identify across methods demonstrate the complex architecture of selection in domesticated cattle.

PMEL p.L18del associates with beef quality of Kumamoto sub-breed of Japanese Brown cattle

Japanese Brown cattle is the second most popular breed among Wagyu breeds and raised mainly in Kumamoto and Kochi Prefectures. Typical coat color of the Kumamoto sub-breed is solid brown, but individuals with diluted coat color are sometimes born. We previously detected four SNPs in PMEL gene and identified p.L18del as the causative polymorphism of this diluted phenotype. The current study examined the association between the SNPs in PMEL gene and carcass traits of the Kumamoto sub-breed. Our association analysis revealed that p.L18del had significant effects on BMS (p = 0.0263), meat brightness (p = 0.0179), meat firmness (p = 0.0102), and meat texture (p = 0.0252) and that del allele of this SNP might be useful to improve these carcass traits.

Genome-Wide Association Analysis Identifies the PMEL Gene Affecting Coat Color and Birth Weight in Simmental × Holstein

Coat color and birth weight, as easily selected traits in cattle, play important roles in cattle breeding. Therefore, we carried out a genome-wide association study on birth weight and coat color to identify loci or potential linkage regions in 233 Simmental × Holstein crossbred beef cattle. The results revealed that nine SNPs were significantly associated with coat color (rs137169378, rs110022687, rs136002689, Hypotrichosis_PMel17, PMEL_1, rs134930689, rs383170073, rs109924971, and rs109146332), and these were in RNF41, ZC3H10, ERBB3, PMEL, and OR10A7 on BTA5. Interestingly, rs137169378, rs110022687, rs136002689, Hypotrichosis_PMel17, and PMEL_1 showed strong linkage disequilibrium (r2 > 0.8) and were significantly associated with coat color. Notably, Hypotrichosis_PMel17 and PMEL_1 were located in the gene PMEL (p = 2.22 × 10-18). Among the five significant SNPs associated with coat color, the birth weight of heterozygous individuals (AB) was greater than that of homozygous individuals (AA). Notably, the birth weight of heterozygous individuals with Hypotrichosis_PMel17 and PMEL_1 genotypes was significantly greater than that of homozygous individuals (0.01 < p < 0.05). Interestingly, the two loci were homozygous in black/white individuals and heterozygous in gray/white individuals, and the birth weight of heterozygous brown/white individuals (43.82 ± 5.25 kg) was greater than that of homozygous individuals (42.58 ± 3.09 kg). The birth weight of calves with the parental color (41.95 ± 3.53 kg) was significantly lower than that of calves with a non-parental color (43.54 ± 4.78 kg) (p < 0.05), and the birth weight of gray/white individuals (49.40 ± 7.11 kg) was the highest. Overall, PMEL appears to be a candidate gene affecting coat color in cattle, and coat color may have a selective effect on birth weight. This study provides a foundation for the breeding of beef cattle through GWAS for coat color and birth weight.

Uncovering genes underlying coat color variation in indigenous cattle breeds through genome-wide positive selection

The identification of candidate genes related to pigmentation and under selective sweep provides insights into the genetic basis of pigmentation and the evolutionary forces that have shaped this variation. The selective sweep events in the genes responsible for normal coat color in Indian cattle groups are still unknown. To find coat color genes displaying signs of selective sweeps in the indigenous cattle, we compiled a list of candidate genes previously investigated for their association with coat color and pigmentation. After that, we performed a genome-wide scan of positive selection signatures using the BovineSNP50K Bead Chip in 187 individuals of seven indigenous breeds. We applied a wide range of methods to find evidence of selection, such as Tajima's D, CLR, iHS, varLD, ROH, and FST. We found a total of sixteen genes under selective sweep, that were involved in coat color and pigmentation physiology. These genes are CRIM1 in Gir, MC1R in Sahiwal, MYO5A, PMEL and POMC in Tharparkar, TYRP1, ERBB2, and ASIP in Red Sindhi, MITF, LOC789175, PAX3 and TYR in Ongole, and IRF2, SDR165 and, KIT in Nelore, ADAMTS19 in Hariana. These genes are related to melanin synthesis, the biology of melanocytes and melanosomes, and the migration and survival of melanocytes during development.

Chromosome-level reference genome for North American bison (Bison bison) and variant database aids in identifying albino mutation

We developed a highly contiguous chromosome-level reference genome for North American bison to provide a platform to evaluate the conservation, ecological, evolutionary, and population genomics of this species. Generated from a F1 hybrid between a North American bison dam and a domestic cattle bull, completeness and contiguity exceed that of other published bison genome assemblies. To demonstrate the utility for genome-wide variant frequency estimation, we compiled a genomic variant database consisting of 3 true albino bison and 44 wild-type pelage color bison. Through the examination of genomic variants fixed in the albino cohort and absent in the controls, we identified a nonsynonymous single nucleotide polymorphism (SNP) mutation on chromosome 29 in exon 3 of the tyrosinase gene (c.1114C>T). A TaqMan SNP Genotyping Assay was developed to genotype this SNP in a total of 283 animals across 29 herds. This assay confirmed the absence of homozygous variants in all animals except 7 true albino bison included in this study. In addition, the only heterozygous animals identified were 2 wild-type pelage color dams of albino offspring. Therefore, we propose that this new high-quality bison genome assembly and incipient variant database provides a highly robust and informative resource for genomics investigations for this iconic North American species.

Analysis of the genetic loci of pigment pattern evolution in vertebrates

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

Genomic Analysis of the Endangered Fonni's Dog Breed: A Comparison of Genomic and Phenotypic Evaluation Scores

The Fonni's dog is an ancient Sardinian breed for livestock and property guarding. In recent years, the number of new registrations to the breeding book has slumped and, thus, this breed risks being lost forever. This work refocuses attention to the Fonni's dog, analysing its genomic makeup and comparing different phenotypical and genetic evaluation scores. Thirty Fonni's dogs were ranked by their general accordance to the breed typicality (typicality score) and to the provisional standard by official judges (judges' score). They were genotyped with a 230K SNP BeadChip and compared with 379 dogs of 24 breeds. Genomically, the Fonni's dogs placed themselves near shepherd dogs and showed a unique genetic signature, which was used to create the genomic score. This score better correlated with typicality (ρ = 0.69, p < 0.0001) than the judges' score (ρ = 0.63, p = 0.0004), which showed little variability among the included dogs. Hair texture or colour were significantly associated in the three scores. The Fonni's dog is confirmed as a well-distinguished breed, despite being selected mainly for its work abilities. The evaluation criteria used during dog expositions can be improved to increase their variability and include elements typical of the breed. The recovery of the Fonni's dog would be possible only with a shared vision between the Italian kennel club and breeders, and the support of regional programs.

Identification of differential methylome signatures of white pigmented skin patches in Nili Ravi buffalo of India

The DNA methylation events mark a major epigenetic change in the genome, reflecting non-genetic disease developments and varied phenotypes. The water buffalo is a dairy production animal with wide agro-climatic distribution in India. Breed-wise the coat color of water buffalo varies from ash-gray to jet black. A typical pigmentation pattern is found in one of the breeds of North India, Nili Ravi, with variedly distributed white patches. The DNA methylation pattern could potentially reveal the epigenetic factors responsible for the pigmentation patterns. To address this question, the DNA isolated from the skin tissues of Nili Ravi with varied white pigmentation and black Murrah buffaloes was subjected to reduced representation bisulfite sequencing. DNA methylation analysis revealed, 68.44%, 63.39%, and 47.94% of the promoter regions were hypermethylated in Nili Ravi over-white versus Murrah, Nili Ravi under-white versus Murrah, and Nili Ravi under-white versus Nili Ravi over-white, respectively. Major genes identified to be differentially methylated among over-white and under-white skin tissues in Nili Ravi included TBX2, SNAI2, HERC2, and CITED1. Overall the results have indicated differential methylation patterns to be potentially involved in hyper or hypopigmentation in Nili Ravi and Murrah buffaloes.

PMEL p.Leu18del dilutes coat color of Kumamoto sub-breed of Japanese Brown cattle

Background

Coat color is important for registration and maintenance of livestock. Standard coat color of Kumamoto sub-breed of Japanese Brown cattle is solid brown, but individuals with diluted coat color have been observed recently. In this study, we attempted to identify polymorphism(s) responsible for coat color dilution by whole genome analysis.

Results

One of the diluted cattle possessed 7302 exonic polymorphisms which could affect genes’ function. Among them, 14 polymorphisms in 10 coat color-related genes were assumed to be specific for the diluted cattle. Subsequent genotyping with three diluted cattle and 74 standard cattle elucidated that PMEL p.Leu18del was the causative polymorphism for coat color dilution in this sub-breed. Individuals with del/del type of this polymorphism showed diluted coat color, but coat color of heterozygotes were intermediate with various dilution rates.

Conclusions

Coat color dilution of Kumamoto sub-breed was caused by PMEL p.Leu18del. The causative del allele has been detected in several genetically distant cattle breeds, suggesting that PMEL p.Leu18del can be used as a DNA marker to control cattle coat color.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08916-8.

Genetic architecture of a composite beef cattle population

Composite breeds are widely used in the beef industry. Composites allow producers to combine desirable traits from the progenitor breeds and simplify herd management, without repeated crossbreeding and maintenance of purebreds. In this study, genomic information was used to evaluate the genetic composition and characteristics of a three-breed beef cattle composite. This composite population referred to as Composite Gene Combination (CGC) consisted of 50% Red Angus, 25% Charolais, 25% Tarentaise. A total of 248 animals were used in this study CGC (n=79), Red Angus (n=61), Charolais (n=79) and Tarentaise (n=29). All animals were genotyped with 777k HD panel. Principal component and ADMIXTURE analyses were carried out to evaluate the genetic structure of CGC animals. The ADMIXTURE revealed the proportion of Tarentaise increased to approximately 57% while Charolais decreased to approximately 5%, and Red Angus decreased to 38% across generations. To evaluate these changes in the genomic composition across different breeds and in CGC across generations runs of homozygosity (ROH) were conducted. This analysis showed Red Angus to have the highest total length of ROH segments per animal with a mean of 349.92 Mb and lowest in CGC with a mean of 141.10 Mb. Furthermore, it showed the formation of new haplotypes in CGC around the sixth generation. Selection signatures were evaluated through Fst and HapFlk analyses. Several selection sweeps in CGC were identified especially in chromosomes 5 and 14 which have previously been reported to be associated with coat color and growth traits. The study supports our previous findings that progenitor combinations are not stable over generations and that either direct or natural selection plays a role in modifying the progenitor proportions. Furthermore, the results showed that Tarentaise contributed useful attributes to the composite in a cool semi-arid environment and suggests a re-exploration of this breed's role may be warranted.

Transcriptome analysis of skin color variation during and after overwintering of Malaysian red tilapia

The commercial value of red tilapia is hampered by variations in skin color during overwintering. In this study, three types of skin of red tilapia, including the skin remained pink color during and after overwintering (P), the skin changed from pink color to black color during overwintering and remained black color after overwintering (P-B), and the skin changed from pink color to black color during overwintering but recovered to pink color when the temperature rose after overwintering (P-B-P), were used to analyze their molecular mechanisms of color variation. The transcriptome results revealed that the P, P-B, and P-B-P libraries had 43, 42, and 43 million clean reads, respectively. The top 10 abundance mRNAs and specific mRNAs (specificity measure SPM > 0.9) were screened. After comparing intergroup gene expression levels, there were 2528, 1924, and 1939 differentially expressed genes (DEGs) between P-B-P and P-B, P-B-P and P, and P-B and P, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of color-related mRNAs showed that a number of DEGs, including tyrp1, tyr, pmel, mitf, mc1r, asip, tat, hpdb, and foxd3, might play a potential role in pigmentation. Additionally, the co-expression patterns of genes were detected within the pigment-related pathways by the PPI network from P-B vs. P group. Furthermore, DEGs from the apoptosis and autophagy pathways, such as baxα, beclin1, and atg7, might be involved in the fading of red tilapia melanocytes. The findings will aid in understanding the molecular mechanism underlying skin color variation in red tilapia during and after overwintering as well as lay a foundation for future research aimed at improving red tilapia skin color characteristics.

Evolutionary conservation genomics reveals recent speciation and local adaptation in threatened takins

Incorrect species delimitation will lead to inappropriate conservation decisions, especially for threatened species. The takin (Budorcas taxicolor) is a large artiodactyl endemic to the Himalayan–Hengduan–Qinling Mountains and is well known for its threatened status and peculiar appearance. However, the speciation, intraspecies taxonomy, evolutionary history, and adaptive evolution of this species still remain unclear, which greatly hampers its scientific conservation. Here, we de novo assembled a high-quality chromosome-level genome of takin and resequenced the genomes of 75 wild takins. Phylogenomics revealed that takin was positioned at the root of Caprinae. Population genomics based on the autosome, X chromosome, and Y chromosome SNPs and mitochondrial genomes consistently revealed the existence of two phylogenetic species and recent speciation in takins: the Himalayan takin (B. taxicolor) and the Chinese takin (B. tibetana), with the support of morphological evidence. Two genetically divergent subspecies were identified in both takin species, rejecting three previously proposed taxonomical viewpoints. Furthermore, their distribution boundaries were determined, suggesting that large rivers play important roles in shaping the genetic partition. Compared with the other subspecies, the Qinling subspecies presented the lowest genomic diversity, higher linkage disequilibrium, inbreeding, and genetic load, thus is in urgent need of genetic management and protection. Moreover, coat color gene (PMEL) variation may be responsible for the adaptive coat color difference between the two species following Gloger’s rule. Our findings provide novel insights into the recent speciation, local adaptation, scientific conservation of takins, and biogeography of the Himalaya–Hengduan biodiversity hotspot.

Cross-population selection signatures in Canchim composite beef cattle

Analyses of livestock genomes have been used to detect selection signatures, which are genomic regions associated with traits under selection leading to a change in allele frequency. The objective of the present study was to characterize selection signatures in Canchim composite beef cattle using cross-population analyses with the founder Nelore and Charolais breeds. High-density single nucleotide polymorphism genotypes were available on 395 Canchim representing the target population, along with genotypes from 809 Nelore and 897 Charolais animals representing the reference populations. Most of the selection signatures were co-located with genes whose functions agree with the expectations of the breeding programs; these genes have previously been reported to associate with meat quality, as well as reproductive traits. Identified genes were related to immunity, adaptation, morphology, as well as behavior, could give new perspectives for understanding the genetic architecture of Canchim. Some selection signatures identified genes that were recently introduced in Canchim, such as the loci related to the polled trait.

Board Invited Review: Crossbreeding beef × dairy cattle for the modern beef production system

Current trends in the United States dairy industry suggest that crossbred beef × dairy calves are replacing a proportion of the calf-fed Holstein steers slaughtered for beef each year. Economic pressures value preweaned beef × dairy calves at a premium over preweaned dairy bull calves; however, there is little modern data to support that intensively fed crossbred calves maintain their premium value over dairy steers across the supply chain. Data from international production systems and from historic research suggests that beef × dairy cattle had greater average daily gains and converted feed to gain more efficiently than dairy steers. Regarding carcass characteristics, across the literature crossbreds consistently yielded heavier carcasses that had lower proportions of trim than dairy steers. Fewer comparisons of beef × dairy and dairy steers exist in the literature for other economically relevant carcass characteristics such as ribeye area, backfat, marbling, tenderness, and eating quality. Existing published data are inconsistent among studies, highlighting the necessity for more research tailored to the United States beef production system.

Shell Biosynthesis and Pigmentation as Revealed by the Expression of Tyrosinase and Tyrosinase-like Protein Genes in Pacific Oyster (Crassostrea gigas) with Different Shell Colors

The widely recognized color polymorphisms of molluscan shell have been appreciated for hundreds of years by collectors and scientists, while molecular mechanisms underlying shell pigmentation are still poorly understood. Tyrosinase is a key rate-limiting enzyme for the biosynthesis of melanin. Here, we performed an extensive multi-omics data mining and identified two tyrosinase genes, including tyrosinase and tyrosinase-like protein 2 (Tyr and Typ-2 respectively), in the Pacific oyster Crassostrea gigas, and investigated the expression patterns of tyrosinase during adults and embryogenesis in black and white shell color C. gigas. Tissue expression analysis showed that two tyrosinase genes were both specifically expressed in the mantle, and the expression levels of Tyr and Typ-2 in the edge mantle were significantly higher than that in the central mantle. Besides, Tyr and Typ-2 genes were black shell-specific compared with white shell oysters. In situ hybridization showed that strong signals for Tyr were detected in the inner surface of the outer fold, whereas positive signals for Typ-2 were mainly localized in the outer surface of the outer fold. In the embryos and larvae, the high expression of Tyr mRNA was detected in eyed-larvae, while Typ-2 mRNA was mainly expressed at the trochophore and early D-veliger. Furthermore, the tyrosinase activity in the edge mantle was significantly higher than that in the central mantle. These findings indicated that Tyr gene may be involved in shell pigmentation, and Typ-2 is more likely to play critical roles not only in the formation of shell prismatic layer but also in shell pigmentation. In particular, Typ-2 gene was likely to involve in the initial non-calcified shell of trochophores. The work provides valuable information for the molecular mechanism study of shell formation and pigmentation in C. gigas.

Exploiting within-breed variability in the autochthonous Reggiana breed identified several candidate genes affecting pigmentation-related traits, stature and udder defects in cattle

Autochthonous cattle breeds constitute important reservoirs of genetic diversity. Reggiana is an Italian local cattle breed reared in the north of Italy for the production of a mono‐breed Parmigiano–Reggiano cheese. Reggiana cattle usually have a classical solid red coat colour and pale muzzle. As part of the strategies designed for the sustainable conservation of this genetic resource, we investigated at the genome‐wise level the within‐breed detected variability of three pigmentation‐related traits (intensity of red coat colour, based on three classes – light/diluted, normal and dark; spotted patterns/piebaldism that sometime emerge in the breed; muzzle colour – pink/pale, grey and black), stature, presence/absence and number of supernumerary teats and teat length. A total of 1776 Reggiana cattle (about two‐thirds of the extant breed population) were genotyped with the GeneSeek GGP Bovine 150k SNP array and single‐marker and haplotype‐based GWASs were carried out. The results indicated that two main groups of genetic factors affect the intensity of red coat colour: darkening genes (including EDN3 and a few other genes) and diluting genes (including PMEL and a few other genes). Muzzle colour was mainly determined by MC1R gene markers. Piebaldism was mainly associated with KIT gene markers. Stature was associated with BTA6 markers upstream of the NCAPG–LCORL genes. Teat defects were associated with TBX3/TBX5, MCC and LGR5 genes. Overall, the identified genomic regions not only can be directly used in selection plans in the Reggiana breed, but also contribute to clarifying the genetic mechanisms involved in determining exterior traits in cattle.

Analysis of whole-genome re-sequencing data of ducks reveals a diverse demographic history and extensive gene flow between Southeast/South Asian and Chinese populations

BACKGROUND

The most prolific duck genetic resource in the world is located in Southeast/South Asia but little is known about the domestication and complex histories of these duck populations.

RESULTS

Based on whole-genome resequencing data of 78 ducks (Anas platyrhynchos) and 31 published whole-genome duck sequences, we detected three geographic distinct genetic groups, including local Chinese, wild, and local Southeast/South Asian populations. We inferred the demographic history of these duck populations with different geographical distributions and found that the Chinese and Southeast/South Asian ducks shared similar demographic features. The Chinese domestic ducks experienced the strongest population bottleneck caused by domestication and the last glacial maximum (LGM) period, whereas the Chinese wild ducks experienced a relatively weak bottleneck caused by domestication only. Furthermore, the bottleneck was more severe in the local Southeast/South Asian populations than in the local Chinese populations, which resulted in a smaller effective population size for the former (7100-11,900). We show that extensive gene flow has occurred between the Southeast/South Asian and Chinese populations, and between the Southeast Asian and South Asian populations. Prolonged gene flow was detected between the Guangxi population from China and its neighboring Southeast/South Asian populations. In addition, based on multiple statistical approaches, we identified a genomic region that included three genes (PNPLA8, THAP5, and DNAJB9) on duck chromosome 1 with a high probability of gene flow between the Guangxi and Southeast/South Asian populations. Finally, we detected strong signatures of selection in genes that are involved in signaling pathways of the nervous system development (e.g., ADCYAP1R1 and PDC) and in genes that are associated with morphological traits such as cell growth (e.g., IGF1R).

CONCLUSIONS

Our findings provide valuable information for a better understanding of the domestication and demographic history of the duck, and of the gene flow between local duck populations from Southeast/South Asia and China.

A conserved haplotype in Wagyu cattle contains RAB4A whose encoded protein regulates glucose trafficking in muscle and fat cells

The Wagyu breed of taurine cattle possess favourable genetics for intramuscular fat (IMF) but genomic loci associated with the trait remain under characterised. Here, we report the identification of a previously unidentified genomic region possessing a particular haplotype structure in Wagyu. Through deployment of a genome-wide haplotype detection analysis that captures regions conserved in a target population but not other populations we screened 100 individual Wagyu and contrasted them with 100 individuals from two independent comparison breeds, Charolais and Angus, using high-density SNPs. An extreme level of Wagyu conservation was assigned to a single genomic window (spanning genomic coordinates BTA28:41 088-300 265 bp). In fact, a five-SNP region spanning 27 096 bp is almost perfectly conserved among the 100 Wagyu individuals assayed and partially overlaps RAB4A. Focussing in, two consecutive SNPs (genomic coordinates 236 949 and 239 950) are apparently fixed within the Wagyu (BB and AA respectively), but at mixed frequencies in the other two breeds. These SNPs are located in the two introns straddling exon 7. In a separate analysis using the 1000 Bulls database, we found that, coincident with exon 7 of RAB4A first allele frequencies were highest in the high IMF Japanese Native (Wagyu) breeds (0.78) and lowest in the low IMF indicine breeds (Nelore and Brahman), with intermediate marbling breeds (Angus and Charolais) assigned intermediate rankings (0.42). RAB4A is known to encode a protein that regulates intracellular trafficking of the insulin-regulated glucose transporter GLUT4. RAB4A can be considered an attractive new positional candidate for IMF development.

A novel SNP of MYO1A gene associated with heat-tolerance in Chinese cattle

With the advent of global climate change, heat-tolerance is becoming more and more important to the sustainability of animal husbandry production systems. Previous studies have shown that MYO1A gene associated with pigmentation may be closely related to heat-tolerance in cattle. In this study, a novel missense mutation (NC_037332.1 g.56390345 A > G) was first detected in MYO1A in 891 individuals of 35 cattle breeds, which transformed the amino acid isoleucine into valine. The purpose of this study was to determine the allele frequencies distribution of this locus in Chinese indigenous cattle and to analyze the relationship between this locus and heat-tolerance. Further analysis showed that frequency of wild allele A decreased gradually from northern cattle to southern cattle, whereas frequency of mutant type allele G showed the opposite pattern, which was consistent with the distribution of various climatic conditions of China. Additionally, association analysis was carried out between genotypes and four climatic conditions (annual mean temperature (T), relative humidity (H), temperature-humidity index (THI) and average annual sunshine hours (100-cloudiness) (SR)). Analysis results showed that genotypes were significantly correlated with climatic conditions. Therefore, our results suggest that the novel SNP (rs209559414) is related to heat-tolerance trait of Chinese indigenous cattle.

Analyses of key genes involved in Arctic adaptation in polar bears suggest selection on both standing variation and de novo mutations played an important role

BACKGROUND

Polar bears are uniquely adapted to an Arctic existence. Since their relatively recent divergence from their closest living relative, brown bears, less than 500,000 years ago, the species has evolved an array of novel traits suited to its Arctic lifestyle. Previous studies sought to uncover the genomic underpinnings of these unique characteristics, and disclosed the genes showing the strongest signal of positive selection in the polar bear lineage. Here, we survey a comprehensive dataset of 109 polar bear and 33 brown bear genomes to investigate the genomic variants within these top genes present in each species. Specifically, we investigate whether fixed homozygous variants in polar bears derived from selection on standing variation in the ancestral gene pool or on de novo mutation in the polar bear lineage.

RESULTS

We find that a large number of sites fixed in polar bears are biallelic in brown bears, suggesting selection on standing variation. Moreover, we uncover sites in which polar bears are fixed for a derived allele while brown bears are fixed for the ancestral allele, which we suggest may be a signal of de novo mutation in the polar bear lineage.

CONCLUSIONS

Our findings suggest that, among other mechanisms, natural selection acting on changes in genes derived from a combination of variation already in the ancestral gene pool, and from de novo missense mutations in the polar bear lineage, may have enabled the rapid adaptation of polar bears to their new Arctic environment.

Indication of Premelanosome Protein (PMEL) Expression Outside of Pigmented Bovine Skin Suggests Functions Beyond Eumelanogenesis

The premelanosome protein (PMEL) is important for fibril formation within melanosomes during vertebrate melanogenesis. Fibrils form a matrix for pigment deposition within pigmented tissues such as skin and hair. PMEL mutations are known to modulate eumelanic pigmentation in vertebrates. However, in bovines, PMEL mutations were also found to alter pheomelanic pigmentation resulting in coat color dilution. Furthermore, epistatic effects of a mutated PMEL allele were detected in the phenotypic expression of the bovine hair defect "rat-tail syndrome" (RTS) characterized by charcoal coat color and hair deformation. Reports about PMEL gene expression in non-pigmented tissues raised the hypothesis that there may be unknown functions of the PMEL protein beyond eumelanin deposition to PMEL fibrils. In our study, we analysed the PMEL protein expression in pigmented skin and non-pigmented bovine tissues (non-pigmented skin, thyroid gland, rumen, liver, kidney, and adrenal gland cortex). We found that a processed form of the bovine PMEL protein is expressed in pigmented as well as in non-pigmented tissues, which is in line with gene expression data from targeted RT-PCR and whole transcriptome RNAseq analysis. The PMEL protein is located in membranes and within the cytosol of epithelial cells. Based on our data from bovine tissues, we concluded that at least in cattle PMEL potentially has additional, yet unexplored functions, which might contribute to effects of PMEL mutations on pheomelanin coat color dilution and charcoal coat color in RTS animals. However, indication of PMEL protein in unpigmented cells and tissues will require further confirmation in the future, because there have been no confirmed reports before, which had detected bovine PMEL protein with specific antibodies either in pigmented or unpigmented tissue.

Selection signatures in tropical cattle are enriched for promoter and coding regions and reveal missense mutations in the damage response gene HELB

Background

Distinct domestication events, adaptation to different climatic zones, and divergent selection in productive traits have shaped the genomic differences between taurine and indicine cattle. In this study, we assessed the impact of artificial selection and environmental adaptation by comparing whole-genome sequences from European taurine and Asian indicine breeds and from African cattle. Next, we studied the impact of divergent selection by exploiting predicted and experimental functional annotation of the bovine genome.

Results

We identified selective sweeps in beef cattle taurine and indicine populations, including a 430-kb selective sweep on indicine cattle chromosome 5 that is located between 47,670,001 and 48,100,000 bp and spans five genes, i.e. HELB, IRAK3, ENSBTAG00000026993, GRIP1 and part of HMGA2. Regions under selection in indicine cattle display significant enrichment for promoters and coding genes. At the nucleotide level, sites that show a strong divergence in allele frequency between European taurine and Asian indicine are enriched for the same functional categories. We identified nine single nucleotide polymorphisms (SNPs) in coding regions that are fixed for different alleles between subspecies, eight of which were located within the DNA helicase B (HELB) gene. By mining information from the 1000 Bull Genomes Project, we found that HELB carries mutations that are specific to indicine cattle but also found in taurine cattle, which are known to have been subject to indicine introgression from breeds, such as N’Dama, Anatolian Red, Marchigiana, Chianina, and Piedmontese. Based on in-house genome sequences, we proved that mutations in HELB segregate independently of the copy number variation HMGA2-CNV, which is located in the same region.

Conclusions

Major genomic sequence differences between Bos taurus and Bos indicus are enriched for promoter and coding regions. We identified a 430-kb selective sweep in Asian indicine cattle located on chromosome 5, which carries SNPs that are fixed in indicine populations and located in the coding sequences of the HELB gene. HELB is involved in the response to DNA damage including exposure to ultra-violet light and is associated with reproductive traits and yearling weight in tropical cattle. Thus, HELB likely contributed to the adaptation of tropical cattle to their harsh environment.

Genome Wide Assessment of Genetic Variation and Population Distinctiveness of the Pig Family in South Africa

Genetic diversity is of great importance and a prerequisite for genetic improvement and conservation programs in pigs and other livestock populations. The present study provides a genome wide analysis of the genetic variability and population structure of pig populations from different production systems in South Africa relative to global populations. A total of 234 pigs sampled in South Africa and consisting of village (n = 91), commercial (n = 60), indigenous (n = 40), Asian (n = 5) and wild (n = 38) populations were genotyped using Porcine SNP60K BeadChip. In addition, 389 genotypes representing village and commercial pigs from America, Europe, and Asia were accessed from a previous study and used to compare population clustering and relationships of South African pigs with global populations. Moderate heterozygosity levels, ranging from 0.204 for Warthogs to 0.371 for village pigs sampled from Capricorn municipality in Eastern Cape province of South Africa were observed. Principal Component Analysis of the South African pigs resulted in four distinct clusters of (i) Duroc; (ii) Vietnamese; (iii) Bush pig and Warthog and (iv) a cluster with the rest of the commercial (SA Large White and Landrace), village, Wild Boar and indigenous breeds of Koelbroek and Windsnyer. The clustering demonstrated alignment with genetic similarities, geographic location and production systems. The PCA with the global populations also resulted in four clusters that where populated with (i) all the village populations, wild boars, SA indigenous and the large white and landraces; (ii) Durocs (iii) Chinese and Vietnamese pigs and (iv) Warthog and Bush pig. K = 10 (The number of population units) was the most probable ADMIXTURE based clustering, which grouped animals according to their populations with the exception of the village pigs that showed presence of admixture. AMOVA reported 19.92%-98.62% of the genetic variation to be within populations. Sub structuring was observed between South African commercial populations as well as between Indigenous and commercial breeds. Population pairwise F ST analysis showed genetic differentiation (P ≤ 0.05) between the village, commercial and wild populations. A per marker per population pairwise F ST analysis revealed SNPs associated with QTLs for traits such as meat quality, cytoskeletal and muscle development, glucose metabolism processes and growth factors between both domestic populations as well as between wild and domestic breeds. Overall, the study provided a baseline understanding of porcine diversity and an important foundation for porcine genomics of South African populations.

Genomic measures of inbreeding coefficients and genome-wide scan for runs of homozygosity islands in Iranian river buffalo, Bubalus bubalis

Background

Consecutive homozygous fragments of a genome inherited by offspring from a common ancestor are known as runs of homozygosity (ROH). ROH can be used to calculate genomic inbreeding and to identify genomic regions that are potentially under historical selection pressure. The dataset of our study consisted of 254 Azeri (AZ) and 115 Khuzestani (KHZ) river buffalo genotyped for ~ 65,000 SNPs for the following two purposes: 1) to estimate and compare inbreeding calculated using ROH (F_ROH), excess of homozygosity (F_HOM), correlation between uniting gametes (F_UNI), and diagonal elements of the genomic relationship matrix (F_GRM); 2) to identify frequently occurring ROH (i.e. ROH islands) for our selection signature and gene enrichment studies.

Results

In this study, 9102 ROH were identified, with an average number of 21.2 ± 13.1 and 33.2 ± 15.9 segments per animal in AZ and KHZ breeds, respectively. On average in AZ, 4.35% (108.8 ± 120.3 Mb), and in KHZ, 5.96% (149.1 ± 107.7 Mb) of the genome was autozygous. The estimated inbreeding values based on F_HOM, F_UNI and F_GRM were higher in AZ than they were in KHZ, which was in contrast to the F_ROH estimates. We identified 11 ROH islands (four in AZ and seven in KHZ). In the KHZ breed, the genes located in ROH islands were enriched for multiple Gene Ontology (GO) terms (P ≤ 0.05). The genes located in ROH islands were associated with diverse biological functions and traits such as body size and muscle development (BMP2), immune response (CYP27B1), milk production and components (MARS, ADRA1A, and KCTD16), coat colour and pigmentation (PMEL and MYO1A), reproductive traits (INHBC, INHBE, STAT6 and PCNA), and bone development (SUOX).

Conclusion

The calculated F_ROH was in line with expected higher inbreeding in KHZ than in AZ because of the smaller effective population size of KHZ. Thus, we find that F_ROH can be used as a robust estimate of genomic inbreeding. Further, the majority of ROH peaks were overlapped with or in close proximity to the previously reported genomic regions with signatures of selection. This tells us that it is likely that the genes in the ROH islands have been subject to artificial or natural selection.

Four Novel SNPs of MYO1A Gene Associated with Heat-Tolerance in Chinese Cattle

Based on the previous GWAS research related to bovine heat tolerance trait, this study aimed to explore the effect of myosin-1a (MYO1A) gene on bovine heat tolerance trait, and find the molecular markers related to the heat tolerance of Chinese cattle. In our study, four novel candidate SNPs highly conserved in B. indicus breeds but barely existed in B. taurus were identified in MYO1A gene according to Bovine Genome Variation Database and Selective Signatures (BGVD). PCR and DNA sequencing were used to genotype 1072 individuals including 34 Chinese indigenous cattle breeds as well as Angus and Indian zebu. Two synonymous mutations (rs208210464 and rs110123931), one missense mutation (rs209999142; Phe172Ser), and one intron mutation (rs135771836) were detected. The frequencies of mutant alleles of the four SNPs gradually increased from northern groups to southern groups of Chinese cattle, which was consistent with the distribution of various climatic conditions of China. Additionally, four SNPs were significantly associated with four climatic conditions including annual mean temperature (T), relative humidity (H), temperature-humidity index (THI), and average annual sunshine hours (100-cloudiness) (SR). Among these, rs209999142 and Hap 1/1 had better performance than others. Our results suggested that rs209999142 was associated with heat-tolerance trait and rs208210464, rs110123931, and rs135771836 showed high phenotypic effect on heat-tolerance trait because of the strong linkage with rs209999142. These SNPs could be used as candidates for marker-assisted selection (MAS) in cattle breeding.

Exploring differentially expressed genes associated with coat color in goat skin using RNA-seq

Fur color in domestic goats is an important, genetically determined characteristic that is associated with economic value. This study was designed to perform a comprehensive expression profiling of genes expressed in the skin tissues from Laiwu Black goat and Lubei White goat. Comparisons of black and white goat skin transcriptomes revealed 102 differentially expressed genes (DEGs), of which 38 were upregulated and 64 downregulated in black skin compared with white skin. Among the DEGs, we identified six genes involved in pigmentation, including agouti signaling protein (ASIP), CAMP responsive element binding protein 3-like 1 (CREB3L1), dopachrome tautomerase (DCT), premelanosome protein (PMEL), transient receptor potential cation channel subfamily M member 1 (TRPM1), and tyrosinase-related protein 1 (TYRP1). Notably, there were no significant differences in the expression of melanocortin 1 receptor, microphthalmia-associated transcription factor, tyrosinase, and KIT proto-oncogene receptor tyrosine kinase between the black and white skin samples, whereas ASIP expression was detected only in white skin. PMEL, TRPM1, TYRP1, and DCT showed higher expression in black goat skin, but ASIP and CREB3L1 had higher expression in white goat skin. Quantitative polymerase chain reaction results for PMEL, TRPM1, DCT, TYRP1, and CREB3L1 expression were consistent with those for RNA-seq. These results will expand our understanding of the complex molecular mechanisms of skin physiology and melanogenesis in goats, and provide a foundation for future studies.

Slc7a11 Modulated by POU2F1 is Involved in Pigmentation in Rabbit

Solute carrier family 7 member 11 (Slc7a11) is a cystine/glutamate xCT transporter that controls the production of pheomelanin pigment to change fur and skin color in animals. Previous studies have found that skin expression levels of Slc7a11 varied significantly with fur color in Rex rabbits. However, the molecular regulation mechanism of Slc7a11 in pigmentation is unknown. Here, rabbit melanocytes were first isolated and identified. The distribution and expression pattern of Slc7a11 was confirmed in skin from rabbits with different fur colors. Slc7a11 affected the expression of pigmentation related genes and thus affected melanogenesis. Meanwhile, Slc7a11 decreased melanocyte apoptosis, but inhibition of Slc7a11 enhanced apoptosis. Furthermore, the POU2F1 protein was found to bind to the −713 to −703 bp region of Slc7a11 promoter to inhibit its activity in a dual-luciferase reporter and site-directed mutagenesis assay. This study reveals the function of the Slc7a11 in melanogenesis and provides in-depth analysis of the mechanism of fur pigmentation.

Comparative transcriptome analysis reveals the genetic basis of coat color variation in Pashmina goat

The genetics of coat color variation remains a classic area. Earlier studies have focused on a limited number of genes involved in color determination; however, the complete set of trait determinants are still not well known. In this study, we used high-throughput sequencing technology to identify and characterize intricate interactions between genes that cause complex coat color variation in Changthangi Pashmina goats, producer of finest and costly commercial animal fiber. We systematically identified differentially expressed mRNAs and lncRNAs from black, brown and white Pashmina goat skin samples by using RNA-sequencing technique. A pairwise comparison of black, white and brown skin samples yielded 2479 significantly dysregulated genes (2422 mRNA and 57 lncRNAs). Differentially expressed genes were enriched in melanin biosynthesis, melanocyte differentiation, developmental pigmentation, melanosome transport activities GO terms. Our analysis suggested the potential role of lncRNAs on color coding mRNAs in cis and trans configuration. We have also developed online data repository as a component of the study to provide a central location for data access, visualization and interpretation accessible through http://pcd.skuastk.org/.

Identification and sequence characterization of melanocortin 1 receptor gene (MC1R) in Bos indicus versus (Bos taurus X Bos indicus)

Melanocortin 1 receptor (MC1R) plays a vital role in melanogenesis and determines coat color of mammals. Polymorphic variants in MC1R, causing coat color variation, were described in few mammals; however, such studies were not done in cattle. The objective of the study was to explore the association of MC1R gene polymorphism within Tharparkar (Bos indicus) and Karan Fries (B. indicus X Bos taurus) cattle. Genomic DNA isolated from blood samples of Tharparkar breed by modified Phenol: Chloroform; Isoamyl alcohol method. Using genomic DNA as template for PCR, MC1R gene was amplified and sequenced. The sequences were analyzed and submitted to Genbank with Acc.No MG373615-MG373644. Comparison of sequence alignment with other bovine species using ClustalW revealed 99-96% similarity. MC1R gene phylogenetic analyses were analyzed using MEGA X. The MC1R gene tree, protein domains and genetic variation of cattle were retrieved from Ensemble Asia Cattle Genome Browser. Eight single nucleotide polymorphisms (SNPs) (c.296T > C, c.583T > C, c.663C > T, c.830T > C, c.853G > A, c.880G > A, c.906C > G, c.927C > T) in CDS reveal high genetic variability. Subsequent to amino acid changes p.L99P, p.F195L, p.F277S, p.A285T and p.D293N, p.R302S, respectively found in seven-transmembrane. Mutations appeared in MC1R of B. taurus with white and black coat color as compared to B. indicus with white coat.

Analysis of Anasplatyrhynchos genome resequencing data reveals genetic signatures of artificial selection

Three artificially selected duck populations (AS), higher lean meat ratios (LTPD), higher fat ratios (FTPD) and higher quality meat (CMD), have been developed in China, providing excellent populations for investigation of artificial selection effects. However, the genetic signatures of artificial selection are unclear. In this study, we sequenced the genome sequences of these three artificially selected populations and their ancestral population (mallard, M). We then compared the genome sequences between AS and M and between LTPD and FTPD using integrated strategies such as anchoring scaffolds to pseudo-chromosomes, mutation detection, selective screening, GO analysis, qRT-PCR, and protein multiple sequences alignment to uncover genetic signatures of selection. We anchored duck scaffolds to pseudo-chromosomes and obtained 28 pseudo-chromosomes, accounting for 84% of duck genome in length. Totally 78 and 99 genes were found to be under selection between AS and M and between LTPD and FTPD. Genes under selection between AS and M mainly involved in pigmentation and heart rates, while genes under selection between LTPD and FTPD involved in muscle development and fat deposition. A heart rate regulator (HCN1), the strongest selected gene between AS and M, harbored a GC deletion in AS and displayed higher mRNA expression level in M than in AS. IGF2R, a regulator of skeletal muscle mass, was found to be under selection between FTPD and LTPD. We also found two nonsynonymous substitutions in IGF2R, which might lead to higher IGF2R mRNA expression level in FTPD than LTPD, indicating the two nonsynonymous substitutions might play a key role for the regulation of duck skeletal muscle mass. Taken together, these results of this study provide valuable insight for the genetic basis of duck artificial selection.

Genome-wide scan reveals divergent selection among taurine and zebu cattle populations from different regions

In this study, to identify genomic signatures of divergent selection, we genotyped 10 cattle breeds/populations (n = 275), representing eight Ethiopian cattle populations (n = 229) and two zebu populations (n = 46) adapted to tropical and sub-tropical environments, using the high-density single-nucleotide polymorphisms (SNPs) derived mainly from Bos indicus breeds, and using five reference taurine breeds (n = 212). Population genetic differentiation (F_ST ) values across sliding windows were estimated between zebu and reference combined taurine breeds. The most differentiated regions (F_ST  ≥ 0.53), representing the top 1% smoothed F_ST values, were considered to represent regions under diversifying selection. In total, 285 and 317 genes were identified in the comparisons of Ethiopian cattle with taurine and Asian zebu with taurine respectively. Some of these genes are involved in stress responses/thermo-tolerance and DNA damage repair (HSPA4, HSF1, CMPK1 and EIF2AK4), pigmentation (ERBB3 and MYO1A), reproduction/fertility (UBE2D3, ID3 and PSPC1), immune response (PIK3CD and AKIRIN2) and body stature and size (MBP2, LYN and NPM1). Additionally, the candidate genes were associated with functional terms (e.g. cellular response to stress, DNA repair, inflammatory response) important for physiological adaptation to environmental stresses. The results of our study may shed light on the influence of artificial and natural selection in shaping the genomic diversity of modern cattle breeds and also may serve as a basis for further genetic investigation of traits of tropical adaptation in cattle.

Feather follicles transcriptome profiles in Bashang long-tailed chickens with different plumage colors

Despite the rich variety in plumage color found in nature, genetic studies on how feather follicles affect pigmentation are often limited to animals that have black and white pigment. To test how gene expression influences plumage color, transcriptomes of chicken feather follicles with white, black, hemp, reed catkins, silvery grey, and landscape plumage colors were generated using Illumina sequencing. We generated six RNA-Seq libraries with over 25 million paired-end clean reads per library with percentage of paired-end clean reads ranging from 96.73 to 96.98%. 78% of the reads mapped to the chicken genome, and approximately 70% of the reads were mapped to exons and 6% mapped to introns. Transcriptomes of feather follicles producing hemp and land plumage were similar, but these two showed moderate differences compared with gray and reed colored plumage. The black and white follicle transcriptomes were most divergent from the other colors. We identified several candidate genes, including GPNMB, PMEL, TYRP1, GPR143, OCA2, SOX10, SLC45A2, KRT75, and TYR. All of these genes are known to induce pigment formation in mice. White feathers result from the lack of pigment formation, and our results suggest that the white chickens due to the recessive insertion mutation of TYR. The formation of black area size and color depth may be due to the expression levels of GPNMB, PMEL, TYRP1, GPR143, OCA2, SOX10, SLC45A2, KRT75, and TYR. The GO analysis of the differentially expressed genes (DEGs) revealed that DEGs in our transcriptome analysis were enriched in cytoskeleton and cell structure related pathways. The black plumage transcriptome showed significant differences in melanogenesis, tyrosine metabolism, and riboflavin metabolism compared with transcriptomes of other plumage colors. The transcriptome profiles of the different chicken plumage colors provide a valuable resource to understand how gene expression influences plumage color, and will be an important resource for identifying candidate genes in breeding programs.

Genetic diversity and divergence among Korean cattle breeds assessed using a BovineHD single-nucleotide polymorphism chip

Objective

In Korea, there are three main cattle breeds, which are distinguished by coat color: Brown Hanwoo (BH), Brindle Hanwoo (BRH), and Jeju Black (JB). In this study, we sought to compare the genetic diversity and divergence among there Korean cattle breeds using a BovineHD chip genotyping array.

Methods

Sample data were collected from 168 cattle in three populations of BH (48 cattle), BRH (96 cattle), and JB (24 cattle). The single-nucleotide polymorphism (SNP) genotyping was performed using the Illumina BovineHD SNP 777K Bead chip.

Results

Heterozygosity, used as a measure of within-breed genetic diversity, was higher in BH (0.293) and BRH (0.296) than in JB (0.266). Linkage disequilibrium decay was more rapid in BH and BRH than in JB, reaching an average r² value of 0.2 before 26 kb in BH and BRH, whereas the corresponding value was reached before 32 kb in JB. Intra-population, inter-population, and Fst analyses were used to identify candidate signatures of positive selection in the genome of a domestic Korean cattle population and 48, 11, and 11 loci were detected in the genomic region of the BRH breed, respectively. A Neighbor-Joining phylogenetic tree showed two main groups: a group comprising BH and BRH on one side and a group containing JB on the other. The runs of homozygosity analysis between Korean breeds indicated that the BRH and JB breeds have high inbreeding within breeds compared with BH. An analysis of differentiation based on a high-density SNP chip showed differences between Korean cattle breeds and the closeness of breeds corresponding to the geographic regions where they are evolving.

Conclusion

Our results indicate that although the Korean cattle breeds have common features, they also show reliable breed diversity.

BREEDING AND GENETICS SYMPOSIUM:Breeding heat tolerant dairy cattle: the case for introgression of the "slick" prolactin receptor variant into dairy breeds

Increasing environmental temperatures are a threat to the sustainability of livestock production and, because of the high metabolic demands of lactation, to dairy production in particular. Summer heat waves in temperate climates reduce feed intake, milk production, and cow comfort. In extreme heat events, there is an increase in cow mortality. In tropical climates, dairy cattle are mostly (zebu) type or zebu crossbred with temperate dairy breeds. Crossbreeding is undertaken to combine the heat tolerance and tick resistance of zebu with the productivity of temperate dairy breeds. In the absence of improved heat tolerance, milk production and fertility of temperate cattle is severely impaired. We have recently identified a key role for the prolactin pathway in regulating heat tolerance. A de novo mutation in prolactin that impairs prolactin activity was discovered in hairy and heat intolerant, New Zealand dairy cattle. The phenotypes produced were remarkably similar to those seen in fescue toxicosis, a syndrome seen in grazing cattle in the U.S. where ingestion of ergovaline, a fungal toxin from infected pasture, inhibits prolactin secretion. Recognition of the role of prolactin in hairy cattle led us to identify a deletion in exon 10 of the long-form of the prolactin receptor in Senepol cattle that causes truncation of the protein and determines the slick coat and heat tolerance traits found in this , beef breed. The short form of the prolactin receptor is predicted to be unaffected by the deletion. Knowledge of this dominant mutation has provided the impetus to begin a crossbreeding program to investigate performance and heat tolerance of temperate dairy cattle carrying the slick, prolactin receptor variant. The perceived opportunity is to introgress this variant into temperate dairy cattle to enable performance and welfare improvement in hot climates. Heat tolerance of cattle with slick coats appears to be mostly associated with coat type although sweating ability may also be enhanced. Further investigation is required of performance traits in cows homozygous for the slick variant because the published data are almost exclusively from heterozygous animals. Combination of the slick mutation with other favorable genes for heat tolerance, especially those for coat color, will be particularly enabled by gene editing technologies, offering opportunities for further improvement in bovine thermotolerance.

Selection signatures in Canchim beef cattle

Background

Recent technological advances in genomics have allowed the genotyping of cattle through single nucleotide polymorphism (SNP) panels. High-density SNP panels possess greater genome coverage and are useful for the identification of conserved regions of the genome due to selection, known as selection signatures (SS). The SS are detectable by different methods, such as the extended haplotype homozygosity (EHH); and the integrated haplotype score (iHS), which is derived from the EHH. The aim of this study was to identify SS regions in Canchim cattle (composite breed), genotyped with high-density SNP panel.

Results

A total of 687,655 SNP markers and 396 samples remained for SS analysis after the genotype quality control. The iHS statistic for each marker was transformed into piHS for better interpretation of the results. Chromosomes BTA5 and BTA14 showed piHS > 5, with 39 and nine statistically significant SNPs (P < 0.00001), respectively. For the candidate selection regions, iHS values were computed across the genome and averaged within non-overlapping windows of 500 Kb. We have identified genes that play an important role in metabolism, melanin biosynthesis (pigmentation), and embryonic and bone development.

Conclusions

The observation of SS indicates that the selection processes performed in Canchim, as well as in the founder breeds (i.e. Charolais), are maintaining specific genomic regions, particularly on BTA5 and BTA14. These selection signatures regions could be associated with Canchim characterization.

Electronic supplementary material

The online version of this article (doi:10.1186/s40104-016-0089-5) contains supplementary material, which is available to authorized users.

Recent development of allele frequencies and exclusion probabilities of microsatellites used for parentage control in the German Holstein Friesian cattle population

Background

Methods for parentage control in cattle have changed since their initial implementation in the late 1950’s from blood group typing to more current single nucleotide polymorphism determination. In the early 1990’s, 12 microsatellites were selected by the International Society for Animal Genetics based on their informativeness and robustness in a variety of different cattle breeds. Since then this panel is used as standard in cattle herd book breeding and its application is accompanied by recurrent international comparison tests ensuring permanent validity for the most common commercial dairy and beef cattle breeds for example Holstein Friesian, Simmental, Angus, and Hereford. Although, nearly every parentage can be resolved using these microsatellites, cases with very close relatives became an emerging resolution problem during recent years. This is mainly due to an increase of monomorphism and a trend to the fixation of alleles, although no direct selection against their variability was applied. Thus other effects must be presumed resulting in a loss of polymorphism information content, heterozygosity, and exclusion probabilities.

Results

To determine changes of allele frequencies and exclusion probabilities, we analyzed the development of these parameters for the 12 microsatellites from 2004 to 2014. One hundred sixty eight thousand recorded Holstein Friesian cattle genotypes were evaluated. During this period certain alleles of nine microsatellites increased significantly (t-values >5). When calculating the exclusion probabilities for 11 microsatellites, reduction was determined for the three situations, i.e. one parent is wrongly identified (p = 0.01), both parents are wrongly identified (p = 0.005), and the genotype of one parent is missing (p = 0.048). With the addition of BM1818 to the marker set in 2009, this development was corrected leading to significant increases in exclusion probabilities. Although, the exclusion probabilities for the three family situations using the 12 microsatellites are >99 %, the clarification of 142 relationships in 40,000 situations where one parent is missing will still be impossible.

Twenty-five sires were identified that are responsible for the most significant microsatellite allele increases in the population. The corresponding alleles are mainly associated with milk protein and fat yield, body weight at birth and weaning, as well as somatic cell score, milk fat percentage, and longissimus muscle area.

Conclusions

Our data show that most of the microsatellites used for parentage control in cattle show directional changes in allele frequencies consistent with the history of artificial selection in the German Holstein population.

An interpretive review of selective sweep studies in Bos taurus cattle populations: identification of unique and shared selection signals across breeds

This review compiles the results of 21 genomic studies of European Bos taurus breeds and thus provides a general picture of the selection signatures in taurine cattle identified by genome-wide selection-mapping scans. By performing a comprehensive summary of the results reported in the literature, we compiled a list of 1049 selection sweeps described across 37 cattle breeds (17 beef breeds, 14 dairy breeds, and 6 dual-purpose breeds), and four different beef-vs.-dairy comparisons, which we subsequently grouped into core selective sweep (CSS) regions, defined as consecutive signals within 1 Mb of each other. We defined a total of 409 CSSs across the 29 bovine autosomes, 232 (57%) of which were associated with a single-breed (Single-breed CSSs), 134 CSSs (33%) were associated with a limited number of breeds (Two-to-Four-breed CSSs) and 39 CSSs (9%) were associated with five or more breeds (Multi-breed CSSs). For each CSS, we performed a candidate gene survey that identified 291 genes within the CSS intervals (from the total list of 5183 BioMart-extracted genes) linked to dairy and meat production, stature, and coat color traits. A complementary functional enrichment analysis of the CSS positional candidates highlighted other genes related to pathways underlying behavior, immune response, and reproductive traits. The Single-breed CSSs revealed an over-representation of genes related to dairy and beef production, this was further supported by over-representation of production-related pathway terms in these regions based on a functional enrichment analysis. Overall, this review provides a comparative map of the selection sweeps reported in European cattle breeds and presents for the first time a characterization of the selection sweeps that are found in individual breeds. Based on their uniqueness, these breed-specific signals could be considered as “divergence signals,” which may be useful in characterizing and protecting livestock genetic diversity.

Genomic regions influencing coat color saturation and facial markings in Fleckvieh cattle

Genomic regions associated with coat color and pigmented areas of the head were identified for Fleckvieh (dual-purpose Simmental), a red-spotted and white-headed cattle breed. Coat color was measured with a chromameter, implementing the CIELAB color space and resulting in numerical representation of lightness, color intensity, red/green and blue/yellow color components, rather than subjective classification. Single marker regression analyses with fixed effects of the sex and barn were applied, and significant regions were determined with the local false discovery rate methodology. The PMEL and ERBB3 genes on chromosome 5 were in the most significant region for the color measurements. In addition to the blue/yellow color component and color intensity, the AP3B2 gene on chromosome 21 was identified. Its function was confirmed for similar traits in a range of model species. The KIT gene on chromosome 6 was found to be strongly associated with the inhibition of circum-ocular pigmentation and pigmented spots on the cheek.

Population genomics reveal recent speciation and rapid evolutionary adaptation in polar bears

Polar bears are uniquely adapted to life in the High Arctic and have undergone drastic physiological changes in response to Arctic climates and a hyper-lipid diet of primarily marine mammal prey. We analyzed 89 complete genomes of polar bear and brown bear using population genomic modeling and show that the species diverged only 479-343 thousand years BP. We find that genes on the polar bear lineage have been under stronger positive selection than in brown bears; nine of the top 16 genes under strong positive selection are associated with cardiomyopathy and vascular disease, implying important reorganization of the cardiovascular system. One of the genes showing the strongest evidence of selection, APOB, encodes the primary lipoprotein component of low-density lipoprotein (LDL); functional mutations in APOB may explain how polar bears are able to cope with life-long elevated LDL levels that are associated with high risk of heart disease in humans.

Selection for complex traits leaves little or no classic signatures of selection

Background

Selection signatures aim to identify genomic regions underlying recent adaptations in populations. However, the effects of selection in the genome are difficult to distinguish from random processes, such as genetic drift. Often associations between selection signatures and selected variants for complex traits is assumed even though this is rarely (if ever) tested. In this paper, we use 8 breeds of domestic cattle under strong artificial selection to investigate if selection signatures are co-located in genomic regions which are likely to be under selection.

Results

Our approaches to identify selection signatures (haplotype heterozygosity, integrated haplotype score and F_ST) identified strong and recent selection near many loci with mutations affecting simple traits under strong selection, such as coat colour. However, there was little evidence for a genome-wide association between strong selection signatures and regions affecting complex traits under selection, such as milk yield in dairy cattle. Even identifying selection signatures near some major loci was hindered by factors including allelic heterogeneity, selection for ancestral alleles and interactions with nearby selected loci.

Conclusions

Selection signatures detect loci with large effects under strong selection. However, the methodology is often assumed to also detect loci affecting complex traits where the selection pressure at an individual locus is weak. We present empirical evidence to suggests little discernible ‘selection signature’ for complex traits in the genome of dairy cattle despite very strong and recent artificial selection.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-246) contains supplementary material, which is available to authorized users.

The challenges and importance of structural variation detection in livestock

Recent studies in humans and other model organisms have demonstrated that structural variants (SVs) comprise a substantial proportion of variation among individuals of each species. Many of these variants have been linked to debilitating diseases in humans, thereby cementing the importance of refining methods for their detection. Despite progress in the field, reliable detection of SVs still remains a problem even for human subjects. Many of the underlying problems that make SVs difficult to detect in humans are amplified in livestock species, whose lower quality genome assemblies and incomplete gene annotation can often give rise to false positive SV discoveries. Regardless of the challenges, SV detection is just as important for livestock researchers as it is for human researchers, given that several productive traits and diseases have been linked to copy number variations (CNVs) in cattle, sheep, and pig. Already, there is evidence that many beneficial SVs have been artificially selected in livestock such as a duplication of the agouti signaling protein gene that causes white coat color in sheep. In this review, we will list current SV and CNV discoveries in livestock and discuss the problems that hinder routine discovery and tracking of these polymorphisms. We will also discuss the impacts of selective breeding on CNV and SV frequencies and mention how SV genotyping could be used in the future to improve genetic selection.

Identification of large selective sweeps associated with major genes in cattle

Selection for new favorable variants can lead to selective sweeps. However, such sweeps might be rare in the evolution of different species for which polygenic adaptation or selection on standing variation might be more common. Still, strong selective sweeps have been described in domestic species such as chicken lines or dog breeds. The goal of our study was to use a panel of individuals from 12 different cattle breeds genotyped at high density (800K SNPs) to perform a whole-genome scan for selective sweeps defined as unexpectedly long stretches of reduced heterozygosity. To that end, we developed a hidden Markov model in which one of the hidden states corresponds to regions of reduced heterozygosity. Some unexpectedly long regions were identified. Among those, six contained genes known to affect traits with simple genetic architecture such as coat color or horn development. However, there was little evidence for sweeps associated with genes underlying production traits.

Colorful DNA polymorphisms in humans

In this review article we summarize current knowledge on how variation on the DNA level influences human pigmentation including color variation of iris, hair, and skin. We review recent progress in the field of human pigmentation genetics by focusing on the genes and DNA polymorphisms discovered to be involved in determining human pigmentation traits, their association with diseases particularly skin cancers, and their power to predict human eye, hair, and skin colors with potential utilization in forensic investigations.

Genes controlling vaccine responses and disease resistance to respiratory viral pathogens in cattle

Farm animals remain at risk of endemic, exotic and newly emerging viruses. Vaccination is often promoted as the best possible solution, and yet for many pathogens, either there are no appropriate vaccines or those that are available are far from ideal. A complementary approach to disease control may be to identify genes and chromosomal regions that underlie genetic variation in disease resistance and response to vaccination. However, identification of the causal polymorphisms is not straightforward as it generally requires large numbers of animals with linked phenotypes and genotypes. Investigation of genes underlying complex traits such as resistance or response to viral pathogens requires several genetic approaches including candidate genes deduced from knowledge about the cellular pathways leading to protection or pathology, or unbiased whole genome scans using markers spread across the genome. Evidence for host genetic variation exists for a number of viral diseases in cattle including bovine respiratory disease and anecdotally, foot and mouth disease virus (FMDV). We immunised and vaccinated a cattle cross herd with a 40-mer peptide derived from FMDV and a vaccine against bovine respiratory syncytial virus (BRSV). Genetic variation has been quantified. A candidate gene approach has grouped high and low antibody and T cell responders by common motifs in the peptide binding pockets of the bovine major histocompatibility complex (BoLA) DRB3 gene. This suggests that vaccines with a minimal number of epitopes that are recognised by most cattle could be designed. Whole genome scans using microsatellite and single nucleotide polymorphism (SNP) markers has revealed many novel quantitative trait loci (QTL) and SNP markers controlling both humoral and cell-mediated immunity, some of which are in genes of known immunological relevance including the toll-like receptors (TLRs). The sequencing, assembly and annotation of livestock genomes and is continuing apace. In addition, provision of high-density SNP chips should make it possible to link phenotypes with genotypes in field populations without the need for structured populations or pedigree information. This will hopefully enable fine mapping of QTL and ultimate identification of the causal gene(s). The research could lead to selection of animals that are more resistant to disease and new ways to improve vaccine efficacy.

Interaction of MC1R and PMEL alleles on solid coat colors in Highland cattle

Six solid colors occur in Highland cattle: black, dun, silver dun and red, yellow, and white. These six coat colors are explained by a non-epistatic interaction of the genotypes at the MC1R and PMEL genes. A three base pair deletion in the PMEL gene leading to the deletion of a leucine from the signal peptide is observed in dilute-colored Highland cattle (c.50_52delTTC, p.Leu18del). The mutant PMEL allele acts in a semi-dominant manner. Dun Galloway cattle also have one copy of the deletion allele, and silver dun Galloway cattle have two copies. The presence of two adjacent leucine residues at the site of this deletion is highly conserved in human, horse, mouse and chicken as well as in cattle with undiluted coat colors. Highland and Galloway cattle thus exhibit a similar dose-dependent dilution effect based on the number of PMEL :c.50_51delTTC alleles, as Charolais cattle with PMEL :c.64G>A alleles. The PMEL :c.64G>A allele was not found in Highland or Galloway cattle.

Information content in genome-wide scans: concordance between patterns of genetic differentiation and linkage mapping associations

BACKGROUND

Scanning the genome with high density SNP markers has become a standard approach for identifying regions of the genome showing substantial between-population genetic differentiation, and thus evidence of diversifying selection. Such regions may contain genes of large phenotypic effect. However, few studies have attempted to address the power or efficacy of such an approach.

RESULTS

In this study, the patterns of allele frequency differences between two cattle breeds based on the Bovine HapMap study were compared with statistical evidence for QTL based on a linkage mapping study of an experimental population formed by a cross between the same breeds. Concordance between the two datasets was seen for chromosomes carrying QTL with strong statistical support, such as BTA5 and BTA18, which carry genes associated with coat color. For these chromosomes, there was a correspondence between the strength of the QTL signal along the chromosome and the degree of genetic differentiation between breeds. However, such an association was not seen in a broader comparison that also included chromosomes carrying QTL with lower significance levels. In addition, other chromosomal regions with substantial QTL effects did not include markers showing extreme between-breed genetic differentiation. Furthermore, the overall consistency between the two studies was weak, with low genome-wide correlation between the statistical values obtained in the linkage mapping study and between-breed genetic differentiation from the HapMap study.

CONCLUSIONS

These results suggest that genomic diversity scans are capable of detecting regions associated with qualitative traits but may be limited in their power to detect regions associated with quantitative phenotypic differences between populations, which may depend on the marker resolution of the study and the level of LD in the populations under investigation.

Quantitative trait loci for variation in immune response to a Foot-and-Mouth Disease virus peptide

Background

Infectious disease of livestock continues to be a cause of substantial economic loss and has adverse welfare consequences in both the developing and developed world. New solutions to control disease are needed and research focused on the genetic loci determining variation in immune-related traits has the potential to deliver solutions. However, identifying selectable markers and the causal genes involved in disease resistance and vaccine response is not straightforward. The aims of this study were to locate regions of the bovine genome that control the immune response post immunisation. 195 F2 and backcross Holstein Charolais cattle were immunised with a 40-mer peptide derived from foot-and-mouth disease virus (FMDV). T cell and antibody (IgG1 and IgG2) responses were measured at several time points post immunisation. All experimental animals (F0, F1 and F2, n = 982) were genotyped with 165 microsatellite markers for the genome scan.

Results

Considerable variability in the immune responses across time was observed and sire, dam and age had significant effects on responses at specific time points. There were significant correlations within traits across time, and between IgG1 and IgG2 traits, also some weak correlations were detected between T cell and IgG2 responses. The whole genome scan detected 77 quantitative trait loci (QTL), on 22 chromosomes, including clusters of QTL on BTA 4, 5, 6, 20, 23 and 25. Two QTL reached 5% genome wide significance (on BTA 6 and 24) and one on BTA 20 reached 1% genome wide significance.

Conclusions

A proportion of the variance in the T cell and antibody response post immunisation with an FDMV peptide has a genetic component. Even though the antigen was relatively simple, the humoral and cell mediated responses were clearly under complex genetic control, with the majority of QTL located outside the MHC locus. The results suggest that there may be specific genes or loci that impact on variation in both the primary and secondary immune responses, whereas other loci may be specifically important for early or later phases of the immune response. Future fine mapping of the QTL clusters identified has the potential to reveal the causal variations underlying the variation in immune response observed.