Comparative analyses by sequencing of transcriptomes during skeletal muscle development between pig breeds differing in muscle growth rate and fatness.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:21637832", "kind": "science-record", "revision": 8, "content_hash": "sha256:9d824bfb1dbeba11ea34be928c7030fa8f00c423fd64a27dc4f3839e5776eb7c", "created_at": "2026-11-01T15:40:01Z", "updated_at": "2026-11-01T22:41:54Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Comparative analyses by sequencing of transcriptomes during skeletal muscle development between pig breeds differing in muscle growth rate and fatness.", "doi": "10.1371/journal.pone.0019774", "pmid": "21637832", "authors": ["Xiao Zhao", "Delin Mo", "Anning Li", "Wen Gong", "Shuqi Xiao", "Yue Zhang", "Limei Qin", "Yuna Niu", "Yunxue Guo", "Xiaohong Liu", "Peiqing Cong", "Zuyong He", "Chong Wang", "Jiaqi Li", "Yaosheng Chen"], "journal": "PloS one", "year": 2011, "abstract": "Understanding the dynamics of muscle transcriptome during development and between breeds differing in muscle growth is necessary to uncover the complex mechanism underlying muscle development. Herein, we present the first transcriptome-wide longissimus dorsi muscle development research concerning Lantang (LT, obese) and Landrace (LR, lean) pig breeds during 10 time-points from 35 days-post-coitus (dpc) to 180 days-post-natum (dpn) using Solexa/Illumina's Genome Analyzer. The data demonstrated that myogenesis was almost completed before 77 dpc, but the muscle phenotypes were still changed from 77 dpc to 28 dpn. Comparative analysis of the two breeds suggested that myogenesis started earlier but progressed more slowly in LT than in LR, the stages ranging from 49 dpc to 77 dpc are critical for formation of different muscle phenotypes. 595 differentially expressed myogenesis genes were identified, and their roles in myogenesis were discussed. Furthermore, GSK3B, IKBKB, ACVR1, ITGA and STMN1 might contribute to later myogenesis and more muscle fibers in LR than LT. Some myogenesis inhibitors (ID1, ID2, CABIN1, MSTN, SMAD4, CTNNA1, NOTCH2, GPC3 and HMOX1) were higher expressed in LT than in LR, which might contribute to more slow muscle differentiation in LT than in LR. We also identified several genes which might contribute to intramuscular adipose differentiation. Most important, we further proposed a novel model in which MyoD and MEF2A controls the balance between intramuscular adipogenesis and myogenesis by regulating CEBP family; Myf5 and MEF2C are essential during the whole myogenesis process while MEF2D affects muscle growth and maturation. The MRFs and MEF2 families are also critical for the phenotypic differences between the two pig breeds. Overall, this study contributes to elucidating the mechanism underlying muscle development, which could provide valuable information for pig meat quality improvement. The raw data have been submitted to Gene Expression Omnibus (GEO) under series GSE25406.", "claims": [{"id": "claim-method", "statement": "Fisher's exact test and a χ2 test were used to classify the GO category, and the false discovery rate (FDR) was calculated to correct the P-value; the smaller the FDR, the smaller the error in judging the p-value.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:21637832", "locator": "sections[39] \"GO analysis\" paragraph 2", "quote": "Fisher's exact test and a χ2 test were used to classify the GO category, and the false discovery rate (FDR) was calculated to correct the P-value; the smaller the FDR, the smaller the error in judging the p-value."}], "state": "complete"}, {"id": "claim-data-access", "statement": "The raw data have been submitted to Gene Expression Omnibus (GEO) under series GSE25406.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:21637832", "locator": "sections[0] \"Abstract\" paragraph 2", "quote": "The raw data have been submitted to Gene Expression Omnibus (GEO) under series GSE25406."}], "state": "complete"}, {"id": "claim-target-result", "statement": "Digital gene expression identified 595 differentially expressed myogenesis genes using TPM-normalized tag counts with significance defined as P<0.005, FDR<0.01, and absolute log2-fold change >0.5, and the raw data were submitted to GEO under series GSE25406.", "kind": "quantitative", "target": {"value": 2.0, "unit": "fold change", "tolerance": null, "source": "To compare the differential expression of genes across samples, the number of raw clean tags in each sample was normalized to Tags per Million (TPM) to obtain normalized gene expression levels. Differential expression of genes or tags across samples was detected according to previously described methods [75]. Genes were deemed significantly differentially expressed with a P-value<0.005, a false discovery rate (FDR) <0.01 and an estimated absolute log2-fold change >0.5 in sequence counts across libraries."}, "evidence": [{"asset": "source:pmid:21637832", "locator": "sections[0] \"Abstract\" paragraph 1", "quote": "595 differentially expressed myogenesis genes were identified, and their roles in myogenesis were discussed."}, {"asset": "source:pmid:21637832", "locator": "sections[0] \"Abstract\" paragraph 2", "quote": "The raw data have been submitted to Gene Expression Omnibus (GEO) under series GSE25406."}, {"asset": "source:pmid:21637832", "locator": "sections[36] \"Differential expression\" paragraph 2", "quote": "To compare the differential expression of genes across samples, the number of raw clean tags in each sample was normalized to Tags per Million (TPM) to obtain normalized gene expression levels. Differential expression of genes or tags across samples was detected according to previously described methods [75]. Genes were deemed significantly differentially expressed with a P-value<0.005, a false discovery rate (FDR) <0.01 and an estimated absolute log2-fold change >0.5 in sequence counts across libraries."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "complete", "objective": "Digital gene expression identified 595 differentially expressed myogenesis genes using TPM-normalized tag counts with significance defined as P<0.005, FDR<0.01, and absolute log2-fold change >0.5, and the raw data were submitted to GEO under series GSE25406.", "steps": [{"id": "read_published_estimate", "order": 1, "name": "Read the published estimate from the source", "kind": "other", "state": "complete", "blocked_on": null, "description": "The claim, method and evidence span published with the source article.", "duration_hours": null, "parameters": {}, "inputs": ["source:pmid:21637832"], "outputs": ["the published estimate and its evidence"], "instruments": [], "reagents": []}, {"id": "obtain_source_data", "order": 2, "name": "Obtain the source data from the named repository", "kind": "data-preparation", "state": "complete", "blocked_on": null, "description": "The GEO series accession resolves and its sample accessions are listed. This lane did not download the dataset or read its licence.", "duration_hours": null, "parameters": {}, "inputs": ["data:geo:GSE25406"], "outputs": ["the source data with its repository provenance"], "instruments": [], "reagents": []}, {"id": "run_analysis", "order": 3, "name": "Run the published analysis", "kind": "analysis", "state": "complete", "blocked_on": null, "description": "Apply the paper's statistical method to the resolved series and compare the result with the published target.", "duration_hours": null, "parameters": {}, "inputs": ["data:geo:GSE25406"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:21637832", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1371/journal.pone.0019774", "sha256": "sha256:faeb3004cfdc6819910d121a122fbfb25afdbc63b7e77781b36260fd1ed56849", "bytes": 88158, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1371/journal.pone.0019774", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC3102668", "source": "europe_pmc_fulltext", "coverage": "whole_body", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE25406", "role": "raw", "format": "geo-series", "media_type": "application/vnd.ncbi.geo-series+json", "uri": "https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE25406", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE25406", "version": null, "note": "The raw data have been submitted to Gene Expression Omnibus (GEO) under series GSE25406.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The raw data have been submitted to Gene Expression Omnibus (GEO) under series GSE25406.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM623882", "GSM623883", "GSM623884", "GSM623885", "GSM623886", "GSM623887", "GSM623888", "GSM623889", "GSM623890", "GSM623891", "GSM623892", "GSM623893", "GSM623894", "GSM623895", "GSM623896", "GSM623897", "GSM623898", "GSM623899", "GSM623900", "GSM623901"], "sample_count": 20, "group_sizes": null, "group_sizes_basis": null, "platform": "9126", "taxon": "Sus scrofa", "title": "Comparative analyses of transcriptome during skeletal muscle development between pig breeds differing in muscle growth rate by deep-sequencing", "resolution_checked_at": "2026-11-01T15:39:54Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}, {"id": "asset:table:pmid:21637832:1", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:9315218656c0cbf5a49932f52b7f152d90040d118598baa6fd33a96804416fef", "bytes": 1461, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[7]/table-wrap[1]#pone-0019774-t001", "version": null, "note": "The numbers of differentially expressed genes between libraries.", "structure": {"pmcid": "PMC3102668", "section_index": 7, "section_title": "Temporal analysis of differentially expressed genes", "row_count": 20, "column_count": 7}, "table": {"label": "Table 1", "caption": "The numbers of differentially expressed genes between libraries.", "rows": [["The comparison between different libraries", "numbers of differentially expressed genes"], ["", "Total DE genes", "Total DE genes-up", "Total DE genes-down", "MD related DE genes", "MD related DE genes-up", "MD related DE genes-down"], ["LR-49 dpc/35 dpc", "2461", "876", "1585", "220", "109", "111"], ["LR-63 dpc/49 dpc", "462", "179", "282", "87", "58", "29"], ["LR-77 dpc/63 dpc", "241", "110", "131", "66", "36", "30"], ["LR-91 dpc/77 dpc", "423", "256", "167", "63", "35", "28"], ["LR-2 dpn/91 dpc", "2006", "771", "1235", "220", "112", "108"], ["LR-28 dpn/2 dpn", "1269", "285", "984", "149", "94", "55"], ["LR-90 dpn/28 dpn", "678", "235", "443", "129", "51", "78"], ["LR-120 dpn/90 dpn", "1373", "808", "565", "178", "56", "122"], ["LR-180 dpn/120 dpn", "702", "247", "455", "85", "44", "41"], ["LT-49 dpc/35 dpc", "1362", "863", "499", "147", "100", "47"], ["LT-63 dpc/49 dpc", "702", "341", "361", "108", "68", "40"], ["LT-77 dpc/63 dpc", "330", "185", "145", "68", "35", "33"], ["LT-91 dpc/77 dpc", "316", "163", "153", "68", "44", "24"], ["LT-2 dpn/91 dpc", "2301", "1067", "694", "271", "93", "178"], ["LT-28 dpn/2 dpn", "1453", "1036", "417", "199", "71", "128"], ["LT-90 dpn/28 dpn", "1036", "573", "463", "134", "66", "68"], ["LT-120 dpn/90 dpn", "642", "357", "285", "97", "57", "40"], ["LT-180 dpn/120 dpn", "403", "124", "279", "90", "51", "39"]], "row_count": 20, "column_count": 7, "numeric_cells": 126, "header_rows": 2, "locator": "section[7]/table-wrap[1]#pone-0019774-t001"}}, {"id": "asset:table:pmid:21637832:2", "role": "table", "format": "jats-table", "media_type": "application/vnd.jats.table+json", "uri": null, "sha256": "sha256:7a3d84e6376d8e558aff93fd529477c2d84d4cc950bce3a0b9f7a71dfcfd316d", "bytes": 2094, "access": "open", "state": "verified", "reader": "read_table", "reader_status": "available", "identifier": "section[9]/table-wrap[1]#pone-0019774-t002", "version": null, "note": "Signaling pathways of DE genes.", "structure": {"pmcid": "PMC3102668", "section_index": 9, "section_title": "Pathway analysis of temporal DE genes", "row_count": 20, "column_count": 19}, "table": {"label": "Table 2", "caption": "Signaling pathways of DE genes.", "rows": [["Pathway", "49/35 dpc", "63/49 dpc", "77/63 dpc", "91/77 dpc", "2 dpn/91 dpc", "28/2 dpn", "90/28 dpn", "120/90 dpn", "180/120 dpn"], ["", "LR", "LT", "LR", "LT", "LR", "LT", "LR", "LT", "LR", "LT", "LR", "LT", "LR", "LT", "LR", "LT", "LR", "LT"], ["PPAR", "↑", "", "", "↑", "", "", "", "", "↑", "↑", "↓", "↑", "↓", "↓", "", "↑", "↑", "↑"], ["MAPK", "↑", "↑", "↑", "", "↑", "↑", "↓", "", "↑", "↑↓", "↓", "↑↓", "↑", "↑↓", "↑↓", "↑↓", "", "↓"], ["Calcium", "↑↓", "↑", "", "", "", "", "", "↑", "↑", "↓", "↓", "↓", "↓", "↑", "↓", "↓", "", ""], ["Chemokine", "↓", "", "", "", "", "", "", "", "", "↓", "", "", "", "", "↓", "↓", "", ""], ["Wnt", "↓", "", "", "", "", "", "", "", "↓", "↓", "↓", "↓", "", "↑", "↓", "↓", "", ""], ["GnRH", "↓", "↑", "↑", "↓", "", "↑", "", "↓", "↑", "↓", "↑", "↓", "↑", "↑", "↓", "↓", "", ""], ["Insulin", "↓", "", "", "↓", "", "", "", "↑", "↑", "↑", "↓", "↓", "↓", "↑↓", "↓", "↓", "↑", "↑"], ["TGF-beta", "↓", "↑↓", "↑", "↓", "↓", "↑", "", "↓", "↓", "↓", "↓", "↓", "", "", "↑↓", "↑↓", "", ""], ["Neurotrophin", "↓", "↑↓", "", "↓", "↓", "", "", "", "↑", "↓", "↓", "↓", "↑", "↑↓", "↓", "↓", "", ""], ["Notch", "↓", "", "", "", "", "↓", "↓", "", "", "↓", "", "↓", "↓", "", "", "", "", ""], ["Adipocytokine", "", "", "", "", "", "", "↑", "", "↑", "", "↓", "", "↓", "", "↓", "", "↑", ""], ["Pentose phosphate", "", "↑", "", "↑", "", "↑", "", "", "↑", "", "↑", "↑", "↑", "↑", "", "", "", "↓"], ["p53", "", "", "", "", "", "", "", "", "↑", "↑", "↓", "↓", "", "↓", "", "", "", ""], ["VEGF", "", "", "", "", "", "", "", "", "↑", "↑", "↑", "", "", "", "↓", "↓", "", ""], ["Jak-STAT", "", "", "", "", "", "", "", "", "↓", "↓", "", "", "↓", "", "", "", "", ""], ["ErbB", "", "", "", "", "", "", "", "", "", "", "", "↓", "↑", "↑", "", "", "↓", ""], ["mTOR", "", "", "", "", "", "↓", "", "", "", "", "", "", "", "", "↓", "↓", "", ""], ["Fatty acid metabolism", "", "↑", "", "↑", "", "", "", "", "", "↑", "", "↑", "", "↑", "", "↓", "", ""]], "row_count": 20, "column_count": 19, "numeric_cells": 10, "header_rows": 2, "locator": "section[9]/table-wrap[1]#pone-0019774-t002"}}], "gaps": ["asset:data:geo:GSE25406 not downloaded", "asset:data:geo:GSE25406 licence not read"]}, "analysis": {"state": "partial", "methods": [{"id": "method-1", "family": "contingency_analysis", "specific": "chi_square", "label": "Categorical (contingency table) analysis", "dsl_operations": ["chi_square_test"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["categorical_cross_tabulation"], "evidence": {"score": 1.0, "band": "verbatim", "quote": "Fisher's exact test and a χ2 test were used to classify the GO category, and the false discovery rate (FDR) was calculated to correct the P-value; the smaller the FDR, the smaller the error in judging the p-value."}}, {"id": "method-2", "family": "correlation_analysis", "specific": "pearson", "label": "Correlation analysis", "dsl_operations": ["pearson_r"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["paired_numeric_columns"], "evidence": {"score": 0.8, "band": "high", "quote": "Pearson's correlation"}}, {"id": "method-3", "family": "contingency_analysis", "specific": "fisher_exact", "label": "Categorical (contingency table) analysis", "dsl_operations": ["contingency_table", "chi_square_test"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["categorical_cross_tabulation"], "evidence": {"score": 0.8, "band": "high", "quote": "Fisher's exact"}}], "targets": [{"id": "target-published-estimate", "claim_id": "claim-target-result", "operation": "chi_square_test", "expected_value": 2.0, "unit": "fold change", "tolerance": null, "state": "complete", "note": "The published estimate as the reader quoted it from the paper. No tolerance is published because the paper states none."}, {"id": "target-published-estimate-p-value", "claim_id": "claim-target-result", "operation": "chi_square_test", "expected_value": 0.005, "unit": "p-value", "tolerance": null, "state": "complete", "note": "The threshold the paper applied to its own result, as the paper states it. A run that reproduces the effect is judged against this threshold."}], "gaps": ["method.reader_differs_from_extractor"]}, "capability": {"state": "complete", "requirements": [{"capability_id": "capability:computational_analysis", "label": "Computational analysis platform", "assay_families": ["meta_analysis", "machine_learning", "survival_analysis", "statistical_regression", "pca", "molecular_docking", "network_analysis", "differential_expression", "gwas"], "availability_state": "owned", "tier_unlocked": "R1", "candidate": "Automated Assays Laravel application and simulated adapters"}], "gaps": []}, "readiness": {"state": "complete", "tier": "R1", "blocking_axis": "none", "blocked_on": "nothing:replicable_now", "axes": {"method": {"tier": "R1", "expressible": true, "required_operations": ["chi_square_test"], "missing_operations": [], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: chi_square."}, "capability": {"tier": "R1", "capability_id": "capability:computational_analysis", "availability_state": "owned", "basis": "The computational analysis is owned. The physical bench the paper used is not, so a physical validation cannot be run here."}, "data": {"tier": "R1", "reachable": true, "accessions": ["GSE25406"], "resolved_accessions": ["GSE25406"], "statements": ["The raw data have been submitted to Gene Expression Omnibus (GEO) under series GSE25406."], "basis": "The paper names GSE25406 in Gene Expression Omnibus (GEO) and the accession resolves in the repository. No dataset was downloaded and no licence was read."}}, "can_run": {"dsl_operations": ["chi_square_test"], "assets": ["source:pmid:21637832", "data:geo:GSE25406", "asset:table:pmid:21637832:1", "asset:table:pmid:21637832:2"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE25406"]}, "gaps": []}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-10-31T18:20:56Z", "independence": "same_provider"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 8, "usage": {"prompt_tokens": 179357, "completion_tokens": 4231, "cached_tokens": 134942}}, "human_critique": null, "stages": [{"id": "extraction", "state": "complete", "reader": "MiniMax-M3.1-Flash", "provider": "minimax"}], "gaps": ["review.independent_reader", "review.machine_critique", "review.human_critique"]}, "execution": {"state": "absent", "status": "not_run", "mode": null, "runs": [], "gaps": ["execution.not_run"]}, "provenance": {"entities": [{"id": "pmid:21637832", "type": "paper", "label": "pmid:21637832"}, {"id": "source:pmid:21637832", "type": "manuscript", "label": "PMID 21637832 full text"}, {"id": "data:geo:GSE25406", "type": "dataset", "label": "GEO GSE25406"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:21637832", "type": "paper", "label": "pmid:21637832"}, {"id": "source:pmid:21637832", "type": "manuscript", "label": "PMID 21637832 full text"}, {"id": "data:geo:GSE25406", "type": "dataset", "label": "GEO GSE25406"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-21637832", "type": "extraction", "used": ["source:pmid:21637832"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE25406", "type": "accession_resolution", "used": ["data:geo:GSE25406"], "agent": "agent:corpus-assay-engine"}], "agents": [{"id": "agent:corpus-assay-engine", "type": "software", "label": "Corpus assay engine"}, {"id": "agent:minimax-reader", "type": "software", "label": "MiniMax literature reader"}], "source_run_id": "reading-sci-2026-11-01", "dsl_catalogue_fingerprint": "343fb30756ba6ec3", "capability_table_fingerprint": "2cb3e42f047d2ecf", "snapshot_id": "ba34b49490d9bfaf"}, "lineage": {"parent_id": "pmid:21637832", "parent_hash": "sha256:6ac0a8ef2ac42c5466f152f691fd0ab21b3114157d28d6e45296887265c337e8", "fork_of": null, "generation": 7}, "gaps": ["asset:data:geo:GSE25406 licence not read", "asset:data:geo:GSE25406 not downloaded", "execution.not_run", "method.reader_differs_from_extractor", "review.human_critique", "review.independent_reader", "review.machine_critique"], "assays": [{"role": "data_generation", "definition_id": "microarray", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.93, "evidence": ["rmation may be lost in similar studies using Microarray or LongSAGE. Tag mapping and gene annotation"], "order_basis": "paper_text_order:7979", "id": "assay:1"}, {"role": "data_generation", "definition_id": "qpcr", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.93, "evidence": ["iased analysis of the transcriptome. And the qPCR confirmation results indicated that the Sole"], "order_basis": "paper_text_order:8268", "id": "assay:2"}, {"role": "analysis", "definition_id": "contingency_analysis", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": ["chi_square_test"], "support_state": "supported", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 1.0, "evidence": ["Fisher's exact test and a χ2 test were used to classify the GO category, and the false discovery rate (FDR) was calculated to correct the P-value; the smaller the FDR, the smaller the error in judging the p-value."], "order_basis": "consumes_earlier_output", "id": "assay:3"}]}