NeutrobodyPlex—monitoring SARS‐CoV‐2 neutralizing immune responses using nanobodies

In light of the COVID‐19 pandemic, there is an ongoing need for diagnostic tools to monitor the immune status of large patient cohorts and the effectiveness of vaccination campaigns. Here, we present 11 unique nanobodies (Nbs) specific for the SARS‐CoV‐2 spike receptor‐binding domain (RBD), of which 8 Nbs potently inhibit the interaction of RBD with angiotensin‐converting enzyme 2 (ACE2) as the major viral docking site. Following detailed epitope mapping and structural analysis, we select two inhibitory Nbs, one of which binds an epitope inside and one of which binds an epitope outside the RBD:ACE2 interface. Based on these, we generate a biparatopic nanobody (bipNb) with viral neutralization efficacy in the picomolar range. Using bipNb as a surrogate, we establish a competitive multiplex binding assay (“NeutrobodyPlex”) for detailed analysis of the presence and performance of neutralizing RBD‐binding antibodies in serum of convalescent or vaccinated patients. We demonstrate that NeutrobodyPlex enables high‐throughput screening and detailed analysis of neutralizing immune responses in infected or vaccinated individuals, to monitor immune status or to guide vaccine design.     

A genotyping‐in‐thousands by sequencing panel to inform invasive deer management using noninvasive fecal and hair samples

Studies in ecology, evolution, and conservation often rely on noninvasive samples, making it challenging to generate large amounts of high‐quality genetic data for many elusive and at‐risk species. We developed and optimized a Genotyping‐in‐Thousands by sequencing (GT‐seq) panel using noninvasive samples to inform the management of invasive Sitka black‐tailed deer (Odocoileus hemionus sitkensis) in Haida Gwaii (Canada). We validated our panel using paired high‐quality tissue and noninvasive fecal and hair samples to simultaneously distinguish individuals, identify sex, and reconstruct kinship among deer sampled across the archipelago, then provided a proof‐of‐concept application using field‐collected feces on SGang Gwaay, an island of high ecological and cultural value. Genotyping success across 244 loci was high (90.3%) and comparable to that of high‐quality tissue samples genotyped using restriction‐site associated DNA sequencing (92.4%), while genotyping discordance between paired high‐quality tissue and noninvasive samples was low (0.50%). The panel will be used to inform future invasive species operations in Haida Gwaii by providing individual and population information to inform management. More broadly, our GT‐seq workflow that includes quality control analyses for targeted SNP selection and a modified protocol may be of wider utility for other studies and systems where noninvasive genetic sampling is employed.     

Modeling population size independent tissue epigenomes by ChIL‐seq with single thin sections

Recent advances in genome‐wide technologies have enabled analyses using small cell numbers of even single cells. However, obtaining tissue epigenomes with cell‐type resolution from large organs and tissues still remains challenging, especially when the available material is limited. Here, we present a ChIL‐based approach for analyzing the diverse cellular dynamics at the tissue level using high‐depth epigenomic data. “ChIL for tissues” allows the analysis of a single tissue section and can reproducibly generate epigenomic profiles from several tissue types, based on the distribution of target epigenomic states, tissue morphology, and number of cells. The proposed method enabled the independent evaluation of changes in cell populations and gene activation in cells from regenerating skeletal muscle tissues, using a statistical model of RNA polymerase II distribution on gene loci. Thus, the integrative analyses performed using ChIL can elucidate in vivo cell‐type dynamics of tissues.     

Endophilin A2 regulates B‐cell endocytosis and is required for germinal center and humoral responses

Antigen‐specific B‐cell responses require endosomal trafficking to regulate antigen uptake and presentation to helper T cells, and to control expression and signaling of immune receptors. However, the molecular composition of B‐cell endosomal trafficking pathways and their specific roles in B‐cell responses have not been systematically investigated. Here, we report high‐throughput identification of genes regulating B‐cell receptor (BCR)‐mediated antigen internalization using genome‐wide functional screens. We show that antigen internalization depends both on constitutive, clathrin‐mediated endocytosis and on antigen‐induced, clathrin‐independent endocytosis mediated by endophilin A2. Although endophilin A2‐mediated endocytosis is dispensable for antigen presentation, it is selectively required for metabolic support of B‐cell proliferation, in part through regulation of iron uptake. Consequently, endophilin A2‐deficient mice show defects in GC B‐cell responses and production of high‐affinity IgG. The requirement for endophilin A2 highlights a unique importance of clathrin‐independent intracellular trafficking in GC B‐cell clonal expansion and antibody responses.     

Body mass and sex,but not breeding condition and season,influence open‐field exploration in the yellow‐necked mouse

Theory predicts that risk taking should be influenced by external (e.g., season) and internal (e.g., breeding condition, sex, and body mass) conditions. We investigated whether these factors are associated with a potentially risky behavior: exploration of a novel environment. We conducted repeated open‐field tests of exploration in a common forest rodent, the yellow‐necked mouse Apodemus flavicollis. Contrary to expectations, the exploration did not vary with the season (spring vs. fall) or the reproductive status of the tested animals. Also unexpectedly, there was an inverted U‐shaped relationship between body mass and exploration: animals with intermediate body mass tended to have the highest exploration tendencies. Males were more exploratory than females. Finally, even after adjusting for the effects of body mass and sex, individuals exhibited consistent, repeatable differences in exploration tendencies (“behavioral types” or “personalities”). The discrepancies between certain broad generalizations and our results suggest that risk taking depends on details of species‐specific biology.     

Method for quick DNA barcode reference library construction

DNA barcoding has become one of the most important techniques in plant species identification. Successful application of this technology is dependent on the availability of reference database of high species coverage. Unfortunately, there are experimental and data processing challenges to construct such a library within a short time. Here, we present our solutions to these challenges. We sequenced six conventional DNA barcode fragments (ITS1, ITS2, matK1, matK2, rbcL1, and rbcL2) of 380 flowering plants on next‐generation sequencing (NGS) platforms (Illumina Hiseq 2500 and Ion Torrent S5) and the Sanger sequencing platform. After comparing the sequencing depths, read lengths, base qualities, and base accuracies, we conclude that Illumina Hiseq2500 PE250 run is suitable for conventional DNA barcoding. We developed a new “Cotu” method to create consensus sequences from NGS reads for longer output sequences and more reliable bases than the other three methods. Step‐by‐step instructions to our method are provided. By using high‐throughput machines (PCR and NGS), labeling PCR, and the Cotu method, it is possible to significantly reduce the cost and labor investments for DNA barcoding. A regional or even global DNA barcoding reference library with high species coverage is likely to be constructed in a few years.     

High‐throughput sequencing outperforms traditional morphological methods in Blue Catfish diet analysis and reveals novel insights into diet ecology

Blue Catfish Ictalurus furcatus are an invasive, yet economically important species in the Chesapeake Bay. However, their impact on the trophic ecology of this system is not well understood. In order to provide in‐depth analysis of predation by Blue Catfish, we identified prey items using high‐throughput DNA sequencing (HTS) of entire gastrointestinal tracts from 134 samples using two genetic markers, mitochondrial cytochrome c oxidase I (COI) and the nuclear 18S ribosomal RNA gene. We compared our HTS results to a more traditional “hybrid” approach that coupled morphological identification with DNA barcoding. The hybrid study was conducted on additional Blue Catfish samples (n = 617 stomachs) collected from the same location and season in the previous year. Taxonomic representation with HTS vastly surpassed that achieved with the hybrid methodology in Blue Catfish. Significantly, our HTS study identified several instances of at‐risk and invasive species consumption not identified using the hybrid method, supporting the hypothesis that previous studies using morphological methods may greatly underestimate consumption of critical species. Finally, we report the novel finding that Blue Catfish diet diversity inversely correlates to daily flow rates, perhaps due to higher mobility and prey‐seeking behaviors exhibited during lower flow.     

NGSpeciesID: DNA barcode and amplicon consensus generation from long‐read sequencing data

Third‐generation sequencing technologies, such as Oxford Nanopore Technologies (ONT) and Pacific Biosciences (PacBio), have gained popularity over the last years. These platforms can generate millions of long‐read sequences. This is not only advantageous for genome sequencing projects, but also advantageous for amplicon‐based high‐throughput sequencing experiments, such as DNA barcoding. However, the relatively high error rates associated with these technologies still pose challenges for generating high‐quality consensus sequences. Here, we present NGSpeciesID, a program which can generate highly accurate consensus sequences from long‐read amplicon sequencing technologies, including ONT and PacBio. The tool includes clustering of the reads to help filter out contaminants or reads with high error rates and employs polishing strategies specific to the appropriate sequencing platform. We show that NGSpeciesID produces consensus sequences with improved usability by minimizing preprocessing and software installation and scalability by enabling rapid processing of hundreds to thousands of samples, while maintaining similar consensus accuracy as current pipelines.     

Predictable fold switching by the SARS‐CoV‐2 protein ORF9b

Extant fold‐switching proteins remodel their secondary structures and change their functions in response to environmental stimuli. These shapeshifting proteins regulate biological processes and are associated with a number of diseases, including tuberculosis, cancer, Alzheimer''s, and autoimmune disorders. Thus, predictive methods are needed to identify more fold‐switching proteins, especially since all naturally occurring instances have been discovered by chance. In response to this need, two high‐throughput predictive methods have recently been developed. Here we test them on ORF9b, a newly discovered fold switcher and potential therapeutic target from the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS‐CoV‐2). Promisingly, both methods correctly indicate that ORF9b switches folds. We then tested the same two methods on ORF9b1, the ORF9b homolog from SARS‐CoV‐1. Again, both methods predict that ORF9b1 switches folds, a finding consistent with experimental binding studies. Together, these results (a) demonstrate that protein fold switching can be predicted using high‐throughput computational approaches and (b) suggest that fold switching might be a general characteristic of ORF9b homologs.     

Prediction of a competing endogenous RNA co‐expression network as a prognostic marker in glioblastoma

Due to its high proliferation capacity and rapid intracranial spread, glioblastoma (GBM) has become one of the least curable malignant cancers. Recently, the competing endogenous RNAs (ceRNAs) hypothesis has become a focus in the researches of molecular biological mechanisms of cancer occurrence and progression. However, there is a lack of correlation studies on GBM, as well as a lack of comprehensive analyses of GBM molecular mechanisms based on high‐throughput sequencing and large‐scale sample sizes. We obtained RNA‐seq data from The Cancer Genome Atlas (TCGA) and Genotype‐Tissue Expression (GTEx) databases. Further, differentially expressed mRNAs were identified from normal brain tissue and GBM tissue. The similarities between the mRNA modules with clinical traits were subjected to weighted correlation network analysis (WGCNA). With the mRNAs from clinical‐related modules, a survival model was constructed by univariate and multivariate Cox proportional hazard regression analyses. Thereafter, we carried out Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Finally, we predicted interactions between lncRNAs, miRNAs and mRNAs by TargetScan, miRDB, miRTarBase and starBase. We identified 2 lncRNAs (NORAD, XIST), 5 miRNAs (hsa‐miR‐3613, hsa‐miR‐371, hsa‐miR‐373, hsa‐miR‐32, hsa‐miR‐92) and 2 mRNAs (LYZ, PIK3AP1) for the construction of a ceRNA network, which might act as a prognostic biomarker of GBM. Combined with previous studies and our enrichment analysis results, we hypothesized that this ceRNA network affects immune activities and tumour microenvironment variations. Our research provides novel aspects to study GBM development and treatment.     

“Protein aggregates” contain RNA and DNA,entrapped by misfolded proteins but largely rescued by slowing translational elongation

All neurodegenerative diseases feature aggregates, which usually contain disease‐specific diagnostic proteins; non‐protein constituents, however, have rarely been explored. Aggregates from SY5Y‐APP_Sw neuroblastoma, a cell model of familial Alzheimer''s disease, were crosslinked and sequences of linked peptides identified. We constructed a normalized “contactome” comprising 11 subnetworks, centered on 24 high‐connectivity hubs. Remarkably, all 24 are nucleic acid‐binding proteins. This led us to isolate and sequence RNA and DNA from Alzheimer''s and control aggregates. RNA fragments were mapped to the human genome by RNA‐seq and DNA by ChIP‐seq. Nearly all aggregate RNA sequences mapped to specific genes, whereas DNA fragments were predominantly intergenic. These nucleic acid mappings are all significantly nonrandom, making an artifactual origin extremely unlikely. RNA (mostly cytoplasmic) exceeded DNA (chiefly nuclear) by twofold to fivefold. RNA fragments recovered from AD tissue were ~1.5‐to 2.5‐fold more abundant than those recovered from control tissue, similar to the increase in protein. Aggregate abundances of specific RNA sequences were strikingly differential between cultured SY5Y‐APP_Sw glioblastoma cells expressing APOE3 vs. APOE4, consistent with APOE4 competition for E‐box/CLEAR motifs. We identified many G‐quadruplex and viral sequences within RNA and DNA of aggregates, suggesting that sequestration of viral genomes may have driven the evolution of disordered nucleic acid‐binding proteins. After RNA‐interference knockdown of the translational‐procession factor EEF2 to suppress translation in SY5Y‐APP_Sw cells, the RNA content of aggregates declined by >90%, while reducing protein content by only 30% and altering DNA content by ≤10%. This implies that cotranslational misfolding of nascent proteins may ensnare polysomes into aggregates, accounting for most of their RNA content.     

High‐altitude adaptation in vertebrates as revealed by mitochondrial genome analyses

The high‐altitude environment may drive vertebrate evolution in a certain way, and vertebrates living in different altitude environments might have different energy requirements. We hypothesized that the high‐altitude environment might impose different influences on vertebrate mitochondrial genomes (mtDNA). We used selection pressure analyses and PIC (phylogenetic independent contrasts) analysis to detect the evolutionary rate of vertebrate mtDNA protein‐coding genes (PCGs) from different altitudes. The results showed that the ratio of nonsynonymous/synonymous substitutions (dN/dS) in the mtDNA PCGs was significantly higher in high‐altitude vertebrates than in low‐altitude vertebrates. The seven rapidly evolving genes were shared by the high‐altitude vertebrates, and only one positive selection gene (ND5 gene) was detected in the high‐altitude vertebrates. Our results suggest the mtDNA evolutionary rate in high‐altitude vertebrates was higher than in low‐altitude vertebrates as their evolution requires more energy in a high‐altitude environment. Our study demonstrates the high‐altitude environment (low atmospheric O₂ levels) drives vertebrate evolution in mtDNA PCGs.     

Scaffold compound L971 exhibits anti‐inflammatory activities through inhibition of JAK/STAT and NFκB signalling pathways

JAK/STAT and NFκB signalling pathways play essential roles in regulating inflammatory responses, which are important pathogenic factors of various serious immune‐related diseases, and function individually or synergistically. To find prodrugs that can treat inflammation, we performed a preliminary high‐throughput screening of 18 840 small molecular compounds and identified scaffold compound L971 which significantly inhibited JAK/STAT and NFκB driven luciferase activities. L971 could inhibit the constitutive and stimuli‐dependent activation of STAT1, STAT3 and IκBα and could significantly down‐regulate the proinflammatory gene expression in mouse peritoneal macrophages stimulated by LPS. Gene expression profiles upon L971 treatment were determined using high‐throughput RNA sequencing, and significant differentially up‐regulated and down‐regulated genes were identified by DESeq analysis. The bioinformatic studies confirmed the anti‐inflammatory effects of L971. Finally, L971 anti‐inflammatory character was further verified in LPS‐induced sepsis shock mouse model in vivo. Taken together, these data indicated that L971 could down‐regulate both JAK/STAT and NFκB signalling activities and has the potential to treat inflammatory diseases such as sepsis shock.     

Oral and long‐acting antipsychotics for relapse prevention in schizophrenia‐spectrum disorders: a network meta‐analysis of 92 randomized trials including 22,645 participants

According to current evidence and guidelines, continued antipsychotic treatment is key for preventing relapse in people with schizophrenia‐spectrum disorders, but evidence‐based recommendations for the choice of the individual antipsychotic for maintenance treatment are lacking. Although oral antipsychotics are often prescribed first line for practical reasons, long‐acting injectable antipsychotics (LAIs) are a valuable resource to tackle adherence issues since the earliest phase of disease. Medline, EMBASE, PsycINFO, CENTRAL and CINAHL databases and online registers were searched to identify randomized controlled trials comparing LAIs or oral antipsychotics head‐to‐head or against placebo, published until June 2021. Relative risks and standardized mean differences were pooled using random‐effects pairwise and network meta‐analysis. The primary outcomes were relapse and dropout due to adverse events. We used the Cochrane Risk of Bias tool to assess study quality, and the CINeMA approach to assess the confidence of pooled estimates. Of 100 eligible trials, 92 (N=22,645) provided usable data for meta‐analyses. Regarding relapse prevention, the vast majority of the 31 included treatments outperformed placebo. Compared to placebo, “high” confidence in the results was found for (in descending order of effect magnitude) amisulpride‐oral (OS), olanzapine‐OS, aripiprazole‐LAI, olanzapine‐LAI, aripiprazole‐OS, paliperidone‐OS, and ziprasidone‐OS. “Moderate” confidence in the results was found for paliperidone‐LAI 1‐monthly, iloperidone‐OS, fluphenazine‐OS, brexpiprazole‐OS, paliperidone‐LAI 1‐monthly, asenapine‐OS, haloperidol‐OS, quetiapine‐OS, cariprazine‐OS, and lurasidone‐OS. Regarding tolerability, none of the antipsychotics was significantly worse than placebo, but confidence was poor, with only aripiprazole (both LAI and OS) showing “moderate” confidence levels. Based on these findings, olanzapine, aripiprazole and paliperidone are the best choices for the maintenance treatment of schizophrenia‐spectrum disorders, considering that both LAI and oral formulations of these antipsychotics are among the best‐performing treatments and have the highest confidence of evidence for relapse prevention. This finding is of particular relevance for low‐ and middle‐income countries and constrained‐resource settings, where few medications may be selected. Results from this network meta‐analysis can inform clinical guidelines and national and international drug regulation policies.     

High‐throughput metabolomics predicts drug–target relationships for eukaryotic proteins

Chemical probes are important tools for understanding biological systems. However, because of the huge combinatorial space of targets and potential compounds, traditional chemical screens cannot be applied systematically to find probes for all possible druggable targets. Here, we demonstrate a novel concept for overcoming this challenge by leveraging high‐throughput metabolomics and overexpression to predict drug–target interactions. The metabolome profiles of yeast treated with 1,280 compounds from a chemical library were collected and compared with those of inducible yeast membrane protein overexpression strains. By matching metabolome profiles, we predicted which small molecules targeted which signaling systems and recovered known interactions. Drug–target predictions were generated across the 86 genes studied, including for difficult to study membrane proteins. A subset of those predictions were tested and validated, including the novel targeting of GPR1 signaling by ibuprofen. These results demonstrate the feasibility of predicting drug–target relationships for eukaryotic proteins using high‐throughput metabolomics.     

Benefits and limitations of a new genome‐based PCR‐RFLP genotyping assay (GB‐RFLP): A SNP‐based detection method for identification of species in extremely young adaptive radiations

High‐throughput DNA sequencing technologies make it possible now to sequence entire genomes relatively easily. Complete genomic information obtained by whole‐genome resequencing (WGS) can aid in identifying and delineating species even if they are extremely young, cryptic, or morphologically difficult to discern and closely related. Yet, for taxonomic or conservation biology purposes, WGS can remain cost‐prohibitive, too time‐consuming, and often constitute a “data overkill.” Rapid and reliable identification of species (and populations) that is also cost‐effective is made possible by species‐specific markers that can be discovered by WGS. Based on WGS data, we designed a PCR restriction fragment length polymorphism (PCR‐RFLP) assay for 19 Neotropical Midas cichlid populations (Amphilophus cf. citrinellus), that includes all 13 described species of this species complex. Our work illustrates that identification of species and populations (i.e., fish from different lakes) can be greatly improved by designing genetic markers using available “high resolution” genomic information. Yet, our work also shows that even in the best‐case scenario, when whole‐genome resequencing information is available, unequivocal assignments remain challenging when species or populations diverged very recently, or gene flow persists. In summary, we provide a comprehensive workflow on how to design RFPL markers based on genome resequencing data, how to test and evaluate their reliability, and discuss the benefits and pitfalls of our approach.     

Degrees of compositional shift in tree communities vary along a gradient of temperature change rates over one decade: Application of an individual‐based temporal beta‐diversity concept

Temporal patterns in communities have gained widespread attention recently, to the extent that temporal changes in community composition are now termed “temporal beta‐diversity.” Previous studies of beta‐diversity have made use of two classes of dissimilarity indices: incidence‐based (e.g., Sørensen and Jaccard dissimilarity) and abundance‐based (e.g., Bray–Curtis and Ružička dissimilarity). However, in the context of temporal beta‐diversity, the persistence of identical individuals and turnover among other individuals within the same species over time have not been considered, despite the fact that both will affect compositional changes in communities. To address this issue, I propose new index concepts for beta‐diversity and the relative speed of compositional shifts in relation to individual turnover based on individual identity information. Individual‐based beta‐diversity indices are novel dissimilarity indices that consider individual identity information to quantitatively evaluate temporal change in individual turnover and community composition. I applied these new indices to individually tracked tree monitoring data in deciduous and evergreen broad‐leaved forests across the Japanese archipelago with the objective of quantifying the effect of climate change trends (i.e., rates of change in both annual mean temperature and annual precipitation) on individual turnover and compositional shifts at each site. A new index explored the relative contributions of mortality and recruitment processes to temporal changes in community composition. Clear patterns emerged showing that an increase in the temperature change rate facilitated the relative contribution of mortality components. The relative speed of compositional shift increased with increasing temperature change rates in deciduous forests but decreased with increasing warming rates in evergreen forests. These new concepts provide a way to identify novel and high‐resolution temporal patterns in communities.     

The formation of “mega‐flocks” depends on vegetation structure in montane coniferous forests of Taiwan

A mixed‐species bird flock is a social assemblage where two or more bird species are moving together while foraging and might benefit from increased foraging efficiency and antipredator vigilance. A “mega‐flock,” which includes flocking species from different vegetation strata, often exhibits high species diversity. Mechanisms for the formation of mega‐flocks have not yet been explored. In this study, we evaluated the influence of vegetation structure and bird species diversity in driving the occurrence of mega‐flocks. We investigated the composition of mixed‐species flocks, local bird communities, and vegetation structure in five vegetation types of two high‐elevation sites in central Taiwan. Mega‐flocks occurred more frequently in pine woodland than later successional stages of coniferous forests. However, species richness/diversity of local bird communities increased along successional stages. Therefore, vegetation variables exhibit more influence on the occurrence of mega‐flocks than local bird communities. Besides foliage height diversity, understory coverage also showed positive effects on flock size of mixed‐species flocks. Our results indicated that pine woodlands with more evenly distributed vegetation layers could facilitate the interactions of canopy and understory flocks and increase the formation of mega‐flocks and thus the complexity of mixed‐species flocks.     

Ancestral ecological regime shapes reaction to food limitation in the Least Killifish,Heterandria formosa

Populations with different densities often show genetically based differences in life histories. The divergent life histories could be driven by several agents of selection, one of which is variation in per‐capita food levels. Its relationship with population density is complex, as it depends on overall food availability, individual metabolic demand, and food‐independent factors potentially affecting density, such as predation intensity. Here, we present a case study of two populations of a small live‐bearing freshwater fish, one characterized by high density, low predation risk, low overall food availability, and presumably low per‐capita food levels, and the other by low density, high predation risk, high overall food availability, and presumably high per‐capita food levels. Using a laboratory experiment, we examined whether fish from these populations respond differently to food limitation, and whether size at birth, a key trait with respect to density variation in this species, is associated with any such differential responses. While at the lower food level growth was slower, body size smaller, maturation delayed, and survival reduced in both populations, these fitness costs were smaller in fish from the high‐density population. At low food, only 15% of high‐density fish died, compared to 75% of low‐density fish. This difference was much smaller at high food (0% vs. 15% mortality). The increased survival of high‐density fish may, at least partly, be due to their larger size at birth. Moreover, being larger at birth enabled fish to mature relatively early even at the lower food level. We demonstrate that sensitivities to food limitation differ between study populations, consistent with selection for a greater ability to tolerate low per‐capita food availability in the high‐density population. While we cannot preclude other agents of selection from operating in these populations simultaneously, our results suggest that variation in per‐capita food levels is one of those agents.     

Mass-transfer based modeling to investigate iodine staining effects for enhanced contrast X-ray computed tomography

Iodine staining combined with X-ray computed tomography (CT) has become a core approach in anatomy, offering three-dimensional and essentially non-destructive imaging of soft tissues. Although there have been rapid advances in methodologies and techniques, the mechanisms underlying diffusible iodine contrast-enhanced CT are not yet fully understood. The protocols for staining samples of differing sizes and tissue types have not yet been justified theoretically. Here we utilize mass transfer modeling to simulate iodine diffusion and predict iodine concentrations within distinct tissue types. We also undertake iodine staining experiments to visualize the detailed anatomy and contrast effects on whole-body avian specimens using different concentrations of iodine solution to compare with model simulation results. The simulations effectively explain most observed concentration changes in differently-sized samples over distinct iodine treatment durations. These results also provide insight into the mechanisms behind the efficacy of solution replenishment for enhancing staining effects. Both consistencies and inconsistencies between our simulation and experimental results regarding iodine concentration in tissues will inform further investigations to optimize iodine staining protocols.