SNPs of MYPN and TTN genes are associated to meat and carcass traits in Italian Large White and Italian Duroc pigs

Genome-wide identification and skin expression of immunoglobulin superfamily in discus fish (Symphysodon aequifasciatus) reveal common genes associated with vertebrate lactation

Discus Symphysodon spp. employs an unusual parental care behavior where fry feed on parental skin mucus after hatching. Studies on discus immunoglobulin superfamily (IgSF) especially during parental care are scarce. Here, a total of 518 IgSF members were identified based on discus genome and clustered into 12 groups, unevenly distributing on 30 linkage groups. A total of 92 pairs of tandem duplication and 40 pairs of segmental duplication that underwent purifying selection were identified. IgSF genes expressed differentially in discus skin during different care stages and between male and female parents. Specifically, the transcription of btn1a1, similar with mammalian lactation, increased after spawning, reached a peak when fry started biting on parents' skin mucus, and then decreased. The expression of btn2a1 and other immune members, e.g., nect4, fcl5 and cd22, were up-regulated when fry stopped biting on mucus. These results suggest the expression differentiation of IgSF genes in skin of discus fish during parental care.

Genomic Uniqueness of Local Sheep Breeds From Morocco

Sheep farming is a major source of meat in Morocco and plays a key role in the country's agriculture. This study aims at characterizing the whole-genome diversity and demographic history of the main Moroccan sheep breeds, as well as to identify selection signatures within and between breeds. Whole genome data from 87 individuals representing the five predominant local breeds were used to estimate their level of neutral genetic diversity and to infer the variation of their effective population size over time. In addition, we used two methods to detect selection signatures: either for detecting selective sweeps within each breed separately or by detecting differentially selected regions by contrasting different breeds. We identified hundreds of genomic regions putatively under selection, which related to several biological terms involved in local adaptation or the expression of zootechnical performances such as Growth, UV protection, Cell maturation or Feeding behavior. The results of this study revealed selection signatures in genes that have an important role in traits of interest and increased our understanding of how genetic diversity is distributed in these local breeds. Thus, Moroccan local sheep breeds exhibit both a high genetic diversity and a large set of adaptive variations, and therefore, represent a valuable genetic resource for the conservation of sheep in the context of climate change.

Estimation of inbreeding and identification of regions under heavy selection based on runs of homozygosity in a Large White pig population

Background

Runs of homozygosity (ROHs) are homozygous segments of the genome where the two haplotypes inherited from the parents are identical. The current availability of genotypes for a very large number of single nucleotide polymorphisms (SNPs) is leading to more accurate characterization of ROHs in the whole genome. Here, we investigated the occurrence and distribution of ROHs in 3,692 Large White pigs and compared estimates of inbreeding coefficients calculated based on ROHs (F_ROH), homozygosity (F_HOM), genomic relationship matrix (F_GRM) and pedigree (F_PED). Furthermore, we identified genomic regions with high ROH frequencies and annotated their candidate genes.

Results

In total, 176,182 ROHs were identified from 3,569 animals, and all individuals displayed at least one ROH longer than 1 Mb. The ROHs identified were unevenly distributed on the autosomes. The highest and lowest coverages of Sus scrofa chromosomes (SSC) by ROH were on SSC14 and SSC13, respectively. The highest pairwise correlation among the different inbreeding coefficient estimates was 0.95 between F_{ROH_total} and F_HOM, while the lowest was − 0.083 between F_GRM and F_PED. The correlations between F_PED and F_ROH using four classes of ROH lengths ranged from 0.18 to 0.37 and increased with increasing ROH length, except for ROH > 10 Mb. Twelve ROH islands were located on four chromosomes (SSC1, 4, 6 and 14). These ROH islands harboured genes associated with reproduction, muscular development, fat deposition and adaptation, such as SIRT1, MYPN, SETDB1 and PSMD4.

Conclusion

F_ROH can be used to accurately assess individual inbreeding levels compared to other inbreeding coefficient estimators. In the absence of pedigree records, F_ROH can provide an alternative to inbreeding estimates. Our findings can be used not only to effectively increase the response to selection by appropriately managing the rate of inbreeding and minimizing the negative effects of inbreeding depression but also to help detect genomic regions with an effect on traits under selection.

Identification of loci affecting sexually dimorphic patterns for height and recurrent laryngeal neuropathy risk in American Belgian Draft Horses

Equine recurrent laryngeal neuropathy (RLN) is a bilateral mononeuropathy with an unknown etiology. In Thoroughbreds (TB), we previously demonstrated that the haplotype association for height (LCORL/NCAPG locus on ECA3, which affects body size) and RLN was coincident. In the present study, we performed a genome-wide association scan (GWAS) for RLN in 458 American Belgian Draft Horses, a breed fixed for the LCORL/NCAPG risk alelle. In this breed, RLN risk is associated with sexually dimorphic differences in height, and we identified a novel locus contributing to height in a sex-specific manner: MYPN (ECA1). Yet this specific locus contributes little to RLN risk, suggesting that other growth traits correlated to height may underlie the correlation to this disease. Controlling for height, we identified a locus on ECA15 contributing to RLN risk specifically in males. These results suggest that loci with sex-specific gene expression play an important role in altering growth traits impacting RLN etiology, but not necessarily adult height. These newly identified genes are promising targets for novel preventative and treatment strategies.

Comparative sequence alignment reveals River Buffalo genomic structural differences compared with cattle

This study sought to characterize differences in gene content, regulation and structure between taurine cattle and river buffalo (one subspecies of domestic water buffalo) using the extensively annotated UMD3.1 cattle reference genome as a basis for comparisons. We identified 127 deletion CNV regions in river buffalo representing 5 annotated cattle genes. We also characterized 583 merged mobile element insertion (MEI) events within the upstream regions of annotated cattle genes. Transcriptome analysis in various tissue types on river buffalo confirmed the absence of four cattle genes. Four genes which may be related to phenotypic differences in meat quality and color, had upstream MEI predictions and were found to have significantly elevated expression in river buffalo compared with cattle. Our comparative alignment approach and gene expression analyses suggested a functional role for many genomic structural variations, which may contribute to the unique phenotypes of river buffalo.

Runs of homozygosity reveal genome-wide autozygosity in Italian sheep breeds

The availability of dense single nucleotide polymorphism (SNP) assays allows for the determination of autozygous segments based on runs of consecutive homozygous genotypes (ROH). The aim of the present study was to investigate the occurrence and distribution of ROH in 21 Italian sheep breeds using medium-density SNP genotypes in order to characterize autozygosity and identify genomic regions that frequently appeared in ROH within individuals, namely ROH islands. After filtering, the final number of animals and SNPs retained for analyses were 502 and 46 277 respectively. A total of 12 302 ROH were identified. The mean number of ROH per breed ranged from 10.58 (Comisana) to 44.54 (Valle del Belice). The average length of ROH across breeds was 4.55 Mb and ranged from 3.85 Mb (Biellese) to 5.51 Mb (Leccese). Valle del Belice showed the highest value of inbreeding on the basis of ROH (F_ROH  = 0.099), whereas Comisana showed the lowest (F_ROH  = 0.016), and high standard deviation values revealed high variability in autozygosity levels within each breed. Differences also existed in the length of ROH. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were <10 Mb in length, with a few long ROH >25 Mb. The levels of ROH that we estimated here reflect the inbreeding history of the investigated sheep breeds. These results also highlight that ancient and recent inbreeding have had an impact on the genome of the Italian sheep breeds and suggest that several animals have experienced recent autozygosity events. Comisana and Bergamasca appeared as the less consanguineous breeds, whereas Barbaresca, Leccese and Valle del Belice showed ROH patterns typically produced by recent inbreeding. Moreover, within the genomic regions most commonly associated with ROH, several candidate genes were detected.

Integrating genome and transcriptome profiling for elucidating the mechanism of muscle growth and lipid deposition in Pekin ducks

Muscle growth and lipid deposition are co-ordinately regulated processes. Cherry Valley Pekin duck is a lean-type duck breed with high growth rate, whereas the native Pekin duck of China has high lipid deposition. Phenotypic analysis showed that native Pekin ducks have smaller fibre diameter and larger density in the breast muscle at 3 weeks of age and higher intramuscular fat content at 6 weeks of age than those in Cherry Valley Pekin ducks. We detected 17 positively selected genes (PSGs) by comparing genes mainly involved with muscle organ development, muscle contraction, peroxisome proliferator activated receptor signalling pathway, and fatty acid metabolism. In all, 52 and 206 differentially expressed genes (DEGs) were identified in transcriptomic comparisons between the two breeds at 3 and 6 weeks of age, respectively, which could potentially affect muscle growth and lipid deposition. Based on the integration of PSGs and DEGs and their functional annotations, we found that 11 and 10 genes were correlated with muscle growth and lipid deposition, respectively. Identification of candidate genes controlling quantitative traits of duck muscle might aid in elucidating the mechanisms of muscle growth and lipid deposition and could help in improving duck breeding.

Association study between single nucleotide polymorphisms in porcine genes and pork quality traits for fresh consumption and processing into Italian dry-cured ham

Single nucleotide polymorphisms (SNPs) of six genes (TTN, PRKAG3, CAST, CTSB, CTSF, and MYPN), known for associations with carcass and meat quality traits, post mortem proteolysis, were screened in a commercial crossed population of 368 heavy pigs (Large White x Landrace)×Duroc, reared according to the rules of Italian Protected Designation of Origin (PDO) dry-cured ham. Carcass, longissimus thoracis et lumborum muscle (LTL), and green ham traits were obtained after slaughtering, main weight losses of dry-cured hams were collected during processing. The results showed the impact of CAST variants on carcass weight, of CTSF on LTL tenderness, ham weight and fatness, of PRKAG3 and TTN on ultimate pH, hamweight. This study, while confirming significant associations between SNPs of genes and qualitative traits of carcass, longissimus and ham, supports CTSF as candidate gene suitable for fresh consumption purpose (tenderness of longissimus at 24h post mortem), and for dry-cured ham processing (higher thickness of ham subcutaneous fat).

Harnessing cross-species alignment to discover SNPs and generate a draft genome sequence of a bighorn sheep (Ovis canadensis)

Background

Whole genome sequences (WGS) have proliferated as sequencing technology continues to improve and costs decline. While many WGS of model or domestic organisms have been produced, a growing number of non-model species are also being sequenced. In the absence of a reference, construction of a genome sequence necessitates de novo assembly which may be beyond the ability of many labs due to the large volumes of raw sequence data and extensive bioinformatics required. In contrast, the presence of a reference WGS allows for alignment which is more tractable than assembly. Recent work has highlighted that the reference need not come from the same species, potentially enabling a wide array of species WGS to be constructed using cross-species alignment. Here we report on the creation a draft WGS from a single bighorn sheep (Ovis canadensis) using alignment to the closely related domestic sheep (Ovis aries).

Results

Two sequencing libraries on SOLiD platforms yielded over 865 million reads, and combined alignment to the domestic sheep reference resulted in a nearly complete sequence (95% coverage of the reference) at an average of 12x read depth (104 SD). From this we discovered over 15 million variants and annotated them relative to the domestic sheep reference. We then conducted an enrichment analysis of those SNPs showing fixed differences between the reference and sequenced individual and found significant differences in a number of gene ontology (GO) terms, including those associated with reproduction, muscle properties, and bone deposition.

Conclusion

Our results demonstrate that cross-species alignment enables the creation of novel WGS for non-model organisms. The bighorn sheep WGS will provide a resource for future resequencing studies or comparative genomics.

Electronic supplementary material

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