An integrated in silico approach for functional and structural impact of non- synonymous SNPs in the MYH1 gene in Jeju Native Pigs

Comparison of differentially expressed genes in longissimus dorsi muscle of Diannan small ears, Wujin and landrace pigs using RNA-seq

Introduction

Pig growth is an important economic trait that involves the co-regulation of multiple genes and related signaling pathways. High-throughput sequencing has become a powerful technology for establishing the transcriptome profiles and can be used to screen genome-wide differentially expressed genes (DEGs). In order to elucidate the molecular mechanism underlying muscle growth, this study adopted RNA sequencing (RNA-seq) to identify and compare DEGs at the genetic level in the longissimus dorsi muscle (LDM) between two indigenous Chinese pig breeds (Diannan small ears [DSE] pig and Wujin pig [WJ]) and one introduced pig breed (Landrace pig [LP]).

Methods

Animals under study were from two Chinese indigenous pig breeds (DSE pig, n = 3; WJ pig, n = 3) and one introduced pig breed (LP, n = 3) were used for RNA sequencing (RNA-seq) to identify and compare the expression levels of DEGs in the LDM. Then, functional annotation, Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, and Protein-Protein Interaction (PPI) network analysis were performed on these DEGs. Then, functional annotation, Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, and Protein-Protein Interaction (PPI) network analysis were performed on these DEGs.

Results

The results revealed that for the DSE, WJ, and LP libraries, more than 66, 65, and 71 million clean reads were generated by transcriptome sequencing, respectively. A total of 11,213 genes were identified in the LDM tissue of these pig breeds, of which 7,127 were co-expressed in the muscle tissue of the three samples. In total, 441 and 339 DEGs were identified between DSE vs. WJ and LP vs. DSE in the study, with 254, 193 up-regulated genes and 187, 193 down-regulated genes in DSE compared to WJ and LP. GO analysis and KEGG signaling pathway analysis showed that DEGs are significantly related to contractile fiber, sarcolemma, and dystrophin-associated glycoprotein complex, myofibril, sarcolemma, and myosin II complex, Glycolysis/Gluconeogenesis, Propanoate metabolism, and Pyruvate metabolism, etc. In combination with functional annotation of DEGs, key genes such as ENO3 and JUN were identified by PPI network analysis.

Discussion

In conclusion, the present study revealed key genes including DES, FLNC, PSMD1, PSMD6, PSME4, PSMB4, RPL11, RPL13A, ROS23, RPS29, MYH1, MYL9, MYL12B, TPM1, TPM4, ENO3, PGK1, PKM2, GPI, and the unannotated new gene ENSSSCG00000020769 and related signaling pathways that influence the difference in muscle growth and could provide a theoretical basis for improving pig muscle growth traits in the future.

In Ovo Injection of CHIR-99021 Promotes Feather Follicle Development via Modulating the Wnt Signaling Pathway and Transcriptome in Goose Embryos (Anser cygnoides)

Feather performs important physiological functions in birds, and it is also one of the economic productions in goose farming. Understanding and modulating feather follicle development during embryogenesis are essential for bird biology and the poultry industry. CHIR-99021 is a potent Wnt/β-catenin signaling pathway activator associated with feather follicle development. In this study, goose embryos (Anser cygnoides) received an in ovo injection of CHIR-9902, which was conducted at the beginning of feather follicle development (E9). The results showed that feather growth and feather follicle development were promoted. The Wnt signaling pathway was activated by the inhibition of GSK-3β. Transcriptomic analyses showed that the transcription changes were related to translation, metabolism, energy transport, and stress in dorsal tissue of embryos that received CHIR-99021, which might be to adapt and coordinate the promoting effects of CHIR-99021 on feather follicle development. This study suggests that in ovo injection of CHIR-99021 is a potential strategy to improve feather follicle development and feather-related traits for goose farming and provides profiling of the Wnt signaling pathway and transcriptome in dorsal tissue of goose embryos for further understanding of feather follicle development.

Integration of multi-omics approaches for functional characterization of muscle related selective sweep genes in Nanchukmacdon

Pig as a food source serves daily dietary demand to a wide population around the world. Preference of meat depends on various factors with muscle play the central role. In this regards, selective breeding abled us to develop "Nanchukmacdon" a pig breeds with an enhanced variety of meat and high fertility rate. To identify genomic regions under selection we performed whole-genome resequencing, transcriptome, and whole-genome bisulfite sequencing from Nanchukmacdon muscles samples and used published data for three other breeds such as Landrace, Duroc, Jeju native pig and analyzed the functional characterization of candidate genes. In this study, we present a comprehensive approach to identify candidate genes by using multi-omics approaches. We performed two different methods XP-EHH, XP-CLR to identify traces of artificial selection for traits of economic importance. Moreover, RNAseq analysis was done to identify differentially expressed genes in the crossed breed population. Several genes (UGT8, ZGRF1, NDUFA10, EBF3, ELN, UBE2L6, NCALD, MELK, SERP2, GDPD5, and FHL2) were identified as selective sweep and differentially expressed in muscles related pathways. Furthermore, nucleotide diversity analysis revealed low genetic diversity in Nanchukmacdon for identified genes in comparison to related breeds and whole-genome bisulfite sequencing data shows the critical role of DNA methylation pattern in identified genes that leads to enhanced variety of meat. This work demonstrates a way to identify the molecular signature and lays a foundation for future genomic enabled pig breeding.

Growth factors improve the proliferation of Jeju black pig muscle cells by regulating myogenic differentiation 1 and growth-related genes

OBJECTIVE

The growth rate of pigs is related to differentiation and proliferation of muscle cells, which are regulated by growth factors and expression of growth-related genes. Thus, the objective of this study was to establish optimal culture conditions for Jeju black pig (JBP) muscle cells and determine the relationship of various factors involved in muscle growth with the proliferation of JBP muscle cells.

METHODS

Muscles were taken from the femur skeletal muscle of JBP embryos. After isolation of the muscle cells, cells were cultured in a 6-well plate under four different culture conditions to optimize culture conditions for JBP muscle cells. To analyze proliferation rate of JBP muscle cells, these muscle cells were seeded into 6-well plates at a density of 1.5×105 cells per well and cultured for 3 days. Western blot and quantitative real-time polymerase chain reaction were applied to verify the myogenic differentiation 1 (MyoD) expression and growth-related gene expression in JBP muscle cells, respectively.

RESULTS

We established a muscle cell line from JBP embryos and optimized its culture conditions. These muscle cells were positive for MyoD, but not for paired box 7. The proliferation rate of these muscle cells was significantly higher in a culture medium containing bFGF and epidermal growth factor + basic fibroblast growth factor (EGF+bFGF) than that without a growth factor or containing EGF alone. Treatment with EGF and bFGF significantly induced the expression of MyoD protein, an important transcription factor in muscle cells. Moreover, we checked the changes of expression of growth-related genes in JBP muscle cells by presence or absence of growth factors. Expression level of collagen type XXI alpha 1 gene was changed only when EGF and bFGF were added together to culture media for JBP muscle cells.

CONCLUSION

Concurrent use of EGF and bFGF increased the expression of MyoD protein, thus regulating the proliferation of JBP muscle cells and the expression of growth-related genes.

Genomic approaches to identify hybrids and estimate admixture times in European wildcat populations

The survival of indigenous European wildcat (Felis silvestris silvestris) populations can be locally threatened by introgressive hybridization with free-ranging domestic cats. Identifying pure wildcats and investigating the ancestry of admixed individuals becomes thus a conservation priority. We analyzed 63k cat Single Nucleotide Polymorphisms (SNPs) with multivariate, Bayesian and gene-search tools to better evaluate admixture levels between domestic and wild cats collected in Europe, timing and ancestry proportions of their hybrids and backcrosses, and track the origin (wild or domestic) of the genomic blocks carried by admixed cats, also looking for possible deviations from neutrality in their inheritance patterns. Small domestic ancestry blocks were detected in the genomes of most admixed cats, which likely originated from hybridization events occurring from 6 to 22 generations in the past. We identified about 1,900 outlier coding genes with excess of wild or domestic ancestry compared to random expectations in the admixed individuals. More than 600 outlier genes were significantly enriched for Gene Ontology (GO) categories mainly related to social behavior, functional and metabolic adaptive processes (wild-like genes), involved in cognition and neural crest development (domestic-like genes), or associated with immune system functions and lipid metabolism (parental-like genes). These kinds of genomic ancestry analyses could be reliably applied to unravel the admixture dynamics in European wildcats, as well as in other hybridizing populations, in order to design more efficient conservation plans.

The first comprehensive description of the expression profile of genes involved in differential body growth and the immune system of the Jeju Native Pig and miniature pig

Sus scrofa provides a major source of animal protein for humans as well as being an excellent biomedical model. This study was carried out to understand, in detail, the genetic and functional variants of Jeju Native Pigs and miniature pigs through differential expression profiling of the genes controlling their immune response, growth performance, and meat quality. The Illumina HiSeq 2000 platform was used for generating 1.3 billion 90 bp paired-end reads, which were mapped to the S. scrofa genome using TopHat2. A total of 2481 and 2768 genes were differentially expressed with 8-log changes in muscle and liver samples, respectively. Five hundred forty-eight genes in muscle and 642 genes in liver samples had BLAST matches within the non-redundant database. GO process and pathway analyses showed enhanced biological processes related to the extracellular structural organization and skeletal muscle cell differentiation in muscle tissue, whereas the liver tissue shares functions related to the inflammatory response. Herein, we identify inflammatory regulatory genes in miniature pigs and growth response genes in Jeju Native Pigs, information which can provide a stronger base for the selection of breeding stock and facilitate further in vitro and in vivo studies for therapeutic purposes.

In silico Design of Laccase Thermostable Mutants From Lacc 6 of Pleurotus Ostreatus

Fungal laccase enzymes have a great biotechnological potential for bioremediation processes due to their ability to degrade compounds such as ρ-diphenol, aminophenols, polyphenols, polyamines, and aryldiamines. These enzymes have activity at different pH and temperature values, however, high temperatures can cause partial or total loss of enzymatic activity, so it is appropriate to do research to modify their secondary and/or tertiary structure to make them more resistant to extreme temperature conditions. In silico, a structure of the Lacc 6 enzyme of Pleurotus ostreatus was constructed using a laccase of Trametes versicolor as a template. From this structure, 16 mutants with possible resistance at high temperature due to ionic interactions, salt bridges and disulfide bonds were also obtained in silico. It was determined that 12 mutants called 4-DB, 3-DB, D233C-T310C, F468P, 3-SB, L132T, N79D, N372D, P203C, P203V, T147E, and W85F, presented the lowest thermodynamic energy. Based on the previous criterion and determining the least flexibility in the protein structures, three mutants (4-DB, 3-DB, and P203C) were selected, which may present high stability at high temperatures without affecting their active site. The obtained results allow the understanding of the molecular base that increase the structural stability of the enzyme Lacc 6 of Pleurotus ostreatus, achieving the in silico generation of mutants, which could have activity at high temperatures.

Transcriptome profile analysis of leg muscle tissues between slow- and fast-growing chickens

Chicken is widely favored by consumers because of some unique features. The leg muscles occupy an important position in the market. However, the specific mechanism for regulating muscle growth speed is not clear. In this experiment, we used Jinghai yellow chickens with different body weights at 300 days as research subjects. The chickens were divided into fast- and slow-growing groups, and we collected leg muscles after slaughtering for use in RNA-seq. After comparing the two groups, 87 differentially expressed genes (DEGs) were identified (fold change ≥ 2 and FDR < 0.05). The fast-growing group had 42 up-regulated genes and 45 down-regulated genes among these DEGs compared to the slow-growing group. Six items were significantly enriched in the biological process: embryo development ending in birth or egg hatching, chordate embryonic development, embryonic skeletal system development, and embryo development as well as responses to ketones and the sulfur compound biosynthetic process. Two significantly enriched pathways were found in the KEGG pathway analysis (P-value < 0.05): the insulin signaling pathway and the adipocytokine signaling pathway. This study provides a theoretical basis for the molecular mechanism of chicken growth and for improving the production of Jinghai yellow chicken.

Genome analysis of Yucatan miniature pigs to assess their potential as biomedical model animals

Objective

Pigs share many physiological, anatomical and genomic similarities with humans, which make them suitable models for biomedical researches. Understanding the genetic status of Yucatan miniature pigs (YMPs) and their association with human diseases will help to assess their potential as biomedical model animals. This study was performed to identify non-synonymous single nucleotide polymorphisms (nsSNPs) in selective sweep regions of the genome of YMPs and present the genetic nsSNP distributions that are potentially associated with disease occurrence in humans.

Methods

nsSNPs in whole genome resequencing data from 12 YMPs were identified and annotated to predict their possible effects on protein function. Sorting intolerant from tolerant (SIFT) and polymorphism phenotyping v2 analyses were used, and gene ontology (GO) network and Kyoto encyclopedia of genes and genomes (KEGG) pathway analyses were performed.

Results

The results showed that 8,462 genes, encompassing 72,067 nsSNPs were identified, and 118 nsSNPs in 46 genes were predicted as deleterious. GO network analysis classified 13 genes into 5 GO terms (p<0.05) that were associated with kidney development and metabolic processes. Seven genes encompassing nsSNPs were classified into the term associated with Alzheimer’s disease by referencing the genetic association database. The KEGG pathway analysis identified only one significantly enriched pathway (p<0.05), hsa04080: Neuroactive ligand-receptor interaction, among the transcripts.

Conclusion

The number of deleterious nsSNPs in YMPs was identified and then these variants-containing genes in YMPs data were adopted as the putative human diseases-related genes. The results revealed that many genes encompassing nsSNPs in YMPs were related to the various human genes which are potentially associated with kidney development and metabolic processes as well as human disease occurrence.

Transformation of animal genomics by next-generation sequencing technologies: a decade of challenges and their impact on genetic architecture

For more than a quarter of a century, sequencing technologies from Sanger's method to next-generation high-throughput techniques have provided fascinating opportunities in the life sciences. The continuing upward trajectory of sequencing technologies will improve livestock research and expedite the development of various new genomic and technological studies with farm animals. The use of high-throughput technologies in livestock research has increased interest in metagenomics, epigenetics, genome-wide association studies, and identification of single nucleotide polymorphisms and copy number variations. Such studies are beginning to provide revolutionary insights into biological and evolutionary processes. Farm animals, such as cattle, swine, and horses, have played a dual role as economically and agriculturally important animals as well as biomedical research models. The first part of this study explores the current state of sequencing methods, many of which are already used in animal genomic studies, and the second part summarizes the state of cattle, swine, horse, and chicken genome sequencing and illustrates its achievements during the last few years. Finally, we describe several high-throughput sequencing approaches for the improved detection of known, unknown, and emerging infectious agents, leading to better diagnosis of infectious diseases. The insights from viral metagenomics and the advancement of next-generation sequencing will strongly support specific and efficient vaccine development and provide strategies for controlling infectious disease transmission among animal populations and/or between animals and humans. However, prospective sequencing technologies will require further research and in-field testing before reaching the marketplace.

Transcriptomic profile of leg muscle during early growth in chicken

The early growth pattern, especially the age of peak growth, of broilers affects the time to market and slaughter weight, which in turn affect the profitability of the poultry industry. However, the underlying mechanisms regulating chicken growth and development have rarely been studied. This study aimed to identify candidate genes involved in chicken growth and investigated the potential regulatory mechanisms of early growth in chicken. RNA sequencing was applied to compare the transcriptomes of chicken muscle tissues at three developmental stages during early growth. In total, 978 differentially expressed genes (DEGs) (fold change ≥ 2; false discovery rate < 0.05) were detected by pairwise comparison. Functional analysis showed that the DEGs are mainly involved in the processes of cell growth, muscle development, and cellular activities (such as junction, migration, assembly, differentiation, and proliferation). Many of the DEGs are well known to be related to chicken growth, such as MYOD1, GH, IGF2BP2, IGFBP3, SMYD1, CEBPB, FGF2, and IGFBP5. KEGG pathway analysis identified that the DEGs were significantly enriched in five pathways (P < 0.1) related to growth and development: extracellular matrix–receptor interaction, focal adhesion, tight junction, insulin signaling pathway, and regulation of the actin cytoskeleton. A total of 42 DEGs assigned to these pathways are potential candidate genes inducing the difference in growth among the three developmental stages, such as MYH10, FGF2, FGF16, FN1, CFL2, MAPK9, IRS1, PHKA1, PHKB, and PHKG1. Thus, our study identified a series of genes and several pathways that may participate in the regulation of early growth in chicken. These results should serve as an important resource revealing the molecular basis of chicken growth and development.

Single nucleotide polymorphism discovery in bovine liver using RNA-seq technology

BACKGROUND

RNA-seq is a useful next-generation sequencing (NGS) technology that has been widely used to understand mammalian transcriptome architecture and function. In this study, a breed-specific RNA-seq experiment was utilized to detect putative single nucleotide polymorphisms (SNPs) in liver tissue of young bulls of the Polish Red, Polish Holstein-Friesian (HF) and Hereford breeds, and to understand the genomic variation in the three cattle breeds that may reflect differences in production traits.

RESULTS

The RNA-seq experiment on bovine liver produced 107,114,4072 raw paired-end reads, with an average of approximately 60 million paired-end reads per library. Breed-wise, a total of 345.06, 290.04 and 436.03 million paired-end reads were obtained from the Polish Red, Polish HF, and Hereford breeds, respectively. Burrows-Wheeler Aligner (BWA) read alignments showed that 81.35%, 82.81% and 84.21% of the mapped sequencing reads were properly paired to the Polish Red, Polish HF, and Hereford breeds, respectively. This study identified 5,641,401 SNPs and insertion and deletion (indel) positions expressed in the bovine liver with an average of 313,411 SNPs and indel per young bull. Following the removal of the indel mutations, a total of 195,3804, 152,7120 and 205,3184 raw SNPs expressed in bovine liver were identified for the Polish Red, Polish HF, and Hereford breeds, respectively. Breed-wise, three highly reliable breed-specific SNP-databases (SNP-dbs) with 31,562, 24,945 and 28,194 SNP records were constructed for the Polish Red, Polish HF, and Hereford breeds, respectively. Using a combination of stringent parameters of a minimum depth of ≥10 mapping reads that support the polymorphic nucleotide base and 100% SNP ratio, 4,368, 3,780 and 3,800 SNP records were detected in the Polish Red, Polish HF, and Hereford breeds, respectively. The SNP detections using RNA-seq data were successfully validated by kompetitive allele-specific PCR (KASPTM) SNP genotyping assay. The comprehensive QTL/CG analysis of 110 QTL/CG with RNA-seq data identified 20 monomorphic SNP hit loci (CARTPT, GAD1, GDF5, GHRH, GHRL, GRB10, IGFBPL1, IGFL1, LEP, LHX4, MC4R, MSTN, NKAIN1, PLAG1, POU1F1, SDR16C5, SH2B2, TOX, UCP3 and WNT10B) in all three cattle breeds. However, six SNP loci (CCSER1, GHR, KCNIP4, MTSS1, EGFR and NSMCE2) were identified as highly polymorphic among the cattle breeds.

CONCLUSIONS

This study identified breed-specific SNPs with greater SNP ratio and excellent mapping coverage, as well as monomorphic and highly polymorphic putative SNP loci within QTL/CGs of bovine liver tissue. A breed-specific SNP-db constructed for bovine liver yielded nearly six million SNPs. In addition, a KASPTM SNP genotyping assay, as a reliable cost-effective method, successfully validated the breed-specific putative SNPs originating from the RNA-seq experiments.