Evidence card: NCBI:GeneID:105318810

Real public-study evidence. Association evidence only; not a validated biomarker.

How to interpret expression as evidence for resilience. Differential expression alone does not establish a resilience marker; higher expression does not necessarily mean greater resilience.

This feature was selected among up to ten displayed findings in at least one study. All available rows for it are shown below, including nonsignificant and exploratory results. Presence in multiple datasets does not establish independent replication: exposures, tissues, life stages, and phenotype meanings differ. No cross-study score or pooled effect is calculated.

Guide to effect sizes, uncertainty and gene descriptions.

Study-specific effects

Study / design RefSeq description Contrast Effect scale Effect Standard error Adjusted p-value Interpretation
PRJEB18614 uncharacterized LOC105318810 alexandrium_minutum_vs_heterocapsa_triquetra_diet log2 fold change -0.3041 0.1829 0.2153 Study-level association; see design caveats
CGIG_HEAT24H_RNASEQ_PRJNA407831 uncharacterized LOC105318810 35C_24h_vs_0h_gill log2 fold change 0.4171 0.08788 0.0001186 Study-level association; see design caveats
CGIG_HEAT6H_RNASEQ_PRJNA407831 uncharacterized LOC105318810 35C_6h_vs_0h_gill log2 fold change 0.2854 0.07038 0.004011 Study-level association; see design caveats
PRJEB86646 uncharacterized LOC105318810 30C_vs_23C_fed_before_oshv1 log2 fold change -0.993 0.1435 1.28e-10 Study-level association; see design caveats
PRJNA516762 uncharacterized LOC105318810 heat_35C_2h_vs_control_12C log2 fold change 0.4041 0.1132 0.07249 Study-level association; see design caveats
CGIG_HEAT_RNASEQ_PRJNA593309 uncharacterized LOC105318810 29C_vs_21C_12h_oshv1_cohabitation log2 fold change -0.9845 0.4039 0.2819 Study-level association; see design caveats
CGIG_HEAT_RNASEQ_PRJNA913164 uncharacterized LOC105318810 30C_seawater_vs_20C_adult_ctenidium log2 fold change -0.4125 0.2039 0.1294 Study-level association; see design caveats
CGIG_MULTISTRESS_RNASEQ_PRJNA913164 uncharacterized LOC105318810 30C_seawater_then_44C_emersion_6h_vs_20C_adult_ctenidium log2 fold change -0.1069 0.1433 0.6766 Study-level association; see design caveats
PRJNA735889 uncharacterized LOC105318810 ph_6.5_6.8_vs_ph_7.4_7.8_23d_juveniles log2 fold change 0.1153 0.239 0.8116 (nominal only) Exploratory; tank clustering unmodeled
PRJNA826964 uncharacterized LOC105318810 ocean_acidification_vs_ambient_control log2 fold change 0.7875 0.2294 0.02319 Study-level association; see design caveats
PRJNA1329250 uncharacterized LOC105318810 oshv1_microvariant_challenge_vs_unchallenged_control log2 fold change -0.2498 0.2253 0.4721 Study-level association; see design caveats
CGIG_OSHV1_RNASEQ_PRJNA593309 uncharacterized LOC105318810 oshv1_cohabitation_48h_vs_0h_21C log2 fold change -1.264 0.6229 0.2799 Study-level association; see design caveats
PRJEB58194 uncharacterized LOC105318810 pesticide_mixture_0_48hpf_vs_control_F0_gastrula_dlarva log2 fold change -0.8066 0.3987 not available Study-level association; see design caveats
PRJNA756710 uncharacterized LOC105318810 low_salinity_10ppt_vs_30ppt_during_vibrio_infection log2 fold change 2.938 0.989 0.02653 Study-level association; see design caveats
PRJEB86618 uncharacterized LOC105318810 starved_vs_fed_ad_libitum_16d_before_oshv1 log2 fold change -2.707 0.1442 1.752e-75 Study-level association; see design caveats
CGIG_TIRELEACHATE_RNASEQ_PRJNA856813 uncharacterized LOC105318810 tire_leachate_high_vs_control_gill_day41 log2 fold change -0.2326 0.3547 0.9985 Study-level association; see design caveats
CGIG_TIREPARTICLE_RNASEQ_PRJNA856813 uncharacterized LOC105318810 tire_microparticles_high_vs_control_gill_day41 log2 fold change -0.05017 0.3397 0.9999 Study-level association; see design caveats
CGIG_VCRASSOSTREAE_RNASEQ_PRJNA515169 uncharacterized LOC105318810 vibrio_crassostreae_J2_9_vs_mock_injection_8h log2 fold change -0.3777 0.142 0.1855 Study-level association; see design caveats
CGIG_VIBRIO12H_RNASEQ_PRJNA852380 uncharacterized LOC105318810 vibrio_12h_vs_control_digestive_gland log2 fold change -0.8235 0.2098 0.002984 Study-level association; see design caveats
CGIG_VIBRIO24H_RNASEQ_PRJNA852380 uncharacterized LOC105318810 vibrio_24h_vs_control_digestive_gland log2 fold change -3.029 0.2856 6.601e-23 Study-level association; see design caveats
CGIG_VTASMANIENSIS_RNASEQ_PRJNA515169 uncharacterized LOC105318810 vibrio_tasmaniensis_LGP32_vs_mock_injection_8h log2 fold change -0.2175 0.156 0.7972 Study-level association; see design caveats

CGIG_HAB_RNASEQ_PRJEB18614

Phenotype interpretation: A. minutum exposure is unfavorable. The publication reports remodelling of the gill cycling transcriptome (loss of oscillations, de novo oscillating transcripts, shortened periods) and toxin accumulation, not survival. AREE records the time-averaged gill transcriptional response to the toxic diet.

Biological replication: 24 pooled RNA-seq libraries in the primary analysis (12 per diet), each pooling gill RNA from 6 oysters. The original 26 libraries cover 13 time-matched pairs; both day 2 ZT 8.5 pools are excluded, leaving 12 pairs. There is no replication within a diet-by-time cell. Each diet used three tanks; sample_size counts analysed pools, not individual oysters or independent tanks.

Evidence classification: stress_response. Tissue: gill; life stage: adult.

Quality flags: raw_reanalysis;harmful_algal_bloom_response;time_averaged_over_diurnal_cycle;timepoint_blocked_contrast;sequencing_pool_partly_confounded.

Study limitations:

  • Each library was sequenced on two lanes; the lanes are technical duplicates of one library, not replicates. Only the lower-numbered lane of each pair is used (26 of the 52 deposited runs, 136.8 GB instead of about 274 GB). The unused run is named in each manifest row.
  • Reads come from two mirrors. 12 runs are the ENA FASTQs, checked against the manifest MD5s. The other 14, listed with NCBI’s size and MD5 in CGIG_HAB_RNASEQ_PRJEB18614_sra_locations.tsv, were converted from full NCBI .sra archives, and their read counts were checked against the manifest. They were moved to NCBI after ENA corrupted multi-connection transfers and then refused requests. Each library comes from one mirror only.
  • Both day 2 ZT 8.5 gill pools (control and A. minutum) carry adductor-muscle RNA from dissection and are excluded from the primary contrast; sample_size counts the 24 remaining pools. The all-26 fit is kept as a comparison.
  • Salmon detected one ENA library (control d1 ZT 12.5, ERR1758944) as ISF instead of the ISR expected for stranded TruSeq, so its read files are swapped relative to the rest. Library type is inferred per sample, so quantification is unaffected.
  • The primary contrast, ~ timepoint + condition, blocks on sampling time and uses 12 time-matched pools per diet after excluding the contaminated pair. The all-26 sensitivity model and publication use 13 per diet. The publication’s central finding is that A. minutum remodels the daily cycling of gill transcripts, so a diet effect averaged over time will miss, or dilute, effects confined to part of the cycle.
  • Diet is partly confounded with sequencing pool. The libraries were multiplexed into four pools, each on its own lane pair of one flowcell. Pool L003-L004 holds 6 A. minutum libraries and no controls, and pool L005-L006 holds 6 controls and 1 A. minutum library; pools L001-L002 and L007-L008 carry both diets. The sensitivity model with sequencing_pool as a covariate keeps the direction of 99.9% of the primary significant genes, including all 30 top genes, which also stay significant. Only 53.8% of genes overall stay significant, so the gene count depends on the model, but the direction does not.
  • Each diet used three tanks, but the publication does not say how the 6 oysters in each pool were drawn from them, so tank effects cannot be modelled.
  • Timepoint labels are Day and Zeitgeber time rounded down: d1_zt04 is Day 1 ZT 4.5 and d3_zt04 is Day 3 ZT 4.5, where ZT 0 is lights on. The methods list the same 13 times as the ArrayExpress SDRF.
  • Condition, time and lane come from the ArrayExpress E-MTAB-5345 SDRF (feeding condition, time and assay name), and every ENA sample alias matches its SDRF source name.
  • The publication aligned reads with STAR to the Ensembl Metazoa C. gigas genome and mitochondrial genome and tested diet with edgeR across the 13 times, reporting 740 differentially expressed transcripts at FDR < 0.01 (89.5% higher under A. minutum). The transcript list is not in the supplement, whose only table lists the 30 most significant cycling transcripts per class, so only the count and direction balance can be compared. AREE finds about 5,070 genes at FDR < 0.01, 59.6% higher under A. minutum. That is far more genes and a weaker up-bias than the published 740 transcripts (89.5% up), with a different aligner, count model and genome assembly.
  • The strongest A. minutum-higher genes form a lysosomal and endocytic programme (dipeptidyl peptidase 1, LAMP1, CD63, VPS4B, syntaxin-8) and include detoxification genes: the ABCB1 multidrug efflux pump, arylacetamide deacetylase and mitochondrial alternative oxidase.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: 185cbdd6f509f43a79442fe2fd9137ea489711c4afe95d0f2a1a21c6c6b4762a.

CGIG_HEAT24H_RNASEQ_PRJNA407831

Phenotype interpretation: Intertidal-origin F1 oysters showed higher survival after acute heat stress (42 C, 1 h) and higher metabolic rates at 35 C than subtidal-origin oysters, with slower growth, and more transcriptional plasticity. AREE records the gill transcriptional response to 35 C averaged over the four populations, not survival per library.

Biological replication: 23 libraries analysed (24 sequenced; one 0 h outlier left out): 11 at 0 h and 12 at 24 h. Each library pools equal gill RNA from five oysters of one culture tank; four populations (two sites x intertidal or subtidal origin) x three tanks, and the same tanks are sampled at each time.

Evidence classification: stress_response. Tissue: gill; life stage: adult.

Quality flags: raw_reanalysis;acute_heat_response;pooled_libraries;population_blocked_contrast;shared_control_group;dependence_group=PRJNA407831_0h_pools;heat_24h;fdr_lt_0.05.

Study limitations:

  • Both AREE records from this study (the other is CGIG_HEAT6H_RNASEQ_PRJNA407831, 6 h) test against the same twelve 0 h pools, so they carry the dependence group PRJNA407831_0h_pools and count as one replication group.
  • Each library pools gill RNA of five oysters, so the error model describes tank-pool variation, not variation between animals. The sister IOCAS study PRJNA826964 deposited resequenced RNA as replicates, so these pools were checked with the expressed-SNP test (scripts/check_replicate_genotypes.py) before analysis.
  • The model blocks on the four source populations (~ population + condition), so the effect is averaged over intertidal and subtidal origins and the two sites; the publication’s main result, that intertidal oysters respond with more plasticity, is not part of this record. A tank-paired model (~ tank + condition) is kept as a sensitivity analysis.
  • Sample aliases code tidal origin as J or X (for example BYQJ01, LTX243); the publication uses I (intertidal) and S (subtidal), and the BioSamples do not record it. J and X are read as the Chinese jian (intertidal) and xia (subtidal), but no AREE result depends on that reading because population is a blocking factor.
  • One 0 h library, h0_LTJ_2, is left out by AREE’s QC rule (median VST correlation below the median minus three robust SDs); it is the library the publication also discarded (LT-I_0h_2). The exclusion lives in scripts/prepare_prjna407831_fastqs.py.
  • ENA served other studies at 0.05-0.35 MB/s in October 2026, so all runs come from NCBI .sra archives checked against NCBI’s MD5s, with read counts checked against the manifest on conversion.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: ba65a504fd1c43fe5df5636afd198fc3f2b95abd6a145604c6ed9ca3caf03421.

CGIG_HEAT6H_RNASEQ_PRJNA407831

Phenotype interpretation: Intertidal-origin F1 oysters showed higher survival after acute heat stress (42 C, 1 h) and higher metabolic rates at 35 C than subtidal-origin oysters, with slower growth, and more transcriptional plasticity. AREE records the gill transcriptional response to 35 C averaged over the four populations, not survival per library.

Biological replication: 23 libraries analysed (24 sequenced; one 0 h outlier left out): 11 at 0 h and 12 at 6 h. Each library pools equal gill RNA from five oysters of one culture tank; four populations (two sites x intertidal or subtidal origin) x three tanks, and the same tanks are sampled at each time.

Evidence classification: stress_response. Tissue: gill; life stage: adult.

Quality flags: raw_reanalysis;acute_heat_response;pooled_libraries;population_blocked_contrast;shared_control_group;dependence_group=PRJNA407831_0h_pools;heat_6h;fdr_lt_0.05.

Study limitations:

  • Both AREE records from this study (the other is CGIG_HEAT24H_RNASEQ_PRJNA407831, 24 h) test against the same twelve 0 h pools, so they carry the dependence group PRJNA407831_0h_pools and count as one replication group.
  • Each library pools gill RNA of five oysters, so the error model describes tank-pool variation, not variation between animals. The sister IOCAS study PRJNA826964 deposited resequenced RNA as replicates, so these pools were checked with the expressed-SNP test (scripts/check_replicate_genotypes.py) before analysis.
  • The model blocks on the four source populations (~ population + condition), so the effect is averaged over intertidal and subtidal origins and the two sites; the publication’s main result, that intertidal oysters respond with more plasticity, is not part of this record. A tank-paired model (~ tank + condition) is kept as a sensitivity analysis.
  • Sample aliases code tidal origin as J or X (for example BYQJ01, LTX243); the publication uses I (intertidal) and S (subtidal), and the BioSamples do not record it. J and X are read as the Chinese jian (intertidal) and xia (subtidal), but no AREE result depends on that reading because population is a blocking factor.
  • One 0 h library, h0_LTJ_2, is left out by AREE’s QC rule (median VST correlation below the median minus three robust SDs); it is the library the publication also discarded (LT-I_0h_2). The exclusion lives in scripts/prepare_prjna407831_fastqs.py.
  • ENA served other studies at 0.05-0.35 MB/s in October 2026, so all runs come from NCBI .sra archives checked against NCBI’s MD5s, with read counts checked against the manifest on conversion.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: 3ff5b26926b5dc5ac2b94196d5b5a8628742dc7a4ebf2460da80389567311c2a.

CGIG_HEAT_RNASEQ_PRJEB86646

Phenotype interpretation: Lower mortality after OsHV-1 mu-Var cohabitation is favorable. The publication reports that 30 C limits POMS development and improves survival, through metabolic reprogramming and stronger baseline and induced antiviral immunity. AREE records the whole-spat transcriptional state at 30 C before infection, not survival per library.

Biological replication: 12 RNA-seq libraries, all sampled at 0 h post-cohabitation: 6 at 23 C and 6 at 30 C, three per temperature in each of the two families. Each library pools the soft tissue of 10 live spat taken from one tank, as in the DECICOMP starvation arm.

Evidence classification: resilience_associated. Tissue: whole_soft_tissue; life stage: spat.

Quality flags: raw_reanalysis;temperature_response;pre_infection_baseline;family_blocked_contrast;shared_control_group_with_PRJEB86618;fdr_lt_0.05;dependence_group=DECICOMP_PRJEB86525_0h_controls.

Study limitations:

  • The six 23 C fed control libraries (PRJEB86525, 0 h) are the same libraries used as controls in CGIG_STARVATION_RNASEQ_PRJEB86618. The two harmonized studies therefore share a reference group, so their effect estimates are correlated and must not be treated as independent in a meta-analysis.
  • The contrast uses only the 0 h libraries, sampled when the OsHV-1 mu-Var donors were added, so it captures temperature without infection. The 36 later libraries (3-72 h) follow infection, which progresses far less at 30 C, and are not used.
  • Each tank held both families, so the F11N and F14R libraries with the same replicate number probably share a tank. The model ~ family + condition treats the six libraries per temperature as independent, although only three tanks per temperature are independent units.
  • The 30 C arm and the 23 C controls were held in separate tanks of one experiment at the Argenton station, so tank is nested in temperature.
  • Condition, family, time and replicate come from the ENA sample aliases and titles (sam_P_<3 fed 23 C | 4 fed 30 C>_), and the two sources agree for every library.
  • The publication reports enhanced baseline antiviral immunity at 30 C. At 0 h this reanalysis finds the opposite for several classic antiviral markers: all 13 significant interferon-induced protein 44 genes, 10 of 12 significant tumor necrosis factor genes and both significant 2’-5’-oligoadenylate synthetase genes are lower at 30 C. The publication is not open access and its gene lists were not compared, so whether this reflects different gene sets, a different time point or a real disagreement is unresolved.
  • Pools at 30 C map about three points less to the transcriptome than pools at 23 C. The shift is systematic but small beside the expression effects.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: abbc33a2dc6e48a0f0a903bbca1d4bce3f90c56927e80ff5c5bc3e88dd6c0c5c.

CGIG_HEAT_RNASEQ_PRJNA516762

Phenotype interpretation: Organism-level thermal-tolerance direction was not measured; molecular heat-response direction only

Biological replication: Direct heat-versus-control contrast: 3 pooled RNA-seq libraries per condition, each pooling 5 individuals (30 individuals total; 6 libraries)

Evidence classification: stress_response. Tissue: gill; life stage: juvenile.

Quality flags: raw_reanalysis;acute_heat_response;pooled_libraries.

Study limitations:

  • This intake represents molecular heat-stress response, not a direct survival or organism-level thermal-tolerance measurement.
  • Each RNA-seq library pooled gill RNA from five oysters, so individual-level variability is unavailable.
  • AREE independently reanalyzed the raw reads against GCF_963853765.1/RS_2024_06; the original study used legacy assembly GCA_000297895.1 and CGI/newGene identifiers, so these results are a current-reference reanalysis rather than an exact reconstruction of the publication pipeline.
  • The direct heat-versus-control contrast yields modest differential expression (99 genes at FDR <= 0.05 of 22,397 tested), and PC1 (26.3% of variance) separates most but not all libraries, with heat_A overlapping the controls.
  • The HSF1 RNAi-only and RNAi-plus-heat contrasts are preserved in the source workbook but excluded from this direct heat-versus-control reanalysis.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: f0bded9f7662170c4985ff78b455a69e7df0b43c114895da85277159a2daba88.

CGIG_HEAT_RNASEQ_PRJNA593309

Phenotype interpretation: Survival 14 days after cohabitation was higher at 29 C (85.7%) than at 21 C (52.4%), and viral load rose more slowly at 29 C. AREE records the whole-oyster transcriptional difference at 12 h, not survival per library.

Biological replication: 5 oysters from 6 libraries: three at 21 C, each from a different tank, and two at 29 C, because two of the three 29 C libraries are the same oyster and were summed. The publication describes one oyster per tank, so the duplicated 29 C oyster means either one tank was sampled twice or one RNA was sequenced twice.

Evidence classification: resilience_associated. Tissue: whole_organism; life stage: juvenile.

Quality flags: raw_reanalysis;temperature_response;during_oshv1_infection;duplicate_animal_libraries_collapsed;one_oyster_per_tank;dependence_group=PRJNA593309_21C_tanks.

Study limitations:

  • Both arms are infected, so the contrast is temperature in the context of early OsHV-1 infection, not temperature alone; before infection the publication finds only nine genes differing between 29 and 21 C.
  • One whole oyster per tank was sequenced at each time point, with three replicate tanks per temperature (45 l tanks), so libraries are tank-level replicates and the contrast has three tanks per arm.
  • The two AREE records from this study use different oysters but the same three 21 C tanks (the 21 C 12 h oysters of CGIG_HEAT_RNASEQ_PRJNA593309 and the 21 C 0 h and 48 h oysters of CGIG_OSHV1_RNASEQ_PRJNA593309), so both carry the dependence group PRJNA593309_21C_tanks and count as one replication group.
  • The BioSamples give a collection date of 2019-05-19, but the publication dates the experiment to May 2016; the experimental groups follow the sample aliases (Cg_<CTR | OsHV>hpc_), which agree with the publication’s design.
  • The study also deposited libraries at 26 C and at 24 and 96 h; the publication does not analyse the 26 C libraries, and AREE uses only the time points of its two contrasts.
  • The publication aligned reads with STAR to the Ensembl oyster genome v9.38 and tested with edgeR; AREE quantifies against GCF_963853765.1/RS_2024_06 with Salmon and DESeq2.
  • Two of the three 29 C libraries at 12 h are the same oyster by expressed-SNP genotype, so this record has two independent 29 C oysters, little power, and fewer significant genes than the publication, whose three-replicate analysis counted that oyster twice.
  • Mapping is about eight points lower at 29 C than at 21 C at 12 h, a systematic shift that may reflect more viral or unspliced RNA at 29 C or a batch difference.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: 539c129666a98b1b71f771636b53a7cd9cbc55ce3de8bb6fe3f1a9a51bfa4c30.

CGIG_HEAT_RNASEQ_PRJNA913164

Phenotype interpretation: Mortality was elevated in the stress treatments relative to the control. AREE records the ctenidium transcriptional response to elevated seawater temperature, not survival per library.

Biological replication: 44 ctenidium libraries analysed, one oyster each (23 controls and 21 treated after QC exclusions), from 48 sequenced: 24 controls and 24 30 C, 12 diploid and 12 triploid in each arm. The 24 control libraries are shared with CGIG_MULTISTRESS_RNASEQ_PRJNA913164 (30 C followed by 44 C emersion).

Evidence classification: stress_response. Tissue: gill; life stage: adult.

Quality flags: raw_reanalysis;quantseq_3prime_tag_seq;ploidy_blocked_contrast;shared_control_group;dependence_group=PRJNA913164_20C_controls;temperature_response.

Study limitations:

  • The BioSamples carry no treatment or ploidy attribute. Both come from the description column of the authors’ NCBI submission sheet (github.com/mattgeorgephd/NOPP-gigas-ploidy-temp, 202107_EXP2/tag-seq/NCBI/attributes.tsv), matched to BioSample by isolate name; sample-name prefixes do not map one-to-one to groups.
  • The study has three arms that share one control group. AREE records each stressor against the 20 C controls as its own study (CGIG_MULTISTRESS_RNASEQ_PRJNA913164 (30 C followed by 44 C emersion) is the other), so the two records share their reference group and their effect estimates are correlated; they must not be treated as independent in a meta-analysis.
  • Diploids and triploids are pooled with ploidy as a blocking factor (~ ploidy + condition), so the effect is averaged over ploidy. The publication’s main finding is that triploids respond differently, especially to the multiple stressor; the ploidy-by-treatment interaction is reported in the quality-control summary, and ploidy-specific effects are not part of this record.
  • QuantSeq 3’ tag-seq yields about one read per transcript, so Salmon was run with –noLengthCorrection and DESeq2 used the counts without the tximport length offset (–tag-seq=yes). Library sizes are small (0.9-8.3 million reads), and isoform-level information is not recoverable.
  • The authors excluded eight libraries as outliers (D54, N56, X44, R53, M43, N54, X42, T62). AREE quantified all 72 and applied its own rule (median VST correlation below the median minus three robust SDs): seven of the authors’ eight fail it, D54 passes and is kept, and R62 and N57 also fail and are left out. The exclusions live in scripts/prepare_prjna913164_fastqs.py and remove one or two libraries from some ploidy-by-treatment cells.
  • The expressed-SNP replicate check (scripts/check_replicate_genotypes.py) could not be run: 3’ tags cover too little of the target transcripts (about 2,000 usable reads per library, even after poly-A and adapter trimming) to call allele frequencies. Each library is a single adult oyster according to the publication, so no pooling or resequencing is expected, but this was not verified genetically.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: f77e0e6f5683b7fe006272a6c53cb16bb293cf3542164602272e155aa0f4bafb.

CGIG_MULTISTRESS_RNASEQ_PRJNA913164

Phenotype interpretation: Mortality after the multiple stressor was higher in triploids than diploids (36.4% versus 14.8%), with metabolic depression and reduced ctenidium Na/K-ATPase activity in triploids. AREE records the ctenidium transcriptional response averaged over ploidy, not survival per library.

Biological replication: 42 ctenidium libraries analysed, one oyster each (23 controls and 19 treated after QC exclusions), from 48 sequenced: 24 controls and 24 30 C + 44 C emersion, 12 diploid and 12 triploid in each arm. The 24 control libraries are shared with CGIG_HEAT_RNASEQ_PRJNA913164 (30 C seawater only).

Evidence classification: stress_response. Tissue: gill; life stage: adult.

Quality flags: raw_reanalysis;quantseq_3prime_tag_seq;ploidy_blocked_contrast;shared_control_group;dependence_group=PRJNA913164_20C_controls;multiple_stressor_response.

Study limitations:

  • The BioSamples carry no treatment or ploidy attribute. Both come from the description column of the authors’ NCBI submission sheet (github.com/mattgeorgephd/NOPP-gigas-ploidy-temp, 202107_EXP2/tag-seq/NCBI/attributes.tsv), matched to BioSample by isolate name; sample-name prefixes do not map one-to-one to groups.
  • The study has three arms that share one control group. AREE records each stressor against the 20 C controls as its own study (CGIG_HEAT_RNASEQ_PRJNA913164 (30 C seawater only) is the other), so the two records share their reference group and their effect estimates are correlated; they must not be treated as independent in a meta-analysis.
  • Diploids and triploids are pooled with ploidy as a blocking factor (~ ploidy + condition), so the effect is averaged over ploidy. The publication’s main finding is that triploids respond differently, especially to the multiple stressor; the ploidy-by-treatment interaction is reported in the quality-control summary, and ploidy-specific effects are not part of this record.
  • QuantSeq 3’ tag-seq yields about one read per transcript, so Salmon was run with –noLengthCorrection and DESeq2 used the counts without the tximport length offset (–tag-seq=yes). Library sizes are small (0.9-8.3 million reads), and isoform-level information is not recoverable.
  • The authors excluded eight libraries as outliers (D54, N56, X44, R53, M43, N54, X42, T62). AREE quantified all 72 and applied its own rule (median VST correlation below the median minus three robust SDs): seven of the authors’ eight fail it, D54 passes and is kept, and R62 and N57 also fail and are left out. The exclusions live in scripts/prepare_prjna913164_fastqs.py and remove one or two libraries from some ploidy-by-treatment cells.
  • The expressed-SNP replicate check (scripts/check_replicate_genotypes.py) could not be run: 3’ tags cover too little of the target transcripts (about 2,000 usable reads per library, even after poly-A and adapter trimming) to call allele frequencies. Each library is a single adult oyster according to the publication, so no pooling or resequencing is expected, but this was not verified genetically.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: d2d6e4a07aa73594d97a4c9c4a593cace73fecddc5c11304a2bbef65277282d6.

CGIG_OA_RNASEQ_PRJNA735889

Phenotype interpretation: Lutier et al. report physiological tipping points at pHT 7.3-6.9 that coincide with a major reshuffling of the transcriptome and membrane lipids; AREE records the whole-animal transcriptional response below the tipping window relative to above it, not a measured tolerance phenotype

Biological replication: 50 oysters, one whole-animal library each: 25 control and 25 acidified, 5 per tank. A 51st library, 13Ebis, is a resequencing of oyster 13E and is not analysed. Each pH level is one tank, so oysters are replicates within tank and the 10 tanks are the replicates of the pH contrast.

Evidence classification: environmental_exposure. Tissue: whole_organism; life stage: juvenile.

Quality flags: raw_reanalysis;ocean_acidification_response;ph_band_contrast_tipping_window_excluded;no_ambient_ph_reference;one_tank_per_ph_level;gametogenesis_blocked_contrast;exploratory_tank_clustering_unmodeled.

Exploratory uncertainty: Tank clustering is unmodeled. Standard errors and adjusted p-values above are nominal only; this row is excluded from pooling and significance rewards.

Study limitations:

  • The study is a 15-level pH gradient (pHT 6.5-7.8, 76 libraries). AREE reduces it to a two-group contrast at the authors’ tipping-point window: the five tanks inside pHT 7.3-6.9 (tanks 1, 6, 4, 11 and 10; 25 libraries) are excluded and not downloaded. Effect sizes are therefore log2 fold changes between two pH bands, not slopes per pH unit, and they mix the response across five pH levels on each side.
  • There is no ambient-pH tank; the reference band (pHT 7.4-7.8) is itself below present-day coastal pH, so the contrast understates the full response to acidification.
  • Each pH level is a single tank, so tank and pH are confounded at the level of a single pH; with five tanks per band the band contrast has tank-level replication, but DESeq2 treats the 50 analysed oysters as independent and does not model tank. Standard errors, p-values and FDR values are exploratory and may overstate precision. A tank-aware reanalysis is required before inferential use; PCA and per-tank residual checks do not correct the error model.
  • Tank 15 is labelled ‘pHT 6.5 bis’ in BioSample: a second tank at pHT 6.5 alongside tank 9.
  • Library 13Ebis (SRR14803132) has a different index from 13E (SRR14803133) and about half its reads. The expressed-SNP check shows it is the same oyster (allele-frequency difference 0.026, below the 0.032 expected for one RNA resequenced; the next closest pair differs by 0.198), so it is kept in the manifest but left out of the analysis.
  • Male gametogenesis is unbalanced between bands: a spermatogenesis module is active in 11 of 25 control and 2 of 25 acidified oysters (Fisher p = 0.008), spread across all five control tanks, which suggests that low pH delays or suppresses the onset of spermatogenesis in 5-month-old juveniles. The primary model blocks on this call (~ gametogenesis + condition) so the evidence records a gametogenesis-adjusted association with pH, without establishing a direct causal response; the unadjusted model is kept as a sensitivity analysis. The call is a transcriptional proxy derived from the same data (scripts/score_prjna735889_gametogenesis.py), not a histological sex or maturation call, and with only two active acidified oysters the gametogenesis effect is estimated mostly from the control band.
  • Condition, tank and pH come from the NCBI BioSample title, treatment and growth_protocol attributes; every library is paired-end. ENA served the FASTQs at about 80 KB/s on 2026-10-02, so all 51 runs are taken from the NCBI Open Data mirror as .sra archives checked against the size and MD5 NCBI publishes (CGIG_OA_RNASEQ_PRJNA735889_sra_locations.tsv), and each converted FASTQ pair must match the manifest read count.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: a8f7485ec9436cc334d22b6819d7137f37dccb753fa688700c49d264cb7c0dda.

CGIG_OA_RNASEQ_PRJNA826964

Phenotype interpretation: The BioProject title frames acidification as compromising energy management; AREE records the hepatopancreas transcriptional response to acidification, not a measured organism-level energetic or survival phenotype

Biological replication: 6 pooled RNA samples, one per arm at each of 7, 28 and 56 days, each sequenced as three libraries (18 libraries); the BioSample records call the libraries biological replicates, but their genotypes show each set of three is one pool of several oysters

Evidence classification: stress_response. Tissue: hepatopancreas; life stage: not_recorded.

Quality flags: raw_reanalysis;ocean_acidification_response;timepoint_blocked_contrast;pooled_rna_samples;possible_genetic_background_confound;fdr_lt_0.05.

Study limitations:

  • No publication is linked to the BioProject, and searches of PubMed, Europe PMC full text and Crossref on 2026-09-28 found none, so there are no published methods, results or effect tables to reconcile the reanalysis against.
  • The acidification intensity (target pH or pCO2), seawater carbonate chemistry, temperature and oyster size or age are not recorded in the public metadata, so exposure intensity and life stage are unresolved.
  • The treatment was probably CO2-driven acidification to pH 7.50. This is inferred, not documented: the same IOCAS group (L. Qiu and colleagues) ran C. gigas CO2 exposures at pH 7.50 sampled from 7 to 56 days (Wang et al. 2021, Ecotoxicol Environ Saf, doi:10.1016/j.ecoenv.2021.112235) and at pH 7.50 for 7 and 60 days (Wang et al. 2020, Sci Total Environ, doi:10.1016/j.scitotenv.2020.140177), and reported hepatopancreas responses to acidification over about two months (Zhang et al. 2021, J Oceanol Limnol, doi:10.1007/s00343-021-0462-x). None of these papers is linked to PRJNA826964, and the closed-access hepatopancreas paper, the most likely companion, has not been checked against it.
  • The primary AREE contrast compares the three acidified pools with the three control pools and blocks on sampling day (~ timepoint + condition); it estimates an acidification effect averaged over 7, 28 and 56 days rather than time-specific or acclimation effects, and it has only two residual degrees of freedom.
  • The three libraries per day-by-condition cell are technical replicates of one pooled RNA sample, not the independent biological replicates the BioSample attributes state. Expressed-SNP allele frequencies (scripts/check_replicate_genotypes.py, 4M read pairs per library) agree within a cell to within resequencing noise and show no single-animal genotype peaks, and within-cell expression dispersion is 2.7% BCV against 10.8-22.4% in the other AREE oyster studies. Treating the 18 libraries as replicates inflates significance roughly twelvefold, so the libraries are summed into 6 pools (–collapse-by=condition+timepoint).
  • Each pool is a different set of animals, and the acidified pools are genetically further from the controls (mean allele-frequency difference 0.24-0.31) than the control pools are from one another (0.085-0.19), so acidified and control oysters may come from different stocks or families. Several of the strongest effects are 4-6 log2 changes in transposon-derived genes (tigger, HARBI1), the pattern expected from presence-absence variation between genetic backgrounds, so part of the contrast may be genetic rather than an acidification response. The evidence is harmonized with every row flagged pooled_rna_samples and possible_genetic_background_confound, so scoring discounts it and evidence cards list the caveat; treat it as provisional acidification evidence until a design with shared genetic backgrounds confirms it.
  • Condition, sampling day and replicate were taken from NCBI BioSample treatment and replicate attributes and agree with the library names (C/O = control/acidified, followed by day and replicate).
  • The BioProject carries no author-deposited processed expression matrix.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: 4a01845d837665c4e977190e7e9111a2d5a7bd3e5f38d89ca473276a96b26f2c.

CGIG_OSHV1_RNASEQ_PRJNA1329250

Phenotype interpretation: Infection up-regulated antiviral defense transcripts, including RIG-I-like receptors and the interferon-like pathway, similarly in both lineages; AREE records molecular antiviral responsiveness, not an organism-level survival difference

Biological replication: 42 RNA-seq libraries: 12 controls (6 Midori, 6 Miyagi) and 30 OsHV-1 microvariant-challenged (15 Midori, 15 Miyagi, split across USA, Australia, and France viral strains at 5 libraries per lineage-strain cell)

Evidence classification: disease_associated. Tissue: whole_soft_tissue; life stage: spat.

Quality flags: raw_reanalysis;pathogen_challenge_response;lineage_balanced_contrast.

Study limitations:

  • This intake represents the molecular antiviral response to experimental OsHV-1 microvariant challenge, not a measured survival or organism-level disease-resistance difference between the two lineages.
  • Challenge dose, exposure duration, and time of sampling after inoculation are not recorded in the public SRA metadata, so exposure intensity and duration are unresolved.
  • The challenged and control arms are unequal in size (30 versus 12 libraries), but lineage is balanced within each arm and within every viral strain: 6 Midori and 6 Miyagi controls, and 5 of each lineage for each of the Australia, France, and USA strains. The lineage effect is therefore not confounded with infection status, and AREE blocks on lineage as ~ population + condition.
  • The publication contrasts the USA and Australia microvariant strains and reports no marked difference between them; the France-strain libraries are present in the BioProject and are pooled into the challenged arm of the AREE contrast rather than analyzed as a separate strain effect.
  • The BioProject carries no author-deposited processed expression matrix, so AREE has no published feature-level effect table to reconcile against.
  • Condition, lineage, and viral-strain assignments were reconstructed from NCBI BioSample attribute tags rather than from a supplementary sample table; the source metadata spells the Miyagi lineage ‘Myagi’, which AREE normalized.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: 3beb15844b8040ff58c64acd452381f853bbbacb730cce2bf434b6ee05d41f2e.

CGIG_OSHV1_RNASEQ_PRJNA593309

Phenotype interpretation: At 21 C survival 14 days after cohabitation was 52.4%, and OsHV-1 DNA had reached high levels in all three oysters sampled at 48 h. AREE records the whole-oyster transcriptional response to infection, not survival per library.

Biological replication: 6 libraries, one whole 8-month-old oyster per tank: three at 0 h and three at 48 h, the same three 21 C tanks sampled at both times.

Evidence classification: disease_associated. Tissue: whole_organism; life stage: juvenile.

Quality flags: raw_reanalysis;pathogen_challenge_response;unpaired_primary_tank_paired_sensitivity;one_oyster_per_tank;dependence_group=PRJNA593309_21C_tanks.

Study limitations:

  • The 0 h and 48 h oysters come from the same three tanks. The primary model is unpaired (~ condition); a tank-paired model (~ tank + condition) is kept as a sensitivity analysis.
  • One whole oyster per tank was sequenced at each time point, with three replicate tanks per temperature (45 l tanks), so libraries are tank-level replicates and the contrast has three tanks per arm.
  • The two AREE records from this study use different oysters but the same three 21 C tanks (the 21 C 12 h oysters of CGIG_HEAT_RNASEQ_PRJNA593309 and the 21 C 0 h and 48 h oysters of CGIG_OSHV1_RNASEQ_PRJNA593309), so both carry the dependence group PRJNA593309_21C_tanks and count as one replication group.
  • The BioSamples give a collection date of 2019-05-19, but the publication dates the experiment to May 2016; the experimental groups follow the sample aliases (Cg_<CTR | OsHV>hpc_), which agree with the publication’s design.
  • The study also deposited libraries at 26 C and at 24 and 96 h; the publication does not analyse the 26 C libraries, and AREE uses only the time points of its two contrasts.
  • The publication aligned reads with STAR to the Ensembl oyster genome v9.38 and tested with edgeR; AREE quantifies against GCF_963853765.1/RS_2024_06 with Salmon and DESeq2.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: cca3b26b0a11f8cd29467b1df45d728788247bf18855e27becda42e3e88e27d8.

CGIG_PESTO_F0_RNASEQ_PRJEB58194

Phenotype interpretation: The publication found no difference in embryo-larval abnormalities at 24 hpf, disturbed swimming, larger and heavier exposed larvae at some times, a lower epinephrine-induced metamorphosis rate and higher field survival, and differential expression of developmental genes (Hox orthologues, cell-cycle regulators) in gastrulae. AREE records the transcriptional response averaged over two stages.

Biological replication: 36 libraries: 18 control and 18 exposed, nine per stage (gastrula, D larva), labelled with group digits 1-3 and replicate letters A-F.

Evidence classification: environmental_exposure. Tissue: whole_organism; life stage: embryo_larva.

Quality flags: raw_reanalysis;pollutant_exposure_response;pesticide_mixture;pooled_libraries;stage_blocked_contrast;direct_exposure;deseq2_significance_unavailable.

Study limitations:

  • One control D-larva sample title (BT2Ed) is attached to two runs (ERR10659105, ERR10659106); both are kept with the run accession appended to the replicate, and the expressed-SNP check tests whether they are one sample.
  • The exposure is a mixture of 18 pesticides and metals at environmentally relevant concentration (nominal sum 2.85 ug/L; measured 3.74 ug/L in F0) during embryo-larval development (0-48 h post fertilization) at the Ifremer Bouin hatchery.
  • Each library is presumably a pool of many embryos or larvae; pool size, culture vessels and the meaning of the F0 group digits (1-3) are not recorded in the public metadata, and pools were checked with the expressed-SNP test (scripts/check_replicate_genotypes.py) before analysis.
  • All runs come from NCBI .sra archives checked against NCBI’s MD5s (11 F0 runs have no ENA FASTQ), with read counts checked against the NCBI spot counts on conversion.
  • The study also deposited F0 pediveligers (sampled after exposure) and F1 offspring exposed in both generations (EE); AREE does not use them.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: 5eb3fa1130b516be240eecebcd8c88a5b405a6739f99d273764f3054814af5bb.

CGIG_SALINITY_RNASEQ_PRJNA756710

Phenotype interpretation: Higher mortality after Vibrio infection is unfavorable; infected oysters reached 100% mortality at 10 ppt versus 33% at 30 ppt over 12 days (log-rank p < 0.01, n = 30 per tank in triplicate). AREE records the gill transcriptional response to low salinity during infection, not survival per library.

Biological replication: 12 RNA-seq libraries used for the contrast: 6 at 30 ppt and 6 at 10 ppt, three per salinity at each of 12 h and 48 h. Each library pools equal gill tissue from three oysters taken from one tank. Two further 0 h libraries (SRR15559848, SRR15559847) are shared pre-transfer baselines with no salinity contrast and are not used.

Evidence classification: disease_associated. Tissue: gill; life stage: not_recorded.

Quality flags: raw_reanalysis;low_salinity_response;infection_background_all_arms;timepoint_blocked_contrast;fdr_lt_0.05.

Study limitations:

  • Every sequenced oyster was infected with Vibrio alginolyticus, so the contrast estimates the effect of low salinity during infection, not of low salinity alone. The uninfected, sham-injected controls held at 30 ppt were used for the survival curve only and were not sequenced. stressor_class is salinity because salinity is the only factor that differs between the contrasted arms.
  • The primary AREE contrast pools 12 h and 48 h and blocks on sampling time (~ timepoint + condition). It estimates a low-salinity effect averaged over the two times, whereas the publication tested each time separately; each time-by-salinity cell has only three libraries.
  • Each salinity arm used three tanks, and a pool of three oysters was taken from each tank at each time, so the 12 h and 48 h libraries with the same replicate number may share a tank. The public metadata do not state that replicate numbers correspond to tanks, so the model does not pair them.
  • Oyster age or size class is not given beyond an average wet weight of 18.2 +/- 2.8 g, so life_stage is not recorded. Temperature during the challenge is not reported.
  • The library names use C23 for the normal-salinity arm, but the publication and BioSample titles give the normal salinity as 30 ppt; C10 matches the 10 ppt arm.
  • Condition, time and replicate come from the NCBI BioSample titles. They agree with the library names (C10/C23 followed by 1-3 for 12 h and 4-6 for 48 h), and the raw read count of every library matches publication Supplementary Table 3.
  • The publication aligned reads with HISAT2 to GCA_902806645.1, counted them with featureCounts and ran DESeq2 per time point (padj < 0.05, |fold change| > 1). It reports 3,038 DEGs at 12 h (1,731 higher, 1,307 lower at 10 ppt) and 3,187 at 48 h (1,656 higher, 1,531 lower) in Supplementary Table 4, keyed by RefSeq LOC identifiers. Only significant genes are listed, so the table cannot be used as a full ranked result.
  • The companion 16S rRNA amplicon data for the digestive gland are in the separate BioProject PRJNA756403 and are outside this record.
  • Against the publication’s per-time DEGs (Supplementary Table 4), 5,206 of 6,225 LOC identifiers map directly to a current GeneID. Of those tested here, 96.0% (12 h) and 95.7% (48 h) change in the same direction in the pooled AREE model, fold changes correlate at r = 0.784 and 0.759, and none is significant in the opposite direction; 44% reach FDR < 0.05 in the pooled model.
  • rRNA genes are over-represented among low-salinity-higher genes (10 observed versus 0.8 expected; adjusted p = 1.2e-10), which points to a difference in residual rRNA between arms rather than a biological response; those loci should not be interpreted.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: 39f1098bfe2737737332f9e0cd5fdcd64159c848ac8f8562e1d8912615e59fe9.

CGIG_STARVATION_RNASEQ_PRJEB86618

Phenotype interpretation: Lower mortality after OsHV-1 mu-Var cohabitation is favorable. Starved spat had lower final mortality than fed spat (F14R 34.8% versus 75.5%; F11N 93.3% versus 100%; log-rank p < 0.05 in both families), and deaths began later (120 h versus 72-96 h post-cohabitation). AREE records the whole-spat transcriptional state produced by starvation before infection, not survival per library.

Biological replication: 12 RNA-seq libraries used for the contrast, all sampled at 0 h post-cohabitation: 6 fed and 6 starved, three per diet in each of the two families. Each library pools the soft tissue of 10 live spat taken from one tank. The experiment used three tanks per diet, and each tank held both families.

Evidence classification: resilience_associated. Tissue: whole_soft_tissue; life stage: spat.

Quality flags: raw_reanalysis;nutritional_limitation_response;pre_infection_baseline;family_blocked_contrast;fdr_lt_0.05;dependence_group=DECICOMP_PRJEB86525_0h_controls.

Study limitations:

  • The contrast uses only the 0 h libraries. They were sampled when the OsHV-1 mu-Var donor oysters were added, so they capture starvation without infection. The other 71 deposited libraries (3-72 h) follow each diet through infection, when the fed spat carry far more virus. A diet effect pooled across that time course would mix starvation with differential infection progression, so they are not used here.
  • Each tank held both families, so the F11N and F14R libraries with the same replicate number probably share a tank. The public metadata do not state that replicate numbers correspond to tanks. The model ~ family + condition treats the six libraries per diet as independent, although only three tanks per diet are independent units, so the error model may be somewhat optimistic.
  • The two families differ strongly in POMS susceptibility and in how much starvation protected them, so a single diet effect averaged over both may hide family-specific responses.
  • The ENA sample titles call every time point ‘hours post-infection’, including 0 h. The publication shows the 0 h pools were sampled at the start of cohabitation, before viral DNA load rose.
  • The publication reports 84 sequenced libraries and 83 are deposited; the starved F14R 6 h replicate 3 is absent. It also excluded one F11N 48 h replicate 2 library with 21% mapping without naming its diet, and both such libraries are deposited. Neither affects the 0 h contrast.
  • The STAR Methods key-resources table swaps the two RNA-seq accessions (PRJEB86618 ‘for the Fed condition’). The data statement, the ENA study titles and every ENA sample title agree that PRJEB86618 is the starved arm and PRJEB86525 the fed arm, and the manifest follows them.
  • Condition, family, time and replicate come from the ENA sample aliases and titles (sam_P_<3 fed | 5 starved>_), and the two sources agree for every library.
  • The publication aligned reads with STAR to the Penaloza et al. 2021 assembly and analysed the 0 h diet contrast with WGCNA, treating the families as replicates, not with a per-gene differential expression test. Supplementary Table S5 lists the 13,238 genes in the three diet-associated modules (6,307 in the fed-associated turquoise module; 3,603 blue and 3,328 brown starved-associated), with module membership scores and, for 9,540 of them, an NCBI GeneID. It supports a direction-level comparison but gives no effect sizes or per-gene significance.
  • The DECICOMP 30 C fed arm (PRJEB86646) shares the same families, design and fed controls, but the publication does not use it and it is outside this record.
  • Against publication Supplementary Table S5, 8,609 genes of the diet-associated WGCNA modules have a GeneID tested here. 98.3% of fed-module genes are lower and 97.5% of starved-module genes are higher in starved pools (99.0% among genes with |module membership| > 0.8), and only 37 are significant in the opposite direction.
  • Nearly half the tested genes are significant because the design has low within-cell variance (pools of 10 full sibs, BCV 9.8%) and 16 days of starvation is a large perturbation. The replicate genotype check rules out resequenced RNA. Pooling still removes individual variation, so the error model describes pool-to-pool variation, not animal-to-animal variation.
  • Starved pools map about two points less to the transcriptome and more to the genome decoy than fed pools. The shift is small beside the expression effects but systematic, and it may reflect more unspliced or intronic RNA under starvation.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: 7f2cdd37eb935232ae95a393eae8b8fd6fc09a04313c6127ab75f9d0422d0d57.

CGIG_TIRELEACHATE_RNASEQ_PRJNA856813

Phenotype interpretation: The publication reports activation of energy metabolism and stress responses under leachate exposure, signs of microbiota dysbiosis at high concentrations, reduced sperm motility and oocyte transcriptional changes, but no significant ecophysiological change in adults and no effect on offspring growth, which it interprets as possible resilience. AREE records the gill transcriptional response, not a measured phenotype per library.

Biological replication: 12 day-41 gill libraries, one oyster each: 6 controls and 6 high-concentration leachate oysters.

Evidence classification: stress_response. Tissue: gill; life stage: adult.

Quality flags: raw_reanalysis;pollutant_exposure_response;shared_control_group;dependence_group=PRJNA856813_gill_controls;tire_leachate.

Study limitations:

  • The publication is not open access and was not read in full; exposure concentrations, tank layout and the number of tanks per treatment are therefore not recorded here. If each treatment was held in a single tank, oyster-level replicates would not account for tank clustering.
  • Both AREE records from this study test against the same six day-41 gill control oysters, so they carry the dependence group PRJNA856813_gill_controls and count as one replication group.
  • AREE uses only the day-41 gill libraries of the control, high-concentration leachate and high-concentration tire microparticle groups. The deposited study also has low concentrations, a natural-particle control, digestive gland (day 44), gonad (day 44, four per group), time-zero baselines and 16S amplicon data, which are not used.
  • Groups come from the SRA sample names (D41_GI_<CTL | LEA_H | MR_H>_; D41 = day 41, GI = gill); the BioSamples record the tissue as gills. The ENA portal was unavailable when the study was staged, so the reads come from NCBI .sra archives checked against NCBI’s MD5s, and each converted FASTQ must match the SRA spot count.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: 9f52f19bfaa110f7d5ba09e6d53e20073cd2ac0f540b2270c2aef8363909af42.

CGIG_TIREPARTICLE_RNASEQ_PRJNA856813

Phenotype interpretation: The publication reports activation of energy metabolism and stress responses under leachate exposure, signs of microbiota dysbiosis at high concentrations, reduced sperm motility and oocyte transcriptional changes, but no significant ecophysiological change in adults and no effect on offspring growth, which it interprets as possible resilience. AREE records the gill transcriptional response, not a measured phenotype per library.

Biological replication: 12 day-41 gill libraries, one oyster each: 6 controls and 6 high-concentration tire-particle oysters.

Evidence classification: stress_response. Tissue: gill; life stage: adult.

Quality flags: raw_reanalysis;pollutant_exposure_response;shared_control_group;dependence_group=PRJNA856813_gill_controls;tire_microparticles.

Study limitations:

  • The publication is not open access and was not read in full; exposure concentrations, tank layout and the number of tanks per treatment are therefore not recorded here. If each treatment was held in a single tank, oyster-level replicates would not account for tank clustering.
  • Both AREE records from this study test against the same six day-41 gill control oysters, so they carry the dependence group PRJNA856813_gill_controls and count as one replication group.
  • AREE uses only the day-41 gill libraries of the control, high-concentration leachate and high-concentration tire microparticle groups. The deposited study also has low concentrations, a natural-particle control, digestive gland (day 44), gonad (day 44, four per group), time-zero baselines and 16S amplicon data, which are not used.
  • Groups come from the SRA sample names (D41_GI_<CTL | LEA_H | MR_H>_; D41 = day 41, GI = gill); the BioSamples record the tissue as gills. The ENA portal was unavailable when the study was staged, so the reads come from NCBI .sra archives checked against NCBI’s MD5s, and each converted FASTQ must match the SRA spot count.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: a4d40ba861f26ac5d8a0c986666287ba5ad8ba6b4e7fd4766c34c1cf11c4e483.

CGIG_VCRASSOSTREAE_RNASEQ_PRJNA515169

Phenotype interpretation: Virulent V. crassostreae J2-9 killed about 60% of oysters within 3 days, through contact-dependent hemocyte lysis requiring the r5.7 gene. AREE records the 8 h whole-spat transcriptional response to infection, not survival per library.

Biological replication: 6 libraries, each a pool of 30 spat: one mock and one V. crassostreae J2-9 pool from each of three independent infection experiments. Pools of 30 animals average out individual variation, so replication is at the level of experiments.

Evidence classification: disease_associated. Tissue: whole_organism; life stage: spat.

Quality flags: raw_reanalysis;pathogen_challenge_response;pooled_libraries;experiment_blocked_contrast;shared_control_group;dependence_group=PRJNA515169_mock_pools;virulent_vibrio_crassostreae.

Study limitations:

  • The two AREE records from this study (the other is CGIG_VTASMANIENSIS_RNASEQ_PRJNA515169, V. tasmaniensis LGP32) test against the same three mock-injected pools, so they carry the dependence group PRJNA515169_mock_pools and count as one replication group.
  • Replicate numbers in the sample aliases are taken to identify the independent infection experiment, as the publication describes one pool per condition per experiment; the model blocks on experiment (~ experiment + condition).
  • The mock control is an injection of sterile seawater, so the contrast isolates the response to the bacteria, not to injection; the study’s untreated oysters and its nonvirulent strains (V. tasmaniensis LMG20012T, V. crassostreae J2-8) are not used.
  • The publication reports 331 oyster genes modulated by Vibrio relative to mock across strains, analysed in the IHPE Galaxy pipeline; AREE quantifies against GCF_963853765.1/RS_2024_06 with Salmon and DESeq2.
  • ENA served the FASTQs at 0.05-0.35 MB/s in October 2026, so the five runs that had not arrived intact (SRR8551080, SRR8551081, SRR8551083, SRR8551087, SRR8551091: one mock, one J2-9 and all three LGP32 pools) were taken from NCBI .sra archives instead (CGIG_VIBRIO_RNASEQ_PRJNA515169_sra_locations.tsv). Each library comes from a single source.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: c3f3c2871db357f202ace0a460515cf823255d294a67864d4284126428ed56d0.

CGIG_VIBRIO12H_RNASEQ_PRJNA852380

Phenotype interpretation: No phenotype is deposited; AREE records the digestive-gland transcriptional response to Vibrio exposure.

Biological replication: 6 digestive-gland libraries: 3 controls and 3 at 12 h after Vibrio exposure. Each library is a pool by expressed-SNP genotype; pool size is not recorded.

Evidence classification: disease_associated. Tissue: digestive_gland; life stage: juvenile.

Quality flags: raw_reanalysis;pathogen_challenge_response;pooled_libraries;vibrio_species_not_stated;no_linked_publication;shared_control_group;dependence_group=PRJNA852380_controls;vibrio_12h;fdr_lt_0.05.

Study limitations:

  • No publication is linked. The Vibrio species, dose and exposure route are not recorded; the submitting group’s 2022 study (Liu et al., Fish Shellfish Immunol 131:1027-1039) used Vibrio mediterranei, and a 2025 study describes digestive-gland data after V. parahaemolyticus at 12 and 24 h, but neither is confirmed to be this dataset.
  • The two AREE records from this study (12 h and 24 h) test against the same three control libraries, so they carry the dependence group PRJNA852380_controls and count as one replication group.
  • Groups come from the SRA sample names (control1-3, Vh_), which agree with the BioSample replicate attributes.
  • All runs come from NCBI .sra archives checked against NCBI’s MD5s, with read counts checked against the NCBI spot counts on conversion.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: 3484a36bd3fc9e7ca5eb367ded0e5847e78e02de5f9a07d6243b43502e0a2db2.

CGIG_VIBRIO24H_RNASEQ_PRJNA852380

Phenotype interpretation: No phenotype is deposited; AREE records the digestive-gland transcriptional response to Vibrio exposure.

Biological replication: 6 digestive-gland libraries: 3 controls and 3 at 24 h after Vibrio exposure. Each library is a pool by expressed-SNP genotype; pool size is not recorded.

Evidence classification: disease_associated. Tissue: digestive_gland; life stage: juvenile.

Quality flags: raw_reanalysis;pathogen_challenge_response;pooled_libraries;vibrio_species_not_stated;no_linked_publication;shared_control_group;dependence_group=PRJNA852380_controls;vibrio_24h;fdr_lt_0.05.

Study limitations:

  • No publication is linked. The Vibrio species, dose and exposure route are not recorded; the submitting group’s 2022 study (Liu et al., Fish Shellfish Immunol 131:1027-1039) used Vibrio mediterranei, and a 2025 study describes digestive-gland data after V. parahaemolyticus at 12 and 24 h, but neither is confirmed to be this dataset.
  • The two AREE records from this study (12 h and 24 h) test against the same three control libraries, so they carry the dependence group PRJNA852380_controls and count as one replication group.
  • Groups come from the SRA sample names (control1-3, Vh_), which agree with the BioSample replicate attributes.
  • All runs come from NCBI .sra archives checked against NCBI’s MD5s, with read counts checked against the NCBI spot counts on conversion.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: 5571ce917cf19997c36f6355b662989aae8ab9a46b93cdcba8c42dd77b580264.

CGIG_VTASMANIENSIS_RNASEQ_PRJNA515169

Phenotype interpretation: Virulent V. tasmaniensis LGP32 killed about 60% of oysters within 3 days, through phagocytosis-dependent hemocyte lysis requiring intracellular T6SS effectors. AREE records the 8 h whole-spat transcriptional response to infection, not survival per library.

Biological replication: 6 libraries, each a pool of 30 spat: one mock and one V. tasmaniensis LGP32 pool from each of three independent infection experiments. Pools of 30 animals average out individual variation, so replication is at the level of experiments.

Evidence classification: disease_associated. Tissue: whole_organism; life stage: spat.

Quality flags: raw_reanalysis;pathogen_challenge_response;pooled_libraries;experiment_blocked_contrast;shared_control_group;dependence_group=PRJNA515169_mock_pools;virulent_vibrio_tasmaniensis.

Study limitations:

  • The two AREE records from this study (the other is CGIG_VCRASSOSTREAE_RNASEQ_PRJNA515169, V. crassostreae J2-9) test against the same three mock-injected pools, so they carry the dependence group PRJNA515169_mock_pools and count as one replication group.
  • Replicate numbers in the sample aliases are taken to identify the independent infection experiment, as the publication describes one pool per condition per experiment; the model blocks on experiment (~ experiment + condition).
  • The mock control is an injection of sterile seawater, so the contrast isolates the response to the bacteria, not to injection; the study’s untreated oysters and its nonvirulent strains (V. tasmaniensis LMG20012T, V. crassostreae J2-8) are not used.
  • The publication reports 331 oyster genes modulated by Vibrio relative to mock across strains, analysed in the IHPE Galaxy pipeline; AREE quantifies against GCF_963853765.1/RS_2024_06 with Salmon and DESeq2.
  • ENA served the FASTQs at 0.05-0.35 MB/s in October 2026, so the five runs that had not arrived intact (SRR8551080, SRR8551081, SRR8551083, SRR8551087, SRR8551091: one mock, one J2-9 and all three LGP32 pools) were taken from NCBI .sra archives instead (CGIG_VIBRIO_RNASEQ_PRJNA515169_sra_locations.tsv). Each library comes from a single source.

Provenance: Reference GCF_963853765.1 / RS_2024_06; mapping exact; workflow 0.1.0.

Harmonized evidence. Processed input. Input SHA-256: 5d73a1e8fde80d460f9fd8139da0c2ff4bf0c58c1a2f0d5e1d0afadc8147c53d.

Next validation step

Independent biological replication with matched exposure and measured phenotype, an appropriate experimental-unit model, and targeted validation in the relevant tissue and life stage.