Fine mapping of a swine quantitative trait locus for number of vertebrae and analysis of an orphan nuclear receptor, germ cell nuclear factor (NR6A1)

Integrating Significant SNPs Identified by GWAS for Genomic Prediction of the Number of Ribs and Carcass Length in Suhuai Pigs

The number of ribs (NRs) and the carcass length (CL) are important economic traits. The traits are usually measured after slaughter. To improve the prediction performance of genomic selection (GS) for NRs and CL, one strategy is to integrate the significant loci identified from whole-genome sequencing (WGS) data by genome-wide association study (GWAS) into the genomic prediction (GP) model. This study investigated the GP of different genomic best linear unbiased prediction (GBLUP) and Bayesian models using chip genotype data, imputed WGS (iWGS) data and modeling significant single-nucleotide polymorphisms (SNPs) in different ways for the GP of NRs and CL in the Suhuai pig population. The prediction accuracy, bias and running time of 15 different GP models were evaluated by 10-fold cross-validation. The prediction accuracy of GBLUP using chip data for NRs and CL was 0.314 ± 0.022 and 0.194 ± 0.040, respectively. For NRs, based on the iWGS data, treating the most significant SNP as fixed effects in the GBLUP model had the highest predictive performance, with a prediction accuracy of 0.528 ± 0.023. For CL, based on the chip data, the model that added all the significant SNPs identified by imputed data by GWAS into the multi-trait GBLUP as the second random additive effect was the highest predictive performance, with a prediction accuracy of 0.305 ± 0.027. This study provides insights into optimizing GP models for small populations with phenotypes that are difficult to measure.

Understanding molecular mechanisms of vertebral number of variations on Mongolian sheep using candidate genes analysis

OBJECTIVE

This study aimed to investigate the genetic link between variations in vertebral number and meat production traits, such as body weight and body measurements (body length, body height, heart girth, and shin width) in Mongolian (Bayantsagaan) sheep. Additionally, we examined the association of single-nucleotide polymorphisms (SNPs) in candidate genes, particularly vertnin (VRTN), nuclear receptor subfamily 6, group A, member 1 (NR6A1), and synapse differentiation-inducing 1-like (SYNDIG1L), with vertebral number variations and their potential impact on meat production traits.

METHODS

The study involved 220 Bayantsagaan sheep from Bayantsagaan soum, Tov province, Mongolia, including 104 sheep with extra vertebrae group and 116 individuals with typical vertebral number as the control group. Morphological data, including body weight and body measurements, were collected, and genetic samples were obtained. The impact of vertebral number on morphological traits was estimated using a general linear model. The SNPs in the VRTN, NR6A1, and SYNDIG1L genes were sequenced, and their association with vertebral number was analyzed using one-way analysis of variance.

RESULTS

Bayantsagaan sheep with extra vertebrae were, on average, 4.45 kg heavier and exhibited higher variability in body size traits compared to the control group. Four polymorphic sites were identified at the VRTN gene, with one polymorphic locus (VRTN1716) showing a significant association with vertebrae number and body size. Sheep with C/C genotype at VRTN1716 locus, had more vertebrae and larger body size compared to other genotypes.

CONCLUSION

The findings suggest that variations in vertebral number and VRTN gene polymorphisms are linked to favorable meat production traits in Bayantsagaan sheep. The identified SNP (VRTN1716) associated with vertebral number and body size offers the potential for marker-assisted selection in breeding programs. These results provide valuable insights into the genetic basis of meat production traits in Bayantsagaan sheep and may contribute to the development of more efficient breeding strategies.

Bioinformatic approach to identifying causative missense polymorphisms in animal genomes

BACKGROUND

Trends in the development of genetic markers for the purposes of genomic and marker-assisted selection primarily focus on identifying causative polymorphisms. Using these polymorphisms as markers enables a more accurate association between genotype and phenotype. Bioinformatic analysis allows calculating the impact of missense polymorphisms on the structural and functional characteristics of proteins, which makes it promising for identifying causative polymorphisms. In this study, a bioinformatic approach is applied to evaluate and differentiate polymorphisms based on their causality in genes that affect the production traits of pigs and cows, which are two important livestock species.

RESULTS

The influence of both known causative and candidate missense polymorphisms in the MC4R, NR6A1, PRKAG3, RYR1, and SYNGR2 genes of pigs, as well as the ABCG2, DGAT1, GHR, and MSTN genes of cows, was assessed. The study included an evaluation of the effect of polymorphisms on protein functions, considering the evolutionary and physicochemical characteristics of amino acids at polymorphic sites. Additionally, it examined the impact of polymorphisms on the stability of tertiary protein structures, including changes in folding, binding of protein monomers, and interaction with ligands.

CONCLUSIONS

The comprehensive bioinformatic analysis used in this study enables the differentiation of polymorphisms into neutral, where both amino acids in the polymorphic site do not significantly affect the structure and function of the protein, and causative, where one of the amino acids significantly impacts the protein's properties. This approach can be employed in future research to screen extensive sets of polymorphisms in animal genomes, identifying the most promising polymorphisms for further investigation in association studies.

Variants in NR6A1 cause a novel oculo-vertebral-renal (OVR) syndrome

Colobomatous microphthalmia is a potentially blinding congenital ocular malformation that can present either in isolation or together with other syndromic features. Despite a strong genetic component to disease, many cases lack a molecular diagnosis. We describe a novel autosomal dominant oculo-vertebral-renal (OVR) syndrome in six independent families characterized by colobomatous microphthalmia, missing vertebrae and congenital kidney abnormalities. Genome sequencing identified six rare variants in the orphan nuclear receptor gene NR6A1 in these families. We performed in silico, cellular and zebrafish experiments to demonstrate the NR6A1 variants were pathogenic or likely pathogenic for OVR syndrome. Knockdown of either or both zebrafish paralogs of NR6A1 results in abnormal eye and somite development, which was rescued by wild-type but not variant NR6A1 mRNA. Illustrating the power of genomic ascertainment in medicine, our study establishes NR6A1 as a critical factor in eye and vertebral development and a pleiotropic gene responsible for OVR syndrome.

Variants in NR6A1 cause a novel oculo-vertebral-renal (OVR) syndrome

Colobomatous microphthalmia is a potentially blinding congenital ocular malformation that can present either in isolation or together with other syndromic features. Despite a strong genetic component to disease, many cases lack a molecular diagnosis. We describe a novel autosomal dominant oculo-vertebral-renal (OVR) syndrome in six independent families characterized by colobomatous microphthalmia, missing vertebrae and congenital kidney abnormalities. Genome sequencing identified six rare variants in the orphan nuclear receptor gene NR6A1 in these families. We performed in silico, cellular and zebrafish experiments to demonstrate the NR6A1 variants were pathogenic or likely pathogenic for OVR syndrome. Knockdown of either or both zebrafish paralogs of NR6A1 results in abnormal eye and somite development, which was rescued by wild-type but not variant NR6A1 mRNA. Illustrating the power of genomic ascertainment in medicine, our study establishes NR6A1 as a critical factor in eye and vertebral development and a pleiotropic gene responsible for OVR syndrome.

Transcriptional Signatures of Domestication Revealed through Meta-Analysis of Pig, Chicken, Wild Boar, and Red Junglefowl Gene Expression Data

Domesticated animals have undergone significant changes in their behavior, morphology, and physiological functions during domestication. To identify the changes in gene expression associated with domestication, we collected the RNA-seq data of pigs, chickens, wild boars, and red junglefowl from public databases and performed a meta-analysis. Gene expression was quantified, and the expression ratio between domesticated animals and their wild ancestors (DW-ratio) was calculated. Genes were classified as "upregulated", "downregulated", or "unchanged" based on their DW-ratio, and the DW-score was calculated for each gene. Gene set enrichment analysis revealed that genes upregulated in pigs were related to defense from viral infection, whereas those upregulated in chickens were associated with aminoglycan and carbohydrate derivative catabolic processes. Genes commonly upregulated in pigs and chickens are involved in the immune response, olfactory learning, epigenetic regulation, cell division, and extracellular matrix. In contrast, genes upregulated in wild boar and red junglefowl are related to stress response, cell proliferation, cardiovascular function, neural regulation, and energy metabolism. These findings provide valuable insights into the genetic basis of the domestication process and highlight potential candidate genes for breeding applications.

Copy number variants selected during pig domestication inferred from whole genome resequencing

Over extended periods of natural and artificial selection, China has developed numerous exceptional pig breeds. Deciphering the germplasm characteristics of these breeds is crucial for their preservation and utilization. While many studies have employed single nucleotide polymorphism (SNP) analysis to investigate the local pig germplasm characteristics, copy number variation (CNV), another significant type of genetic variation, has been less explored in understanding pig resources. In this study, we examined the CNVs of 18 Wanbei pigs (WBP) using whole genome resequencing data with an average depth of 12.61. We identified a total of 8,783 CNVs (~30.07 Mb, 1.20% of the pig genome) in WBP, including 8,427 deletions and 356 duplications. Utilizing fixation index (Fst), we determined that 164 CNVs were within the top 1% of the Fst value and defined as under selection. Functional enrichment analyses of the genes associated with these selected CNVs revealed genes linked to reproduction (SPATA6, CFAP43, CFTR, BPTF), growth and development (NR6A1, SMYD3, VIPR2), and immunity (PARD3, FYB2). This study enhances our understanding of the genomic characteristics of the Wanbei pig and offers a theoretical foundation for the future breeding of this breed.

Unlocking the Transcriptional Control of NCAPG in Bovine Myoblasts: CREB1 and MYOD1 as Key Players

Muscle formation directly determines meat production and quality. The non-SMC condensin I complex subunit G (NCAPG) is strongly linked to the growth features of domestic animals because it is essential in controlling muscle growth and development. This study aims to elucidate the tissue expression level of the bovine NCAPG gene, and determine the key transcription factors for regulating the bovine NCAPG gene. In this study, we observed that the bovine NCAPG gene exhibited high expression levels in longissimus dorsi and spleen tissues. Subsequently, we cloned and characterized the promoter region of the bovine NCAPG gene, consisting of a 2039 bp sequence, through constructing the deletion fragment double-luciferase reporter vector and site-directed mutation-identifying core promoter region with its key transcription factor binding site. In addition, the key transcription factors of the core promoter sequence of the bovine NCAPG gene were analyzed and predicted using online software. Furthermore, by integrating overexpression experiments and the electrophoretic mobility shift assay (EMSA), we have shown that cAMP response element binding protein 1 (CREB1) and myogenic differentiation 1 (MYOD1) bind to the core promoter region (-598/+87), activating transcription activity in the bovine NCAPG gene. In conclusion, these findings shed important light on the regulatory network mechanism that underlies the expression of the NCAPG gene throughout the development of the muscles in beef cattle.

The expanding roles of Nr6a1 in development and evolution

The Nuclear Receptor (NR) family of transcriptional regulators possess the ability to sense signalling molecules and directly couple that to a transcriptional response. While this large class of proteins are united by sequence and structural homology, individual NR functional output varies greatly depending on their expression, ligand selectivity and DNA binding sequence specificity. Many NRs have remained somewhat enigmatic, with the absence of a defined ligand categorising them as orphan nuclear receptors. One example is Nuclear Receptor subfamily 6 group A member 1 (Nr6a1), an orphan nuclear receptor that has no close evolutionary homologs and thus is alone in subfamily 6. Nonetheless, Nr6a1 has emerged as an important player in the regulation of key pluripotency and developmental genes, as functionally critical for mid-gestational developmental progression and as a possible molecular target for driving evolutionary change in animal body plan. Here, we review the current knowledge on this enigmatic nuclear receptor and how it impacts development and evolution.

Advancements in Genetic Marker Exploration for Livestock Vertebral Traits with a Focus on China

In livestock breeding, the number of vertebrae has gained significant attention due to its impact on carcass quality and quantity. Variations in vertebral traits have been observed across different animal species and breeds, with a strong correlation to growth and meat production. Furthermore, vertebral traits are classified as quantitative characteristics. Molecular marker techniques, such as marker-assisted selection (MAS), have emerged as efficient tools to identify genetic markers associated with vertebral traits. In the current review, we highlight some key potential genes and their polymorphisms that play pivotal roles in controlling vertebral traits (development, length, and number) in various livestock species, including pigs, donkeys, and sheep. Specific genetic variants within these genes have been linked to vertebral development, number, and length, offering valuable insights into the genetic mechanisms governing vertebral traits. This knowledge has significant implications for selective breeding strategies to enhance structural characteristics and meat quantity and quality in livestock, ultimately improving the efficiency and quality of the animal husbandry industry.

Candidate SNPs for meat quality and carcass composition in free-range Iberian pigs

Several genetic markers, previously associated with meat quality traits, have been proposed to be included in Iberian pig breeding programs. However, before being implemented, effects of these candidate SNPs on premium cuts' yield should be evaluated to avoid potential undesirable antagonistic effects. Therefore, the main goals of this study were to evaluate the effects of a set of 26 polymorphisms on premium cuts weights and to corroborate their effects on meat quality in a larger population. Phenotypic data of approximately 1550 Iberian pigs were recorded. The PRKAG3_rs319678464C and PRKAG3_rs1108399077G alleles were associated with an increase of shear force, water losses and color values and a reduction of average daily gain (ADG). The CAST_rs196949783G > A and ADIPOQ_rs3476515794T > G SNPs mainly affected IMF content, with this last SNP being the only one with significant effects on any of the carcass measures, specifically Longissimus thoracis et lumborum (LTL) weight. The ELOVL6_rs3473714672A, FASN_rs331694510A, MTTP_rs335896411C and ACACA_rs340781986C alleles were linked to a higher percentage of oleic acid and monounsaturated FA and a decrease in palmitic, palmitoleic and saturated FA. Besides, suggestive effects were observed between ELOVL6_rs3473714672A > G and ham and shoulder weights, and between MTTP_rs335896411T > C and LTL muscle weight. Finally, the NR6A1_rs326780270T allele was associated with a significant increase in ADG and tended to reduce the ham weight. The SNPs mapped to PRKAG3 (rs319678464G > C), ACACA, FASN and CAST genes could be implemented in breeding programs to improve meat quality traits without undesirable effects on carcass composition. The SNPs mapped to PRKAG3 (rs1108399077G > A) ELOVL6, MTTP and NR6A1 should continue to be tested in a larger number of animals.

The control of transitions along the main body axis

Although vertebrates display a large variety of forms and sizes, the mechanisms controlling the layout of the basic body plan are substantially conserved throughout the clade. Following gastrulation, head, trunk, and tail are sequentially generated through the continuous addition of tissue at the caudal embryonic end. Development of each of these major embryonic regions is regulated by a distinct genetic network. The transitions from head-to-trunk and from trunk-to-tail development thus involve major changes in regulatory mechanisms, requiring proper coordination to guarantee smooth progression of embryonic development. In this review, we will discuss the key cellular and embryological events associated with those transitions giving particular attention to their regulation, aiming to provide a cohesive outlook of this important component of vertebrate development.

Signatures of Admixture and Genetic Uniqueness in the Autochthonous Greek Black Pig Breed Deduced from Gene Polymorphisms Affecting Domestication-Derived Traits

The Greek Black Pig (or Greek Pig) is the only recognized autochthonous pig breed raised in Greece, usually in extensive or semi-extensive production systems. According to its name, the characteristic breed coat color is solid black. In this study, with the aim to start a systematic genetic characterization of the Greek Black Pig breed, we investigated polymorphisms in major genes well known to affect exterior and production traits (MC1R, KIT, NR6A1, VRTN and IGF2) and compared these data with population genetic information available in other Mediterranean and Western Balkan pig breeds and wild boars. None of the investigated gene markers were fixed for one allele, suggesting that, in the past, this breed experienced introgression from wild boars and admixture from cosmopolitan pig breeds, enriching the breed genetic pool that should be further investigated to design appropriate conservation genetic strategies. We identified a new MC1R allele, containing two missense mutations already reported in two other independent alleles, but here present in the same haplotype. This allele might be useful to disclose biological information that can lead to better understanding the cascade transmission of signals to produce melanin pigments. This study demonstrated that autochthonous genetic resources can be an interesting reservoir of unexpected genetic variants.

Integrated Analysis of Transcriptome Expression Profiles Reveals miRNA-326-NKX3.2-Regulated Porcine Chondrocyte Differentiation

The porcine body length trait is an essential factor affecting meat production and reproductive performance. It is evident that the development/lengthening of individual vertebrae is one of the main reasons for increases in body length; however, the underlying molecular mechanism remains unclear. In this study, RNA-seq analysis was used to profile the transcriptome (lncRNA, mRNA, and miRNA) of the thoracic intervertebral cartilage (TIC) at two time points (1 and 4 months) during vertebral column development in Yorkshire (Y) and Wuzhishan pigs (W). There were four groups: 1- (Y1) and 4-month-old (Y4) Yorkshire pigs and 1- (W1) and 4-month-old (W4) Wuzhishan pigs. In total, 161, 275, 86, and 126 differentially expressed (DE) lncRNAs, 1478, 2643, 404, and 750 DE genes (DEGs), and 74,51, 34, and 23 DE miRNAs (DE miRNAs) were identified in the Y4 vs. Y1, W4 vs. W1, Y4 vs. W4, and Y1 vs. W1 comparisons, respectively. Functional analysis of these DE transcripts (DETs) demonstrated that they had participated in various biological processes, such as cellular component organization or biogenesis, the developmental process, the metabolic process, bone development, and cartilage development. The crucial bone development-related candidate genes NK3 Homeobox 2 (NKX3.2), Wnt ligand secretion mediator (WLS), gremlin 1 (GREM1), fibroblast growth factor receptor 3 (FGFR3), hematopoietically expressed homeobox (HHEX), (collagen type XI alpha 1 chain (COL11A1), and Wnt Family Member 16 (WNT16)) were further identified by functional analysis. Moreover, lncRNA, miRNA, and gene interaction networks were constructed; a total of 55 lncRNAs, 6 miRNAs, and 7 genes formed lncRNA-gene, miRNA-gene, and lncRNA-miRNA-gene pairs, respectively. The aim was to demonstrate that coding and non-coding genes may co-regulate porcine spine development through interaction networks. NKX3.2 was identified as being specifically expressed in cartilage tissues, and it delayed chondrocyte differentiation. miRNA-326 regulated chondrocyte differentiation by targeting NKX3.2. The present study provides the first non-coding RNA and gene expression profiles in the porcine TIC, constructs the lncRNA-miRNA-gene interaction networks, and confirms the function of NKX3.2 in vertebral column development. These findings contribute to the understanding of the potential molecular mechanisms regulating pig vertebral column development. They expand our knowledge about the differences in body length between different pig species and provide a foundation for future studies.

Identification of breed-specific genomic variants in Colombian Creole pig breeds by whole-genome sequencing

Dissecting genetic variation of local breeds is important for the success of conservation. In this research, we investigated the genomic variation of Colombian Creole (CR) pigs, with a focus on the breed-specific variants in the exonic region of 34 genes with reported effects on adaptive and economic traits. Seven individuals of each of the three CR breeds (CM, Casco de Mula; SP, San Pedreño; and ZU, Zungo) were whole-genome sequenced along with 7 Iberian (IB) pigs and 7 pigs of each of the four most used cosmopolitan (CP) breeds (Duroc, Landrace × Large White, and Pietrain). Molecular variability in CR (6,451,218 variants; from 3,919,242, in SP, to 4,648,069, in CM) was comparable to that in CP, but higher than in IB. For the investigated genes, SP pigs displayed less exonic variants (178) than ZU (254), CM (263), IB (200), and the individual CP genetic types (201 to 335). Sequence variation in these genes confirmed the resemblance of CR to IB and indicates that CR pigs, particularly ZU and CM, are not exempt from selective introgression of other breeds. A total of 50 exonic variants were identified as being potentially specific to CR, including a high-impact deletion in the intron between exons 15 and 16 of the leptin receptor gene, which was only found in CM and ZU. The identification of breed-specific variants in genes related to adaptive and economical traits can bolster the understanding of the role of gene-environment interactions on local adaptation and points the way for effective breeding and conservation of CR pigs.

Integration of Selection Signatures and Protein Interactions Reveals NR6A1, PAPPA2, and PIK3C2B as the Promising Candidate Genes Underlying the Characteristics of Licha Black Pig

Licha black (LI) pig has the specific characteristics of larger body length and appropriate fat deposition among Chinese indigenous pigs. Body length is one of the external traits that affect production performance, and fat deposition influences meat quality. However, the genetic characteristics of LI pigs have not yet been systematically uncovered. Here, the genomic information from 891 individuals of LI pigs, commercial pigs, and other Chinese indigenous pigs was used to analyze the breed characteristics of the LI pig with runs of homozygosity, haplotype, and FST selection signatures. The results showed the growth traits-related genes (i.e., NR6A1 and PAPPA2) and the fatness traits-related gene (i.e., PIK3C2B) were the promising candidate genes that closely related to the characteristics of LI pigs. In addition, the protein–protein interaction network revealed the potential interactions between the promising candidate genes and the FASN gene. The RNA expression data from FarmGTEx indicated that the RNA expression levels of NR6A1, PAPPA2, PIK3C2B, and FASN were highly correlated in the ileum. This study provides valuable molecular insights into the mechanisms that affect pig body length and fat deposition, which can be used in the further breeding process to improve meat quality and commercial profitability.

Genome-wide association study reveals new QTL and functional candidate genes for the number of ribs and carcass length in pigs

The number of ribs (NR) and carcass length (CL) are important economic traits in pig breeding programs. Pigs with a higher NR and longer CL produce greater pork yields. In the present study, Suhuai pigs with NR and CL phenotypes were genotyped using the Neogen^® GGP Porcine 80 K SNP array to identify the QTL affecting NR and CL and dissect the candidate genes for the two traits. The SNP-chip data was imputed to the whole-genome sequence (iWGS) to increase the probability of identifying causal variants. Through genome-wide association studies (GWAS) based on both chip and iWGS data, significant SNPs were detected on Sus scrofa chromosome (SSC) 1, SSC4 and SSC7 for NR and on SSC5, SSC16 and SSC17 for CL. Moreover, two SNPs (H3GA0022644 and WU_10.2_7_103460706) on SSC7 detected in chip-based GWAS were significantly associated with both NR and CL. Through Bayes fine mapping, one reported QTL for NR on SSC7 and two reported QTL for CL on SSC17 were verified, and two new QTL (SSC1: 14.05-15.84 Mb and SSC4: 64.83-66.59 Mb) affecting NR and two new QTL (SSC5: 58.31-59.84 Mb and SSC16: 22.98-23.43 Mb) affecting CL were detected. According to the biological functions of genes, MTHFD1L on SSC1 and SULF1 on SSC4 are novel functional candidate genes for NR, and EMP1 on SSC5 and EGFLAM on SSC16 are novel functional candidate genes for CL. Overall, our findings provide a basis for identifying new causal genes and mutations affecting NR and CL.

Transcriptomics and Selection Pressure Analysis Reveals the Influence Mechanism of PLIN1 Protein on the Development of Small Size in Min Pigs

Body size is an important biological phenotypic trait that has attracted substantial attention. Small domestic pigs can serve as excellent animal models for biomedicine and also help meet sacrificial culture needs in human societies. Although the mechanisms underlying vertebral development regulating body size variation in domestic pigs during the embryonic period have been well described, few studies have examined the genetic basis of body size variation in post embryonic developmental stages. In this study, seven candidate genes-PLIN1, LIPE, PNPLA1, SCD, FABP5, KRT10 and IVL-significantly associated with body size were identified in Min pigs, on the basis of weighted gene co-expression network analysis (WGCNA), and most of their functions were found to be associated with lipid deposition. Six candidate genes except for IVL were found to have been subjected to purifying selection. PLIN1 had the lowest ω value (0.139) and showed heterogeneous selective pressure among domestic pig lineages with different body sizes (p < 0.05). These results suggested that PLIN1 is an important genetic factor regulating lipid deposition and consequently affecting body size variation in pigs. The culture of whole pig sacrifice in Manchu during the Qing Dynasty in China might have contributed to the strong artificial domestication and selection of Hebao pigs.

The selected genes NR6A1, RSAD2-CMPK2, and COL3A1 contribute to body size variation in Meishan pigs through different patterns

The high-fertility Meishan pig is currently categorized into medium sized (MMS) and small sized (SMS) based on body size. To identify causal genes responsible for the variation in body size within the two categories, we sequenced individuals representing the entire consanguinity of the existing Meishan pig. This enabled us to conduct genome selective signal analysis. Our findings revealed the genomes of MMS and SMS are stratified, with selective sweep regions formed by differential genomic intervals between the two categories enriched in multiple pig body size related quantitative trait loci (QTLs). Furthermore, the missense mutation c.575T > C of candidate causal gene NR6A1, accounting for the variation in lumbar vertebrae number in pigs, was positively selected in MMS only, leading to an increase in body length of MMS at 6 months of age. To precisely identify causal genes accounting for body size variation through multi-omics, we collected femoral cartilage and liver transcription data from MMS and SMS respectively, and re-sequencing data from pig breeds exhibiting varying body sizes. We found that two selected regions where the RSAD2-CMPK2 and COL3A1 genes are located, respectively, showed different haplotypes in pig breeds of varying body size, and was associated with body or carcass length in hybridized Suhuai pig. Additionally, the above three hub genes, were significantly greater expressed in SMS femoral cartilage and liver tissues compared to MMS. These three genes could strengthen the pathways related to bone resorption and metabolism in SMS, potentially hindering bone and skeletal development and resulting in a smaller body size in SMS. These findings provide valuable insights into the genetic mechanism of body size variation in Meishan pig population.

Nr6a1 controls Hox expression dynamics and is a master regulator of vertebrate trunk development

The vertebrate main-body axis is laid down during embryonic stages in an anterior-to-posterior (head-to-tail) direction, driven and supplied by posteriorly located progenitors. Whilst posterior expansion and segmentation appears broadly uniform along the axis, there is developmental and evolutionary support for at least two discrete modules controlling processes within different axial regions: a trunk and a tail module. Here, we identify Nuclear receptor subfamily 6 group A member 1 (Nr6a1) as a master regulator of trunk development in the mouse. Specifically, Nr6a1 was found to control vertebral number and segmentation of the trunk region, autonomously from other axial regions. Moreover, Nr6a1 was essential for the timely progression of Hox signatures, and neural versus mesodermal cell fate choice, within axial progenitors. Collectively, Nr6a1 has an axially-restricted role in all major cellular and tissue-level events required for vertebral column formation, supporting the view that changes in Nr6a1 levels may underlie evolutionary changes in axial formulae.

Identification of LTBP2 gene polymorphisms and their association with thoracolumbar vertebrae number, body size, and carcass traits in Dezhou donkeys

The number of thoracolumbar vertebrae in Dezhou donkeys varies from 22 to 24 and is associated with body size and carcass traits. In mammals, the latent transforming growth factor beta binding protein 2 (LTBP2) has been found to have some functions in the development of thoracolumbar vertebrae. The relationship between LTBP2 and TLN (the number of thoracolumbar vertebrae) of Dezhou donkeys is yet to be reported. The purposes of this study are as follows: 1) to quantify the effect of thoracolumbar vertebrae number variation of Dezhou donkeys on body size and carcass trait; 2) to study the distribution of single nucleotide variants (SNVs) in the LTBP2 gene of Dezhou donkeys; and 3) to explore whether these SNVs can be used as candidate sites to study the mechanism of Dezhou donkey muti-thoracolumbar vertebrae development. The TLN, body size, and carcass traits of 392 individuals from a Dezhou donkey breed were recorded. All animals were sequenced for LTBP2 using GBTS liquid chip and 16 SNVs were used for further analysis. We then analyzed the relationship between these SNVs with TLN, body size, and carcass traits. The results showed that: 1) c.5547 + 860 C > T, c.5251 + 281 A > C, c.3769 + 40 C > T, and c.2782 + 3975 A > G were complete genetic linkages and significantly associated with thoracic vertebrae number (TN) (p < 0.05) (wild-type homozygotes had more TN than heterozygotes); 2) c.1381 + 768 T > G and c.1381 + 763 G > T were significantly associated with lumber vertebrae number (LN) (p < 0.05); 3) c.1003 + 704 C > T, c.1003 + 651 C > T, c.1003 + 626 A > G, and c.812 + 22526 T > G were significantly associated with chest circumference (CHC), front carcass weight (CWF), after carcass weight (CWA), and carcass weight (CW) (p < 0.05) (wild-type homozygotes were larger than other genotypes in CHC, CWF, CWA, and CW); and 4) the effect of variation is not consistent in c.565 + 11921 A > G, c.565 + 6840 A > G, c.565 + 3453 C > T, and c.494 + 5808 C > T. These results provide useful information that the polymorphism of LTBP2 is significantly associated with TLN, body size, and carcass traits in Dezhou donkeys, which can serve as a molecule marker to improve donkey production performance.

Establishment of a Prognostic Model for Hepatocellular Carcinoma Based on Bioinformatics and the Role of NR6A1 in the Progression of HCC

BACKGROUND AND AIMS

Generally acceptable prognostic models for hepatocellular carcinoma (HCC) are not available. This study aimed to establish a prognostic model for HCC by identifying immune-related differentially expressed genes (IR-DEGs) and to investigate the potential role of NR6A1 in the progression of HCC.

METHODS

Bioinformatics analysis using The Cancer Genome Atlas and ImmPort databases was used to identify IR-DEGs. Lasso Cox regression and multivariate Cox regression analysis were used to establish a prognostic model of HCC. Kaplan-Meier analysis and the receiver operating characteristic (ROC) curves were used to evaluate the performance of the prognostic model, which was further verified in the International Cancer Genome Consortium (ICGC) database. Gene set enrichment analysis was used to explore the potential pathways of NR6A1. Cell counting kit 8, colony formation, wound healing, and Transwell migration assays using Huh7 cells, and tumor formation models in nude mice were conducted.

RESULTS

A prognostic model established based on ten identified IR-DEGs including HSPA4, FABP6, MAPT, NDRG1, APLN, IL17D, LHB, SPP1, GLP1R, and NR6A1, effectively predicted the prognosis of HCC patients, was confirmed by the ROC curves and verified in ICGC database. NR6A1 expression was significantly up-regulated in HCC patients, and NR6A1 was significantly associated with a low survival rate. Gene set enrichment analysis showed the enrichment of cell cycle, mTOR, WNT, and ERBB signaling pathways in patients with high NR6A1 expression. NR6A1 promoted cell proliferation, invasiveness, migration, and malignant tumor formation and growth in vitro and in vivo.

CONCLUSIONS

An effective prognostic model for HCC, based on a novel signature of 10 immune-related genes, was established. NR6A1 was up-regulated in HCC and was associated with a poor prognosis of HCC. NR6A1 promoted cell proliferation, migration, and growth of HCC, most likely through the cell cycle, mTOR, WNT, and ERBB signaling pathways.

Current Analytical Methods and Research Trends Are Used to Identify Domestic Pig and Wild Boar DNA in Meat and Meat Products

The pig, one of the most important livestock species, is a meaningful source of global meat production. It is necessary, however, to prove whether a food product that a discerning customer selects in a store is actually made from pork or venison, or does not contain it at all. The problem of food authenticity is widespread worldwide, and cases of meat adulteration have accelerated the development of food and the identification methods of feed species. It is worth noting that several different molecular biology techniques can identify a porcine component. However, the precise differentiation between wild boar and a domestic pig in meat products is still challenging. This paper presents the current state of knowledge concerning the species identification of the domestic pig and wild boar DNA in meat and its products.

Natural adaptation and human selection of northeast African sheep genomes

African sheep manifest diverse but distinct physio-anatomical traits, which are the outcomes of natural- and human-driven selection. Here, we generated 34.8 million variants from 150 indigenous northeast African sheep genomes sequenced at an average depth of ∼54× for 130 samples (Ethiopia, Libya) and ∼20× for 20 samples (Sudan). These represented sheep from diverse environments, tail morphology and post-Neolithic introductions to Africa. Phylogenetic and model-based admixture analysis provided evidence of four genetic groups corresponding to altitudinal geographic origins, tail morphotypes and possible historical introduction and dispersal of the species into and across the continent. Running admixture at higher levels of K (6 ≤ K ≤ 25), revealed cryptic levels of genome intermixing as well as distinct genetic backgrounds in some populations. Comparative genomic analysis identified targets of selection that spanned conserved haplotype structures overlapping clusters of genes and gene families. These were related to hypoxia responses, ear morphology, caudal vertebrae and tail skeleton length, and tail fat-depot structures. Our findings provide novel insights underpinning morphological variation and response to human-driven selection and environmental adaptation in African indigenous sheep.

Multi-Thoracolumbar Variations and NR6A1 Gene Polymorphisms Potentially Associated with Body Size and Carcass Traits of Dezhou Donkey

The number of thoracolumbar vertebrae is a quantitative trait positively correlated with the economic traits of livestock. More thoracolumbar vertebrae individuals could genetically be used to improve the livestock population, as more thoracolumbar vertebrae means a longer carcass, which could bring more meat production. Nuclear receptor subfamily 6 group A member 1 (NR6A1) is considered a strong candidate gene for effecting the number of vertebrae in livestock. The purposes of this study are as follows: (a) Analyzing the effect of TLN variation on body size and carcass traits of Dezhou donkey; (b) Studying the distribution of seven single nucleotide variants (SNVs) in NR6A1 gene of Dezhou donkey; (c) Exploring the relationship between latent SNVs and TLN, the body size and carcass traits. We examined the thoracic and lumbar vertebrae number and seven SNVs in NR6A1 gene of 455 Dezhou donkeys, and analyzed the relationships between them. Five types of thoracolumbar combinations (T17L5 (individual with 17 thoracic and five lumbar vertebrae) 2.4%, T18L5 75.8%, T19L5 1.1%, T17L6 11.9%, and T18L6 8.8%) of Dezhou donkeys were detected in this study. For one thoracolumbar vertebra added, the body length of Dezhou donkey increases by 3 cm and the carcass weight increases by 6 kg. Seven SNVs (g.18093100G > T, g.18094587G > T, g.18106043G > T, g.18108764G > T, g.18110615T > G, g.18112000C > T and g.18114954T > G) of the NR6A1 gene were found to have a significant association with the TLN, body size and carcass traits of Dezhou donkey (p < 0.05), respectively. For instance, g.18114954C > T is significantly associated with lumber vertebrae number, the total number of thoracolumbar, and carcass weight, and individuals with TT genotype had significantly larger value than CC genotype (p < 0.05). Using these 7SNVs, 16 different haplotypes were estimated. Compared to Hap3Hap3, individuals homozygous for Hap2Hap2 showed significantly longer length in one thoracic spine (STL), the total thoracic vertebrae and one thoracolumbar spine. Our study will not only extend the understanding of genetic variation in the NR6A1 gene of Dezhou donkey, but also provide useful information for marker assisted selection in donkey breeding program.

Regulation of posterior Hox genes by sex steroids explains vertebral variation in inbred mouse strains

A series of elegant embryo transfer experiments in the 1950s demonstrated that the uterine environment could alter vertebral patterning in inbred mouse strains. In the intervening decades, attention has tended to focus on the technical achievements involved and neglected the underlying biological question: how can genetically homogenous individuals have a heterogenous number of vertebrae? Here I revisit these experiments and, with the benefit of knowledge of the molecular-level processes of vertebral patterning gained over the intervening decades, suggest a novel hypothesis for homeotic transformation of the last lumbar vertebra to the adjacent sacral type through regulation of Hox genes by sex steroids. Hox genes are involved in both axial patterning and development of male and female reproductive systems and have been shown to be sensitive to sex steroids in vitro and in vivo. Regulation of these genes by sex steroids and resulting alterations to vertebral patterning may hint at a deep evolutionary link between the ribless lumbar region of mammals and the switch from egg-laying to embryo implantation. An appreciation of the impact of sex steroids on Hox genes may explain some puzzling aspects of human disease, and highlights the spine as a neglected target for in utero exposure to endocrine disruptors.

Whole-genome resequencing to investigate the determinants of the multi-lumbar vertebrae trait in sheep

Multi-lumbar vertebrae trait is a beneficial mutation that can significantly improve livestock meat production. However, the genetic basis of the multi-lumbar vertebrae in sheep is still unclear. Here, we analysed the number of lumbar vertebrae of Duolang sheep and found three different traits of lumbar vertebrae number. Compared with the normal sheep, the length and weight of animal carcass from the multi-lumbar vertebrae sheep increased by 2.21 cm and 0.78 kg, respectively. We performed high-throughput genome resequencing on multi-lumbar vertebrae (n = 18) and normal (n = 11) Duolang sheep and obtained a total of more than 528.87 GB data. We found that the most significantly selective region were located in the 49.68-49.74 MB of chromosome 4 by selective-sweep analysis. We annotated this region and found that it contains SFRP4 which is known to regulate bone development. We further used the PCR-SSCP technology to detect the single nucleotide polymorphism (SNP) of the putative candidate SFRP4 and found that the two SNPs (rs600370085:C > T and rs415133338: A > G) of this gene were significantly associated with the multi-lumbar vertebrae of Duolang sheep. Our study indicates that the SFRP4 may be a potential major gene that affects the number of lumbar vertebrae in Duolang sheep, and has the potential to be utilized for sheep breeding in the future.

Identification and validation of a regulatory mutation upstream of the BMP2 gene associated with carcass length in pigs

Background

Carcass length is very important for body size and meat production for swine, thus understanding the genetic mechanisms that underly this trait is of great significance in genetic improvement programs for pigs. Although many quantitative trait loci (QTL) have been detected in pigs, very few have been fine-mapped to the level of the causal mutations. The aim of this study was to identify potential causal single nucleotide polymorphisms (SNPs) for carcass length by integrating a genome-wide association study (GWAS) and functional assays.

Results

Here, we present a GWAS in a commercial Duroc × (Landrace × Yorkshire) (DLY) population that reveals a prominent association signal (P = 4.49E−07) on pig chromosome 17 for carcass length, which was further validated in two other DLY populations. Within the detected 1 Mb region, the BMP2 gene stood out as the most likely causal candidate because of its functions in bone growth and development. Whole-genome gene expression studies showed that the BMP2 gene was differentially expressed in the cartilage tissues of pigs with extreme carcass length. Then, we genotyped an additional 267 SNPs in 500 selected DLY pigs, followed by further whole-genome SNP imputation, combined with deep genome resequencing data on multiple pig breeds. Reassociation analyses using genotyped and imputed SNP data revealed that the rs320706814 SNP, located approximately 123 kb upstream of the BMP2 gene, was the strongest candidate causal mutation, with a large association with carcass length, with a ~ 4.2 cm difference in length across all three DLY populations (N = 1501; P = 3.66E−29). This SNP segregated in all parental lines of the DLY (Duroc, Large White and Landrace) and was also associated with a significant effect on body length in 299 pure Yorkshire pigs (P = 9.2E−4), which indicates that it has a major value for commercial breeding. Functional assays showed that this SNP is likely located within an enhancer and may affect the binding affinity of transcription factors, thereby regulating BMP2 gene expression.

Conclusions

Taken together, these results suggest that the rs320706814 SNP on pig chromosome 17 is a putative causal mutation for carcass length in the widely used DLY pigs and has great value in breeding for body size in pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00689-0.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Analysis of the eighth intron polymorphism of NR6A1 gene in sheep and its correlation with lumbar spine number

For revealing molecular markers related to the traits of multiple lumbar vertebrae in sheep, we analyze the relationship between NR6A1 gene polymorphism and lumbar vertebrae number traits in Xinjiang Kazakh sheep. Lumbar muscle tissues were collected from 6-lumbar spine (L6) Kazak sheep and 7-lumbar spine (L7) Kazak sheep and the intron-8 of NR6A1 gene was amplified by PCR. The SNP locus was detected by the PCR-SSCP method. One-Way ANOVA and an Independent Chi-square Test is adopted to analyze the genotype association with lumbar spine number variation. There were two SNP loci in the intron-8 of the NR6A1 gene: IVS8-188 and IVS8-281. One-Way ANOVA and Independent Chi-square Test indicated a significant association between IVS8-281 and lumbar spine number. The SNP locus of NR6A1 gene intron 8 (IVS8-281G > A) could play a certain role in the variation of lumbar spine number in Xinjiang Kazakh sheep and demonstrates potential to accelerate the sheep breeding of selection process.

Genome Divergence and Dynamics in the Thin-Tailed Desert Sheep From Sudan

With climate change bound to affect food and feed production, emphasis will shift to resilient and adapted indigenous livestock to sustain animal production. However, indigenous livestock comprise several varieties, strains and ecotypes whose genomes are poorly characterized. Here, we investigated genomic variation in an African thin-tailed Desert Sheep sampled in Sudan, using 600K genotype data generated from 92 individuals representing five ecotypes. We included data from 18 fat-tailed and 45 thin-tailed sheep from China, to investigate shared ancestry and perform comparative genomic analysis. We observed a clear genomic differentiation between the African thin-tailed Desert Sheep and the Chinese thin-tailed and fat-tailed sheep, suggesting a broad genetic structure between the fat-tailed and thin-tailed sheep in general, and that at least two autosomal gene pools comprise the genome profile of the thin-tailed sheep. Further analysis detected two distinct genetic clusters in both the African thin-tailed Desert Sheep and the Chinese thin-tailed sheep, suggesting a fine-scale and complex genome architecture in thin-tailed sheep. Selection signature analysis suggested differences in adaptation, production, reproduction and morphology likely underly the fine-scale genetic structure in the African thin-tailed Desert Sheep. This may need to be considered in designing breeding programs and genome-wide association studies.

Differentiating Pigs from Wild Boars Based on NR6A1 and MC1R Gene Polymorphisms

Simple Summary

Wild boar meat is much more expensive than pork. Therefore, there are cases when pork is added to wild boar meat products, but this information is not included on the product label. Currently, there is no fully reliable method that would allow the identification of wild boar and domestic swine products. In this study, we tested the possibility of distinguishing two subspecies using polymorphisms within the MC1R and NR6A1 genes. For this purpose, we used two techniques commonly used in molecular biology, PCR-RFLP and Real-time PCR.

Abstract

This preliminary study aimed to differentiate domestic pigs from wild boars based on MC1R and NR6A1 polymorphisms and to identify admixture between these genomes. We studied samples obtained from wild boars from two regions of Poland and five pig breeds: Polish Landrace, Polish Large White, Złotnicka White, Pulawska and Duroc. Along the MC1R gene sequence, we identified four polymorphic loci comprising three codons. The “wild type” allele was primarily found in wild boar but also in the Duroc and Złotnicka White breeds. Non-wild type alleles were identified in the vast majority of domestic pig samples and in two wild boar samples. Based on MC1R profiles, we conducted a population study, and revealed admixture between both genomes using STRUCTURE and NETWORK Software. Interestingly, an allelic discrimination assay with NR6A1 g.748C > T TaqMan probes revealed a clear separation of samples into two groups: wild boar samples representing the C allele and domestic breeds representing the T allele. Based on the obtained results, we conclude that NR6A1 g.748C > T is an effective marker for differentiating between wild boars and domestic pigs, where this is supported by MC1R data, to identify admixed profiles. We recommend that a larger sample of genomes is studied to verify this method.

Genome scanning reveals novel candidate genes for vertebral and teat number in the Beijing Black Pig

Porcine vertebral and teat numbers are variable and important economic traits in pig production. However, the quantitative trait loci (QTL) and candidate genes for both of these traits in the Beijing Black Pig are not yet known. In the present study, number of vertebrae and number of teats were obtained for 891 individuals of the Beijing Black Pig and genotyped using the Illumina Porcine 50 K BeadChip. Genome scanning was performed to detect associated variants and candidate genes for both traits using a genome-wide association study by tassel software. For vertebral number, 15 significant SNPs were located on SSC7. According to linkage disequilibrium analysis on SSC7, a haplotype block of 221 kb from 97.4 to 97.6 Mb was shown to contain a good candidate gene ABCD4. Interestingly, on SSC12, we recorded a novel QTL containing three significant SNPs and 34 annotated genes from 24.0 to 25.7 Mb for vertebral number. Of the 34 genes, nine Hox family genes (HOXB 1-7, 9, and 13) were found to be good candidate genes. Using the 34 genes, a gene ontology analysis was performed to detect enrichment of anterior/posterior pattern specification. For teat number, a novel chromosome-wide significant QTL was identified on SSC10. In this QTL region, one significant SNP was identified. The nearest gene, NTRK2, was regarded as a candidate gene. The present results expanded the QTL for vertebral and teat numbers and provided useful molecular markers for breeding in the Beijing Black Pig population.

A further survey of the quantitative trait loci affecting swine body size and carcass traits in five related pig populations

Breeding for good meat quality performance while maintaining large body size and desirable carcass traits has been the major challenge for modern swine selective breeding. To address this goal, in the present work we studied five related populations produced by two commercial breeds (Berkshire and Duroc) and two Chinese breeds (Licha black pig and Lulai black pig). A single-trait GWAS performed on 20 body size and carcass traits using a self-developed China Chip-1 porcine SNP50K BeadChip identified 11 genome-wide significant QTL on nine chromosomes and 22 suggestive QTL on 15 chromosomes. For the 11 genome-wide significant QTL, eight were detected in at least two populations, and the rest were population-specific and only mapped in Shanxia black pig. Most of the genome-wide significant QTL were pleiotropic; for example, the QTL around 75.65 Mb on SSC4 was associated with four traits at genome-wide significance level. After screening the genes within 50 kb of the top SNP for each genome-wide significant QTL, NR6A1 and VRTN were chosen as candidate genes for vertebrae number; PLAG1 and BMP2 were identified as candidate genes for body size; and MC4R was the strong candidate gene for body weight. The four genes have been reported as candidates for thoracic vertebrae number, lumbar vertebrae number, carcass length and body weight respectively in previous studies. The effects of VRTN on thoracic vertebrae number, carcass length and body length have been verified in Shanxia black pig. Therefore, the VRTN genotype could be used in gene-assisted selection, and this could accelerate genetic improvement of body size and carcass traits in Shanxia black pig.

Genomic Analysis Reveals Human-Mediated Introgression From European Commercial Pigs to Henan Indigenous Pigs

Introgression of genetic features from European pigs into Chinese pigs was reported possibly contributing to improvements in productivity traits, such as feed conversion efficiency and body size. However, the genomic differences from European pigs and the potential role of introgression in Henan indigenous pigs remains unclear. In this study, we found significant introgression from European pigs into the genome of Chinese indigenous pigs, especially in Henan indigenous pigs. The introgression in Henan indigenous pigs, particularly in the Nanyang black pig, was mainly derived from Duroc pigs. Most importantly, we found that the NR6A1, GPD2, and CSRNP3 genes were introgressed and reshaped by artificial selection, and these may have contributed to increases in pig body size and feed conversion efficiency. Our results suggest that human-mediated introgression and selection have reshaped the genome of Henan pigs and improved several of their desired traits. These findings contribute to our understanding of the history of Henan indigenous pigs and provide insights into the genetic mechanisms affecting economically important traits in pig populations.

Single-marker and haplotype-based genome-wide association studies for the number of teats in two heavy pig breeds

The number of teats is a reproductive‐related trait of great economic relevance as it affects the mothering ability of the sows and thus the number of properly weaned piglets. Moreover, genetic improvement of this trait is fundamental to parallelly help the selection for increased litter size. We present the results of single‐marker and haplotypes‐based genome‐wide association studies for the number of teats in two large cohorts of heavy pig breeds (Italian Large White and Italian Landrace) including 3990 animals genotyped with the 70K GGP Porcine BeadChip and other 1927 animals genotyped with the Illumina PorcineSNP60 BeadChip. In the Italian Large White population, genome scans identified three genome regions (SSC7, SSC10, and SSC12) that confirmed the involvement of the VRTN gene (as we previously reported) and highlighted additional loci known to affect teat counts, including the FRMD4A and HOXB1 gene regions. A different picture emerged in the Italian Landrace population, with a total of 12 genome regions in eight chromosomes (SSC3, SSC6, SSC8, SSC11, SSC13, SSC14, SSC15, and SSC16) mainly detected via the haplotype‐based genome scan. The most relevant QTL was close to the ARL4C gene on SSC15. Markers in the VRTN gene region were not significant in the Italian Landrace breed. The use of both single‐marker and haplotype‐based genome‐wide association analyses can be helpful to exploit and dissect the genome of the pigs of different populations. Overall, the obtained results supported the polygenic nature of the investigated trait and better elucidated its genetic architecture in Italian heavy pigs.

Redefinition of the Mora Romagnola Pig Breed Herd Book Standard Based on DNA Markers Useful to Authenticate Its "Mono-Breed" Products: An Example of Sustainable Conservation of a Livestock Genetic Resource

Simple Summary

Autochthonous breeds are, in general, well adapted to their production systems in which they have been constituted but they are usually less efficient than commercial breeds. Therefore, conservation strategies of livestock genetic resources should be designed to assure profitability to the farmers. The development of “mono-breed” brand products is one of the most effective actions towards this aim. These products are usually sold at a higher price compared to undifferentiated ones, as the consumers consider positively the link between these breeds and the perceived quality of their products. The premium price, however, also attracts fraudsters that unscrupulously see an economic advantage by selling mis-labelled products to obtain an unjustified additional economic gain. These frauds undermine the whole strategy designed to support a sustainable conservation of autochthonous genetic resources. Mora Romagnola is a local pig breed raised in the north of Italy. Mono-breed pork products derived from this breed are part of an important niche value chain that is intrinsically linked to the conservation of this local genetic resource. In this study we present how the Mora Romagnola Herd Book standard integrated information of DNA markers of two genes (MC1R and NR6A1), affecting morphological traits, to allow the authentication of mono-breed products of this breed. This is one of the first examples of sustainable conservation of a pig genetic resource designed starting from the genotype of the animals registered to the breed herd book, with the specific purpose to combat frauds.

Abstract

Mora Romagnola is an autochthonous pig breed, raised in the north of Italy. Mono-breed pork products of this breed are part of important niche value chain that is intrinsically linked to the conservation of this local genetic resources that can only survive due to the premium price that these products can obtain on the market. However, the added value attracts fraudsters that unscrupulously sell mis-labelled Mora Romagnola products, causing consumer distrust that, in turn, undermines the conservation strategy of this breed. To monitor and better characterise this local breed, we phenotyped 826 Mora Romagnola pigs for three breed-specific traits. Then, we genotyped almost all living sows and boars registered to the Herd Book (n. = 357 animals) for polymorphisms in the MC1R and NR6A1 genes (affecting coat colour and vertebral number, respectively). The results were used to re-define the breed descriptors of the Mora Romagnala breed that included information on the allowed genotypes at these two genes. A few pigs that did not carry the allowed genotypes were excluded from its Herd Book. Finally, we evaluated the usefulness of these DNA markers to authenticate Mora Romagnola meat against meat derived from other 11 pig breeds and wild boars. To our knowledge, the Mora Romagnola Herd Book is one of the first examples that established a direct link between a genetic standard of a breed with the possibility to authenticate mono-breed products using DNA markers with the specific purpose to combat frauds and, indirectly, support the conservation of a livestock genetic resource.

Runs of homozygosity provide a genome landscape picture of inbreeding and genetic history of European autochthonous and commercial pig breeds

ROHs are long stretches of DNA homozygous at each polymorphic position. The proportion of genome covered by ROHs and their length are indicators of the level and origin of inbreeding. Frequent common ROHs within the same population define ROH islands and indicate hotspots of selection. In this work, we investigated ROHs in a total of 1131 pigs from 20 European local pig breeds and in three cosmopolitan breeds, genotyped with the GGP Porcine HD Genomic Profiler. plink software was used to identify ROHs. Size classes and genomic inbreeding parameters were evaluated. ROH islands were defined by evaluating different thresholds of homozygous SNP frequency. A functional overview of breed-specific ROH islands was obtained via over-representation analyses of GO biological processes. Mora Romagnola and Turopolje breeds had the largest proportions of genome covered with ROH (~1003 and ~955 Mb respectively), whereas Nero Siciliano and Sarda breeds had the lowest proportions (~207 and 247 Mb respectively). The highest proportion of long ROH (>16 Mb) was in Apulo-Calabrese, Mora Romagnola and Casertana. The largest number of ROH islands was identified in the Italian Landrace (n = 32), Cinta Senese (n = 26) and Lithuanian White Old Type (n = 22) breeds. Several ROH islands were in regions encompassing genes known to affect morphological traits. Comparative ROH structure analysis among breeds indicated the similar genetic structure of local breeds across Europe. This study contributed to understanding of the genetic history of the investigated pig breeds and provided information to manage these pig genetic resources.

A genome-wide association study for the number of teats in European rabbits (Oryctolagus cuniculus) identifies several candidate genes affecting this trait

In the European rabbit (Oryctolagus cuniculus), a polytocous livestock species, the number of teats indirectly impacts the doe reproduction efficiency and, in turn, the sustainable production of rabbit meat. In this study, we carried out a genome-wide association study (GWAS) for the total number of teats in 247 Italian White does included in the Italian White rabbit breed selection program, by applying a selective genotyping approach. Does had either 8 (n = 121) or 10 teats (n = 126). All rabbits were genotyped with the Affymetrix Axiom OrcunSNP Array. Genomic data from the two extreme groups of rabbits were also analysed with the single-marker fixation index statistic and combined with the GWAS results. The GWAS identified 50 significant SNPs and the fixation index analysis identified a total of 20 SNPs that trespassed the 99.98th percentile threshold, 19 of which confirmed the GWAS results. The most significant SNP (P = 4.31 × 10^-11 ) was located on OCU1, close to the NUDT2 gene, a breast carcinoma cells proliferation promoter. Another significant SNP identified as candidate gene NR6A1, which is well known to play an important role in affecting the correlated number of vertebrae in pigs. Other significant markers were close to candidate genes involved in determining body length in mice. Markers associated with increased number of teats could be included in selection programmes to speed up the improvement for this trait in rabbit lines that need to increase maternal performances.

NR6A1 Allelic Frequencies as an Index for both Miniaturizing and Increasing Pig Body Size

BACKGROUND/AIM

The number of vertebrae in swine varies from 19 to 23 and is associated with body size. Nuclear receptor subfamily 6 group A member 1 (NR6A1) is considered a strong candidate for affecting the number of vertebrae in swine. Wild boars, which uniformly have 19 vertebrae, have the wild type allele while multi-vertebrae European commercial pigs have the mutated allele. Our aim was to confirm the factor of the miniaturization.

MATERIALS AND METHODS

We examined vertebrae number and NR6A1 polymorphism in the Microminipig and three domestic breeds that vary in body size.

RESULTS

The Microminipig had 19 or less vertebrae and a wild type NR6A1 genotype. Three domestic breeds had more than 21 vertebrae while the largest vertebrae number was observed in multi-vertebrae-fixed Large White. Heterozygous genotypes were observed in the middle-sized indigenous pig while homozygous NR6A1 mutations were observed in European commercial breeds.

CONCLUSION

NR6A1 could be a useful index for both miniaturizing and increasing pig body size.

Selective sweep analysis reveals extensive parallel selection traits between large white and Duroc pigs

In the process of pig genetic improvement, different commercial breeds have been bred for the same purpose, improving meat production. Most of the economic traits, such as growth and fertility, have been selected similarly despite the discrepant selection pressure, which is known as parallel selection. Here, 28 whole-genome sequencing data of Danish large white pigs with an approximately 25-fold depth each were generated, resulting in about 12 million high-quality SNPs for each individual. Combined with the sequencing data of 27 Duroc and 23 European wild boars, we investigated the parallel selection of Danish large white and Duroc pigs using two complementary methods, Fst and iHS. In total, 67 candidate regions were identified as the signatures of parallel selection. The genes in candidate regions of parallel selection were mainly associated with sensory perception, growth rate, and body size. Further functional annotation suggested that the striking consistency of the terms may be caused by the polygenetic basis of quantitative traits, and revealing the complex genetic basis of parallel selection. Besides, some unique terms were enriched in population-specific selection regions, such as the limb development-related terms enriched in Duroc-specific selection regions, suggesting unique selections of breed-specific selected traits. These results will help us better understand the parallel selection process of different breeds. Moreover, we identified several potential causal SNPs that may contribute to the pig genetic breeding process.

Population genome of the newly discovered Jinchuan yak to understand its adaptive evolution in extreme environments and generation mechanism of the multirib trait

The adaptation and diversity of animals to the extreme environments of the Qinghai-Tibet Plateau (QTP) are typical materials to study adaptive evolution. The recently discovered Jinchuan yak population has many individuals with multiple ribs. However, little is known about this yak's origin, evolution, and the genetic mechanisms that formed its unique multirib trait. Here, we report a valuable population genome resource of the Jinchuan yak by resequencing the whole genome of 150 individuals. Population genetic polymorphism and structure analysis reveal that Jinchuan yak can be differentiated as a unique and original yak population among the domestic yak. Combined with geological change, the Jinchuan yak's evolutionary origin is speculated to be about 6290 years ago, which may be related to the unique geographical environment of the eastern edge of the QTP during this period. Compared with other domestic yaks, this new population has 280 positively selected genes. The genes related to skeletal function hold a considerable and remarkable proportion, suggesting that the specific skeletal characteristics have been enhanced in the adaptive evolution of Jinchuan yak in the extreme plateau environment. The genome-wide association study has revealed that TUBA8 and TUBA4A, the genes that regulate the cytoskeleton, are potential genes associated with the multirib trait. Our findings provide a basis to further understand the generation mechanism of the adaptive evolution of this new population in high-altitude extreme environments and the multivertebrate trait of domestic animals.

Whole-genome sequencing of European autochthonous and commercial pig breeds allows the detection of signatures of selection for adaptation of genetic resources to different breeding and production systems

Background

Natural and artificial directional selection in cosmopolitan and autochthonous pig breeds and wild boars have shaped their genomes and resulted in a reservoir of animal genetic diversity. Signatures of selection are the result of these selection events that have contributed to the adaptation of breeds to different environments and production systems. In this study, we analysed the genome variability of 19 European autochthonous pig breeds (Alentejana, Bísara, Majorcan Black, Basque, Gascon, Apulo-Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano, Sarda, Krškopolje pig, Black Slavonian, Turopolje, Moravka, Swallow-Bellied Mangalitsa, Schwäbisch-Hällisches Schwein, Lithuanian indigenous wattle and Lithuanian White old type) from nine countries, three European commercial breeds (Italian Large White, Italian Landrace and Italian Duroc), and European wild boars, by mining whole-genome sequencing data obtained by using a DNA-pool sequencing approach. Signatures of selection were identified by using a single-breed approach with two statistics [within-breed pooled heterozygosity (H_P) and fixation index (F_ST)] and group-based F_ST approaches, which compare groups of breeds defined according to external traits and use/specialization/type.

Results

We detected more than 22 million single nucleotide polymorphisms (SNPs) across the 23 compared populations and identified 359 chromosome regions showing signatures of selection. These regions harbour genes that are already known or new genes that are under selection and relevant for the domestication process in this species, and that affect several morphological and physiological traits (e.g. coat colours and patterns, body size, number of vertebrae and teats, ear size and conformation, reproductive traits, growth and fat deposition traits). Wild boar related signatures of selection were detected across all the genome of several autochthonous breeds, which suggests that crossbreeding (accidental or deliberate) occurred with wild boars.

Conclusions

Our findings provide a catalogue of genetic variants of many European pig populations and identify genome regions that can explain, at least in part, the phenotypic diversity of these genetic resources.

Polymorphism of VRTN Gene g.20311_20312ins291 Was Associated with the Number of Ribs, Carcass Diagonal Length and Cannon Bone Circumference in Suhuai Pigs

Simple Summary

An increase in the number of ribs (RIB) could improve carcass length (CL) and body size. Cannon bone circumference (CBC) is a pivotal body size trait, and a large CBC could enhance the capacity to bear excessive body weight, vigorous exercise, and resistance to injuries. Several researchers showed that the vertnin (VRTN) gene g.20311_20312ins291 (NC_010449.5 7: g.20311_20312ins291) is an important variant that is related to RIB and CL of Western pigs. However, it is unknown whether this variant could affect the CBC of pigs. Our study showed that this variant was significantly associated with RIB, carcass diagonal length (CDL), and CBC in Suhuai pigs; therefore, it could be used as a potential molecular marker for improving RIB, CDL, and CBC in this breed.

Abstract

The vertnin (VRTN) gene g.20311_20312ins291 was reported as an important variant related to the number of ribs (RIB), and the ins/ins genotype was advantageous for improving RIB of Western pigs. The purpose of this study was to determine whether the VRTN gene g.20311_20312ins291 influences RIB, carcass traits, and body size traits, including cannon bone circumference (CBC) in Chinese Suhuai pigs. We found that the VRTN gene g.20311_20312ins291 was polymorphic in Suhuai fattening pigs and gilts. The polymorphism of g.20311_20312ins291 was significantly associated with RIB and CDL in Suhuai fattening pigs (p < 0.01), whereas this variant had no influence on carcass weight (CWT). There was a tendency of association between this variant and carcass straight length (CSL) in Suhuai fattening pigs (p = 0.06). The polymorphism of g.20311_20312ins291 was also significantly associated with CBC in Suhuai gilts (p = 0.04). Furthermore, CBC was positively genetically correlated with body length (0.22, p < 0.01) and body weight (0.15, p < 0.01). Our results indicated that the VRTN gene g.20311_20312ins291 could be used as a potential marker for improving RIB, CDL, and CBC in Suhuai pigs.

Matching STR and SNP genotyping to discriminate between wild boar, domestic pigs and their recent hybrids for forensic purposes

The genetic discrimination between phylogenetically close taxa can be challenging if their gene pools are not differentiated and there are many shared polymorphisms. The gene flow between wild boar (Sus scrofa) and domestic pig (S. s. domesticus) has never been interrupted from domestication onwards, due to non-stop natural and human-mediated crossbreeding. To date there are no individual genetic markers that are able to distinguish between the two forms, nor even to identify effectively their hybrids. We developed a combined molecular protocol based on multiplex porcine-specific STR-profiling system and new real time PCR-based assays of single polymorphisms in the NR6A1 and MC1R genes to gain high diagnostic power in the differentiation of wild boar, pig and hybrids for forensic purposes. The combined approach correctly assigned individuals to one or the other parental gene pool and identified admixed genotypes. Evidence was found for substantial reduction of false negative results by using multiple marker systems jointly, compared to their use individually. Our protocol is a powerful and cost-effective diagnostic tool that can easily be adopted by most forensic laboratories to assist authorities contrast food adulteration, assure veterinary public health and fight against wildlife crimes, like poaching and illegal detention of wild animals.

A Novel 13 bp Deletion within the NR6A1 Gene Is Significantly Associated with Growth Traits in Donkeys

Simple Summary

The detection of genes potentially associated with economic traits and identification of effective variants can provide a basis for molecular marker-assisted selection of livestock. NR6A1 is a member of the nuclear receptor family and is an important candidate gene related to body size traits. Previous studies showed that NR6A1 gene was associated with body size traits in pigs and other livestock, however, it has not yet been observed in donkeys. In the current study, a 13 bp deletion in NR6A1 gene was firstly identified in donkeys. Analysis showed that this deletion had significant associations with body size traits.

Abstract

Nuclear receptor subfamily 6, group A, member 1 (NR6A1), as an important member of the nuclear receptor family, plays an important role in regulating growth, metabolism, and differentiation of embryonic stem cells. For this reason, the NR6A1 gene is considered to be a promising candidate for economic traits and was found to be associated with body size traits in many livestock. However, no studies have been conducted on NR6A1 in donkeys so far. Thus, in this research, we focused on donkeys and identified a 13 bp deletion in intron-1 of the NR6A1 gene among 408 individuals from Guanzhong and Dezhou donkeys using polyacrylamide gel electrophoresis. Three genotypes were identified, namely II, ID, and DD. The association analysis indicated that the body lengths and body heights5f genotype II individuals were significantly different to those of genotype ID in Dezhou donkeys. Conclusively, the 13 bp deletion was associated with growth traits in both Guanzhong donkeys and Dezhou donkeys, indicating that the NR6A1 gene could be a possible candidate gene in marker-assisted selection for donkey breeding programs.

Generalized glycogenosis in Brahman-derived breeds: diagnosis and prevalence in Argentina

Generalized glycogenosis is a lethal autosomal recessive disease caused by a deficient activity of the acidic 1,4-α-glucosidase enzyme and characterized by an accumulation of glycogen within lysosomes. Three mutations in the GAA gene causing bovine generalized glycogenosis have been identified in two cattle breeds, Brahman and Shorthorn. The objective of this study was to evaluate the prevalence of carriers of the E7 mutation in the GAA gene in Argentinean Brahman-derived herds. A total of 930 Braford, 94 Brangus, and 8 Brahman samples were analyzed. The genotyping was done by polymerase chain reaction and restriction fragment length polymorphism (PCR/RFLP). We found that 12.02% (95% CI 12.00-12.04) of the total number of samples received were heterozygous (i.e., carriers) for the E7 mutation, while 12.58% (95% CI 12.56-12.60) of the Braford, 6.38% (95% CI 6.26-6.51) of the Brangus, and 12.50% (95% CI 9.82-15.18) of the Brahman samples were carriers of this loss-of-function allele. Neither breed nor sex were significantly associated to the presence of the mutation. The prevalence informed in this study is similar to the average prevalence reported for Australian Brahmans. The finding of heterozygous animals suggests that breeders and insemination centers should continue screening their herds to minimize the dissemination of this deleterious allele.