DNA sequence polymorphisms in a panel of eight candidate bovine imprinted genes and their association with performance traits in Irish Holstein-Friesian cattle

New insight into copy number variations of goat SMAD2 gene and their associations with litter size and semen quality

Copy number variations (CNV) contribute significantly to genetic variations. Numerous studies have shown that CNV affects phenotypic traits in livestock. The SMAD family member 2 (SMAD2) is a leading candidate gene in reproduction and has a crucial effect on litter size. Additionally, SMAD2 is also required for male reproduction and influences male germ cell development. However, there are no reports on investigating the effect of CNVs in the SMAD2 gene on reproductive traits in goat. Therefore, the goal of this study was to explore associations between CNV of the SMAD2 gene and litter size and semen quality in Shaanbei white cashmere (SBWC) goats. In this study, two CNVs within the SMAD2 were identified in 352 SBWC goats (50 males and 302 females). The association analysis revealed that only CNV2 was significantly associated with female goat first-born litter size (P = 3.59 × 10^-4), male semen concentration (P < 0.01), ejaculation volume, live sperm count, and sperm deformity rate (P < 0.05). In terms of phenotypic performance, the individuals with loss genotypes outperformed those with other genotypes. CNV1 and CNV2 genotype combinations containing their dominant genotypes were also associated with goat litter size (P = 1.7 × 10^-5), but no differences in semen quality were found. In summary, CNV2 of the SMAD2 gene is useful for molecular marker-assisted selection breeding, as it is associated with essential goat reproductive traits.

Genome-wide association study revealed suggestive QTLs for production and reproduction traits in Indian Murrah buffalo

The present study was aimed to identify the genome-wide SNPs associated with production and reproduction traits in 96 Indian Murrah buffalo genotyped based on ddRAD approach using Genome-Wide Association Study (GWAS) along with phenotypes of contemporary animals using mixed linear model for production and reproduction traits. A total of 27,735 SNPs identified using ddRAD approach in 96 Indian Murrah buffaloes were used for GWAS. A total of 28 SNPs were found to be associated with production and reproductive traits. Among these, 14 SNPs were present in the intronic region of AK5, BACH2, DIRC2, ECPAS, MPZL1, MYO16, QRFPR, RASGRF1, SLC9A4, TANC1, and TRIM67 genes and one SNP in long non-coding region of LOC102414911. Out of these 28 SNPs, 9 SNPs were found to have pleiotropic effect over milk production traits and were present in chromosome number BBU 1, 2, 4, 6, 9, 10, 12, 19, and 20. SNPs in the intronic region of AK5, TRIM67 genes were found to be associated with milk production traits. Eleven and five SNPs in the intergenic region were associated with milk production and reproduction traits respectively. The above genomic information may be used for selection of Murrah animals for genetic improvement.

An insight into the runs of homozygosity distribution and breed differentiation in Mangalitsa pigs

Mangalitsa pigs exhibit three distinct coat color patterns based on which they are described as Red, Blond, and Swallow-bellied. The current study investigated genome-wide diversity and selection signatures in the three breeds using fixation index, runs of homozygosity and population structure analyses. The analyses were originally based on quality-controlled data on 77 Mangalitsa animals from Germany, including 23 Blond, 30 Swallow-bellied and 24 Red Mangalitsa genotyped with a customized version of the ProcineSNP60 v2 Genotyping Bead Chip. Also, 20 Hungarian Mangalitsa genotypes were included as outgroup data for comparison. Estimates of observed heterozygosity were 0.27, 0.28, and 0.29, and inbreeding coefficients estimated based on runs of homozygosity were 24.11%, 20.82%, and 16.34% for Blond, Swallow-bellied and Red Mangalitsa, respectively. ROH islands were detected in all breeds, however, none of these were shared amongst them. The KIF16B gene previously reported to play a role in synaptic signaling was found in a ROH island (SSC17: 16–26) in Swallow-bellied Mangalitsa. The same gene was found to harbor a significantly differentiated SNP (MARC0032380) while contrasting either Blond or Red to Swallow-belied Mangalitsa. In the Red Mangalitsa, some ROH islands were associated with genes that play a role in meat quality traits, i.e., ABCA12, VIL1, PLSCR5, and USP37. Our population structure analysis highlighted a separation of the three breeds, but also showed the closest relatedness between Red and Blond Mangalitsa pigs. Findings of this study improve our understanding of the diversity in the three breeds of Mangalitsa pigs.

Detection of Genomic Imprinting for Carcass Traits in Cattle Using Imputed High-Density Genotype Data

Genomic imprinting is an epigenetic phenomenon defined as the silencing of an allele, at least partially, at a given locus based on the sex of the transmitting parent. The objective of the present study was to detect the presence of SNP-phenotype imprinting associations for carcass weight (CW), carcass conformation (CC) and carcass fat (CF) in cattle. The data used comprised carcass data, along with imputed, high-density genotype data on 618,837 single nucleotide polymorphisms (SNPs) from 23,687 cattle; all animal genotypes were phased with respect to parent of origin. Based on the phased genotypes and a series of single-locus linear models, 24, 339, and 316 SNPs demonstrated imprinting associations with CW, CC, and CF, respectively. Regardless of the trait in question, no known imprinted gene was located within 0.5 Mb of the SNPs demonstrating imprinting associations in the present study. Since all imprinting associations detected herein were at novel loci, further investigation of these regions may be warranted. Nonetheless, knowledge of these associations might be useful for improving the accuracy of genomic evaluations for these traits, as well as mate allocations systems to exploit the effects of genomic imprinting.

Gene Banks as Reservoirs to Detect Recent Selection: The Example of the Asturiana de los Valles Bovine Breed

Gene banks, framed within the efforts for conserving animal genetic resources to ensure the adaptability of livestock production systems to population growth, income, and climate change challenges, have emerged as invaluable resources for biodiversity and scientific research. Allele frequency trajectories over the few last generations contain rich information about the selection history of populations, which cannot be obtained from classical selection scan approaches based on present time data only. Here we apply a new statistical approach taking advantage of genomic time series and a state of the art statistic (nSL) based on present time data to disentangle both old and recent signatures of selection in the Asturiana de los Valles cattle breed. This local Spanish originally multipurpose breed native to Asturias has been selected for beef production over the last few generations. With the use of SNP chip and whole-genome sequencing (WGS) data, we detect candidate regions under selection reflecting the effort of breeders to produce economically valuable beef individuals, e.g., by improving carcass and meat traits with genes such as MSTN, FLRT2, CRABP2, ZNF215, RBPMS2, OAZ2, or ZNF609, while maintaining the ability to thrive under a semi-intensive production system, with the selection of immune (GIMAP7, GIMAP4, GIMAP8, and TICAM1) or olfactory receptor (OR2D2, OR2D3, OR10A4, and 0R6A2) genes. This kind of information will allow us to take advantage of the invaluable resources provided by gene bank collections from local less competitive breeds, enabling the livestock industry to exploit the different mechanisms fine-tuned by natural and human-driven selection on different populations to improve productivity.

Selection signatures in two oldest Russian native cattle breeds revealed using high-density single nucleotide polymorphism analysis

Native cattle breeds can carry specific signatures of selection reflecting their adaptation to the local environmental conditions and response to the breeding strategy used. In this study, we comprehensively analysed high-density single nucleotide polymorphism (SNP) genotypes to characterise the population structure and detect the selection signatures in Russian native Yaroslavl and Kholmogor dairy cattle breeds, which have been little influenced by introgression with transboundary breeds. Fifty-six samples of pedigree-recorded purebred animals, originating from different breeding farms and representing different sire lines, of the two studied breeds were genotyped using a genome-wide bovine genotyping array (Bovine HD BeadChip). Three statistical analyses—calculation of fixation index (F_ST) for each SNP for the comparison of the pairs of breeds, hapFLK analysis, and estimation of the runs of homozygosity (ROH) islands shared in more than 50% of animals—were combined for detecting the selection signatures in the genome of the studied cattle breeds. We confirmed nine and six known regions under putative selection in the genomes of Yaroslavl and Kholmogor cattle, respectively; the flanking positions of most of these regions were elucidated. Only two of the selected regions (localised on BTA 14 at 24.4–25.1 Mbp and on BTA 16 at 42.5–43.5 Mb) overlapped in Yaroslavl, Kholmogor and Holstein breeds. In addition, we detected three novel selection sweeps in the genome of Yaroslavl (BTA 4 at 4.74–5.36 Mbp, BTA 15 at 17.80–18.77 Mbp, and BTA 17 at 45.59–45.61 Mbp) and Kholmogor breeds (BTA 12 at 82.40–81.69 Mbp, BTA 15 at 16.04–16.62 Mbp, and BTA 18 at 0.19–1.46 Mbp) by using at least two of the above-mentioned methods. We expanded the list of candidate genes associated with the selected genomic regions and performed their functional annotation. We discussed the possible involvement of the identified candidate genes in artificial selection in connection with the origin and development of the breeds. Our findings on the Yaroslavl and Kholmogor breeds obtained using high-density SNP genotyping and three different statistical methods allowed the detection of novel putative genomic regions and candidate genes that might be under selection. These results might be useful for the sustainable development and conservation of these two oldest Russian native cattle breeds.

Genetic Architecture of Carcass and Meat Quality Traits in Montana Tropical® Composite Beef Cattle

The Montana Tropical^® Composite is a recently developed beef cattle population that is rapidly expanding in Brazil and other tropical countries. This is mainly due to its improved meat quality and adaptation to tropical climate conditions compared to Zebu and Taurine cattle breeds, respectively. This study aimed to investigate the genetic architecture of ultrasound-based carcass and meat quality traits in Montana Tropical^® Composite beef cattle. Therefore, we estimated variance components and genetic parameters and performed genome-wide association studies using the weighted single-step Genomic Best Linear Unbiased Prediction (GBLUP) approach. A pedigree dataset containing 28,480 animals was used, in which 1,436 were genotyped using a moderate-density Single Nucleotide Polymorphism panel (30K; 30,105 SNPs). A total of 9,358, 5,768, 7,996, and 1,972 phenotypic records for the traits Longissimus muscle area (LMA), backfat thickness (BFT), rump fat thickness (RFT), and for marbling score (MARB), respectively, were used for the analyses. Moderate to high heritability estimates were obtained and ranged from 0.16 ± 0.03 (RFT) to 0.33 ± 0.05 (MARB). A high genetic correlation was observed between BFT and RFT (0.97 ± 0.02), suggesting that a similar set of genes affects both traits. The most relevant genomic regions associated with LMA, BFT, RFT, and MARB were found on BTA10 (5.4–5.8 Mb), BTA27 (25.2–25.5 Mb), BTA18 (60.6–61.0 Mb), and BTA21 (14.8–15.4 Mb). Two overlapping genomic regions were identified for RFT and MARB (BTA13:47.9–48.1 Mb) and for BFT and RFT (BTA13:61.5–62.3 Mb). Candidate genes identified in this study, including PLAG1, LYN, WWOX, and PLAGL2, were previously reported to be associated with growth, stature, skeletal muscle growth, fat thickness, and fatty acid composition. Our results indicate that ultrasound-based carcass and meat quality traits in the Montana Tropical^® Composite beef cattle are heritable, and therefore, can be improved through selective breeding. In addition, various novel and already known genomic regions related to these traits were identified, which contribute to a better understanding of the underlying genetic background of LMA, BFT, RFT, and MARB in the Montana Tropical Composite population.

Metabolism-associated genome-wide epigenetic changes in bovine oocytes during early lactation

Dietary intake in early lactating cows is outmatched by milk production. These cows experience a negative energy balance, resulting in a distinct blood metabolism and poor reproductive function due to impaired ovulation and increased embryo loss. We hypothesize that oocytes from lactating cows undergoing transient metabolic stress exhibit a different epigenetic profile crucial for developmental competence. To investigate this, we collected oocytes from metabolically-profiled cows at early- and mid-postpartum stages and characterized their epigenetic landscape compared with control heifers using whole-genome bisulfite sequencing. Early-postpartum cows were metabolically deficient with a significantly lower energy balance and significantly higher concentrations of non-esterified fatty acids and beta-hydroxybutyrate than mid-postpartum animals and control heifers. Accordingly, 32,990 early-postpartum-specific differentially methylated regions (DMRs) were found in genes involved in metabolic pathways, carbon metabolism, and fatty acid metabolism, likely descriptive of the epigenetic regulation of metabolism in early-postpartum oocytes. DMRs found overlapping CpG islands and exons of imprinted genes such as MEST and GNAS in early-postpartum oocytes suggest that early lactation metabolic stress may affect imprint acquisition, which could explain the embryo loss. This whole-genome approach introduces potential candidate genes governing the link between metabolic stress and the reproductive outcome of oocytes.

Copy number variation detection in Chinese indigenous cattle by whole genome sequencing

Copy number variation (CNV) refers to a kind of structural variation, having functional and evolutionary effects on phenotypes. Thus far, further elucidation of the CNVs in different Chinese indigenous cattle breeds by whole genome sequencing have yet not been done. In this study, a comprehensive genomic analysis was performed on 75 cattle individuals including six Chinese indigenous cattle breeds and two non-native specialized beef cattle breeds. Based on the 11,486 CNVRs discovered, population analysis was performed, showed that all the cattle breeds clustered in to three clades, consistent with their lineages Bos taurus, Bos taurus × Bos indicus and Bos indicus. Importantly, a set of CNVRs related genes were found to be associated with the traits of interest, which include meat production or quality (CAST, ACTC1, etc.), adaption (BLA-DQB, EGLN2, etc.) and coat color (KIT, MITF, etc.). These results provide valuable full genome variation resources for Chinese bovine genome research and would be helpful for cattle breeding and selection programs in the future.

Genetic and epigenetic architecture of paternal origin contribute to gestation length in cattle

The length of gestation can affect offspring health and performance. Both maternal and fetal effects contribute to gestation length; however, paternal contributions to gestation length remain elusive. Using genome-wide association study (GWAS) in 27,214 Holstein bulls with millions of gestation records, here we identify nine paternal genomic loci associated with cattle gestation length. We demonstrate that these GWAS signals are enriched in pathways relevant to embryonic development, and in differentially methylated regions between sperm samples with long and short gestation length. We reveal that gestation length shares genetic and epigenetic architecture in sperm with calving ability, body depth, and conception rate. While several candidate genes are detected in our fine-mapping analysis, we provide evidence indicating ZNF613 as a promising candidate for cattle gestation length. Collectively, our findings support that the paternal genome and epigenome can impact gestation length potentially through regulation of the embryonic development.

Lingzhao Fang et al. studied the paternal genetic variants that affect gestational length in cattle. They found that paternal genes from pathways involved in embryonic development were associated with gestation length, and that these were often found in differentially methylated regions of the genome.

Scanning the genomes of parents for imprinted loci acting in their un-genotyped progeny

Depending on their parental origin, alleles at imprinted loci are fully or partially inactivated through epigenetic mechanisms. Their effects contribute to the broader class of parent-of-origin effects. Standard methodology for mapping imprinted quantitative trait loci in association studies requires phenotypes and parental origin of marker alleles (ordered genotypes) to be simultaneously known for each individual. As such, many phenotypes are known from un-genotyped offspring in ongoing breeding programmes (e.g. meat animals), while their parents have known genotypes but no phenotypes. By theoretical considerations and simulations, we showed that the limitations of standard methodology can be overcome in such situations. This is achieved by first estimating parent-of-origin effects, which then serve as dependent variables in association analyses, in which only imprinted loci give a signal. As a theoretical foundation, the regression of parent-of-origin effects on the number of B-alleles at a biallelic locus — representing the un-ordered genotype — equals the imprinting effect. The applicability to real data was demonstrated for about 1800 genotyped Brown Swiss bulls and their un-genotyped fattening progeny. Thus, this approach unlocks vast data resources in various species for imprinting analyses and offers valuable clues as to what extent imprinted loci contribute to genetic variability.

Probe-based association analysis identifies several deletions associated with average daily gain in beef cattle

BACKGROUND

Average daily gain (ADG) is an important trait that contributes to the production efficiency and economic benefits in the beef cattle industry. The molecular mechanisms of ADG have not yet been fully explored because most recent association studies for ADG are based on SNPs or haplotypes. We reported a systematic CNV discovery and association analysis for ADG in Chinese Simmental beef cattle.

RESULTS

Our study identified 4912 nonredundant CNVRs with a total length of ~ 248.7 Mb, corresponding to ~ 8.9% of the cattle genome. Using probe-based CNV association, we identified 24 and 12 significant SNP probes within five deletions and two duplications for ADG, respectively. Among them, we found one common deletion with 89 kb imbedded in LHFPL Tetraspan Subfamily Member 6 (LHFPL6) at 22.9 Mb on BTA12, which has high frequency (12.9%) dispersing across population. CNV selection test using V_ST statistic suggested this common deletion may be under positive selection in Chinese Simmental cattle. Moreover, this deletion was not overlapped with any candidate SNP for ADG compared with previous SNPs-based association studies, suggesting its important role for ADG. In addition, we identified one rare deletion near gene Growth Factor Receptor-bound Protein 10 (GRB10) at 5.1 Mb on BTA4 for ADG using both probe-based association and region-based approaches.

CONCLUSIONS

Our results provided some valuable insights to elucidate the genetic basis of ADG in beef cattle, and these findings offer an alternative perspective to understand the genetic mechanism of complex traits in terms of copy number variations in farm animals.

Human-Mediated Introgression of Haplotypes in a Modern Dairy Cattle Breed

Domestic animals can serve as model systems of adaptive introgression and their genomic signatures. In part, their usefulness as model systems is due to their well-known histories. Different breeding strategies such as introgression and artificial selection have generated numerous desirable phenotypes and superior performance in domestic animals. The modern Danish Red Dairy Cattle is studied as an example of an introgressed population. It originates from crossing the traditional Danish Red Dairy Cattle with the Holstein and Brown Swiss breeds, both known for high milk production. This crossing happened, among other things due to changes in the production system, to raise milk production and overall performance. The genomes of modern Danish Red Dairy Cattle are heavily influenced by regions introgressed from the Holstein and Brown Swiss breeds and under subsequent selection in the admixed population. The introgressed proportion of the genome was found to be highly variable across the genome. Haplotypes introgressed from Holstein and Brown Swiss contained or overlapped known genes affecting milk production, as well as protein and fat content (CD14, ZNF215, BCL2L12, and THRSP for Holstein origin and ITPR2, BCAT1, LAP3, and MED28 for Brown Swiss origin). Genomic regions with high introgression signals also contained genes and enriched QTL associated with calving traits, body confirmation, feed efficiency, carcass, and fertility traits. These introgressed signals with relative identity-by-descent scores larger than the median showing Holstein or Brown Swiss introgression are mostly significantly correlated with the corresponding test statistics from signatures of selection analyses in modern Danish Red Dairy Cattle. Meanwhile, the putative significant introgressed signals have a significant dependency with the putative significant signals from signatures of selection analyses. Artificial selection has played an important role in the genomic footprints of introgression in the genome of modern Danish Red Dairy Cattle. Our study on a modern cattle breed contributes to an understanding of genomic consequences of selective introgression by demonstrating the extent to which adaptive effects contribute to shape the specific genomic consequences of introgression.

Review: Genomics of bull fertility

Fertility is one of the most economically important traits in both beef and dairy cattle production; however, only female fertility is typically subjected to selection. Male and female fertility have only a small positive genetic correlation which is likely due to the existence of a relatively small number of genetic variants within each breed that cause embryonic and developmental losses. Genomic tools have been developed that allow the identification of lethal recessive loci based upon marker haplotypes. Selection against haplotypes harbouring lethal alleles in conjunction with selection to improve female fertility will result in an improvement in male fertility. Genomic selection has resulted in a two to fourfold increase in the rate of genetic improvement of most dairy traits in US Holstein cattle, including female fertility. Considering the rapidly increasing rate of adoption of high-throughput single nucleotide polymorphism genotyping in both the US dairy and beef industries, genomic selection should be the most effective of all currently available approaches to improve male fertility. However, male fertility phenotypes are not routinely recorded in natural service mating systems and when artificial insemination is used, semen doses may be titrated to lower post-thaw progressively motile sperm numbers for high-merit and high-demand bulls. Standardization of sperm dosages across bull studs for semen distributed from young bulls would allow the capture of sire conception rate phenotypes for young bulls that could be used to generate predictions of genetic merit for male fertility in both males and females. These data would allow genomic selection to be implemented for male fertility in addition to female fertility within the US dairy industry. While the rate of use of artificial insemination is much lower within the US beef industry, the adoption of sexed semen in the dairy industry has allowed dairy herds to select cows from which heifer replacements are produced and cows that are used to produce terminal crossbred bull calves sired by beef breed bulls. Capture of sire conception rate phenotypes in dairy herds utilizing sexed semen will contribute data enabling genomic selection for male fertility in beef cattle breeds. As the commercial sector of the beef industry increasingly adopts fixed-time artificial insemination, sire conception rate phenotypes can be captured to facilitate the development of estimates of genetic merit for male fertility within US beef breeds.

Regulation of DNA methylation on EEF1D and RPL8 expression in cattle

Dynamic changes to the epigenome play a critical role in a variety of biology processes and complex traits. Many important candidate genes have been identified through our previous genome wide association study (GWAS) on milk production traits in dairy cattle. However, the underlying mechanism of candidate genes have not yet been clearly understood. In this study, we analyzed the methylation variation of the candidate genes, EEF1D and RPL8, which were identified to be strongly associated with milk production traits in dairy cattle in our previous studies, and its effect on protein and mRNA expression. We compared DNA methylation profiles and gene expression levels of EEF1D and RPL8 in five different tissues (heart, liver, mammary gland, ovary and muscle) of three cows. Both genes showed the highest expression level in mammary gland. For RPL8, there was no difference in the DNA methylation pattern in the five tissues, suggesting no effect of DNA methylation on gene expression. For EEF1D, the DNA methylation levels of its first CpG island differed in the five tissues and were negatively correlated with the gene expression levels. To further investigate the function of DNA methylation on the expression of EEF1D, we collected blood samples of three cows at early stage of lactation and in dry period and analyzed its expression and the methylation status of the first CpG island in blood. As a result, the mRNA expression of EEF1D in the dry period was higher than that at the early stage of lactation, while the DNA methylation level in the dry period was lower than that at the early stage of lactation. Our result suggests that the DNA methylation of EEF1D plays an important role in the spatial and temporal regulation of its expression and possibly have an effect on the milk production traits.

Correlating single nucleotide polymorphisms in the myostatin gene with performance traits in rabbit

The Myostatin (MSTN), or Growth and Differentiation Factor 8 (GDF8), gene has been implicated in the double muscling phenomenon, in which a series of mutations render the gene inactive and unable to properly regulate muscle fibre deposition. Single nucleotide polymorphisms (SNPs) in the MSTN gene have been correlated to production traits, making it a candidate target gene to enhance livestock and fowl productivity. This study aimed to assess any association of three SNPs in the rabbit MSTN gene (c.713T&gt;A in exon 2, c.747+34C&gt;T in intron 2, and c.*194A&gt;G in 3’-untranslated region) and their combinations, with carcass, production and reproductive traits. The investigated traits included individual body weight, daily body weight gain, carcass traits and reproductive traits. The 3 SNPs were screened using PCR-restriction fragment length polymorphism (RFLP)-based analysis and the effects of the different SNP genotypes and their combinations were estimated in a rabbit population. Additionally, additive and dominance effects were estimated for significant traits. The results found no significant association between the c.713 T&gt;A SNP and all the examined traits. Allele T at the c.747+34C&gt;T SNP was only significantly associated (P&lt;0.05) with increased body weight at 12 wk of age. However, for the SNP residing in the 3’ untranslated region (c.*194A&gt;G), allele G was significantly associated (P&lt;0.05) with increased body weight and high growth rate. Genotype GG at the c.*194A&gt;G SNP also had positive effects on most carcass traits. The estimated additive genetic effect for the c.*194A&gt;G SNP was significant (P&lt;0.05) with most body weight, daily gain and carcass traits. No significant association was obtained between any MSTN SNPs and reproductive traits. In the combinations analysis, regardless of the genotypes of SNPs at c.713T&gt;A and c.747+34C&gt;T, GG at the c.*194A&gt;G SNP correlated with highest values in body weight and daily weight gain. In conclusion, the ‘G’ allele at the c.*194A&gt;G SNP had positive effects on growth and carcass traits and so could be used as a favourable allele in planning rabbit selection. Further population-wide studies are necessary to test the association of the c.*194A&gt;G SNP with carcass traits. We also recommend evaluation of the potential effects of the c.*194A&gt;G SNP on MSTN gene expression.

Use of canonical discriminant analysis to study signatures of selection in cattle

BACKGROUND

Cattle include a large number of breeds that are characterized by marked phenotypic differences and thus constitute a valuable model to study genome evolution in response to processes such as selection and domestication. Detection of "signatures of selection" is a useful approach to study the evolutionary pressures experienced throughout history. In the present study, signatures of selection were investigated in five cattle breeds farmed in Italy using a multivariate approach.

METHODS

A total of 4094 bulls from five breeds with different production aptitudes (two dairy breeds: Italian Holstein and Italian Brown Swiss; two beef breeds: Piemontese and Marchigiana; and one dual purpose breed: Italian Simmental) were genotyped using the Illumina BovineSNP50 v.1 beadchip. Canonical discriminant analysis was carried out on the matrix of single nucleotide polymorphisms (SNP) genotyping data, separately for each chromosome. Scores for each canonical variable were calculated and then plotted in the canonical space to quantify the distance between breeds. SNPs for which the correlation with the canonical variable was in the 99th percentile for a specific chromosome were considered to be significantly associated with that variable. Results were compared with those obtained using an FST-based approach.

RESULTS

Based on the results of the canonical discriminant analysis, a large number of signatures of selection were detected, among which several had strong signals in genomic regions that harbour genes known to have an impact on production and morphological bovine traits, including MSTN, LCT, GHR, SCD, NCAPG, KIT, and ASIP. Moreover, new putative candidate genes were identified, such as GCK, B3GALNT1, MGAT1, GALNTL1, PRNP, and PRND. Similar results were obtained with the FST-based approach.

CONCLUSIONS

The use of canonical discriminant analysis on 50 K SNP genotypes allowed the extraction of new variables that maximize the separation between breeds. This approach is quite straightforward, it can compare more than two groups simultaneously, and relative distances between breeds can be visualized. The genes that were highlighted in the canonical discriminant analysis were in concordance with those obtained using the FST index.

Epigenetics and developmental programming of welfare and production traits in farm animals

The concept that postnatal health and development can be influenced by events that occur in utero originated from epidemiological studies in humans supported by numerous mechanistic (including epigenetic) studies in a variety of model species. Referred to as the 'developmental origins of health and disease' or 'DOHaD' hypothesis, the primary focus of large-animal studies until quite recently had been biomedical. Attention has since turned towards traits of commercial importance in farm animals. Herein we review the evidence that prenatal risk factors, including suboptimal parental nutrition, gestational stress, exposure to environmental chemicals and advanced breeding technologies, can determine traits such as postnatal growth, feed efficiency, milk yield, carcass composition, animal welfare and reproductive potential. We consider the role of epigenetic and cytoplasmic mechanisms of inheritance, and discuss implications for livestock production and future research endeavours. We conclude that although the concept is proven for several traits, issues relating to effect size, and hence commercial importance, remain. Studies have also invariably been conducted under controlled experimental conditions, frequently assessing single risk factors, thereby limiting their translational value for livestock production. We propose concerted international research efforts that consider multiple, concurrent stressors to better represent effects of contemporary animal production systems.

Plasma exosome profiles from dairy cows with divergent fertility phenotypes

Cell-to-cell communication in physiological and pathological conditions may be influenced by neighboring cells, distant tissues, or local environmental factors. Exosomes are specific subsets of extracellular vesicles that internalize and deliver their content to near and distant sites. Exosomes may play a role in the maternal-embryo crosstalk vital for the recognition and maintenance of a pregnancy; however, their role in dairy cow reproduction has not been established. This study aimed to characterize the exosome profile in the plasma of 2 strains of dairy cow with divergent fertility phenotypes. Plasma was obtained and characterized on the basis of genetic ancestry as fertile (FERT; <23% North American genetics, New Zealand Holstein-Friesian strain, n=8) or subfertile (SUBFERT; >92% North American genetics, North American Holstein-Friesian strain, n=8). Exosomes were isolated by differential and buoyant density centrifugation and characterized by size distribution (nanoparticle tracking analysis, NanoSight NS500, NanoSight Ltd., Amesbury, UK), the presence of CD63 (Western blot), and their morphology (electron microscopy). The total number of exosomes was determined by quantifying the immunoreactive CD63 (ExoELISA kit, System Biosciences), and the protein content established by mass spectrometry. Enriched exosome fractions were identified as cup-shape vesicles with diameters around 100 nm and positive for the CD63 marker. The concentration of exosomes was 50% greater in FERT cows. Mass spectrometry identified 104 and 117 proteins in FERT and SUBFERT cows, of which 23 and 36 were unique, respectively. Gene ontology analysis revealed enrichment for proteins involved in immunomodulatory processes and cell-to-cell communication. Although the role of exosomes in dairy cow reproduction remains to be elucidated, their quantification and content in models with divergent fertility phenotypes could provide novel information to support both physiological and genetic approaches to improving dairy cow fertility.

Genome-wide association between single nucleotide polymorphisms with beef fatty acid profile in Nellore cattle using the single step procedure

Background

Saturated fatty acids can be detrimental to human health and have received considerable attention in recent years. Several studies using taurine breeds showed the existence of genetic variability and thus the possibility of genetic improvement of the fatty acid profile in beef. This study identified the regions of the genome associated with saturated, mono- and polyunsaturated fatty acids, and n-6 to n-3 ratios in the Longissimus thoracis of Nellore finished in feedlot, using the single-step method.

Results

The results showed that 115 windows explain more than 1 % of the additive genetic variance for the 22 studied fatty acids. Thirty-one genomic regions that explain more than 1 % of the additive genetic variance were observed for total saturated fatty acids, C12:0, C14:0, C16:0 and C18:0. Nineteen genomic regions, distributed in sixteen different chromosomes accounted for more than 1 % of the additive genetic variance for the monounsaturated fatty acids, such as the sum of monounsaturated fatty acids, C14:1 cis-9, C18:1 trans-11, C18:1 cis-9, and C18:1 trans-9. Forty genomic regions explained more than 1 % of the additive variance for the polyunsaturated fatty acids group, which are related to the total polyunsaturated fatty acids, C20:4 n-6, C18:2 cis-9 cis12 n-6, C18:3 n-3, C18:3 n-6, C22:6 n-3 and C20:3 n-6 cis-8 cis-11 cis-14. Twenty-one genomic regions accounted for more than 1 % of the genetic variance for the group of omega-3, omega-6 and the n-6:n-3 ratio.

Conclusions

The identification of such regions and the respective candidate genes, such as ELOVL5, ESSRG, PCYT1A and genes of the ABC group (ABC5, ABC6 and ABC10), should contribute to form a genetic basis of the fatty acid profile of Nellore (Bos indicus) beef, contributing to better selection of the traits associated with improving human health.

Electronic supplementary material

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

Use of Genomic Tools to Improve Cattle Health in the Context of Infectious Diseases

Although infectious diseases impose a heavy economic burden on the cattle industry, the etiology of many disorders that affect livestock is not fully elucidated, and effective countermeasures are often lacking. The main tools available until now have been vaccines, antibiotics and antiparasitic drugs. Although these have been very successful in some cases, the appearance of parasite and microbial resistance to these treatments is a cause of concern. Next-generation sequencing provides important opportunities to tackle problems associated with pathogenic illnesses. This review describes the rapid gains achieved to track disease progression, identify the pathogens involved, and map pathogen interactions with the host. Use of novel genomic tools subsequently aids in treatment development, as well as successful creation of breeding programs aimed toward less susceptible livestock. These may be important tools for mitigating the long term effects of combating infection and helping reduce the reliance on antibiotic treatment.

Prospects and challenges for the conservation of farm animal genomic resources, 2015-2025

Livestock conservation practice is changing rapidly in light of policy developments, climate change and diversifying market demands. The last decade has seen a step change in technology and analytical approaches available to define, manage and conserve Farm Animal Genomic Resources (FAnGR). However, these rapid changes pose challenges for FAnGR conservation in terms of technological continuity, analytical capacity and integrative methodologies needed to fully exploit new, multidimensional data. The final conference of the ESF Genomic Resources program aimed to address these interdisciplinary problems in an attempt to contribute to the agenda for research and policy development directions during the coming decade. By 2020, according to the Convention on Biodiversity's Aichi Target 13, signatories should ensure that “…the genetic diversity of …farmed and domesticated animals and of wild relatives …is maintained, and strategies have been developed and implemented for minimizing genetic erosion and safeguarding their genetic diversity.” However, the real extent of genetic erosion is very difficult to measure using current data. Therefore, this challenging target demands better coverage, understanding and utilization of genomic and environmental data, the development of optimized ways to integrate these data with social and other sciences and policy analysis to enable more flexible, evidence-based models to underpin FAnGR conservation. At the conference, we attempted to identify the most important problems for effective livestock genomic resource conservation during the next decade. Twenty priority questions were identified that could be broadly categorized into challenges related to methodology, analytical approaches, data management and conservation. It should be acknowledged here that while the focus of our meeting was predominantly around genetics, genomics and animal science, many of the practical challenges facing conservation of genomic resources are societal in origin and are predicated on the value (e.g., socio-economic and cultural) of these resources to farmers, rural communities and society as a whole. The overall conclusion is that despite the fact that the livestock sector has been relatively well-organized in the application of genetic methodologies to date, there is still a large gap between the current state-of-the-art in the use of tools to characterize genomic resources and its application to many non-commercial and local breeds, hampering the consistent utilization of genetic and genomic data as indicators of genetic erosion and diversity. The livestock genomic sector therefore needs to make a concerted effort in the coming decade to enable to the democratization of the powerful tools that are now at its disposal, and to ensure that they are applied in the context of breed conservation as well as development.

Genomic imprinting effects on complex traits in domesticated animal species

Monoallelically expressed genes that exert their phenotypic effect in a parent-of-origin specific manner are considered to be subject to genomic imprinting, the most well understood form of epigenetic regulation of gene expression in mammals. The observed differences in allele specific gene expression for imprinted genes are not attributable to differences in DNA sequence information, but to specific chemical modifications of DNA and chromatin proteins. Since the discovery of genomic imprinting some three decades ago, over 100 imprinted mammalian genes have been identified and considerable advances have been made in uncovering the molecular mechanisms regulating imprinted gene expression. While most genomic imprinting studies have focused on mouse models and human biomedical disorders, recent work has highlighted the contributions of imprinted genes to complex trait variation in domestic livestock species. Consequently, greater understanding of genomic imprinting and its effect on agriculturally important traits is predicted to have major implications for the future of animal breeding and husbandry. In this review, we discuss genomic imprinting in mammals with particular emphasis on domestic livestock species and consider how this information can be used in animal breeding research and genetic improvement programs.

Alternative parameterizations of relatedness in whole genome association analysis of pre-weaning traits of Nelore-Angus calves

Gestation length, birth weight, and weaning weight of F₂ Nelore-Angus calves (n = 737) with designed extensive full-sibling and half-sibling relatedness were evaluated for association with 34,957 SNP markers. In analyses of birth weight, random relatedness was modeled three ways: 1) none, 2) random animal, pedigree-based relationship matrix, or 3) random animal, genomic relationship matrix. Detected birth weight-SNP associations were 1,200, 735, and 31 for those parameterizations respectively; each additional model refinement removed associations that apparently were a result of the built-in stratification by relatedness. Subsequent analyses of gestation length and weaning weight modeled genomic relatedness; there were 40 and 26 trait-marker associations detected for those traits, respectively. Birth weight associations were on BTA14 except for a single marker on BTA5. Gestation length associations included 37 SNP on BTA21, 2 on BTA27 and one on BTA3. Weaning weight associations were on BTA14 except for a single marker on BTA10. Twenty-one SNP markers on BTA14 were detected in both birth and weaning weight analyses.

Genome wide CNV analysis reveals additional variants associated with milk production traits in Holsteins

Background

Milk production is an economically important sector of global agriculture. Much attention has been paid to the identification of quantitative trait loci (QTL) associated with milk, fat, and protein yield and the genetic and molecular mechanisms underlying them. Copy number variation (CNV) is an emerging class of variants which may be associated with complex traits.

Results

In this study, we performed a genome-wide association between CNVs and milk production traits in 26,362 Holstein bulls and cows. A total of 99 candidate CNVs were identified using Illumina BovineSNP50 array data, and association tests for each production trait were performed using a linear regression analysis with PCA correlation. A total of 34 CNVs on 22 chromosomes were significantly associated with at least one milk production trait after false discovery rate (FDR) correction. Some of those CNVs were located within or near known QTL for milk production traits. We further investigated the relationship between associated CNVs with neighboring SNPs. For all 82 combinations of traits and CNVs (less than 400 kb in length), we found 17 cases where CNVs directly overlapped with tag SNPs and 40 cases where CNVs were adjacent to tag SNPs. In 5 cases, CNVs located were in strong linkage disequilibrium with tag SNPs, either within or adjacent to the same haplotype block. There were an additional 20 cases where CNVs did not have a significant association with SNPs, suggesting that the effects of those CNVs were probably not captured by tag SNPs.

Conclusion

We conclude that combining CNV with SNP analyses reveals more genetic variations underlying milk production traits than those revealed by SNPs alone.

Electronic supplementary material

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

Imprinted loci in domestic livestock species as epigenomic targets for artificial selection of complex traits

The phenomenon of genomic imprinting, whereby a subset of mammalian genes display parent-of-origin-specific monoallelic expression, is one of the most active areas of epigenetics research. Over the past two decades, more than 100 imprinted mammalian genes have been identified, while considerable advances have been made in elucidating the molecular mechanisms governing imprinting. These studies have helped to unravel the epigenome--a separate layer of regulatory information contained in eukaryotic chromosomes that influences gene expression and phenotypes without involving changes to the underlying DNA sequence. Although most studies of genomic imprinting in mammals have focussed on mouse models or human biomedical disorders, there is burgeoning interest in the phenotypic effects of imprinted genes in domestic livestock species. In particular, research has focused on imprinted genes influencing foetal growth and development, which are associated with economically important production traits in cattle, sheep and pigs. These findings, when coupled with the data emerging from the various different livestock genome projects, have major implications for the future of animal breeding, health and management. Here, we review current scientific knowledge regarding genomic imprinting in livestock species and evaluate how this information can be used in modern livestock improvement programmes.

Comparative epigenetics: relevance to the regulation of production and health traits in cattle

With the development of genomic, transcriptomic and bioinformatic tools, recent advances in molecular technologies have significantly impacted bovine bioscience research and are revolutionising animal selection and breeding. Integration of epigenetic information represents yet another challenging molecular frontier. Epigenetics is the study of biochemical modifications to DNA and to histones, the proteins that provide stability to DNA. These epigenetic changes are induced by environmental stimuli; they alter gene expression and are potentially heritable. Epigenetics research holds the key to understanding how environmental factors contribute to phenotypic variation in traits of economic importance in cattle including development, nutrition, behaviour and health. In this review, we discuss the potential applications of epigenetics in bovine research, using breakthroughs in human and murine research to signpost the way.

Genomic imprinting and parent-of-origin effects on complex traits

Parent-of-origin effects occur when the phenotypic effect of an allele depends on whether it is inherited from the mother or the father. Several phenomena can cause parent-of-origin effects, but the best characterized is parent-of-origin-dependent gene expression associated with genomic imprinting. The development of new mapping approaches applied to the growing abundance of genomic data has demonstrated that imprinted genes can be important contributors to complex trait variation. Therefore, to understand the genetic architecture and evolution of complex traits, including complex diseases and traits of agricultural importance, it is crucial to account for these parent-of-origin effects. Here, we discuss patterns of phenotypic variation associated with imprinting, evidence supporting its role in complex trait variation and approaches for identifying its molecular signatures.

Maternal Diet during Pregnancy Induces Gene Expression and DNA Methylation Changes in Fetal Tissues in Sheep

Studies in rats and mice have established that maternal nutrition induces epigenetic modifications, sometimes permanently, that alter gene expression in the fetus, which in turn leads to phenotypic changes. However, limited data is available on the influence of maternal diet on epigenetic modifications and gene expression in sheep. Therefore, the objectives of this study were to investigate the impact of different maternal dietary energy sources on the expression of imprinted genes in fetuses in sheep. Ewes were naturally bred to a single sire and from days 67 ± 3 of gestation until necropsy (days 130 ± 1), they were fed one of three diets of alfalfa haylage (HY; fiber), corn (CN; starch), or dried corn distiller’s grains (DG; fiber plus protein plus fat). A total of 26 fetuses were removed from the dams and longissimus dorsi, semitendinosus, perirenal adipose depot, and subcutaneous adipose depot tissues were collected for expression and DNA methylation analyses. Expression analysis of nine imprinted genes and three DNA methyltransferase (DNMTs) genes showed significant effects of the different maternal diets on the expression of these genes. The methylation levels of CpG islands of both IGF2R and H19 were higher in HY and DG than CN fetuses in both males and females. This result is consistent with the low amino acid content of the CN diet, a source of methyl group donors, compared to HY and DG diets. Thus, results of this study provide evidence of association between maternal nutrition during pregnancy and transcriptomic and epigenomic alterations of imprinted genes and DNMTs in the fetal tissues.

Association of single nucleotide polymorphic sites in candidate genes with aggressiveness and deoxynivalenol production in Fusarium graminearum causing wheat head blight

BACKGROUND

Fusarium graminearum sensu stricto (s.s.) is an ubiquitous pathogen of cereals. The economic impact of Fusarium head blight (FHB) is characterized by crop losses and mycotoxin contamination. Our objective was to associate SNP diversity within candidate genes with phenotypic traits. A total of 77 F. graminearum s.s. isolates was tested for severity of fungal infection (= aggressiveness) and deoxynivalenol (DON) production in an inoculated field experiment at two locations in each of two years. For seven genes known to control fungal growth (MetAP1, Erf2) or DON production (TRI1, TRI5, TRI6 TRI10 and TRI14) single nucleotides polymorphic sites (SNPs) were determined and evaluated for the extent of linkage disequilibrium (LD). Associations of SNPs with both phenotypic traits were tested using linear mixed models.

RESULTS

Decay of LD was in most instances fast. Two neighboring SNPs in MetAP1 and one SNP in Erf2 were significantly (P < 0.05) associated with aggressiveness explaining proportions of genotypic variance (pG) of 25.6%, 0.5%, and 13.1%, respectively. One SNP in TRI1 was significantly associated with DON production (pG = 4.4).

CONCLUSIONS

We argue that using the published sequence information of Fusarium graminearum as a template to amplify comparative sequence parts of candidate genes is an effective method to detect quantitative trait loci. Our findings underline the potential of candidate gene association mapping approaches to identify functional SNPs underlying aggressiveness and DON production for F. graminearum s.s populations.

Gene polymorphisms: the keys for marker assisted selection and unraveling core regulatory pathways for mastitis resistance

One of the most frequent mammary diseases impacting lactating animals is mastitis, an inflammation of the mammary gland most commonly caused by bacterial infection. The severity of mastitis is greatly influenced by the invading organism and the subsequent immune response which must recognize the foreign organism, recruit immune cells, eliminate the invading pathogen, and resolve the inflammatory response. The speed, strength, and duration of this response and subsequent disease susceptibility are critically tied to the genetic background of an animal. However, the genetic contribution has been difficult to identify due to the complex interactions that must occur for effective disease resistance. Recent studies have utilized polymorphisms to better define the genes and chromosomal regions that contribute to mastitis resistance. This review will examine these studies with primary emphasis in bovine systems, as the most work regarding mastitis has been conducted in this species.

Genome Scan for Parent-of-Origin QTL Effects on Bovine Growth and Carcass Traits

Parent-of-origin effects (POE) such as genomic imprinting influence growth and body composition in livestock, rodents, and humans. Here, we report the results of a genome scan to detect quantitative trait loci (QTL) with POE on growth and carcass traits in Angus × Brahman cattle crossbreds. We identified 24 POE–QTL on 15 Bos taurus autosomes (BTAs) of which six were significant at 5% genome-wide (GW) level and 18 at the 5% chromosome-wide (CW) significance level. Six QTL were paternally expressed while 15 were maternally expressed. Three QTL influencing post-weaning growth map to the proximal end of BTA2 (linkage region of 0–9 cM; genomic region of 5.0–10.8 Mb), for which only one imprinted ortholog is known so far in the human and mouse genomes, and therefore may potentially represent a novel imprinted region. The detected QTL individually explained 1.4 ∼ 5.1% of each trait’s phenotypic variance. Comparative in silico analysis of bovine genomic locations show that 32 out of 1,442 known mammalian imprinted genes from human and mouse homologs map to the identified QTL regions. Although several of the 32 genes have been associated with quantitative traits in cattle, only two (GNAS and PEG3) have experimental proof of being imprinted in cattle. These results lend additional support to recent reports that POE on quantitative traits in mammals may be more common than previously thought, and strengthen the need to identify and experimentally validate cattle orthologs of imprinted genes so as to investigate their effects on quantitative traits.

QTLMAS 2010: simulated dataset

Background

Objective was to simulate the data for the QTLMAS 2010 Workshop under a model that includes major additive, epistatic and parent-of-origin effects.

Results

Data were simulated for 3226 individuals in 5 generations. Genomic data for 5 chromosomes were simulated using coalescent model. In total, the data included 10,031 SNPs, 30 additive QTLs, 2 interacting QTL pairs, and 3 imprinted loci. The density was 20 SNPs/1Mb, whereas mean linkage disequilibrium between adjacent SNPs was 0.1. One quantitative and one binary trait were simulated with heritability of 0.39-0.52 and additive correlation of 0.59. The data can be used as a benchmark for comparison of QTL mapping methods and models for genomic breeding value estimation under complex genetic architecture.