Search

Inbreeding and loss of genetic diversity are serious problems in local Chinese pig breeds. Runs of homozygosity (ROHs) are contiguous lengths of homozygous genotypes that can provide information on inbreeding levels, mating schemes, selection pressure, and genetic events. The Anqing six-end-white (AQ) pig, an important autochthonous pig breed bred in the Anhui province of China, plays a key role in the local pig industry. In recent decades, AQ pigs have become a closed breed with a small population size in conservation farms, raising the issue of inbreeding decline. In this study, we used 10× resequencing to detect the whole genome of 24 AQ pigs and six Asian wild boars (AWBs). We used the Plink software to analyze the homozygosity of the genome and the distribution of ROHs in the genome based on the resequencing data. We obtained 935.04 Gb of raw data from the resequencing results and more than 45 822 239 SNPs. Additionally, we identified 28 702 ROHs. Compared with AWB, AQ pigs had higher ROH numbers, longer ROH rates, and higher F

(1) Background: Piglet diarrhea is one of the most serious diseases in pigs and has brought great economic losses to the pig industry. Alteration of the gut microbiota is an important factor in the etiology of piglet diarrhea. Therefore, this study aimed to analyze the differences in the gut microbial structures and fecal metabolic profile between post-weaning diarrhea and healthy Chinese Wannan Black pigs. (2) Methods: An integrated approach of 16S rRNA gene sequencing combined with LC/MS-based metabolomics was employed in this study. (3) Results: We found an increase in the relative abundance of the bacterial genus Campylobacter and a decrease in phylum Bacteroidetes and the species Streptococcus gallolyticus subsp. macedonicus. (S. macedonicus) in piglet diarrhea. Meanwhile, obvious changes in the fecal metabolic profile of diarrheic piglets were also detected, particularly higher levels of polyamines (spermine and spermidine). Moreover, there were substantial associations between the disturbed gut microbiota and the altered fecal metabolites, especially a strong positive relationship between spermidine and Campylobacter. (4) Conclusions: These observations may provide novel insights into potential etiologies related to post-weaning diarrhea and further enhance our understanding of the role of gut microbiota in host homeostasis and in modulating gut microbial structure.

Copy number variation (CNV) has been widely used to study the evolution of different species. We first discovered different CNVs in 24 Anqingliubai pigs and 6 Asian wild boars using next-generation sequencing at the whole-genome level with 10× depth to understand the relationship between genetic evolution and production traits in wild boars and domestic pigs. A total of 97,489 CNVs were identified and divided into 10,429 copy number variation regions (CNVRs), occupying 32.06% of the porcine genome. Chromosome 1 had the most CNVRs, and chromosome 18 had the least. Ninety-six CNVRs were selected using VST 1% based on the signatures of all CNVRs, and sixty-five genes were identified in the selected regions. These genes were strongly correlated with traits distinguishing groups by enrichment in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways, such as growth (CD36), reproduction (CIT, RLN), detoxification (CYP3A29), and fatty acid metabolism (ELOVL6). The QTL overlapping regions were associated with meat traits, growth, and immunity, which was consistent with CNV analysis. Our findings increase the understanding of evolved genome structural variations between wild boars and domestic pigs, and provide new molecular biomarkers to guide breeding and the efficient use of available genetic resources.

As a distinguished Chinese indigenous pig breed that exhibits disease resistance and high meat quality, the Anqing six-end-white (AQ) pig represents a valuable germplasm resource for improving the quality of the pig breeding industry. In this study, 24 AQ pigs that were distantly blood-related and 6 Asian Wild Boar (AWB) were selected for 10× deep-genome resequencing. The signatures of the selection were analyzed to explore the genetic basis of their germplasm characteristics and to identify excellent germplasm-related functional genes based on NGS data. A total of 49,289,052 SNPs and 6,186,123 indels were detected across the genome in 30 pigs. Most of the genetic variations were synonym mutations and existed in the intergenic region. We identified 275 selected regions (top 1%) harboring 85 genes by applying a crossover approach based on genetic differentiation (FST) and polymorphism levels (π ratio). Some genes were found to be positively selected in AQ pigs’ breeding. The SMPD4 and DDX18 genes were involved in the immune response to pseudorabies virus (PRV) and porcine reproductive and respiratory syndrome virus (PRRSV). The BCL6 and P2RX6 genes were involved in biological regulation of immune T cells and phagocytes. The SLC7A4 and SPACA4 genes were related to reproductive performance. The MSTN and HIF1A genes were related to fat deposition and muscle development. Moreover, 138 overlapping regions were detected in selected regions and ROH islands of AQ pigs. Additionally, we found that the QTLs with the most overlapping regions were related to back fat thickness, meat color, pH value, fatty acid content, immune cells, parasitic immunity, and bacterial immunity. Based on functional enrichment analysis and QTLs mapping, we conducted further research on the molecular genetic basis of germplasm traits (disease resistance and excellent meat quality). These results are a reliable resource for conserving germplasm resources and exploiting molecular markers of AQ pigs.

Intramuscular fat (IMF) content is vital for pork quality, serving an important role in economic performance in pig industry. Non-coding RNAs, with mRNAs, are involved in IMF deposition; however, their functions and regulatory mechanisms in porcine IMF remain elusive. This study assessed the whole transcriptome expression profiles of the Longissimus dorsi muscle of pigs with high (H) and low (L) IMF content to identify genes implicated in porcine IMF adipogenesis and their regulatory functions. Hundreds of differentially expressed RNAs were found to be involved in fatty acid metabolic processes, lipid metabolism, and fat cell differentiation. Furthermore, combing co-differential expression analyses, we constructed competing endogenous RNAs (ceRNA) regulatory networks, showing crosstalk among 30 lncRNAs and 61 mRNAs through 20 miRNAs, five circRNAs and 11 mRNAs through four miRNAs, and potential IMF deposition-related ceRNA subnetworks. Functional lncRNAs and circRNAs (such as MSTRG.12440.1, ENSSSCT00000066779, novel_circ_011355, novel_circ_011355) were found to act as ceRNAs of important lipid metabolism-related mRNAs (LEP, IP6K1, FFAR4, CEBPA, etc.) by sponging functional miRNAs (such as ssc-miR-196a, ssc-miR-200b, ssc-miR10391, miR486-y). These findings provide potential regulators and molecular regulatory networks that can be utilized for research on IMF traits in pigs, which would aid in marker-assisted selection to improve pork quality.

Background Porcine epidemic diarrhea virus (PEDV) is a causative agent of serious viral enteric disease in suckling pigs. Such diseases cause considerable economic losses in the global swine industry. Enhancing our knowledge of PEDV-induced transcriptomic responses in host cells is imperative to understanding the molecular mechanisms involved in the immune response. Here, we analyzed the transcriptomic profile of intestinal porcine epithelial cell line J2 (IPEC-J2) after infection with a classical strain of PEDV to explore the host response. Results In total, 854 genes were significantly differentially expressed after PEDV infection, including 716 upregulated and 138 downregulated genes. Functional annotation analysis revealed that the differentially expressed genes were mainly enriched in the influenza A, TNF signaling, inflammatory response, cytokine receptor interaction, and other immune-related pathways. Next, the putative promoter regions of the 854 differentially expressed genes were examined for the presence of transcription factor binding sites using the MEME tool. As a result, 504 sequences (59.02%) were identified as possessing at least one binding site of signal transducer and activator of transcription (STAT), and five STAT transcription factors were significantly induced by PEDV infection. Furthermore, we revealed the regulatory network induced by STAT members in the process of PEDV infection. Conclusion Our transcriptomic analysis described the host genetic response to PEDV infection in detail in IPEC-J2 cells, and suggested that STAT transcription factors may serve as key regulators in the response to PEDV infection. These results further our understanding of the pathogenesis of PEDV.

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 Toll signaling pathway was initially identified for its involvement in the control of early embryogenesis. It was later shown to be also part of a major innate immune pathway controlling the expression of anti-microbial peptides in many eukaryotes including humans; cactus, the essential negative regulator of this pathway in flies, was found to be induced in parallel to the Toll-dependent activation process during immune defenses. We were interested in the mechanisms of this dual effect and provide here evidence that upon pathogenic stimuli, dorsal, one of the transcription factors of the fly Toll pathway, can induce the expression of the E3 ligase Bre1. We further show that Bre1 complexes with the E2 Rad6 to mono-ubiquitinate histone H2B and to promote the transcription of cactus to achieve homeostasis of the Toll immune response. Our studies characterize a Toll signal-dependent regulatory machinery in governing the Toll pathway in Drosophila.

Today, many state-of-the-art automatic speech recognition (ASR) systems apply all-neural models that map audio to word sequences trained end-to-end along one global optimisation criterion in a fully data driven fashion. These models allow high precision ASR for domains and words represented in the training material but have difficulties recognising words that are rarely or not at all represented during training, i.e. trending words and new named entities. In this paper, we use a text-to-speech (TTS) engine to provide synthetic audio for out-of-vocabulary (OOV) words. We aim to boost the recognition accuracy of a recurrent neural network transducer (RNN-T) on OOV words by using the extra audio-text pairs, while maintaining the performance on the non-OOV words. Different regularisation techniques are explored and the best performance is achieved by fine-tuning the RNN-T on both original training data and extra synthetic data with elastic weight consolidation (EWC) applied on the encoder. This yields a 57% relative word error rate (WER) reduction on utterances containing OOV words without any degradation on the whole test set., To appear in Proc. ICASSP2021, June 06-11, 2021, Toronto, Ontario, Canada

Alternative splicing (AS) is a process during gene expression that results in a single gene coding for different protein variants. AS contributes to transcriptome and proteome diversity. In order to characterize AS in pigs, genome-wide transcripts and AS events were detected using RNA sequencing of 34 different tissues in Duroc pigs. In total, 138,403 AS events and 29,270 expressed genes were identified. An alternative donor site was the most common AS form and accounted for 44% of the total AS events. The percentage of the other three AS forms (exon skipping, alternative acceptor site, and intron retention) was approximately 19%. The results showed that the most common AS events involving alternative donor sites could produce different transcripts or proteins that affect the biological processes. The expression of genes with tissue-specific AS events showed that gene functions were consistent with tissue functions. AS increased proteome diversity and resulted in novel proteins that gained or lost important functional domains. In summary, these findings extend porcine genome annotation and highlight roles that AS could play in determining tissue identity.

The intestine, a high-turnover tissue, plays a critical role in regulating aging and health in both vertebrates and invertebrates. Maintaining the epithelial barrier function of the intestine by preserving innate immune homeostasis significantly delays aging and prevents mortality. In an effort to explore effective chemicals and materials that can improve intestinal integrity, we performed a nonbiased screen utilizing Drosophila as an animal model. We showed that long-term uptake of aspirin markedly prevented age-onset gut leakage, the over-proliferation of intestinal stem cells, and the dysbiosis of commensal microbiota in fruit flies. Mechanistically, aspirin efficiently downregulated chronic activation of intestinal immune deficiency signaling during aging. Furthermore, our in vivo and in vitro biochemical analyses indicated that aspirin is a negative modulator in control of the K63-linked ubiquitination of Imd. Our findings uncover a novel regulatory mechanism by which aspirin positively modulates intestinal homeostasis, thus delaying aging, in Drosophila.

Long noncoding RNA (lncRNA) play a critical role in mammary development and breast cancer biology. Despite their important role in the mammary gland, little is known of the roles of lncRNA in bovine lactation, particularly regarding the molecular processes underlying it. To characterize the role of lncRNA in bovine lactation, 4 samples of Holstein cow mammary gland tissue at peak and late lactation stages were examined after biopsy. We then profiled the transcriptome of the mammary gland using RNA sequencing technology. Further, functional lncRNA-mRNA coexpression pairs were constructed to infer the function of lncRNA using a generalized linear model, followed by gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. More than 1,000 putative lncRNA were identified, 117 of which were differentially expressed between peak and late lactation stages. Bovine lncRNA were shorter, with fewer exon numbers, and expressed at significantly lower levels than protein-coding genes. Seventy-two differentially expressed (DE) lncRNA were coexpressed with 340 different protein-coding genes. The KEGG pathway analysis showed that target mRNA for DE lncRNA were mainly related to lipid and glucose metabolism, including the peroxisome proliferator-activated receptors and 5' adenosine monophosphate-activated protein kinase signaling pathways. Further bioinformatics and integrative analyses revealed that 12 DE lncRNA potentially played important roles in bovine lactation. Our findings provide a valuable resource for future bovine transcriptome studies, facilitate the understanding of bovine lactation biology, and offer functional information for cattle lactation.

Background Alternative splicing (AS) is a process that mRNA precursor splices intron to form the mature mRNA. AS plays important roles in contributing to transcriptome and proteome divert. However, to date there is no research about pig AS in genome-wide level by RNA sequencing. Results To characterize the AS in pigs, herein we detected genome-wide transcripts and events by RNA sequencing technology (RNA-seq) 34 different tissues in Duroc pigs. In total, we identified 138, 403 AS events and 29, 270 expressed genes. We found alternative donor site was the most common AS form, which is accounted for 44% of the total AS events. The percentage of the other 3 AS forms (Exon skipping, Alternative acceptor site and Intron retention) are all around 19%. The results showed that the most common AS events (alternative donor site) can produce different transcripts or different proteins which affect the biological process. Among these AS events, 109, 483 were novel AS events, and the number of alternative donor splice site has increased the most (Accounting for 44% of the novel AS events). Conclusions The expression of gene with tissue specific AS events showed that the functions of these genes were consistent with the tissue function. AS increased proteome diversity and resulted in novel proteins that gained and lost important functional domains. In summary, these findings extend genome annotation and highlight roles that AS acts in tissue identity in pig.

Ribosomal protein L8 (RPL8) was considered as a promising candidate gene for the milk fat percentage trait in dairy cattle in our previous genome-wide association studies, but the mechanism remains to be determined. Here we investigated the molecular mechanism underlying the effect of bovine RPL8 on milk fat percentage. We demonstrated that RPL8 silencing in bovine mammary epithelial cells affected the expression of genes encoding fat-related enzymes (ACACA, FASN, ACSS1, FABP3, SREBP-1, DGAT1, GPAM, PLIN2, PLIN5 and CIDEA). Furthermore, we showed here that a single nucleotide polymorphism, g.-931G > T (chr14:1508300, UMD3.1) in the putative RPL8 promoter region significantly reduced its promoter activity. Interestingly, this decrease in activity was paralleled by lower RPL8 expression in mammary gland tissues of dairy cattle with the homozygous TT genotype compared to that of cattle with the wild-type homozygous GG genotype. Importantly, we found g.-931G > T added a paired box 6 (Pax6)-binding site and this mutation located in the presumed Pax6-binding site. EMSA and co-immunoprecipitation (Co-IP) assays confirmed the interaction between RPL8 and Pax6 and the T allele exhibited a higher affinity of DNA/protein interactions than G allele, suggesting that Pax6 is an important transcription factor for RPL8 expression. In addition, lactating cows with the GG and GT genotypes presented a significant decrease in milk fat percentage compared to cows with TT genotypes. Altogether, our study indicated that g.-931G > T at RPL8 promoter altered its expression by affecting the interplay between Pax6 and RPL8, which may account for the association with milk fat traits. Findings herein first elucidated the biological function of RPL8 gene in milk fat and the identified SNP g.-931G > T may be considered as genetic makers for breeding in dairy cattle.

The intestine, a high-turnover tissue, plays a critical role in regulating aging and health in both vertebrates and invertebrates. Maintaining the epithelial barrier function of the intestine by preserving innate immune homeostasis significantly delays aging and prevents mortality. In an effort to explore effective chemicals and materials that can improve intestinal integrity, we performed a nonbiased screen utilizing Drosophila as an animal model. We showed that long-term uptake of aspirin markedly prevented age-onset gut leakage, the over-proliferation of intestinal stem cells, and the dysbiosis of commensal microbiota in fruit flies. Mechanistically, aspirin efficiently downregulated chronic activation of intestinal immune deficiency signaling during aging. Furthermore, our in vivo and in vitro biochemical analyses indicated that aspirin is a negative modulator in control of the K63-linked ubiquitination of Imd. Our findings uncover a novel regulatory mechanism by which aspirin positively modulates intestinal homeostasis, thus delaying aging, in Drosophila.

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.

Figure S1. The predicted secondary structures of the mitochondrial tRNAs. (A) Comparison of tRNA-Phe secondary structures between the reference structure of NC_000845.1 (left) and the structure with mutation sites (right). The green circle referred to the reference base, while the yellow was the mutation base in our study. (B-H) Comparison of secondary structures of tRNA-Val, tRNA-Ala, tRNA-Cys, tRNA-Gly, tRNA-Leu, tRNA-Glu and tRNA-Thr, respectively. Figure S2. The predicted tertiary structures of the mitochondrial tRNAs. (A) Comparison of tRNA-Phe tertiary structures between the reference structure of NC_000845.1 (left) and the structure with mutation sites (right). (B-H) Comparison of tertiary structures of tRNA-Val, tRNA-Ala, tRNA-Cys, tRNA-Gly, tRNA-Leu, tRNA-Glu and tRNA-Thr, respectively. (DOCX 1582 kb)

Aim: To construct a comprehensive pig noncoding transcriptome and further enhance porcine noncoding genome annotation. Materials & methods: We performed a tissue-based long noncoding RNA (lncRNA) profiling via exploiting 32,212 nonredundant lncRNA isoforms corresponding to 18,676 lncRNA loci across 34 normal pig tissues using high-throughput sequencing. Furthermore, the potential relationship between our identified lncRNAs and known protein-coding genes were globally assessed via a comprehensive computation-based strategy, developing a genome-wide lncRNA-targeted genome draft for further functional studies on noncoding genes. Results & conclusion: Among these lncRNAs, ubiquitously expressed lncRNA appeared at a higher level compared with tissue-specific one. Findings herein provide insight into comprehensive knowledge of porcine noncoding RNAs and further enhance pig noncoding annotation. For ease of accessing the information of the identified lncRNAs, we deposited those with high confidence in the publicly available NONCODE database, providing a valuable resource for facilitating pig noncoding genomic studies.

Mammary gland is an important organ for milk synthesis and secretion. It undergoes dramatic physiological changes to adapt the shift from peak to late lactation stage. Protein plays a final very vital role in many life functions, and the protein changes during different lactation stages potentially reflect the biology of lactation and the functions of mammary gland in cows. In current study, we adopted tandem mass tags label-based quantitative analysis technique and to investigate proteome changes occurring in bovine mammary gland from peak to late lactation stages. A total of 3753 proteins from mammary tissues taken at two lactation points from four individual cows by biopsy were quantified, out of which 179 proteins were expressed differentially between two stages. We observed five new DEPs (AACS, DHCR7, GSTM3, SFRP1 and SFRP4) and nine functional well-studies known proteins (PLIN2, LPIN1, PLIN3, GSN, CD74, MMP2, SOD1, SOD3 and GPX3) related to milk performance and mammary morphology. Bioinformatics analyses of the DEPs showed a majority of the up-regulated proteins during late lactation stage were related to apoptosis and immune process, while the downregulated proteins were mainly involved in localization, lipid metabolic and transport process. This suggests that the mammary gland can adapt to different molecular functions according to the biological need of the animal. From the integrated analysis of the differentially expressed proteins with known quantitative trait loci and genome-wide association study data, we identified 95 proteins may potentially affect milking performance. We expect findings in this study could be a valuable resource for future studies investigating the bovine proteome and functional studies.

DNA methylation changes play essential roles in regulating the activities of genes involved in immune responses. Understanding of variable DNA methylation linked to immune responses may contribute to identifying biologically promising epigenetic markers for pathogenesis of diseases. Here, we generated genome-wide DNA methylation and transcriptomic profiles of six pairs of polyinosinic-polycytidylic acid-treated pig peripheral blood mononuclear cell (PBMC) samples and corresponding controls using methylated DNA immunoprecipitation sequencing and RNA sequencing. Comparative methylome analyses identified 5,827 differentially methylated regions and 615 genes showing differential expression between the two groups. Integrative analyses revealed inverse associations between DNA methylation around transcriptional start site and gene expression levels. Furthermore, 70 differentially methylated and expressed genes were identified such as TNFRSF9, IDO1 and EBI3. Functional annotation revealed the enriched categories including positive regulation of immune system process and regulation of leukocyte activation. These findings demonstrated DNA methylation changes occurring in immune responses of PBMCs to poly I:C stimulation and a subset of genes potentially regulated by DNA methylation in the immune responses. The PBMC DNA methylome provides an epigenetic overview of this physiological system in response to viral infection, and we expect it to constitute a valuable resource for future epigenetic epidemiology studies in pigs.