14:I[71133,["4507","static/chunks/e2799680-446c489506865058.js","7665","static/chunks/7665-81c5407bc4489ff9.js","8500","static/chunks/8500-e4dd1cb1bb626827.js","9753","static/chunks/9753-028c0a26aa71f2d4.js","8437","static/chunks/8437-a09d0b83fcec76af.js","4557","static/chunks/4557-64c5efa221bc528f.js","947","static/chunks/947-410627998426400f.js","5357","static/chunks/5357-6142894e88cb7da2.js","8382","static/chunks/app/dataset/%5Bid%5D/data/view/page-177ea8adea7aed86.js"],"default"] 15:T53c,

This dataset is one of many generated from an iron incubation experiment conducted as part of the PUPCYCLE II cruise in May and June 2023. The different analyses are listed here, with links to other datasets from this study in the Related Datasets section.

Analyses include:

Physiology (dataset 994607)
- Dissolved inorganic nutrients (nitrate, phosphate, silicic acid)
- Size-fractionated chlorophyll a
- Size-fractionated isotope uptake rates (nitrate and DIC)
- Flow cytometry
- FlowCAM analysis
- FLP experiments for mixotroph cell-specific grazing rates
RNA taxonomy data for normalized reads (dataset 1003525)
Differential expression for +Fe versus DFB treatments at T1 and T2: Pfam
- Mixotroph Pfam (this dataset 1004659)
- Diatom Pfam (dataset 1004674)
Differential expression for +Fe versus DFB treatments at T1 and T2 for mixotrophs and diatoms: KEGG gene (dataset 1004689)
Differential expression for +Fe versus DFB treatments at T1 and T2 for mixotrophs and diatoms: KEGG pathways (this dataset 1004704)
Dissolved iron concentrations (dataset 1004007)

See Related Datasets section below for links to above mentioned datasets.

16:T12c7,

Data collection took place on the R/V Sally Ride from May 29th to June 10th, 2023. To simulate upwelling conditions under different iron treatments, an onboard incubation experiment was conducted. Seawater for the incubation experiment was collected within the northern CCS, off the southern coast of Oregon at 43°02'42.7"N, 124°33'07.2"W. The collected seawater was deemed as freshly upwelled water and taken from a depth of 55 meters – which corresponds to a depth slightly below the euphotic zone receiving less than 1% irradiance – using trace-metal clean techniques on May 30th, 2023, 13:30 GMT. The seawater was pumped and homogenized using trace metal clean techniques, then transferred into a total of thirty 20 L low-density polyethylene cubitainers. Three cubitainers were immediately harvested for the initial timepoint (T0). The remaining twenty-seven cubitainers were assigned treatments: nine were unamended (Ctrl), nine were amended with 5 nM FeCl₂(+Fe), and nine were amended with 200 nM desferrioxamine B, a strong iron chelator that inhibits dissolved iron uptake (DFB).The cubitainers were placed in an on-deck incubator covered with two layers of neutral density screening to simulate 26% of incident irradiance. The incubator was supplied with flow-through surface seawater to maintain ambient surface temperature. Three cubitainers from each treatment were harvested for each of the three subsequent timepoints: 48 hours (T1), 168 hours (T2), and 264 hours (T3) after incubation.

RNA Collection and Metatranscriptomics. Cubitainers were harvested for RNA by collecting and filtering approximately 2.5 to 4 L of each sample onto 0.8 μm Pall Supor filters (142 mm) using a peristaltic pump; samples were then immediately flash frozen in liquid nitrogen and stored at -80°C. RNA was extracted using the RNAqueous-4PCR kit, per instructions from the manufacturer with a few modifications – filters were cut up due to their large size, 200 μL of glass beads were added, and 3 mL of lysis buffer was added. RNA samples were sent to GENEWIZ for library preparation and sequencing with PolyA tail selection. Sequencing was performed on an Illumina HiSeq 4000 with a 2x150 bp configuration. GENEWIZ provided raw paired-end read sequences for each sample.

Raw reads were trimmed using Trim Galore v0.6.10 (Krueger, 2015), which uses Cutadapt (Martin, 2011), and quality control was assessed through FastQC (Andrews, 2010). A de novo metatranscriptome assembly was conducted using rnaSPAdes v3.15.5 for individual assemblies (Bushmanova et al., 2019) and CD-HIT v4.8.1 for a grand assembly (Li & Godzik, 2006). TransDecoder v5.7.1 was used to determine the most likely open reading frame (ORF) for each contig (Haas, 2023). ORFs were annotated using the Marine Functional Eukaryotic Reference Taxa (MarFERReT) database v1.1.1 (e-value < 1e^-06) (Groussman et al., 2023) for NCBI taxonomies (Federhen, 2011), PR² taxonomies (Guillou et al., 2012) and Pfam 34.0 functions (Mistry et al., 2021). To remove any potential bacterial contamination, ORFs were also annotated against the PhyloDB v1076 database using EUKulele (e-value < 1e^-06) (Krinos et al., 2020). Any ORF with a higher bitscore for a bacteria or virus compared to its MarFERReT annotation was removed (resulting in < 3% of ORFs removed). To maximize functional annotation coverage, ORFs were annotated using eggNOG-mapper v2.1.12 (Cantalapiedra et al., 2021) to obtain annotations from the Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa, 2000). Trimmed samples were then aligned to the grand assembly using Salmon v1.10.3 (Patro et al., 2017) and exported into a comprehensive counts table using tximport (Soneson et al., 2016). Reads not mapping to a protist group as defined by MarFERReT annotations were removed before downstream analysis.

DESeq2 was used to obtain normalized counts and conduct differential expression analyses (Love et al., 2014). Differential expression analysis for Pfam protein domains was performed for reads mapping to mixotrophs. Reads not mapping to mixotrophs were removed prior to normalization and statistical testing. Mixotrophs were identified based on constitutive mixotroph genera within the Mixoplankton Database (MDB) (Mitra et al., 2023). Reads were summed based on annotation ID prior to entering DESeq2 analysis; if a contig annotated to multiple IDs, its reads were evenly distributed between the mapped annotations. The differential expression analyses only focused on comparing +Fe vs DFB at T1 and T2. Normalized counts for each Pfam protein domain were also obtained to determine relative contributions of mixotroph genera to the annotation ID.

17:T7f5,

The metatranscriptomic pipeline code used for this project is publicly available at https://github.com/emilyspeciale/Speciale-Metatranscriptomics.

Software: Trim Galore v0.6.10; FastQC v0.12.1; rnaSPAdes v3.15.5; CD-HIT v4.8.1; TransDecoder v5.7.1; EUKulele v2.0.9; eggNOG-mapper v2.1.12; Salmon v1.10.3; tximport v1.3.9; DESeq2 v1.42.1.

Databases/reference resources: MarFERReT v1.1.1; NCBI Taxonomy; PR2; Pfam 34.0; PhyloDB v1076; Kyoto Encyclopedia of Genes and Genomes (KEGG); Mixoplankton Database (MDB).

1. Raw reads were trimmed with Trim Galore v0.6.10 (Krueger, 2015)
2. Read quality was assessed with FastQC v0.12.1 (Andrews, 2010)
3. De novo metatranscriptome assemblies were generated using rnaSPAdes v3.15.5 (Bushmanova et al., 2019) for individual assemblies and CD-HIT v4.8.1 (Li & Godzik, 2006) for the grand assembly
4. ORFs were identified with TransDecoder v5.7.1 (Haas, 2023).
5. Taxonomy and functional annotation via Marine Functional Eukaryotic Reference Taxa (MarFERReT) database v1.1.1 (Groussman et al., 2023); includes NCBI taxonomies (Federhen, 2011), PR2 taxonomies (Guillou et al., 2012), and Pfam 34.0 functions (Mistry et al., 2021).
6. PhyloDB v1076 filtering using EUKulele v2.0.9 (Krinos et al., 2020).
7. eggNOG-mapper v2.1.12 (Cantalapiedra et al., 2021) / KEGG annotation (Kanehisa, 2000)
8. Reads were quantified with Salmon v1.10.3 (Patro et al., 2017)
9. Reads were imported with tximport v1.3.9 (Soneson et al., 2016)
10. Taxonomic filtering where reads outside the target taxonomic group were removed
11. Annotation-ID read summation
12. DESeq2 v1.42.1 (Love et al., 2014) was used to compare +Fe versus DFB at T1 and T2

18:T605,This dataset includes RNA-seq differential expression results of Pfam domains within mixotrophs from an iron incubation experiment using upwelled waters sampled in the California Current System (CCS) during the PUPCYCLE II cruise with Chief Scientist Adrian Marchetti. PUPCYCLE II (Phytoplankton response to the UPwelling CYCLE) took place onboard the R/V Sally Ride from May 29th to June 10th, 2023. Freshly upwelled seawater for the incubation experiment was collected within the northern CCS and placed into designated cubitainers. Three cubitainers were immediately harvested for the initial timepoint (T0). The remaining twenty-seven cubitainers were assigned treatments: nine were unamended (Ctrl), nine were amended with 5 nM FeCl2 (+Fe), and nine were amended with 200 nM Desferrioxamine B, a strong iron chelator (DFB). Three cubitainers from each treatment were harvested for each of the three subsequent timepoints: 48 hours (T1), 168 hours (T2), and 264 hours (T3) after incubation. Approximately 2.5 to 4 liters of each sample was harvested for RNA-seq of the natural phytoplankton community, or metatranscriptomics. The dataset provides differential expression results of Pfam domains within mixotrophs (and which genera they map to), specifically comparing expression levels between +Fe and DFB at T1 and T2. These data include domains that significantly differ in expression between treatments in mixotrophs, and results were used to investigate phagotrophic versus phototrophic gene expression based on iron bioavailability.19:Tc97,NSF Award Abstract: Upwelling zones are hotspots of photosynthesis that are very dynamic in space and time. Microscopic algae, known as phytoplankton, bloom when deep, nutrient-rich waters are upwelled into sunlit surface layers of the ocean, providing nourishment that supports productive food webs and draws down carbon dioxide (CO2) from the atmosphere to the deep ocean. Photosynthetic microbes in these regions must constantly adapt to changes in their chemical and physical environments. For example, subsurface populations respond to changes in light as they approach the surface. When upwelled waters move offshore, cells sink out of the illuminated zone, establishing seed populations that remain inactive until the next upwelling event. This process is called the upwelling conveyor belt cycle (UCBC). How phytoplankton respond to these changes in environmental conditions and how they may influence their nutrient requirements remains unknown. With future ocean changes predicted to alter seawater chemistry, including ocean acidification and decreased iron availability, some phytoplankton groups may be more vulnerable than others. Accompanying educational activities provide learning experiences to enhance understanding and awareness of marine microbes. The development of a research hub at UNC aims to provide infrastructure and support for scientists and students conducting research on environmental genomics. A laboratory component for an upper-level undergraduate course focused on marine phytoplankton is being developed. Educational outreach activities to broader communities include creation of a lesson plan on phytoplankton in upwelling zones and a virtual research cruise experience for middle-school students, as well as a hands-on lab activity for a local museum focused on marine phytoplankton and the important roles they play in shaping our planet. The project examines how phytoplankton respond at the molecular and physiological level to the different UCBC stages, which seed populations (i.e., surface versus subsurface) contribute most to phytoplankton blooms during upwelling events of varying intensity, how phytoplankton elemental compositions are altered throughout UCBC stages, and how future predicted ocean conditions will affect the phytoplankton responses to UCBC conditions. This project contains both laboratory and fieldwork. In the laboratory, phytoplankton isolates recently obtained from upwelling regions are exposed to simulated UCBC conditions to examine changes in gene expression, growth and photosynthetic characteristics and elemental composition. Cultures are subjected to both current and future ocean conditions, including reduced iron availability and higher CO2. In the field, research cruises within upwelling regions study the dynamics of natural phytoplankton communities (both surface and subsurface) experiencing upwelling and relaxation and within simulated upwelling incubation experiments. Knowledge of how phytoplankton are affected by UCBC conditions at an integrated molecular, physiological and elemental level under both current and future scenarios is imperative for the proper conservation and management of these critically important ecosystems.6:["$","$L14",null,{"dataset":{"_id":"1004659","description":"$15","problem":"Analyses focused on +Fe vs. DFB treatments at T1 and T2 to compare an iron-replete versus iron-limited environment; Ctrl and T3 were excluded due to large variability among biological replicates, possibly caused by low read and annotation mappings.","acquisition":"$16","processing":"$17","bcodmoProcessing":"CURATION ACTIONS PERFORMED ON DATA\r\n- Loaded source file \"BCO_DMO_PUPCYCLEII_MixoPfamDE.csv\" into the BCO-DMO system converting missing data identifier NA to empty cell\r\n- Added columns for latitude and longitude of sampling station.\r\n- Set data types for mixotroph columns, coordinates, and differential expression statistics to numeric, and identifier/categorical columns to string.\r\n- Exported final file as \"1004659_v1_mixotroph_pfam.csv\"\r\n\r\nCURATION ACTIONS PERFORMED ON METADATA\r\n- Changed in-line citations to have format of (author, year) replacing number references.\r\n\r\nISSUES POTENTIALLY IMPACTING REUSE\r\n- N/A\r\n\r\n","title":"Pfam differential expression data for mixotrophs from an iron incubation experiment aboard the R/V Sally Ride during the PUPCYCLE II cruise within the California Current System in May-June 2023","nickname":"Mixotroph Pfam differential expression (PUPCYCLE II)","dataUrl":"","currentState":"Final no updates expected","currentStateId":"7","validated":true,"restricted":false,"version":"1","version_date":"2026-09-11T04:00:00.000Z","abstract":"$18","citation":"Speciale, E., Lim, P., Jeong, Y., Cook, C., Cohen, N., McClure, W., Schnetzer, A., Marchetti, A. (2026) Pfam differential expression data for mixotrophs from an iron incubation experiment aboard the R/V Sally Ride during the PUPCYCLE II cruise within the California Current System in May-June 2023. Biological and Chemical Oceanography Data Management Office (BCO-DMO). (Version 1) Version Date 2026-09-11 [if applicable, indicate subset used]. doi:10.26008/1912/bco-dmo.1004659.1 [access date]","citationParts":{"identifier":"doi:10.26008/1912/bco-dmo.1004659.1","versionDate":"2026-09-11","versionId":"1","year":"2026","title":"Pfam differential expression data for mixotrophs from an iron incubation experiment aboard the R/V Sally Ride during the PUPCYCLE II cruise within the California Current System in May-June 2023","authors":"Speciale, E., Lim, P., Jeong, Y., Cook, C., Cohen, N., McClure, W., Schnetzer, A., Marchetti, A.","publisher":"Biological and Chemical Oceanography Data Management Office (BCO-DMO)."},"doi":"10.26008/1912/bco-dmo.1004659.1","release_date":null,"field_dataset_location":"California Current System, off the southern coast of Oregon at 43°02'42.7\"N, 124°33'07.2\"W, depth 55m","temporal_extent":{"start":"2023-05-29T04:00:00.000Z","end":"2023-06-10T04:00:00.000Z"},"temporal":{"resolution":"time:unitDay","type":"time:Interval","start":{"value":"1685332800","label":"2023-05-29"},"end":{"value":"1686369600","label":"2023-06-10"}},"spatial_extent":{"shape":"point","center":{"latitude":43.045194,"longitude":-124.552},"wkt":" POINT (-124.552 43.045194)","north":0,"east":0,"south":0,"west":0},"dataset_type":["Cruise Results","experimental"],"keywords":["biota","oceans"],"data_files":[{"id":"1007693","title":"1004659_v1_mixotroph_pfam.csv","url":"https://datadocs.bco-dmo.org/dataset/1004659/file/nGEJXBohK9G72M/1004659_v1_mixotroph_pfam.csv","showTableUrl":"https://www.bco-dmo.org/dataset/1004659/data/view/1007693","description":"Differential expression results for Pfam domains in mixotrophs from an iron incubation experiment using upwelled waters sampled in the California Current System (CCS) during the PUPCYCLE II cruises ; Primary data file for dataset ID 1004659, version 1","filename":"1004659_v1_mixotroph_pfam.csv","bytesize":2991353,"checksum":{"value":"45d86ccefc8a5ce250ec142ed5947eb8","Algorithm":{"_id":"190","name":"MD5"}},"mimetype":"text/csv","fileType":{"_id":"192","name":"Comma Separated Values (.csv)"},"isErddap":"","date":"","isPublicDuringEmbargo":false,"primaryDataFile":true,"showTable":true,"restricted":false}],"supplemental_files":[],"roles":[{"role":{"_id":"138","name":"Principal Investigator"},"performedBy":{"_id":"527120","name":"Adrian Marchetti"},"performedAtOrganization":{"_id":"20","name":"University of North Carolina at Chapel Hill"}},{"role":{"_id":"130","name":"Scientist"},"performedBy":{"_id":"663696","name":"Natalie Cohen"},"performedAtOrganization":{"_id":"86","name":"University of Georgia"}},{"role":{"_id":"130","name":"Scientist"},"performedBy":{"_id":"51599","name":"Astrid Schnetzer"},"performedAtOrganization":{"_id":"542130","name":"North Carolina State University"}},{"role":{"_id":"128","name":"Student"},"performedBy":{"_id":"994624","name":"Claire Cook"},"performedAtOrganization":{"_id":"86","name":"University of Georgia"}},{"role":{"_id":"128","name":"Student"},"performedBy":{"_id":"994622","name":"Yubeen Jeong"},"performedAtOrganization":{"_id":"20","name":"University of North Carolina at Chapel Hill"}},{"role":{"_id":"128","name":"Student"},"performedBy":{"_id":"994619","name":"Prisca Lim"},"performedAtOrganization":{"_id":"20","name":"University of North Carolina at Chapel Hill"}},{"role":{"_id":"128","name":"Student"},"performedBy":{"_id":"994628","name":"Will McClure"},"performedAtOrganization":{"_id":"542130","name":"North Carolina State University"}},{"role":{"_id":"128","name":"Student"},"performedBy":{"_id":"994615","name":"Emily Speciale"},"performedAtOrganization":{"_id":"20","name":"University of North Carolina at Chapel Hill"}},{"role":{"_id":"145","name":"BCO-DMO Data Manager"},"performedBy":{"_id":"816518","name":"Dana Stuart Gerlach"},"performedAtOrganization":{"_id":"191","name":"Woods Hole Oceanographic Institution"}}],"deployments":[{"_id":"994610","name":"SR2311","synonyms":["PUPCYCLE II cruise"],"dates":{"start":"2023-05-29T00:00:00.000Z","end":"2023-06-10T00:00:00.000Z"},"platform":{"_id":"777407","name":"R/V Sally Ride"},"chiefScientists":[{"id":"527120","name":"Adrian Marchetti","role":"Chief Scientist"}]}],"funding":[{"_id":"768005","name":"OCE-1751805","fundingSource":{"_id":"355","name":"NSF Division of Ocean Sciences"},"url":"https://www.nsf.gov/awardsearch/show-award?AWD_ID=1751805"}],"instruments":[{"_id":"1004673","name":"Illumina HiSeq 4000","description":"Sequencing was performed on an Illumina HiSeq 4000 with a 2x150 bp configuration.","type_id":"649","type_name":"Automated DNA Sequencer","type_desc":"A DNA sequencer is an instrument that determines the order of deoxynucleotides in deoxyribonucleic acid sequences.","type_acronym":"Automated Sequencer"},{"_id":"1004672","name":"Niskin bottle","description":"Niskin bottles were used to collect the seawater used for metatranscriptomics.\n ","type_id":"413","type_name":"Niskin bottle","type_desc":"A Niskin bottle (a next generation water sampler based on the Nansen bottle) is a cylindrical, non-metallic water collection device with stoppers at both ends. The bottles can be attached individually on a hydrowire or deployed in 12, 24, or 36 bottle Rosette systems mounted on a frame and combined with a CTD. Niskin bottles are used to collect discrete water samples for a range of measurements including pigments, nutrients, plankton, etc.","type_acronym":"Niskin bottle"}],"parameters":[{"_id":"1006011","name":"Taxa_Group","format":"","frictionlessType":"string","description":"

Taxa group being targeted for differential expression analysis. All reads not mapping to the targeted taxa group were removed prior to analysis.

","units":{"_id":"","name":"unitless"},"type":{"_id":"994","name":"taxon","description":"

taxonomic group or entity. This may be a family, class, genus, species, etc.; usually this parameter will contain a mixture of taxonomic entities.

"}},{"_id":"1006012","name":"Timepoint","format":"","frictionlessType":"string","description":"

Timepoint identification (T1 or T2)

","units":{"_id":"","name":"unitless"},"type":{"_id":"472026","name":"time_point","description":"

Arbitrary time point of sampling in a sampling sequence (T0, T1,....Tn)

"}},{"_id":"1006013","name":"DE_Comparison","format":"","frictionlessType":"string","description":"

Treatments being targeted for differential expression analysis. This will always be +FevsDFB, representing the iron-replete vs iron-limited treatments.

","units":{"_id":"","name":"unitless"},"type":{"_id":"489066","name":"treatment","description":"

Experimental conditions applied to experimental units. In comparative experiments, members of the complementary group, the control group, receive either no treatment or a standard treatment.

"}},{"_id":"1006014","name":"Pfam_ID","format":"","frictionlessType":"string","description":"

Number ID of Pfam domain.

","units":{"_id":"","name":"unitless"},"type":{"_id":"1971","name":"accession number","description":"

Database identifier assigned by repository and linked to GenBank or other repository.

"}},{"_id":"1006015","name":"Pfam_Name","format":"","frictionlessType":"string","description":"

Short name of Pfam domain.

","units":{"_id":"","name":"unitless"},"type":{"_id":"1005958","name":"Pfam_domain_name","description":"

Domain name from the Pfam (Protein families) database

"}},{"_id":"1006016","name":"Pfam_Description","format":"","frictionlessType":"string","description":"

Long name of Pfam domain.

","units":{"_id":"","name":"unitless"},"type":{"_id":"1005958","name":"Pfam_domain_name","description":"

Domain name from the Pfam (Protein families) database

"}},{"_id":"1006017","name":"baseMean","format":"","frictionlessType":"number","description":"

The average expression level of the Pfam domain across all samples (normalized).

","units":{"_id":"","name":"number of normalized reads"},"type":{"_id":"491108","name":"mean","description":"

A central value of a discrete set of numbers: specifically, the sum of the values divided by the number of values

"}},{"_id":"1006018","name":"log2FoldChange","format":"","frictionlessType":"number","description":"

DESeq2-estimated log2 fold change in Pfam-domain expression for +Fe (iron-amended treatment) relative to DFB (iron-limited treatment); positive values indicate higher expression in +Fe and negative values indicate higher expression in DFB

","units":{"_id":"","name":"unitless"},"type":{"_id":"1005960","name":"fold_change","description":"

Fold change is a number describing how much a quantity changes going from an initial to a final value or from one condition to another condition (e.g., from 10 to 20 is a 2-fold increase, from 10 to 50 is a five-fold increase).

"}},{"_id":"1006019","name":"lfcSE","format":"","frictionlessType":"number","description":"

Standard error of log2FoldChange estimate.

","units":{"_id":"","name":"unitless"},"type":{"_id":"1975","name":"standard error","description":"

Generic standard error taking on the dimensions of the parameter it is tied to.

"}},{"_id":"1006020","name":"stat","format":"","frictionlessType":"number","description":"

Test statistic used in Wald test to determine significance; log2FoldChange divided by lfcSE, which is used to generate a two-tailed pvalue via the Wald Test

","units":{"_id":"","name":"unitless"},"type":{"_id":"1006000","name":"Wald_z_stat","description":"

The Wald statistic is a statistic used during a Wald test, a test of significance of the regression coefficient.

"}},{"_id":"1006021","name":"pvalue","format":"","frictionlessType":"number","description":"

Raw p-value from Wald test.

","units":{"_id":"","name":"unitless"},"type":{"_id":"887009","name":"p-value","description":"

A quantitative confidence value that represents the probability of obtaining a result at least as extreme as that actually obtained, assuming that the actual value was the result of chance alone.

"}},{"_id":"1006022","name":"padj","format":"","frictionlessType":"number","description":"

The False Discovery Rate adjusted p-value using the Benjamini-Hochberg correction.

","units":{"_id":"","name":"unitless"},"type":{"_id":"887009","name":"p-value","description":"

A quantitative confidence value that represents the probability of obtaining a result at least as extreme as that actually obtained, assuming that the actual value was the result of chance alone.

"}},{"_id":"1006023","name":"Gonyaulax","format":"","frictionlessType":"number","description":"

The number of normalized reads that the mixotroph genera Gonyaulax contributes to the Pfam domain.

","units":{"_id":"","name":"number of normalized reads"},"type":{"_id":"1005010","name":"norm_seq_read_count","description":"

Normalized sequence read count; Sequence read count that has been normalized to account for technical differences among samples, such as sequencing depth or library size.