Comprehensive variation discovery and recovery of missing sequence in the pig genome using multiple de novo assemblies

Porcine low-density lipoprotein receptor plays an important role in classical swine fever virus infection

Several cellular factors have been reported to be required for replication of classical swine fever virus (CSFV), a member of the genus Pestivirus within the family Flaviviridae. However, many steps of its replication cycle are still poorly understood. The low-density lipoprotein receptor (LDLR) is involved in cell entry and post-entry processes of different viruses including other members of the Flaviviridae. In this study, the relevance of LDLR in replication of CSFV and another porcine pestivirus, Bungowannah pestivirus (BuPV), was investigated by antibody-mediated blocking of LDLR and genetically engineered porcine cell lines providing altered LDLR expression levels. An LDLR-specific antibody largely blocked infection with CSFV, but had only a minor impact on BuPV. Infections of the genetically modified cells confirmed an LDLR-dependent replication of CSFV. Compared to wild type cells, lower and higher expression of LDLR resulted in a 3.5-fold decrease or increase in viral titers already 20 h post infection. Viral titers were 25-fold increased in LDLR-overexpressing cells compared to cells with reduced LDLR expression at 72 h post infection. The varying LDLR expression levels had no clear effect on permissivity to BuPV. A decoy receptor assay using recombinant soluble LDLR provided no evidence that LDLR may function as a receptor for CSFV or BuPV. Differences in their dependency on LDLR suggest that CSFV and BuPV likely use different mechanisms to interact with their host cells. Moreover, this study reveals similarities in the replication cycles of CSFV and other members of the family Flaviviridae that are dependent on LDLR.

Effect of DHCR7 on adipocyte differentiation in goats

Cholesterol is regarded as a signaling molecule in regulating the metabolism and function of fat cells, in which 7-Dehydrocholesterol reductase (DHCR7) is a key enzyme that catalyzes the conversion of 7-dehydrocholesterol to cholesterol, however, the exact function of DHCR7 in goat adipocytes remains unknown. Here, the effect of DHCR7 on the formation of subcutaneous and intramuscular fat in goats was investigated in vitro, and the result indicated that the mRNA level of DHCR7 showed a gradual downward trend in subcutaneous adipogenesis, but an opposite trend in intramuscular adipogenesis. In the process of subcutaneous preadipocytes differentiation, overexpression of DHCR7 inhibited the expression of adipocytes differentiation marker genes (CEBP/α, CEBP/β, SREBP1 and AP2), lipid metabolism-related genes (AGPAT6, FASN, SCD1 and LPL), and the lipid accumulation. However, in intramuscular preadipocyte differentiation, DHCR7 overexpression showed a promoting effect on adipocyte differentiation marker genes (CEBP/α, CEBP/β, PPARγ and SREBP1) and lipid metabolism-related genes (GPAM, AGPAT6, DGAT1 and SCD1) expression, and on lipid accumulation. In summary, our work demonstrated that DHCR7 played an important role in regulating adipogenic differentiation and lipid metabolism in preadipocytes in goats, which is of great significance for uncovering the underlying molecular mechanism of adipocyte differentiation and improving goat meat quality.

Whole-genome resequencing of Hu sheep identifies candidate genes associated with agronomic traits

The phenotypic diversity resulting from artificial or natural selection of sheep has made a significant contribution to human civilization. Hu sheep are a local sheep breed unique to China with high reproductive rates and rapid growth. Genomic selection signatures have been widely used to investigate the genetic mechanisms underlying phenotypic variation in livestock. Here, we conduct whole-genome sequencing of 207 Hu sheep and compare them with the wild ancestors of domestic sheep (Asiatic mouflon) to investigate the genetic characteristics and selection signatures of Hu sheep. Based on six signatures of selection approaches, we detect genomic regions containing genes related to reproduction (BMPR1B, BMP2, PGFS, CYP19, CAMK4, GGT5, and GNAQ), vision (ALDH1A2, SAG, and PDE6B), nervous system (NAV1), and immune response (GPR35, SH2B2, PIK3R3, and HRAS). Association analysis with a population of 1299 Hu sheep reveals that those missense mutations in the GPR35 (GPR35 g.952651 A>G; GPR35 g.952496 C>T) and NAV1 (NAV1 g.84216190 C>T; NAV1 g.84227412 G>A) genes are significantly associated (P < 0.05) with immune and growth traits in Hu sheep, respectively. This research offers unique insights into the selection characteristics of Hu sheep and facilitates further genetic improvement and molecular investigations.

Uncovering the Interaction between TRAF1 and MAVS in the RIG-I Pathway to Enhance the Upregulation of IRF1/ISG15 during Classical Swine Fever Virus Infection

Classical swine fever (CSF) is caused by the classical swine fever virus (CSFV), which poses a threat to swine production. The activation of host innate immunity through linker proteins such as tumor necrosis factor receptor (TNF-R)-associated factor (TRAF) is crucial for the induction of the NF-κB pathway. Recent research has revealed the involvement of mitochondrial antiviral-signaling protein (MAVS) in the interaction with TRAF2, 3, 5, and 6 to activate both the NF-κB and IRF3 pathways. This study revealed that CSFV infection led to the upregulation of TRAF1 mRNA and protein levels; moreover, TRAF1 overexpression inhibited CSFV replication, while TRAF1 knockdown promoted replication, highlighting its importance in the host response to CSFV infection. Additionally, the expression of RIG-I, MAVS, TRAF1, IRF1, and ISG15 were detected in PK-15 cells infected with CSFV, revealing that TRAF1 plays a role in regulating IRF1 and ISG15 within the RIG-I pathway. Furthermore, Co-IP, GST pull-down, and IFA analyses demonstrated that TRAF1 interacted with MAVS and co-localized in the cytoplasm during CSFV infection. Ultimately, TRAF1 acted as a novel member of the TRAF family, bound to MAVS as a linker molecule, and functioned as a mediator downstream of MAVS in the RIG-I/MAVS pathway against CSFV replication.

Genomic, transcriptomic and epigenomic analysis towards the understanding of porcine semen quality traits. Past, current and future trends

The importance of boar reproductive traits, including semen quality, in the sustainability of pig production system is increasingly being acknowledged by academic and industrial sectors. Research is needed to understand the biology and genetic components underlying these traits so that they can be incorporated into selection schemes and managerial decisions. This article reviews our current understanding of genome biology and technologies for genome, transcriptome and epigenome analysis which now facilitate the identification of causal variants affecting phenotypes more than ever before. Genetic and transcriptomic analysis of candidate genes, Genome-Wide Association Studies, expression microarrays, RNA-Seq of coding and noncoding genes and epigenomic evaluations have been conducted to profile the molecular makeups of pig sperm. These studies have provided insightful information for a several semen-related parameters. Nonetheless, this research is still incipient. The spermatozoon harbors a reduced transcriptome and highly modified epigenome, and it is assumed to be transcriptionally silent for nuclear gene expression. For this reason, the extent to which the sperm's RNA and epigenome recapitulate sperm biology and function is unclear. Hence, we anticipate that single-cell level analyses of the testicle and other male reproductive organs, which can reveal active transcription and epigenomic profiles in cells influencing sperm quality, will gain popularity and markedly advance our understanding of sperm-related traits. Future research will delve deeper into sperm fertility, boar resilience to environmental changes or harsh conditions, especially in the context of global warming, and also in transgenerational inheritance and how the environment influences the sperm transcriptome and epigenome.

A chromosome-level genome of Chenghua pig provides new insights into the domestication and local adaptation of pigs

As a country with abundant genetic resources of pigs, the domestication history of pigs in China and the adaptive evolution of Chinese pig breeds at different latitudes have rarely been elucidated at the genome-wide level. To fill this gap, we first assembled a high-quality chromosome-level genome of the Chenghua pig and used it as a benchmark to analyse the genomes of 272 samples from three genera of three continents. The divergence of the three species belonging to three genera, Phacochoerus africanus, Potamochoerus porcus, and Sus scrofa, was assessed. The introgression of pig breeds redefined that the migration routes were basically from southern China to central and southwestern China, then spread to eastern China, arrived in northern China, and finally reached Europe. The domestication of pigs in China occurred ∼12,000 years ago, earlier than the available Chinese archaeological domestication evidence. In addition, FBN1 and NR6A1 were identified in our study as candidate genes related to extreme skin thickness differences in Eurasian pig breeds and adaptive evolution at different latitudes in Chinese pig breeds, respectively. Our study provides a new resource for the pig genomic pool and refines our understanding of pig genetic diversity, domestication, migration, and adaptive evolution at different latitudes.

Potential Role of Pig UCP3 in Modulating Adipocyte Browning via the Beta-Adrenergic Receptor Signaling Pathway

Adipose tissue plays an important role in regulating body temperature and metabolism, with white adipocytes serving as storage units for energy. Recent research focused on the browning of white adipocytes (beige adipocytes), causing thermogenesis and lipolysis. The process of browning is linked to the activation of uncoupling protein (UCP) expression, which can be mediated by the β3 adrenergic receptor pathway. Transcriptional factors, such as peroxisome proliferator activated receptor γ (PPARγ) and PPARγ coactivator 1 alpha, play vital roles in cell fate determination for fat cells. Beige adipocytes have metabolic therapeutic potential to combat diseases such as obesity, diabetes mellitus, and dyslipidemia, owing to their significant impact on metabolic functions. However, the molecular mechanisms that cause the induction of browning are unclear. Therefore, research using animal models and primary culture is essential to provide an understanding of browning for further application in human metabolic studies. Pigs have physiological similarities to humans; hence, they are valuable models for research on adipose tissue. This study demonstrates the browning potential of pig white adipocytes through primary culture experiments. The results show that upregulation of UCP3 gene expression and fragmentation of lipid droplets into smaller particles occur due to isoproterenol stimulation, which activates beta-adrenergic receptor signaling. Furthermore, PPARγ and PGC-1α were found to activate the UCP3 promoter region, similar to that of UCP1. These findings suggest that pigs undergo metabolic changes that induce browning in white adipocytes, providing a promising approach for metabolic research with potential implications for human health. This study offers valuable insights into the mechanism of adipocyte browning using pig primary culture that can enhance our understanding of human metabolism, leading to cures for commonly occurring diseases.

Constructing eRNA-mediated gene regulatory networks to explore the genetic basis of muscle and fat-relevant traits in pigs

BACKGROUND

Enhancer RNAs (eRNAs) play a crucial role in transcriptional regulation. While significant progress has been made in understanding epigenetic regulation mediated by eRNAs, research on the construction of eRNA-mediated gene regulatory networks (eGRN) and the identification of critical network components that influence complex traits is lacking.

RESULTS

Here, employing the pig as a model, we conducted a comprehensive study using H3K27ac histone ChIP-seq and RNA-seq data to construct eRNA expression profiles from multiple tissues of two distinct pig breeds, namely Enshi Black (ES) and Duroc. In addition to revealing the regulatory landscape of eRNAs at the tissue level, we developed an innovative network construction and refinement method by integrating RNA-seq, ChIP-seq, genome-wide association study (GWAS) signals and enhancer-modulating effects of single nucleotide polymorphisms (SNPs) measured by self-transcribing active regulatory region sequencing (STARR-seq) experiments. Using this approach, we unraveled eGRN that significantly influence the growth and development of muscle and fat tissues, and identified several novel genes that affect adipocyte differentiation in a cell line model.

CONCLUSIONS

Our work not only provides novel insights into the genetic basis of economic pig traits, but also offers a generalizable approach to elucidate the eRNA-mediated transcriptional regulation underlying a wide spectrum of complex traits for diverse organisms.

Haplotype-resolved 3D chromatin architecture of the hybrid pig

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

Characterization of the complete mitochondrial genome and identification of signature sequence of Indian wild pig

Mitochondrial DNA (mtDNA) serves as a valuable molecular marker for constructing matrilineal genealogies and tracing the evolutionary history of animals. This study aimed to characterize the complete mitochondrial genome of the Indian wild pig (IWB) (Sus scrofa cristatus) and identify IWB-specific DNA sequences that could be used as genomic signatures to differentiate IWB from domestic Indian pigs (IDP) in forensic cases. For the purpose, three wild IWB from a rescue centre were used for the characterization of the mitochondrial genome of the IWB. The mitochondrial genome was sequenced by the primer walking technique using 30 overlapping primers. The mitochondrial genome of the IWB was found to be 16,689 bp long containing 37 genes coding for 2 rRNAs, 22 tRNAs, 13 protein coding genes, and 1 D-loop region similar to the mitogenome of other pigs. Sequence analysis of the D-loop of IWB with other IDP indicated some signature sequence for IWB like duplication and transition event from 1090th to 1099th position, deletion of a 10 bp sequence at the 755th position, insertion of (CA) at the 137th position, and substitution of AT to GA at the 638th position. These variations specially the duplication along with transition event causes creation of unique signature sequence (-ACACAAACCT-) in the IWB that could serve as signature sequences for the IWB and be used as markers for differentiation of IWB from IDP breeds in academic as well as forensic or vetero-legal cases. Overall, a total of 36 polymorphic positions were identified in the IWB, with 29 sites being unique to the IWB only and seven being common to the Doom and HDK75 pig breeds. None of the common polymorphic sites were identified in prevailing domestic pig populations. Phylogenetic analysis of the mitochondrial genome revealed the distinct separation of the IWB from IDP. The results of genetic distance evaluation showed that the Doom pig breed was the closest to the IWB. This study provides valuable insights into the mitogenome characterisation, signature sequence and genetic distance analysis of the IWB and establishes a foundation for future studies on the conservation of this protected species.

Assembly of novel sequences for Chinese domestic pigs reveals new genes and regulatory variants providing new insights into their diversity

There is an increasing understanding that a reference genome representing an individual cannot capture all the gene repertoire of a species. Here, we conduct a population-scale missing sequences detection of Chinese domestic pigs using whole-genome sequencing data from 534 individuals. We identify 132.41 Mb of sequences absent in the reference assembly, including eight novel genes. In particular, the breeds spread in Chinese high-altitude regions perform significantly different frequencies of new sequences in promoters than other breeds. Furthermore, we dissect the role of non-coding variants and identify a novel sequence inserted in the 3'UTR of the FMO3 gene, which may be associated with the intramuscular fat phenotype. This novel sequence could be a candidate marker for meat quality. Our study provides a comprehensive overview of the missing sequences in Chinese domestic pigs and indicates that this dataset is a valuable resource for understanding the diversity and biology of pigs.

Pig pangenome graph reveals functional features of non-reference sequences

BACKGROUND

The reliance on a solitary linear reference genome has imposed a significant constraint on our comprehensive understanding of genetic variation in animals. This constraint is particularly pronounced for non-reference sequences (NRSs), which have not been extensively studied.

RESULTS

In this study, we constructed a pig pangenome graph using 21 pig assemblies and identified 23,831 NRSs with a total length of 105 Mb. Our findings revealed that NRSs were more prevalent in breeds exhibiting greater genetic divergence from the reference genome. Furthermore, we observed that NRSs were rarely found within coding sequences, while NRS insertions were enriched in immune-related Gene Ontology terms. Notably, our investigation also unveiled a close association between novel genes and the immune capacity of pigs. We observed substantial differences in terms of frequencies of NRSs between Eastern and Western pigs, and the heat-resistant pigs exhibited a substantial number of NRS insertions in an 11.6 Mb interval on chromosome X. Additionally, we discovered a 665 bp insertion in the fourth intron of the TNFRSF19 gene that may be associated with the ability of heat tolerance in Southern Chinese pigs.

CONCLUSIONS

Our findings demonstrate the potential of a graph genome approach to reveal important functional features of NRSs in pig populations.

When Livestock Genomes Meet Third-Generation Sequencing Technology: From Opportunities to Applications

Third-generation sequencing technology has found widespread application in the genomic, transcriptomic, and epigenetic research of both human and livestock genetics. This technology offers significant advantages in the sequencing of complex genomic regions, the identification of intricate structural variations, and the production of high-quality genomes. Its attributes, including long sequencing reads, obviation of PCR amplification, and direct determination of DNA/RNA, contribute to its efficacy. This review presents a comprehensive overview of third-generation sequencing technologies, exemplified by single-molecule real-time sequencing (SMRT) and Oxford Nanopore Technology (ONT). Emphasizing the research advancements in livestock genomics, the review delves into genome assembly, structural variation detection, transcriptome sequencing, and epigenetic investigations enabled by third-generation sequencing. A comprehensive analysis is conducted on the application and potential challenges of third-generation sequencing technology for genome detection in livestock. Beyond providing valuable insights into genome structure analysis and the identification of rare genes in livestock, the review ventures into an exploration of the genetic mechanisms underpinning exemplary traits. This review not only contributes to our understanding of the genomic landscape in livestock but also provides fresh perspectives for the advancement of research in this domain.

An Analysis of a Transposable Element Expression Atlas during 27 Developmental Stages in Porcine Skeletal Muscle: Unveiling Molecular Insights into Pork Production Traits

The development and growth of porcine skeletal muscle determine pork quality and yield. While genetic regulation of porcine skeletal muscle development has been extensively studied using various omics data, the role of transposable elements (TEs) in this context has been less explored. To bridge this gap, we constructed a comprehensive atlas of TE expression throughout the developmental stages of porcine skeletal muscle. This was achieved by integrating porcine TE genomic coordinates with whole-transcriptome RNA-Seq data from 27 developmental stages. We discovered that in pig skeletal muscle, active Tes are closely associated with active epigenomic marks, including low levels of DNA methylation, high levels of chromatin accessibility, and active histone modifications. Moreover, these TEs include 6074 self-expressed TEs that are significantly enriched in terms of muscle cell development and myofibril assembly. Using the TE expression data, we conducted a weighted gene co-expression network analysis (WGCNA) and identified a module that is significantly associated with muscle tissue development as well as genome-wide association studies (GWAS) of the signals of pig meat and carcass traits. Within this module, we constructed a TE-mediated gene regulatory network by adopting a unique multi-omics integration approach. This network highlighted several established candidate genes associated with muscle-relevant traits, including HES6, CHRNG, ACTC1, CHRND, MAMSTR, and PER2, as well as novel genes like ENSSSCG00000005518, ENSSSCG00000033601, and PIEZO2. These novel genes hold promise for regulating muscle-related traits in pigs. In summary, our research not only enhances the TE-centered dissection of the genetic basis underlying pork production traits, but also offers a general approach for constructing TE-mediated regulatory networks to study complex traits or diseases.

The jigsaw puzzle of pedigree: whole-genome resequencing reveals genetic diversity and ancestral lineage in Sunong black pigs

The Sunong black pig is a new composite breed under development generated from Chinese indigenous pig breeds (i.e., Taihu and Huai) and intensive pig breeds (i.e., Landrace and Berkshire), which is an important genetic material for studying breeding mechanisms. However, there is currently limited knowledge about the genetic structure and germplasm characteristics of Sunong black pigs. To comprehensively understand their genetic composition and ancestry proportions, we performed population structure and local ancestry inference analysis based on whole-genome sequencing information. The results showed that Sunong black pigs could be clustered independently into a group, whose pedigree was intermediate between indigenous and commercial pig breeds, but closer to commercial pigs. Furthermore, local ancestry inference analysis revealed that Sunong black pigs inherited immune and reproductive traits from indigenous pig breeds, including CC and CXC chemokine family, Toll-like receptor family, IFN gene family, ESR1, AREG and EREG gene, while growth and development-related traits were inherited from commercial pig breeds, including IGF1 and GSY2 gene. Overall, Sunong black pigs have formed a relatively stable genome structure with some advantageous traits inherited from their ancestral breeds. This study deepened the understanding of the breeding mechanism of Sunong black pigs and provided a reference for cross-breeding programmes in livestock.

In Silico Investigation of the Clinical Translatability of Competitive Clearance Glucose-Responsive Insulins

The glucose-responsive insulin (GRI) MK-2640 from Merck was a pioneer in its class to enter the clinical stage, having demonstrated promising responsiveness in in vitro and preclinical studies via a novel competitive clearance mechanism (CCM). The smaller pharmacokinetic response in humans motivates the development of new predictive, computational tools that can improve the design of therapeutics such as GRIs. Herein, we develop and use a new computational model, IM3PACT, based on the intersection of human and animal model glucoregulatory systems, to investigate the clinical translatability of CCM GRIs based on existing preclinical and clinical data of MK-2640 and regular human insulin (RHI). Simulated multi-glycemic clamps not only validated the earlier hypothesis of insufficient glucose-responsive clearance capacity in humans but also uncovered an equally important mismatch between the in vivo competitiveness profile and the physiological glycemic range, which was not observed in animals. Removing the inter-species gap increases the glucose-dependent GRI clearance from 13.0% to beyond 20% for humans and up to 33.3% when both factors were corrected. The intrinsic clearance rate, potency, and distribution volume did not apparently compromise the translation. The analysis also confirms a responsive pharmacokinetics local to the liver. By scanning a large design space for CCM GRIs, we found that the mannose receptor physiology in humans remains limiting even for the most optimally designed candidate. Overall, we show that this computational approach is able to extract quantitative and mechanistic information of value from a posteriori analysis of preclinical and clinical data to assist future therapeutic discovery and development.

Genetic variation and domestication of horses revealed by 10 chromosome-level genomes and whole-genome resequencing

Understanding the genetic variations of the horse (Equus caballus) genome will improve breeding conservation and welfare. However, genetic variations in long segments, such as structural variants (SVs), remain understudied. We de novo assembled 10 chromosome-level three-dimensional horse genomes, each representing a distinct breed, and analysed horse SVs using a multi-assembly approach. Our findings suggest that SVs with the accumulation of mammalian-wide interspersed repeats related to long interspersed nuclear elements might be a horse-specific mechanism to modulate genome-wide gene regulatory networks. We found that olfactory receptors were commonly loss and accumulated deleterious mutations, but no purge of deleterious mutations occurred during horse domestication. We examined the potential effects of SVs on the spatial structure of chromatin via topologically associating domains (TADs). Breed-specific TADs were significantly enriched by breed-specific SVs. We identified 4199 unique breakpoint-resolved novel insertions across all chromosomes that account for 2.84 Mb sequences missing from the reference genome. Several novel insertions might have potential functional consequences, as 519 appeared to reside within 449 gene bodies. These genes are primarily involved in pathogen recognition, innate immune responses and drug metabolism. Moreover, 37 diverse horses were resequenced. Combining this with public data, we analysed 97 horses through a comparative population genomics approach to identify the genetic basis underlying breed characteristics using Thoroughbreds as a case study. We provide new scientific evidence for horse domestication, an understanding of the genetic mechanism underlying the phenotypic evolution of horses, and a comprehensive genetic variation resource for further genetic studies of horses.

The pig pangenome provides insights into the roles of coding structural variations in genetic diversity and adaptation

Structural variations have emerged as an important driving force for genome evolution and phenotypic variation in various organisms, yet their contributions to genetic diversity and adaptation in domesticated animals remain largely unknown. Here we constructed a pangenome based on 250 sequenced individuals from 32 pig breeds in Eurasia and systematically characterized coding sequence presence/absence variations (PAVs) within pigs. We identified 308.3-Mb nonreference sequences and 3438 novel genes absent from the current reference genome. Gene PAV analysis showed that 16.8% of the genes in the pangene catalog undergo PAV. A number of newly identified dispensable genes showed close associations with adaptation. For instance, several novel swine leukocyte antigen (SLA) genes discovered in nonreference sequences potentially participate in immune responses to productive and respiratory syndrome virus (PRRSV) infection. We delineated previously unidentified features of the pig mobilome that contained 490,480 transposable element insertion polymorphisms (TIPs) resulting from recent mobilization of 970 TE families, and investigated their population dynamics along with influences on population differentiation and gene expression. In addition, several candidate adaptive TE insertions were detected to be co-opted into genes responsible for responses to hypoxia, skeletal development, regulation of heart contraction, and neuronal cell development, likely contributing to local adaptation of Tibetan wild boars. These findings enhance our understanding on hidden layers of the genetic diversity in pigs and provide novel insights into the role of SVs in the evolutionary adaptation of mammals.

Analysis of polygenic selection in purebred and crossbred pig genomes using generation proxy selection mapping

BACKGROUND

Artificial selection on quantitative traits using breeding values and selection indices in commercial livestock breeding populations causes changes in allele frequency over time at hundreds or thousands of causal loci and the surrounding genomic regions. In population genetics, this type of selection is called polygenic selection. Researchers and managers of pig breeding programs are motivated to understand the genetic basis of phenotypic diversity across genetic lines, breeds, and populations using selection mapping analyses. Here, we applied generation proxy selection mapping (GPSM), a genome-wide association analysis of single nucleotide polymorphism (SNP) genotypes (38,294-46,458 markers) of birth date, in four pig populations (15,457, 15,772, 16,595 and 8447 pigs per population) to identify loci responding to artificial selection over a period of five to ten years. Gene-drop simulation analyses were conducted to provide context for the GPSM results. Selected loci within and across each population of pigs were compared in the context of swine breeding objectives.

RESULTS

The GPSM identified 49 to 854 loci as under selection (Q-values less than 0.10) across 15 subsets of pigs based on combinations of populations. The number of significant associations increased when data were pooled across populations. In addition, several significant associations were identified in more than one population. These results indicate concurrent selection objectives, similar genetic architectures, and shared causal variants responding to selection across these pig populations. Negligible error rates (less than or equal to 0.02%) of false-positive associations were found when testing GPSM on gene-drop simulated genotypes, suggesting that GPSM distinguishes selection from random genetic drift in actual pig populations.

CONCLUSIONS

This work confirms the efficacy and the negligible error rates of the GPSM method in detecting selected loci in commercial pig populations. Our results suggest shared selection objectives and genetic architectures across swine populations. The identified polygenic selection highlights loci that are important to swine production.

AgAnimalGenomes: browsers for viewing and manually annotating farm animal genomes

Current genome sequencing technologies have made it possible to generate highly contiguous genome assemblies for non-model animal species. Despite advances in genome assembly methods, there is still room for improvement in the delineation of specific gene features in the genomes. Here we present genome visualization and annotation tools to support seven livestock species (bovine, chicken, goat, horse, pig, sheep, and water buffalo), available in a new resource called AgAnimalGenomes. In addition to supporting the manual refinement of gene models, these browsers provide visualization tracks for hundreds of RNAseq experiments, as well as data generated by the Functional Annotation of Animal Genomes (FAANG) Consortium. For species with predicted gene sets from both Ensembl and RefSeq, the browsers provide special tracks showing the thousands of protein-coding genes that disagree across the two gene sources, serving as a valuable resource to alert researchers to gene model issues that may affect data interpretation. We describe the data and search methods available in the new genome browsers and how to use the provided tools to edit and create new gene models.

Young SINEs in pig genomes impact gene regulation, genetic diversity, and complex traits

Transposable elements (TEs) are a major source of genetic polymorphisms and play a role in chromatin architecture, gene regulatory networks, and genomic evolution. However, their functional role in pigs and contributions to complex traits are largely unknown. We created a catalog of TEs (n = 3,087,929) in pigs and found that young SINEs were predominantly silenced by histone modifications, DNA methylation, and decreased accessibility. However, some transcripts from active young SINEs showed high tissue-specificity, as confirmed by analyzing 3570 RNA-seq samples. We also detected 211,067 dimorphic SINEs in 374 individuals, including 340 population-specific ones associated with local adaptation. Mapping these dimorphic SINEs to genome-wide associations of 97 complex traits in pigs, we found 54 candidate genes (e.g., ANK2 and VRTN) that might be mediated by TEs. Our findings highlight the important roles of young SINEs and provide a supplement for genotype-to-phenotype associations and modern breeding in pigs.

Population analysis reveals the roles of DNA methylation in tomato domestication and metabolic diversity

DNA methylation is an important epigenetic marker, yet its diversity and consequences in tomato breeding at the population level are largely unknown. We performed whole-genome bisulfite sequencing (WGBS), RNA sequencing, and metabolic profiling on a population comprising wild tomatoes, landraces, and cultivars. A total of 8,375 differentially methylated regions (DMRs) were identified, with methylation levels progressively decreasing from domestication to improvement. We found that over 20% of DMRs overlapped with selective sweeps. Moreover, more than 80% of DMRs in tomato were not significantly associated with single-nucleotide polymorphisms (SNPs), and DMRs had strong linkages with adjacent SNPs. We additionally profiled 339 metabolites from 364 diverse accessions and further performed a metabolic association study based on SNPs and DMRs. We detected 971 and 711 large-effect loci via SNP and DMR markers, respectively. Combined with multi-omics, we identified 13 candidate genes and updated the polyphenol biosynthetic pathway. Our results showed that DNA methylation variants could complement SNP profiling of metabolite diversity. Our study thus provides a DNA methylome map across diverse accessions and suggests that DNA methylation variation can be the genetic basis of metabolic diversity in plants.

Taming transposable elements in livestock and poultry: a review of their roles and applications

Livestock and poultry play a significant role in human nutrition by converting agricultural by-products into high-quality proteins. To meet the growing demand for safe animal protein, genetic improvement of livestock must be done sustainably while minimizing negative environmental impacts. Transposable elements (TE) are important components of livestock and poultry genomes, contributing to their genetic diversity, chromatin states, gene regulatory networks, and complex traits of economic value. However, compared to other species, research on TE in livestock and poultry is still in its early stages. In this review, we analyze 72 studies published in the past 20 years, summarize the TE composition in livestock and poultry genomes, and focus on their potential roles in functional genomics. We also discuss bioinformatic tools and strategies for integrating multi-omics data with TE, and explore future directions, feasibility, and challenges of TE research in livestock and poultry. In addition, we suggest strategies to apply TE in basic biological research and animal breeding. Our goal is to provide a new perspective on the importance of TE in livestock and poultry genomes.

Cut site preference allows influenza A virus PA-X to discriminate between host and viral mRNAs

Many viruses block host gene expression to take over the infected cell. This process, termed host shutoff, is thought to promote viral replication by preventing antiviral responses and redirecting cellular resources to viral processes. Several viruses from divergent families accomplish host shutoff through RNA degradation by endoribonucleases. However, viruses also need to ensure expression of their own genes. The influenza A virus endoribonuclease PA-X solves this problem by sparing viral mRNAs and some host RNAs necessary for viral replication. To understand how PA-X distinguishes between RNAs, we characterized PA-X cut sites transcriptome-wide using 5' rapid amplification of complementary DNA ends coupled to high-throughput sequencing. This analysis, along with RNA structure predictions and validation experiments using reporters, shows that PA-Xs from multiple influenza strains preferentially cleave RNAs at GCUG tetramers in hairpin loops. Importantly, GCUG tetramers are enriched in the human but not the influenza transcriptome. Moreover, optimal PA-X cut sites inserted in the influenza A virus genome are quickly selected against during viral replication in cells. This finding suggests that PA-X evolved these cleavage characteristics to preferentially target host over viral mRNAs in a manner reminiscent of cellular self versus non-self discrimination.

A review of the pangenome: how it affects our understanding of genomic variation, selection and breeding in domestic animals?

As large-scale genomic studies have progressed, it has been revealed that a single reference genome pattern cannot represent genetic diversity at the species level. While domestic animals tend to have complex routes of origin and migration, suggesting a possible omission of some population-specific sequences in the current reference genome. Conversely, the pangenome is a collection of all DNA sequences of a species that contains sequences shared by all individuals (core genome) and is also able to display sequence information unique to each individual (variable genome). The progress of pangenome research in humans, plants and domestic animals has proved that the missing genetic components and the identification of large structural variants (SVs) can be explored through pangenomic studies. Many individual specific sequences have been shown to be related to biological adaptability, phenotype and important economic traits. The maturity of technologies and methods such as third-generation sequencing, Telomere-to-telomere genomes, graphic genomes, and reference-free assembly will further promote the development of pangenome. In the future, pangenome combined with long-read data and multi-omics will help to resolve large SVs and their relationship with the main economic traits of interest in domesticated animals, providing better insights into animal domestication, evolution and breeding. In this review, we mainly discuss how pangenome analysis reveals genetic variations in domestic animals (sheep, cattle, pigs, chickens) and their impacts on phenotypes and how this can contribute to the understanding of species diversity. Additionally, we also go through potential issues and the future perspectives of pangenome research in livestock and poultry.

Genome-Wide Association Study on Reproductive Traits Using Imputation-Based Whole-Genome Sequence Data in Yorkshire Pigs

Reproductive traits have a key impact on production efficiency in the pig industry. It is necessary to identify the genetic structure of potential genes that influence reproductive traits. In this study, a genome-wide association study (GWAS) based on chip and imputed data of five reproductive traits, namely, total number born (TNB), number born alive (NBA), litter birth weight (LBW), gestation length (GL), and number of weaned (NW), was performed in Yorkshire pigs. In total, 272 of 2844 pigs with reproductive records were genotyped using KPS Porcine Breeding SNP Chips, and then chip data were imputed to sequencing data using two online software programs: the Pig Haplotype Reference Panel (PHARP v2) and Swine Imputation Server (SWIM 1.0). After quality control, we performed GWAS based on chip data and the two different imputation databases by using fixed and random model circulating probability unification (FarmCPU) models. We discovered 71 genome-wide significant SNPs and 25 potential candidate genes (e.g., SMAD4, RPS6KA2, CAMK2A, NDST1, and ADCY5). Functional enrichment analysis revealed that these genes are mainly enriched in the calcium signaling pathway, ovarian steroidogenesis, and GnRH signaling pathways. In conclusion, our results help to clarify the genetic basis of porcine reproductive traits and provide molecular markers for genomic selection in pig breeding.

Pangenome obtained by long-read sequencing of 11 genomes reveal hidden functional structural variants in pigs

Long-read sequencing (LRS) facilitates both the genome assembly and the discovery of structural variants (SVs). Here, we built a graph-based pig pangenome by incorporating 11 LRS genomes with an average of 94.01% BUSCO completeness score, revealing 206-Mb novel sequences. We discovered 183,352 nonredundant SVs (63% novel), representing 12.12% of the reference genome. By genotyping SVs in an additional 196 short-read sequencing samples, we identified thousands of population stratified SVs. Particularly, we detected 7,568 Tibetan specific SVs, some of which demonstrate significant population differentiation between Tibetan and low-altitude pigs, which might be associated with the high-altitude hypoxia adaptation in Tibetan pigs. Further integrating functional genomic data, the most promising candidate genes within the SVs that might contribute to the high-altitude hypoxia adaptation were discovered. Overall, our study generates a benchmark pangenome resource for illustrating the important roles of SVs in adaptive evolution, domestication, and genetic improvement of agronomic traits in pigs.

Integrative Analysis of Blood Transcriptomics and Metabolomics Reveals Molecular Regulation of Backfat Thickness in Qinchuan Cattle

A crucial goal of reducing backfat thickness (BFT) is to indirectly improve feed conversion efficiency. This phenotype has been reported in certain papers; however, the molecular mechanism has yet to be fully revealed. Two extreme BFT groups, consisting of four Qinchuan cattle, were chosen for this study. We performed metabolite and transcriptome analyses of blood from cattle with a high BFT (H-BFT with average = 1.19) and from those with a low BFT (L-BFT with average = 0.39). In total, 1106 differentially expressed genes (DEGs) and 86 differentially expressed metabolites (DEMs) were identified in the extreme trait. In addition, serum ceramide was strongly correlated with BFT and could be used as a potential biomarker. Moreover, the most notable finding was that the functional genes (SMPD3 and CERS1) and metabolite (sphingosine 1-phosphate (S1P)) were filtered out and significantly enriched in the processes related to the sphingolipid metabolism. This investigation contributed to a better understanding of the subcutaneous fat depots in cattle. In general, our results indicated that the sphingolipid metabolism, involving major metabolites (serum ceramide and S1P) and key genes (SMPD3 and CERS1), could regulate BFT through blood circulation.

Distinct traces of mixed ancestry in western commercial pig genomes following gene flow from Chinese indigenous breeds

Studying gene flow between different livestock breeds will benefit the discovery of genes related to production traits and provide insight into human historical breeding. Chinese pigs have played an indispensable role in the breeding of Western commercial pigs. However, the differences in the timing and volume of the contribution of pigs from different Chinese regions to Western pigs are not yet apparent. In this paper, we combine the whole-genome sequencing data of 592 pigs from different studies and illustrate patterns of gene flow from Chinese pigs into Western commercial pigs. We describe introgression patterns from four distinct Chinese indigenous groups into five Western commercial groups. There were considerable differences in the number and length of the putative introgressed segments from Chinese pig groups that contributed to Western commercial pig breeds. The contribution of pigs from different Chinese geographical locations to a given western commercial breed varied more than that from a specific Chinese pig group to different Western commercial breeds, implying admixture within Europe after introgression. Within different Western commercial lines from the same breed, the introgression patterns from a given Chinese pig group seemed highly conserved, suggesting that introgression of Chinese pigs into Western commercial pig breeds mainly occurred at an early stage of breed formation. Finally, based on analyses of introgression signals, allele frequencies, and selection footprints, we identified a ∼2.65 Mb Chinese-derived haplotype under selection in Duroc pigs (CHR14: 95.68-98.33 Mb). Functional and phenotypic studies demonstrate that this PRKG1 haplotype is related to backfat and loin depth in Duroc pigs. Overall, we demonstrate that the introgression history of domestic pigs is complex and that Western commercial pigs contain distinct traces of mixed ancestry, likely derived from various Chinese pig breeds.

Detection of Candidate Genes Associated with Fecundity through Genome-Wide Selection Signatures of Katahdin Ewes

One of the strategies to genetically improve reproductive traits, despite their low inheritability, has been the identification of candidate genes. Therefore, the objective of this study was to detect candidate genes associated with fecundity through the fixation index (F_ST) and runs of homozygosity (ROH) of selection signatures in Katahdin ewes. Productive and reproductive records from three years were used and the genotypes (OvineSNP50K) of 48 Katahdin ewes. Two groups of ewes were identified to carry out the genetic comparison: with high fecundity (1.3 ± 0.03) and with low fecundity (1.1 ± 0.06). This study shows for the first time evidence of the influence of the CNOT11, GLUD1, GRID1, MAPK8, and CCL28 genes in the fecundity of Katahdin ewes; in addition, new candidate genes were detected for fecundity that were not reported previously in ewes but that were detected for other species: ANK2 (sow), ARHGAP22 (cow and buffalo cow), GHITM (cow), HERC6 (cow), DPF2 (cow), and TRNAC-GCA (buffalo cow, bull). These new candidate genes in ewes seem to have a high expression in reproduction. Therefore, future studies are needed focused on describing the physiological basis of changes in the reproductive behavior influenced by these genes.

LncPLAAT3-AS Regulates PLAAT3-Mediated Adipocyte Differentiation and Lipogenesis in Pigs through miR-503-5p

Animal fat deposition has a significant impact on meat flavor and texture. However, the molecular mechanisms of fat deposition are not well understood. LncPLAAT3-AS is a naturally occurring transcript that is abundant in porcine adipose tissue. Here, we focus on the regulatory role of lncPLAAT3-AS in promoting preadipocyte proliferation and adipocyte differentiation. By overexpressing or repressing lncPLAAT3 expression, we found that lncPLAAT3-AS promoted the transcription of its host gene PLAAT3, a regulator of adipocyte differentiation. In addition, we predicted the region of lncPLAAT3-AS that binds to miR-503-5p and showed by dual luciferase assay that lncPLAAT3-AS acts as a sponge to absorb miR-503-5p. Interestingly, miR-503-5p also targets and represses PLAAT3 expression and helps regulate porcine preadipocyte proliferation and differentiation. Taken together, these results show that lncPLAAT3-AS upregulates PLAAT3 expression by absorbing miR-503-5p, suggesting a potential regulatory mechanism based on competing endogenous RNAs. Finally, we explored lncPLAAT3-AS and PLAAT3 expression in adipose tissue and found that both molecules were expressed at significantly higher levels in fatty pig breeds compared to lean pig breeds. In summary, we identified the mechanism by which lncPLAAT3-AS regulates porcine preadipocyte proliferation and differentiation, contributing to our understanding of the molecular mechanisms of lipid deposition in pigs.

Whole-genome sequence association study identifies cyclin dependent kinase 8 as a key gene for the number of mummified piglets

OBJECTIVE

Pigs, an ideal biomedical model for human diseases, suffer from about 50% early embryonic and fetal death, a major cause of fertility loss worldwide. However, identifying the causal variant remains a huge challenge. This study aimed to detect single nucleotide polymorphisms (SNPs) and candidate genes for the number of mummified (NM) piglets using the imputed whole-genome sequence (WGS) and validate the potential candidate genes.

METHODS

The imputed WGS was introduced from genotyping-by-sequencing (GBS) using a multi-breed reference population. We performed genome-wide association studies (GWAS) for NM piglets at birth from a Landrace pig populatiGWAS peak located on SSC11: 0.10 to 7.11 Mbp (Top SNP, SSC11:1,889,658 bp; p = 9.98E-13) was identified in cyclin dependent kinase on. A total of 300 Landrace pigs were genotyped by GBS. The whole-genome variants were imputed, and 4,252,858 SNPs were obtained. Various molecular experiments were conducted to determine how the genes affected NM in pigs.

RESULTS

A strong GWAS peak located on SSC11: 0.10 to 7.11 Mbp (Top SNP, SSC11:1,889,658 bp; p = 9.98E-13) was identified in cyclin dependent kinase 8 (CDK8) gene, which plays a crucial role in embryonic retardation and lethality. Based on the molecular experiments, we found that Y-box binding protein 1 (YBX1) was a crucial transcription factor for CDK8, which mediated the effect of CDK8 in the proliferation of porcine ovarian granulosa cells via transforming growth factor beta/small mother against decapentaplegic signaling pathway, and, as a consequence, affected embryo quality, indicating that this pathway may be contributing to mummified fetal in pigs.

CONCLUSION

A powerful imputation-based association study was performed to identify genes associated with NM in pigs. CDK8 was suggested as a functional gene for the proliferation of porcine ovarian granulosa cells, but further studies are required to determine causative mutations and the effect of loci on NM in pigs.

An effect of large-scale deletions and duplications on transcript expression

Since copy number variants (CNVs) have been recognized as an important source of genetic and transcriptomic variation, we aimed to characterize the impact of CNVs located within coding, intergenic, upstream, and downstream gene regions on the expression of transcripts. Regions in which deletions occurred most often were introns, while duplications in coding regions. The transcript expression was lower for deleted coding (P = 0.008) and intronic regions (P = 1.355 × 10^-10), but it was not changed in the case of upstream and downstream gene regions (P = 0.085). Moreover, the expression was decreased if duplication occurred in the coding region (P = 8.318 × 10^-5). Furthermore, a negative correlation (r =  - 0.27) between transcript length and its expression was observed. The correlation between the percent of deleted/duplicated transcript and transcript expression level was not significant for all concerned genomic regions in five out of six animals. The exceptions were deletions in coding regions (P = 0.004) and duplications in introns (P = 0.01) in one individual. CNVs in coding (deletions, duplications) and intronic (deletions) regions are important modulators of transcripts by reducing their expression level. We hypothesize that deletions imply severe consequences by interrupting genes. The negative correlation between the size of the transcript and its expression level found in this study is consistent with the hypothesis that selection favours shorter introns and a moderate number of exons in highly expressed genes. This may explain the transcript expression reduction by duplications. We did not find the correlation between the size of deletions/duplications and transcript expression level suggesting that expression is modulated by CNVs regardless of their size.

Introgression of a Complex Genomic Structural Variation Causes Hybrid Male Sterility in GJ Rice (Oryza sativa L.) Subspecies

Hybrids between different subspecies of rice Oryza sativa L. commonly show hybrid sterility. Here we show that a widely planted commercial japonica/GJ variety, DHX2, exhibited hybrid sterility when crossing with other GJ varieties. Using the high-quality genome assembly, we identified three copies of the Sc gene in DHX2, whereas Nipponbare (Nip) had only one copy of Sc. Knocking out the extra copies of Sc in DHX2 significantly improved the pollen fertility of the F1 plant of DHX2/Nip cross. The population structure analysis revealed that a slight introgression from Basmati1 might occur in the genome of DHX2. We demonstrated that both DHX2 and Basmati1 harbored three copies of Sc. Moreover, the introgression of GS3 and BADH2/fgr from Basmati1 confers the slender and fragrance grain of DHX2. These results add to our understanding of the hybrid sterility of inter-subspecies and intra-subspecies and may provide a novel strategy for hybrid breeding.

Integration of multi-omics data reveals cis-regulatory variants that are associated with phenotypic differentiation of eastern from western pigs

Background

The genetic mechanisms that underlie phenotypic differentiation in breeding animals have important implications in evolutionary biology and agriculture. However, the contribution of cis-regulatory variants to pig phenotypes is poorly understood. Therefore, our aim was to elucidate the molecular mechanisms by which non-coding variants cause phenotypic differences in pigs by combining evolutionary biology analyses and functional genomics.

Results

We obtained a high-resolution phased chromosome-scale reference genome with a contig N50 of 18.03 Mb for the Luchuan pig breed (a representative eastern breed) and profiled potential selective sweeps in eastern and western pigs by resequencing the genomes of 234 pigs. Multi-tissue transcriptome and chromatin accessibility analyses of these regions suggest that tissue-specific selection pressure is mediated by promoters and distal cis-regulatory elements. Promoter variants that are associated with increased expression of the lysozyme (LYZ) gene in the small intestine might enhance the immunity of the gastrointestinal tract and roughage tolerance in pigs. In skeletal muscle, an enhancer-modulating single-nucleotide polymorphism that is associated with up-regulation of the expression of the troponin C1, slow skeletal and cardiac type (TNNC1) gene might increase the proportion of slow muscle fibers and affect meat quality.

Conclusions

Our work sheds light on the molecular mechanisms by which non-coding variants shape phenotypic differences in pigs and provides valuable resources and novel perspectives to dissect the role of gene regulatory evolution in animal domestication and breeding.

A transcriptional cross species map of pancreatic islet cells

OBJECTIVE

Pancreatic islets of Langerhans secrete hormones to regulate systemic glucose levels. Emerging evidence suggests that islet cells are functionally heterogeneous to allow a fine-tuned and efficient endocrine response to physiological changes. A precise description of the molecular basis of this heterogeneity, in particular linking animal models to human islets, is an important step towards identifying the factors critical for endocrine cell function in physiological and pathophysiological conditions.

METHODS

In this study, we used single-cell RNA sequencing to profile more than 50'000 endocrine cells isolated from healthy human, pig and mouse pancreatic islets and characterize transcriptional heterogeneity and evolutionary conservation of those cells across the three species. We systematically delineated endocrine cell types and α- and β-cell heterogeneity through prior knowledge- and data-driven gene sets shared across species, which altogether capture common and differential cellular properties, transcriptional dynamics and putative driving factors of state transitions.

RESULTS

We showed that global endocrine expression profiles correlate, and that critical identity and functional markers are shared between species, while only approximately 20% of cell type enriched expression is conserved. We resolved distinct human α- and β-cell states that form continuous transcriptional landscapes. These states differentially activate maturation and hormone secretion programs, which are related to regulatory hormone receptor expression, signaling pathways and different types of cellular stress responses. Finally, we mapped mouse and pig cells to the human reference and observed that the spectrum of human α- and β-cell heterogeneity and aspects of such functional gene expression are better recapitulated in the pig than mouse data.

CONCLUSIONS

Here, we provide a high-resolution transcriptional map of healthy human islet cells and their murine and porcine counterparts, which is easily queryable via an online interface. This comprehensive resource informs future efforts that focus on pancreatic endocrine function, failure and regeneration, and enables to assess molecular conservation in islet biology across species for translational purposes.

Expression and Genetic Effects of GLI Pathogenesis-Related 1 Gene on Backfat Thickness in Pigs

Backfat thickness (BFT) is an important carcass composition trait and regarded as a breeding focus. Our initial transcriptome analysis of pig BFT identified GLI pathogenesis-related 1 (GLIPR1) as one of the promising candidate genes. This study was conducted to identify the expression profiles, polymorphisms, and genetic effects of the GLIPR1 gene on BFT in pigs. The expression of the GLIPR1 gene existed in every detected tissue, and there was a significantly higher expression in spleen and adipose tissue than others (p < 0.05). At the different ages of pig, the expression of the GLIPR1 gene was low at an early age, increased with growth, and reached the highest level at 180 days. Genetic polymorphism analysis was detected in 553 individuals of the Large White × Minzhu F2 population. Four SNPs in the promoter significantly associated with 6−7 rib BFT (p < 0.05) were predicted to alter the transcription factor binding sites (TFBS), and the mutations of g.38758089 T>G and g.38758114 G>C were predicted to change the TFs associated with the regulation of adipogenesis. Haplotypes were formed by the detected SNPs, and one block showed a strong association with BFT (p < 0.05). In summary, our results indicate that the expression profiles and genetic variants of GLIPR1 affected the BFT of pigs. To our knowledge, this study is the first to demonstrate the biological function and genetic effects of the GLIPR1 gene on the BFT of pig and provide genetic markers to optimize breeding for BFT in pigs.

A missense mutation in ISPD contributes to maintain muscle fiber stability

Background

Results

Conclusion

Loss of Monoallelic Expression of IGF2 in the Adult Liver Via Alternative Promoter Usage and Chromatin Reorganization

In mammals, genomic imprinting operates via gene silencing mechanisms. Although conservation of the imprinting mechanism at the H19/IGF2 locus has been generally described in pigs, tissue-specific imprinting at the transcript level, monoallelic-to-biallelic conversion, and spatio-temporal chromatin reorganization remain largely uninvestigated. Here, we delineate spatially regulated imprinting of IGF2 transcripts, age-dependent hepatic mono- to biallelic conversion, and reorganization of topologically associating domains at the porcine H19/IGF2 locus for better translation to human and animal research. Whole-genome bisulfite sequencing (WGBS) and RNA sequencing (RNA-seq) of normal and parthenogenetic porcine embryos revealed the paternally hypermethylated H19 differentially methylated region and paternal expression of IGF2. Using a polymorphism-based approach and omics datasets from chromatin immunoprecipitation sequencing (ChIP–seq), whole-genome sequencing (WGS), RNA-seq, and Hi-C, regulation of IGF2 during development was analyzed. Regulatory elements in the liver were distinguished from those in the muscle where the porcine IGF2 transcript was monoallelically expressed. The IGF2 transcript from the liver was biallelically expressed at later developmental stages in both pigs and humans. Chromatin interaction was less frequent in the adult liver compared to the fetal liver and skeletal muscle. The duration of genomic imprinting effects within the H19/IGF2 locus might be reduced in the liver with biallelic conversion through alternative promoter usage and chromatin remodeling. Our integrative omics analyses of genome, epigenome, and transcriptome provided a comprehensive view of imprinting status at the H19/IGF2 cluster.

Detection of non-reference porcine endogenous retrovirus loci in the Vietnamese native pig genome

The Vietnamese native pig (VnP)-a porcine breed with a small body-has proven suitable as a biomedical animal model. Here, we demonstrate that, compared to other breeds, VnPs have fewer copies of porcine endogenous retroviruses (PERVs), which pose a risk for xenotransplantation of pig organs to humans. More specifically, we sought to characterize non-reference PERVs (nrPERVs) that were previously unidentified in the reference genome. To this end, we used whole-genome sequencing data to identify nrPERV loci with long terminal repeat (LTR) sequences in VnPs. RetroSeq was used to estimate nrPERV loci based on the most current porcine reference genome (Sscrofa11.1). LTRs were detected using de novo sequencing read assembly near the loci containing the target site duplication sequences in the inferred regions. A total of 21 non-reference LTR loci were identified and separated into two subtypes based on phylogenetic analysis. Moreover, PERVs within the detected LTR loci were identified, the presence of which was confirmed using conventional PCR and Sanger sequencing. These novel loci represent previously unknown PERVs as they have not been identified in the porcine reference genome. Thus, our RetroSeq method accurately detects novel PERV loci, and can be applied for development of a useful biomedical model.

Genome‐wide analysis of hybridization in wild boar populations reveals adaptive introgression from domestic pig

The admixture of domestic pig into French wild boar populations has been monitored since the 1980s thanks to the existence of a cytogenetic difference between the two sub‐species. The number of chromosomes is 2n = 36 in wild boar and 2n = 38 in pig, respectively. This difference makes it possible to assign the “hybrid” status to wild boar individuals controlled with 37 or 38 chromosomes. However, it does not make it possible to determine the timing of the hybridization(s), nor to guarantee the absence of domestic admixture in an animal with 2n = 36 chromosomes. In order to analyze hybridization in greater detail and to avoid the inherent limitations of the cytogenetic approach, 362 wild boars (WB) recently collected in different French geographical areas and in different environments (farms, free ranging in protected or unprotected areas, animals with 2n = 36, 37 or 38 chromosomes) were genotyped on a 70K SNP chip. Principal component analyses allowed the identification of 13 “outliers” (3.6%), for which the proportion of the genome of “domestic” origin was greater than 40% (Admixture analyses). These animals were probably recent hybrids, having Asian domestic pig ancestry for most of them. For the remaining 349 animals studied, the proportion of the genome of “wild” origin varied between 83% and 100% (median: 94%). This proportion varied significantly depending on how the wild boar populations were managed. Local ancestry analyses revealed adaptive introgression from domestic pig, suggesting a critical role of genetic admixture in improving the fitness and population growth of WB. Overall, our results show that the methods used to monitor the domestic genetic contributions to wild boar populations should evolve in order to limit the level of admixture between the two gene pools.

The Visayan warty pig (Sus cebifrons) genome provides insight into chromosome evolution and sensory adaptation in pigs

It is largely unknown how mammalian genomes evolve under rapid speciation and environmental adaptation. An excellent model for understanding fast evolution is provided by the genus Sus, which diverged relatively recently and lacks post-zygotic isolation. Here, we present a high-quality reference genome of the Visayan warty pig, which is specialized to a tropical island environment. Comparing the genome sequences and chromatin contact maps of the Visayan warty pig (Sus cebifrons) and domestic pig (Sus scrofa), we characterized the dynamics of chromosomal structure evolution during Sus speciation, revealing the similar chromosome conformation as the potential biological mechanism of frequent post-divergence hybridization among Suidae. We further investigated the different signatures of adaptive selection and domestication in Visayan warty pig and domestic pig with specific emphasize on the evolution of olfactory and gustatory genes, elucidating higher olfactory diversity in Visayan warty pig and positive and relaxed evolution of bitter and fat taste receptors, respectively, in domestic pig. Our comprehensive evolutionary and comparative genome analyses provide insight into the dynamics of genomes and how these change over relative short evolutionary times, as well as how these genomic differences encode for differences in the phenotypes.

Machine-Learning Prospects for Detecting Selection Signatures Using Population Genomics Data

Natural selection has been given a lot of attention because it relates to the adaptation of populations to their environments, both biotic and abiotic. An allele is selected when it is favored by natural selection. Consequently, the favored allele increases in frequency in the population and neighboring linked variation diminishes, causing so-called selective sweeps. A high-throughput genomic sequence allows one to disentangle the evolutionary forces at play in populations. With the development of high-throughput genome sequencing technologies, it has become easier to detect these selective sweeps/selection signatures. Various methods can be used to detect selective sweeps, from simple implementations using summary statistics to complex statistical approaches. One of the important problems of these statistical models is the potential to provide inaccurate results when their assumptions are violated. The use of machine learning (ML) in population genetics has been introduced as an alternative method of detecting selection by treating the problem of detecting selection signatures as a classification problem. Since the availability of population genomics data is increasing, researchers may incorporate ML into these statistical models to infer signatures of selection with higher predictive accuracy and better resolution. This article describes how ML can be used to aid in detecting and studying natural selection patterns using population genomic data.

Genome-wide detection of genetic structure and runs of homozygosity analysis in Anhui indigenous and Western commercial pig breeds using PorcineSNP80k data

Background

Runs of homozygosity (ROH) are continuous homozygous regions typically located in the DNA sequence of diploid organisms. Identifications of ROH that lead to reduced performance can provide valuable insight into the genetic architecture of complex traits. Here, we systematically investigated the population genetic structure of five Anhui indigenous pig breeds (AHIPs), and compared them to those of five Western commercial pig breeds (WECPs). Furthermore, we examined the occurrence and distribution of ROHs in the five AHIPs and estimated the inbreeding coefficients based on the ROHs (F_ROH) and homozygosity (F_HOM). Finally, we identified genomic regions with high frequencies of ROHs and annotated candidate genes contained therein.

Results

The WECPs and AHIPs were clearly differentiated into two separate clades consistent with their geographical origins, as revealed by the population structure and principal component analysis. We identified 13,530 ROHs across all individuals, of which 4,555 and 8,975 ROHs were unique to AHIPs and WECPs, respectively. Most ROHs identified in our study were short (< 10 Mb) or medium (10–20 Mb) in length. WECPs had significantly higher numbers of short ROHs, and AHIPs generally had longer ROHs. F_ROH values were significantly lower in AHIPs than in WECPs, indicating that breed improvement and conservation programmes were successful in AHIPs. On average, F_ROH and F_HOM values were highly correlated (0.952–0.991) in AHIPs and WECPs. A total of 27 regions had a high frequency of ROHs and contained 17 key candidate genes associated with economically important traits in pigs. Among these, nine candidate genes (CCNT2, EGR2, MYL3, CDH13, PROX1, FLVCR1, SETD2, FGF18, and FGF20) found in WECPs were related to muscular and skeletal development, whereas eight candidate genes (CSN1S1, SULT1E1, TJP1, ZNF366, LIPC, MCEE, STAP1, and DUSP) found in AHIPs were associated with health, reproduction, and fatness traits.

Conclusion

Our findings provide a useful reference for the selection and assortative mating of pig breeds, laying the groundwork for future research on the population genetic structures of AHIPs, ultimately helping protect these local varieties.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08583-9.

The de novo assembly of a European wild boar genome revealed unique patterns of chromosomal structural variations and segmental duplications

The rapid progress of sequencing technology has greatly facilitated the de novo genome assembly of pig breeds. However, the assembly of the wild boar genome is still lacking, hampering our understanding of chromosomal and genomic evolution during domestication from wild boars into domestic pigs. Here, we sequenced and de novo assembled a European wild boar genome (ASM2165605v1) using the long‐range information provided by 10× Linked‐Reads sequencing. We achieved a high‐quality assembly with contig N50 of 26.09 Mb. Additionally, 1.64% of the contigs (222) with lengths from 107.65 kb to 75.36 Mb covered 90.3% of the total genome size of ASM2165605v1 (~2.5 Gb). Mapping analysis revealed that the contigs can fill 24.73% (93/376) of the gaps present in the orthologous regions of the updated pig reference genome (Sscrofa11.1). We further improved the contigs into chromosome level with a reference‐assistant scaffolding method. Using the ‘assembly‐to‐assembly’ approach, we identified intra‐chromosomal large structural variations (SVs, length >1 kb) between ASM2165605v1 and Sscrofa11.1 assemblies. Interestingly, we found that the number of SV events on the X chromosome deviated significantly from the linear models fitting autosomes (R² > 0.64, p < 0.001). Specifically, deletions and insertions were deficient on the X chromosome by 66.14 and 58.41% respectively, whereas duplications and inversions were excessive on the X chromosome by 71.96 and 107.61% respectively. We further used the large segmental duplications (SDs, >1 kb) events as a proxy to understand the large‐scale inter‐chromosomal evolution, by resolving parental‐derived relationships for SD pairs. We revealed a significant excess of SD movements from the X chromosome to autosomes (p < 0.001), consistent with the expectation of meiotic sex chromosome inactivation. Enrichment analyses indicated that the genes within derived SD copies on autosomes were significantly related to biological processes involving nervous system, lipid biosynthesis and sperm motility (p < 0.01). Together, our analyses of the de novo assembly of ASM2165605v1 provides insight into the SVs between European wild boar and domestic pig, in addition to the ongoing process of meiotic sex chromosome inactivation in driving inter‐chromosomal interaction between the sex chromosome and autosomes.

Genome-wide annotation of protein-coding genes in pig

BACKGROUND

There is a need for functional genome-wide annotation of the protein-coding genes to get a deeper understanding of mammalian biology. Here, a new annotation strategy is introduced based on dimensionality reduction and density-based clustering of whole-body co-expression patterns. This strategy has been used to explore the gene expression landscape in pig, and we present a whole-body map of all protein-coding genes in all major pig tissues and organs.

RESULTS

An open-access pig expression map ( www.rnaatlas.org ) is presented based on the expression of 350 samples across 98 well-defined pig tissues divided into 44 tissue groups. A new UMAP-based classification scheme is introduced, in which all protein-coding genes are stratified into tissue expression clusters based on body-wide expression profiles. The distribution and tissue specificity of all 22,342 protein-coding pig genes are presented.

CONCLUSIONS

Here, we present a new genome-wide annotation strategy based on dimensionality reduction and density-based clustering. A genome-wide resource of the transcriptome map across all major tissues and organs in pig is presented, and the data is available as an open-access resource ( www.rnaatlas.org ), including a comparison to the expression of human orthologs.

The genome variation and developmental transcriptome maps reveal genetic differentiation of skeletal muscle in pigs

Natural and artificial directional selections have resulted in significantly genetic and phenotypic differences across breeds in domestic animals. However, the molecular regulation of skeletal muscle diversity remains largely unknown. Here, we conducted transcriptome profiling of skeletal muscle across 27 time points, and performed whole-genome re-sequencing in Landrace (lean-type) and Tongcheng (obese-type) pigs. The transcription activity decreased with development, and the high-resolution transcriptome precisely captured the characterizations of skeletal muscle with distinct biological events in four developmental phases: Embryonic, Fetal, Neonatal, and Adult. A divergence in the developmental timing and asynchronous development between the two breeds was observed; Landrace showed a developmental lag and stronger abilities of myoblast proliferation and cell migration, whereas Tongcheng had higher ATP synthase activity in postnatal periods. The miR-24-3p driven network targeting insulin signaling pathway regulated glucose metabolism. Notably, integrated analysis suggested SATB2 and XLOC_036765 contributed to skeletal muscle diversity via regulating the myoblast migration and proliferation, respectively. Overall, our results provide insights into the molecular regulation of skeletal muscle development and diversity in mammals.

Long-read assembly of the Chinese indigenous Ningxiang pig genome and identification of genetic variations in fat metabolism among different breeds

Advances in long-read sequencing technology and genome assembly provide an opportunity to improve the pig genome and reveal the full range of structural variations (SVs) between local Chinese and European pigs. To date, little is known about the genomes of some unique Chinese indigenous breeds, such as the Ningxiang pig. Here, we report the sequencing and assembly of a highly contiguous Ningxiang pig genome (NX) via an integration of PacBio single-molecule real-time sequencing, Illumina next-generation sequencing, BioNano optical mapping and Hi-C (chromosome conformation capture) approaches. The assembled genome comprises 2.44 Gb with a contig N50 of 26.1 Mb and 418 contigs in total. These contigs are organized into 121 scaffolds with a scaffold N50 of 139.0 Mb. More than 99.1% of the assembled sequence could be localized to 19 pseudochromosomes and is annotated with 20,914 protein-coding genes and 34.04% repetitive sequences. Comparisons between the NX and European Duroc assemblies revealed many SVs in genes involved in the immune system, nervous system, lipid metabolism and environmental adaptation. The genetic variants include 47 Chinese domestic pig-specific SVs and the associated 74 genes may contribute to the differences in domestic traits compared to European pigs. Moreover, single nucleotide polymorphisms (SNPs) identified from whole genome resequencing data of 73 Chinese pigs, representing 17 geographically isolated breeds, showed their specific genetic variations, population structure and evolutionary patterns. Finally, we explore transcriptional regulation in the first intron of the MYL4 gene, as the genomic SV (281-bp deletion) in Ningxiang pig promotes its subcutaneous fat compared to European pig breeds. This work identifies a set of Asian-specific SVs and SNPs, which will be important resources for modern pig breeding and genetic conservation.