Pub Date : 2025-07-26DOI: 10.1101/2023.06.06.543983
Gabriel Birzu, Harihara Subrahmaniam Muralidharan, Danielle Goudeau, Rex R Malmstrom, Daniel S Fisher, Devaki Bhaya
Bacterial species often undergo rampant recombination yet maintain cohesive genomic identity. Ecological differences can generate recombination barriers between species and sustain genomic clusters in the short term. But can these forces prevent genomic mixing during long-term coevolution? Cyanobacteria in Yellowstone hot springs comprise several diverse species that have coevolved for hundreds of thousands of years, providing a rare natural experiment. By analyzing more than 300 single-cell genomes, we show that despite each species forming a distinct genomic cluster, much of the diversity within species is the result of hybridization driven by selection, which has mixed their ancestral genotypes. This widespread mixing is contrary to the prevailing view that ecological barriers can maintain cohesive bacterial species and highlights the importance of hybridization as a source of genomic diversity.
{"title":"Hybridization breaks species barriers in long-term coevolution of a cyanobacterial population.","authors":"Gabriel Birzu, Harihara Subrahmaniam Muralidharan, Danielle Goudeau, Rex R Malmstrom, Daniel S Fisher, Devaki Bhaya","doi":"10.1101/2023.06.06.543983","DOIUrl":"10.1101/2023.06.06.543983","url":null,"abstract":"<p><p>Bacterial species often undergo rampant recombination yet maintain cohesive genomic identity. Ecological differences can generate recombination barriers between species and sustain genomic clusters in the short term. But can these forces prevent genomic mixing during long-term coevolution? Cyanobacteria in Yellowstone hot springs comprise several diverse species that have coevolved for hundreds of thousands of years, providing a rare natural experiment. By analyzing more than 300 single-cell genomes, we show that despite each species forming a distinct genomic cluster, much of the diversity within species is the result of hybridization driven by selection, which has mixed their ancestral genotypes. This widespread mixing is contrary to the prevailing view that ecological barriers can maintain cohesive bacterial species and highlights the importance of hybridization as a source of genomic diversity.</p>","PeriodicalId":72407,"journal":{"name":"bioRxiv : the preprint server for biology","volume":" ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2025-07-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://ftp.ncbi.nlm.nih.gov/pub/pmc/oa_pdf/24/fb/nihpp-2023.06.06.543983v1.PMC10274767.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"10302167","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-07-26DOI: 10.1101/2023.09.19.558262
Lutz Menzel, Hengbo Zhou, James W Baish, Meghan J O'Melia, Laurel B Darragh, Derek N Effiom, Emma Specht, Juliane Czapla, Pin-Ji Lei, Johanna J Rajotte, Lingshan Liu, Mohammad R Nikmaneshi, Mohammad S Razavi, Matthew G Vander Heiden, Jessalyn M Ubellacker, Lance L Munn, Sana D Karam, Genevieve M Boland, Sonia Cohen, Timothy P Padera
Lymphatic transport facilitates the presentation of cancer antigens in tumor-draining lymph nodes (tdLNs), leading to T cell activation and the generation of systemic anti-cancer immune surveillance. Surgical removal of tdLNs to control cancer progression is routine in clinical practice. However, whether removing tdLNs impairs immune checkpoint blockade (ICB) is still controversial. Our analysis demonstrates that melanoma patients remain responsive to PD-1 checkpoint blockade after regional LN dissection. We were able to recapitulate the persistent response to ICB after regional LN resection in murine melanoma and mammary carcinoma models. Mechanistically, soluble antigen is diverted to distant LNs after tdLN dissection. Consistently, robust ICB responses in patients with head and neck cancer after primary tumor and tdLN resection correlated with the presence of reactive LNs in distant sites. These findings indicate that distant LNs sufficiently compensate for the removal of direct tdLNs and sustain the response to ICB.
{"title":"Distant lymph nodes compensate for resected tumor-draining lymph nodes during cancer immunotherapy.","authors":"Lutz Menzel, Hengbo Zhou, James W Baish, Meghan J O'Melia, Laurel B Darragh, Derek N Effiom, Emma Specht, Juliane Czapla, Pin-Ji Lei, Johanna J Rajotte, Lingshan Liu, Mohammad R Nikmaneshi, Mohammad S Razavi, Matthew G Vander Heiden, Jessalyn M Ubellacker, Lance L Munn, Sana D Karam, Genevieve M Boland, Sonia Cohen, Timothy P Padera","doi":"10.1101/2023.09.19.558262","DOIUrl":"10.1101/2023.09.19.558262","url":null,"abstract":"<p><p>Lymphatic transport facilitates the presentation of cancer antigens in tumor-draining lymph nodes (tdLNs), leading to T cell activation and the generation of systemic anti-cancer immune surveillance. Surgical removal of tdLNs to control cancer progression is routine in clinical practice. However, whether removing tdLNs impairs immune checkpoint blockade (ICB) is still controversial. Our analysis demonstrates that melanoma patients remain responsive to PD-1 checkpoint blockade after regional LN dissection. We were able to recapitulate the persistent response to ICB after regional LN resection in murine melanoma and mammary carcinoma models. Mechanistically, soluble antigen is diverted to distant LNs after tdLN dissection. Consistently, robust ICB responses in patients with head and neck cancer after primary tumor and tdLN resection correlated with the presence of reactive LNs in distant sites. These findings indicate that distant LNs sufficiently compensate for the removal of direct tdLNs and sustain the response to ICB.</p>","PeriodicalId":72407,"journal":{"name":"bioRxiv : the preprint server for biology","volume":" ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2025-07-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10541098/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41179676","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-07-24DOI: 10.1101/2021.12.29.474455
Duygu Payzin-Dogru, Tim Froitzheim, Steven J Blair, Siddhartha G Jena, Hani Singer, Julia C Paoli, Ryan T Kim, Emil Kriukov, Sarah E Wilson, Renzhi Hou, Aaron M Savage, Victor Cat, Louis V Cammarata, S Y Celeste Wu, Vivien Bothe, Burcu Erdogan, Shifa Hossain, Noah Lopez, Julia Losner, Juan Velazquez Matos, Sangwon Min, Sebastian Böhm, Anthony E Striker, Kelly E Dooling, Adam H Freedman, Bobby Groves, Benjamin Tajer, Glory Kalu, Eric Wynn, Alan Y L Wong, Nadia Fröbisch, Petr Baranov, Maksim V Plikus, Jason D Buenrostro, Brian J Haas, Isaac M Chiu, Timothy B Sackton, Jessica L Whited
Many species regenerate lost body parts following amputation. Most limb regeneration research has focused on the immediate injury site. Meanwhile, body-wide injury responses remain largely unexplored but may be critical for regeneration. Here, we discovered a role for the sympathetic nervous system in stimulating a body-wide stem cell activation response to amputation that drives enhanced limb regeneration in axolotls. This response is mediated by adrenergic signaling, which coordinates distant cellular activation responses via the α2A-adrenergic receptor, and local regeneration responses via β-adrenergic receptors. Both α2A- and β-adrenergic signaling act upstream of mTOR signaling. Notably, systemically-activated axolotls regenerate limbs faster than naïve animals, suggesting a potential selective advantage in environments where injury from cannibalism or predation is common. This work challenges the predominant view that cellular responses underlying regeneration are confined to the injury site and argues instead for body-wide cellular priming as a foundational step that enables localized tissue regrowth.
{"title":"Adrenergic signaling coordinates distant and local responses to amputation in axolotl.","authors":"Duygu Payzin-Dogru, Tim Froitzheim, Steven J Blair, Siddhartha G Jena, Hani Singer, Julia C Paoli, Ryan T Kim, Emil Kriukov, Sarah E Wilson, Renzhi Hou, Aaron M Savage, Victor Cat, Louis V Cammarata, S Y Celeste Wu, Vivien Bothe, Burcu Erdogan, Shifa Hossain, Noah Lopez, Julia Losner, Juan Velazquez Matos, Sangwon Min, Sebastian Böhm, Anthony E Striker, Kelly E Dooling, Adam H Freedman, Bobby Groves, Benjamin Tajer, Glory Kalu, Eric Wynn, Alan Y L Wong, Nadia Fröbisch, Petr Baranov, Maksim V Plikus, Jason D Buenrostro, Brian J Haas, Isaac M Chiu, Timothy B Sackton, Jessica L Whited","doi":"10.1101/2021.12.29.474455","DOIUrl":"10.1101/2021.12.29.474455","url":null,"abstract":"<p><p>Many species regenerate lost body parts following amputation. Most limb regeneration research has focused on the immediate injury site. Meanwhile, body-wide injury responses remain largely unexplored but may be critical for regeneration. Here, we discovered a role for the sympathetic nervous system in stimulating a body-wide stem cell activation response to amputation that drives enhanced limb regeneration in axolotls. This response is mediated by adrenergic signaling, which coordinates distant cellular activation responses via the α<sub>2A</sub>-adrenergic receptor, and local regeneration responses via β-adrenergic receptors. Both α<sub>2A</sub>- and β-adrenergic signaling act upstream of mTOR signaling. Notably, systemically-activated axolotls regenerate limbs faster than naïve animals, suggesting a potential selective advantage in environments where injury from cannibalism or predation is common. This work challenges the predominant view that cellular responses underlying regeneration are confined to the injury site and argues instead for body-wide cellular priming as a foundational step that enables localized tissue regrowth.</p>","PeriodicalId":72407,"journal":{"name":"bioRxiv : the preprint server for biology","volume":"10 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2025-07-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12330692/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"82563427","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-07-23DOI: 10.1101/2022.10.14.512297
Vishaka Gopalan, Chun Wai Wong, Rotem Leshem, Luke Owen, Tuulia Vallius, Yingxiao Shi, Yuhong Jiang, Eva Pérez-Guijarro, Emily Wu, Sung Chin, Jessica Ebersole, Cari Smith, Antonella Sassano, Maira Alves Constantino, Michael J Haley, Ferenc Livak, R Mark Simpson, Chi-Ping Day, Adam Hurlstone, Sridhar Hannenhalli, Glenn Merlino, Kerrie L Marie
Melanoma plasticity, driven by phenotype state switching, underlies clinically relevant traits such as metastasis and therapy resistance. As melanoma progression is thought to recapitulate aspects of neural crest cell (NCC) development, understanding embryonic melanocyte specification and lineage fate decisions of closely related NCCs may illuminate the pathways co-opted during disease evolution. Here, we use a mouse model to isolate and sequence Dopachrome tautomerase (Dct) expressing NCCs, the precursors of melanocytes, at two key developmental stages. We classify these lineages and devise a Developmental Gene Module (DGM) scoring system to interrogate lineage state switching in melanoma samples. In bulk transcriptomes, activation of DGMs representing embryonic Schwann Cell Precursors (SCPs)-multipotent stem cells-in patient tumors predicts poor response to immune checkpoint inhibitors (ICI). Co-activation of SCP and Mesenchymal-like (Mes.) modules further correlates with resistance to MAPK inhibitors. Notably, single-cell analyses reveal that melanoma cells can simultaneously express multiple DGMs, forming "hybrid" states. Cells in a hybrid Neural/SCP state are enriched in early metastasis and ICI-resistant tumors and are insensitive to inflammatory stimuli. We demonstrate that targeting Hdac2, a histone deacetylase associated with this Neural/SCP hybrid state, promotes a mesenchymal-like state switch, remodels the tumor microenvironment, and sensitizes melanoma cells to TNFα and tumors to ICI therapy. Our methodology thus reveals dynamic patterns of lineage state switching correlated with melanoma tumor evolution to drive insight into new therapeutic targets.
{"title":"Transitory Schwann Cell Precursor and hybrid states underpin melanoma therapy resistance and metastasis.","authors":"Vishaka Gopalan, Chun Wai Wong, Rotem Leshem, Luke Owen, Tuulia Vallius, Yingxiao Shi, Yuhong Jiang, Eva Pérez-Guijarro, Emily Wu, Sung Chin, Jessica Ebersole, Cari Smith, Antonella Sassano, Maira Alves Constantino, Michael J Haley, Ferenc Livak, R Mark Simpson, Chi-Ping Day, Adam Hurlstone, Sridhar Hannenhalli, Glenn Merlino, Kerrie L Marie","doi":"10.1101/2022.10.14.512297","DOIUrl":"10.1101/2022.10.14.512297","url":null,"abstract":"<p><p>Melanoma plasticity, driven by phenotype state switching, underlies clinically relevant traits such as metastasis and therapy resistance. As melanoma progression is thought to recapitulate aspects of neural crest cell (NCC) development, understanding embryonic melanocyte specification and lineage fate decisions of closely related NCCs may illuminate the pathways co-opted during disease evolution. Here, we use a mouse model to isolate and sequence Dopachrome tautomerase (Dct) expressing NCCs, the precursors of melanocytes, at two key developmental stages. We classify these lineages and devise a Developmental Gene Module (DGM) scoring system to interrogate lineage state switching in melanoma samples. In bulk transcriptomes, activation of DGMs representing embryonic Schwann Cell Precursors (SCPs)-multipotent stem cells-in patient tumors predicts poor response to immune checkpoint inhibitors (ICI). Co-activation of SCP and Mesenchymal-like (Mes.) modules further correlates with resistance to MAPK inhibitors. Notably, single-cell analyses reveal that melanoma cells can simultaneously express multiple DGMs, forming \"hybrid\" states. Cells in a hybrid Neural/SCP state are enriched in early metastasis and ICI-resistant tumors and are insensitive to inflammatory stimuli. We demonstrate that targeting <i>Hdac2</i>, a histone deacetylase associated with this Neural/SCP hybrid state, promotes a mesenchymal-like state switch, remodels the tumor microenvironment, and sensitizes melanoma cells to TNFα and tumors to ICI therapy. Our methodology thus reveals dynamic patterns of lineage state switching correlated with melanoma tumor evolution to drive insight into new therapeutic targets.</p>","PeriodicalId":72407,"journal":{"name":"bioRxiv : the preprint server for biology","volume":"9 4","pages":""},"PeriodicalIF":0.0,"publicationDate":"2025-07-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12330744/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91463601","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-07-22DOI: 10.1101/2023.03.04.531124
Marisa Carrasco, Caroline Myers, Mariel Roberts
For human adults, visual perception varies with polar angle at isoeccentric locations from the center of gaze. The same visual information yields better performance along the horizontal than vertical meridian (horizontal vertical anisotropy, HVA) and along the lower than upper vertical meridian (vertical meridian asymmetry, VMA). For children, performance is better along the horizontal than vertical meridian (HVA) but does not differ between the lower and upper vertical meridian. Here, we investigated whether the extent of the HVA varies and whether the VMA emerges and develops during adolescence, or whether the VMA only emerges in adulthood. We found that both the HVA and VMA develop gradually throughout adolescence and become as pronounced as those of adults only in late adolescence.
{"title":"Visual field asymmetries develop throughout adolescence.","authors":"Marisa Carrasco, Caroline Myers, Mariel Roberts","doi":"10.1101/2023.03.04.531124","DOIUrl":"10.1101/2023.03.04.531124","url":null,"abstract":"<p><p>For human adults, visual perception varies with polar angle at isoeccentric locations from the center of gaze. The same visual information yields better performance along the horizontal than vertical meridian (horizontal vertical anisotropy, HVA) and along the lower than upper vertical meridian (vertical meridian asymmetry, VMA). For children, performance is better along the horizontal than vertical meridian (HVA) but does not differ between the lower and upper vertical meridian. Here, we investigated whether the extent of the HVA varies and whether the VMA emerges and develops during adolescence, or whether the VMA only emerges in adulthood. We found that both the HVA and VMA develop gradually throughout adolescence and become as pronounced as those of adults only in late adolescence.</p>","PeriodicalId":72407,"journal":{"name":"bioRxiv : the preprint server for biology","volume":" ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2025-07-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10028823/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"9201952","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-07-21DOI: 10.1101/2023.07.13.548817
Giacomo Bignardi, Michel G Nivard, H Lina Schaare, Boris C Bernhardt, Richard A I Bethlehem, Simon E Fisher, Sofie L Valk
In humans, many neurobiological features of the cortex-including gene expression patterns, microstructure, and functional connectivity-vary systematically along a sensorimotor-association (S-A) axis of brain organisation. To date, it is still poorly understood whether inter-individual differences in patterns of S-A axis capture these robust spatial relationships across neurobiological properties observed at the group-level. Here, we examine inter-individual differences in structural and functional properties of the S-A axis, namely cortical microstructure, geodesic distances, and the functional gradient, in a sample of young adults from the Human Connectome Project (N = 992, including 328 twins). We quantified heritable variation associated with inter-individual differences in the S-A axis, and assessed whether structural and functional properties that are highly spatially correlated at the group-level also share genetic underpinnings. To consider measurement errors in resting-state functional connectivity data and their impact on properties of the S-A axis, we used a multivariate twin design capable of disentangling individual-level variation in both intra- and inter-individual differences. After accounting for some of the intra-individual variation, we found average heritable individual differences in both the functional gradient , cortical microstructure , and geodesic distances . However, these genetic influences were mostly distinct and deviated from group-level patterns. In particular, we found no significant genetic correlation between the functional gradient and microstructure, while we found both positive and negative genetic associations between the functional gradient and geodesic distances. Our approach highlights the complexity of genetic contributions to brain organisation and may have potential implications for understanding cognitive variability within the S-A axis framework.
{"title":"Distinct genetic underpinnings of inter-individual differences in the sensorimotor-association axis of cortical organisation.","authors":"Giacomo Bignardi, Michel G Nivard, H Lina Schaare, Boris C Bernhardt, Richard A I Bethlehem, Simon E Fisher, Sofie L Valk","doi":"10.1101/2023.07.13.548817","DOIUrl":"10.1101/2023.07.13.548817","url":null,"abstract":"<p><p>In humans, many neurobiological features of the cortex-including gene expression patterns, microstructure, and functional connectivity-vary systematically along a sensorimotor-association (S-A) axis of brain organisation. To date, it is still poorly understood whether inter-individual differences in patterns of S-A axis capture these robust spatial relationships across neurobiological properties observed at the group-level. Here, we examine inter-individual differences in structural and functional properties of the S-A axis, namely cortical microstructure, geodesic distances, and the functional gradient, in a sample of young adults from the Human Connectome Project (N = 992, including 328 twins). We quantified heritable variation associated with inter-individual differences in the S-A axis, and assessed whether structural and functional properties that are highly spatially correlated at the group-level also share genetic underpinnings. To consider measurement errors in resting-state functional connectivity data and their impact on properties of the S-A axis, we used a multivariate twin design capable of disentangling individual-level variation in both intra- and inter-individual differences. After accounting for some of the intra-individual variation, we found average heritable individual differences in both the functional gradient <math> <mrow> <mfenced> <mrow> <msup><msub><mi>h</mi> <mrow><mrow><mtext>twin</mtext></mrow> </mrow> </msub> <mn>2</mn></msup> <mo>=</mo> <mn>57</mn> <mtext>%</mtext></mrow> </mfenced> </mrow> </math> , cortical microstructure <math> <mrow> <mfenced> <mrow> <msup><msub><mi>h</mi> <mrow><mrow><mtext>twin</mtext></mrow> </mrow> </msub> <mn>2</mn></msup> <mo>=</mo> <mn>43</mn> <mtext>%</mtext></mrow> </mfenced> </mrow> </math> , and geodesic distances <math> <mrow> <mfenced> <mrow> <msup><msub><mi>h</mi> <mrow><mrow><mtext>twin</mtext></mrow> </mrow> </msub> <mn>2</mn></msup> <mo>=</mo> <mn>34</mn> <mtext>%</mtext></mrow> </mfenced> </mrow> </math> . However, these genetic influences were mostly distinct and deviated from group-level patterns. In particular, we found no significant genetic correlation between the functional gradient and microstructure, while we found both positive and negative genetic associations between the functional gradient and geodesic distances. Our approach highlights the complexity of genetic contributions to brain organisation and may have potential implications for understanding cognitive variability within the S-A axis framework.</p>","PeriodicalId":72407,"journal":{"name":"bioRxiv : the preprint server for biology","volume":"15 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2025-07-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12330504/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"84324053","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-07-21DOI: 10.1101/2022.02.24.481854
Robert M Flight, Praneeth S Bhatt, Hunter Nb Moseley
Background: Almost all correlation measures currently available are unable to directly handle missing values. Typically, missing values are either ignored completely by removing them or are imputed and used in the calculation of the correlation coefficient. In either case, the correlation value will be impacted based on a perspective that the missing data represents no useful information. However, missing values occur in real data sets for a variety of reasons. In omics data sets that are derived from analytical measurements, a major reason for missing values is that a specific measurable phenomenon falls below the detection limits of the analytical instrumentation (left-censored values). These missing data are not missing at random, but represent useful information by virtue of their "missingness" at one end of the data distribution.
Results: To include this information due to left-censorship missingness, we propose the information-content-informed Kendall-tau (ICI-Kt) methodology. We show how left-censored missing values can be included within the definition of the Kendall-tau correlation coefficient, and how that inclusion leads to an interpretation of information being added to the correlation. We also implement calculations for additional measures of theoretical maxima and pairwise completeness that add further layers of information interpretation in the methodology. Using both simulated and real data sets from RNA-seq, metabolomics, and lipidomics experiments, we demonstrate that the ICI-Kt methodology allows for the inclusion of left-censored missing data values as interpretable information, enabling both improved determination of outlier samples and improved feature-feature network construction. We provide explicitly parallel implementations in both R and Python that allow fast calculations of all the variables used when applying the ICI-Kt methodology on large numbers of samples.
Conclusions: The ICI-Kt methods are available as an R package and Python module on GitHub at https://github.com/moseleyBioinformaticsLab/ICIKendallTau and https://github.com/moseleyBioinformaticsLab/icikt, respectively.
{"title":"Information-Content-Informed Kendall-tau Correlation Methodology: Interpreting Missing Values as Useful Information.","authors":"Robert M Flight, Praneeth S Bhatt, Hunter Nb Moseley","doi":"10.1101/2022.02.24.481854","DOIUrl":"10.1101/2022.02.24.481854","url":null,"abstract":"<p><strong>Background: </strong>Almost all correlation measures currently available are unable to directly handle missing values. Typically, missing values are either ignored completely by removing them or are imputed and used in the calculation of the correlation coefficient. In either case, the correlation value will be impacted based on a perspective that the missing data represents no useful information. However, missing values occur in real data sets for a variety of reasons. In omics data sets that are derived from analytical measurements, a major reason for missing values is that a specific measurable phenomenon falls below the detection limits of the analytical instrumentation (left-censored values). These missing data are not missing at random, but represent useful information by virtue of their \"missingness\" at one end of the data distribution.</p><p><strong>Results: </strong>To include this information due to left-censorship missingness, we propose the information-content-informed Kendall-tau (ICI-Kt) methodology. We show how left-censored missing values can be included within the definition of the Kendall-tau correlation coefficient, and how that inclusion leads to an interpretation of information being added to the correlation. We also implement calculations for additional measures of theoretical maxima and pairwise completeness that add further layers of information interpretation in the methodology. Using both simulated and real data sets from RNA-seq, metabolomics, and lipidomics experiments, we demonstrate that the ICI-Kt methodology allows for the inclusion of left-censored missing data values as interpretable information, enabling both improved determination of outlier samples and improved feature-feature network construction. We provide explicitly parallel implementations in both R and Python that allow fast calculations of all the variables used when applying the ICI-Kt methodology on large numbers of samples.</p><p><strong>Conclusions: </strong>The ICI-Kt methods are available as an R package and Python module on GitHub at https://github.com/moseleyBioinformaticsLab/ICIKendallTau and https://github.com/moseleyBioinformaticsLab/icikt, respectively.</p>","PeriodicalId":72407,"journal":{"name":"bioRxiv : the preprint server for biology","volume":"38 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2025-07-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12330630/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"75088567","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-07-21DOI: 10.1101/2023.08.19.553999
Zitong Lu, Julie D Golomb
Remarkably, human brains have the ability to accurately perceive and process the real-world size of objects, despite vast differences in distance and perspective. While previous studies have delved into this phenomenon, distinguishing the processing of real-world size from other visual properties, like depth, has been challenging. Using the THINGS EEG2 dataset with human EEG recordings and more ecologically valid naturalistic stimuli, our study combines human EEG and representational similarity analysis to disentangle neural representations of object real-world size from retinal size and perceived depth, leveraging recent datasets and modeling approaches to address challenges not fully resolved in previous work. We report a representational timeline of visual object processing: object real-world depth processed first, then retinal size, and finally, real-world size. Additionally, we input both these naturalistic images and object-only images without natural background into artificial neural networks. Consistent with the human EEG findings, we also successfully disentangled representation of object real-world size from retinal size and real-world depth in all three types of artificial neural networks (visual-only ResNet, visual-language CLIP, and language-only Word2Vec). Moreover, our multi-modal representational comparison framework across human EEG and artificial neural networks reveals real-world size as a stable and higher-level dimension in object space incorporating both visual and semantic information. Our research provides a temporally resolved characterization of how certain key object properties - such as object real-world size, depth, and retinal size - are represented in the brain, which offers further advances and insights into our understanding of object space and the construction of more brain-like visual models.
{"title":"Human EEG and artificial neural networks reveal disentangled representations and processing timelines of object real-world size and depth in natural images.","authors":"Zitong Lu, Julie D Golomb","doi":"10.1101/2023.08.19.553999","DOIUrl":"10.1101/2023.08.19.553999","url":null,"abstract":"<p><p>Remarkably, human brains have the ability to accurately perceive and process the real-world size of objects, despite vast differences in distance and perspective. While previous studies have delved into this phenomenon, distinguishing the processing of real-world size from other visual properties, like depth, has been challenging. Using the THINGS EEG2 dataset with human EEG recordings and more ecologically valid naturalistic stimuli, our study combines human EEG and representational similarity analysis to disentangle neural representations of object real-world size from retinal size and perceived depth, leveraging recent datasets and modeling approaches to address challenges not fully resolved in previous work. We report a representational timeline of visual object processing: object real-world depth processed first, then retinal size, and finally, real-world size. Additionally, we input both these naturalistic images and object-only images without natural background into artificial neural networks. Consistent with the human EEG findings, we also successfully disentangled representation of object real-world size from retinal size and real-world depth in all three types of artificial neural networks (visual-only ResNet, visual-language CLIP, and language-only Word2Vec). Moreover, our multi-modal representational comparison framework across human EEG and artificial neural networks reveals real-world size as a stable and higher-level dimension in object space incorporating both visual and semantic information. Our research provides a temporally resolved characterization of how certain key object properties - such as object real-world size, depth, and retinal size - are represented in the brain, which offers further advances and insights into our understanding of object space and the construction of more brain-like visual models.</p>","PeriodicalId":72407,"journal":{"name":"bioRxiv : the preprint server for biology","volume":" ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2025-07-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://ftp.ncbi.nlm.nih.gov/pub/pmc/oa_pdf/56/01/nihpp-2023.08.19.553999v1.PMC10473678.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"10165932","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-07-16DOI: 10.1101/2023.11.19.567585
Katie E Holmes, David VanInsberghe, Lucas M Ferreri, Baptiste Elie, Ketaki Ganti, Chung-Young Lee, Anice C Lowen
For many viruses, narrow bottlenecks acting during transmission sharply reduce genetic diversity in a recipient host relative to the donor. Since genetic diversity represents adaptive potential, such losses of diversity are thought to limit the opportunity for viral populations to undergo antigenic change and other adaptive processes. Thus, a detailed picture of evolutionary dynamics during transmission is critical to understanding the forces driving viral evolution at an epidemiologic scale. To advance this understanding, we used a barcoded virus library and a guinea pig model of transmission to decipher where in the transmission process influenza A virus populations lose diversity. In inoculated guinea pigs, we show that a high level of viral barcode diversity is maintained. Within-host continuity in the barcodes detected across time furthermore indicates that stochastic effects are not pronounced within the inoculated hosts. Importantly, in both aerosol-exposed and direct contact animals, we observed many barcodes at the earliest time point(s) positive for infectious virus, indicating robust transfer of diversity through the environment. This high viral diversity is short-lived, however, with a sharp decline seen 1-2 days after initiation of infection. Although major losses of diversity at transmission are well described for influenza A virus, our data indicate that events that occur following viral transfer and during the earliest stages of natural infection have a central role in this process. This finding suggests that host factors, such as immune effectors, may have greater opportunity to impose selection during influenza A virus transmission than previously recognized.
{"title":"Viral expansion after transfer is a primary driver of influenza A virus transmission bottlenecks.","authors":"Katie E Holmes, David VanInsberghe, Lucas M Ferreri, Baptiste Elie, Ketaki Ganti, Chung-Young Lee, Anice C Lowen","doi":"10.1101/2023.11.19.567585","DOIUrl":"10.1101/2023.11.19.567585","url":null,"abstract":"<p><p>For many viruses, narrow bottlenecks acting during transmission sharply reduce genetic diversity in a recipient host relative to the donor. Since genetic diversity represents adaptive potential, such losses of diversity are thought to limit the opportunity for viral populations to undergo antigenic change and other adaptive processes. Thus, a detailed picture of evolutionary dynamics during transmission is critical to understanding the forces driving viral evolution at an epidemiologic scale. To advance this understanding, we used a barcoded virus library and a guinea pig model of transmission to decipher where in the transmission process influenza A virus populations lose diversity. In inoculated guinea pigs, we show that a high level of viral barcode diversity is maintained. Within-host continuity in the barcodes detected across time furthermore indicates that stochastic effects are not pronounced within the inoculated hosts. Importantly, in both aerosol-exposed and direct contact animals, we observed many barcodes at the earliest time point(s) positive for infectious virus, indicating robust transfer of diversity through the environment. This high viral diversity is short-lived, however, with a sharp decline seen 1-2 days after initiation of infection. Although major losses of diversity at transmission are well described for influenza A virus, our data indicate that events that occur following viral transfer and during the earliest stages of natural infection have a central role in this process. This finding suggests that host factors, such as immune effectors, may have greater opportunity to impose selection during influenza A virus transmission than previously recognized.</p>","PeriodicalId":72407,"journal":{"name":"bioRxiv : the preprint server for biology","volume":" ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2025-07-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10680852/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138447445","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-07-11DOI: 10.1101/2023.11.01.565061
John W Benning, Jedidiah Carlson, Olivia S Smith, Ruth G Shaw, Arbel Harpak
The scientific literature has seen a resurgence of interest in genetic influences on human behavior and socioeconomic outcomes. Such studies face the central difficulty of distinguishing possible causal influences, in particular genetic and non-genetic ones. When confounding between possible influences is not rigorously addressed, it invites over- and misinterpretation of data. We illustrate the breadth of this problem through a discussion of the literature and a reanalysis of two examples. Clark (2023) suggested that patterns of similarity in social status between relatives indicate that social status is largely determined by one's DNA. We show that the paper's conclusions are based on the conflation of genetic and non-genetic transmission (for example, of wealth) within families. Song & Zhang (2024) posited that genetic variants underlying bisexual behavior are maintained in the population because they also affect risk-taking behavior, thereby conferring an evolutionary fitness advantage through increased sexual promiscuity. In this case, too, we show that possible explanations cannot be distinguished, but only one is chosen and presented as a conclusion. We discuss how issues of confounding apply more broadly to studies that claim to establish genetic underpinnings to human behavior and societal outcomes.
{"title":"Confounding Fuels Misinterpretation in Human Genetics.","authors":"John W Benning, Jedidiah Carlson, Olivia S Smith, Ruth G Shaw, Arbel Harpak","doi":"10.1101/2023.11.01.565061","DOIUrl":"10.1101/2023.11.01.565061","url":null,"abstract":"<p><p>The scientific literature has seen a resurgence of interest in genetic influences on human behavior and socioeconomic outcomes. Such studies face the central difficulty of distinguishing possible causal influences, in particular genetic and non-genetic ones. When confounding between possible influences is not rigorously addressed, it invites over- and misinterpretation of data. We illustrate the breadth of this problem through a discussion of the literature and a reanalysis of two examples. Clark (2023) suggested that patterns of similarity in social status between relatives indicate that social status is largely determined by one's DNA. We show that the paper's conclusions are based on the conflation of genetic and non-genetic transmission (for example, of wealth) within families. Song & Zhang (2024) posited that genetic variants underlying bisexual behavior are maintained in the population because they also affect risk-taking behavior, thereby conferring an evolutionary fitness advantage through increased sexual promiscuity. In this case, too, we show that possible explanations cannot be distinguished, but only one is chosen and presented as a conclusion. We discuss how issues of confounding apply more broadly to studies that claim to establish genetic underpinnings to human behavior and societal outcomes.</p>","PeriodicalId":72407,"journal":{"name":"bioRxiv : the preprint server for biology","volume":" ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2025-07-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10635045/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"92157654","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}