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

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

Comprehensive Atlas of Alternative Splicing Reveals NSRP1 Promoting Adipogenesis through CCDC18.

Alternative splicing (AS) plays a crucial role in regulating gene expression, function, and diversity. However, limited reports exist on the identification and comparison of AS in Eastern and Western pigs. Here, we analyzed 243 transcriptome data from eight tissues, integrating information on transcription factors (TFs), selection signals, splicing factors (SFs), and quantitative trait loci (QTL) to comprehensively study alternative splicing events (ASEs) in pigs. Five ASE types were identified, with Mutually Exclusive Exon (MXE) and Skipped Exon (SE) ASEs being the most prevalent. A significant portion of genes with ASEs (ASGs) showed conservation across all eight tissues (63.21-76.13% per tissue). Differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) exhibited tissue specificity, with blood and adipose tissues having more DASGs. Functional enrichment analysis revealed coDASG_DEGs in adipose were enriched in pathways associated with adipose deposition and immune inflammation, while coDASG_DEGs in blood were enriched in pathways related to immune inflammation and metabolism. Adipose deposition in Eastern pigs might be linked to the down-regulation of immune-inflammation-related pathways and reduced insulin resistance. The TFs, selection signals, and SFs appeared to regulate ASEs. Notably, ARID4A (TF), NSRP1 (SF), ANKRD12, IFT74, KIAA2026, CCDC18, NEXN, PPIG, and ROCK1 genes in adipose tissue showed potential regulatory effects on adipose-deposition traits. NSRP1 could promote adipogenesis by regulating alternative splicing and expression of CCDC18. Conducting an in-depth investigation into AS, this study has successfully identified key marker genes essential for pig genetic breeding and the enhancement of meat quality, which will play important roles in promoting the diversity of pork quality and meeting market demand.

S-dinotefuran affects the social behavior of honeybees (Apis mellifera)and increases their risk in the colony.

The toxic effects of neonicotinoid pesticides on honeybees is a global concern, whereas little is known about the effect of stereoisomeric pesticides among honeybee social behavior. In this study, we investigated the effects of stereoisomeric dinotefuran on honeybee social behavior. We found that honeybees exhibit a preference for consuming food containing S-dinotefuran, actively engage in trophallaxis with S-dinotefuran-consuming peers, and consequently acquire higher levels of S-dinotefuran compared with R-dinotefuran. In comparison to R-dinotefuran, S-dinotefuran stimulates honeybees to elevate their body temperature, thereby attracting more peers for trophallaxis. Transcriptome analysis revealed a significant enrichment of thermogenesis pathways due to S-dinotefuran exposure. Additionally, metabolome data indicated that S-dinotefuran may enhance body temperature by promoting lipid synthesis in the lysine degradation pathway. Consequently, body temperature emerges as a key factor influencing honeybee social behavior. Our study is the first to highlight the propensity of S-dinotefuran to raise honeybee body temperature, which prompts honeybee to preferentially engage in trophallaxis with peers exhibiting higher body temperatures. This preference may lead honeybees to collect more dinotefuran-contaminated food in the wild, significantly accelerating dinotefuran transmission within a population. Proactive trophallaxis further amplifies the risk of neonicotinoid pesticide transmission within a population, making honeybees that have consumed S-dinotefuran particularly favored within their colonies. These findings may contribute to our understanding of the higher risk associated with neonicotinoid use compared with other pesticides.

Revisiting the Pig IGHC Gene Locus in Different Breeds Uncovers Nine Distinct IGHG Genes.

IgG subclass diversification is common in placental mammals. It has been well documented in humans and mice that different IgG subclasses, with diversified functions, synergistically regulate humoral immunity. However, our knowledge on the genomic and functional diversification of IgG subclasses in the pig, a mammalian species with high agricultural and biomedical importance, is incomplete. Using bacterial artificial chromosome sequencing and newly assembled genomes generated by the PacBio sequencing approach, we characterized and mapped the IgH C region gene locus in three indigenous Chinese breeds (Erhualian, Xiang, and Luchuan) and compared them to that of Duroc. Our data revealed that IGHG genes in Chinese pigs differ from the Duroc, whereas the IGHM, IGHD, IGHA, and IGHE genes were all single copy and highly conserved in the pig breeds examined. Most striking were differences in numbers of IGHG genes: there are seven genes in Erhualian pigs, six in the Duroc, but only five in Xiang pigs. Phylogenetic analysis suggested that all reported porcine IGHG genes could be classified into nine subclasses: IGHG1, IGHG2a, IGHG2b, IGHG2c, IGHG3, IGHG4, IGHG5a, IGHG5b, and IGHG5c. Using sequence information, we developed a mouse mAb specific for IgG3. This study offers a starting point to investigate the structure-function relationship of IgG subclasses in pigs.

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

BACKGROUND: The genetic mechanisms that underlie phenotypic differentiation in breeding animals have important implications in evolutionary biology and agriculture. However, the contribution of cis-regulatory variants to pig phenotypes is poorly understood. Therefore, our aim was to elucidate the molecular mechanisms by which non-coding variants cause phenotypic differences in pigs by combining evolutionary biology analyses and functional genomics. RESULTS: We obtained a high-resolution phased chromosome-scale reference genome with a contig N50 of 18.03 Mb for the Luchuan pig breed (a representative eastern breed) and profiled potential selective sweeps in eastern and western pigs by resequencing the genomes of 234 pigs. Multi-tissue transcriptome and chromatin accessibility analyses of these regions suggest that tissue-specific selection pressure is mediated by promoters and distal cis-regulatory elements. Promoter variants that are associated with increased expression of the lysozyme (LYZ) gene in the small intestine might enhance the immunity of the gastrointestinal tract and roughage tolerance in pigs. In skeletal muscle, an enhancer-modulating single-nucleotide polymorphism that is associated with up-regulation of the expression of the troponin C1, slow skeletal and cardiac type (TNNC1) gene might increase the proportion of slow muscle fibers and affect meat quality. CONCLUSIONS: Our work sheds light on the molecular mechanisms by which non-coding variants shape phenotypic differences in pigs and provides valuable resources and novel perspectives to dissect the role of gene regulatory evolution in animal domestication and breeding.

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

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

Differential day-night expression of tight junction components in murine retinal pigment epithelium.

Strong communication and interaction between the retinal pigment epithelium (RPE) and the photoreceptor (PR) cells is essential for vision. RPE cells are essential for supporting and maintaining PR cells by transporting nutrients, waste products and ions, and phagocytosing photoreceptor outer segments (POS). POS phagocytosis follows a circadian pattern, taking place in the morning in human, mice and other organisms. However, it remains unknown whether other RPE processes follow a daily rhythm. To study the daily rhythm of RPE cells, we isolated murine RPE cells at six different time points during a 24 h period, after which RNA was isolated and sequenced. Murine RPE flatmounts were isolated at four different time points to study daily rhythm in protein abundance and localisation. EnrichR pathway analysis resulted in 13 significantly-enriched KEGG pathways (p < 0.01) of which seven showed a large number of overlapping genes. Several genes were involved in intracellular trafficking, possibly playing a role in nutrient transport, POS phagocytosis or membrane protein trafficking, with different expression patterns during the day-night cycle. Other genes were involved in actin cytoskeleton building, remodelling and crosslinking and showed a high expression in the morning, suggesting actin cytoskeleton remodelling at this time point. Finally, tight junction proteins Cldn2 and Cldn4 showed a difference in RNA and protein expression and tight junction localisation over time. Our study suggests that several important processes in the RPE follow a day-night rhythm, including intracellular trafficking, and processes involving the actin cytoskeleton and tight junctions. The differential protein localisation of Cldn2 in the RPE during the day-night cycle suggest that Cldn2 may facilitate paracellular water and sodium transport during the day.

Genome-wide association study revealed a promising region and candidate genes for eggshell quality in an F2 resource population.

BACKGROUND: Eggshell is subject to quality loss with aging process of laying hens, and damaged eggshells result in economic losses of eggs. However, the genetic architecture underlying the dynamic eggshell quality remains elusive. Here, we measured eggshell quality traits, including eggshell weight (ESW), eggshell thickness (EST) and eggshell strength (ESS) at 11 time points from onset of laying to 72 weeks of age and conducted comprehensive genome-wide association studies (GWAS) in 1534 F2 hens derived from reciprocal crosses between White Leghorn (WL) and Dongxiang chickens (DX). RESULTS: ESWs at all ages exhibited moderate SNP-based heritability estimates (0.30 ~ 0.46), while the estimates for EST (0.21 ~ 0.31) and ESS (0.20 ~ 0.27) were relatively low. Eleven independent univariate genome-wide screens for each trait totally identified 1059, 1026 and 1356 significant associations with ESW, EST and ESS, respectively. Most significant loci were in a region spanning from 57.3 to 71.4 Mb of chromosome 1 (GGA1), which together account for 8.4 ~ 16.5% of the phenotypic variance for ESW from 32 to 72 weeks of age, 4.1 ~ 6.9% and 2.95 ~ 16.1% for EST and ESS from 40 to 72 weeks of age. According to linkage disequilibrium (LD) and conditional analysis, the significant SNPs in this region were in extremely strong linkage disequilibrium status. Ultimately, two missense SNPs in GGA1 and one in GGA4 were considered as promising loci on three independent genes including ITPR2, PIK3C2G, and NCAPG. The homozygotes of advantageously effective alleles on PIK3C2G and ITPR2 possessed the best eggshell quality and could partly counteract the negative effect of aging process. NCAPG had certain effect on eggshell quality for young hens. CONCLUSIONS: Identification of the promising region as well as potential candidate genes will greatly advance our understanding of the genetic basis underlying dynamic eggshell quality and has the practical significance in breeding program for the improvement of eggshell quality, especially at the later part of laying cycle.

Genome-wide association study dissects genetic architecture underlying longitudinal egg weights in chickens.

BACKGROUND: As a major economic trait in chickens, egg weight (EW) receives widespread interests in breeding, production and consumption. However, limited information is available for underlying genetic architecture of longitudinal trend in EW. Herein, we measured EWs at nine time points from onset of laying to 60 week of age, and conducted comprehensive genome-wide association studies (GWAS) in 1,534 F2 hens derived from reciprocal crosses between White Leghorn and Dongxiang chickens. RESULTS: Egg weights at all ages except the first egg weight (FEW) exhibited high SNP-based heritability estimates (0.47~0.60). Strong pair-wise genetic correlations (0.77~1.00) were found among all EWs. Nine separate univariate genome-wide screens suggested 73 signals showing significant associations with longitudinal EWs. After multivariate and conditional analyses, four variants on three chromosomes remained independent contributions. The minor alleles at two loci exerted consistent and positive substitution effects on EWs, and other two were negative. The four loci together accounted for 3.84 % of the phenotypic variance for FEW and 7.29~11.06 % for EWs from 32 to 60 week of age. We obtained five candidate genes, of which NCAPG harbors a non-synonymous SNP (rs14491030) causing a valine-to-alanine amino-acid substitution. Genome partitioning analysis indicated a strong linear correlation between the variance explained by each chromosome and its length, which provided evidence that EW follows a highly polygenic nature of inheritance. CONCLUSIONS: Identification of significant genetic causes that together implicate EWs at different ages will greatly advance our understanding of the genetic basis behind longitudinal EWs, and would be helpful to illuminate the future breeding direction on how to select desired egg size.

Genome-wide association studies for feed intake and efficiency in two laying periods of chickens.

BACKGROUND: Feed contributes to over 60 % of the total production costs in the poultry industry. Increasing feed costs prompt geneticists to include feed intake and efficiency as selection goals in breeding programs. In the present study, we used an F2 chicken population in a genome-wide association study (GWAS) to detect potential genetic variants and candidate genes associated with daily feed intake (FI) and feed efficiency, including residual feed intake (RFI) and feed conversion ratio (FCR). METHODS: A total of 1534 F2 hens from a White Leghorn and Dongxiang reciprocal cross were phenotyped for feed intake and efficiency between 37 and 40 weeks (FI1, RFI1, and FCR1) and between 57 and 60 weeks (FI2, RFI2, and FCR2), and genotyped using the chicken 600 K single nucleotide polymorphism (SNP) genotyping array. Univariate, bivariate, and conditional genome-wide association studies (GWAS) were performed with GEMMA, a genome-wide efficient mixed model association algorithm. The statistical significance threshold for association was inferred by the simpleM method. RESULTS: We identified eight genomic regions that each contained at least one genetic variant that showed a significant association with FI. Genomic regions on Gallus gallus (GGA) chromosome 4 coincided with known quantitative trait loci (QTL) that affect feed intake of layers. Of particular interest, eight SNPs on GGA1 in the region between 169.23 and 171.55 Mb were consistently associated with FI in both univariate and bivariate GWAS, which explained 3.72 and 2.57 % of the phenotypic variance of FI1 and FI2, respectively. The CAB39L gene can be considered as a promising candidate for FI1. For RFI, a haplotype block on GGA27 harbored a significant SNP associated with RFI2. The major allele of rs315135692 was favorable for a lower RFI, with a phenotypic difference of 3.35 g/day between opposite homozygous genotypes. Strong signals on GGA1 were detected in the bivariate GWAS for FCR. CONCLUSIONS: The results demonstrated the polygenic nature of feed intake. GWAS identified novel variants and confirmed a QTL that was previously reported for feed intake in chickens. Genetic variants associated with feed efficiency may be used in genomic breeding programs to select more efficient layers.

Genetic parameters of feed efficiency traits in laying period of chickens.

Laying records on 1,534 F2 hens, derived from a reciprocal cross between White Leghorns and Dongxiang blue-shelled chickens, were used to estimate genetic parameters for residual feed intake (RFI), feed conversion ratio (FCR), daily feed intake (FI), metabolic BW (MBW), BW gain (BWG), and daily egg mass (EM) at 37 to 40 (T1) and 57 to 60 wk age (T2), respectively. Genetic analysis was subsequently conducted with the AI-REML method using an animal model. Estimates for heritability of RFI, FCR, and FI were 0.21, 0.19, and 0.20 in T1, and 0.29, 0.13, and 0.26 in T2, respectively. In T1 and T2, RFI showed high and positive genetic correlations with FCR (0.51, 0.43) and FI (0.72, 0.84), whereas the genetic correlation between FI and FCR was very low (-0.09, 0.11). Genetically, negative correlations were found between RFI and its component traits (-0.01 to -0.47). In addition, high genetic correlations, from 0.76 to 0.94, were observed between T1 and T2 for RFI, FCR, and FI, suggesting that feed efficiency traits in the 2 stages had a similar genetic background. The results indicate that selection for low RFI could reduce FI without significant changes in EM, while selection on FCR will increase EM. The present study lays the foundation for genetic improvement of feed efficiency during the laying period of chickens.

Genome-wide patterns of copy number variation in the diversified chicken genomes using next-generation sequencing.

BACKGROUND: Copy number variation (CNV) is important and widespread in the genome, and is a major cause of disease and phenotypic diversity. Herein, we performed a genome-wide CNV analysis in 12 diversified chicken genomes based on whole genome sequencing. RESULTS: A total of 8,840 CNV regions (CNVRs) covering 98.2 Mb and representing 9.4% of the chicken genome were identified, ranging in size from 1.1 to 268.8 kb with an average of 11.1 kb. Sequencing-based predictions were confirmed at a high validation rate by two independent approaches, including array comparative genomic hybridization (aCGH) and quantitative PCR (qPCR). The Pearson's correlation coefficients between sequencing and aCGH results ranged from 0.435 to 0.755, and qPCR experiments revealed a positive validation rate of 91.71% and a false negative rate of 22.43%. In total, 2,214 (25.0%) predicted CNVRs span 2,216 (36.4%) RefSeq genes associated with specific biological functions. Besides two previously reported copy number variable genes EDN3 and PRLR, we also found some promising genes with potential in phenotypic variation. Two genes, FZD6 and LIMS1, related to disease susceptibility/resistance are covered by CNVRs. The highly duplicated SOCS2 may lead to higher bone mineral density. Entire or partial duplication of some genes like POPDC3 may have great economic importance in poultry breeding. CONCLUSIONS: Our results based on extensive genetic diversity provide a more refined chicken CNV map and genome-wide gene copy number estimates, and warrant future CNV association studies for important traits in chickens.

Genetic analysis for dynamic changes of egg weight in 2 chicken lines.

One of the main concerns for poultry producers is how to maintain egg uniformity and stability in size and weight following the rapid growth during the early laying period. In this study, we aimed to investigate the increase in egg weight with advancing hen age, and to estimate genetic parameters of these increment traits in 2 pure lines of chickens (i.e., 2,010 White Leghorn and 1,200 brown-egg dwarf hens), using the restricted maximum likelihood method with the DMU procedure. We collected age at first egg (AFE), first egg weight (FEW) and kept records of egg weight per 10 wk from 30 to 60 wk of age. Meanwhile, the increments of egg weight were calculated for the evaluation of age-dependent dynamic changes. The increment of egg weight gained dramatically before 30 wk of age and became slower with the advance of age. Heritability estimates of AFE were larger than 0.32, and the low to moderate genetic correlations between AFE and FEW were observed in the 2 lines. The FEW showed high variation level compared with egg weights at later ages in the 2 lines, and had moderate heritability estimates in White Leghorns (0.20) and dwarf hens (0.33). Egg weights at different ages were highly heritable in the 2 lines (h(2) ≥ 0.35), and had strong genetic and phenotypic correlations among different ages. The estimates of heritability for most increment traits were low to moderate, especially those increments for 10-wk intervals ranging from 0.00 to 0.14. The genetic correlations among 3 consecutive egg weight increments for 10-wk intervals were low to moderate. Our results in the 2 lines should provide important insights into the genetic architecture of increment traits and offer some suggestions for producing uniform and stable eggs in response to advancing age.

Quercetin inhibits the epithelial-mesenchymal transition and reverses CDK4/6 inhibitor resistance in breast cancer by regulating circHIAT1/miR-19a-3p/CADM2 axis.

Breast cancer (BC) cells have a high risk of metastasis due to epithelial-mesenchymal transition (EMT). Palbociclib (CDK4/6 inhibitor) is an approved drug for BC treatment. However, the drug resistance and metastasis can impair the treatment outcome of Palbociclib. Understanding the mechanisms of EMT and Palbociclib drug resistance in BC is conducive to the formulation of novel therapeutic strategy. Here, we investigated the role of circHIAT1/miR-19a-3p/CADM2 axis in modulating EMT and Palbociclib resistance in BC. circHIAT1 and CADM2 were down-regulated in BC tissues and cell lines, and miR-19a-3p showed an up-regulation. circHIAT1 could interact with miR-19a-3p and suppress its activity, while miR-19a-3p functioned to negatively regulate CADM2. Forced over-expression of circHIAT1 could impaired the EMT status and migratory ability of BC cells, and this effect was inhibited by miR-19a-3p mimic. In addition, we also generated Palbociclib resistant BC cells, and showed that circHIAT1 and CADM2 were down-regulated in the resistant BC cells while miR-19a-3p showed an up-regulation. Forced circHIAT1 over-expression re-sensitized BC cells to Palbociclib treatment. Quercetin, a bioactive flavonoid, could suppressed the migration and invasion of BC cells, and re-sensitized BC cells to Palbociclib. The anti-cancer effect of quercetin could be attributed to its regulatory effect on circHIAT1/miR-19a-3p/CADM2 axis. In vivo tumorigenesis experiment further revealed that quercetin administration enhanced the anti-cancer effect of Palbociclib, an effect was dependent on the up-regulation of circHIAT1 by quercetin. In summary, this study identified quercetin as a potential anti-cancer compound to reverse Palbociclib resistance and impair EMT in BC cells by targeting circHIAT1/miR-19a-3p/CADM2 axis.

Whole-genome analysis reveals distinct adaptation signatures to diverse environments in Chinese domestic pigs.

BACKGROUND: Long-term natural and artificial selection has resulted in many genetic footprints within the genomes of pig breeds across distinct agroecological zones. Nevertheless, the mechanisms by which these signatures contribute to phenotypic diversity and facilitate environmental adaptation remain unclear. RESULTS: Here, we leveraged whole-genome sequencing data from 82 individuals from 6 domestic pig breeds originating in tropical, high-altitude, and frigid regions. Population genetic analysis suggested that habitat isolation significantly shaped the genetic diversity and contributed to population stratification in local Chinese pig breeds. Analysis of selection signals revealed regions under selection for adaptation in tropical (55.5 Mb), high-altitude (43.6 Mb), and frigid (17.72 Mb) regions. The potential functions of the selective sweep regions were linked to certain complex traits that might play critical roles in different geographic environments, including fat coverage in frigid environments and blood indicators in tropical and high-altitude environments. Candidate genes under selection were significantly enriched in biological pathways involved in environmental adaptation. These pathways included blood circulation, protein degradation, and inflammation for adaptation to tropical environments; heart and lung development, hypoxia response, and DNA damage repair for high-altitude adaptation; and thermogenesis, cold-induced vasodilation (CIVD), and the cell cycle for adaptation to frigid environments. By examining the chromatin state of the selection signatures, we identified the lung and ileum as two candidate functional tissues for environmental adaptation. Finally, we identified a mutation (chr1: G246,175,129A) in the cis-regulatory region of ABCA1 as a plausible promising variant for adaptation to tropical environments. CONCLUSIONS: In this study, we conducted a genome-wide exploration of the genetic mechanisms underlying the adaptability of local Chinese pig breeds to tropical, high-altitude, and frigid environments. Our findings shed light on the prominent role of cis-regulatory elements in environmental adaptation in pigs and may serve as a valuable biological model of human plateau-related disorders and cardiovascular diseases.

An intronic SNP affects skeletal muscle development by regulating the expression of TP63.

Background: Porcine skeletal muscle development is pivotal for improving meat production. TP63, a transcription factor, regulates vital cellular processes, yet its role in skeletal muscle proliferation is unclear. Methods: The effects of TP63 on skeletal muscle cell viability and proliferation were investigated using both mouse and porcine skeletal muscle myoblasts. Selective sweep analysis in Western pigs identified TP63 as a potential candidate gene for skeletal muscle development. The correlation between TP63 overexpression and cell proliferation was assessed using quantitative real-time PCR (RT-qPCR) and 5-ethynyl-2'-deoxyuridine (EDU). Results: The study revealed a positive correlation between TP63 overexpression and skeletal muscle cell proliferation. Bioinformatics analysis predicted an interaction between MEF2A, another transcription factor, and the mutation site of TP63. Experimental validation through dual-luciferase assays confirmed that a candidate enhancer SNP could influence MEF2A binding, subsequently regulating TP63 expression and promoting skeletal muscle cell proliferation. Conclusion: These findings offer experimental evidence for further exploration of skeletal muscle development mechanisms and the advancement of genetic breeding strategies aimed at improving meat production traits.

Single-nucleus transcriptional and chromatin accessible profiles reveal critical cell types and molecular architecture underlying chicken sex determination.

INTRODUCTION: Understanding the sex determination mechanisms in birds has great significance for the biological sciences and production in the poultry industry. Sex determination in chickens is a complex process that involves fate decisions of supporting cells such as granulosa or Sertoli cells. However, a systematic understanding of the genetic regulation and cell commitment process underlying sex determination in chickens is still lacking. OBJECTIVES: We aimed to dissect the molecular characteristics associated with sex determination in the gonads of chicken embryos. METHODS: Single-nucleus RNA-seq (snRNA-seq) and ATAC-seq (snATAC-seq) analysis were conducted on the gonads of female and male chickens at embryonic day 3.5 (E3.5), E4.5, and E5.5. RESULTS: Here, we provided a time-course transcriptional and chromatin accessible profiling of gonads during chicken sex determination at single-cell resolution. We uncovered differences in cell composition and developmental trajectories between female and male gonads and found that the divergence of transcription and accessibility in gonadal cells first emerged at E5.5. Furthermore, we revealed key cell-type-specific transcription factors (TFs) and regulatory networks that drive lineage commitment. Sex determination signaling pathways, dominated by BMP signaling, are preferentially activated in males during gonadal development. Further pseudotime analysis of the supporting cells indicated that granulosa cells were regulated mainly by the TEAD gene family and that Sertoli cells were driven by the DMRT1 regulons. Cross-species analysis suggested high conservation of both cell types and cell-lineage-specific TFs across the six vertebrates. CONCLUSIONS: Overall, our study will contribute to accelerating the development of sex manipulation technology in the poultry industry and the application of chickens as a unique model for studying cell fate decisions.

Uncovering Evolutionary Adaptations in Common Warthogs through Genomic Analyses.

In the Suidae family, warthogs show significant survival adaptability and trait specificity. This study offers a comparative genomic analysis between the warthog and other Suidae species, including the Luchuan pig, Duroc pig, and Red River hog. By integrating the four genomes with sequences from the other four species, we identified 8868 single-copy orthologous genes. Based on 8868 orthologous protein sequences, phylogenetic assessments highlighted divergence timelines and unique evolutionary branches within suid species. Warthogs exist on different evolutionary branches compared to DRCs and LCs, with a divergence time preceding that of DRC and LC. Contraction and expansion analyses of warthog gene families have been conducted to elucidate the mechanisms of their evolutionary adaptations. Using GO, KEGG, and MGI databases, warthogs showed a preference for expansion in sensory genes and contraction in metabolic genes, underscoring phenotypic diversity and adaptive evolution direction. Associating genes with the QTLdb-pigSS11 database revealed links between gene families and immunity traits. The overlap of olfactory genes in immune-related QTL regions highlighted their importance in evolutionary adaptations. This work highlights the unique evolutionary strategies and adaptive mechanisms of warthogs, guiding future research into the distinct adaptability and disease resistance in pigs, particularly focusing on traits such as resistance to African Swine Fever Virus.

The sodium channel gene family is specifically expressed in hen uterus and associated with eggshell quality traits.

BACKGROUND: Eggshell quality is important for the poultry industry. During eggshell formation a mass of inorganic minerals is deposited. The Sodium Channel (SCNN1) gene family plays an essential role in cation transportation. The objective of this study was to investigate the pattern of expression of members of the SCNN1 gene family, their variation and their effects on eggshell quality. RESULT: The highest expression of SCNN1a, SCNN1b, and SCNN1g genes were in the active uterus during eggshell mineralization, while SCNN1d showed its highest expression level in the quiescent uterus (no egg present). Nineteen candidate SNPs from the four genes were genotyped in a population of 338 White Leghorn layers. Association analysis between SNPs (haplotypes/diplotypes) and eggshell traits was performed. Among seven significant SNPs, five SNPs were associated with eggshell strength, eggshell thickness, eggshell percentage or/and egg weight, while the other two SNPs within SCNN1d were only associated with eggshell percentage. These SNPs had a 0.25-6.99% contribution to phenotypic variance, depending on the trait. In haplotype analysis, SCNN1b and SCNN1d were associated with egg weight. The SCNN1b and SCNN1g were significantly associated with eggshell weight while only SCNN1g explained 2.04% of phenotypic variance. All the alleles of the members of SCNN1 gene family were associated with eggshell percentage and eggshell thickness, and others members had an association with eggshell strength except for SCNN1a. The contribution of different haplotypes of the SCNN1 gene family to eggshell phenotypic variance ranged from 0.09% to 5.74%. CONCLUSIONS: Our study indicated that the SCNN1 gene family showed tissue expression specificity and was significantly associated with eggshell traits in chicken. This study provides evidence that genetic variation in members of the sodium channel can influence eggshell quality.