Phylogenetic Analysis of Wall-Associated Kinase Genes in Triticum Species and Characterization of TaWAK7 Involved in Wheat Powdery Mildew Resistance

Wall-associated kinases (WAKs), a group of receptor-like kinases, have been found to play important roles in defending against pathogens and in various developmental processes. However, the importance of this family in wheat remains largely unknown. Wheat powdery mildew is caused by Blumeria graminis f. sp. tritici (Bgt), which initiates infection on the cell surface and forms haustoria inside the cell; therefore, the defense to Bgt involves extracellular and subsequently intracellular signals. In this study, WAKs were identified genome-wide and analyzed phylogenetically, and then a transmembrane WAK gene that putatively participated in pathogen-associated molecular pattern-triggered immunity and effector-triggered immunity to Bgt was functionally and evolutionarily investigated. In total, 1,193 WAKs were identified from wheat and its Gramineae relatives. Phylogenetic analysis indicated that WAKs expanded through tandem duplication or segment duplication. TaWAK7, from chromosome 2A, was identified as a Bgt-inducible gene both in susceptible and resistant materials, but it showed distinct responsive patterns. Functional analysis showed that TaWAK7 was involved in both the basal and resistance gene-mediated resistances. The specific gene structures and protein characteristics of TaWAK7, along with its orthologs, were characterized both in subgenomes of Triticum spp. and in the A genome of multiple wheat accessions, which revealed that TaWAK7 orthologs underwent complex evolution with frequent gene fusion and domain deletion. In addition, three cytoplasmic proteins interacting with TaWAK7 were indicated by yeast two-hybrid and bimolecular fluorescence complementation assays. Binding of TaWAK7 with these proteins could change its subcellular localization from the plasma membrane to the cytoplasm. This study provides a better understanding of the evolution of WAKs at the genomic level and TaWAK7 at the gene level and provides useful clues for further investigation of how WAKs transmit the extracellular signals to the cytoplasm to activate defense responses.

MRI radiomics predicts the efficacy of EGFR-TKI in EGFR-mutant non-small-cell lung cancer with brain metastasis

AIM

To develop and validate models based on magnetic resonance imaging (MRI) radiomics for predicting the efficacy of epidermal growth factor receptor tyrosine kinase inhibitor (EGFR-TKI) in EGFR-mutant non-small-cell lung cancer (NSCLC) patients with brain metastases.

MATERIALS AND METHODS

117 EGFR-mutant NSCLC patients with brain metastases who received EGFR-TKI treatment were included in this study from January 1, 2014 to December 31, 2021. Patients were randomly divided into training and validation cohorts in a ratio of 2:1. Radiomics features extracted from brain MRI were screened by least absolute shrinkage and selection operator (LASSO) algorithm. Logistic regression analysis and Cox proportional hazard regression analysis were used to screen clinical risk factors. Clinical (C), radiomics (R), and combined (C + R) nomograms were constructed in models predicting short-term efficacy and intracranial progression-free survival (iPFS), respectively. Calibration curves, Harrell's concordance index (C-index), and decision curve analysis (DCA) were used to evaluate the performance of models.

RESULTS

Overall response rate (ORR) was 57.3% and median iPFS was 12.67 months. The C + R nomograms were more effective. In the short-term efficacy model, the C-indexes of C + R nomograms in training cohort and validation cohort were 0.860 (0.820-0.901, 95%CI) and 0.843 (0.783-0.904, 95%CI). In iPFS model, the C-indexes of C + R nomograms in training cohort and validation cohort were 0.837 (0.751-0.923, 95%CI) and 0.850 (0.763-0.937, 95%CI).

CONCLUSION

The C + R nomograms were more effective in predicting EGFR-TKI efficacy of EGFR-mutant NSCLC patients with brain metastases than single clinical or radiomics nomograms.

Introgression of an all-stage and broad-spectrum powdery mildew resistance gene Pm3VS from Dasypyrum villosum chromosome 3V into wheat

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is a serious disease tothat threatens wheat production globally. It is imperative to explore novel resistance genes in order to control this disease throughby developing and planting resistant varieties. Here, we identified a wheat-Dasypyrum villosum 3V (3D) disomic substitution line, NAU3815 (2n=42), with a high level of powdery mildew resistance at both the seedling and adult-plant stages. Subsequently, NAU3815 was used to generate recombination between chromosomes 3V and 3D. Through genomic in situ hybridization (GISH), fluorescence in situ hybridization (FISH)GISH/FISH and 3VS, 3VL-specific markers analysis, four introgression lines were developed from the selfing progenies of 3V and 3D double monosomic line NAU3816, which was derived from the F1 hybrids of NAU3815/NAU0686. There were t3VS (3D) ditelosomic substitution line NAU3817, t3VL (3D) ditelosomic substitution line NAU3818, homozygous T3DL·3VS translocation line NAU3819, and homozygous T3DS·3VL translocation line NAU3820. Powdery mildew tests of these lines confirmed the presence of an all-stage and broad-spectrum powdery mildew resistance gene, Pm3VS, located on chromosome arm 3VS. When compared with the recurrent parent NAU0686 plants, the T3DL·3VS translocation line NAU3819 showed no obvious negative effect on yield-related traits. However, the introduction of the T3DL·3VS translocated chromosome had a strong effect on reducing the flag-leaf length. Consequently, the T3DL·3VS translocation line NAU3819 provides a new germplasm in breeding for both resistance and plant architecture.

Gold(I) ion and the phosphine ligand are necessary for the anti-Toxoplasma gondii activity of auranofin

Auranofin, an FDA-approved drug for rheumatoid arthritis, has emerged as a promising antiparasitic medication in recent years. The gold(I) ion in auranofin is postulated to be responsible for its antiparasitic activity. Notably, aurothiomalate and aurothioglucose also contain gold(I), and, like auranofin, they were previously used to treat rheumatoid arthritis. Whether they have antiparasitic activity remains to be elucidated. Herein, we demonstrated that auranofin and similar derivatives, but not aurothiomalate and aurothioglucose, inhibited the growth of Toxoplasma gondii in vitro. We found that auranofin affected the T. gondii biological cycle (lytic cycle) by inhibiting T. gondii's invasion and triggering its egress from the host cell. However, auranofin could not prevent parasite replication once T. gondii resided within the host. Auranofin treatment induced apoptosis in T. gondii parasites, as demonstrated by its reduced size and elevated phosphatidylserine externalization (PS). Notably, the gold from auranofin enters the cytoplasm of T. gondii, as demonstrated by scanning transmission electron microscopy-energy dispersive X-ray spectroscopy (STEM-EDS) and Inductively Coupled Plasma-Mass Spectrometry (ICP-MS).IMPORTANCEToxoplasmosis, caused by Toxoplasma gondii, is a devastating disease affecting the brain and the eyes, frequently affecting immunocompromised individuals. Approximately 60 million people in the United States are already infected with T. gondii, representing a population at-risk of developing toxoplasmosis. Recent advances in treating cancer, autoimmune diseases, and organ transplants have contributed to this at-risk population's exponential growth. Paradoxically, treatments for toxoplasmosis have remained the same for more than 60 years, relying on medications well-known for their bone marrow toxicity and allergic reactions. Discovering new therapies is a priority, and repurposing FDA-approved drugs is an alternative approach to speed up drug discovery. Herein, we report the effect of auranofin, an FDA-approved drug, on the biological cycle of T. gondii and how both the phosphine ligand and the gold molecule determine the anti-parasitic activity of auranofin and other gold compounds. Our studies would contribute to the pipeline of candidate anti-T. gondii agents.

Wheat Pm55 alleles exhibit distinct interactions with an inhibitor to cause different powdery mildew resistance

Powdery mildew poses a significant threat to wheat crops worldwide, emphasizing the need for durable disease control strategies. The wheat-Dasypyrum villosum T5AL·5 V#4 S and T5DL·5 V#4 S translocation lines carrying powdery mildew resistant gene Pm55 shows developmental-stage and tissue-specific resistance, whereas T5DL·5 V#5 S line carrying Pm5V confers resistance at all stages. Here, we clone Pm55 and Pm5V, and reveal that they are allelic and renamed as Pm55a and Pm55b, respectively. The two Pm55 alleles encode coiled-coil, nucleotide-binding site-leucine-rich repeat (CNL) proteins, conferring broad-spectrum resistance to powdery mildew. However, they interact differently with a linked inhibitor gene, SuPm55 to cause different resistance to wheat powdery mildew. Notably, Pm55 and SuPm55 encode unrelated CNL proteins, and the inactivation of SuPm55 significantly reduces plant fitness. Combining SuPm55/Pm55a and Pm55b in wheat does not result in allele suppression or yield penalty. Our results provide not only insights into the suppression of resistance in wheat, but also a strategy for breeding durable resistance.

Remnant cholesterol is associated with hip BMD and low bone mass in young and middle-aged men: a cross-sectional study

PURPOSE

Remnant cholesterol (RC) is a contributor to cardiovascular diseases, obesity, diabetes, and metabolic syndrome. However, the specific relationship between RC and bone metabolism remains unexplored. Therefore, we aimed to investigate the relationships of RC with hip bone mineral density (BMD) and the risk of low bone mass.

METHODS

Physical examination data was collected from men aged < 60 years as part of the Kailuan Study between 2014 and 2018. The characteristics of the participants were compared between RC quartile groups. A generalized linear regression model was used to evaluate the relationship between RC and hip BMD and a logistic regression model was used to calculate odds ratios (ORs) and 95% confidence intervals (CIs) for low bone mass. Additional analyses were performed after stratification by body mass index (BMI) (≥ or < 24 kg/m2). Sensitivity analyses were performed by excluding individuals who were taking lipid-lowering therapy or had cancer, cardiovascular diseases, or diabetes.

RESULTS

Data from a total of 7,053 participants were included in the analysis. After adjustment for confounding factors, RC negatively correlated with hip BMD (β =  - 0.0079, 95% CI: - 0.0133, - 0.0025). The risk of low bone mass increased from the lowest to the highest RC quartile, with ORs of 1 (reference), 1.09 (95% CI: (0.82, 1.44), 1.35 (95%CI: 1.02, 1.77), and 1.43 (95% CI: 1.09, 1.89) for Q1, Q2, Q3, and Q4, respectively (P for trend = 0.004) in the fully adjusted model. Compared to RC < 0.80 mmol/l group, the risk of low bone mass increased 39% in RC ≥ 0.80 mmol/l group (P < 0.001). The correlation between RC and hip BMD was stronger in participants with BMI ≥ 24 kg/m2 group (β =  - 0.0159, 95% CI: - 0.0289, - 0.0029). The results of sensitivity analyses were consistent with the main results.

CONCLUSION

We have identified a negative correlation between serum RC and hip BMD, and a higher RC concentration was found to be associated with a greater risk of low bone mass in young and middle-aged men.

Comparison of different dosages of propofol combined with its equivalent alfentanil in outpatient abortion: a prospective, double-blinded, randomized trial

OBJECTIVE

This study aimed at determining the optimal dose combination of alfentanil and propofol for outpatient abortion anesthesia.

PATIENTS AND METHODS

The study was separated into two parts. In the first part, patients were to determine the median effective dose (ED50) and the 95% effective dose (ED95) of alfentanil in combination with 2.5 mg·kg-1 propofol to inhibit body movements during the abortion using the Dixon up-and-down sequential allocation method. In the second part, 170 patients were randomly divided into group C (2.0 mg·kg-1 propofol with alfentanil 12.16 μg·kg-1) and group E (2.5 mg·kg-1 propofol with its ED95) to compare the anesthetic effect. The primary outcome was the sedation level during general anesthesia. The secondary outcomes were circulation, respiratory complications, and postoperative recovery quality.

RESULTS

The ED50 and the ED95 values of alfentanil were 3.37 μg·kg-1 (95% CI: 2.58-3.97 μg·kg-1) and 4.68 μg·kg-1 (95% CI: 4.04-9.32 μg·kg-1). The frequency of deep sedation in group E was significantly higher than in group C (76.5% vs. 60%). Patients in group C showed more wakefulness even during the surgery (14.3% vs. 4.4%). The results of our exploratory analyses did not reveal differences in respiratory depression, circulatory depression, postoperative side effects, or recovery outcomes.

CONCLUSIONS

The combination of 2.5 mg·kg-1 propofol and 4.68 μg·kg-1 alfentanil produces a better sedative effect than the combination of 2.0 mg·kg-1 propofol and 12.16 μg·kg-1 alfentanil without increasing additional risks associated with anesthesia.

A new strategy for enhanced electrochemically activation of persulfate on B and Co co-modified carbon felt in flow-through system for efficient organic pollutants degradation

Electrochemical activation of persulfate (EA-PS) is gradually attracting attention as an emerging method for wastewater treatment. In this study, a novelty flow-through EA-PS system was first attempted for pollutants degradation using boron and cobalt co-doping carbon felt (B, Co-CF) as the cathode. SEM images, XPS and XRD spectra of B, Co-CF were investigated. The optimal doping ration between B and Co was 1:2. Increasing current density, PS concentration and flow rate, decreasing initial pH accelerated the removal of AO7. The mechanism involved in EA-PS were the comprehensive effect of DET, •OH and SO4•-. B, Co-CF cathode for flow-through system was stable with five cycles efficient AO7 decay performance. EA-PS in flow-through system was an efficient method with low cost and efficient pollutants degradation. This work provides a feasible strategy for synergistically enhancing PS activation and promoting the degradation of organic pollutants.

Identification and map-based cloning of long glume mutant gene lgm1 in barley

The spikes of gramineous plants are composed of specialized units called spikelets. Two bracts at the spikelet bases are known as glumes. The spikelet glumes in barley are degenerated into threadlike structures. Here, we report a long glume mutant, lgm1, similar in appearance to a lemma with a long awn at the apex. Map-based cloning showed that the mutant lgm1 allele has an approximate 1.27 Mb deletion of in chromosome 2H. The deleted segment contains five putative high-confidence genes, among which HORVU.MOREX.r3.2HG0170820 encodes a C2H2 zinc finger protein, an ortholog of rice NSG1/LRG1 and an important candidate for the Lgm1 allele. Line GA01 with a long glume and short awn was obtained in progenies of crosses involving the lgm1 mutant. Interestingly, lsg1, a mutant with long glumes on lateral spikelets, was obtained in the progenies of the lgm1 mutant. The long glume variant increased the weight of kernels in the lateral spikelets and increased kernel uniformity across the entire spike, greatly improving the potential of six-rowed barley for malting.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01448-x.

A Self-Attention-Based Neural Network for Predicting Immune Checkpoint Inhibitors Response

PURPOSE/OBJECTIVE(S)

Cancer cells evade immune system by negatively regulating T cells via immune checkpoints (e.g., PD-1). By blocking these checkpoints, the ability of immune system to recognize and kill cancer cells restores. Individual response rate of checkpoint blockade varies among patients, with 50%-80% in specific types of cancer such as melanoma, while only 15%-30% in most other tumors. Yet it is still an open question what is the set of biomarkers that are crucial to the response to immune checkpoint inhibitors (ICI). The overall goal of this study is to develop and validate a biologically-aware interpretable deep learning model to identify the biomarkers that can predict the survival outcome to ICI treatment.

MATERIALS/METHODS

The self-attention mechanism could yield interpretable results where important biomarkers may have more "attention". However, in classical self-attention mechanism, the prior biological knowledge of protein interactions (PPI) and gene pathways are not incorporated. In this study, we propose a weighted biologically-aware attention score, where it is weighted against the gene centrality and pathway length. The genes that are closely connected to mutated genes receive 'high attention', while the genes that are far away from mutated genes along the pathway receive "lower attention". We then train, validate and test our model using 1,660 patients of nine types of cancer. To validate the prediction, 1. We evaluate the accuracy via concordance index. 2. We identified the genes that receive high attention and verify their functions in existed literature. 3. We perform sanity check by removing these genes from the data, re-training and predicting again, and comparing the prediction accuracy.

RESULTS

Our framework has achieved an average accuracy (measured via c-index) of 0.60 ± 0.06 for NSCLC and 0.58 ± 0.07 for melanoma, which is superior to both the gold standard COX-PH model (0.57 ± 0.06 for NSCLC and 0.53 ± 0.03 for melanoma) and DeepSurv (0.54 ± 0.05 for NSCLC and 0.51 ± 0.10 for melanoma). Genes that receive high attention have been validated by supporting literature, which provides an additional means of verifying the prediction in comparison to "black box" deep learning models, where there is no way to comprehend the reason behind predictions. Removing the top 8% high-attention genes (∼25 genes) from the data while using the remaining 92% for making predictions resulted in a drop in accuracy to 0.55 ± 0.073 for NSCLC and 0.56 ± 0.03 for melanoma, underscoring the significance of these genes. Patient stratification is also performed by dividing patients into responders and non-responders based on prediction score.

CONCLUSION

In this study, we propose and validate a biologically-aware self-attention based deep learning model which outperforms commonly-used survival models. Additionally, this tool has the potential to identify key biomarkers while assist in clinical decision-making, which demonstrates a promising step for immunotherapy response prediction.

Accelerating Volumetric CT and MRI Imaging by Reference-Free Deep Learning Transformation from Low-Resolution to High-Resolution

PURPOSE/OBJECTIVE(S)

High-resolution (HR) images are important in precision radiation oncology. However, acquiring HR volumetric CT and MRI images is often time consuming; also, the resolution in some direction(s) (e.g., z-direction in the case of CT) is often limited by imaging hardware or fundamental imaging principle. Super-resolution (SR) imaging, i.e., the low-resolution (LR) to HR image transformation, is widely used to improve image resolution. Data-driven deep learning (DL) methods have achieved great success in SR imaging, yet they can hardly be applied to medical imaging as they require large amount of LR-HR image pairs to train the model. We therefore propose a reference-free DL method to increase resolutions of volumetric medical images in an efficient way.

MATERIALS/METHODS

We propose a maximum likelihood estimation (MLE)-based implicit neural representation (INR) network for SR imaging. The INR network aims to represent an image as a continuous function parameterized by a coordinate-based multi-layer perceptron. The INR network takes image coordinates as input and outputs corresponding pixel intensities. To train the network without using any HR images, we use a MLE framework to model LR observations' statistics and their relation to the latent HR image. The predicted HR image from the INR's output is transformed to LR images based on the MLE, and the network parameters are then optimized by minimizing the distance between the transformed LR images and actual LR observations. We demonstrate the efficacy of the proposed method on CT and MRI images for 2x, 4x, and 8x SR using only one or two LR image(s). The performance is compared with a conventional SR method named plain MLE, in terms of visual quality and numerical qualities of PSNR and SSIM.

RESULTS

Our method outperformed the plain MLE method in the experiment. Table 1 reports the numerical improvements of our method over the compared plain MLE method. For 2x SR with a single LR image, our method achieved significant improvements in both PSNR and SSIM. When using two LR images, the better structural restoration capability of our method became more obvious with higher SR magnifications, as indicated by the increased SSIM differences. Better noise suppression capability of our method is observed in all our studies, as indicated by the PSNR values. In visual quality evaluation, we observed sharper image details with less noise in SR images generated by the proposed method, compared with the plain MLE method.

CONCLUSION

The proposed novel reference-free DL method can efficiently provide high-quality HR images with only one or two LR images for CT and MRI imaging. This method can be easily generalized to many other radiation therapy related applications without the requirement for HR reference images.

Reinforcement Learning Powered Station Parameter Optimized Radiation Therapy (SPORT): A Novel Treatment Planning and Beam Delivery Technique

PURPOSE/OBJECTIVE(S)

Conventional intensity modulated radiation therapy (IMRT) with a typical 5-20 fixed beams often does not provide sufficient angular sampling required for conformal dose shaping, whereas current volumetric modulated arc therapy (VMAT) discretizes the angular space into equally spaced control points without considering the differential need for intensity modulation of different angles, leading to undersampling at some angles while oversampling at some other angles. Our goal is to develop a node or station parameter optimized radiation therapy (SPORT) strategy with simultaneously optimized angular sampling and beam modulation by leveraging state-of-the-art reinforcement learning and the unique capability of modern digital LINACs in dose delivery through programmable nodal points.

MATERIALS/METHODS

We developed a SPORT optimization framework, in which, the process of programming control points (or station parameters) was formulated as a stochastic dynamic programming problem, which was solved by a reinforcement learning-based algorithm. On-policy reinforcement learning method, namely, state-action-reward-state-action (SARSA) was integrated with deep convolutional neural network to predict station parameters by utilizing the patient's anatomical structures meanwhile considering the delivery capability of a typical digital LINAC machine. Here, the deep convolutional neural network estimated the state-action value by using the quality of the plan with current station parameters when a next potential station parameter was selected. The state-action value was then updated by SARSA learning. The quality of the plan was quantified by dosimetry constraints. The model was assessed by a retrospective study on a cohort of patients underwent head-and-neck radiation therapy. Dosimetric analysis and delivery efficiency comparisons were used to evaluate the performance of the proposed framework.

RESULTS

Our model was used to generate 16 plans unseen in the original training set. All the plans predicted by our model achieved better dose distributions without violating clinical planning constraints. Moreover, instead of using 4 full standard arcs in the original clinically used plans obtained via manual optimization, the predicted plans only used one full standard arc (about 178 control points) plus boost from a few sub-arcs (less than 30 degrees of gantry angles), which significantly improved the efficiency of the beam delivery. We are in the process of integrating the sub-arcs into the full arc by considering the programmable capability of modern LINACs.

CONCLUSION

We demonstrated that a machine learning-based SPORT framework capable of optimizing the spatial sampling and beam modulation simultaneously for modern radiation therapy. The framework not only significantly improves the quality and efficiency of beam delivery, but also has the potential to be incorporated into current clinical workflow to improve the efficiency of dose planning and delivery.

Automated Small Tumor Segmentation by a Template-Based Global Hierarchical Attention Method

PURPOSE/OBJECTIVE(S)

Accurate segmentation of tumors is significant for radiation therapy treatment planning and clinical decision-making. While deep convolutional neural network-based methods have found valuable applications in automatic medical segmentation, tumor segmentation, especially small tumor segmentation, remains challenging due to deficiencies of current deep learning in convolutional and pooling operations, which often results in the loss of small object information. This research proposes a global hierarchical attention-based method for accurate and automated segmentation of small tumors by exploiting the associations between small tumors and the feature maps of large tumors.

MATERIALS/METHODS

This study included 131 patients with liver cancer. The in-plane resolution of the patients' CTs is from 0.55 mm to 1.0 mm and slice spacing from 0.45 mm to 6.0 mm. We randomly selected 100 CT scans as the training set and others as the testing set. Each CT slice of the testing set was separated into groups according to tumor size as follows: 0.1-2.0, 2.1-5.0, 5.1-10.0, and 10.1-20.0 cm. The CT slice without tumor or tumor size > 20 cm were excluded. This work presents a tumor template-based hierarchical attention method to quantify the relation between small and large tumors by computing their feature maps. The relation of small-large tumors can compensate for the information loss of small tumors during the convolutional and pooling operations and improve the performance of small tumor segmentation.

RESULTS

Among 20,693 CT slices of the 31 testing patients, 3.0% CT slices with tumors ≤2 cm, 6.7% ≤5 cm, 10.6% ≤10 cm, and 13.4%≤20 cm. We compared our method with six widely used segmentation models. The results show our model outperforms other methods on all sizes of liver tumors, especially for small size tumors: For the 0.1-2.0 cm liver tumor, it achieved 8.4%, 10.0%, 11.3%, 9.1%, 10.9%, and 9.6% improvement compared to Unet, PAN, DeepLabV3, FPN, LinkNet, and PSPNet, respectively.

CONCLUSION

We found that the small-large tumors relation can significantly improve small tumor segmentation, which is valuable for treatment planning, and clinical decision-making. Our experimental results show that our method can significantly improve the accuracy of segmenting small liver tumors compared to existing deep-learning-based models. The method is quite general and can be extended to other types of tumor detection and segmentation. We discovered that the relationship between small and large tumors can significantly enhance the segmentation of small tumors, which has significant value for treatment planning and clinical decision-making. Our experiments demonstrate that our approach significantly improves the accuracy of small liver tumor segmentation compared to existing deep learning-based models. Our method is quite versatile and can be extended to other types of tumor detection and segmentation.

An Efficient Transformer Model for Synthesizing Dual Energy CT from Single Energy Scanner

PURPOSE/OBJECTIVE(S)

Dual-energy CT can be used to optimize radiation treatment. Recently, deep learning has been demonstrated to synthesize high-energy CT images from low-energy ones for dose reduction and lower CT system burden. As the state-of-the-art deep learning architecture, the computation burden of Transformer increases quadratically with the feature size, making the model training resource-demanding or even infeasible. Here, we introduce an efficient transformer for the balance between CT image synthesis quality and computational burden.

MATERIALS/METHODS

The model is a U-shape deep neural network with encoders and decoders built by Transformer blocks. The model input is low-energy 100kVp CT image and the output is high-energy 140kVp one. Each block has a Self Channel Correlation Unit (SCCU) and a Self Spatial Attention Unit (SSAU). Local shortcuts are applied for both units. Under-sampling operation achieved by pixel shuffling is used to obtain multi-scale feature maps, and the transformer block is applied on each feature scale. Symmetric skip connection sending features from shallow layers to deep layers, thus an additional 1 × 1 convolution is used for feature fusion in each decoder. In a SCCU, the feature is first mapped to one Query, one Key, and one Value. Then the Query and the Key tensors perform matrix multiplication to compute cross covariance of feature channels. The channel correlation score can be obtained by normalization of the covariance, and it is used to weight the Value tensor. As a result, the model complexity only increases linearly with the feature size. Besides the channel weighting, we enhance spatial information using SSAU, where the feature is mapped to two tensors. One tensor after activation is used to point-wisely calibrate another tensor. Additional Transformer blocks are cascaded to the last decoder for feature refinement. Because of the structure similarity of low- and high-energy CT images, a global shortcut is used to ease model training. Clinical iodine contrast-enhanced dual energy CT image datasets of 19 patients are used in this study. The dual-energy scanning is performed by a SOMATOM Definition Flash DECT scanner. We split the datasets into training dataset of 15 patients, validation dataset of 1 patient, and testing dataset of 3 patients. The image size is 512 × 512 with pixel size 0.5 × 0.5 mm2.

RESULTS

The U-Net model with 1.95M parameters and 44.87G FLOPS achieved the averaged PSNR value of 44.55 dB (s.t.d. 1.34) and averaged RMSE value of 0.0060 (s.t.d. 0.001). In comparison, our efficient Transformer with 1.408M parameters and 31.375G FLOPS achieved the averaged PSNR value of 44.78 dB (s.t.d. 1.37) and RMSE value of 0.0059 (s.t.d. 0.001), demonstrating our model has better performance with small model size and less computation.

CONCLUSION

The efficient Transformer model allows high-resolution CT image synthesis with small model scale and computation burden from low-energy CT image.

The Role of CT-Based Radiomics Nomogram in Differential Diagnosis of Immune Checkpoint Inhibitor-Related Pneumonitis from Radiation Pneumonitis for Patients with ESCC

PURPOSE/OBJECTIVE(S)

The combination of immunotherapy and chemoradiotherapy has widely used for patients with esophageal squamous cell carcinoma (ESCC) and induced treatment-related adverse effects, particularly immune checkpoint inhibitor-related pneumonitis (CIP) and radiation pneumonitis (RP). The aim of this study is to differentiate between CIP and RP by the CT radiomics and clinical or radiological parameters.

MATERIALS/METHODS

A total of 76 ESCC patients with pneumonitis were enrolled in this retrospective study and divided into training dataset (n = 53) and validation dataset (n = 23). A total of 837 radiomics features were extracted from regions of interest (ROIs) based on the lung parenchyma window of CT images. A radiomics signature was constructed on the basis of the predictive features by the least absolute shrinkage and selection operator (LASSO). A logistic regression was applied to develop radiomics nomogram. Receiver operating characteristics (ROC) curve and area under the curve (AUC) were applied to evaluate the performance of pneumonitis etiology identification.

RESULTS

No significant difference was detected between training dataset and validation dataset. The radiomics signature which was made up of four radiomics features shown a favorable performance on differentiating between CIP and RP with the α-binormal-based and empirical AUC = 0.831 and 0.843. Patients with RP had a close relationship with location (p = 0.003) and shape of lesions (p = 0.002). The nomogram that combined with radiomics signature and clinical factors improved the classifying performance on discrimination in the training dataset (AUCαbin = 0.963 and AUCemp = 0.964). The results were verified in the validation dataset with AUC = 0.967 and 0.964.

CONCLUSION

CT-based radiomics features have potential values for differentiating between patients with CIP and RP. Addition of bilateral changes and sharp border produced superior model performance on classifying, which could be a useful method to improve related clinical decision-making.

Deep Learning-Based Single-View Fluorescence Dose Reconstruction for 3D Dosimetry

PURPOSE/OBJECTIVE(S)

3D dose distribution measurement is crucial for precise radiotherapy. Radiation-excited fluorescence imaging has potential for the 3D dosimetry with high spatial resolution, but multiple fluorescence images from different view-angles are required for analytical reconstruction techniques. Furthermore, the imaging data are contaminated by anisotropic Cherenkov light emission and statistical noise. This project aims to establish a novel deep learning-based model to predict 3D dose distributions from a single-view 2D fluorescence image while simultaneously removing the adverse effects of Cherenkov signals and other noises.

MATERIALS/METHODS

A total of 124 single-aperture static photon beams were delivered to an acrylic tank containing 1 g/L quinine hemisulfate water solution with varying aperture shapes and collimator angle. The emitted optical signals were detected by a low-cost CMOS camera for 20 seconds, and image pre-processing was performed to obtain input 2D fluorescence images with 0.3 × 0.3 mm spatial resolution. 3D back-projected dose distribution images were also calculated from the input fluorescence images. Ground-truth of 3D dose distributions and 2D field map images were obtained from a clinical treatment planning system with 1.4 × 1.4 × 1.4 mm spatial resolution. The proposed deep learning-based dose reconstruction method involved 3 steps. First, 2D fluence map images at the bottom plane of the tank were predicted from the fluorescence images by using a customized convolutional neural network (CNN). Second, the predicted fluence map images were transformed into the 2D field map images on the isocenter plane by applying perspective transformation. Finally, 2D dose distributions at a given radiological depth were calculated by using the predicted field map images, the back-projected dose distribution images, and the radiological depth value as inputs of a shallow CNN. Both CNN models were trained separately, and the 3D dose distributions were predicted by concatenating the output 2D dose distributions at various radiological depths.

RESULTS

The proposed CNN model yielded accurate 2D field map images. Averaged Dice similarity coefficient and mean absolute error of the field maps in the test data was 92.0% ± 4.6% and 0.0132 ± 0.0113, respectively. Moreover, our deep learning-based approach was able to predict accurate 3D dose distributions from the 2D fluorescence images. Mean squared error and averaged 3D gamma passing ratio (3%/3mm) were 9.55 mGy ± 6.8 mGy and 86.3% ± 9.86%, respectively.

CONCLUSION

Theproposed deep learning-based method calculated accurate 3D dose distributions from a single-view 2D fluorescence image. Since this technique require only a single CMOS camera image and fluorescent material, it can be readily used for any external radiotherapy modalities, including SRS/SBRT with small fields. This method is useful for acquiring 3D dose distribution data for precise dose verification within a few seconds.

Clinical Implementation of an Automated IMRT/VMAT Treatment Planning Tool

PURPOSE/OBJECTIVE(S)

To create an in-house automated treatment planning tool for IMRT/VMAT treatments and evaluate the dosimetric plan quality against manually generated plans.

MATERIALS/METHODS

A scripting application programming interface is employed to interact with a commercial treatment planning system (TPS) to implement automatic plan evaluation and update optimization parameters by mimicking the human planning process. The automated planning performs in an iterative fashion until reaching an acceptable tradeoff among target coverage/dose homogeneity and sparing of critical organs at risk. In each iteration, the dose constraints, priorities, and optimization structures for are automatically updated based on the results of the current iteration. Twenty previously treated plans (10 prostate and 10 head and neck), were preliminarily used to evaluate the performance of the automated planning tool. The differences in target and organ-at-risk metrics from the manually generated clinical plans were analyzed using paired t-test to evaluate clinical acceptability of tour automated planning tool. The current in-house-developed automated planning solution is able to create plans for different disease sites, including head & neck, prostate, pelvis, and lung. So far, the VMAT plans for more than 150 different cases have been generated with the tool. The results for these were also evaluated.

RESULTS

Compared to the manually generated clinical head and neck plans, all auto plans achieved PTV D95% coverage and critical organs at risk sparing without statistically significant change in average global Dmax (107.4% for manual vs 107.3% for automated plans). The auto-planning solution provided reduced maximum doses to brainstem and spinal cord (average reductions with standard deviations of 5.1 ± 2.6 Gy and 2.9 ± 1.4 Gy, respectively, all p <0.03), reduced average mean doses to contralateral parotid, ipsilateral parotid, contralateral submandibular gland, pharynx, esophagus, cochleae (reductions of 2.2 ± 2.9 Gy, 4.8 ± 4.7 Gy, 3.6 ± 5.2 Gy, 2.0 ± 7.1 Gy, 3.9 ± 2.6 Gy, 3.8 ± 5.0 Gy, respectively, all p < 0.045). Similar results were observed for the prostate plans. With the same PTV coverage and without statistically significant change in average global Dmax (106.5% for manual vs 106.8% for automated plans), the automated solution provided superior sparing for both bladder and rectum. Bladder V75, V70, V65 were reduced by 0.6% ± 2.1%, 0.8% ± 2.5%, and 0.9% ± 2.9% (all p <0.04), respectively. Rectum V75, V70, V65, V60 were reduced by 1.0% ± 2.3%, 1.2% ± 2.8%, 1.3% ± 3.2%, 1.6% ± 3.6% (all p < 0.01), respectively.

CONCLUSION

Our automated treatment planning solution is capable of efficiently generating VMAT plans for different disease sites with superior dosimetric indices compared to manually generated plans. Our tool is integrated within a commercial TPS platform, so it has the advantage of seamless adoption into the standard workflow to improve plan quality and treatment planning efficiency in our clinic.

Patient-Adaptive Automated Segmentation in Daily kVCT Images for Radiotherapy of Head and Neck and Prostate Cancer

PURPOSE/OBJECTIVE(S)

The purpose of this study was to examine the use of transfer learning in deep learning-based auto-segmentation of daily kilovoltage computed tomography (kVCT) images for patient-specific adaptive radiotherapy. Using data from the first cohort of patients treated with the innovative BgRT system, the objective of this study was to evaluate the potential benefits of this approach in facilitating efficient and effective adaptive radiotherapy.

MATERIALS/METHODS

For the head and neck (HaN) site and pelvic site, we first trained a deep convolutional segmentation network using a population dataset, consisting of 67 and 56 patient cases, respectively. This population network was then fine-tuned for a specific patient using a transfer learning approach, adapting the network weights. The auto-segmentation network utilized in this study was a 23-layer U-Net with batch normalization, a dropout rate of 0.5, and four skip connections between the encoder and decoder at different levels. We used initial planning CT and 5-26 sets of daily kVCT scans with a total of 8,039 images for patient-specific learning in the 6 HaN cases and 4 pelvic cases, particularly analyzing the relationship between the number of sequential patient-specific training data and the performance of the auto-segmentation. We compared the performance of the patient-specific network with the population network and the clinical rigid registration method, using the Dice similarity coefficient (DSC) as the evaluation metric. Additionally, we investigated the corresponding dosimetric impacts of the different auto-segmentation and registration methods.

RESULTS

The patient-specific network showed improved mean DSC scores of 0.88 and 0.90 for three HaN organs at risk (OARs) and eight pelvic targets and OARs, respectively, compared to the population network (0.70 and 0.63) and the registration method (0.72 and 0.72). The DSC of the patient-specific network steadily improved as the number of longitudinal training cases increased, reaching near saturation after 6 training cases. The use of the patient-specific auto-segmentation resulted in a reduction of the mean discrepancy in target and OAR doses between delivery and planning from 5.5% with the clinical rigid registration to 1.1%.

CONCLUSION

The use of patient-specific transfer learning in auto-segmenting kVCT images showed higher accuracy compared to a conventional population network and clinical registration-based method. This approach holds promise for enhancing dose evaluation accuracy in adaptive radiotherapy.

Real Time Volumetric MRI for MR-Guided 3D Motion Tracking via Sparse Prior-Augmented Neural Representation Learning

PURPOSE/OBJECTIVE(S)

To reconstruct volumetric MRI from orthogonal cine acquisition aided by sparse priors of 2 static 3D MRI through implicit neural representation (NeRP) learning, with the goal of eliminating large-scale training datasets for data-driven sparse MRI reconstruction and supporting clinical workflow of real time 3D motion tracking during MR-guided radiotherapy.

MATERIALS/METHODS

A multi-layer perceptron network was trained to learn the NeRP of a patient-specific MRI dataset, where the network takes 4D data coordinates of voxel locations and motion states as inputs and outputs corresponding voxel intensities. By first learning the NeRP of 2 static 3D MRI with different breathing motion states, prior knowledge of patient breathing motion was embedded into network weights through optimization. The prior knowledge was then augmented from 2 to 31 motion states by querying the optimized network at interpolated/extrapolated motion state coordinates. Starting from the prior-augmented network as an initialization point, the network was further trained using sparse samples of 2 orthogonal cine slices. The final volumetric reconstruction was obtained by querying the trained network at desired 3D spatial locations. We evaluated the proposed method using 5-minute volumetric MRI time series with 340 ms temporal resolution collected from 7 liver carcinoma patients. The time series was acquired using golden-angle radial MRI sequence and reconstructed through retrospective sorting. Two MRI with inhale and exhale states respectively were selected from the first 30 sec of the time series for prior embedding and augmentation. The remaining 4.5-min time series was used for volumetric reconstruction evaluation, where we retrospectively subsampled each MRI to 2 orthogonal slices and compared network-reconstructed images to ground truth images in terms of image quality and the capability of supporting 3D target motion tracking.

RESULTS

Across the 7 patients evaluated, the peak signal to noise ratio between model reconstruction and ground truth was 54.66 ± 6.16 dB and the structural similarity index measure was 0.99 ± 0.01. Gross tumor volume (GTV) contours estimated by deforming a reference state MRI to model-reconstructed and ground truth MRI showed good consistency. The 95-percentile Hausdorff distance between GTV contours was 1.89 ± 1.13 mm, which is less than the voxel dimension. The mean GTV centroid position difference between ground truth and model estimation was less than 1 mm in all 3 orthogonal directions.

CONCLUSION

Volumetric MRI from orthogonal cine acquisition with sparse priors is feasible by modeling prior knowledge through implicit neural representation learning. The model-reconstructed images showed sufficient accuracy in supporting 3D motion tracking of abdominal targets. By eliminating the need for large scale training datasets, the method promises to enable clinical implementation of 3D motion tracking for precision radiation therapy.

Analyzing the Predictive Effects of PD-L1 Expression, Early Changes of bTMB and Circulated CD8+T Cells during Treatment for Responses of RT Combined with ICI in NSCLC

PURPOSE/OBJECTIVE(S)

The beneficial role of immunotherapy and the clinical relevance of current biomarkers remain inconclusive; thus, appropriate strategies and reliable predictors need further definition. A rational combination of biomarkers is needed. Here, we estimated potential predictive factors for responses of radiotherapy (RT) combined with immune checkpoint inhibitor (ICI) in a phase II trial to determine the efficacy and safety of combination of moderate hypofractionated RT with ICI in patients with oligometastatic NSCLC (NCT03557411).

MATERIALS/METHODS

Pretreatment tumor tissue samples and longitudinal blood were collected for immune and tumor biomarker analysis. We examined pre-treatment (pre-ICI) PD-L1 expression in tumor cells. Circulating tumor cell (CTC), PD-L1+CTC, blood tumor mutation burden (bTMB), CD8+T cells, CD4+T cells, NK cells, B cells in circulation were acquired pre-ICI and 1 month after ICI starting (1-mth). In addition, early changes of CTC (CTC), PD-L1+CTC (PD-L1+CTC), bTMB (bTMB), CD8+T cells (CD8+T cells), CD4+T cells (CD4+T cells), NK cells (NK cells), B cells (B cells) were also analyzed to estimate the predictive effects for treatment.

RESULTS

High pre-ICI bTMB and increased CD8+T cells at 1 month was associated with better PFS (p = 0.016; p = 0.006). Interaction analyses revealed that each combination of two markers in the 5 markers including PD-L1, pre-ICI bTMB, 1-mth bTMB, 1-mth CD8+T cells and CD8+T cells was significantly associated with PFS, except for CTC, PD-L1+CTC, CD4+T cells, NK cells and B cells in circulation due to low power. Unsupervised cluster analysis based on these markers revealed three sub-cohorts. Cohort-1 was overrepresented by patients with progressive disease (81%) of whom were negative for 3-4 of the 5 biomarkers. Cohort-3 was overrepresented by patients with partial response (70%) of whom were positive for 3-4 of the 5 biomarkers. Survival analyses of the 3 cohorts indicated a significant association with PFS (p = 0.017).

CONCLUSION

This study suggests that a combination of PD-L1 expression, early changes of bTMB and circulated CD8+T cells as a better predictive biomarker for response to RT combined with ICI. Consequently, refinement of this set of biomarkers and validation in a larger set of patients is warranted.

Does Epigallocatechin Gallate as a Radiation Protective Agent Reduce the Anti-Tumor Effect of Radiotherapy in Postoperative Breast Cancer Radiotherapy?

PURPOSE/OBJECTIVE(S)

Based on the previous encouraging results, we further explored whether EGCG would have a protective effect on potential tumor lesions, that is, reduce the efficacy of radiotherapy. We selected patients with stage III breast cancer with or without EGCG. The local control rate, distant metastasis rate, DFS and OS were compared between the two groups.

MATERIALS/METHODS

Patients with stage III breast cancer who were treated with EGCG and radiotherapy was selected from a phase II clinical study (ClinicalTrials.gov, No. NCT02580279). Each patient was matched with one control patient without EGCG From the medical database of our hospital matching for age and stage. The control group of stage-and age-matched patients was selected at random from the medical database of our hospital RESULTS: We identified 43 EGCG patients and 43 matched controls. The median age was 45 years (range: 26-67). Between the two groups, there were no obvious differences in the baseline demographic or clinical features. When compared to the placebo group, the mean radiation-induced dermatitis index (RIDI) in the EGCG group was substantially lower (2.56±1.14 vs 3.36±1.16 T = -3.232, P = 0.002). Repeated measures ANOVA indicated the significant differences in the RTOG score during the course of radiotherapy between the two groups (F = 9.611 P = 0.003). The patients mostly experienced RID two or three weeks after starting radiotherapy, although in the EGCG group, it appeared later (3.19±0.91 weeks) than it did in the placebo group (2.67±0.84 weeks), P = 0.008. The median follow-up for patients in the EGCG group at the time of data collection was 50.6 months with 95% confidence interval (CI) from 43.9 to 57.3. While it was 48.6 months (95% CI: 43.4-53.9) for patients in the control group. There was no significant difference in overall survival (OS), disease free-survival (DFS) and freedom from locoregional (LRF) and distant failure (DMF) (P > 0.05). At the data cut-off (December 2021), the 4-year DFS with EGCG was 71.4% compared to 65.4% with conventional therapy, and the 4-year OS was 77.0% compared to 80.3%.

CONCLUSION

The prophylactic use of EGCG solution reduced the RID score of stage III breast cancer patients without negatively impacting the therapeutic effect of radiotherapy on the tumor. EGCG is safe and feasible choice for RID for breast cancer during radiotherapy.

Prospective Clinical Adoption of Artificial Intelligence for Organ Contouring in Head and Neck Radiation Treatment Planning

PURPOSE/OBJECTIVE(S)

Patients that undergo head and neck (H&N) radiation therapy (RT) require laborious delineation of organs-at-risk (OARs) on computed tomography (CT) scans in a treatment planning system (TPS) to minimize radiation to normal tissue. This task can be completed rapidly and accurately with recently developed artificial intelligence-based semantic segmentation models. The current study aims to deploy and evaluate a strategy for improving clinical practice with this technology.

MATERIALS/METHODS

Deep learning models were trained and tested with CT scans and OAR contours from previous H&N RT cases at our clinic. Two medical physicists vetted the models and selected a 2.5D U-Net for further implementation. The model was embedded in a dedicated server at the hospital, programmed to read H&N CT scans staged for import into the TPS, generate auto-contours, and write them into a TPS-compatible format made available alongside the scan. In the pilot implementation, the auto-contouring service was utilized for more than 60 cases, prospectively. The auto-contours were quantitatively evaluated against the treatment-approved contours to determine how much modification was performed by the clinical team.

RESULTS

The 2.5D U-Net selected for clinical integration segments 21 OARs in less than 3 minutes per scan. Across all the prospective cases, the mean Dice score and mean 95th percentile Hausdorff distance (mm) between the auto-contour and treatment-approved contour for each of the 21 OARs were as follows, respectively: brainstem (0.93, 1.94), optic chiasm (0.70, 2.96), left cochlea (0.69, 2.37), right cochlea (0.68, 2.44), esophagus (0.88, 2.46), left globe (0.93, 1.50), right globe (0.93, 1.63), glottis (0.91, 2.13), larynx (0.93, 2.76), mandible (0.90, 4.86), left optic nerve (0.78, 1.64), right optic nerve (0.82, 1.65), oral cavity (0.86, 8.46), left parotid gland (0.91, 2.78), right parotid gland (0.91, 2.39), pharynx (0.85, 2.39), spinal cord (0.87, 2.27), left submandibular gland (0.85, 3.46), right submandibular gland (0.83, 3.69), left temporal lobe (0.94, 2.20), and right temporal lobe (0.95, 2.09). The auto-contours for the optic chiasm, optic nerves, cochleas, and submandibular glands differed substantially from the final contours, a finding corroborated by the clinical team; the rest were clinically acceptable with minor or no edits necessary.

CONCLUSION

The proposed strategy provides a sophisticated starting point for treatment planning that has garnered overall favorable feedback from the participating radiation oncologists and dosimetrists. Consequently, the technique is being extended to other treatment sites.

Mechanisms of circRNA/lncRNA-miRNA interactions and applications in disease and drug research

In recent years, breakthroughs in bioinformatics have been made with the discovery of many functionally significant non-coding RNAs (ncRNAs). The discovery of these ncRNAs has further demonstrated the multi-level characteristics of intracellular gene expression regulation, which plays an important role in assisting diagnosis, guiding clinical drug use and determining prognosis in the treatment process of various diseases. microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) are the three major types of ncRNAs that interact with each other. Studies have shown that lncRNAs and circRNAs can sponge miRNAs, thereby influencing normal physiological processes and regulating mRNA expression and, thus, the physiological state of cells. This paper summarizes the mechanism of action and research progress of the three ncRNA and seven types of modalities. This summary is intended to provide new ideas for diagnosing and treating diseases and researching and developing new drugs.

Identification and transfer of a new Pm21 haplotype with high genetic diversity and a special molecular resistance mechanism

A new functional Pm21 haplotype, Pm21(8#), was cloned from the new wheat-H. villosa translocation line T6VS(8#)·6DL, which confers the same strong resistance to powdery mildew through a different resistance mechanism. Broad-spectrum disease resistance genes are desirable in crop breeding for conferring stable, durable resistance in field production. Pm21(4#) is a gene introduced from wild Haynaldia villosa into wheat that confers broad-spectrum resistance to wheat powdery mildew and has been widely used in wheat production for approximately 30 years. The discovery and transfer of new functional haplotypes of Pm21 into wheat will expand its genetic diversity in production and avoid the breakdown of resistance conferred by a single gene on a large scale. Pm21(4#) previously found from T6VS(4#)·6AL has been cloned. In this study, a new wheat-H. villosa translocation, T6VS(8#)·6DL, was identified. A new functional Pm21 haplotype, designated Pm21(8#), was cloned and characterized. The genomic structures and the splicing patterns of Pm21(4#) and Pm21(8#) were different, and widespread sequence diversity was observed in the gene coding region and the promoter region. In the field, Pm21(8#) conferred resistance to Blumeria graminis f. sp. tritici (Bgt), similar to Pm21(4#), indicating that Pm21(8#) was also a resistance gene. However, Bgt development during the infection stage was obviously different between Pm21(4#)- and Pm21(8#)-containing materials under the microscopic observation. Pm21(4#) inhibited the formation of haustoria and the development of hyphae in the initial infection stage, while Pm21(8#) limited the growth of hyphae and inhibited the formation of conidiophores in the late infection stage. Therefore, Pm21(8#) is a new functional Pm21 haplotype that provides a new gene resource for wheat breeding.

Particulate suspension coating of capillary tubes

The displacement of a suspension of particles by an immiscible fluid in a capillary tube or in porous media is a canonical configuration that finds application in a large number of natural and industrial applications, including water purification, dispersion of colloids and microplastics, coating and functionalization of tubings. The influence of particles dispersed in the fluid on the interfacial dynamics and on the properties of the liquid film left behind remain poorly understood. Here, we study the deposition of a coating film on the walls of a capillary tube induced by the translation of a suspension plug pushed by air. We identify the different deposition regimes as a function of the translation speed of the plug, the particle size, and the volume fraction of the suspension. The thickness of the coating film is characterized, and we show that similarly to dip coating, three coating regimes are observed, liquid only, heterogeneous, and thick films. We also show that, at first order, the thickness of films thicker than the particle diameter can be predicted using the effective viscosity of the suspension. Nevertheless, we also report that for large particles and concentrated suspensions, a shear-induced migration mechanism leads to local variations in volume fraction and modifies the deposited film thickness and composition.

Fine mapping of powdery mildew and stripe rust resistance genes Pm5V/Yr5V transferred from Dasypyrum villosum into wheat without yield penalty

The novel wheat powdery mildew and stripe rust resistance genes Pm5V/Yr5V are introgressed from Dasypyrum villosum and fine mapped to a narrowed region in 5VS, and their effects on yield-related traits were characterized. The powdery mildew and stripe rust seriously threaten wheat production worldwide. Dasypyrum villosum (2n = 2x = 14, VV), a relative of wheat, is a valuable resource of resistance genes for wheat improvement. Here, we describe a platform for rapid introgression of the resistance genes from D. villosum into the wheat D genome. A complete set of new wheat-D. villosum V (D) disomic substitution lines and 11 D/V Robertsonian translocation lines are developed and characterized by molecular cytogenetic method. A new T5DL·5V#5S line NAU1908 shows resistance to both powdery mildew and stripe rust, and the resistances associated with 5VS are confirmed to be conferred by seedling resistance gene Pm5V and adult-plant resistance gene Yr5V, respectively. We flow-sort chromosome arm 5VS and sequence it using the Illumina NovaSeq 6000 system that allows us to generate 5VS-specific markers for genetic mapping of Pm5V/Yr5V. Fine mapping shows that Pm5V and Yr5V are closely linked and the location is narrowed to an approximately 0.9 Mb region referencing the sequence of Chinese Spring 5DS. In this region, a NLR gene in scaffold 24,874 of 5VS orthologous to TraesCS5D02G044300 is the most likely candidate gene for Pm5V. Soft- and hard-grained T5DL·5V#5S introgressions confer resistance to both powdery mildew and stripe rust in diverse wheat genetic backgrounds without yield penalty. Meanwhile, significant decrease in plant height and increase in yield were observed in NIL-5DL·5V#5S compared with that in NIL-5DL·5DS. These results indicate that Pm5V/Yr5V lines might have the potential value to facilitate wheat breeding for disease resistance.

Aged Callus Skeletal Stem/Progenitor Cells Contain an Inflammatory Osteogenic Population With Increased IRF and NF-κB Pathways and Reduced Osteogenic Potential

Skeletal stem/progenitor cells (SSPCs) are critical for fracture repair by providing osteo-chondro precursors in the callus, which is impaired in aging. However, the molecular signatures of callus SSPCs during aging are not known. Herein, we performed single-cell RNA sequencing on 11,957 CD45-CD31-Ter119- SSPCs isolated from young and aged mouse calluses. Combining unsupervised clustering, putative makers, and DEGs/pathway analyses, major SSPC clusters were annotated as osteogenic, proliferating, and adipogenic populations. The proliferating cluster had a differentiating potential into osteogenic and adipogenic lineages by trajectory analysis. The osteoblastic/adipogenic/proliferating potential of individual clusters was further evidenced by elevated expression of genes related to osteoblasts, adipocytes, or proliferation. The osteogenic cluster was sub-clustered into house-keeping and inflammatory osteogenic populations that were decreased and increased in aged callus, respectively. The majority of master regulators for the inflammatory osteogenic population belong to IRF and NF-κB families, which was confirmed by immunostaining, RT-qPCR, and Western blot analysis. Furthermore, cells in the inflammatory osteogenic sub-cluster had reduced osteoblast differentiation capacity. In conclusion, we identified 3 major clusters in callus SSPCs, confirming their heterogeneity and, importantly, increased IRF/NF-κB-mediated inflammatory osteogenic population with decreased osteogenic potential in aged cells.

A Survey on the Willingness of Ganzhou Residents to Participate in "Internet + Nursing Services" and Associated Factors

Objective

To investigate the willingness of Ganzhou residents to participate in “Internet + Nursing services” and analyse the relevant influencing factors.

Methods

From May to June 2021, 426 Ganzhou residents were surveyed using an Internet + Nursing services questionnaire and the relevant influencing factors were analysed. The questionnaire comprised two parts: demographic characteristic section and a questionnaire on residents’ willingness to participate in Internet + Nursing services including for dimensions (awareness, participation, trust and need), a 5-point Likert scale was used.

Results

A total of 397 valid questionnaires were recovered, and the total willingness of Ganzhou residents to participate in the service was derived as 11.59 ± 2.14. The results of multiple linear regression analyses showed that the presence of family members with a chronic disease or mobility difficulties, and an awareness and trust of Internet + Nursing services were influencing factors of residents’ participation willingness (P < 0.05).

Conclusion

The participation willingness of Ganzhou residents in Internet + Nursing services is modestly low, and the reasons for participation varied. It is suggested that the government and pilot hospitals strengthen the publicity surrounding these services, improve safety measures, strengthen team training, and develop products suitable for the elderly to increase residents’ participation willingness.

Genome-wide identification of the NLR gene family in Haynaldia villosa by SMRT-RenSeq

Background

Nucleotide-binding and leucine-rich repeat (NLR) genes have attracted wide attention due to their crucial role in protecting plants from pathogens. SMRT-RenSeq, combining PacBio sequencing after resistance gene enrichment sequencing (RenSeq), is a powerful method for selectively capturing and sequencing full-length NLRs. Haynaldia villosa, a wild grass species with a proven potential for wheat improvement, confers resistance to multiple diseases. So, genome-wide identification of the NLR gene family in Haynaldia villosa by SMRT-RenSeq can facilitate disease resistance genes exploration.

Results

In this study, SMRT-RenSeq was performed to identify the genome-wide NLR complement of H. villosa. In total, 1320 NLRs were annotated in 1169 contigs, including 772 complete NLRs. All the complete NLRs were phylogenetically analyzed and 11 main clades with special characteristics were derived. NLRs could be captured with high efficiency when aligned with cloned R genes, and cluster expansion in some specific gene loci was observed. The physical location of NLRs to individual chromosomes in H. villosa showed a perfect homoeologous relationship with group 1, 2, 3, 5 and 6 of other Triticeae species, however, NLRs physically located on 4VL were largely in silico predicted to be located on the homoeologous group 7. Fifteen types of integrated domains (IDs) were integrated in 52 NLRs, and Kelch and B3 NLR-IDs were found to have expanded in H. villosa, while DUF948, NAM-associated and PRT_C were detected as unique integrated domains implying the new emergence of NLR-IDs after H. villosa diverged from other species.

Conclusion

SMRT-RenSeq is a powerful tool to identify NLR genes from wild species using the baits of the evolutionary related species with reference sequences. The availability of the NLRs from H. villosa provide a valuable library for R gene mining and transfer of disease resistance into wheat.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08334-w.

Expression pattern of CRYAB and CTGF genes in two pig breeds at different altitudes

ABSTRACT Tibetan pigs are characterized by significant phenotypic differences relative to lowland pigs. Our previous study demonstrated that the genes CRYAB and CTGF were differentially expressed in heart tissues between Tibetan (highland breed) and Yorkshire (lowland breed) pigs, indicating that they might participate in hypoxia adaptation. CRYAB (ɑB-crystallin) and CTGF (connective tissue growth factor) have also been reported to be associated with lung development. However, the expression patterns of CRYAB and CTGF in lung tissues at different altitudes and their genetic characterization are not well understood. In this study, qRT-PCR and western blot of lung tissue revealed higher CRYAB expression levels in highland and middle-highland Tibetan and Yorkshire pigs than in their lowland counterparts. With an increase in altitude, the expression level of CTGF increased in Tibetan pigs, whereas it decreased in Yorkshire pigs. Furthermore, two novel single-nucleotide polymorphism were identified in the 5′ flanking region of CRYAB (g.39644482C>T and g.39644132T>C) and CTGF (g.31671748A>G and g.31671773T>G). The polymorphism may partially contribute to the differences in expression levels between groups at the same altitude. These findings provide novel insights into the high-altitude hypoxia adaptations of Tibetan pigs.

[Clinicopathological characteristics of SMARCB1(INI1)-deficient sinonasal carcinoma]

Objective: To investigate the clinicopathological characteristics, diagnosis, differential diagnosis and prognostic factors of SMARCB1 (INI1)-deficient sinonasal carcinoma (SDSC). Methods: Sixteen cases of SDSC diagnosed in the Department of Pathology, Beijing Tongren Hospital from January 2016 to September 2020 were enrolled. Ninety-nine cases of small round cell malignant tumors of the head and neck were selected as the control, including poorly-differentiated squamous cell carcinoma (n=10), poorly-differentiated adenocarcinoma (n=5), undifferentiated carcinoma (SNUC, n=4), NUT carcinoma (n=5), neuroendocrine carcinoma (n=10), and other non-epithelial tumors [olfactory neuroblastoma (n=10), rhabdomyosarcoma (n=10), NK/T-cell lymphoma (n=10), malignant melanoma (n=10), Ewing's sarcoma/primitive neuroectodermal tumor (EWS/PNET, n=5)] and non-keratinizing undifferentiated nasopharyngeal carcinoma (n=20). The clinical and pathologic characteristics of SDSC, and immunohistochemical (IHC) expression of broad-spectrum CKpan, CK7, CK8/18, CK5/6, p63, p40, p16, INI1, NUT and neuroendocrine markers (Syn, CgA, CD56) were evaluated. In situ hybridization (ISH) was used to detect EBER and fluorescence in situ hybridization (FISH) to detect INI1 gene deletion. Results: The 16 cases of SDSC accounted for 1.3% (16/1 218) of all malignant sinonasal tumors in the author's unit during this time period, and 2.4% (16/657) of all malignant epithelial tumors. Microscopically, there was no clear squamous and adenomatous differentiation, but "rhabdoid-like" cells, are often seen. All SDSC cases were positive for CKpan and CK8/18, negative for INI1; Epstein-Barr virus was not detected by ISH; and INI1 gene deletion was observed in all 11 SDSC patients with FISH. Twelve cases were followed up for 3-47 months. One died of tumor-related diseases half a year after diagnosis, and the remaining patients were alive with tumor, the longest survival time was 47 months. Conclusion: SDSC should be differentiated from a variety of poorly-differentiated tumors in the sinonasal area. Histologically, SDSC has no clear differentiation, but the tumor cells are characteristically basal-like or rhabdoid-like, with non-specific vacuoles, translucent or vacuolar nuclei, prominent nucleoli and necrotic foci. They are negative for INI1 IHC staining, and FISH demonstrates INI1 gene deletion. The clinical prognosis is still unclear, further studies on its biologic behavior and treatment methods are warranted.