Targeting Mitochondria-Located circRNA SCAR Alleviates NASH via Reducing mROS Output.

Mitochondria, which play central roles in immunometabolic diseases, have their own genome. However, the functions of mitochondria-located noncoding RNAs are largely unknown due to the absence of a specific delivery system. By circular RNA (circRNA) expression profile analysis of liver fibroblasts from patients with nonalcoholic steatohepatitis (NASH), we observe that mitochondrial circRNAs account for a considerable fraction of downregulated circRNAs in NASH fibroblasts. By constructing mitochondria-targeting nanoparticles, we observe that Steatohepatitis-associated circRNA ATP5B Regulator (SCAR), which is located in mitochondria, inhibits mitochondrial ROS (mROS) output and fibroblast activation. circRNA SCAR, mediated by PGC-1α, binds to ATP5B and shuts down mPTP by blocking CypD-mPTP interaction. Lipid overload inhibits PGC-1α by endoplasmic reticulum (ER) stress-induced CHOP. In vivo, targeting circRNA SCAR alleviates high fat diet-induced cirrhosis and insulin resistance. Clinically, circRNA SCAR is associated with steatosis-to-NASH progression. Collectively, we identify a mitochondrial circRNA that drives metaflammation and serves as a therapeutic target for NASH.

Complement Signals Determine Opposite Effects of B Cells in Chemotherapy-Induced Immunity.

Understanding molecular mechanisms that dictate B cell diversity is important for targeting B cells as anti-cancer treatment. Through the single-cell dissection of B cell heterogeneity in longitudinal samples of patients with breast cancer before and after neoadjuvant chemotherapy, we revealed that an ICOSL+ B cell subset emerges after chemotherapy. Using three immunocompetent mouse models, we recapitulated the subset switch of human tumor-infiltrating B cells during chemotherapy. By employing B-cell-specific deletion mice, we showed that ICOSL in B cells boosts anti-tumor immunity by enhancing the effector to regulatory T cell ratio. The signature of ICOSL+ B cells is imprinted by complement-CR2 signaling, which is triggered by immunogenic cell death. Moreover, we identified that CD55, a complement inhibitory protein, determines the opposite roles of B cells in chemotherapy. Collectively, we demonstrated a critical role of the B cell subset switch in chemotherapy response, which has implications in designing novel anti-cancer therapies. VIDEO ABSTRACT.

Noncoding RNAs in skeletal development and disorders.

Protein-encoding genes only constitute less than 2% of total human genomic sequences, and 98% of genetic information was previously referred to as "junk DNA". Meanwhile, non-coding RNAs (ncRNAs) consist of approximately 60% of the transcriptional output of human cells. Thousands of ncRNAs have been identified in recent decades, and their essential roles in the regulation of gene expression in diverse cellular pathways associated with fundamental cell processes, including proliferation, differentiation, apoptosis, and metabolism, have been extensively investigated. Furthermore, the gene regulation networks they form modulate gene expression in normal development and under pathological conditions. In this review, we integrate current information about the classification, biogenesis, and function of ncRNAs and how these ncRNAs support skeletal development through their regulation of critical genes and signaling pathways in vivo. We also summarize the updated knowledge of ncRNAs involved in common skeletal diseases and disorders, including but not limited to osteoporosis, osteoarthritis, rheumatoid arthritis, scoliosis, and intervertebral disc degeneration, by highlighting their roles established from in vivo, in vitro, and ex vivo studies.

EuRBPDB: a comprehensive resource for annotation, functional and oncological investigation of eukaryotic RNA binding proteins (RBPs).

RNA binding proteins (RBPs) are a large protein family that plays important roles at almost all levels of gene regulation through interacting with RNAs, and contributes to numerous biological processes. However, the complete list of eukaryotic RBPs including human is still unavailable. Here, we systematically identified RBPs in 162 eukaryotic species based on both computational analysis of RNA binding domains (RBDs) and large-scale RNA binding proteomic data, and established a comprehensive eukaryotic RBP database, EuRBPDB (http://EuRBPDB.syshospital.org). We identified a total of 311 571 RBPs with RBDs (corresponding to 6368 ortholog groups) and 3,651 non-canonical RBPs without known RBDs. EuRBPDB provides detailed annotations for each RBP, including basic information and functional annotation. Moreover, we systematically investigated RBPs in the context of cancer biology based on published literatures, PPI-network and large-scale omics data. To facilitate the exploration of the clinical relevance of RBPs, we additionally designed a cancer web interface to systematically and interactively display the biological features of RBPs in various types of cancers. EuRBPDB has a user-friendly web interface with browse and search functions, as well as data downloading function. We expect that EuRBPDB will be a widely-used resource and platform for both the communities of RNA biology and cancer biology.

The subunit of RNA N6-methyladenosine methyltransferase OsFIP regulates early degeneration of microspores in rice.

N6-Methyladenosine (m6A) RNA methylation plays important roles during development in different species. However, knowledge of m6A RNA methylation in monocots remains limited. In this study, we reported that OsFIP and OsMTA2 are the components of m6A RNA methyltransferase complex in rice and uncovered a previously unknown function of m6A RNA methylation in regulation of plant sporogenesis. Importantly, OsFIP is essential for rice male gametogenesis. Knocking out of OsFIP results in early degeneration of microspores at the vacuolated pollen stage and simultaneously causes abnormal meiosis in prophase I. We further analyzed the profile of rice m6A modification during sporogenesis in both WT and OsFIP loss-of-function plants, and identified a rice panicle specific m6A modification motif "UGWAMH". Interestingly, we found that OsFIP directly mediates the m6A methylation of a set of threonine protease and NTPase mRNAs and is essential for their expression and/or splicing, which in turn regulates the progress of sporogenesis. Our findings revealed for the first time that OsFIP plays an indispensable role in plant early sporogenesis. This study also provides evidence for the different functions of the m6A RNA methyltransferase complex between rice and Arabidopsis.

Reprogramming of m6A epitranscriptome is crucial for shaping of transcriptome and proteome in response to hypoxia.

Hypoxia causes a series of responses supporting cells to survive in harsh environments. Substantial post-transcriptional and translational regulation during hypoxia has been observed. However, detailed regulatory mechanism in response to hypoxia is still far from complete. RNA m6A modification has been proven to govern the life cycle of RNAs. Here, we reported that total m6A level of mRNAs was decreased during hypoxia, which might be mediated by the induction of m6A eraser, ALKBH5. Meanwhile, expression levels of most YTH family members of m6A readers were systematically down-regulated. Transcriptome-wide analysis of m6A revealed a drastic reprogramming of m6A epitranscriptome during cellular hypoxia. Integration of m6A epitranscriptome with either RNA-seq based transcriptome analysis or mass spectrometry (LC-MS/MS) based proteome analysis of cells upon hypoxic stress revealed that reprogramming of m6A epitranscriptome reshaped the transcriptome and proteome, thereby supporting efficient generation of energy for adaption to hypoxia. Moreover, ATP production was blocked when silencing an m6A eraser, ALKBH5, under hypoxic condition, demonstrating that m6A pathway is an important regulator during hypoxic response. Collectively, our studies indicate that crosstalk between m6A and HIF1 pathway is essential for cellular response to hypoxia, providing insights into the underlying molecular mechanisms during hypoxia.

Overexpression of lncRNAs with endogenous lengths and functions using a lncRNA delivery system based on transposon.

BACKGROUND: Long noncoding RNAs (lncRNAs) play important roles in many physiological and pathological processes, this indicates that lncRNAs can serve as potential targets for gene therapy. Stable expression is a fundamental technology in the study of lncRNAs. The lentivirus is one of the most widely used delivery systems for stable expression. However, it was initially designed for mRNAs, and the applicability of lentiviral vectors for lncRNAs is largely unknown. RESULTS: We found that the lentiviral vector produces lncRNAs with improper termination, appending an extra fragment of ~ 2 kb to the 3'-end. Consequently, the secondary structures were changed, the RNA-protein interactions were blocked, and the functions were impaired in certain lncRNAs, which indicated that lentiviral vectors are not ideal delivery systems of lncRNAs. Here, we developed a novel lncRNA delivery method called the Expression of LncRNAs with Endogenous Characteristics using the Transposon System (ELECTS). By inserting a termination signal after the lncRNA sequence, ELECTS produces transcripts without 3'-flanking sequences and retains the native features and function of lncRNAs, which cannot be achieved by lentiviral vectors. Moreover, ELECTS presents no potential risk of infection for the operators and it takes much less time. ELECTS provides a reliable, convenient, safe, and efficient delivery method for stable expression of lncRNAs. CONCLUSIONS: Our study demonstrated that improper transcriptional termination from lentiviral vectors have fundamental effects on molecular action and cellular function of lncRNAs. The ELECTS system developed in this study will provide a convenient and reliable method for the lncRNA study.

CIRDES: an efficient genome-wide method for in vivo RNA-RNA interactome analysis.

Complex RNA-RNA interactions underlie fundamental biological processes. However, a large number of RNA-RNA interactions remain unknown. Most existing methods used to map RNA-RNA interactions are based on proximity ligation, but these strategies also capture a huge amount of intramolecular RNA secondary structures, making it almost impossible to detect most RNA-RNA interactions. To overcome this limitation, we developed an efficient, genome-wide method, Capture Interacting RNA and Deep Sequencing (CIRDES) for in vivo capturing of the RNA interactome. We designed multiple 20-nt CIRDES probes tiling the whole RNA sequence of interest. This strategy obtained high selectivity and low background noise proved by qRT-PCR data. CIRDES enriched target RNA and its interacting RNAs from cells crosslinked by formaldehyde in high efficiency. After hybridization and purification, the captured RNAs were converted to the cDNA library after a highly efficient ligation to a 3' end infrared-dye-conjugated RNA adapter based on adapter ligation library construction. Using CIRDES, we detected highly abundant known interacting RNA, as well as a large number of novel targets of U6 snRNA. The enrichment of U4 snRNA, which interacts with U6, confirmed the robustness of the identification of the RNA-RNA interaction by CIRDES. These results suggest that the CIRDES is an efficient strategy for genome-wide RNA-RNA interactome analysis.

Deep sequencing reveals a global reprogramming of lncRNA transcriptome during EMT.

Several studies have shown that long non-coding RNAs (lncRNAs) may play an essential role in Epithelial-Mesenchymal Transition (EMT), which is an important step in tumor metastasis; however, little is known about the global change of lncRNA transcriptome during EMT. To investigate how lncRNA transcriptome alterations contribute to EMT progression regulation, we deep-sequenced the whole-transcriptome of MCF10A as the cells underwent TGF-ß-induced EMT. RESULTS: Deep-sequencing results showed that the long RNA transcriptome of MCF10A had undergone global changes as early as 8h after treatment with TGF-ß. The expression of 3403 known and novel lncRNAs, and 570 known and novel circRNAs were altered during EMT. To identify the key lncRNA-regulator, we constructed the co-expression network and found all junction nodes in the network are lncRNAs. One junction node, RP6-65G23.5, was further verified as a key regulator of EMT. Intriguingly, we identified 216 clusters containing lncRNAs which were located in "gene desert" regions. The expressions of all lncRNAs in these clusters changed concurrently during EMT, strongly suggesting that these clusters might play important roles in EMT. Our study reveals a global reprogramming of lncRNAs transcriptome during EMT and provides clues for the future study of the molecular mechanism of EMT.

Both endo-siRNAs and tRNA-derived small RNAs are involved in the differentiation of primitive eukaryote Giardia lamblia.

Small RNAs (sRNAs), including microRNAs and endogenous siRNAs (endo-siRNAs), regulate most important biologic processes in eukaryotes, such as cell division and differentiation. Although sRNAs have been extensively studied in various eukaryotes, the role of sRNAs in the early emergence of eukaryotes is unclear. To address these questions, we deep sequenced the sRNA transcriptome of four different stages in the differentiation of Giardia lamblia, one of the most primitive eukaryotes. We identified a large number of endo-siRNAs in this fascinating parasitic protozoan and found that they were produced from live telomeric retrotransposons and three genomic regions (i.e., endo-siRNA generating regions [eSGRs]). eSGR-derived endo-siRNAs were proven to target mRNAs in trans. Gradual up-regulation of endo-siRNAs in the differentiation of Giardia suggested that they might be involved in the regulation of this process. This hypothesis was supported by the impairment of the differentiation ability of Giardia when GLDICER, essential for the biogenesis of endo-siRNAs, was knocked down. Endo-siRNAs are not the only sRNA regulators in Giardia differentiation, because a great number of tRNAs-derived sRNAs showed more dramatic expression changes than endo-siRNAs in this process. We totally identified five novel kinds of tRNAs-derived sRNAs and found that the biogenesis in four of them might be correlated with that of stress-induced tRNA-derived RNA (sitRNA), which was discovered in our previous studies. Our studies reveal an unexpected complex panorama of sRNA in G. lamblia and shed light on the origin and functional evolution of eukaryotic sRNAs.

Comparative transcriptome analysis of small noncoding RNAs in different stages of Trypanosoma brucei.

Trypanosoma brucei, a pathogen of human and domestic animals, is an early evolved parasitic protozoan with a complex life cycle. Most genes of this parasite are post-transcriptionally regulated. However, the mechanisms and the molecules involved remain largely unknown. We have deep-sequenced the small RNAs of two life stages of this parasite--the bloodstream form and the procyclic form. Our results show that the small RNAs of T. brucei could derive from multiple sources, including NATs (natural antisense transcripts), tRNAs, and rRNAs. Most of these small RNAs in the two stages were found to share uniform characteristics. However, our results demonstrate that their variety and expression show significant differences between different stages, indicating possible functional differentiation. Dicer-knockdown evidence further proved that some of the small interfering RNAs (siRNAs) could regulate the expression of genes. Based on the genome-wide analysis of the small RNAs in the two stages of T. brucei, our results not only provide evidence to study their differentiation but also shed light on questions regarding the origins and evolution of small RNA-based mechanisms in early eukaryotes.

Pseudogene-derived small interference RNAs regulate gene expression in African Trypanosoma brucei.

Pseudogenes have been shown to acquire unique regulatory roles from more and more organisms. We report the observation of a cluster of siRNAs derived from pseudogenes of African Trypanosoma brucei using high through-put analysis. We show that these pseudogene-derived siRNAs suppress gene expression through RNA interference. The discovery that siRNAs may originate from pseudogenes and regulate gene expression in a unicellular eukaryote provides insights into the functional roles of pseudogenes and into the origin of noncoding small RNAs.

Glycolysis and tumor progression promoted by the m6A writer VIRMA via m6A-dependent upregulation of STRA6 in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is one of the most aggressive and lethal malignancies, highlighting the urgent need to elucidate the underlying oncogenic mechanisms. VIRMA is a classic isoform of methyltransferases that participates in epigenetic transcriptomic modification in eukaryotic mRNAs. However, the exact roles of VIRMA in PDAC remain unclear. Here, we identified that VIRMA is highly expressed in PDAC, and histone modifications of the promoter may partly account for this dysregulation. Moreover, VIRMA is closely related to glycolysis and poor prognosis in PDAC. We further determined that STRA6 is a direct downstream target of VIRMA in PDAC by RNA sequencing (RNA-seq) and m6A sequencing (m6A-seq). VIRMA is involved in gene expression regulation via 3' UTR targeting of STRA6 mRNA. Furthermore, the m6A reader IGF2BP2 was shown to critically contribute to the stability of STRA6 mRNA. We describe the role of VIRMA in promoting signaling via the STRA6/STAT3 axis, which results in increased levels of HIF-1α, a key activator of glycolysis. In vivo and in vitro experiments reveal that the VIRMA-STRA6-STAT3-HIF-1α axis plays an instrumental role in glycolysis and tumor progression in PDAC. In conclusion, we demonstrate that VIRMA can increase glycolysis in PDAC by upregulating STRA6, a cell surface membrane protein that stimulates the STAT3 pathway, thereby activating HIF-1α and leading to pancreatic cancer malignancy. Overall, our data strongly suggest that the VIRMA-STRA6-STAT3-HIF-1α axis is a viable therapeutic target in PDAC.

Niche stiffness sustains cancer stemness via TAZ and NANOG phase separation.

Emerging evidence shows that the biomechanical environment is required to support cancer stem cells (CSCs), which play a crucial role in drug resistance. However, how mechanotransduction signals regulate CSCs and its clinical significance has remained unclear. Using clinical-practice ultrasound elastography for patients' lesions and atomic force microscopy for surgical samples, we reveal that increased matrix stiffness is associated with poor responses to neoadjuvant chemotherapy, worse prognosis, and CSC enrichment in patients with breast cancer. Mechanically, TAZ activated by biomechanics enhances CSC properties via phase separation with NANOG. TAZ-NANOG phase separation, which is dependent on acidic residues in the N-terminal activation domain of NANOG, promotes the transcription of SOX2 and OCT4. Therapeutically, targeting NANOG or TAZ reduces CSCs and enhances the chemosensitivity in vivo. Collectively, this study demonstrated that the phase separation of a pluripotency transcription factor links mechanical cues in the niche to the fate of CSCs.

Drastic expression change of transposon-derived piRNA-like RNAs and microRNAs in early stages of chicken embryos implies a role in gastrulation.

Recent studies have shown that endogenous small RNAs regulate a variety of biological processes during vertebrate development; however, little is known about the role of small RNAs in regulating developmental signaling pathways during early embryogenesis. In this study, we applied Illumina sequencing to characterize an unexpected endogenous small RNA catalog and demonstrated a dramatic transition from transposon-derived piRNA-like small RNAs (pilRNAs) to microRNAs (miRNAs) in pre- and post-gastrula chicken embryos. The comprehensive expression profile of chicken miRNAs at the pre- and post-gastrula stages revealed that most known and new miRNAs were dynamically regulated during development. In addition to embryonic stem cell-related miRNAs, Gene Ontology (GO) analysis showed that miRNAs enriched in early stage chicken embryos targeted multiple signal transduction pathways associated with the reproductive process and embryogenesis, including Wnt and TGF-ß, which specifies the neural fate of blastodermal cells. Intriguingly, a large cohort of pilRNAs primarily derived from the active and most abundant transposable elements (TEs) were enriched in chicken stage X blastoderms. Within stage X blastoderms, pilRNAs were specifically localized to the primordial germ cells (PGCs), indicating their post-zygotic origin. Together, these findings imply a role for small RNAs in gastrulation in early stage chicken embryos.

Genome-wide evolutionary analysis of the noncoding RNA genes and noncoding DNA of Paramecium tetraurelia.

The compact genome of the unicellular eukaryote Paramecium tetraurelia contains noncoding DNA (ncDNA) distributed into >39,000 intergenic sequences and >90,000 introns of 390 base pairs (bp) and 25 bp on average, respectively. Here we analyzed the molecular features of the ncRNA genes, introns, and intergenic sequences of this genome. We mainly used computational programs and comparative genomics possible because the P. tetraurelia genome had formed throughout whole-genome duplications (WGDs). We characterized 417 5S rRNA, snRNA, snoRNA, SRP RNA, and tRNA putative genes, 415 of which map within intergenic sequences, and two, within introns. The evolution of these ncRNA genes appears to have mainly involved purifying selection and gene deletion. We then compared the introns that interrupt the protein-coding gene duplicates arisen from the recent WGD and identified a population of a few thousands of introns having evolved under most stringent constraints (>95% of identity). We also showed that low nucleotide substitution levels characterize the 50 and 80-115 base pairs flanking, respectively, the stop and start codons of the protein-coding genes. Lower substitution levels mark the base pairs flanking the highly transcribed genes, or the start codons of the genes of the sets with a high number of WGD-related sequences. Finally, adjacent to protein-coding genes, we characterized 32 DNA motifs able to encode stable and evolutionary conserved RNA secondary structures and defining putative expression controlling elements. Fourteen DNA motifs with similar properties map distant from protein-coding genes and may encode regulatory ncRNAs.

mirExplorer: detecting microRNAs from genome and next generation sequencing data using the AdaBoost method with transition probability matrix and combined features.

microRNAs (miRNAs) represent an abundant group of small regulatory non-coding RNAs in eukaryotes. The emergence of Next-generation sequencing (NGS) technologies has allowed the systematic detection of small RNAs (sRNAs) and de novo sequencing of genomes quickly and with low cost. As a result, there is an increased need to develop fast miRNA prediction tools to annotate miRNAs from various organisms with a high level of accuracy, using the genome sequence or the NGS data. Several miRNA predictors have been proposed to achieve this purpose. However, the accuracy and fitness for multiple species of existing predictors needed to be improved. Here, we present a novel prediction tool called mirExplorer, which is based on an integrated adaptive boosting method and contains two modules. The first module named mirExplorer-genome was designed to de novo predict pre-miRNAs from genome, and the second module named mirExplorer-NGS was used to discover miRNAs from NGS data. A set of novel features of pre-miRNA secondary structure and miRNA biogenesis has been extracted to distinguish real pre-miRNAs from pseudo ones. We used outer-ten-fold cross-validation to verify the mirExplorer-genome computation, which obtained a specificity of 95.03% and a sensitivity of 93.71% on human data. This computation was made on test data from 16 species, and it achieved an overall accuracy of 95.53%. Systematic outer-ten-fold cross-validation of the mirExplorer-NGS model achieved a specificity of 98.3% and a sensitivity of 97.72%. We found that the good performance of the mirExplorer-NGS model was upheld across species from vertebrates to plants in test datasets. The mirExplorer is available as both web server and software package at http://biocenter.sysu.edu.cn/mir/.

Dynamic m6A mRNA Methylation Reveals the Role of METTL3/14-m6A-MNK2-ERK Signaling Axis in Skeletal Muscle Differentiation and Regeneration.

N6-methyladenosine (m6A) RNA methylation has emerged as an important factor in various biological processes by regulating gene expression. However, the dynamic profile, function and underlying molecular mechanism of m6A modification during skeletal myogenesis remain elusive. Here, we report that members of the m6A core methyltransferase complex, METTL3 and METTL14, are downregulated during skeletal muscle development. Overexpression of either METTL3 or METTL14 dramatically blocks myotubes formation. Correspondingly, knockdown of METTL3 or METTL14 accelerates the differentiation of skeletal muscle cells. Genome-wide transcriptome analysis suggests ERK/MAPK is the downstream signaling pathway that is regulated to the greatest extent by METTL3/METTL14. Indeed, METTL3/METTL14 expression facilitates ERK/MAPK signaling. Via MeRIP-seq, we found that MNK2, a critical regulator of ERK/MAPK signaling, is m6A modified and is a direct target of METTL3/METTL14. We further revealed that YTHDF1 is a potential reader of m6A on MNK2, regulating MNK2 protein levels without affecting mRNA levels. Furthermore, we discovered that METTL3/14-MNK2 axis was up-regulated notably after acute skeletal muscle injury. Collectively, our studies revealed that the m6A writers METTL3/METTL14 and the m6A reader YTHDF1 orchestrate MNK2 expression posttranscriptionally and thus control ERK signaling, which is required for the maintenance of muscle myogenesis and may contribute to regeneration.

Stress-induced tRNA-derived RNAs: a novel class of small RNAs in the primitive eukaryote Giardia lamblia.

Giardia lamblia is an early diverging and evolutionarily successful protozoan as it can enter into a dormant cyst stage from a vegetative trophozoite. During dormant stage, its metabolic rate decreases dramatically. However, to date, the regulatory molecules participating in the initiation and maintenance of this process have not been fully investigated. In this study, we have identified a class of abundant small RNAs named sitRNAs, which are approximately 46 nucleotides in length and accumulate in G. lamblia encysting cultures. Remarkably, they are derived from the 3' portion of fully matured tRNAs by cleavage of the anticodon left arm, with the 3' terminal CCA triplex still connected. During differentiation, only a limited portion of mature tRNAs is cleaved, but this cleavage occurs almost in the entire tRNA family. sitRNAs begin to accumulate as early as 3 h after initiation of encystation and are maintained at a relatively stable level during the whole process, exhibiting an expression peak at around 24 hr. Our studies further show that sitRNAs can be induced by several other stress factors, and in the case of serum deprivation, both tRNAs and sitRNAs degrade rapidly, with the accumulation of tRNA being halved. Our results may provide new insight into a novel mechanism for stressed G. lamblia to regulate gene expression globally.

Protocol for Single-Cell Analysis of Tumor-Infiltrating B Cells Isolated from Human Breast Cancer Tissue Before and After Neo-adjuvant Chemotherapy.

Single-cell analysis of tumor-infiltrating lymphocytes obtained before and after preoperative therapy reflects the dynamic interplay of the tumor and immune system during treatment. Here, we present a protocol to implement single-cell analysis of tumor-infiltrating B cells, which were isolated from paired human breast cancers before and after neo-adjuvant chemotherapy. This protocol also facilitates isolation and single-cell analysis of other tumor-infiltrating lymphocytes. For complete information on the generation and use of this protocol, please refer to Lu et al. (2020).