Search

The mammalian Y chromosome offers a unique perspective on the male reproduction and paternal evolutionary histories. However, further understanding of the Y chromosome biology for most mammals is hindered by the lack of a Y chromosome assembly. This study presents an integrated in silico strategy for identifying and assembling the goat Y‐linked scaffolds using existing data. A total of 11.5 Mb Y‐linked sequences were clustered into 33 scaffolds, and 187 protein‐coding genes were annotated. We also identified high abundance of repetitive elements. A 5.84 Mb subset was further ordered into an assembly with the evidence from the goat radiation hybrid map (RH map). The existing whole‐genome resequencing data of 96 goats (worldwide distribution) were utilized to exploit the paternal relationships among bezoars and domestic goats. Goat paternal lineages were clearly divided into two clades (Y1 and Y2), predating the goat domestication. Demographic history analyses indicated that maternal lineages experienced a bottleneck effect around 2,000 YBP (years before present), after which goats belonging to the A haplogroup spread worldwide from the Near East. As opposed to this, paternal lineages experienced a population decline around the 10,000 YBP. The evidence from the Y chromosome suggests that male goats were not affected by the A haplogroup worldwide transmission, which implies sexually unbalanced contribution to the goat trade and population expansion in post‐Neolithic period., In this paper, we first assembled goat Y chromosome into 33 scaffolds totaling 11.5 Mb, and a subset of 5.84 Mb was further ordered into an assembly with the validation of goat radiation hybrid map. Then, we identified SNPs of 96 goat samples worldwide to exploit the paternal relationships among wild bezoars and domestic goats. The paternal lineages showed a clear division of bezoars and goats into two clades, and they were independent of goat domestication. Further demographical history of paternal and maternal lineages showed evident differences, with a bottleneck effect in mtDNA around 2,000 YBP, which may be related to the worldwide dispersal of goat in the post‐Neolithic period, implied the possibility of gender preference in goat domestication. Our research provided a new perspective of goat evolution history.

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., To date, little is known about the regulatory landscape of the pig genome. Here, the authors characterize cis-regulatory elements in the pig genome using RNA-seq, ChIP-Seq and ATAC-seq, finding a higher degree of usage conservation between pig and human than mouse and human.

Background The feed conversion ratio (FCR) is an important productive trait that greatly affects profits in the pig industry. Elucidating the genetic mechanisms underpinning FCR may promote more efficient improvement of FCR through artificial selection. In this study, we integrated a genome-wide association study (GWAS) with transcriptome analyses of different tissues in Yorkshire pigs (YY) with the aim of identifying key genes and signalling pathways associated with FCR. Results A total of 61 significant single nucleotide polymorphisms (SNPs) were detected by GWAS in YY. All of these SNPs were located on porcine chromosome (SSC) 5, and the covered region was considered a quantitative trait locus (QTL) region for FCR. Some genes distributed around these significant SNPs were considered as candidates for regulating FCR, including TPH2, FAR2, IRAK3, YARS2, GRIP1, FRS2, CNOT2 and TRHDE. According to transcriptome analyses in the hypothalamus, TPH2 exhibits the potential to regulate intestinal motility through serotonergic synapse and oxytocin signalling pathways. In addition, GRIP1 may be involved in glutamatergic and GABAergic signalling pathways, which regulate FCR by affecting appetite in pigs. Moreover, GRIP1, FRS2, CNOT2, and TRHDE may regulate metabolism in various tissues through a thyroid hormone signalling pathway. Conclusions Based on the results from GWAS and transcriptome analyses, the TPH2, GRIP1, FRS2, TRHDE, and CNOT2 genes were considered candidate genes for regulating FCR in Yorkshire pigs. These findings improve the understanding of the genetic mechanisms of FCR and may help optimize the design of breeding schemes.

Transgenic pigs play an important role in biomedicine and agriculture. The "safe harbor" locus maintains consistent foreign gene expression in cells and is important for transgenic pig generation. However, as only several safe harbor loci(Rosa26, pH11 and Pifs501) have been identified in pigs, meeting the needs of the insertion of various foreign genes is difficult. In this study, we develop a novel strategy for the efficient knock-in of gene-of-interest fragments into endogenous beta-actin(ACTB) gene via CRISPR/Cas9 mediated homologous recombination with normal expression of ACTB. Thus, we provide an alternative strategy to integrate exogenous genes into the pig genome that can be applied to agricultural breeding and biomedical models.

The relationship between the conceptus and the maternal uterine environment is crucial for the successful establishment and maintenance of pregnancy in cattle. Gene expression analysis of the conceptus and maternal reproductive tissues is a favourable method to assess the embryonic maternal interaction. The reliability of the commonly used method reverse transcription-quantitative polymerase chain reaction (RT-qPCR) depends on proper normalization to stable reference genes (RGs). The objective of this study was to determine the expression stability of 10 potential RGs in maternal reproductive tissues and foetal tissues, and to analyse the effect of RG selection on the calculation of the relative expression of target genes. The expression stability of 10 potential RGs was analysed in eight different tissues from three pregnant dairy cows. Three programs-GeNorm, NormFinder and Bestkeeper-were used to identify the best RGs. According to all three programs, the most stable RG was CNOT11, whereas the least stable RGs were GAPDH and HPRT1. GeNorm analysis showed that a combination of five RGs (SDHA, PPIA, CNOT11, RPS9 and RPL19) was necessary for appropriate data normalization. However, NormFinder analysis indicated that the combination of CNOT11 and PPIA was the most suitable. When target genes were normalized to these RGs, the relative expression of the Radical S-adenosyl methionine domain containing 2 gene was not affected by the choice of RGs, whereas a large difference was observed in the expression profile of the Nuclear erythroid2-related factor 2 gene between the most stable and least stable RGs. The results indicate that careful selection of RGs is crucial under different conditions, especially for target genes with relatively small fold changes. Furthermore, the results provide useful information for the selection of RGs for evaluating genes affecting bovine reproduction.

Advances in high-throughput sequencing technologies have reduced the cost of genotyping dramatically and led to genomic prediction being widely used in animal and plant breeding, and increasingly in human genetics. Inspired by the efficient computing of linear mixed model and the accurate prediction of Bayesian methods, we propose a machine learning-based method incorporating cross-validation, multiple regression, grid search, and bisection algorithms named KAML that aims to combine the advantages of prediction accuracy with computing efficiency. KAML exhibits higher prediction accuracy than existing methods, and it is available at https://github.com/YinLiLin/KAML.

Semen traits are crucial in commercial pig production since semen from boars is widely used in artificial insemination for both purebred and crossbred pig production. Revealing the genetic architecture of semen traits potentially promotes the efficiencies of improving semen traits through artificial selection. This study is aimed to identify candidate genes related to the semen traits in Duroc boars. First, we identified the genes that were significantly associated with three semen traits, including sperm motility (MO), sperm concentration (CON), and semen volume (VOL) in a Duroc boar population through a genome-wide association study (GWAS). Second, we performed a weighted gene co-expression network analysis (WGCNA). A total of 2, 3, and 20 single-nucleotide polymorphisms were found to be significantly associated with MO, CON, and VOL, respectively. Based on the haplotype block analysis, we identified one genetic region associated with MO, which explained 6.15% of the genetic trait variance. ENSSSCG00000018823 located within this region was considered as the candidate gene for regulating MO. Another genetic region explaining 1.95% of CON genetic variance was identified, and, in this region, B9D2, PAFAH1B3, TMEM145, and CIC were detected as the CON-related candidate genes. Two genetic regions that accounted for 2.23% and 2.48% of VOL genetic variance were identified, and, in these two regions, WWC2, CDKN2AIP, ING2, TRAPPC11, STOX2, and PELO were identified as VOL-related candidate genes. WGCNA analysis showed that, among these candidate genes, B9D2, TMEM145, WWC2, CDKN2AIP, TRAPPC11, and PELO were located within the most significant module eigengenes, confirming these candidate genes’ role in regulating semen traits in Duroc boars. The identification of these candidate genes can help to better understand the genetic architecture of semen traits in boars. Our findings can be applied for semen traits improvement in Duroc boars.

Cartilage dysplasia is one of the important reasons for the weakness of pig limbs and hooves. Porcine rickets with weak limbs and hooves bring huge economic losses to the pig industry. However, research on the development of pig cartilage is lacking. This study investigated the key genes and molecular mechanisms involved in cartilage development via an RNA-seq technique. Samples of proximal tibia cartilage were collected from three normal piglets with 1 day, 14 days, and 28 days of age, respectively, and then these samples were divided into two comparison groups (1-day vs. 14-day group, 14-day vs. 28-day group). Through the transcriptome analysis, 108 differentially expressed genes (DEGs), such as FORL2, were obtained from 1-day vs. 14-day comparison group, and 3602 DEGs were obtained from 14-day vs. 28-day comparison group, including SOX9, BMP6, and MMP13. The gene ontology (GO) functional and KEGG pathway enrichment revealed that many functions of DEGs were related to bone development. The pathways of DEGs from Day 1 vs. Day 14 were mainly enriched in mineral absorption, but the DEGs of Day 14 vs. Day 28 were enriched in osteoclast differentiation. Then, the expression patterns of six candidate genes were verified via qPCR. In conclusion, candidate genes affecting cartilage development in Yorkshire pigs were obtained by transcriptome analysis, and the clues showed that Day 14 to Day 28 is a more active and extensive period in cartilage developments, which played a key role in revealing the molecular mechanism of pig cartilage development basis, also compensating for vacancies in cartilage research.

Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 technology is effective for genome editing and now widely used in life science research. However, the key factors determining its editing efficiency and off-target cleavage activity for single-guide RNA (sgRNA) are poorly documented. Here, we systematically evaluated the effects of sgRNA length on genome editing efficiency and specificity. Results showed that sgRNA 5'-end lengths can alter genome editing activity. Although the number of predicted off-target sites significantly increased after sgRNA length truncation, sgRNAs with different lengths were highly specific. Because only a few predicted off-targets had detectable cleavage activity as determined by Target capture sequencing (TargetSeq). Interestingly, > 20% of the predicted off-targets contained microsatellites for selected sgRNAs targeting the dystrophin gene, which can produce genomic instability and interfere with accurate assessment of off-target cleavage activity. We found that sgRNA activity and specificity can be sensitively detected by TargetSeq in combination with in silico prediction. Checking whether the on- and off-targets contain microsatellites is necessary to improve the accuracy of analyzing the efficiency of genome editing. Our research provides new features and novel strategies for the accurate assessment of CRISPR sgRNA activity and specificity.

Our previous association study revealed the mutation in candidate immune genes (NOD1 & NLRP9) was significantly associated with FEC of Haemonchus contortus infection in Yichang white goats, but the relative expression of mRNA of those genes associated with resistance to H. contortus was not investigated. Aim of the current experiment was to evaluate the susceptible and resistant individuals to nematode infection within the population of Yichang white goat (YWG) and assess the differential level of mRNA expression of those candidate genes in the abomasal tissues of susceptible and resistant goats. Fecal egg count (FEC) was determined using a modified McMaster technique, and the hematological parameter was measured by Mindray auto hematology analyzer. Phenotype data were collected and analyzed using a generalized linear model with SAS statistical program. Field investigate revealed that the prevalence (76%) with maxium parasite load (734.34±84.21epg) of H. contortus occurred in August within the experimental flock FEC in resistant group (103.38±1.20epg) and susceptible group (1180.25±43.53epg) group demonstrated the presence of two distinct goat populations within this breed. Four resistant and four susceptible goats were selected from each group. The parasite infection was established by artificially challenge with 5000 infective L3 larvae of H. contortus. Abomasal tissues were collected from all experimental goats after 42 days of post-infection. FEC, Body weight, packed cell volume, and hemoglobin value were significantly different (P < 0.01) between resistant and susceptible group of goats. Quantitative real-time PCR in abomasal tissues revealed that the expression level of mRNA for NOD1 (P < 0.00001), IFNG (P < 0.0001), NLRP9, TLR8, IL32, (P < 0.001) and IGF1 (P < 0.01) was higher in resistant goat compared to susceptible, except SFTPA1. These findings revealed the presence of genetic resistant individuals to H. contortus within the goat breed and expression of NOD1 and NLRP9 genes proved the positive finding of our previous study. Presence of genetic resistant individuals in Yichang white goat YWG breed could be a good candidate for selective breeding and highly expressed genes related to resistant could be used as biomarkers to develop H. contortus resistant goat population. Progressive Agriculture 30 (2): 194-208, 2019

N6-methyladenosine (m6A) is considered as a reversible RNA modification occurring more frequently on the GAC than AAC context in vivo, which regulates post-transcriptional gene expression in mammal...

Homeobox A10 (Hoxa10) gene is one of the most important candidate genes associated with the reproductive performance of humans and mice. Overexpression of Hoxa10 in mouse endometrium can increase litter size. Moreover, Hoxa10 plays a key role in regulating the embryo implantation of sows. This study aimed to generate transgenic pigs using Hoxa10 via somatic cell nuclear transfer (SCNT). We established seven Hoxa10-transgenic cell lines, and two of the cell lines were selected as nuclear donors for the transfer. A total of 1270 cloned embryos were generated and transferred to five surrogate mothers (Landrace×Yorkshire). Eight cloned male piglets were produced including one with cryptorchidism. Six transgenic piglets grew up healthy and produced 56 offspring. Finally, we obtained six transgenic male pigs and 26 transgenic positive offspring that can be used to further study the regulatory mechanism of Hoxa10 on the reproductive performance of pigs.

Early diagnosis of pregnancy is important in livestock production, but there is no reliable technology used for pregnancy diagnosis within the first three weeks after insemination. During early pregnancy, the expression of interferon-stimulating genes (ISGs) in peripheral blood leukocytes (PBL) is significantly increased. However, due to different strains, detection sample types, detection methods, threshold value, etc. the specific effectiveness of early pregnancy diagnosis using ISGs is worth further study. The purpose of this study was to test interferon-stimulated protein 15 (ISG15), 2'–5'-oligoadenylate synthetase 1 (OAS1) and radical S-adenosyl methionine domain containing 2 (RSAD2) for early pregnancy diagnosis in dairy cows. The expression of ISG15, OAS1, and RSAD2 in PBL of pregnant and non-pregnant heifers on days 0, 14, 18, 21 and 28 after artificial insemination (AI) was detected by fluorescence quantitative polymerase chain reaction (PCR). The sensitivity and specificity of the pregnancy diagnosis was analyzed using expression of these three genes separately or in combination by receiver operating characteristic curve. The combination with the highest accuracy used probe primers and duplex fluorescence quantitative PCR. The single quantitative PCR results showed that expression of ISG15, OAS1 and RSAD2 on day 18 after AI was significantly higher in pregnant than in non-pregnant cows. When these three genes were used separately, or in combination, for early pregnancy diagnosis, the sensitivity for the RSAD2 gene was 100%, and the combination of ISG15 with RSAD2 was 94.7%. The duplex quantitative PCR showed that, although the sensitivity of ISG15 alone was 100%, its specificity was only 88.2% (cut-off value 1.402). The sensitivity of RSAD2 alone was 89.5%, and the specificity was 88.2%; however, when the two genes were used in combination, the sensitivity, specificity and diagnostic cut-off value were consistent with the results of single quantitative PCR. These results indicated that a duplex quantitative PCR assay system for early pregnancy diagnosis in cows using ISG15 and RSAD2 was established. Simultaneous detection of expression of ISG15 and RSAD2 by duplex quantitative PCR can effectively improve the diagnostic accuracy for dairy cows.

Mycobacterium tuberculosis (M.tb) secretes numerous proteins to interfere with host immune response for its long-term survival. As one of the top abundant M.tb secreted proteins, Rv3722c was found to be essential for bacilli growth. However, it remains elusive how this protein interferes with the host immune response and regulates M.tb survival. Here, we confirmed that Rv3722c interacted with host TRAF3 to promote M.tb replication in macrophages. Knock-down of TRAF3 attenuated the effect of Rv3722c on the intracellular M.tb survival. The interaction between Rv3722c and TRAF3 hampered MAPK and NF-κB pathways, resulting in a significant increase of IFN-β expression and decrease of IL-1β, IL-6, IL-12p40, and TNF-α expression. Our study revealed that Rv3722c interacted with TRAF3 and interrupted its downstream pathways to promote M.tb survival in macrophages. These findings facilitate further understanding of the mechanism of M.tb secreted proteins in regulating the host cell immune response and promoting its intracellular survival.

Background Piglet mortality is an economically important complex trait that impacts sow prolificacy in the pig industry. The genetic parameters estimations and genome-wide association studies will help us to better understand the genetic fundamentals of piglet mortality. However, compared with other economically important traits, a little breakthrough in the genetic analyses of the trait has been achieved. Results In this study, we used multi-breed data sets from Yorkshire, Landrace, and Duroc sows and characterized the genetic and genomic properties of mortality rate at birth by treating each parity as a different trait. The heritability of mortality rate from parity I to III were estimated to be 0.0630, 0.1031, and 0.1140, respectively. The phenotypic and genetic correlations with its component traits were all positive with ranges from 0.0897 to 0.9054, and 0.2388 to 0.9999, respectively. Integrating the results, we identified 21 loci that were detected at least by two tools from standard MLM, FarmCPU, BLINK and mrMLM, and these loci were annotated to 22 genes. The annotations revealed that the gene expressions were associated with the reproductive system, nervous system, digestive system, and embryonic development, which are reasonably related to the piglet mortality. Conclusions In brief, the genetic properties of piglet mortality at birth were reported. These findings are expected to provide much information for understanding the genetic and genomic fundamentals of farrowing mortality and also identify candidate molecular markers for breeding practice.

Attachment of conceptus to the endometrial luminal epithelium (LE) is a critical event for early placentation in Eutheria. Since the attachment occurs at a particular site within the uterus, a coordinated communication between three spatially distinct compartments (conceptus and endometrial LE from two anatomical regions of the uterus to which conceptus attaches and does not attach) is essential but remains to be fully characterized. Using the laser capture microdissection (LCM) technique, we firstly developed an approach that can allow us to pair the pig conceptus sample with its nearby endometrial epithelium sample without losing the native spatial information. Then, a comprehensive spatio-temporal transcriptomic profile without losing the original conceptus-endometrium coordinates was constructed. The analysis shows that an apparent difference in transcriptional responses to the conceptus exists between the endometrial LE from the two anatomically distinct regions in the uterus. In addition, we identified the communication pathways that link the conceptus and endometrial LE and found that these pathways have important roles in conceptus attachment. Furthermore, a number of genes whose expression is spatially restricted in the two different anatomical regions within the uterus were characterized for the first time and two of them (SULT2A1 and MEP1B) may cooperatively contribute to establish conceptus attachment in pigs. The results from our study have implications in understanding of conceptus/embryo attachment in pigs and other large polytocous species.

Background Genetic analysis of gene expression level is a promising approach for characterizing candidate genes that are involved in complex economic traits such as meat quality. In the present study, we conducted expression quantitative trait loci (eQTL) and allele-specific expression (ASE) analyses based on RNA-sequencing (RNAseq) data from the longissimus muscle of 189 Duroc × Luchuan crossed pigs in order to identify some candidate genes for meat quality traits. Results Using a genome-wide association study based on a mixed linear model, we identified 7192 cis-eQTL corresponding to 2098 cis-genes (p ≤ 1.33e-3, FDR ≤ 0.05) and 6400 trans-eQTL corresponding to 863 trans-genes (p ≤ 1.13e-6, FDR ≤ 0.05). ASE analysis using RNAseq SNPs identified 9815 significant ASE-SNPs in 2253 unique genes. Integrative analysis between the cis-eQTL and ASE target genes identified 540 common genes, including 33 genes with expression levels that were correlated with at least one meat quality trait. Among these 540 common genes, 63 have been reported previously as candidate genes for meat quality traits, such as PHKG1 (q-value = 1.67e-6 for the leading SNP in the cis-eQTL analysis), NUDT7 (q-value = 5.67e-13), FADS2 (q-value = 8.44e-5), and DGAT2 (q-value = 1.24e-3). Conclusions The present study confirmed several previously published candidate genes and identified some novel candidate genes for meat quality traits via eQTL and ASE analyses, which will be useful to prioritize candidate genes in further studies.

Background The CD4 protein is an important surface marker of T lymphocytes, which can mediate the antigen presentation process by interacting with MHC II and TCR molecules in human and mouse. Results In this study, two haplotypes (A and B) of the CD4 gene were found within Chinese indigenous and Western commercial pig breeds. These two haplotypes were defined by 22 fully linked SNPs in the CDS region of the CD4 gene. The expression level and localization of the CD4 protein were significantly different between haplotypes A and B. Transcriptome analysis revealed that the immune response-related genes and signaling pathways were down-regulated in genotype AA. Finally, three linked functional SNPs were identified, which affected the expression level and membrane localization of the CD4 protein in pigs. These three SNPs led to the replacements of two amino acids in the IgV1 domain of the CD4 protein, and related to the function of the CD4 protein in the immune response. Conclusion These three linked SNPs were the key functional mutation sites in the CD4 gene, which played important roles in the immune response, and could be utilized as new molecular markers in breeding for disease resistance in pigs.

Background: Feed conversion ratio (FCR) is an important productive trait that largely affects profits in pig industry. Elucidating the genetic mechanisms underpinning the FCR potentially promote the efficiencies of improving FCR through artificial selection. In this study, we integrated a genome-wide association study (GWAS) with transcriptome analyses in different tissues in Yorkshire pigs (YY), aimed at identifying key genes and signaling pathways significantly associated with FCR.Results: A total of 61 significant single nucleotide polymorphism (SNPs) were detected by GWAS in YY. All of these SNPs are located on porcine chromosome (SSC) 5 and the covered region was considered as a quantitative trait locus (QTL) region for FCR. Some genes that distributed around these significant SNPs were considered as the candidates for regulating FCR, including TPH2, FAR2, IRAK3, YARS2, GRIP1, FRS2, CNOT2 and TRHDE. According to the transcriptome analyses in hypothalamus, TPH2 exhibits abilities of regulating the intestinal motility by a serotonergic synapse and an oxytocin signaling pathway. In addition, GRIP1 is involved in a glutamatergic and GABAergic signaling pathway, which regulates FCR through affecting the appetite in pigs. Moreover, GRIP1, FRS2, CNOT2, TRHDE regulates the metabolism in various tissues by a thyroid hormone signaling pathway.Conclusions: Synthesizes results from GWAS and transcriptome analyses, TPH2, GRIP1, FRS2, TRHDE, CNOT2 genes were considered as candidate genes for regulating FCR in Yorkshire pigs. These findings help to improve the understandings of the genetic mechanism of FCR and potentially optimize the design of breeding schemes.

Pigs, as one of the most common livestock species worldwide, are expected to have a fast growth rate and lower subcutaneous fatness but higher intramuscular fat (“marbling meat”). Nowadays, it is believed that not only host genetics but also its gut microbiomes can modulate farm animal phenotypes, however, many of the mechanisms remain elusive. We measured the body weight (BW), average daily gain (ADG), backfat thickness (BFT), and intramuscular fatness (IMF) of 91 Enshi pigs at 260 days of age, then genotyped each one individually using a 50K single nucleotide polymorphism array and performed 16S ribosomal RNA gene sequencing on 455 microbial samples from the jejunum, ileum, cecum, colon, and rectum. The microbial diversity showed notable spatial variation across the entire intestinal tract, with the cecum and colon having the highest α-diversity. The cecal and colonic microbiotas made greater contributions to BW and ADG and accounted for 22–37% of the phenotypic variance. The jejunal and cecal microbiotas contributed more (13–31%) to the BFT and IMF than the other segments. Finally, from cecum, colon, and jejunum, we identified eight microbial taxa that were significantly correlated with the target traits. The genera Alloprevotella and Ruminococcaceae UCG-005 were highly positively correlated with BW and ADG. The genera Prevotellaceae UCG-001 and Alistipes in the cecum and Clostridium sensu stricto 1 in the jejunum were highly positively correlated with BFT and IMF. The genera Stenotrophomonas, Sphaerochaeta, and Desulfovibrio were negatively associated with the mentioned traits. These findings could aid in developing strategies for manipulating the gut microbiota to alter production performance in pigs.

The analyses of multi-omics data have revealed candidate genes for objective traits. However, they are integrated poorly, especially in non-model organisms, and they pose a great challenge for prioritizing candidate genes for follow-up experimental verification. Here, we present a general convolutional neural network model that integrates multi-omics information to prioritize the candidate genes of objective traits. By applying this model to Sus scrofa, which is a non-model organism, but one of the most important livestock animals, the model precision was 72.9%, recall 73.5%, and F1-Measure 73.4%, demonstrating a good prediction performance compared with previous studies in Arabidopsis thaliana and Oryza sativa. Additionally, to facilitate the use of the model, we present ISwine (http://iswine.iomics.pro/), which is an online comprehensive knowledgebase in which we incorporated almost all the published swine multi-omics data. Overall, the results suggest that the deep learning strategy will greatly facilitate analyses of multi-omics integration in the future., Yuhua Fu et al. develop a CNN model that integrates multi-omics information to prioritize candidate genes of objective traits. Their model performs well when applied to important livestock non-model animals like Sus scrofa. Finally, the authors present ISwine, an online comprehensive knowledgebase which includes all published swine omics data to facilitate the integration of heterogeneous data.

African swine fever (ASF) is one of the most severe diseases of pigs. In this study, a CRISPR-Cas12a (also known as Cpf1) system coupled with nucleic acid amplification was optimized for the detection of ASF virus (ASFV). Two novel single-stranded DNA-fluorophore-quencher (ssDNA-FQ) reporters were developed to increase the brightness of the fluorescent signal for the visualization of nucleic acid detection. The CRISPR-Cas12a system was used to simultaneously cleave the polymerase chain reaction (PCR) or loop-mediated isothermal amplification (LAMP) amplicons and the newly developed ssDNA-FQ reporter, resulting in fluorescence that could be easily detected in multiple platforms, especially on cheap and portable blue or UV light transilluminators. This specific cleavage with fluorescence reveals the presence of the amplicon and confirms its identity, thereby preventing false-positive test results from nonspecific amplicons. This method is also uninterfered by the presence of large amounts of irrelevant background DNA and displays no cross-reactivity with other porcine DNA or RNA viruses. When coupled with LAMP, the Cas12a platform can detect a plasmid containing p72 with as few as 2 copies/μL reaction. Our results indicate that the CRISPR-Cas12a enhanced fluorescence assay coupled with nucleic acid amplification is robust, convenient, specific, confirmatory, affordable, and potentially adaptable for ASF diagnosis.

Colostrum quality is a vital factor in mortality and growth performance for piglets. Lactoferrin is an immuno-active milk protein that contributes to the formation of a protective layer above intestinal mucosa, possesses the antibacterial and antiviral activities that are favorable for piglet development. However, there is a notable reduction in lactoferrin in sow milk during lactation after the first few days, which causes many piglets to fail to ingest enough colostrum thereby leading to an increase in piglet mortality. In this study, we successfully constructed genome-edited Large-White pigs with marker-free site-specific knock-in of lactoferrin gene in the 3′-end of Casein alpha-s1 via CRISPR/Cas9 mediated homologous recombination. Thus, the lactoferrin protein can be expressed in the mammary gland in the control of Casein alpha-s1 promoter. As expected, the lactoferrin protein in genetically modified pigs sustained high expression in both colostrum and milk when compared with wild-type pigs. Moreover, the bacterial plate assay indicated that the milk from genetically modified pigs showed bacteriostatic effects when compared with control pigs. Taken together, our study demonstrated that the milk from genetically modified pigs had antibacterial activity which may reduce the costs of veterinary drug and improve the surviving rate of piglets, which is promising for pig breeding.

It is well known that the estrogen receptor alpha gene (ESR1) affects the reproductive traits of pigs; however, the immune role of ESR1 gene has not yet been resolved. Here, we characterized the pleiotropic aspects of ESR1 gene in immunity using the pig model. Tissue expression profile showed that the ESR1 gene had a broad ectopic expression in multiple reproductive and immune-related tissues/organs, which provided the tissue-level spatial fundamental of ESR1 gene that might function as a pleiotropic immune regulator. Using the peripheral blood cell model, a coupling transcriptome analytical strategy was proposed and verified that there existed strong positive or negative correlations of ESR1 gene with hundreds of differentially expressed genes that were involved in the immune regulation, indicating that the ESR1 gene might affect or be affected by, directly or indirectly, dozens of immune-related genes in the peripheral blood cells. Furthermore, the results of genetic association analysis showed that the SmaI-polymorphism of ESR1 gene had significant or highly significant associations with multiple immune traits, including platelet (PLT), hematocrit (HCT), the number of CD4-CD8-CD3- cells, plateletcrit (PCT), mean corpuscular volume (MCV), and mean corpuscular hemoglobin concentration (MCHC). Multiple evidences supported the immune pleiotropic roles of ESR1 gene in pigs. The study advances our understanding of the cross-species immune pleiotropic landscape of ESR1 gene and also provides a potential pleiotropic molecular marker for disease-resistant breeding in pigs.

Despite the broad variety of available microRNA (miRNA) research tools and methods, their application to the identification, annotation, and target prediction of miRNAs in nonmodel organisms is still limited. In this study, we collected nearly all public sRNA-seq data to improve the annotation for known miRNAs and identify novel miRNAs that have not been annotated in pigs (Sus scrofa). We newly annotated 210 mature sequences in known miRNAs and found that 43 of the known miRNA precursors were problematic due to redundant/missing annotations or incorrect sequences. We also predicted 811 novel miRNAs with high confidence, which was twice the current number of known miRNAs for pigs in miRBase. In addition, we proposed a correlation-based strategy to predict target genes for miRNAs by using a large amount of sRNA-seq and RNA-seq data. We found that the correlation-based strategy provided additional evidence of expression compared with traditional target prediction methods. The correlation-based strategy also identified the regulatory pairs that were controlled by nonbinding sites with a particular pattern, which provided abundant complementarity for studying the mechanism of miRNAs that regulate gene expression. In summary, our study improved the annotation of known miRNAs, identified a large number of novel miRNAs, and predicted target genes for all pig miRNAs by using massive public data. This large data-based strategy is also applicable for other nonmodel organisms with incomplete annotation information.

As an emerging swine enteropathogenic coronavirus, porcine deltacoronavirus (PDCoV) not only causes serious diarrhea in suckling piglets but also possesses the potential for cross-species transmission, which has sparked growing interest when studying this emerging virus. We previously identified a novel accessory protein NS7a encoded by PDCoV; however, the function of NS7a was not resolved. In this study, we demonstrated that PDCoV NS7a is an interferon antagonist. Overexpression of NS7a notably inhibited Sendai virus (SeV)-induced interferon-β (IFN-β) production and the activation of IRF3 rather than NF-κB. NS7a also inhibited IFN-β promoter activity induced by RIG-I, MDA5, MAVS, TBK1, and IKKe, which are key components of the RIG-I-like receptor (RLR) signaling pathway but not IRF3, the transcription factor downstream of TBK1/IKKe. Surprisingly, NS7a specifically interacts with IKKe but not with the closely related TBK1. Furthermore, NS7a interacts simultaneously with the kinase domain (KD) and the scaffold dimerization domain (SDD) of IKKe, competing with TRAF3, and IRF3 for binding to IKKe, leading to the reduction of RLR-mediated IFN-β production. The interactions of TRAF3-IKKe and IKKe-IRF3 are also attenuated in PDCoV-infected cells. Taken together, our results demonstrate that PDCoV NS7a inhibits IFN-β production by disrupting the association of IKKe with both TRAF3 and IRF3, revealing a new mechanism utilized by a PDCoV accessory protein to evade the host antiviral innate immune response.

Investigating the patterns of homozygosity, linkage disequilibrium, effective population size and inbreeding coefficients in livestock contributes to our understanding of the genetic diversity and evolutionary history. Here we used Illumina PorcineSNP50 Bead Chip to identify the runs of homozygosity (ROH) and estimate the linkage disequilibrium (LD) across the whole genome, and then predict the effective population size. In addition, we calculated the inbreeding coefficients based on ROH in 305 Pi&eacute, train pigs and compared its effect with the other two types of inbreeding coefficients obtained by different calculation methods. A total of 23,434 ROHs were detected, and the average length of ROH per individual was about 507.27 Mb. There was no regularity on how those runs of homozygosity distributed in genome. The comparisons of different categories suggested that the formation of long ROH was probably related with recent inbreeding events. Although the density of genes located in ROH core regions is lower than that in the other genomic regions, most of them are related with Pi&eacute, train commercial traits like meat qualities. Overall, the results provide insight into the way in which ROH is produced and the identified ROH core regions can be used to map the genes associated with commercial traits in domestic animals.

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

Birth weight of pigs is an important economic factor in the livestock industry. The identification of the genes and variants that underlie birth weight is of great importance. In this study, we integrated two genotyping methods, single nucleotide polymorphism (SNP) chip analysis and restriction site associated DNA sequencing (RAD-seq) to genotype genome-wide SNPs. In total, 45,175 and 139,634 SNPs were detected with the SNP chip and RAD-seq, respectively. The genome-wide association study (GWAS) of the combined SNP panels identified two significant loci located at chr1: 97,745,041 and chr4: 112,031,589, that explained 6.36% and 4.25% of the phenotypic variance respectively. To reduce interval containing causal variants, we imputed sequence-level SNPs in the GWAS identified regions and fine-mapped the causative variants into two narrower genomic intervals: a ∼100 kb interval containing 71 SNPs and a broader ∼870 kb interval with 432 SNPs. This fine-mapping highlighted four promising candidate genes, SKOR2, SMAD2, VAV3, and NTNG1. Additionally, the functional genes, SLC25A24, PRMT6 and STXBP3, are also located near the fine-mapping region. These results suggest that these candidate genes may have contribute substantially to the birth weight of pigs.

miRNAs reportedly participate in various biological processes, such as skeletal muscle proliferation and differentiation. However, the regulation of differentially expressed (DE) miRNAs and their function in myogenesis remain unclear. Herein, miRNA expression profiles and regulation during C2C12 differentiation were analyzed in relation to chromatin states by RNA-seq, ATAC-seq, and ChIP-seq. We identified 19 known and nine novel differentially expressed miRNAs at days 0, 1, 2, and 4. The expression of the differentially expressed miRNAs was related to the chromatin states of the 113 surrounding open chromatin regions defined by ATAC-seq peaks. Of these open chromatin regions, 44.25% were colocalized with MyoD/MyoG binding sites. The remainder of the above open chromatin regions were enriched with motifs of the myoblast-expressed AP-1 family, Ctcf, and Bach2 transcription factors (TFs). Additionally, the target genes of the above differentially expressed miRNAs were enriched primarily in muscle growth and development pathways, especially the Hippo signaling pathway. Moreover, via combining a loss-of-function assay with Q-PCR, western blotting, and immunofluorescence, we confirmed that the Hippo signaling pathway was responsible for C2C12 myoblast differentiation. Thus, our results showed that these differentially expressed miRNAs were regulated by chromatin states and affected muscle differentiation through the Hippo signaling pathway. Our findings provide new insights into the function of these differentially expressed miRNAs and the regulation of their expression during myoblast differentiation.

Objective This study aims to clarify the mechanisms underlying transcriptional regulation and regulatory roles of lncRNAs in skeletal muscle cell differentiation. Methods We analysed the expression patterns of lncRNAs via time‐course RNA‐seq. Then, we further combined the ATAC‐seq and ChIP‐seq to investigate the governing mechanisms of transcriptional regulation of differentially expressed (DE) lncRNAs. Weighted correlation network analysis and GO analysis were conducted to identify the transcription factor (TF)‐lncRNA pairs related to skeletal muscle cell differentiation. Results We identified 385 DE lncRNAs during C2C12 differentiation, the transcription of which is determined by chromatin states around their transcriptional start sites. The TF‐lncRNA correlation network showed substantially concordant changes in DE lncRNAs between C2C12 differentiation and satellite cell rapid growth stages. Moreover, the up‐regulated lncRNAs showed a significant decrease following the differentiation capacity of satellite cells, which gradually declines during skeletal muscle development. Notably, inhibition of the lncRNA Atcayos and Trp53cor1 led to the delayed differentiation of satellite cells. Those lncRNAs were significantly up‐regulated during the rapid growth stage of satellite cells (4‐6 weeks) and down‐regulated with reduced differentiation capacity (8‐12 weeks). It confirms that these lncRNAs are positively associated with myogenic differentiation of satellite cells during skeletal muscle development. Conclusions This study extends the understanding of mechanisms governing transcriptional regulation of lncRNAs and provides a foundation for exploring their functions in skeletal muscle cell differentiation., Schematic of TF (MyoD, MyoG, Mef2c and Heyl) and lncRNA (Atcayos and Trp53cor1) regulatory interactions in differentiating C2C12 and ageing satellite cells.

Embryonic and neonatal skeletal muscles grow via the proliferation and fusion of myogenic cells, whereas adult skeletal muscle adapts largely by remodelling pre-existing myofibers and optimizing metabolic balance. It has been reported that miRNAs played key roles during skeletal muscle development through targeting different genes at post-transcriptional level. In this study, we show that a single miRNA (miR-208b) can modulate both the myogenesis and homoeostasis of skeletal muscle by distinct targets. As results, miR-208b accelerates the proliferation and inhibits the differentiation of myogenic cells by targeting the E-protein family member transcription factor 12 (TCF12). Also, miR-208b can stimulate fast-to-slow fibre conversion and oxidative metabolism programme through targeting folliculin interacting protein 1 (FNIP1) but not TCF12 gene. Further, miR-208b could active the AMPK/PGC-1a signalling and mitochondrial biogenesis through targeting FNIP1. Thus, miR-208b could mediate skeletal muscle development and homoeostasis through specifically targeting of TCF12 and FNIP1.

Somatic cloning has had a significant impact on the life sciences and is important in a variety of processes, including medical research and animal production. However, the application of somatic cloning has been limited due to its low success rate. Therefore, potential epigenetic variations between cloned and donor animals are still unclear. DNA methylation, one of the factors which is responsible for phenotypic differences in animals, is a commonly researched topic in epigenetic studies of mammals. To investigate the epigenetic variations between cloned and donor animals, we selected blood and ear fibroblasts of a donor pig and a cloned pig to perform whole-genome bisulfite sequencing (WGBS). A total of 215 and 707 differential methylation genes (DMGs) were identified in blood and ear fibroblasts, respectively. Functional annotation revealed that DMGs are enriched in many pathways, including T/B or natural killer (NK) cell differentiation, oocyte maturation, embryonic development, and reproductive hormone secretion. Moreover, 22 DMGs in the blood and 75 in the ear were associated with immune responses (e.g., CD244, CDK6, CD5, CD2, CD83, and CDC7). We also found that 18 DMGs in blood and 53 in ear fibroblasts were involved in reproduction. Understanding the expression patterns of DMGs, especially in relation to immune responses and reproduction, will reveal insights that will aid the advancement of future somatic cloning techniques in swine.

Improving meat quality is a crucial purpose of commercial production and breeding systems. In this study, multiomics techniques were used to investigate the molecular mechanisms that impact the excessive diversity of meat quality in Enshi black pigs. The results suggest that 120 differentially expressed genes (DEGs) and 171 significantly changed metabolites (SCMs) contribute to the content of intramuscular fat (IMF) through the processes of fat accumulation and regulation of lipolysis. A total of 141 DEGs and 47 SCMs may regulate meat color through the processes of nicotinate and nicotinamide metabolism. Herein, we found some candidate genes associated with IMF and meat color. We also presented a series of metabolites that are potentially available biological indicators to measure meat quality. This research provides further insight into the detection of intramuscular fat accumulation and meat color variation and provides a reference for molecular mechanisms in the regulation of IMF and meat color.

Prostate cancer (PCa) has a predominantly luminal phenotype. Basal cells were previously identified as a cell of origin for PCa, but increasing evidence implicates luminal cells as a preferred cell of origin for PCa, as well as key drivers of tumor development and progression. Prostate luminal cells are understudied compared with basal cells. In this review, we describe the contribution of prostate luminal progenitor (LP) cells to luminal cell development and their role in prostate development, androgen-mediated regeneration of castrated prostate, and tumorigenesis. We also discuss the potential value of LP transcriptomics to identify new targets and therapies to treat aggressive PCa. Finally, we propose future research directions focusing on molecular mechanisms underlying LP cell biology and heterogeneity in normal and diseased prostate.

Duration of fertility, (DF) is an important functional trait in poultry production and lncRNAs have emerged as important regulators of various process including fertility. In this study we applied a genome-guided strategy to reconstruct the uterovaginal junction (UVJ) transcriptome of 14 egg-laying birds with long- and short-DF (n = 7); and sought to uncover key lncRNAs related to duration of fertility traits by RNA-sequencing technology. Examination of RNA-seq data revealed a total of 9977 lncRNAs including 2576 novel lncRNAs. Differential expression (DE) analysis of lncRNA identified 223 lncRNAs differentially expressed between the two groups. DE-lncRNA target genes prediction uncovered over 200 lncRNA target genes and functional enrichment tests predict a potential function of DE-lncRNAs. Gene ontology classification and pathway analysis revealed 8 DE-lncRNAs, with the majority of their target genes enriched in biological functions such as reproductive structure development, developmental process involved in reproduction, response to cytokine, carbohydrate binding, chromatin organization, and immune pathways. Differential expression of lncRNAs and target genes were confirmed by qPCR. Together, these results significantly expand the utility of the UVJ transcriptome and our analysis identification of key lncRNAs and their target genes regulating DF will form the baseline for understanding the molecular functions of lncRNAs regulating DF.

Alternative splicing (AS) and fusion transcripts produce a vast expansion of transcriptomes and proteomes diversity. However, the reliability of these events and the extend of epigenetic mechanisms have not been adequately addressed due to its limitation of uncertainties about the complete structure of mRNA. Here we combined single-molecule real-time sequencing, Illumina RNA-seq and DNA methylation data to characterize the landscapes of DNA methylation on AS, fusion isoforms formation and lncRNA feature and further to unveil the transcriptome complexity of pig. Our analysis identified an unprecedented scale of high-quality full-length isoforms with over 28,127 novel isoforms from 26,881 novel genes. More than 92,000 novel AS events were detected and intron retention predominated in AS model, followed by exon skipping. Interestingly, we found that DNA methylation played an important role in generating various AS isoforms by regulating splicing sites, promoter regions and first exons. Furthermore, we identified a large of fusion transcripts and novel lncRNAs, and found that DNA methylation of the promoter and gene body could regulate lncRNA expression. Our results significantly improved existed gene models of pig and unveiled that pig AS and epigenetic modify were more complex than previously thought.

Up to date, the scarcity of publicly available complete mitochondrial sequences for European wild pigs hampers deeper understanding about the genetic changes following domestication. Here, we have assembled 26 de novo mtDNA sequences of European wild boars from next generation sequencing (NGS) data and downloaded 174 complete mtDNA sequences to assess the genetic relationship, nucleotide diversity, and selection. The Bayesian consensus tree reveals the clear divergence between the European and Asian clade and a very small portion (10 out of 200 samples) of maternal introgression. The overall nucleotides diversities of the mtDNA sequences have been reduced following domestication. Interestingly, the selection efficiencies in both European and Asian domestic pigs are reduced, probably caused by changes in both selection constraints and maternal population size following domestication. This study suggests that de novo assembled mitogenomes can be a great boon to uncover the genetic turnover following domestication. Further investigation is warranted to include more samples from the ever-increasing amounts of NGS data to help us to better understand the process of domestication.

Aim: This study aims to couple DNA methylation changes and evolution of retrogenes. Materials & methods: A new two-step strategy was developed to screen retrogenes. Further, reduced representation bisulfite sequencing and RNA-seq data of eight tissues were used to analyze retrogenes. Results: A total of 964 retrocopies were identified and new retrocopies were available for the synthesis of glycans and lipids corresponding to pig phenotypic traits. Retrogenes were consistently hypermethylated. Hypomethylation of parental genes presented more susceptibility to retroposition. Promoter DNA methylation of retrogenes was negatively correlated with evolutionary time and played important roles in regulating retrogene tissue-specific expression pattern. Conclusion: A two-step procedure is effective and necessary for identifying retrogenes. DNA methylation drives origination, survival, evolution and expression of retrogenes.

Lipid metabolism has important roles in animal growth, development, and reproduction. As a regulator of lipid metabolism, CIDEc promotes unilocular development of lipid droplets and stimulates intracellular lipid deposition, and has two isoforms, CIDEc-l and CIDEc-s. CIDEc-l has ten more N-terminal amino acids than CIDEc-s. However, the functions of two isoforms are largely unknown. In this study, the expression profiles of two isoforms in Bama pigs differed, with cidec-l dominant in the liver and small intestine, and cidec-s dominant in muscle and adipose tissue. Fasting and consuming a high-fat diet resulted in changes in the expression of the two isoforms that were closely related to changes in blood and muscle triglyceride (TG) concentrations. Comparison of gene expression and TG concentration suggested that CIDEc-l accelerated lipid secretion and that CIDEc-s promoted lipid deposition, implying that the two isoforms had different functions. Study In vitro confirmed that CIDEc-s stimulated lipid deposition in C2C12 muscle cells and CIDEc-l promoted lipid secretion in HepG2 liver cells. The results showed that two tissue-specific CIDEc isoforms had different roles in lipid deposition and secretion. They may be potential targets for regulation of fat content.

It is one of the important events that trophoblast cells within the placental folds differentiate into two types that differ in cell shape during placental development in pigs. This study showed that all the trophoblast cells were of similar shape between Yorkshire and Chinese Meishan pigs on day 26 of gestation; thereafter, the trophoblast cells located at the top of the placental folds became high columnar, while those cells at the base of the placental folds were cuboidal on day 50 of gestation. Additionally, on day 95 of gestation, all the trophoblast cells in Meishan pigs became cuboidal, but the trophoblast cells located at the top of the placental folds in Yorkshire pigs still remained columnar. The membranous E-cadherin and β-catenin were strongly co-expressed by the high columnar trophoblast cells but very weakly expressed by those cuboidal cells. Consistently, the expression pattern of ZEB2, the E-cadherin repressor, was inversely correlated with that of E-cadherin in the two types of trophoblast cells in the two breeds. Furthermore, electrophoretic mobility shift assays demonstrated the binding of ZEB2 to the E-cadherin promoter in nuclear extracts from porcine placental tissue. These findings suggest a ZEB2-dependent mechanism of trophoblast cell differentiation during placental development in pigs.

Intensive oxidative stress occurs during high-fat-diet-induced hepatic fat deposition, suggesting a critical role for redox signaling in liver metabolism. Intriguingly, evidence shows that fasting could also result in redox-profile changes largely through reduced oxidant or increased antioxidant levels. However, a comprehensive landscape of redox-modified hepatic substrates is lacking, thereby hindering our understanding of liver metabolic homeostasis. We employed a proteomic approach combining iodoacetyl tandem mass tag and nanoliquid chromatography tandem mass spectrometry to quantitatively probe the effects of high-fat feeding and fasting on in vivo redox-based cysteine modifications. Compared with control groups, ∼60% of cysteine residues exhibited downregulated oxidation ratios by fasting, whereas ∼94% of these ratios were upregulated by high-fat feeding. Importantly, in fasted livers, proteins exhibiting diminished cysteine oxidation were annotated in pathways associated with fatty acid metabolism, carbohydrate metabolism, insulin, peroxisome proliferator-activated receptors, and oxidative respiratory chain signaling, suggesting that fasting-induced redox changes targeted major metabolic pathways and consequently resulted in hepatic lipid accumulation.

Neutrophils are vital components of defense mechanisms against invading pathogens and are closely linked with the individual antiviral capacity of pigs and other mammals. Neutrophilia is a well-known clinical characteristic of viral and bacterial infections. Using Affymetrix porcine genome microarrays, we investigated the gene expression profiles associated with neutrophil variation in porcine peripheral blood before and after polyriboinosinic-polyribocytidylic acid stimulation. Transcriptomic analysis showed 796 differentially expressed genes (DEGs) in extreme response (ER) pigs and 192 DEGs in moderate response (MR) pigs. Most DEGs were related to immune responses, included MXD1, CXCR4, CREG1, MyD88, CD14, TLR2, TLR4, IRF3 and IRF7. Gene ontology analysis indicated that the DEGs of both ER and MR pigs were involved in common biological processes, such as cell proliferation, growth regulation, immune response, inflammatory response and cell activation. The ER and MR groups also showed differences in DEGs involved in biological processes. DEGs involved in cell division and cell cycle were specifically found in the ER pigs, whereas DEGs involved in cell migration were specifically found in the MR pigs. The study provides a basic understanding of the molecular basis for the antiviral capacity of pigs and other mammals.

An important characteristic of intrauterine growth restricted (IUGR) neonate is the impaired intestinal barrier function. With the use of a pig model, this study was conducted to identify the responsible microRNA (miRNA) for the intestinal damage in IUGR neonates through comparing the miRNA profile of IUGR and normal porcine neonates and to investigate the regulation mechanism. Compared with the normal ones, we identified 83 upregulated and 76 downregulated miRNAs in the jejunum of IUGR pigs. Notably, IUGR is associated with profoundly increasesd miR-29 family and decreased expression of extracellular matrix (ECM) and tight junction (TJ) proteins in the jejunum. Furthermore, in vitro study using theporcine intestinal epithelial cell line (IPEC-1) showed that inhibition of miR-29a expression could improve the monolayer integrity by increasing cell proliferation and transepithelial resistance. Also, overexpression/inhibition of miR-29a in IPEC-1 cells can suppress/increase the expression of integrin-β1, collagen I, collagen IV, fibronectin, and claudin 1, both at transcriptional and translational levels. Subsequent luciferase reporter assay confirmed a direct interaction between miR-29a and the 3′-untranslated regions of these genes. In conclusion, this study reveals that IUGR-impaired intestinal barrier function is associated with downregulated ECM and TJ protein expression mediated by the upregulation of miR-29a. NEW & NOTEWORTHY Intrauterine growth restricted (IUGR) remains a major problem for both human health and animal production due to its association with high rates of preweaning morbidity and mortality. We have identified the abnormal expression of microRNA-29a (miR-29a) in the small intestine of IUGR neonates, as well as its targets and mechanisms. These results provide new information about biological characteristics of IUGR-affected intestinal dysfunction and can lead to the development of potentially solution for preventing and treating IUGR in the future.

Haemonchus contortus is a major gastrointestinal nematode affecting goat in China and Bangladesh like many countries in the world. The aim of this study was to evaluate the susceptibility of Haemonchus contortus infection in different goat breeds/populations in China and Bangladesh under the natural grazing condition. Fecal Egg of Haemonchus contortus was counted from 430 goats in both countries using McMaster’s technique. Body weight, Pack Cell Volume and Hemoglobin value were measured for Bangladeshi goats. Statistical analysis was done after transforming data into log10 (n+1), where n is the number of egg per gram feces and analysis of variance was done by using Generalized Linear Model procedures of computer package of SAS. Prevalence of parasite and Fecal Egg Count varied significantly (P0.05). Black Bengal goat in the hilly region of Bangladesh was the least susceptible to Haemonchus contortus infection in term of Fecal Egg Count. Our study suggests that Chinese goat breeds were more susceptible to Haemonchus contortus infection as compared to Black Bengal goat of Bangladesh under natural grazing condition. Further studies on immune-genetics aspects of goats might be helpful to find out the actual causes of such differences and to develop disease resistant breed of goat.Progressive Agriculture 27 (4): 473-481, 2016

Determination of genetic diversity and population structure plays an important role in supporting genetic improvement programs and future conservation plans. In this study, 352 individuals representing nine Chinese indigenous goat populations distributed in China were genotyped at 15 microsatellite loci. The mean number of alleles (MNA) per breed ranged from 3.867 (Matou goat, MT) to 5.400 (Ujumqin white goat, UW), the expected heterozygosity (He) varied from 0.482 (Hainan black goat, HNB) to 0.659 (UW). Allelic diversity and heterozygosity measures in the studied populations were much lower than Chinese Cashmere and meat type goats. Global F statistics revealed 9.8% of total variance explained among breeds while 90.2% of variance was due to diversities within breeds. Three major population clusters were observed broadly conforming to geographical locations of different goat populations. Three-dimensional scatterplot derived from three largest principal components supported the observed phylogeny based on genetic distance estimations. Goats from Northern China and Island region were distinct while strong admixture was observed among goat populations from Central China. The study revealed market orientation and geographical distances among populations might have contributed to the genetic structure and population sub-division among Chinese indigenous goats. Our study provided a new insight into understanding the genetic diversity and structure of Chinese indigenous goat breeds, and will be helpful to determine the strategies for breeding and conservation programs.

Designing efficient and specific CRISPR single-guide RNAs (sgRNAs) is vital for the successful application of CRISPR technology. Currently, a growing number of new RNA-guided endonucleases with a different protospacer adjacent motif (PAM) have been discovered, suggesting the necessity to develop a versatile tool for designing sgRNA to meet the requirement of different RNA-guided DNA endonucleases. Here, we report the development of a flexible sgRNA design program named "CRISPR-offinder". Support for user-defined PAM and sgRNA length was provided to increase the targeting range and specificity. Additionally, evaluation of on- and off-target scoring algorithms was integrated into the CRISPR-offinder. The CRISPR-offinder has provided the bench biologist a rapid and efficient tool for identification of high quality target sites, and it is freely available at https://sourceforge.net/projects/crispr-offinder-v1-2/ or http://www.biootools.com.

Japanese encephalitis virus (JEV) is a mosquito-borne zoonotic flavivirus that causes encephalitis and reproductive disorders in mammalian species. However, the host factors critical for its entry, replication, and assembly are poorly understood. Here, we design a porcine genome-scale CRISPR/Cas9 knockout (PigGeCKO) library containing 85,674 single guide RNAs targeting 17,743 protein-coding genes, 11,053 long ncRNAs, and 551 microRNAs. Subsequently, we use the PigGeCKO library to identify key host factors facilitating JEV infection in porcine cells. Several previously unreported genes required for JEV infection are highly enriched post-JEV selection. We conduct follow-up studies to verify the dependency of JEV on these genes, and identify functional contributions for six of the many candidate JEV-related host genes, including EMC3 and CALR. Additionally, we identify that four genes associated with heparan sulfate proteoglycans (HSPGs) metabolism, specifically those responsible for HSPGs sulfurylation, facilitate JEV entry into porcine cells. Thus, beyond our development of the largest CRISPR-based functional genomic screening platform for pig research to date, this study identifies multiple potentially vulnerable targets for the development of medical and breeding technologies to treat and prevent diseases caused by JEV., Here the authors report the construction of a genome-scale porcine CRISPR/Cas9 library, called PigGeCKO, for screening and analyses of host resistance genes and factors associated with Japanese encephalitis virus replication.

Japanese encephalitis virus (JEV) is a mosquito-borne zoonotic flavivirus that causes encephalitis and reproductive disorders in mammalian species. However, key host genes involved in the JEV life cycle and cell death are poorly understood. Here, we designed 85,674 single guide RNAs (sgRNAs) targeting 17,743 protein-coding genes, 11,053 long ncRNAs (lncRNAs), and 551 microRNAs (miRNAs) in the porcine genome, and subsequently developed a porcine sgRNA library and genome-scale CRISPR/Cas9 knockout (PigGeCKO) system. These sgRNAs were delivered into porcine kidney-15 (PK-15) cells that constitutively express Cas9, positive selection screening of the resulting PigGeCKO cell collection for resistance to JEV-induced cell death led to the identification of several previously unreported genes required for JEV infection. We conducted follow-up studies to verify the dependency of JEV on these genes, and identified functional contributions for six of the many candidate JEV-related host genes, including an endoplasmic reticulum membrane protein complex subunit 3 (EMC3) and calreticulin (CALR). Additionally, we identified that four genes associated with heparan sulfate proteoglycans (HSPGs) metabolism, specifically those responsible for HSPG sulfurylation, facilitated JEV entry into PK-15 cells. Thus, beyond our development of the largest CRISPR-based functional genomic screening platform for pig research to date, this work deepens our basic understanding of flavivirus infection and identifies multiple potentially vulnerable targets for the development of medical and breeding technologies to prevent and treat diseases caused by Japanese encephalitis virus.