Live-attenuated virus vaccine defective in RNAi suppression induces rapid protection in neonatal and adult mice lacking mature B and T cells.

Global control of infectious diseases depends on the continuous development and deployment of diverse vaccination strategies. Currently available live-attenuated and killed virus vaccines typically take a week or longer to activate specific protection by the adaptive immunity. The mosquito-transmitted Nodamura virus (NoV) is attenuated in mice by mutations that prevent expression of the B2 viral suppressor of RNA interference (VSR) and consequently, drastically enhance in vivo production of the virus-targeting small-interfering RNAs. We reported recently that 2 d after immunization with live-attenuated VSR-disabled NoV (NoVΔB2), neonatal mice become fully protected against lethal NoV challenge and develop no detectable infection. Using Rag1-/- mice that produce no mature B and T lymphocytes as a model, here we examined the hypothesis that adaptive immunity is dispensable for the RNAi-based protective immunity activated by NoVΔB2 immunization. We show that immunization of both neonatal and adult Rag1-/- mice with live but not killed NoVΔB2 induces full protection against NoV challenge at 2 or 14 d postimmunization. Moreover, NoVΔB2-induced protective antiviral immunity is virus-specific and remains effective in adult Rag1-/- mice 42 and 90 d after a single-shot immunization. We conclude that immunization with the live-attenuated VSR-disabled RNA virus vaccine activates rapid and long-lasting protective immunity against lethal challenges by a distinct mechanism independent of the adaptive immunity mediated by B and T cells. Future studies are warranted to determine whether additional animal and human viruses attenuated by VSR inactivation induce similar protective immunity in healthy and adaptive immunity-compromised individuals.

Ubiquitin-induced RNF168 condensation promotes DNA double-strand break repair.

Rapid accumulation of repair factors at DNA double-strand breaks (DSBs) is essential for DSB repair. Several factors involved in DSB repair have been found undergoing liquid-liquid phase separation (LLPS) at DSB sites to facilitate DNA repair. RNF168, a RING-type E3 ubiquitin ligase, catalyzes H2A.X ubiquitination for recruiting DNA repair factors. Yet, whether RNF168 undergoes LLPS at DSB sites remains unclear. Here, we identified K63-linked polyubiquitin-triggered RNF168 condensation which further promoted RNF168-mediated DSB repair. RNF168 formed liquid-like condensates upon irradiation in the nucleus while purified RNF168 protein also condensed in vitro. An intrinsically disordered region containing amino acids 460-550 was identified as the essential domain for RNF168 condensation. Interestingly, LLPS of RNF168 was significantly enhanced by K63-linked polyubiquitin chains, and LLPS largely enhanced the RNF168-mediated H2A.X ubiquitination, suggesting a positive feedback loop to facilitate RNF168 rapid accumulation and its catalytic activity. Functionally, LLPS deficiency of RNF168 resulted in delayed recruitment of 53BP1 and BRCA1 and subsequent impairment in DSB repair. Taken together, our finding demonstrates the pivotal effect of LLPS in RNF168-mediated DSB repair.

Mammalian viral suppressors of RNA interference.

The antiviral defense directed by the RNAi pathway employs distinct specificity and effector mechanisms compared with other immune responses. The specificity of antiviral RNAi is programmed by siRNAs processed from virus-derived double-stranded RNA by Dicer endonuclease. Argonaute-containing RNA-induced silencing complex loaded with the viral siRNAs acts as the effector to mediate specific virus clearance by RNAi. Recent studies have provided evidence for the production and antiviral function of virus-derived siRNAs in both undifferentiated and differentiated mammalian cells infected with a range of RNA viruses when the cognate virus-encoded suppressor of RNAi (VSR) is rendered nonfunctional. In this review, we discuss the function, mechanism, and evolutionary origin of the validated mammalian VSRs and cell culture assays for their identification.

Leveraging the local genetic structure for trans-ancestry association mapping.

Over the past two decades, genome-wide association studies (GWASs) have successfully advanced our understanding of the genetic basis of complex traits. Despite the fruitful discovery of GWASs, most GWAS samples are collected from European populations, and these GWASs are often criticized for their lack of ancestry diversity. Trans-ancestry association mapping (TRAM) offers an exciting opportunity to fill the gap of disparities in genetic studies between non-Europeans and Europeans. Here, we propose a statistical method, LOG-TRAM, to leverage the local genetic architecture for TRAM. By using biobank-scale datasets, we showed that LOG-TRAM can greatly improve the statistical power of identifying risk variants in under-represented populations while producing well-calibrated p values. We applied LOG-TRAM to the GWAS summary statistics of various complex traits/diseases from BioBank Japan, UK Biobank, and African populations. We obtained substantial gains in power and achieved effective correction of confounding biases in TRAM. Finally, we showed that LOG-TRAM can be successfully applied to identify ancestry-specific loci and the LOG-TRAM output can be further used for construction of more accurate polygenic risk scores in under-represented populations.

RAP80 phase separation at DNA double-strand break promotes BRCA1 recruitment.

RAP80 has been characterized as a component of the BRCA1-A complex and is responsible for the recruitment of BRCA1 to DNA double-strand breaks (DSBs). However, we and others found that the recruitment of RAP80 and BRCA1 were not absolutely temporally synchronized, indicating that other mechanisms, apart from physical interaction, might be implicated. Recently, liquid-liquid phase separation (LLPS) has been characterized as a novel mechanism for the organization of key signaling molecules to drive their particular cellular functions. Here, we characterized that RAP80 LLPS at DSB was required for RAP80-mediated BRCA1 recruitment. Both cellular and in vitro experiments showed that RAP80 phase separated at DSB, which was ascribed to a highly disordered region (IDR) at its N-terminal. Meanwhile, the Lys63-linked poly-ubiquitin chains that quickly formed after DSBs occur, strongly enhanced RAP80 phase separation and were responsible for the induction of RAP80 condensation at the DSB site. Most importantly, abolishing the condensation of RAP80 significantly suppressed the formation of BRCA1 foci, encovering a pivotal role of RAP80 condensates in BRCA1 recruitment and radiosensitivity. Together, our study disclosed a new mechanism underlying RAP80-mediated BRCA1 recruitment, which provided new insight into the role of phase separation in DSB repair.

Mendelian randomization for causal inference accounting for pleiotropy and sample structure using genome-wide summary statistics.

Mendelian randomization (MR) is a valuable tool for inferring causal relationships among a wide range of traits using summary statistics from genome-wide association studies (GWASs). Existing summary-level MR methods often rely on strong assumptions, resulting in many false-positive findings. To relax MR assumptions, ongoing research has been primarily focused on accounting for confounding due to pleiotropy. Here, we show that sample structure is another major confounding factor, including population stratification, cryptic relatedness, and sample overlap. We propose a unified MR approach, MR-APSS, which 1) accounts for pleiotropy and sample structure simultaneously by leveraging genome-wide information; and 2) allows the inclusion of more genetic variants with moderate effects as instrument variables (IVs) to improve statistical power without inflating type I errors. We first evaluated MR-APSS using comprehensive simulations and negative controls and then applied MR-APSS to study the causal relationships among a collection of diverse complex traits. The results suggest that MR-APSS can better identify plausible causal relationships with high reliability. In particular, MR-APSS can perform well for highly polygenic traits, where the IV strengths tend to be relatively weak and existing summary-level MR methods for causal inference are vulnerable to confounding effects.

Ni/Photoredox-Catalyzed Enantioselective Acylation of α-Bromobenzoates with Aldehydes: A Formal Approach to Aldehyde-Aldehyde Cross-Coupling.

The catalytic cross-coupling of identical or similar functional groups is a cornerstone strategy for carbon-carbon bond formation, as exemplified by renowned methods, such as olefin cross-metathesis, Kolbe electrolysis, and various cross-electrophile couplings. However, similar methodologies for coupling aldehydes─fundamental building blocks in organic synthesis─remain underdeveloped. While the benzoin-type condensation, first reported in 1832, offers a reliable route for aldehyde dimerization, the chemo- and enantioselective cross-coupling of nonidentical yet similar aldehydes remains an unsolved challenge. Herein, we report a unified platform enabling highly chemo- and enantioselective cross-coupling of aldehydes. By leveraging nickel photoredox catalysis in tandem with discrete activation strategies for each aldehyde, this mechanistically distinct approach facilitates the enantioselective union of an aldehyde-derived α-oxy radical with an acyl radical, photocatalytically generated from a distinct aldehyde. This novel strategy enables modular access to enantioenriched α-oxygenated ketones with two minimally differentiated aliphatic substituents, a feat not achievable with existing chemocatalytic or biocatalytic techniques. The synthetic utility of this method is demonstrated by its application in the streamlined asymmetric synthesis of various medicinally relevant molecules. Additionally, mechanistic investigations rationalize the versatility of nickel photoredox catalysis to exploit new pathways for addressing long-standing synthetic challenges.

A unified framework for cross-population trait prediction by leveraging the genetic correlation of polygenic traits.

The development of polygenic risk scores (PRSs) has proved useful to stratify the general European population into different risk groups. However, PRSs are less accurate in non-European populations due to genetic differences across different populations. To improve the prediction accuracy in non-European populations, we propose a cross-population analysis framework for PRS construction with both individual-level (XPA) and summary-level (XPASS) GWAS data. By leveraging trans-ancestry genetic correlation, our methods can borrow information from the Biobank-scale European population data to improve risk prediction in the non-European populations. Our framework can also incorporate population-specific effects to further improve construction of PRS. With innovations in data structure and algorithm design, our methods provide a substantial saving in computational time and memory usage. Through comprehensive simulation studies, we show that our framework provides accurate, efficient, and robust PRS construction across a range of genetic architectures. In a Chinese cohort, our methods achieved 7.3%-198.0% accuracy gain for height and 19.5%-313.3% accuracy gain for body mass index (BMI) in terms of predictive R2 compared to existing PRS approaches. We also show that XPA and XPASS can achieve substantial improvement for construction of height PRSs in the African population, suggesting the generality of our framework across global populations.

FIRM: Flexible integration of single-cell RNA-sequencing data for large-scale multi-tissue cell atlas datasets.

Single-cell RNA-sequencing (scRNA-seq) is being used extensively to measure the mRNA expression of individual cells from deconstructed tissues, organs and even entire organisms to generate cell atlas references, leading to discoveries of novel cell types and deeper insight into biological trajectories. These massive datasets are usually collected from many samples using different scRNA-seq technology platforms, including the popular SMART-Seq2 (SS2) and 10X platforms. Inherent heterogeneities between platforms, tissues and other batch effects make scRNA-seq data difficult to compare and integrate, especially in large-scale cell atlas efforts; yet, accurate integration is essential for gaining deeper insights into cell biology. We present FIRM, a re-scaling algorithm which accounts for the effects of cell type compositions, and achieve accurate integration of scRNA-seq datasets across multiple tissue types, platforms and experimental batches. Compared with existing state-of-the-art integration methods, FIRM provides accurate mixing of shared cell type identities and superior preservation of original structure without overcorrection, generating robust integrated datasets for downstream exploration and analysis. FIRM is also a facile way to transfer cell type labels and annotations from one dataset to another, making it a reliable and versatile tool for scRNA-seq analysis, especially for cell atlas data integration.

PALM: a powerful and adaptive latent model for prioritizing risk variants with functional annotations.

MOTIVATION: The findings from genome-wide association studies (GWASs) have greatly helped us to understand the genetic basis of human complex traits and diseases. Despite the tremendous progress, much effects are still needed to address several major challenges arising in GWAS. First, most GWAS hits are located in the non-coding region of human genome, and thus their biological functions largely remain unknown. Second, due to the polygenicity of human complex traits and diseases, many genetic risk variants with weak or moderate effects have not been identified yet. RESULTS: To address the above challenges, we propose a powerful and adaptive latent model (PALM) to integrate cell-type/tissue-specific functional annotations with GWAS summary statistics. Unlike existing methods, which are mainly based on linear models, PALM leverages a tree ensemble to adaptively characterize non-linear relationship between functional annotations and the association status of genetic variants. To make PALM scalable to millions of variants and hundreds of functional annotations, we develop a functional gradient-based expectation-maximization algorithm, to fit the tree-based non-linear model in a stable manner. Through comprehensive simulation studies, we show that PALM not only controls false discovery rate well, but also improves statistical power of identifying risk variants. We also apply PALM to integrate summary statistics of 30 GWASs with 127 cell type/tissue-specific functional annotations. The results indicate that PALM can identify more risk variants as well as rank the importance of functional annotations, yielding better interpretation of GWAS results. AVAILABILITY AND IMPLEMENTATION: The source code is available at https://github.com/YangLabHKUST/PALM. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Epitope delimitation: A new method for defining epitopes of human IgG-reactive antigenic peptides based on rabbit-recognized epitope motifs.

The use of precise epitope peptides as antigens is essential for accurate serological diagnosis of viral-infected individuals, but now it remains an unsolvable problem for mapping precise B cell epitopes (BCEs) recognized by human serum. To address this challenge, we propose a novel epitope delimitation (ED) method to uncover BCEs in the delineated human IgG-reactive (HR) antigenic peptides (APs). Specifically, the method based on the rationale of similarities in humoral immune responses between mammalian species consists of a pair of elements: experimentally delineated HR-AP and rabbit-recognized (RR) BCE motif and corresponding pair of sequence alignment analysis. As a result of using the ED approach, after decoding four RR-epitomes of human papillomavirus types 16/18-E6 and E7 proteins utilizing rabbit serum against each recombinant protein and sequence alignment analysis of HR-APs and RR-BCEs, 19 fine BCEs in 17 of 22 known HR-APs were defined based on each corresponding RR-BCE motifs, including the type-specificity of each delimited BCE in homologous proteins. The test with 22 known 16/20mer HR-APs demonstrated that the ED method is effective and efficient, indicating that it can be used as an alternative method to the conventional identification of fine BCEs using overlapping 8mer peptides.

Development of a universal method for vibrational analysis of the terminal alkyne C≡C stretch.

The terminal alkyne C≡C stretch has a large Raman scattering cross section in the "silent" region for biomolecules. This has led to many Raman tag and probe studies using this moiety to study biomolecular systems. A computational investigation of these systems is vital to aid in the interpretation of these results. In this work, we develop a method for computing terminal alkyne vibrational frequencies and isotropic transition polarizabilities that can easily and accurately be applied to any terminal alkyne molecule. We apply the discrete variable representation method to a localized version of the C≡C stretch normal mode. The errors of (1) vibrational localization to the terminal alkyne moiety, (2) anharmonic normal mode isolation, and (3) discretization of the Born-Oppenheimer potential energy surface are quantified and found to be generally small and cancel each other. This results in a method with low error compared to other anharmonic vibrational methods like second-order vibrational perturbation theory and to experiments. Several density functionals are tested using the method, and TPSS-D3, an inexpensive nonempirical density functional with dispersion corrections, is found to perform surprisingly well. Diffuse basis functions are found to be important for the accuracy of computed frequencies. Finally, the computation of vibrational properties like isotropic transition polarizabilities and the universality of the localized normal mode for terminal alkynes are demonstrated.

Associations between non-insulin-based insulin resistance indices and heart failure prevalence in overweight/obesity adults without diabetes mellitus: evidence from the NHANES 2001-2018.

BACKGROUND: The triglyceride glucose (TyG) index and triglyceride-to-high-density lipoprotein cholesterol (TG/HDL-C) ratio are recognized as simple non-insulin-based insulin resistance indices. Our study aimed to explore the relationship between these two indicators and heart failure (HF) in overweight or obesity individuals without diabetes. METHODS: This cross-sectional study selected 13,473 participants from the National Health and Nutrition Examination Survey (NHANES) 2001-2018 dataset. Weighted multivariable logistic regression and subgroup analysis were employed to evaluate the relationships between TyG index, TG/HDL-C ratio, and HF prevalence, respectively. Additionally, smooth curve fitting was utilized to analyze the dose-response relationships. RESULTS: A total of 13,473 obesity or overweight people without diabetes were included in this study through screening, among whom 291 (2.16%) had comorbid HF. The results of multivariable logistic regression suggested that the highest TyG index (OR = 2.4, 95% CI = 1.4-4.2, p = 0.002) and the highest TG/HDL-C ratio (OR = 1.2, 95% CI = 1.1-1.3, p < 0.001) both increased the prevalence of HF, especially in the non-Hispanic population. Dose-response relationships suggested nonlinear relationships between these two indicators and HF. CONCLUSION: Our study demonstrated that elevated TyG index and TG/HDL-C ratio were closely associated with the prevalence of HF, and both exhibited nonlinear relationships with HF prevalence in overweight/obesity adults without diabetes. Based on these findings, additional prospective studies are needed for further validation.

AggMapNet: enhanced and explainable low-sample omics deep learning with feature-aggregated multi-channel networks.

Omics-based biomedical learning frequently relies on data of high-dimensions (up to thousands) and low-sample sizes (dozens to hundreds), which challenges efficient deep learning (DL) algorithms, particularly for low-sample omics investigations. Here, an unsupervised novel feature aggregation tool AggMap was developed to Aggregate and Map omics features into multi-channel 2D spatial-correlated image-like feature maps (Fmaps) based on their intrinsic correlations. AggMap exhibits strong feature reconstruction capabilities on a randomized benchmark dataset, outperforming existing methods. With AggMap multi-channel Fmaps as inputs, newly-developed multi-channel DL AggMapNet models outperformed the state-of-the-art machine learning models on 18 low-sample omics benchmark tasks. AggMapNet exhibited better robustness in learning noisy data and disease classification. The AggMapNet explainable module Simply-explainer identified key metabolites and proteins for COVID-19 detections and severity predictions. The unsupervised AggMap algorithm of good feature restructuring abilities combined with supervised explainable AggMapNet architecture establish a pipeline for enhanced learning and interpretability of low-sample omics data.

Antibiotic-free production of sucrose isomerase in Bacillus subtilis by genome integration.

Sucrose isomerase (SIase) catalyzes the hydrolysis and isomerization of sucrose to form isomaltulose, a valuable functional sugar widely used in the food industry. However, the lack of safe and efficient heterologous expression systems hinders SIase production and application. In this study, we achieved antibiotic-free SIase expression in Bacillus subtilis through genome integration. Using CRISPR/Cas9 system, SIase expression cassettes were integrated into various genomic loci, including amyE and ctc, both individually and in combination, resulting in single-copy and muti-copy integration strains. Engineered strains with a maltose-inducible promoter effectively expressed and secreted SIase. Notably, multi-copy strain exhibited enhanced SIase production, achieving 4.4 U/mL extracellular activity in shake flask cultivations. Furthermore, crude enzyme solution from engineered strain transformed high concentrations sucrose into high yields of isomaltulose, reaching a maximum yield of 94.6%. These findings demonstrate antibiotic-free SIase production in B. subtilis via genome integration, laying the foundation for its industrial production and application.

Neoadjuvant Chemotherapy With CAPOX Versus Chemoradiation for Locally Advanced Rectal Cancer With Uninvolved Mesorectal Fascia (CONVERT): Initial Results of a Phase III Trial.

OBJECTIVE: To compare neoadjuvant chemotherapy (nCT) with CAPOX alone versus neoadjuvant chemoradiotherapy (nCRT) with capecitabine in locally advanced rectal cancer (LARC) with uninvolved mesorectal fascia (MRF). BACKGROUND DATA: nCRT is associated with higher surgical complications, worse long-term functional outcomes, and questionable survival benefits. Comparatively, nCT alone seems a promising alternative treatment in lower-risk LARC patients with uninvolved MRF. METHODS: Patients between June 2014 and October 2020 with LARC within 12 cm from the anal verge and uninvolved MRF were randomly assigned to nCT group with 4 cycles of CAPOX (Oxaliplatin 130 mg/m2 IV day 1 and Capecitabine 1000 mg/m2 twice daily for 14 d. Repeat every 3 wk) or nCRT group with Capecitabine 825 mg/m² twice daily administered orally and concurrently with radiation therapy (50 Gy/25 fractions) for 5 days per week. The primary end point is local-regional recurrence-free survival. Here we reported the results of secondary end points: histopathologic response, surgical events, and toxicity. RESULTS: Of the 663 initially enrolled patients, 589 received the allocated treatment (nCT, n=300; nCRT, n=289). Pathologic complete response rate was 11.0% (95% CI, 7.8-15.3%) in the nCT arm and 13.8% (95% CI, 10.1-18.5%) in the nCRT arm ( P =0.33). The downstaging (ypStage 0 to 1) rate was 40.8% (95% CI, 35.1-46.7%) in the nCT arm and 45.6% (95% CI, 39.7-51.7%) in the nCRT arm ( P =0.27). nCT was associated with lower perioperative distant metastases rate (0.7% vs. 3.1%, P =0.03) and preventive ileostomy rate (52.2% vs. 63.6%, P =0.008) compared with nCRT. Four patients in the nCT arm received salvage nCRT because of local disease progression after nCT. Two patients in the nCT arm and 5 in the nCRT arm achieved complete clinical response and were treated with a nonsurgical approach. Similar results were observed in subgroup analysis. CONCLUSIONS: nCT achieved similar pCR and downstaging rates with lower incidence of perioperative distant metastasis and preventive ileostomy compared with nCRT. CAPOX could be an effective alternative to neoadjuvant therapy in LARC with uninvolved MRF. Long-term follow-up is needed to confirm these results.

XPXP: improving polygenic prediction by cross-population and cross-phenotype analysis.

MOTIVATION: As increasing sample sizes from genome-wide association studies (GWASs), polygenic risk scores (PRSs) have shown great potential in personalized medicine with disease risk prediction, prevention and treatment. However, the PRS constructed using European samples becomes less accurate when it is applied to individuals from non-European populations. It is an urgent task to improve the accuracy of PRSs in under-represented populations, such as African populations and East Asian populations. RESULTS: In this article, we propose a cross-population and cross-phenotype (XPXP) method for construction of PRSs in under-represented populations. XPXP can construct accurate PRSs by leveraging biobank-scale datasets in European populations and multiple GWASs of genetically correlated phenotypes. XPXP also allows to incorporate population-specific and phenotype-specific effects, and thus further improves the accuracy of PRS. Through comprehensive simulation studies and real data analysis, we demonstrated that our XPXP outperformed existing PRS approaches. We showed that the height PRSs constructed by XPXP achieved 9% and 18% improvement over the runner-up method in terms of predicted R2 in East Asian and African populations, respectively. We also showed that XPXP substantially improved the stratification ability in identifying individuals at high genetic risk of type 2 diabetes. AVAILABILITY AND IMPLEMENTATION: The XPXP software and all analysis code are available at github.com/YangLabHKUST/XPXP. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Efficient Dicer processing of virus-derived double-stranded RNAs and its modulation by RIG-I-like receptor LGP2.

The interferon-regulated antiviral responses are essential for the induction of both innate and adaptive immunity in mammals. Production of virus-derived small-interfering RNAs (vsiRNAs) to restrict virus infection by RNA interference (RNAi) is a recently identified mammalian immune response to several RNA viruses, which cause important human diseases such as influenza and Zika virus. However, little is known about Dicer processing of viral double-stranded RNA replicative intermediates (dsRNA-vRIs) in mammalian somatic cells. Here we show that infected somatic cells produced more influenza vsiRNAs than cellular microRNAs when both were produced by human Dicer expressed de novo, indicating that dsRNA-vRIs are not poor Dicer substrates as previously proposed according to in vitro Dicer processing of synthetic long dsRNA. We report the first evidence both for canonical vsiRNA production during wild-type Nodamura virus infection and direct vsiRNA sequestration by its RNAi suppressor protein B2 in two strains of suckling mice. Moreover, Sindbis virus (SINV) accumulation in vivo was decreased by prior production of SINV-targeting vsiRNAs triggered by infection and increased by heterologous expression of B2 in cis from SINV genome, indicating an antiviral function for the induced RNAi response. These findings reveal that unlike artificial long dsRNA, dsRNA-vRIs made during authentic infection of mature somatic cells are efficiently processed by Dicer into vsiRNAs to direct antiviral RNAi. Interestingly, Dicer processing of dsRNA-vRIs into vsiRNAs was inhibited by LGP2 (laboratory of genetics and physiology 2), which was encoded by an interferon-stimulated gene (ISG) shown recently to inhibit Dicer processing of artificial long dsRNA in cell culture. Our work thus further suggests negative modulation of antiviral RNAi by a known ISG from the interferon response.

Effects of neoadjuvant radiochemotherapy for anorectal function in locally advanced rectal cancer patients: a study protocol for a prospective, observational, controlled, multicentre study.

BACKGROUND: Neoadjuvant chemoradiotherapy (NCRT) and total mesorectal excision are standard treatment regimen for patients with locally advanced rectal cancer (LARC). This sphincter-saving treatment strategy may be accompanied by a series of anorectal functional disorders. Yet, prospective studies that dynamically evaluating the respective roles of radiotherapy, chemotherapy and surgery on anorectal function are lacking. PATIENTS/DESIGN: The study is a prospective, observational, controlled, multicentre study. After screening for eligibility and obtaining informed consent, a total of 402 LARC patients undergoing NCRT followed by surgery, or neoadjuvant chemotherapy followed by surgery, or surgery only would be included in the trial. The primary outcome measure is the average resting pressure of anal sphincter. The secondary outcome measures are maximum anal sphincter contraction pressure, Wexner continence score and low anterior resection syndrome (LARS) score. Evaluations will be carried out at the following stages: baseline (T1), after radiotherapy or chemotherapy (before surgery, T2), after surgery (before closing the temporary stoma, T3), and at follow-up visits (every 3 to 6 months, T4, T5……). Follow-up for each patient will be at least 2 years. DISCUSSION: We expect the program to provide more information of neoadjuvant radiotherapy and/or chemotherapy on anorectal function, and to optimize the treatment strategy to reduce anorectal dysfunction for LARC patients. TRIAL REGISTRATION: ClinicalTrials.gov (NCT05671809). Registered on 26 December 2022.

Characterization of the chloroplast genome of medicinal herb Polygonatum cyrtonema and identification of molecular markers by comparative analysis.

Polygonatum cyrtonema Hua is a traditional Chinese herb medicine, and it is widely distributed in China. The intrageneric taxonomy and phylogenetic relationships within Polygonatum have long been controversial due to their morphological similarity and lacking special DNA barcodes. In this paper, the complete chloroplast genome is a relatively conserved quadripartite structure including a large single copy region of 84 711 bp, a small single copy region of 18 210 bp, and a pair of inverted repeats region of 26 142 bp. A total of 342 simple sequence repeats were identified, and most of them were found to be composed of A/T, including 126 mono-nucleotides and 179 di-nucleotides. Nucleotide diversity was analyzed and eight highly variable regions (psbl∼trnT-CGU, atpF∼atpH, trnT-GGU∼psbD, psaJ∼rps20, trnL-UAG∼ndhD, ndhG∼ndhl, ndhA, and rpl32∼ccsA) were identified as potential molecular markers. Phylogenetic analysis based on the whole chloroplast genome showed that P. cyrtonema, within the family Asparagaceae, is closely related to Polygonatum sibiricum and Polygonatum kingianum. The sequence matK, trnT-GGU & ccsA, and ndhG∼ndhA were identified as three DNA barcodes. The assembly and comparative analysis of P. cyrtonema complete chloroplast genome will provide essential molecular information about the evolution and molecular biology for further study.