Haplotype-resolved 3D chromatin architecture of the hybrid pig.

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

Spatio-temporal transcriptome dynamics coordinate rapid transition of core crop functions in 'lactating' pigeon.

Pigeons (Columba livia) are among a select few avian species that have developed a specialized reproductive mode wherein the parents produce a 'milk' in their crop to feed newborn squabs. Nonetheless, the transcriptomic dynamics and role in the rapid transition of core crop functions during 'lactation' remain largely unexplored. Here, we generated a de novo pigeon genome assembly to construct a high resolution spatio-temporal transcriptomic landscape of the crop epithelium across the entire breeding stage. This multi-omics analysis identified a set of 'lactation'-related genes involved in lipid and protein metabolism, which contribute to the rapid functional transitions in the crop. Analysis of in situ high-throughput chromatin conformation capture (Hi-C) sequencing revealed extensive reorganization of promoter-enhancer interactions linked to the dynamic expression of these 'lactation'-related genes between stages. Moreover, their expression is spatially localized in specific epithelial layers, and can be correlated with phenotypic changes in the crop. These results illustrate the preferential de novo synthesis of 'milk' lipids and proteins in the crop, and provides candidate enhancer loci for further investigation of the regulatory elements controlling pigeon 'lactation'.

Comparative three-dimensional genome architectures of adipose tissues provide insight into human-specific regulation of metabolic homeostasis.

Elucidating the regulatory mechanisms of human adipose tissues (ATs) evolution is essential for understanding human-specific metabolic regulation, but the functional importance and evolutionary dynamics of three-dimensional (3D) genome organizations of ATs are not well defined. Here, we compared the 3D genome architectures of anatomically distinct ATs from humans and six representative mammalian models. We recognized evolutionarily conserved and human-specific chromatin conformation in ATs at multiple scales, including compartmentalization, topologically associating domain (TAD), and promoter-enhancer interactions (PEI), which have not been described previously. We found PEI are much more evolutionarily dynamic with respect to compartmentalization and topologically associating domain. Compared to conserved PEIs, human-specific PEIs are enriched for human-specific sequence, and the binding motifs of their potential mediators (transcription factors) are less conserved. Our data also demonstrated that genes involved in the evolutionary dynamics of chromatin organization have weaker transcriptional conservation than those associated with conserved chromatin organization. Furthermore, the genes involved in energy metabolism and the maintenance of metabolic homeostasis are enriched in human-specific chromatin organization, while housekeeping genes, health-related genes, and genetic variations are enriched in evolutionarily conserved compared to human-specific chromatin organization. Finally, we showed extensively divergent human-specific 3D genome organizations among one subcutaneous and three visceral ATs. Together, these findings provide a global overview of 3D genome architecture dynamics between ATs from human and mammalian models and new insights into understanding the regulatory evolution of human ATs.

A mutator-derived prognostic eRNA signature provides insight into the pathogenesis of breast cancer.

Abundant evidence suggests that enhancer RNA (eRNA) is closely related to tumorigenesis, and the role of eRNA transcription in promoting genomic instability in cancers is gradually unveiled. However, research on the evaluation of the prognostic value and molecular mechanisms of genomic instability associated eRNAs in breast cancer is long overdue. Here, we integratively analyzed eRNA expression and somatic mutation profiles in breast cancer genome. We identified genomic instability associated eRNAs and developed a prognostic signature based on these eRNAs with the area under the curve (AUC) around 0.8 at 9-year survival. We further found the prognostic value of this signature is independent of common clinical factors and is better than TP53 status. Higher expression of genomic instability associated genes in the high-risk group was observed, suggesting that this eRNA signature may serve as an indicator of genomic instability in breast cancer. We found prognostic eRNA co-expressed genes are mainly enriched in Gene set 'Breast Cancer 8P12-P11 Amplicon', Gene set 'Metabolism of lipids' and GO process 'Ubiquitin protein ligase binding'. Furthermore, 11 eRNA-signature regulated genes are identified by assessing promoter-enhancer interaction. Among these genes, F11R, BHLHE40, and NECTIN4 are previously reported oncogenes and EGOT is a tumor suppressor gene, indicating the direct roles of eRNAs in tumorigenesis.

Exploring high-resolution chromatin interaction changes and functional enhancers of myogenic marker genes during myogenic differentiation.

Skeletal muscle differentiation (myogenesis) is a complex and highly coordinated biological process regulated by a series of myogenic marker genes. Chromatin interactions between gene's promoters and their enhancers have an important role in transcriptional control. However, the high-resolution chromatin interactions of myogenic genes and their functional enhancers during myogenesis remain largely unclear. Here, we used circularized chromosome conformation capture coupled with next generation sequencing (4C-seq) to investigate eight myogenic marker genes in C2C12 myoblasts (C2C12-MBs) and C2C12 myotubes (C2C12-MTs). We revealed dynamic chromatin interactions of these marker genes during differentiation and identified 163 and 314 significant interaction sites (SISs) in C2C12-MBs and C2C12-MTs, respectively. The interacting genes of SISs in C2C12-MTs were mainly involved in muscle development, and histone modifications of the SISs changed during differentiation. Through functional genomic screening, we also identified 25 and 41 putative active enhancers in C2C12-MBs and C2C12-MTs, respectively. Using luciferase reporter assays for putative enhancers of Myog and Myh3, we identified eight activating enhancers. Furthermore, dCas9-KRAB epigenome editing and RNA-Seq revealed a role for Myog enhancers in the regulation of Myog expression and myogenic differentiation in the native genomic context. Taken together, this study lays the groundwork for understanding 3D chromatin interaction changes of myogenic genes during myogenesis and provides insights that contribute to our understanding of the role of enhancers in regulating myogenesis.

Milletdb: a multi-omics database to accelerate the research of functional genomics and molecular breeding of millets.

Millets are a class of nutrient-rich coarse cereals with high resistance to abiotic stress; thus, they guarantee food security for people living in areas with extreme climatic conditions and provide stress-related genetic resources for other crops. However, no platform is available to provide a comprehensive and systematic multi-omics analysis for millets, which seriously hinders the mining of stress-related genes and the molecular breeding of millets. Here, a free, web-accessible, user-friendly millets multi-omics database platform (Milletdb, http://milletdb.novogene.com) has been developed. The Milletdb contains six millets and their one related species genomes, graph-based pan-genomics of pearl millet, and stress-related multi-omics data, which enable Milletdb to be the most complete millets multi-omics database available. We stored GWAS (genome-wide association study) results of 20 yield-related trait data obtained under three environmental conditions [field (no stress), early drought and late drought] for 2 years in the database, allowing users to identify stress-related genes that support yield improvement. Milletdb can simplify the functional genomics analysis of millets by providing users with 20 different tools (e.g., 'Gene mapping', 'Co-expression', 'KEGG/GO Enrichment' analysis, etc.). On the Milletdb platform, a gene PMA1G03779.1 was identified through 'GWAS', which has the potential to modulate yield and respond to different environmental stresses. Using the tools provided by Milletdb, we found that the stress-related PLATZs TFs (transcription factors) family expands in 87.5% of millet accessions and contributes to vegetative growth and abiotic stress responses. Milletdb can effectively serve researchers in the mining of key genes, genome editing and molecular breeding of millets.

Comparative 3D genome architecture in vertebrates.

BACKGROUND: The three-dimensional (3D) architecture of the genome has a highly ordered and hierarchical nature, which influences the regulation of essential nuclear processes at the basis of gene expression, such as gene transcription. While the hierarchical organization of heterochromatin and euchromatin can underlie differences in gene expression that determine evolutionary differences among species, the way 3D genome architecture is affected by evolutionary forces within major lineages remains unclear. Here, we report a comprehensive comparison of 3D genomes, using high resolution Hi-C data in fibroblast cells of fish, chickens, and 10 mammalian species. RESULTS: This analysis shows a correlation between genome size and chromosome length that affects chromosome territory (CT) organization in the upper hierarchy of genome architecture, whereas lower hierarchical features, including local transcriptional availability of DNA, are selected through the evolution of vertebrates. Furthermore, conservation of topologically associating domains (TADs) appears strongly associated with the modularity of expression profiles across species. Additionally, LINE and SINE transposable elements likely contribute to heterochromatin and euchromatin organization, respectively, during the evolution of genome architecture. CONCLUSIONS: Our analysis uncovers organizational features that appear to determine the conservation and transcriptional regulation of functional genes across species. These findings can guide ongoing investigations of genome evolution by extending our understanding of the mechanisms shaping genome architecture.

Revisiting Assessment of Computational Methods for Hi-C Data Analysis.

The performances of algorithms for Hi-C data preprocessing, the identification of topologically associating domains, and the detection of chromatin interactions and promoter-enhancer interactions have been mostly evaluated using semi-quantitative or synthetic data approaches, without utilizing the most recent methods, since 2017. In this study, we comprehensively evaluated 24 popular state-of-the-art methods for the complete end-to-end pipeline of Hi-C data analysis, using manually curated or experimentally validated benchmark datasets, including a CRISPR dataset for promoter-enhancer interaction validation. Our results indicate that, although no single method exhibited superior performance in all situations, HiC-Pro, DomainCaller, and Fit-Hi-C2 showed relatively balanced performances of most evaluation metrics for preprocessing, topologically associating domain identification, and chromatin interaction/promoter-enhancer interaction detection, respectively. The comprehensive comparison presented in this manuscript provides a reference for researchers to choose Hi-C analysis tools that best suit their needs.

Identification of enhancers responsible for the coordinated expression of myosin heavy chain isoforms in skeletal muscle.

BACKGROUND: Skeletal muscles consist of fibers of differing contractility and metabolic properties, which are primarily determined by the content of myosin heavy chain (MYH) isoforms (MYH7, MYH2, MYH1, and MYH4). The regulation of Myh genes transcription depends on three-dimensional chromatin conformation interaction, but the mechanistic details remain to be determined. RESULTS: In this study, we characterized the interaction profiles of Myh genes using 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing). The interaction profile of Myh genes changed between fast quadriceps and slow soleus muscles. Combining chromatin immunoprecipitation-sequencing (ChIP-seq) and transposase accessible chromatin with high-throughput sequencing (ATAC-seq), we found that a 38 kb intergenic region interacting simultaneously with fast Myh genes promoters controlled the coordinated expression of fast Myh genes. We also identified four active enhancers of Myh7, and revealed that binding of MYOG and MYOD increased the activity of Myh7 enhancers. CONCLUSIONS: This study provides new insight into the chromatin interactions that regulate Myh genes expression.

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

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

Temporal microRNA expression profile of pig peripheral blood during postnatal development.

Gene expression profiles of blood can reflect the physiopathologic status of the immune system. The dynamic microRNA (miRNA) expression profiles of peripheral blood from pigs at different developmental stages, and how differential expression of miRNAs might relate to immune system development, are unknown. In this study, peripheral blood samples taken at five developmental stages were used to construct 15 miRNA libraries (three biological replicates/stage): 0 days (newborn), 30 days (weaning), 60 days (weaned), and 180 and 360 days (puberty). We identified 295 known mature miRNAs. Hierarchical clustering of the miRNA expression profile showed significant differences between individuals at the neonatal and postnatal stages. Functional enrichment analysis revealed that miRNAs differentially expressed between pairwise comparisons of the developmental stages were over-represented in immune-related pathways such as toll-like receptor signaling. The time-course of expression of the over-representated miRNAs exhibited a pattern of steady decline over time, for both the complete miRNA compendium and immune-related miRNAs. We identified six marker miRNAs that were highly negatively correlated with chronologic age and enriched for genes involved in immune-related pathways. This study of a peripheral blood miRNA transcriptome offers insight into immune system development in swine and provides a resource for pig genome annotation.

Current Advances on the Important Roles of Enhancer RNAs in Molecular Pathways of Cancer.

Enhancers are critical genomic elements that can cooperate with promoters to regulate gene transcription in both normal and cancer cells. Recent studies reveal that enhancer regions are transcribed to produce a class of noncoding RNAs referred to as enhancer RNAs (eRNAs). Emerging evidence shows that eRNAs play important roles in enhancer activation and enhancer-driven gene regulation, and the expression of eRNAs may be a critical factor in tumorigenesis. The important roles of eRNAs in cancer signaling pathways are also gradually unveiled, providing a new insight into cancer therapy. Here, we review the roles of eRNAs in regulating cancer signaling pathways and discuss the potential of eRNA-targeted therapy for human cancers.

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

Uncovering genetic variation through resequencing is limited by the fact that only sequences with similarity to the reference genome are examined. Reference genomes are often incomplete and cannot represent the full range of genetic diversity as a result of geographical divergence and independent demographic events. To more comprehensively characterize genetic variation of pigs (Sus scrofa), we generated de novo assemblies of nine geographically and phenotypically representative pigs from Eurasia. By comparing them to the reference pig assembly, we uncovered a substantial number of novel SNPs and structural variants, as well as 137.02-Mb sequences harboring 1737 protein-coding genes that were absent in the reference assembly, revealing variants left by selection. Our results illustrate the power of whole-genome de novo sequencing relative to resequencing and provide valuable genetic resources that enable effective use of pigs in both agricultural production and biomedical research.

MicroRNA-351-5p mediates skeletal myogenesis by directly targeting lactamase-ß and is regulated by lnc-mg.

Skeletal muscle is an important and complex organ with a variety of functions in humans and animals. Skeletal myogenesis is a multistep and complex process, and increasing evidence suggests that microRNAs (miRNAs) and long noncoding RNAs (lncRNAs) play critical roles in skeletal myogenesis. In this study the expression of miR-351-5p is dynamically regulated during skeletal myogenesis in vitro and in vivo. Cell-counting kit-8, qRT-PCR, and EdU immunofluorescence analysis showed that miR-351-5p overexpression promoted the proliferation and inhibited the differentiation of C2C12 myoblast, whereas inhibition of miR-351-5p had the opposite effect. In addition, miR-351-5p mediated the regulation of muscle fiber type transition in vivo. In vitro, loss of miR-351-5p in muscle tissues promoted muscle hypertrophy and increased slow-twitch fibers in the gastrocnemius muscles of mice. Luciferase reporter assay and functional analyses demonstrated that lactamase ß ( LACTB) is a direct target of miR-351-5p involved in the regulation of skeletal myogenesis. Expression levels of a myogenesis-associated lncRNA ( lnc-mg) correlated negatively with miR-351-5p and positively with LACTB during C2C12 myoblast proliferation and differentiation. Further analyses showed that lnc-mg acted as a molecular sponge for miR-351-5p, demonstrating its involvement in the negative regulation of LACTB by miR-351-5p during skeletal myogenesis. These findings indicate that miRNA-351-5p functions in skeletal myogenesis by targeting LACTB and is regulated by lnc-mg, supporting the role of the competing endogenous RNA network in skeletal myogenesis.-Du, J., Zhang, P., Zhao, X., He, J., Xu, Y., Zou, Q., Luo, J., Shen, L., Gu, H., Tang, Q., Li, M., Jiang, Y., Tang, G., Bai, L., Li, X., Wang, J., Zhang, S., Zhu, L. MicroRNA-351-5p mediates skeletal myogenesis by directly targeting lactamase ß and is regulated by lnc-mg.

XBP1 promotes triple-negative breast cancer by controlling the HIF1α pathway.

Cancer cells induce a set of adaptive response pathways to survive in the face of stressors due to inadequate vascularization. One such adaptive pathway is the unfolded protein (UPR) or endoplasmic reticulum (ER) stress response mediated in part by the ER-localized transmembrane sensor IRE1 (ref. 2) and its substrate XBP1 (ref. 3). Previous studies report UPR activation in various human tumours, but the role of XBP1 in cancer progression in mammary epithelial cells is largely unknown. Triple-negative breast cancer (TNBC)--a form of breast cancer in which tumour cells do not express the genes for oestrogen receptor, progesterone receptor and HER2 (also called ERBB2 or NEU)--is a highly aggressive malignancy with limited treatment options. Here we report that XBP1 is activated in TNBC and has a pivotal role in the tumorigenicity and progression of this human breast cancer subtype. In breast cancer cell line models, depletion of XBP1 inhibited tumour growth and tumour relapse and reduced the CD44(high)CD24(low) population. Hypoxia-inducing factor 1α (HIF1α) is known to be hyperactivated in TNBCs. Genome-wide mapping of the XBP1 transcriptional regulatory network revealed that XBP1 drives TNBC tumorigenicity by assembling a transcriptional complex with HIF1α that regulates the expression of HIF1α targets via the recruitment of RNA polymerase II. Analysis of independent cohorts of patients with TNBC revealed a specific XBP1 gene expression signature that was highly correlated with HIF1α and hypoxia-driven signatures and that strongly associated with poor prognosis. Our findings reveal a key function for the XBP1 branch of the UPR in TNBC and indicate that targeting this pathway may offer alternative treatment strategies for this aggressive subtype of breast cancer.

Transcriptomic analysis between Normal and high-intake feeding geese provides insight into adipose deposition and susceptibility to fatty liver in migratory birds.

BACKGROUND: Dysregulation of adipogenesis causes metabolic diseases, like obesity and fatty liver. Migratory birds such as geese have a high tolerance of massive energy intake and exhibit little pathological development. Domesticated goose breeds, derivatives of the wild greyleg goose (Anser anser) or swan goose (Anser cygnoides), have high tolerance of energy intake resembling their ancestor species. Thus, goose is potentially a model species to study mechanisms associated with adipogenesis. RESULTS: Phenotypically, goose liver exhibited higher fat accumulation than adipose tissues during fattening (liver increased by 3.35 fold than 1.65 fold in adipose), showing a priority of fat accumulation in liver. We found the number of differentially expressed genes in liver (13.97%) was nearly twice the number of that in adipose (6.60%). These differentially expressed genes in liver function in several important lipid metabolism pathways, immune response, regulation of cancer, while in adipose, terms closely related to protein binding, gluconeogenesis were enriched. Typically, genes like MDH2 and SCD, which have key roles in glycolysis and fatty acids metabolism, had higher fold change in liver than in adipose tissues. Three hundred two differentially expressed long noncoding RNAs involved in regulation of metabolism in liver were also identified. For example, lncRNA XLOC_292762, which was 5.7 kb downstream of FERMT2, a gene involved phosphatidylinositol-3,4,5-trisphosphate binding, was significantly down-regulated after the high-intake feeding period. Further investigation of documented obesity-related orthologous genes in goose suggested that understanding the evolutionary split from mammals in adipogenesis will make goose fatty liver a better resource for future research. CONCLUSIONS: Our research reveals that goose uses liver as the major tissue to regulate a distinct lipid synthesis and degradation flux and the dynamic expression network analyses showed numerous layers of positive responses to both massive energy intake and possible pathological development. Our results offer insights into goose adipogenesis and provide a new perspective for research in human metabolic dysregulation.

Population genomics identifies patterns of genetic diversity and selection in chicken.

BACKGROUND: There are hundreds of phenotypically distinguishable domestic chicken breeds or lines with highly specialized traits worldwide, which provide a unique opportunity to illustrate how selection shapes patterns of genetic variation. There are many local chicken breeds in China. RESULTS: Here, we provide a population genome landscape of genetic variations in 86 domestic chickens representing 10 phenotypically diverse breeds. Genome-wide analysis indicated that sex chromosomes have less genetic diversity and are under stronger selection than autosomes during domestication and local adaptation. We found an evidence of admixture between Tibetan chickens and other domestic population. We further identified strong signatures of selection affecting genomic regions that harbor genes underlying economic traits (typically related to feathers, skin color, growth, reproduction and aggressiveness) and local adaptation (to high altitude). By comparing the genomes of the Tibetan and lowland fowls, we identified genes associated with high-altitude adaptation in Tibetan chickens were mainly involved in energy metabolism, body size maintenance and available food sources. CONCLUSIONS: The work provides crucial insights into the distinct evolutionary scenarios occurring under artificial selection for agricultural production and under natural selection for success at high altitudes in chicken. Several genes were identified as candidates for chicken economic traits and other phenotypic traits.

miR-27a-5p Attenuates Hypoxia-induced Rat Cardiomyocyte Injury by Inhibiting Atg7.

Acute myocardial infarction (AMI) is an ischemic heart disease with high mortality worldwide. AMI triggers a hypoxic microenvironment and induces extensive myocardial injury, including autophagy and apoptosis. MiRNAs, which are a class of posttranscriptional regulators, have been shown to be involved in the development of ischemic heart diseases. We have previously reported that hypoxia significantly alters the miRNA transcriptome in rat cardiomyoblast cells (H9c2), including miR-27a-5p. In the present study, we further investigated the potential function of miR-27a-5p in the cardiomyocyte response to hypoxia, and showed that miR-27a-5p expression was downregulated in the H9c2 cells at different hypoxia-exposed timepoints and the myocardium of a rat AMI model. Follow-up experiments revealed that miR-27a-5p attenuated hypoxia-induced cardiomyocyte injury by regulating autophagy and apoptosis via Atg7, which partly elucidated the anti-hypoxic injury effects of miR-27a-5p. Taken together, this study shows that miR-27a-5p has a cardioprotective effect on hypoxia-induced H9c2 cell injury, suggesting it may be a novel target for the treatment of hypoxia-related heart diseases.

Comprehensive Analysis of lncRNAs and circRNAs Reveals the Metabolic Specialization in Oxidative and Glycolytic Skeletal Muscles.

The biochemical and functional differences between oxidative and glycolytic muscles could affect human muscle health and animal meat quality. However, present understanding of the epigenetic regulation with respect to lncRNAs and circRNAs is rudimentary. Here, porcine oxidative and glycolytic skeletal muscles, which were at the growth curve inflection point, were sampled to survey variant global expression of lncRNAs and circRNAs using RNA-seq. A total of 4046 lncRNAs were identified, including 911 differentially expressed lncRNAs (p < 0.05). The cis-regulatory analysis identified target genes that were enriched for specific GO terms and pathways (p < 0.05), including the oxidation-reduction process, glycolytic process, and fatty acid metabolic. All these were closely related to different phenotypes between oxidative and glycolytic muscles. Additionally, 810 circRNAs were identified, of which 137 were differentially expressed (p < 0.05). Interestingly, some circRNA-miRNA-mRNA networks were found, which were closely linked to muscle fiber-type switching and mitochondria biogenesis in muscles. Furthermore, 44.69%, 39.19%, and 54.01% of differentially expressed mRNAs, lncRNAs, and circRNAs respectively were significantly enriched in pig quantitative trait loci (QTL) regions for growth and meat quality traits. This study reveals a mass of candidate lncRNAs and circRNAs involved in muscle physiological functions, which may improve understanding of muscle metabolism and development from an epigenetic perspective.

Identifying SNPs and candidate genes for three litter traits using single-step GWAS across six parities in Landrace and Large White pigs.

Total number born (TNB), number born alive (NBA), and litter weight born alive (LWB) are critically important traits in pig production. The sow's parity is one of the major factors influencing litter traits. Because of monogenic or polygenic contributions and the presence of temporal gene effects in different sows' parities, it is difficult to clarify the biological and genetic background. To systematically explore the genetic mechanism of litter traits, we conducted 18 GWASs using single-step GWAS (ssGWAS) based on two breeds (908 Landrace and 1,130 Large White sow litter records) for each litter trait in different parities. A total of 300 Landrace and 300 Large White sows were genotyped by sequencing (GBS). ssGWAS was performed separately for each breed and each parity due to population stratification and temporal gene effect. In summary, we identified 80 (15 for Landrace and 65 for Large White), 227 (52 for Landrace, 175 for Large White), and 187 (34 for Landrace, 153 for Large White) single nucleotide polymorphisms (SNPs) affecting TNB, NBA, and LWB, respectively. Of them, we suggest that a total of 22 loci (SSC1: 125098202, SSC1: 117560058, SSC14: 147794697, SSC8: 84823302, SSC9: 143554876, and SSC9: 138766097 for Landrace; SSC1: 4023577, SSC1: 3859573, SSC1: 4891063, SSC16: 5197665, SSC10: 32050819, SSC13: 13552924, SSC13: 92819, SSC17: 3579607, SSC13: 196698221, SSC7: 30918403, SSC16: 46221484, SSC16: 46169204, SSC2: 41988642, SSC2: 44475457, SSC2: 42521875, and SSC7: 58411951 for Large White) are shared by TNB, NBA, and LWB. These results indicate the existence of gene temporal effect in each parity. Furthermore, our findings suggest four interesting candidate genes (FBXL7, ALDH1A2, LEPR, and DDX1) associated with litter traits in different parities that have a major effect on embryonic development progression. In conclusion, 22 crucial SNPs and four interesting candidate genes were identified for three litter traits across six parities. These findings advance our understanding of the genetic architecture of litter traits and confirm the presence of temporal gene effects in different parities. Importantly, functional validation studies for findings of particular interest are recommended in litter traits.