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

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.

pyRSD-CoEv: A python package for selective sweep detection and co-evolutionary gene cluster identification

Genes undergo distinct selective sweeps, and also interact and coevolve, forming the bases of complex phenotypic traits. Therefore, the identification of genes that coevolve or are under artificial selective sweeps is of great importance. However, previous computational methods have been designed for either populations of closely related breeds or individuals of distinct species. Approaches intended specifically for closely related individuals without replicate (i.e. each breed/strain is represented by only one individual) are long overdue. We present a free, powerful, open source package, pyRSD-CoEv, that allows the identification of genes undergoing coevolution and/or selection-based sweeps. pyRSD-CoEv includes two main analysis workflows for genomic variant data: (i) the identification of selective sweeps using relative homozygous single nucleotide variant density (RSD); and (ii) the identification of coevolutionary gene clusters based on correlated evolutionary rates. The python package pyRSD-CoEv is written using python 3.7 and is freely available from the github website at https://github.com/QianZiTang/pyRSD-CoEv. It runs on Linux.

Human-Mediated Admixture and Selection Shape the Diversity on the Modern Swine (Sus scrofa) Y Chromosomes

Throughout its distribution across Eurasia, domestic pig (Sus scrofa) populations have acquired differences through natural and artificial selection, and have often interbred. We resequenced 80 Eurasian pigs from nine different Asian and European breeds; we identify 42,288 reliable SNPs on the Y chromosome in a panel of 103 males, among which 96.1% are newly detected. Based on these new data, we elucidate the evolutionary history of pigs through the lens of the Y chromosome. We identify two highly divergent haplogroups: one present only in Asia and one fixed in Europe but present in some Asian populations. Analyzing the European haplotypes present in Asian populations, we find evidence of three independent waves of introgression from Europe to Asia in last 200 years, agreeing well with the literature and historical records. The diverse European lineages were brought in China by humans and left significant imprints not only on the autosomes but also on the Y chromosome of geographically and genetically distinct Chinese pig breeds. We also find a general excess of European ancestry on Y chromosomes relative to autosomes in Chinese pigs, an observation that cannot be explained solely by sex-biased migration and genetic drift. The European Y haplotype is associated with leaner meat production, and we hypothesize that the European Y chromosome increased in frequency in Chinese populations due to artificial selection. We find evidence of Y chromosomal gene flow between Sumatran wild boar and Chinese pigs. Our results demonstrate how human-mediated admixture and selection shaped the distribution of modern swine Y chromosomes.

Landscape of tissue-specific RNA Editome provides insight into co-regulated and altered gene expression in pigs (Sus-scrofa)

RNA editing generates genetic diversity in mammals by altering amino acid sequences, miRNA targeting site sequences, influencing the stability of targeted RNAs, and causing changes in gene expression. However, the extent to which RNA editing affect gene expression via modifying miRNA binding site remains unexplored. Here, we first profiled the dynamic A-to-I RNA editome across tissues of Duroc and Luchuan pigs. The RNA editing events at the miRNA binding sites were generated. The biological function of the differentially edited gene in skeletal muscle was further characterized in pig muscle-derived satellite cells. RNA editome analysis revealed a total of 171,909 A-to-I RNA editing sites (RESs), and examination of its features showed that these A-to-I editing sites were mainly located in SINE retrotransposons PRE-1/Pre0_SS element. Analysis of differentially edited sites (DESs) revealed a total of 4,552 DESs across tissues between Duroc and Luchuan pigs, and functional category enrichment analysis of differentially edited gene (DEG) sets highlighted a significant association and enrichment of tissue-developmental pathways including TGF-beta, PI3K-Akt, AMPK, and Wnt signaling pathways. Moreover, we found that RNA editing events at the miRNA binding sites in the 3'-UTR of HSPA12B mRNA could prevent the miRNA-mediated mRNA downregulation of HSPA12B in the muscle-derived satellite (MDS) cell, consistent with the results obtained from the Luchuan skeletal muscle. This study represents the most systematic attempt to characterize the significance of RNA editing in regulating gene expression, particularly in skeletal muscle, constituting a new layer of regulation to understand the genetic mechanisms behind phenotype variance in animals.Abbreviations: A-to-I: Adenosine-to-inosine; ADAR: Adenosine deaminase acting on RNA; RES: RNA editing site; DEG: Differentially edited gene; DES: Differentially edited site; FDR: False discovery rate; GO: Gene Ontology; KEGG: Kyoto Encyclopaedia of Genes and Genomes; MDS cell: musclederived satellite cell; RPKM: Reads per kilobase of exon model in a gene per million mapped reads; UTR: Untranslated coding regions.

Breed-specific reference sequence optimized mapping accuracy of NGS analyses for pigs

Background

Reference sequences play a vital role in next-generation sequencing (NGS), impacting mapping quality during genome analyses. However, reference genomes usually do not represent the full range of genetic diversity of a species as a result of geographical divergence and independent demographic events of different populations. For the mitochondrial genome (mitogenome), which occurs in high copy numbers in cells and is strictly maternally inherited, an optimal reference sequence has the potential to make mitogenome alignment both more accurate and more efficient. In this study, we used three different types of reference sequences for mitogenome mapping, i.e., the commonly used reference sequence (CU-ref), the breed-specific reference sequence (BS-ref) and the sample-specific reference sequence (SS-ref), respectively, and compared the accuracy of mitogenome alignment and SNP calling among them, for the purpose of proposing the optimal reference sequence for mitochondrial DNA (mtDNA) analyses of specific populations

Results

Four pigs, representing three different breeds, were high-throughput sequenced, subsequently mapping reads to the reference sequences mentioned above, resulting in a largest mapping ratio and a deepest coverage without increased running time when aligning reads to a BS-ref. Next, single nucleotide polymorphism (SNP) calling was carried out by 18 detection strategies with the three tools SAMtools, VarScan and GATK with different parameters, using the bam results mapping to BS-ref. The results showed that all eighteen strategies achieved the same high specificity and sensitivity, which suggested a high accuracy of mitogenome alignment by the BS-ref because of a low requirement for SNP calling tools and parameter choices.

Conclusions

This study showed that different reference sequences representing different genetic relationships to sample reads influenced mitogenome alignment, with the breed-specific reference sequences being optimal for mitogenome analyses, which provides a refined processing perspective for NGS data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08030-1.

Antiviral Activity of Canine RIG-I against Canine Influenza Virus and Interactions between Canine RIG-I and CIV

RIG-I functions as a virus sensor that induces a cellular antiviral response. Although it has been investigated in other species, there have been no further studies to date on canine RIG-I against canine influenza virus (CIV). In the present study, we cloned the RIG-I gene of beagle dogs and characterized its expression, subcellular localization, antiviral response, and interactions with CIV proteins. RIG-I was highly expressed and mainly localized in the cytoplasm, with low levels detected in the nucleus. The results revealed that overexpression of the CARD domain of RIG-I and knockdown of RIG-I showed its ability to activate the RLR pathway and induced the expression of downstream interferon-stimulated genes. Moreover, overexpression of canine RIG-I suppressed the replication of CIV. The association between RIG-I and CIV was evaluated with the luciferase assay and by indirect immunofluorescence and bimolecular fluorescence complementation analyses. The results showed that CIV nonstructural protein 1 (NS1) can strongly suppress the RIG-I–mediated innate immune response, and the novel interactions between CIV matrix proteins (M1 and M2) and canine RIG-I were disclosed. These findings provide a basis for investigating the antiviral mechanism of canine RIG-I against CIV, which can lead to effective strategies for preventing CIV infection in dogs.

PRE-1 Revealed Previous Unknown Introgression Events in Eurasian Boars during the Middle Pleistocene

Introgression events and population admixture occurred among Sus species across the Eurasian mainland in the Middle Pleistocene, which reflects the local adaption of different populations and contributes to evolutionary novelty. Previous findings on these population introgressions were largely based on extensive genome-wide single-nucleotide polymorphism information, ignoring structural variants (SVs) as an important alternative resource of genetic variations. Here, we profiled the genome-wide SVs and explored the formation of pattern-related SVs, indicating that PRE1-SS is a recently active subfamily that was strongly associated with introgression events in multiple Asian and European pig populations. As reflected by the three different combination haplotypes from two specific patterns and known phylogenetic relationships in Eurasian boars, we identified the Asian Northern wild pigs as having experienced introgression from European wild boars around 0.5–0.2 Ma and having received latitude-related selection. During further exploration of the influence of pattern-related SVs on gene functions, we found substantial sequence changes in 199 intron regions of 54 genes and 3 exon regions of 3 genes (HDX, TRO, and SMIM1), implying that the pattern-related SVs were highly related to positive selection and adaption of pigs. Our findings revealed novel introgression events in Eurasian wild boars, providing a timeline of population admixture and divergence across the Eurasian mainland in the Middle Pleistocene.

The chicken pan-genome reveals gene content variation and a promoter region deletion in IGF2BP1 affecting body size

Domestication and breeding have reshaped the genomic architecture of chicken, but the retention and loss of genomic elements during these evolutionary processes remain unclear. We present the first chicken pan-genome constructed using 664 individuals, which identified an additional ∼66.5 Mb sequences that are absent from the reference genome (GRCg6a). The constructed pan-genome encoded 20,491 predicated protein-coding genes, of which higher expression level are observed in conserved genes relative to dispensable genes. Presence/absence variation (PAV) analyses demonstrated that gene PAV in chicken was shaped by selection, genetic drift, and hybridization. PAV-based GWAS identified numerous candidate mutations related to growth, carcass composition, meat quality, or physiological traits. Among them, a deletion in the promoter region of IGF2BP1 affecting chicken body size is reported, which is supported by functional studies and extra samples. This is the first time to report the causal variant of chicken body size QTL located at chromosome 27 which was repeatedly reported. Therefore, the chicken pan-genome is a useful resource for biological discovery and breeding. It improves our understanding of chicken genome diversity and provides materials to unveil the evolution history of chicken domestication.

Eight soybean reference genome resources from varying latitudes and agronomic traits

Comparative analysis of multiple reference genomes representing diverse genetic backgrounds is critical for understanding the role of key alleles important in domestication and genetic breeding of important crops such as soybean. To enrich the genetic resources for soybean, we describe the generation, technical assessment, and preliminary genomic variation analysis of eight de novo reference-grade soybean genome assemblies from wild and cultivated accessions. These resources represent soybeans cultured at different latitudes and exhibiting different agronomical traits. Of these eight soybeans, five are from new accessions that have not been sequenced before. We demonstrate the usage of these genomes to identify small and large genomic variations affecting known genes as well as screening for genic PAV regions for identifying candidates for further functional studies.

Novel functional sequences uncovered through a bovine multiassembly graph

Many genomic analyses start by aligning sequencing reads to a linear reference genome. However, linear reference genomes are imperfect, lacking millions of bases of unknown relevance and are unable to reflect the genetic diversity of populations. This makes reference-guided methods susceptible to reference-allele bias. To overcome such limitations, we build a pangenome from six reference-quality assemblies from taurine and indicine cattle as well as yak. The pangenome contains an additional 70,329,827 bases compared to the Bos taurus reference genome. Our multiassembly approach reveals 30 and 10.1 million bases private to yak and indicine cattle, respectively, and between 3.3 and 4.4 million bases unique to each taurine assembly. Utilizing transcriptomes from 56 cattle, we show that these nonreference sequences encode transcripts that hitherto remained undetected from the B. taurus reference genome. We uncover genes, primarily encoding proteins contributing to immune response and pathogen-mediated immunomodulation, differentially expressed between Mycobacterium bovis-infected and noninfected cattle that are also undetectable in the B. taurus reference genome. Using whole-genome sequencing data of cattle from five breeds, we show that reads which were previously misaligned against the Bos taurus reference genome now align accurately to the pangenome sequences. This enables us to discover 83,250 polymorphic sites that segregate within and between breeds of cattle and capture genetic differentiation across breeds. Our work makes a so-far unused source of variation amenable to genetic investigations and provides methods and a framework for establishing and exploiting a more diverse reference genome.

Accessing Livestock Resources in Ensembl

Genome assembly is cheaper, more accurate and more automated than it has ever been. This is due to a combination of more cost-efficient chemistries, new sequencing technologies and better algorithms. The livestock community has been at the forefront of this new wave of genome assembly, generating some of the highest quality vertebrate genome sequences. Ensembl's goal is to add functional and comparative annotation to these genomes, through our gene annotation, genomic alignments, gene trees, regulatory, and variation data. We run computationally complex analyses in a high throughput and consistent manner to help accelerate downstream science. Our livestock resources are continuously growing in both breadth and depth. We annotate reference genome assemblies for newly sequenced species and regularly update annotation for existing genomes. We are the only major resource to support the annotation of breeds and other non-reference assemblies. We currently provide resources for 13 pig breeds, maternal and paternal haplotypes for hybrid cattle and various other non-reference or wild type assemblies for livestock species. Here, we describe the livestock data present in Ensembl and provide protocols for how to view data in our genome browser, download via it our FTP site, manipulate it via our tools and interact with it programmatically via our REST API.

Classical, Molecular, and Genomic Cytogenetics of the Pig, a Clinical Perspective

The chromosomes of the domestic pig (Sus scrofa domesticus) are known to be prone to reciprocal chromosome translocations and other balanced chromosome rearrangements with concomitant fertility impairment of carriers. In response to the remarkable prevalence of chromosome rearrangements in swine herds, clinical cytogenetics laboratories have been established in several countries in order to screen young boars for chromosome rearrangements prior to service. At present, clinical cytogenetics laboratories typically apply classical cytogenetics techniques such as giemsa-trypsin (GTG)-banding to produce high-quality karyotypes and reveal large-scale chromosome ectopic exchanges. Further refinements to clinical cytogenetics practices have led to the implementation of molecular cytogenetics techniques such as fluorescent in-situ hybridization (FISH), allowing for rearrangements to be visualized and breakpoints refined using fluorescently labelled painting probes. The next-generation of clinical cytogenetics include the implementation of DNA microarrays, and next-generation sequencing (NGS) technologies such as DNA sequencing to better explore tentative genome architecture changes. The implementation of these cytogenomics techniques allow the genomes of rearrangement carriers to be deciphered at the highest resolution, allowing rearrangements to be detected; breakpoints to be delineated; and, most importantly, potential gene implications of those chromosome rearrangements to be interrogated. Clinical cytogenetics has become an integral tool in the livestock industry, identifying rearrangements and allowing breeders to make informed breeding decisions.

The Meishan pig genome reveals structural variation-mediated gene expression and phenotypic divergence underlying Asian pig domestication

There are wide genomic and phenotypic differences between Asian and European pig breeds, yet the current reference genome is the European Duroc pig genome. A high-quality pig genome is lacking for genetic analysis of agricultural traits in Asian pigs. Here, using a hybrid approach, a high-quality reference genome (MSCAAS v1) for the Asian Meishan breed is assembled with a contig N50 size of 48.05 Mb. MSCAAS v1 outperforms the Duroc genome as a reference genome for Asian breeds. Genomic comparison reveals 49,103 structural variations (SVs) between Meishan and Duroc, 4.02% of which are Asian-specific SVs (AP-SVs). Notably, a 30-Mb hotspot for AP-SVs on chromosome X enriched for genes associated with Asian-pig-specific phenotypes is present in Asian domestic pig breeds, but absent in Asian wild boars, suggesting that Asian domestic breeds share a common ancestor. Interbreed transcriptomics reveals transcriptional suppression roles of AP-SVs in multiple tissues. Finally, transcriptional regulation in the intron of IGF2R is reported, as genomic SV (274-bp deletion) in Tibetan pig limits its growth compared to domestic pig breeds. In summary, this study provides insights regarding the genetic changes underlying pig domestication and presents a benchmark-setting resource for the utilization of agricultural valuable loci in Asian pigs.

A compendium and comparative epigenomics analysis of cis-regulatory elements in the pig genome

Although major advances in genomics have initiated an exciting new era of research, a lack of information regarding cis-regulatory elements has limited the genetic improvement or manipulation of pigs as a meat source and biomedical model. Here, we systematically characterize cis-regulatory elements and their functions in 12 diverse tissues from four pig breeds by adopting similar strategies as the ENCODE and Roadmap Epigenomics projects, which include RNA-seq, ATAC-seq, and ChIP-seq. In total, we generate 199 datasets and identify more than 220,000 cis-regulatory elements in the pig genome. Surprisingly, we find higher conservation of cis-regulatory elements between human and pig genomes than those between human and mouse genomes. Furthermore, the differences of topologically associating domains between the pig and human genomes are associated with morphological evolution of the head and face. Beyond generating a major new benchmark resource for pig epigenetics, our study provides basic comparative epigenetic data relevant to using pigs as models in human biomedical research.

New Gene Markers Expressed in Porcine Oviductal Epithelial Cells Cultured Primary In Vitro Are Involved in Ontological Groups Representing Physiological Processes of Porcine Oocytes

Changes that occur within oviducts after fertilization are dependent on post-ovulation events, including oocyte-oviduct interactions. Although general processes are well-defined, the molecular basis are poorly understood. Recently, new marker genes involved in ‘cell development’, ‘cell growth’, ‘cell differentiation’ and ‘cell maturation’ processes have been identified in porcine oocytes. The aim of the study was to assess the expression profile of genes in primary in vitro cultured oviductal epithelial cells (OECs), clustered in Gene Ontology groups which enveloped markers also identified in porcine oocytes. OECs (from 45 gilts) were surgically removed and cultured in vitro for ≤ 30 days, and then subjected to molecular analyses. The transcriptomic and proteomic profiles of cells cultured during 7, 15 and 30 days were investigated. Additionally, morphological/histochemical analyzes were performed. The results of genes expression profiles were validated after using RT-qPCR. The results showed a significant upregulation of UNC45B, NOX4, VLDLR, ITGB3, FMOD, SGCE, COL1A2, LOX, LIPG, THY1 and downregulation of SERPINB2, CD274, TXNIP, CELA1, DDX60, CRABP2, SLC5A1, IDO1, ANPEP, FST. Detailed knowledge of the molecular pathways occurring in the OECs and the gametes that contact them may contribute both to developments of basic science of physiology, and new possibilities in advanced biotechnology of assisted reproduction.

Analysis of Structural Variants Reveal Novel Selective Regions in the Genome of Meishan Pigs by Whole Genome Sequencing

Structural variants (SVs) represent essential forms of genetic variation, and they are associated with various phenotypic traits in a wide range of important livestock species. However, the distribution of SVs in the pig genome has not been fully characterized, and the function of SVs in the economic traits of pig has rarely been studied, especially for most domestic pig breeds. Meishan pig is one of the most famous Chinese domestic pig breeds, with excellent reproductive performance. Here, to explore the genome characters of Meishan pig, we construct an SV map of porcine using whole-genome sequencing data and report 33,698 SVs in 305 individuals of 55 globally distributed pig breeds. We perform selective signature analysis using these SVs, and a number of candidate variants are successfully identified. Especially for the Meishan pig, 64 novel significant selection regions are detected in its genome. A 140-bp deletion in the Indoleamine 2,3-Dioxygenase 2 (IDO2) gene, is shown to be associated with reproduction traits in Meishan pig. In addition, we detect two duplications only existing in Meishan pig. Moreover, the two duplications are separately located in cytochrome P450 family 2 subfamily J member 2 (CYP2J2) gene and phospholipase A2 group IVA (PLA2G4A) gene, which are related to the reproduction trait. Our study provides new insights into the role of selection in SVs' evolution and how SVs contribute to phenotypic variation in pigs.

Integrated hybrid de novo assembly technologies to obtain high-quality pig genome using short and long reads

With the rapid progress of sequencing technologies, various types of sequencing reads and assembly algorithms have been designed to construct genome assemblies. Although recent studies have attempted to evaluate the appropriate type of sequencing reads and algorithms for assembling high-quality genomes, it is still a challenge to set the correct combination for constructing animal genomes. Here, we present a comparative performance assessment of 14 assembly combinations-9 software programs with different short and long reads of Duroc pig. Based on the results of the optimization process for genome construction, we designed an integrated hybrid de novo assembly pipeline, HSCG, and constructed a draft genome for Duroc pig. Comparison between the new genome and Sus scrofa 11.1 revealed important breakpoints in two S. scrofa 11.1 genes. Our findings may provide new insights into the pan-genome analysis studies of agricultural animals, and the integrated assembly pipeline may serve as a guide for the assembly of other animal genomes.

Cytosine Base Editor (hA3A-BE3-NG)-Mediated Multiple Gene Editing for Pyramid Breeding in Pigs

Pig is an important agricultural economic animal, providing large amount of meat products. With the development of functional genomics and bioinformatics, lots of genes and functional single nucleotide polymorphisms (SNPs) related to disease resistance and (or) economic traits in pigs have been identified, which provides the targets for genetic improvement by genome editing. Base editors (BEs), combining Cas9 nickase and cytidine or adenine deaminase, achieve all four possible transition mutations (C-to-T, A-to-G, T-to-C, and G-to-A) efficiently and accurately without double strand breaks (DSBs) under the protospacer adjacent motif (PAM) sequence of NGG. However, the NGG PAM in canonical CRISPR-Cas9 can only cover approximately 8.27% in the whole genome which limits its broad application. In the current study, hA3A-BE3-NG system was constructed with the fusion of SpCas9-NG variant and hA3A-BE3 to create C-to-T conversion at NGN PAM sites efficiently. The editing efficiency and scope of hA3A-BE3-NG were confirmed in HEK293T cells and porcine fetal fibroblast (PFF) cells. Results showed that the efficiency of hA3A-BE3-NG was much higher than that of hA3A-BE3 on NGH (H = A, C, or T) PAM sites (21.27 vs. 2.81% at average). Further, nonsense and missense mutations were introduced efficiently and precisely via hA3A-BE3-NG in multiple pig economic trait-related genes (CD163, APN, MSTN, and MC4R) in PFF cells by one transfection. The current work indicates the potential applications of hA3A-BE3-NG for pyramid breeding studies in livestock.

Evolutionary and Functional Analysis of Korean Native Pig Using Single Nucleotide Polymorphisms

Time and cost-effective production of next-generation sequencing data has enabled the performance of population-scale comparative and evolutionary studies for various species, which are essential for obtaining the comprehensive insight into molecular mechanisms underlying species- or breed-specific traits. In this study, the evolutionary and functional analysis of Korean native pig (KNP) was performed using single nucleotide polymorphism (SNP) data by comparative and population genomic approaches with six different mammalian species and five pig breeds. We examined the evolutionary history of KNP SNPs, and the specific genes of KNP based on the uniqueness of non-synonymous SNPs among the used species and pig breeds. We discovered the evolutionary trajectory of KNP SNPs within the used mammalian species as well as pig breeds. We also found olfaction-associated functions that have been characterized and diversified during evolution, and quantitative trait loci associated with the unique traits of KNP. Our study provides new insight into the evolution of KNP and serves as a good example for a better understanding of domestic animals in terms of evolution and domestication using the combined approaches of comparative and population genomics.

Pan-genomic open reading frames: A potential supplement of single nucleotide polymorphisms in estimation of heritability and genomic prediction

Pan-genomic open reading frames (ORFs) potentially carry protein-coding gene or coding variant information in a population. In this study, we suggest that pan-genomic ORFs are promising to be utilized in estimation of heritability and genomic prediction. A Saccharomyces cerevisiae dataset with whole-genome SNPs, pan-genomic ORFs, and the copy numbers of those ORFs is used to test the effectiveness of ORF data as a predictor in three prediction models for 35 traits. Our results show that the ORF-based heritability can capture more genetic effects than SNP-based heritability for all traits. Compared to SNP-based genomic prediction (GBLUP), pan-genomic ORF-based genomic prediction (OBLUP) is distinctly more accurate for all traits, and the predictive abilities on average are more than doubled across all traits. For four traits, the copy number of ORF-based prediction(CBLUP) is more accurate than OBLUP. When using different numbers of isolates in training sets in ORF-based prediction, the predictive abilities for all traits increased as more isolates are added in the training sets, suggesting that with very large training sets the prediction accuracy will be in the range of the square root of the heritability. We conclude that pan-genomic ORFs have the potential to be a supplement of single nucleotide polymorphisms in estimation of heritability and genomic prediction.

Pan-Genome of Wild and Cultivated Soybeans

Soybean is one of the most important vegetable oil and protein feed crops. To capture the entire genomic diversity, it is needed to construct a complete high-quality pan-genome from diverse soybean accessions. In this study, we performed individual de novo genome assemblies for 26 representative soybeans that were selected from 2,898 deeply sequenced accessions. Using these assembled genomes together with three previously reported genomes, we constructed a graph-based genome and performed pan-genome analysis, which identified numerous genetic variations that cannot be detected by direct mapping of short sequence reads onto a single reference genome. The structural variations from the 2,898 accessions that were genotyped based on the graph-based genome and the RNA sequencing (RNA-seq) data from the representative 26 accessions helped to link genetic variations to candidate genes that are responsible for important traits. This pan-genome resource will promote evolutionary and functional genomics studies in soybean.

An improved pig reference genome sequence to enable pig genetics and genomics research

Background

The domestic pig (Sus scrofa) is important both as a food source and as a biomedical model given its similarity in size, anatomy, physiology, metabolism, pathology, and pharmacology to humans. The draft reference genome (Sscrofa10.2) of a purebred Duroc female pig established using older clone-based sequencing methods was incomplete, and unresolved redundancies, short-range order and orientation errors, and associated misassembled genes limited its utility.

Results

We present 2 annotated highly contiguous chromosome-level genome assemblies created with more recent long-read technologies and a whole-genome shotgun strategy, 1 for the same Duroc female (Sscrofa11.1) and 1 for an outbred, composite-breed male (USMARCv1.0). Both assemblies are of substantially higher (>90-fold) continuity and accuracy than Sscrofa10.2.

Conclusions

These highly contiguous assemblies plus annotation of a further 11 short-read assemblies provide an unprecedented view of the genetic make-up of this important agricultural and biomedical model species. We propose that the improved Duroc assembly (Sscrofa11.1) become the reference genome for genomic research in pigs.

Diffusion Kurtosis Imaging of Leptin Intervention in Early Hypoxic-Ischemic Brain Edema

The role of leptin in neuroprotection has recently been recognized. However, there are few reports on the use of imaging methods to dynamically evaluate the neuroprotection role of leptin. Diffusion kurtosis imaging (DKI), which is a method used to measure non-Gaussian water diffusion, can reflect the real water diffusion in brain tissues. In this study, a newborn piglet model was used to dynamically evaluate the leptin intervention in early hypoxic-ischemic brain edema via DKI. Thirty-two Yorkshire newborn piglets were divided into three groups: the hypoxic-ischemic encephalopathy (HIE) group, the leptin group, and the control group. DKI scanning was performed at time points of 3, 6, 9, 12, 16, and 24 h after hypoxic-ischemic exposure. After scanning, arterial blood was extracted from all piglets to measure NSE and S100β levels. Then, the brain was completely extracted for pathological examination. In the lesion areas, the MK, Ka, and Kr values in the leptin group were significantly lower than those in the HIE group, the MD, Da, and Dr values showed an opposite trend. The lesion areas in the leptin group were smaller than those of in the HIE group. In addition, the pathological results showed that less cell and organelle injury occurred in the leptin group. Our findings indicate that leptin can effectively reduce hypoxic-ischemic brain edema, and DKI can be more sensitive than conventional diffusivity metrics for visualizing the microstructural changes of HIE. This provides a new clue for the treatment and evaluation of HIE.

Pangenomics Comes of Age: From Bacteria to Plant and Animal Applications

The pangenome refers to a collection of genomic sequence found in the entire species or population rather than in a single individual; the sequence can be core, present in all individuals, or accessory (variable or dispensable), found in a subset of individuals only. While pangenomic studies were first undertaken in bacterial species, developments in genome sequencing and assembly approaches have allowed construction of pangenomes for eukaryotic organisms, fungi, plants, and animals, including two large-scale human pangenome projects. Analysis of the these pangenomes revealed key differences, most likely stemming from divergent evolutionary histories, but also surprising similarities.

Genomic analyses provide insights into breed-of-origin effects from purebreds on three-way crossbred pigs

Crossbreeding is widely used aimed at improving crossbred performance for poultry and livestock. Alleles that are specific to different purebreds will yield a large number of heterozygous single-nucleotide polymorphisms (SNPs) in crossbred individuals, which are supposed to have the power to alter gene function or regulate gene expression. For pork production, a classic three-way crossbreeding system of Duroc × (Landrace × Yorkshire) (DLY) is generally used to produce terminal crossbred pigs with stable and prominent performance. Nonetheless, little is known about the breed-of-origin effects from purebreds on DLY pigs. In this study, we first estimated the distribution of heterozygous SNPs in three kinds of three-way crossbred pigs via whole genome sequencing data originated from three purebreds. The result suggested that DLY is a more effective strategy for three-way crossbreeding as it could yield more stably inherited heterozygous SNPs. We then sequenced a DLY pig family and identified 95, 79, 132 and 42 allele-specific expression (ASE) genes in adipose, heart, liver and skeletal muscle, respectively. Principal component analysis and unrestricted clustering analyses revealed the tissue-specific pattern of ASE genes, indicating the potential roles of ASE genes for development of DLY pigs. In summary, our findings provided a lot of candidate SNP markers and ASE genes for DLY three-way crossbreeding system, which may be valuable for pig breeding and production in the future.

‘Cell cycle’ and ‘cell death’- related genes are differentially expressed during long – term in vitro real-time cultivation of porcine oviductal epithelial cells

Alterations in cells depend on their genetic material, its activation and translation of the products. The genes responsible for the cell cycle processes and apoptosis of porcine oviductal cells have been presented in our study. The processes occurring in the reproductive system of females are extremely complex and require in-depth knowledge. Thanks to in vitro studies on the fallopian tube epithelium cells, we can get closer to understanding the biochemical and morphological changes occurring in mammalian organisms. Our research was conducted on fallopian tubes obtained from commercially bred pigs and its aim was to assess the expression profile of genes responsible for the most important processes of cellular life. Cell cultures were carried out for 30 days, with the obtained cells subjected to molecular analysis. We have shown significant regulation of “cell death” and “cell cycle” genes, some of which are related to the reproductive system. The alterations in transcriptomic profile and mutual relations between the genes were analyzed and related to the literature findings. The knowledge gained could help in identifying new potential markers of the in vitro occurrence of processes described by the ontology groups of interest.

Running title: pig, oocytes, microarray assays, in vitro maturation (IVM)

Development and Genome Sequencing of a Laboratory-Inbred Miniature Pig Facilitates Study of Human Diabetic Disease

Pig has been proved to be a valuable large animal model used for research on diabetic disease. However, their translational value is limited given their distinct anatomy and physiology. For the last 30 years, we have been developing a laboratory Asian miniature pig inbred line (Bama miniature pig [BM]) from the primitive Bama xiang pig via long-term selective inbreeding. Here, we assembled a BM reference genome at full chromosome-scale resolution with a total length of 2.49 Gb. Comparative and evolutionary genomic analyses identified numerous variations between the BM and commercial pig (Duroc), particularly those in the genetic loci associated with the features advantageous to diabetes studies. Resequencing analyses revealed many differentiated gene loci associated with inbreeding and other selective forces. These together with transcriptome analyses of diabetic pig models provide a comprehensive genetic basis for resistance to diabetogenic environment, especially related to energy metabolism.

Building a sequence map of the pig pan-genome from multiple de novo assemblies and Hi-C data

Pigs were domesticated independently in the Near East and China, indicating that a single reference genome from one individual is unable to represent the full spectrum of divergent sequences in pigs worldwide. Therefore, 12 de novo pig assemblies from Eurasia were compared in this study to identify the missing sequences from the reference genome. As a result, 72.5 Mb of non-redundant sequences (∼3% of the genome) were found to be absent from the reference genome (Sscrofa11.1) and were defined as pan-sequences. Of the pan-sequences, 9.0 Mb were dominant in Chinese pigs, in contrast with their low frequency in European pigs. One sequence dominant in Chinese pigs contained the complete genic region of the tazarotene-induced gene 3 (TIG3) gene which is involved in fatty acid metabolism. Using flanking sequences and Hi-C based methods, 27.7% of the sequences could be anchored to the reference genome. The supplementation of these sequences could contribute to the accurate interpretation of the 3D chromatin structure. A web-based pan-genome database was further provided to serve as a primary resource for exploration of genetic diversity and promote pig breeding and biomedical research.

A Western-Type Dietary Pattern Induces an Atherogenic Gene Expression Profile in the Coronary Arteries of the Ossabaw Pig

BACKGROUND

Current cardiovascular risk reduction guidance focuses on shifts in dietary patterns, rather than single foods or nutrients. Experimental studies are needed to identify the mechanisms by which food-based diets affect the development and progression of atherosclerosis.

OBJECTIVES

The aim of this study was to investigate the effect of 2 food-based dietary patterns and statin therapy on the transcriptome of the left anterior descending coronary artery of the Ossabaw pig.

METHODS

Pigs were randomly assigned to 1 of 4 groups and fed isocaloric diets for 6 mo; Heart Healthy-style diet (HHD) (high in unsaturated fat, unrefined grain, fruits/vegetables) or Western-style diet (WD) (high in saturated fat, cholesterol, refined grain), with or without atorvastatin. A 2-factor edge R analysis was used to determine differential gene expression in the left anterior descending coronary artery.

RESULTS

Relative to the HHD, the WD resulted in the differential expression of 143 genes, of which 139 genes were upregulated and 4 genes were downregulated (all log fold change ≥0.6, false discovery rate <0.10). The WD, compared with the HHD, resulted in the statistically significant upregulation of 8 atherosclerosis-associated pathways implicated in immune and inflammatory processes. There were no genes with significant differential expression attributable to statin therapy.

CONCLUSIONS

These data suggest that a WD induces alterations in the transcriptome of the coronary artery consistent with an inflammatory atherogenic phenotype in the Ossabaw pig with no significant modification by concurrent statin therapy.

Deep Genome Resequencing Reveals Artificial and Natural Selection for Visual Deterioration, Plateau Adaptability and High Prolificacy in Chinese Domestic Sheep

Sheep were one of the earliest domesticated animals. Both artificial and natural selection during domestication has resulted in remarkable changes in behavioral, physiological, and morphological phenotypes; however, the genetic mechanisms underpinning these changes remain unclear, particularly for indigenous Chinese sheep. In the present study, we performed pooled whole-genome resequencing of 338 sheep from five breeds representative of indigenous Chinese breeds and compared them to the wild ancestors of domestic sheep (Asian mouflon, Ovis orientalis) for detection of genome-wide selective sweeps. Comparative genomic analysis between domestic sheep and Asian mouflon showed that selected regions were enriched for genes involved in bone morphogenesis, growth regulation, and embryonic and neural development in domestic sheep. Moreover, we identified several vision-associated genes with funtional mutations, such as PDE6B (c.G2994C/p.A982P and c.C2284A/p.L762M mutations), PANK2, and FOXC1/GMSD in all five Chinese native breeds. Breed-specific selected regions were determined including genes such as CYP17 for hypoxia adaptability in Tibetan sheep and DNAJB5 for heat tolerance in Duolang sheep. Our findings provide insights into the genetic mechanisms underlying important phenotypic changes that have occurred during sheep domestication and subsequent selection.

Resequencing core accessions of a pedigree identifies derivation of genomic segments and key agronomic trait loci during cotton improvement

Upland cotton (Gossypium hirsutum) is the world's largest source of natural fibre and dominates the global textile industry. Hybrid cotton varieties exhibit strong heterosis that confers high fibre yields, yet the genome-wide effects of artificial selection that have influenced Upland cotton during its breeding history are poorly understood. Here, we resequenced Upland cotton genomes and constructed a variation map of an intact breeding pedigree comprising seven elite and 19 backbone parents. Compared to wild accessions, the 26 pedigree accessions underwent strong artificial selection during domestication that has resulted in reduced genetic diversity but stronger linkage disequilibrium and higher extents of selective sweeps. In contrast to the backbone parents, the elite parents have acquired significantly improved agronomic traits, with an especially pronounced increase in the lint percentage. Notably, identify by descent (IBD) tracking revealed that the elite parents inherited abundant beneficial trait segments and loci from the backbone parents and our combined analyses led to the identification of a core genomic segment which was inherited in the elite lines from the parents Zhong 7263 and Ejing 1 and that was strongly associated with lint percentage. Additionally, SNP correlation analysis of this core segment showed that a non-synonymous SNP (A-to-G) site in a gene encoding the cell wall-associated receptor-like kinase 3 (GhWAKL3) protein was highly correlated with increased lint percentage. Our results substantially increase the valuable genomics resources available for future genetic and functional genomics studies of cotton and reveal insights that will facilitate yield increases in the molecular breeding of cotton.

Absence of a human ortholog of rodent Kupffer cell galactose-binding receptor encoded by the CLEC4f gene

The murine CLEC4f gene encodes the Kupffer cell receptor, a galactose-binding receptor containing a C-type carbohydrate-recognition domain. Orthologs have been identified in nearly 100 species. The receptors from rat and mouse have previously been characterized and data presented here show that functional CLEC4f protein is expressed in domestic cattle (Bos taurus). However, the human CLEC4f gene does not encode a functional receptor because a mutation in the splice acceptor site of the final exon prevents appropriate splicing and a missense mutation disrupts the sugar-binding site. Transcriptomic and PCR analysis of transcripts confirms the absence of a spliced transcript containing the final exon and only background levels of transcripts are detected in human tissues. These mutations are also present in the CLEC4f gene in Neanderthals. In contrast to humans, closely related species, including chimpanzees, do have CLEC4f genes that encode full-length receptors. Affinity chromatography and glycan array results demonstrate that the chimpanzee, bovine and murine proteins all bind to galactose, but they show preferences for different subsets of galactose-containing glycans. In non-human primates, the receptor is expressed in spleen rather than in liver. The results indicate that the CLEC4f protein probably has distinct functions in different species. Absence of the receptor precludes using it for targeting of glycoconjugates to cells in human liver. The fact that CLEC4f protein is expressed in spleen in non-human primates and the close evolutionary relationship of the CLEC4f protein to langerin (CD207) suggest that it may function in the immune system, possibly as a pathogen receptor.

The genomic organization and expression pattern of the low-affinity Fc gamma receptors (FcγR) in the Göttingen minipig

Safety and efficacy of therapeutic antibodies are often dependent on their interaction with Fc receptors for IgG (FcγRs). The Göttingen minipig represents a valuable species for biomedical research but its use in preclinical studies with therapeutic antibodies is hampered by the lack of knowledge about the porcine FcγRs. Genome analysis and sequencing now enabled the localization of the previously described FcγRIIIa in the orthologous location to human FCGR3A. In addition, we identified nearby the gene coding for the hitherto undescribed putative porcine FcγRIIa. The 1′241 bp long FCGR2A cDNA translates to a 274aa transmembrane protein containing an extracellular region with high similarity to human and cattle FcγRIIa. Like in cattle, the intracellular part does not contain an immunoreceptor tyrosine-based activation motif (ITAM) as in human FcγRIIa. Flow cytometry of the whole blood and single-cell RNA sequencing of peripheral blood mononuclear cells (PBMCs) of Göttingen minipigs revealed the expression profile of all porcine FcγRs which is compared to human and mouse. The new FcγRIIa is mainly expressed on platelets making the minipig a good model to study IgG-mediated platelet activation and aggregation. In contrast to humans, minipig blood monocytes were found to express inhibitory FcγRIIb that could lead to the underestimation of FcγR-mediated effects of monocytes observed in minipig studies with therapeutic antibodies.

Age-associated changes in DNA methylation and expression of the TNFα gene in pigs

DNA methylation is an important mediator of gene expression regulation and has been shown to be closely linked to aging. Immune-related genes tend to be influenced by DNA methylation at different ages. To explore DNA methylation changes in the porcine TNFα gene and analyze their potential effects on gene expression, we measured the methylation level of the TNFα promoter and TNFα mRNA expression in the spleen of Meishan piglets at six developmental stages (1, 7, 14, 21, 28 and 35 days old) by bisulfite sequencing PCR and quantitative PCR. The results revealed a trend for TNFα promoter methylation level to increase and mRNA expression to decrease with age. Correlation analysis showed a significant negative association between methylation level and mRNA expression (Pearson's r = -0.775, P = 4.87E-07). In addition, the transcription factor Sp1 was revealed to bind with the TNFα promoter and regulate TNFα expression. DNA methylation in the TNFα promoter was found to decrease the promoter's activity, and methylation inhibition could enhance the expression level of TNFα, providing functional evidence that promoter methylation controls TNFα expression. Together, our data provide insights into age-associated changes in promoter methylation of the TNFα gene in the spleen and contribute to our understanding of regulatory mechanisms for TNFα expression in the immune system of pigs.

Evidence of evolutionary history and selective sweeps in the genome of Meishan pig reveals its genetic and phenotypic characterization

Background

Meishan is a pig breed indigenous to China and famous for its high fecundity. The traits of Meishan are strongly associated with its distinct evolutionary history and domestication. However, the genomic evidence linking the domestication of Meishan pigs with its unique features is still poorly understood. The goal of this study is to investigate the genomic signatures and evolutionary evidence related to the phenotypic traits of Meishan via large-scale sequencing.

Results

We found that the unique domestication of Meishan pigs occurred in the Taihu Basin area between the Majiabang and Liangzhu Cultures, during which 300 protein-coding genes have underwent positive selection. Notably, enrichment of the FoxO signaling pathway with significant enrichment signal and the harbored gene IGF1R were likely associated with the high fertility of Meishan pigs. Moreover, NFKB1 exhibited strong selective sweep signals and positively participated in hyaluronan biosynthesis as the key gene of NF-kB signaling, which may have resulted in the wrinkled skin and face of Meishan pigs. Particularly, three population-specific synonymous single-nucleotide variants occurred in PYROXD1, MC1R, and FAM83G genes; the T305C substitution in the MCIR gene explained the black coat of the Meishan pigs well. In addition, the shared haplotypes between Meishan and Duroc breeds confirmed the previous Asian-derived introgression and demonstrated the specific contribution of Meishan pigs.

Conclusions

These findings will help us explain the unique genetic and phenotypic characteristics of Meishan pigs and offer a plausible method for their utilization of Meishan pigs as valuable genetic resources in pig breeding and as an animal model for human wrinkled skin disease research.

Comparative transcriptomics of 5 high-altitude vertebrates and their low-altitude relatives

Background

Species living at high altitude are subject to strong selective pressures due to inhospitable environments (e.g., hypoxia, low temperature, high solar radiation, and lack of biological production), making these species valuable models for comparative analyses of local adaptation. Studies that have examined high-altitude adaptation have identified a vast array of rapidly evolving genes that characterize the dramatic phenotypic changes in high-altitude animals. However, how high-altitude environment shapes gene expression programs remains largely unknown.

Findings

We generated a total of 910 Gb of high-quality RNA-seq data for 180 samples derived from 6 tissues of 5 agriculturally important high-altitude vertebrates (Tibetan chicken, Tibetan pig, Tibetan sheep, Tibetan goat, and yak) and their cross-fertile relatives living in geographically neighboring low-altitude regions. Of these, ∼75% reads could be aligned to their respective reference genomes, and on average ∼60% of annotated protein coding genes in each organism showed FPKM expression values greater than 0.5. We observed a general concordance in topological relationships between the nucleotide alignments and gene expression-based trees. Tissue and species accounted for markedly more variance than altitude based on either the expression or the alternative splicing patterns. Cross-species clustering analyses showed a tissue-dominated pattern of gene expression and a species-dominated pattern for alternative splicing. We also identified numerous differentially expressed genes that could potentially be involved in phenotypic divergence shaped by high-altitude adaptation.

Conclusions

These data serve as a valuable resource for examining the convergence and divergence of gene expression changes between species as they adapt or acclimatize to high-altitude environments.

Genome-wide detection of selection signatures in Chinese indigenous Laiwu pigs revealed candidate genes regulating fat deposition in muscle

Background

Currently, genome-wide scans for positive selection signatures in commercial breed have been investigated. However, few studies have focused on selection footprints of indigenous breeds. Laiwu pig is an invaluable Chinese indigenous pig breed with extremely high proportion of intramuscular fat (IMF), and an excellent model to detect footprint as the result of natural and artificial selection for fat deposition in muscle.

Result

In this study, based on GeneSeek Genomic profiler Porcine HD data, three complementary methods, F_ST, iHS (integrated haplotype homozygosity score) and CLR (composite likelihood ratio), were implemented to detect selection signatures in the whole genome of Laiwu pigs. Totally, 175 candidate selected regions were obtained by at least two of the three methods, which covered 43.75 Mb genomic regions and corresponded to 1.79% of the genome sequence. Gene annotation of the selected regions revealed a list of functionally important genes for feed intake and fat deposition, reproduction, and immune response. Especially, in accordance to the phenotypic features of Laiwu pigs, among the candidate genes, we identified several genes, NPY1R, NPY5R, PIK3R1 and JAKMIP1, involved in the actions of two sets of neurons, which are central regulators in maintaining the balance between food intake and energy expenditure.

Conclusions

Our results identified a number of regions showing signatures of selection, as well as a list of functionally candidate genes with potential effect on phenotypic traits, especially fat deposition in muscle. Our findings provide insights into the mechanisms of artificial selection of fat deposition and further facilitate follow-up functional studies.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0622-y) contains supplementary material, which is available to authorized users.

Long non-coding RNAs and mRNAs profiling during spleen development in pig

Genome-wide transcriptomic studies in humans and mice have become extensive and mature. However, a comprehensive and systematic understanding of protein-coding genes and long non-coding RNAs (lncRNAs) expressed during pig spleen development has not been achieved. LncRNAs are known to participate in regulatory networks for an array of biological processes. Here, we constructed 18 RNA libraries from developing fetal pig spleen (55 days before birth), postnatal pig spleens (0, 30, 180 days and 2 years after birth), and the samples from the 2-year-old Wild Boar. A total of 15,040 lncRNA transcripts were identified among these samples. We found that the temporal expression pattern of lncRNAs was more restricted than observed for protein-coding genes. Time-series analysis showed two large modules for protein-coding genes and lncRNAs. The up-regulated module was enriched for genes related to immune and inflammatory function, while the down-regulated module was enriched for cell proliferation processes such as cell division and DNA replication. Co-expression networks indicated the functional relatedness between protein-coding genes and lncRNAs, which were enriched for similar functions over the series of time points examined. We identified numerous differentially expressed protein-coding genes and lncRNAs in all five developmental stages. Notably, ceruloplasmin precursor (CP), a protein-coding gene participating in antioxidant and iron transport processes, was differentially expressed in all stages. This study provides the first catalog of the developing pig spleen, and contributes to a fuller understanding of the molecular mechanisms underpinning mammalian spleen development.

The porcine translational research database: a manually curated, genomics and proteomics-based research resource

BACKGROUND

The use of swine in biomedical research has increased dramatically in the last decade. Diverse genomic- and proteomic databases have been developed to facilitate research using human and rodent models. Current porcine gene databases, however, lack the robust annotation to study pig models that are relevant to human studies and for comparative evaluation with rodent models. Furthermore, they contain a significant number of errors due to their primary reliance on machine-based annotation. To address these deficiencies, a comprehensive literature-based survey was conducted to identify certain selected genes that have demonstrated function in humans, mice or pigs.

RESULTS

The process identified 13,054 candidate human, bovine, mouse or rat genes/proteins used to select potential porcine homologs by searching multiple online sources of porcine gene information. The data in the Porcine Translational Research Database (( http://www.ars.usda.gov/Services/docs.htm?docid=6065 ) is supported by >5800 references, and contains 65 data fields for each entry, including >9700 full length (5' and 3') unambiguous pig sequences, >2400 real time PCR assays and reactivity information on >1700 antibodies. It also contains gene and/or protein expression data for >2200 genes and identifies and corrects 8187 errors (gene duplications artifacts, mis-assemblies, mis-annotations, and incorrect species assignments) for 5337 porcine genes.

CONCLUSIONS

This database is the largest manually curated database for any single veterinary species and is unique among porcine gene databases in regard to linking gene expression to gene function, identifying related gene pathways, and connecting data with other porcine gene databases. This database provides the first comprehensive description of three major Super-families or functionally related groups of proteins (Cluster of Differentiation (CD) Marker genes, Solute Carrier Superfamily, ATP binding Cassette Superfamily), and a comparative description of porcine microRNAs.

Origin and dispersal of early domestic pigs in northern China

It is widely accepted that modern pigs were domesticated independently at least twice, and Chinese native pigs are deemed as direct descendants of the first domesticated pigs in the corresponding domestication centers. By analyzing mitochondrial DNA sequences of an extensive sample set spanning 10,000 years, we find that the earliest pigs from the middle Yellow River region already carried the maternal lineages that are dominant in both younger archaeological populations and modern Chinese pigs. Our data set also supports early Neolithic pig utilization and a long-term in situ origin for northeastern Chinese pigs during 8,000–3,500 BP, suggesting a possibly independent domestication in northeast China. Additionally, we observe a genetic replacement in ancient northeast Chinese pigs since 3,500 BP. The results not only provide increasing evidence for pig origin in the middle Yellow River region but also depict an outline for the process of early pig domestication in northeast China.

Characterization of porcine simple sequence repeat variation on a population scale with genome resequencing data

Simple sequence repeats (SSRs) are used as polymorphic molecular markers in many species. They contribute very important functional variations in a range of complex traits; however, little is known about the variation of most SSRs in pig populations. Here, using genome resequencing data, we identified ~0.63 million polymorphic SSR loci from more than 100 individuals. Through intensive analysis of this dataset, we found that the SSR motif composition, motif length, total length of alleles and distribution of alleles all contribute to SSR variability. Furthermore, we found that CG-containing SSRs displayed significantly lower polymorphism and higher cross-species conservation. With a rigorous filter procedure, we provided a catalogue of 16,527 high-quality polymorphic SSRs, which displayed reliable results for the analysis of phylogenetic relationships and provided valuable summary statistics for 30 individuals equally selected from eight local Chinese pig breeds, six commercial lean pig breeds and Chinese wild boars. In addition, from the high-quality polymorphic SSR catalogue, we identified four loci with potential loss-of-function alleles. Overall, these analyses provide a valuable catalogue of polymorphic SSRs to the existing pig genetic variation database, and we believe this catalogue could be used for future genome-wide genetic analysis.