Spatiotemporal evolution of the distribution of Chronic bee paralysis virus (CBPV) in honey bee colonies.

Chronic bee paralysis virus (CBPV) is a Apis mellifera viral infectious disease, exhibiting dark and hairless abdomen in workers with tremors and ataxita. Clinical signs are also typically linked to adverse weather conditions and overcrowding in the hive. The disease occurs in spring but recently it has been observed cases increase and seasonality loss of the disease incidence. This study analyses the evolution of CBPV in Italy, through data collected from 2009 to 2023 within three monitoring projects comprising nationwide extended detection networks, aimed to investigate the evolution of the CBPV spatial distribution, identifying high-risk areas for the virus spread. This study highlights an increased risk over years. Prevalence increased from 4.3% during 2009-2010 to 84.7% during 2021-2023 monitoring years. CBPV outbreaks were irregular between investigated seasons, highlighting Spring and Autumn as the most susceptible seasons. Risk of CBPV infection has increased, reaching high-risk in last years of monitoring. Sequence analysis showed a high similarity to other isolated Italian CBPVs. The study offers an epidemiological insight into the aetiology of this disease. CBPV distribution is a prerequisite to predict its future spread and factors involved in its propagation not only in honey bees but also in other pollinators and environments.

Spatial transcriptomics reveals organized and distinct immune activation in cutaneous granulomatous disorders.

BACKGROUND: Non-infectious (inflammatory) cutaneous granulomatous disorders include cutaneous sarcoidosis (CS), granuloma annulare (GA), necrobiosis lipoidica (NL), and necrobiotic xanthogranuloma (NXG). These disorders share macrophage predominant inflammation histologically, but the inflammatory architecture and the pattern of extracellular matrix alteration varies. The underlying molecular explanations for these differences remain unclear. OBJECTIVE: To understand spatial gene expression characteristics in these disorders. METHODS: We performed spatial transcriptomics in cases of CS, GA, NL, and NXG to compare patterns of immune activation and other molecular features in a spatially resolved fashion. RESULTS: CS is characterized by a polarized, spatially organized T helper (Th) 1 predominant response with classical macrophage activation. GA is characterized by a mixed, but spatially organized pattern of Th1 and Th2 polarization with both classical and alternative macrophage activation. NL showed concomitant activation of Th1, Th2, and Th17 immunity with a mixed pattern of macrophage activation. Activation of type 1 immunity was shared among, CS, GA, and NL and included upregulation of IL-32. NXG showed upregulation of CXCR4-CXCL12/14 chemokine signaling and exaggerated alternative macrophage polarization. Histologic alteration of extracellular matrix correlated with hypoxia and glycolysis programs and type 2 immune activation. CONCLUSIONS: Inflammatory cutaneous granulomatous disorders show distinct and spatially organized immune activation that correlate with hallmark histologic changes.

Improving spatial exposure data for everyone - lifecourse social context and ascertaining residential history.

Measuring age-specific, contextual exposures is crucial for lifecourse epidemiology research. Longitudinal residential data offers a "golden ticket" to cumulative exposure metrics and can enhance our understanding of health disparities. Residential history can be linked to myriad spatiotemporal databases to characterize environmental, socioeconomic, and policy contexts that a person experienced throughout life. However, obtaining accurate residential history is challenging in the United States due to the limitations of administrative registries and self-reports. Xu et al. (Am J Epidemiol. 2024; 193(2):348-359) detail an approach to linking residential history sourced from LexisNexis ® Accurint ® to a Wisconsin-based research cohort, offering insights into challenges with residential history collection. Researchers must analyze the magnitude of selection and misclassification biases inherent to ascertaining residential history from cohort data. A lifecourse framework can provide insights into why the frequency and distance of moves is patterned by age, birth cohort, racial/ethnic identity, socioeconomic status, and urbanicity. Historic and contemporary migration patterns of marginalized people seeking economic and political opportunities must guide interpretations of residential history data. We outline methodologic priorities for use of residential history in health disparities research, including contextualizing residential history data with determinants of residential moves, triangulating spatial exposure assessment methods, and transparently quantifying measurement error.

Mapping of on-field soil nutrient variabilities as a guiding force for smart farming: a case study from FarmerZone sentinel-1 from three potato agroecological zones of India.

Mapping of soil nutrient parameters using experimental measurements and geostatistical approaches to assist site-specific fertiliser advisories is anticipated to play a significant role in Smart Agriculture. FarmerZone is a cloud service envisioned by the Department of Biotechnology, Government of India, to provide advisories to assist smallholder farmers in India in enhancing their overall farm production. As a part of the project, we evaluated the soil spatial variability of three potato agroecological zones in India and provided soil health cards along with field-specific fertiliser recommendations for potato cultivation to farmers. Specifically, 705 surface samples were collected from three representative potato-growing districts of Indian states (Meerut, UP; Jalandhar, Punjab and Lahaul and Spiti, HP) and analysed for soil parameters such as organic carbon, macronutrients (NPK), micronutrients (Zn, Fe, Mn, and Cu), pH, and EC. The soil parameters were integrated into a geodatabase and subjected to kriging interpolation to create spatial soil maps of the targeted potato agroecological zones through best-fit experimental semivariograms. The spatial distribution showed a deficiency of soil organic carbon in two studied zones and available nitrogen among all studied zones. The available phosphorus and potassium varied among the agroecological zones. The micronutrient levels were largely sufficient in all the zones except at a few specific sites where nutrient advisories are recommended to replenish. The general management strategies were recommended based on the nutrient status in the studied area. This study clearly supports the significance of site-specific soil analytics and interpolated spatial soil mapping over any targeted agroecological zones as a promising strategy to deliver reliable advisories of fertiliser recommendations for smart farming.

Open-source, high-throughput targeted in-situ transcriptomics for developmental and tissue biology.

Multiplexed spatial profiling of mRNAs has recently gained traction as a tool to explore the cellular diversity and the architecture of tissues. We propose a sensitive, open-source, simple and flexible method for the generation of in-situ expression maps of hundreds of genes. We exploit direct ligation of padlock probes on mRNAs, coupled with rolling circle amplification and hybridization-based in situ combinatorial barcoding, to achieve high detection efficiency, high throughput and large multiplexing. We validate the method across a number of species, and show its use in combination with orthogonal methods such as antibody staining, highlighting its potential value for developmental and tissue biology studies. Finally, we provide an end-to-end computational workflow that covers the steps of probe design, image processing, data extraction, cell segmentation, clustering and annotation of cell types. By enabling easier access to high-throughput spatially resolved transcriptomics, we hope to encourage a diversity of applications and the exploration of a wide range of biological questions.

Achieving Image Encryption Quantum Dot-Functionalized Encryption Camera with Designed Films.

The risk of information leaks increases as images become a crucial medium for information sharing. There is a great need to further develop the versatility of image encryption technology to protect confidential and sensitive information. Herein, using high spatial redundancy (strong correlation of neighboring pixels) of the image and the in situ encryption function of a quantum dot functionalized encryption camera, in situ image encryption is achieved by designing quantum dot films (size, color, and full width at half maximum) to modify the correlation and reduce spatial redundancy of the captured image during encryption processing. The correlation coefficients of simulated encrypted image closely apporach to 0. High-quality decrypted images are achieved with a PSNR of more than 35 dB by a convolutional neural network-based algorithm that meets the resolution requirements of human visual perception. Compared with the traditional image encryption algorithms, chaotic image encryption algorithms and neural network-based encryption algorithms described previously, it provides a universal, efficient and effective in situ image encryption method.

Visualization of Pulpal Structures by SWIR in Endodontic Access Preparation.

Endodontic access preparation is one of the initial steps in root canal treatments and can be hindered by the obliteration of pulp canals and formation of tertiary dentin. Until now, methods for direct intraoperative visualization of the 3-dimensional anatomy of teeth have been missing. Here, we evaluate the use of shortwave infrared radiation (SWIR) for navigation during stepwise access preparation. Nine teeth (3 anteriors, 3 premolars, and 3 molars) were explanted en bloc with intact periodontium including alveolar bone and mucosa from the upper or lower jaw of human body donors. Analysis was performed at baseline as well as at preparation depths of 5 mm, 7 mm, and 9 mm, respectively. For reflection, SWIR was used at a wavelength of 1,550 nm from the occlusal direction, whereas for transillumination, SWIR was passed through each sample at the marginal gingiva from the buccal as well as oral side at a wavelength of 1,300 nm. Pulpal structures could be identified as darker areas approximately 2 mm before reaching the pulp chamber using SWIR transillumination, although they were indistinguishable under normal circumstances. Furcation areas in molars appeared with higher intensity than areas with canals. The location of pulpal structures was confirmed by superimposition of segmented micro-computed tomography (µCT) images. By radiomic analysis, significant differences between pulpal and parapulpal areas could be detected in image features. With hierarchical cluster analysis, both segments could be confirmed and associated with specific clusters. The local thickness of µCTs was calculated and correlated with SWIR transillumination images, by which a linear dependency of thickness and intensity could be demonstrated. Lastly, by in silico simulations of light propagation, dentin tubules were shown to be a crucial factor for understanding the visibility of the pulp. In conclusion, SWIR transillumination may allow direct clinical live navigation during endodontic access preparation.

Drivers of change behind the spatial distribution and fate of typical trace organic pollutants in fresh waste leachate across China.

There have been growing concerns regarding the health and environmental impacts of trace organic pollutants (TOPs). However, fresh leachate from municipal solid waste (MSW) has been overlooked as a potential reservoir of TOPs. Therefore, we investigated 90 legacy and emerging TOPs in fresh leachate from 14 provinces and municipalities in China. Additionally, the fate and final discharge impacts of TOPs in 14 leachate treatment systems were analyzed. The results revealed that the detection rate of 90 TOPs was over 50 % in all samples. Notably, polychlorinated biphenyls, banned for 40 years, were frequently detected in fresh leachate. The concentration of pseudo-persistent TOPs (105-107 ng/L) is significantly higher than that of persistent TOPs (102-104 ng/L). Spatial distribution patterns of TOPs in fresh leachate suggest that economy, population, climate, and policies impact TOPs discharge from MSW. For example, economically developed and densely populated areas displayed higher TOPs concentrations, whereas warmer climates facilitate TOPs leaching from MSW. We confirmed that waste classification policies were a key driver of the decline in multiple TOPs in leachate. Mass balance analysis shows that the final effluent and sludge from current dominant leachate treatment systems contain refractory TOPs, especially perfluoroalkyl acids, which must be prioritized for control. This paper was the first comprehensive investigation of multiple TOPs in fresh leachate at a large geographic scale. The factors affecting the occurrence, spatial distribution, and fate of TOPs in fresh leachate were revealed. It provides a valuable reference for the establishment of policies for the management of TOPs in MSW and the associated leachate.

Spatial transcriptomics: Integrating morphology and molecular mechanisms of kidney diseases.

The recent arrival of high-resolution spatial transcriptomics (ST) technologies is generating a veritable revolution in life sciences, enabling biomolecules to be measured in their native spatial context. By integrating morphology and molecular biology, ST technologies offer the potential of improving the understanding of tissue biology and disease and may also provide meaningful clinical insights. In this review, we describe the main spatial transcriptomics technologies currently available, the computational analysis for data interpretation and visualization, and illustrate their scientific and potential medical interest in the context of kidney disease. Finally, we discuss the perspectives and challenges of these booming new technologies.

The impact of response type on affordance and spatial compatibility effects in human and object interactions.

Everyday social interactions or goal-directed actions towards objects activate action plans appropriate to their affordances. The spatial compatibility of a stimulus and a response might interfere with the activation of these action plans. In the present study, we examined how framing of interactions affects the interplay between affordance and spatial compatibility effects towards humans and objects in two separate experiments. In a motor priming task designed to simultaneously assess these two effects, participants were presented with interactive hand gestures and objects with a single handle. Participants responded either with their left or right hand according to the colour mask of the stimulus, regardless of the spatial position or the affordance-related orientation of the stimulus. In Experiment 1, when responses were given by keypresses, we found independent affordance and spatial compatibility effects towards objects. Surprisingly, interactive hand gestures induced a reversed affordance effect, that is, imitative action tendencies. Changing the responses from keypresses to the performance of grasping actions in Experiment 2 drastically altered these findings, resulting in the enhancement of affordance and the elimination of spatial compatibility effects for both human and object interactions. These findings highlight the importance of contextual influences on the emergence of automatic action tendencies.

Changes in urban green space configuration and connectivity using spatial graph-based metrics in Ardabil developing city, Iran.

Urban planning is essential for managing the diverse impacts of urban green spaces, such as public access, stormwater control, urban life quality, and landscape aesthetics, promoting sustainable urban development and urban residents' well-being by integrating green space considerations into city planning. The aim of this study is to use graph-based metrics to calculate the connectivity of UGS across the main municipal zones of Ardabil city over consecutive periods under different population growth rates. Another objective of this study is to compare the connectivity values of UGS in the four municipal zones and to evaluate changes in the connectivity indices at various distance thresholds of UGS patches. After identifying UGS in different periods, the changes in graph-based connectivity indices at various distance thresholds of UGS patches were analyzed. Additionally, the changes in connectivity indices over different periods and across various municipal zones were compared and analyzed. The findings reveal that UGS areas were larger in the past but have recently had smaller patch sizes. Connectivity between UGS nodes (dNL) decreased at various distances over the study years, showing a declining trend in different connectivity indices. UGS connectivity decreased in municipal zones 1, 2, and 3 but increased in recent years after a decline until 2012 across all four zones of Ardabil city. Zone 4 had the highest UGS connectivity due to newly developed urban areas and well-allocated UGSs. Integrating the ecological impacts of UGS connectivity in urban development and design will enhance trade-offs between conservation, public health, and social equity. New urban areas should allocate sufficient land for UGS and parks, ensuring accessibility to support health and leisure through municipal planning. The study highlights the need for sustainable urban development policies that prioritize the allocation and maintenance of UGSs.

A possible origin of the inverted vertebrate retina revealed by physical modeling.

The evolutionary origin of the inverted retina in the vertebrate eye is unknown. This paper explores a hypothetical evolutionary scenario that explains the unique orientation of the photoreceptors in the vertebrate retina. The proposed scenario follows the scientific accepted scenario for eye evolution and gradually builds up towards an eye prototype by considering light direction detection and increase in achievable spatial resolution as the driving forces. It suggests that eye retinas developed along two different morphological processes, an evagination process that results in the inverted retina in vertebrate eyes and an invagination process that results in a verted retina in cephalopod eyes. The development of the inverted vertebrate retina and eye prototype morphology is strongly substantiated by physics of vision. The proposed evolutionary sequence for vertebrate eye development is simple and has the full potential to explain the origin of the inverted retina and leads to an eye prototype enabling visual detection and orientation. It allows the emergence of eye structures like, extraocular muscles, tapetum lucidum, biconvex lens, cornea, and pupil. This study supports the suggestion that a primitive inverted retina in the predecessor of vertebrates is of ectodermal origin and available before neurulation occurred.

SARC-UNet: A coronary artery segmentation method based on spatial attention and residual convolution.

BACKGROUND AND OBJECTIVE: Coronary artery segmentation is a pivotal field that has received increasing attention in recent years. However, this task remains challenging because of the inhomogeneous distributions of the contrast agent and dim light, resulting in noise, vascular breakages and small vessel losses in the obtained segmentation results. METHODS: To acquire better automatic blood vessel segmentation results for coronary angiography images, a UNet-based segmentation network (SARC-UNet) is constructed for coronary artery segmentation; this approach is based on residual convolution and spatial attention. First, we use the low-light image enhancement (LIME) approach to increase the contrast and clarity levels of coronary angiography images. Then, we design two residual convolution fusion modules (RCFM1 and RCFM2) that can successfully fuse the local and global information of coronary images while also capturing the characteristics of finer-grained blood vessels, hence preventing the loss of tiny blood vessels in the segmentation findings. Finally, using a cascaded waterfall structure, we create a new location-enhanced spatial attention (LESA) mechanism that can efficiently improve the long-distance dependencies between coronary vascular pixel features, eradicating vascular ruptures and noise in the segmentation results. RESULTS: This article subjectively and objectively evaluates the experimental results. This method has performed well on five general indicators. Furthermore, it outperforms the connectivity indicators proposed in this article. This method can effectively segment blood vessels and obtain higher accuracy results. CONCLUSIONS: Numerous experiments have shown that the suggested method outperforms the state-of-the-art approaches, particularly in terms of vessel connectivity and small blood vessel segmentation.

aKNNO: single-cell and spatial transcriptomics clustering with an optimized adaptive k-nearest neighbor graph.

Typical clustering methods for single-cell and spatial transcriptomics struggle to identify rare cell types, while approaches tailored to detect rare cell types gain this ability at the cost of poorer performance for grouping abundant ones. Here, we develop aKNNO to simultaneously identify abundant and rare cell types based on an adaptive k-nearest neighbor graph with optimization. Benchmarking on 38 simulated and 20 single-cell and spatial transcriptomics datasets demonstrates that aKNNO identifies both abundant and rare cell types more accurately than general and specialized methods. Using only gene expression aKNNO maps abundant and rare cells more precisely compared to integrative approaches.

Spatial effects of infiltrating T cells on neighbouring cancer cells and prognosis in stage III CRC patients.

Colorectal cancer (CRC) is one of the most frequently occurring cancers, but prognostic biomarkers identifying patients at risk of recurrence are still lacking. In this study, we aimed to investigate in more detail the spatial relationship between intratumoural T cells, cancer cells, and cancer cell hallmarks as prognostic biomarkers in stage III colorectal cancer patients. We conducted multiplexed imaging of 56 protein markers at single-cell resolution on resected fixed tissue from stage III CRC patients who received adjuvant 5-fluorouracil (5FU)-based chemotherapy. Images underwent segmentation for tumour, stroma, and immune cells, and cancer cell 'state' protein marker expression was quantified at a cellular level. We developed a Python package for estimation of spatial proximity, nearest neighbour analysis focusing on cancer cell-T-cell interactions at single-cell level. In our discovery cohort (Memorial Sloan Kettering samples), we processed 462 core samples (total number of cells: 1,669,228) from 221 adjuvant 5FU-treated stage III patients. The validation cohort (Huntsville Clearview Cancer Center samples) consisted of 272 samples (total number of cells: 853,398) from 98 stage III CRC patients. While there were trends for an association between the percentage of cytotoxic T cells (across the whole cancer core), it did not reach significance (discovery cohort: p = 0.07; validation cohort: p = 0.19). We next utilised our region-based nearest neighbour approach to determine the spatial relationships between cytotoxic T cells, helper T cells, and cancer cell clusters. In both cohorts, we found that shorter distance between cytotoxic T cells, T helper cells, and cancer cells was significantly associated with increased disease-free survival. An unsupervised trained model that clustered patients based on the median distance between immune cells and cancer cells, as well as protein expression profiles, successfully classified patients into low-risk and high-risk groups (discovery cohort: p = 0.01; validation cohort: p = 0.003). © 2024 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Understanding air pollution reduction from the perspective of synergy with carbon mitigation in China from 2008 to 2017.

In response to the dual challenges of air pollution control and carbon mitigation, China has strategically shifted its focus towards the synergistic reduction of air pollutants and CO2 emissions. This study identifies the potential areas and specific air pollutant species (including CO, NOx, and SO2) for co-reduction with carbon mitigation. We also reveal the driving forces behind the emissions of each air pollutant at both the national and regional scales. Our findings are as follows: (1) The potential for synergistic reduction of CO and SO2 with CO2 emissions has diminished in economically developed areas. There is a significant opportunity for co-reduction of SO2 and CO2 in the western and northern regions of China, particularly within Heilongjiang Province. (2) NOx is the key species for synergistic reduction with CO2 emissions across China, especially in the Chengyu Plain. (3) Cleaner production and the synergistic reduction effect are the primary contributors to national air pollutant reduction in China from 2008 to 2017. Conversely, efforts in economic development and energy efficiency have led to emission increases. Energy and industrial structures have only made limited contributions to emission reductions, and carbon mitigation shows an inhibition effect on emission reductions. These results offer valuable insights for developing targeted regional strategies for deeper air pollution control, considering the specific characteristics and needs of each region. Additionally, our findings highlight the importance of addressing policy misalignments and strengthening mutual-influence mechanisms between air pollution control and carbon mitigation, ensuring that policies for carbon reduction also effectively contribute to air quality improvements.

A single-cell transcriptomic map of the developing Atoh1 lineage identifies neural fate decisions and neuronal diversity in the hindbrain.

Proneural transcription factors establish molecular cascades to orchestrate neuronal diversity. One such transcription factor, Atonal homolog 1 (Atoh1), gives rise to cerebellar excitatory neurons and over 30 distinct nuclei in the brainstem critical for hearing, breathing, and balance. Although Atoh1 lineage neurons have been qualitatively described, the transcriptional programs that drive their fate decisions and the full extent of their diversity remain unknown. Here, we analyzed single-cell RNA sequencing and ATOH1 DNA binding in Atoh1 lineage neurons of the developing mouse hindbrain. This high-resolution dataset identified markers for specific brainstem nuclei and demonstrated that transcriptionally heterogeneous progenitors require ATOH1 for proper migration. Moreover, we identified a sizable population of proliferating unipolar brush cell progenitors in the mouse Atoh1 lineage, previously described in humans as the origin of one medulloblastoma subtype. Collectively, our data provide insights into the developing mouse hindbrain and markers for functional assessment of understudied neuronal populations.

The FXR1 network acts as a signaling scaffold for actomyosin remodeling.

It is currently not known whether mRNAs fulfill structural roles in the cytoplasm. Here, we report the fragile X-related protein 1 (FXR1) network, an mRNA-protein (mRNP) network present throughout the cytoplasm, formed by FXR1-mediated packaging of exceptionally long mRNAs. These mRNAs serve as an underlying condensate scaffold and concentrate FXR1 molecules. The FXR1 network contains multiple protein binding sites and functions as a signaling scaffold for interacting proteins. We show that it is necessary for RhoA signaling-induced actomyosin reorganization to provide spatial proximity between kinases and their substrates. Point mutations in FXR1, found in its homolog FMR1, where they cause fragile X syndrome, disrupt the network. FXR1 network disruption prevents actomyosin remodeling-an essential and ubiquitous process for the regulation of cell shape, migration, and synaptic function. Our findings uncover a structural role for cytoplasmic mRNA and show how the FXR1 RNA-binding protein as part of the FXR1 network acts as an organizer of signaling reactions.

Suitability of rainwater harvesting in saline and arsenic affected areas of Bangladesh.

A major portion of Bangladesh is currently experiencing a scarcity of safe drinking water because of arsenic contamination, high salinity and human-induced pollution. The objectives of this study were to identify locations with a high scarcity of drinking water and suitability of harvesting rainwater. Kriging interpolation algorithms of Geographical Information System (GIS) was employed to identify the probable water scarce zones as well as suitable zones of harvesting rain water from the available data of secondary sources. Statistical methods were employed to cluster, correlate, and regress variables such as rainfall, salinity, and As. The results showed that groundwater quality in the southwestern parts of Bangladesh is saline with high concentration (>10000 µS/cm). On the other hand, the northeastern and southwestern parts of Bangladesh are also vulnerable to arsenic contamination (60 %-97 % of tubewells), compared to other regions. The rainfall zonation map, covering the years 1951-2022, indicated that the Sylhet division had the highest potential for rainfall (ranging from 2600 to 3900 mm). From this study it was demonstrated that Sylhet, Noakhali, Bhola, Barishall, Patuakhali, Bagerhat, and Khulna were identified as suitable places for sustainable rainwater harvesting (RWH). The findings of this study may play significant role towards achieving sustainable potable water supply in vulnerable zones, if they receive attention from policymakers.

Deciphering the microbial succession and color formation mechanism of "green-covering and red-heart" Guanyin Tuqu.

"Green-covering and red-heart" Guanyin Tuqu (GRTQ), as a type of special fermentation starter, is characterized by the "green-covering" formed on the surface of Guanyin Tuqu (SQ) and the "red-heart" in the center of Guanyin Tuqu (CQ). However, the mechanisms that promote temporal succession in the GRTQ microbial ecology and the formation of "green-covering and red-heart" characteristics remain unclear. Herein, we correlated the temporal profiles of microbial community succession with the main environmental variables (temperature, moisture, and acidity) and spatial position (center and surface) in GRTQ throughout fermentation. According to the results of high-throughput sequencing and culture-dependent methods, the microbial communities in the CQ and SQ demonstrated functional complementarity. For instance, the bacterial richness index of the CQ was greater than that of SQ, and the fungal richness index of the SQ was greater than that of CQ at the later stage of fermentation. Furthermore, Saccharomycopsis, Saccharomyces, Aspergillus, Monascus, Lactobacillus, Bacillus, Rhodanobacter, and Chitinophaga were identified as the dominant microorganisms in the center, while the surface was represented by Saccharomycopsis, Aspergillus, Monascus, Lactobacillus, Acetobacter, and Weissella. By revealing the physiological characteristics of core microorganisms at different spatial positions of GRTQ, such as Aspergillus clavatus and Monascus purpureus, as well as their interactions with environmental factors, we elucidated the color formation mechanism behind the phenomenon of "green" outside and "red" inside. This study provides fundamental information support for optimizing the production process of GRTQ.