Agave REVEILLE1 regulates the onset and release of seasonal dormancy in Populus

Deciduous woody plants like poplar (Populus spp.) have seasonal bud dormancy. It has been challenging to simultaneously delay the onset of bud dormancy in the fall and advance bud break in the spring, as bud dormancy, and bud break were thought to be controlled by different genetic factors. Here, we demonstrate that heterologous expression of the REVEILLE1 gene (named AaRVE1) from Agave (Agave americana) not only delays the onset of bud dormancy but also accelerates bud break in poplar in field trials. AaRVE1 heterologous expression increases poplar biomass yield by 166% in the greenhouse. Furthermore, we reveal that heterologous expression of AaRVE1 increases cytokinin contents, represses multiple dormancy-related genes, and up-regulates bud break-related genes, and that AaRVE1 functions as a transcriptional repressor and regulates the activity of the DORMANCY-ASSOCIATED PROTEIN 1 (DRM1) promoter. Our findings demonstrate that AaRVE1 appears to function as a regulator of bud dormancy and bud break, which has important implications for extending the growing season of deciduous trees in frost-free temperate and subtropical regions to increase crop yield.

AlphaFold2 models indicate that protein sequence determines both structure and dynamics

AlphaFold 2 (AF2) has placed Molecular Biology in a new era where we can visualize, analyze and interpret the structures and functions of all proteins solely from their primary sequences. We performed AF2 structure predictions for various protein systems, including globular proteins, a multi-domain protein, an intrinsically disordered protein (IDP), a randomized protein, two larger proteins (> 1000 AA), a heterodimer and a homodimer protein complex. Our results show that along with the three dimensional (3D) structures, AF2 also decodes protein sequences into residue flexibilities via both the predicted local distance difference test (pLDDT) scores of the models, and the predicted aligned error (PAE) maps. We show that PAE maps from AF2 are correlated with the distance variation (DV) matrices from molecular dynamics (MD) simulations, which reveals that the PAE maps can predict the dynamical nature of protein residues. Here, we introduce the AF2-scores, which are simply derived from pLDDT scores and are in the range of [0, 1]. We found that for most protein models, including large proteins and protein complexes, the AF2-scores are highly correlated with the root mean square fluctuations (RMSF) calculated from MD simulations. However, for an IDP and a randomized protein, the AF2-scores do not correlate with the RMSF from MD, especially for the IDP. Our results indicate that the protein structures predicted by AF2 also convey information of the residue flexibility, i.e., protein dynamics.

Exosome-mediated transfer of MIF confers temozolomide resistance by regulating TIMP3/PI3K/AKT axis in gliomas

Temozolomide (TMZ) resistance is an important cause of clinical treatment failure and poor prognosis in gliomas. Increasing evidence indicates that cancer-derived exosomes contribute to chemoresistance; however, the specific contribution of glioma-derived exosomes remains unclear. The aim of this study was to explore the role and underlying mechanisms of exosomal macrophage migration inhibitory factor (MIF) on TMZ resistance in gliomas. We first demonstrated that MIF was upregulated in the exosomes of TMZ-resistant cells, engendering the transfer of TMZ resistance to sensitive cells. Our results indicated that exosomal MIF conferred TMZ resistance to sensitive cells through the enhancement of cell proliferation and the repression of cell apoptosis upon TMZ exposure. MIF knockdown enhanced TMZ sensitivity in resistant glioma cells by upregulating Metalloproteinase Inhibitor 3 (TIMP3) and subsequently suppressing the PI3K/AKT signaling pathway. Additionally, exosomal MIF promoted tumor growth and TMZ resistance of glioma cells in vivo, while IOS-1 (MIF inhibitor) promotes glioma TMZ sensitive in vivo. Taken together, our study demonstrated that exosome-mediated transfer of MIF enhanced TMZ resistance in glioma through downregulating TIMP3 and further activating the PI3K/AKT signaling pathway, highlighting a prognostic biomarker and promising therapeutic target for TMZ treatment in gliomas.

A quantitative yeast aging proteomics analysis reveals novel aging regulators

Calorie restriction (CR) is the most robust longevity intervention, extending lifespan from yeast to mammals. Numerous conserved pathways regulating aging and mediating CR have been identified; however, the overall proteomic changes during these conditions remain largely unexplored. We compared proteomes between young and replicatively aged yeast cells under normal and CR conditions using the Stable-Isotope Labeling by Amino acids in Cell culture (SILAC) quantitative proteomics and discovered distinct signatures in the aging proteome. We found remarkable proteomic similarities between aged and CR cells, including induction of stress response pathways, providing evidence that CR pathways are engaged in aged cells. These observations also uncovered aberrant changes in mitochondria membrane proteins as well as a proteolytic cellular state in old cells. These proteomics analyses help identify potential genes and pathways that have causal effects on longevity.

Protein interaction potential landscapes for yeast replicative aging

We proposed a novel interaction potential landscape approach to map the systems-level profile changes of gene networks during replicative aging in Saccharomyces cerevisiae. This approach enabled us to apply quasi-potentials, the negative logarithm of the probabilities, to calibrate the elevation of the interaction landscapes with young cells as a reference state. Our approach detected opposite landscape changes based on protein abundances from transcript levels, especially for intra-essential gene interactions. We showed that essential proteins play different roles from hub proteins on the age-dependent interaction potential landscapes. We verified that hub proteins tend to avoid other hub proteins, but essential proteins prefer to interact with other essential proteins. Overall, we showed that the interaction potential landscape is promising for inferring network profile change during aging and that the essential hub proteins may play an important role in the uncoupling between protein and transcript levels during replicative aging.

Complementary performances of convolutional and capsule neural networks on classifying microfluidic images of dividing yeast cells

Microfluidic-based assays have become effective high-throughput approaches to examining replicative aging of budding yeast cells. Deep learning may offer an efficient way to analyze a large number of images collected from microfluidic experiments. Here, we compare three deep learning architectures to classify microfluidic time-lapse images of dividing yeast cells into categories that represent different stages in the yeast replicative aging process. We found that convolutional neural networks outperformed capsule networks in terms of accuracy, precision, and recall. The capsule networks had the most robust performance in detecting one specific category of cell images. An ensemble of three best-fitted single-architecture models achieves the highest overall accuracy, precision, and recall due to complementary performances. In addition, extending classification classes and data augmentation of the training dataset can improve the predictions of the biological categories in our study. This work lays a useful framework for sophisticated deep-learning processing of microfluidic-based assays of yeast replicative aging.

Overexpression of an Agave Phosphoenolpyruvate Carboxylase Improves Plant Growth and Stress Tolerance

It has been challenging to simultaneously improve photosynthesis and stress tolerance in plants. Crassulacean acid metabolism (CAM) is a CO2-concentrating mechanism that facilitates plant adaptation to water-limited environments. We hypothesized that the ectopic expression of a CAM-specific phosphoenolpyruvate carboxylase (PEPC), an enzyme that catalyzes primary CO2 fixation in CAM plants, would enhance both photosynthesis and abiotic stress tolerance. To test this hypothesis, we engineered a CAM-specific PEPC gene (named AaPEPC1) from Agave americana into tobacco. In comparison with wild-type and empty vector controls, transgenic tobacco plants constitutively expressing AaPEPC1 showed a higher photosynthetic rate and biomass production under normal conditions, along with significant carbon metabolism changes in malate accumulation, the carbon isotope ratio δ13C, and the expression of multiple orthologs of CAM-related genes. Furthermore, AaPEPC1 overexpression enhanced proline biosynthesis, and improved salt and drought tolerance in the transgenic plants. Under salt and drought stress conditions, the dry weight of transgenic tobacco plants overexpressing AaPEPC1 was increased by up to 81.8% and 37.2%, respectively, in comparison with wild-type plants. Our findings open a new door to the simultaneous improvement of photosynthesis and stress tolerance in plants.

Atomistic probing of aptameric binding of CD19 outer membrane domain reveals an "aptamer walking" mechanism

We propose an integrated structural approach to search potential aptamer molecules for targeting cancer receptor proteins. We used the outer cellular domain of the B-lymphocyte antigen, CD19, as the target for this study. First, using available protein-aptamer structures deposited in the protein data bank as resources, structural annotation was performed to seek the most probable binding aptamer and its potential initial configuration to the CD19 structure. Using this initial structure, molecular dynamics (MD) simulations were performed for adjustment of the aptamer-binding. During this process, we observed an "aptamer walking" mechanism of the binding of the single-stranded RNA-aptamer to CD19: the aptamer molecule gradually adjusts its configurations and shifts toward favorable binding positions. However, the target molecule CD19 maintained a relatively stable conformation during this process. The interface area between the RNA-aptamer and CD19 increased from less than 8 nm² to over 12 nm² during a 2-μs MD simulation. Using a stable binding pose as the starting structure, we manually mutated the RNA-aptamer to a DNA-aptamer and found that the interface area was further increased to over 16 nm² , indicating a stronger affinity compared to the RNA-aptamer. The RNA- and DNA-aptamers and their stable binding-poses to the CD19 molecule may be used as templates in designing potential aptamer molecules that target the B-cell marker molecule CD19 with enhanced specificity and stability.

[Numerical simulation study of type B aortic dissection using patient-specific reverse engineering and fluid-structure interaction]

Objective: To construct computational fluid model of type B aortic dissection using patient-specific reverse engineering and fluid-structure interaction, and evaluate the application of computational fluid model on aortic remodeling of type B aortic dissection. Methods: Consecutive computed tomographic angiograph data was acquired from a patient with type B aortic dissection at initial diagnosis, 1 week and 6 years after endovascular repair of primary tear entry and 3 months after endovascular repair of distal tear erosion. Three-dimensional model of aortic dissection was reversely reconstructed by Mimics, and then the model was smoothened by Geomagic. Computational fluid dynamic numerical simulation was performed in ANSYS by the means of two-way fluid-structure interaction, and the relation between blood dynamic characteristic and thrombosed remodeling of type B aortic dissection was evaluated. Results: The computational fluid model of type B aortic dissection using patient-specific reverse engineering and fluid-structure interaction method was successfully constructed. Local peak of blood pressure on the convex surface of junction at aortic arch and descending aorta was found. The wall stress was much higher at the false lumen than that at the true lumen, and the peak of wall stress converged on the edge and tear entry of false lumen. After the exclusion of proximal tear entry, the blood streamline was decreased significantly and flowed reversely. Blood flow in the remaining false lumen was retrograded from the entry at left iliac artery and formed turbulence at the top of false lumen, which was benefit for dissection thrombus remodeling. The higher pressure at the false lumen was associated with the new formation of aortic aneurysm at the distal tear. Conclusion: The computational fluid model of aortic dissection based on patient-specific reverse engineering and fluid-structure interaction method can successfully reveal the relatively truly blood dynamic and wall pressure characteristic of type B aortic dissection.

A New Calmodulin-Binding Protein Expresses in the Context of Secondary Cell Wall Biosynthesis and Impacts Biomass Properties in Populus

A greater understanding of biosynthesis, signaling and regulatory pathways involved in determining stem growth and secondary cell wall chemistry is important for enabling pathway engineering and genetic optimization of biomass properties. The present study describes a new functional role of PdIQD10, a Populus gene belonging to the IQ67-Domain1 family of IQD genes, in impacting biomass formation and chemistry. Expression studies showed that PdIQD10 has enhanced expression in developing xylem and tension-stressed tissues in Populus deltoides. Molecular dynamics simulation and yeast two-hybrid interaction experiments suggest interactions with two calmodulin proteins, CaM247 and CaM014, supporting the sequence-predicted functional role of the PdIQD10 as a calmodulin-binding protein. PdIQD10 was found to interact with specific Populus isoforms of the Kinesin Light Chain protein family, shown previously to function as microtubule-guided, cargo binding and delivery proteins in Arabidopsis. Subcellular localization studies showed that PdIQD10 localizes in the nucleus and plasma membrane regions. Promoter-binding assays suggest that a known master transcriptional regulator of secondary cell wall biosynthesis (PdWND1B) may be upstream of an HD-ZIP III gene that is in turn upstream of PdIQD10 gene in the transcriptional network. RNAi-mediated downregulation of PdIQD10 expression resulted in plants with altered biomass properties including higher cellulose, wall glucose content and greater biomass quantity. These results present evidence in support of a new functional role for an IQD gene family member, PdIQD10, in secondary cell wall biosynthesis and biomass formation in Populus.

Diel rewiring and positive selection of ancient plant proteins enabled evolution of CAM photosynthesis in Agave

Background

Crassulacean acid metabolism (CAM) enhances plant water-use efficiency through an inverse day/night pattern of stomatal closure/opening that facilitates nocturnal CO₂ uptake. CAM has evolved independently in over 35 plant lineages, accounting for ~ 6% of all higher plants. Agave species are highly heat- and drought-tolerant, and have been domesticated as model CAM crops for beverage, fiber, and biofuel production in semi-arid and arid regions. However, the genomic basis of evolutionary innovation of CAM in genus Agave is largely unknown.

Results

Using an approach that integrated genomics, gene co-expression networks, comparative genomics and protein structure analyses, we investigated the molecular evolution of CAM as exemplified in Agave. Comparative genomics analyses among C₃, C₄ and CAM species revealed that core metabolic components required for CAM have ancient genomic origins traceable to non-vascular plants while regulatory proteins required for diel re-programming of metabolism have a more recent origin shared among C₃, C₄ and CAM species. We showed that accelerated evolution of key functional domains in proteins responsible for primary metabolism and signaling, together with a diel re-programming of the transcription of genes involved in carbon fixation, carbohydrate processing, redox homeostasis, and circadian control is required for the evolution of CAM in Agave. Furthermore, we highlighted the potential candidates contributing to the adaptation of CAM functional modules.

Conclusions

This work provides evidence of adaptive evolution of CAM related pathways. We showed that the core metabolic components required for CAM are shared by non-vascular plants, but regulatory proteins involved in re-reprogramming of carbon fixation and metabolite transportation appeared more recently. We propose that the accelerated evolution of key proteins together with a diel re-programming of gene expression were required for CAM evolution from C₃ ancestors in Agave.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4964-7) contains supplementary material, which is available to authorized users.

A 5-Enolpyruvylshikimate 3-Phosphate Synthase Functions as a Transcriptional Repressor in Populus

Long-lived perennial plants, with distinctive habits of inter-annual growth, defense, and physiology, are of great economic and ecological importance. However, some biological mechanisms resulting from genome duplication and functional divergence of genes in these systems remain poorly studied. Here, we discovered an association between a poplar (Populus trichocarpa) 5-enolpyruvylshikimate 3-phosphate synthase gene (PtrEPSP) and lignin biosynthesis. Functional characterization of PtrEPSP revealed that this isoform possesses a helix-turn-helix motif in the N terminus and can function as a transcriptional repressor that regulates expression of genes in the phenylpropanoid pathway in addition to performing its canonical biosynthesis function in the shikimate pathway. We demonstrated that this isoform can localize in the nucleus and specifically binds to the promoter and represses the expression of a SLEEPER-like transcriptional regulator, which itself specifically binds to the promoter and represses the expression of PtrMYB021 (known as MYB46 in Arabidopsis thaliana), a master regulator of the phenylpropanoid pathway and lignin biosynthesis. Analyses of overexpression and RNAi lines targeting PtrEPSP confirmed the predicted changes in PtrMYB021 expression patterns. These results demonstrate that PtrEPSP in its regulatory form and PtrhAT form a transcriptional hierarchy regulating phenylpropanoid pathway and lignin biosynthesis in Populus.

The Kalanchoë genome provides insights into convergent evolution and building blocks of crassulacean acid metabolism

Crassulacean acid metabolism (CAM) is a water-use efficient adaptation of photosynthesis that has evolved independently many times in diverse lineages of flowering plants. We hypothesize that convergent evolution of protein sequence and temporal gene expression underpins the independent emergences of CAM from C3 photosynthesis. To test this hypothesis, we generate a de novo genome assembly and genome-wide transcript expression data for Kalanchoë fedtschenkoi, an obligate CAM species within the core eudicots with a relatively small genome (~260 Mb). Our comparative analyses identify signatures of convergence in protein sequence and re-scheduling of diel transcript expression of genes involved in nocturnal CO2 fixation, stomatal movement, heat tolerance, circadian clock, and carbohydrate metabolism in K. fedtschenkoi and other CAM species in comparison with non-CAM species. These findings provide new insights into molecular convergence and building blocks of CAM and will facilitate CAM-into-C3 photosynthesis engineering to enhance water-use efficiency in crops.

Populus trichocarpa encodes small, effector-like secreted proteins that are highly induced during mutualistic symbiosis

During symbiosis, organisms use a range of metabolic and protein-based signals to communicate. Of these protein signals, one class is defined as ‘effectors’, i.e., small secreted proteins (SSPs) that cause phenotypical and physiological changes in another organism. To date, protein-based effectors have been described in aphids, nematodes, fungi and bacteria. Using RNA sequencing of Populus trichocarpa roots in mutualistic symbiosis with the ectomycorrhizal fungus Laccaria bicolor, we sought to determine if host plants also contain genes encoding effector-like proteins. We identified 417 plant-encoded putative SSPs that were significantly regulated during this interaction, including 161 SSPs specific to P. trichocarpa and 15 SSPs exhibiting expansion in Populus and closely related lineages. We demonstrate that a subset of these SSPs can enter L. bicolor hyphae, localize to the nucleus and affect hyphal growth and morphology. We conclude that plants encode proteins that appear to function as effector proteins that may regulate symbiotic associations.

Overexpression of a Domain of Unknown Function 266-containing protein results in high cellulose content, reduced recalcitrance, and enhanced plant growth in the bioenergy crop Populus

BACKGROUND

Domain of Unknown Function 266 (DUF266) is a plant-specific domain. DUF266-containing proteins (DUF266 proteins) have been categorized as 'not classified glycosyltransferases (GTnc)' due to amino acid similarity with GTs. However, little is known about the function of DUF266 proteins.

RESULTS

Phylogenetic analysis revealed that DUF266 proteins are only present in the land plants including moss and lycophyte. We report the functional characterization of one member of DUF266 proteins in Populus, PdDUF266A. PdDUF266A was ubiquitously expressed with high abundance in the xylem. In Populus transgenic plants overexpressing PdDUF266A (OXPdDUF266A), the glucose and cellulose contents were significantly higher, while the lignin content was lower than that in the wild type. Degree of polymerization of cellulose in OXPdDUF266A transgenic plants was also higher, whereas cellulose crystallinity index remained unchanged. Gene expression analysis indicated that cellulose biosynthesis-related genes such as CESA and SUSY were upregulated in mature leaf and xylem of OXPdDUF266A transgenic plants. Moreover, PdDUF266A overexpression resulted in an increase of biomass production. Their glucose contents and biomass phenotypes were further validated via heterologous expression of PdDUF266A in Arabidopsis. Results from saccharification treatment demonstrated that the rate of sugar release was increased by approximately 38% in the OXPdDUF266A transgenic plants.

CONCLUSIONS

These results suggest that the overexpression of PdDUF266A can increase cellulose content, reduce recalcitrance, and enhance biomass production, and that PdDUF266A is a promising target for genetic manipulation for biofuel production.

Pinning effect of reactive elements on adhesion energy and adhesive strength of incoherent Al2O3/NiAl interface

The profound effects of reactive elements (REs) on the adhesion energy and adhesive strength of the α-Al2O3/β-NiAl interface in thermal barrier coating (TBC) systems have attracted increasing attention because RE-doping has played a significant role in improving the thermal cycling lifetime of TBCs. However, the fundamental mechanism is, so far, not well understood due to the experimental difficulty and theoretical complexity in interface modelling. For this purpose, in the present study we have performed comprehensive density functional theory calculations and information targeted experiments to underline the origin of the surprising enhancement of interface adhesion, stability and mechanical strength of the α-Al2O3/β-NiAl interface by different RE doping levels. Our results suggest that the interface failure firstly appears within the NiAl layer adjacent to the Al-terminated oxide under mechanical loading, while the formation of O-RE-Ni bond pairs at the interface can effectively hinder the interface de-cohesion, providing a higher mechanical strength. By comparing several typical REs, it is observed that Hf can emerge not only with the highest interface adhesion energy, but also the highest mechanical strength; in agreement with our experimental results. By continuously increasing the dopant concentration, the strengthening effect may increase correspondingly, but is limited by the solute solubility. These results shed light into the effect of REs on the stability and strength of the α-Al2O3/β-NiAl interface, providing theoretical guidance for interface design via a combinational analysis of bond topology and electronic structure.

Orchestration of carbohydrate processing for crassulacean acid metabolism

The production of phosphoenolpyruvate as a substrate for nocturnal CO2 uptake represents a significant sink for carbohydrate in CAM plants which has to be balanced with the provisioning of carbohydrate for growth and maintenance. In starch-storing CAM species, diversification in chloroplast metabolite transporters, and the deployment of both phosphorolytic and hydrolytic routes of starch degradation accommodate a division of labour in directing C-skeletons towards nocturnal carboxylation or production of sucrose for growth. In soluble-sugar storing CAM plants, the vacuole plays a central role in managing carbon homeostasis. The molecular identities of various types of vacuolar sugar transporters have only been identified for C3 species within the last 10 years. The recent availability of CAM genomes enables the identification of putative orthologues of vacuolar sugar transporters which represent strategic targets for orchestrating the diel provisioning of substrate for nocturnal carboxylation and growth.

The pineapple genome and the evolution of CAM photosynthesis

Pineapple (Ananas comosus (L.) Merr.) is the most economically valuable crop possessing crassulacean acid metabolism (CAM), a photosynthetic carbon assimilation pathway with high water-use efficiency, and the second most important tropical fruit. We sequenced the genomes of pineapple varieties F153 and MD2 and a wild pineapple relative, Ananas bracteatus accession CB5. The pineapple genome has one fewer ancient whole-genome duplication event than sequenced grass genomes and a conserved karyotype with seven chromosomes from before the ρ duplication event. The pineapple lineage has transitioned from C3 photosynthesis to CAM, with CAM-related genes exhibiting a diel expression pattern in photosynthetic tissues. CAM pathway genes were enriched with cis-regulatory elements associated with the regulation of circadian clock genes, providing the first cis-regulatory link between CAM and circadian clock regulation. Pineapple CAM photosynthesis evolved by the reconfiguration of pathways in C3 plants, through the regulatory neofunctionalization of preexisting genes and not through the acquisition of neofunctionalized genes via whole-genome or tandem gene duplication.

A roadmap for research on crassulacean acid metabolism (CAM) to enhance sustainable food and bioenergy production in a hotter, drier world

Crassulacean acid metabolism (CAM) is a specialized mode of photosynthesis that features nocturnal CO2 uptake, facilitates increased water-use efficiency (WUE), and enables CAM plants to inhabit water-limited environments such as semi-arid deserts or seasonally dry forests. Human population growth and global climate change now present challenges for agricultural production systems to increase food, feed, forage, fiber, and fuel production. One approach to meet these challenges is to increase reliance on CAM crops, such as Agave and Opuntia, for biomass production on semi-arid, abandoned, marginal, or degraded agricultural lands. Major research efforts are now underway to assess the productivity of CAM crop species and to harness the WUE of CAM by engineering this pathway into existing food, feed, and bioenergy crops. An improved understanding of CAM has potential for high returns on research investment. To exploit the potential of CAM crops and CAM bioengineering, it will be necessary to elucidate the evolution, genomic features, and regulatory mechanisms of CAM. Field trials and predictive models will be required to assess the productivity of CAM crops, while new synthetic biology approaches need to be developed for CAM engineering. Infrastructure will be needed for CAM model systems, field trials, mutant collections, and data management.

X-ray structure of a Hg2+ complex of mercuric reductase (MerA) and quantum mechanical/molecular mechanical study of Hg2+ transfer between the C-terminal and buried catalytic site cysteine pairs

Mercuric reductase, MerA, is a key enzyme in bacterial mercury resistance. This homodimeric enzyme captures and reduces toxic Hg²⁺ to Hg⁰, which is relatively unreactive and can exit the cell passively. Prior to reduction, the Hg²⁺ is transferred from a pair of cysteines (C558′ and C559′ using Tn501 numbering) at the C-terminus of one monomer to another pair of cysteines (C136 and C141) in the catalytic site of the other monomer. Here, we present the X-ray structure of the C-terminal Hg²⁺ complex of the C136A/C141A double mutant of the Tn501 MerA catalytic core and explore the molecular mechanism of this Hg transfer with quantum mechanical/molecular mechanical (QM/MM) calculations. The transfer is found to be nearly thermoneutral and to pass through a stable tricoordinated intermediate that is marginally less stable than the two end states. For the overall process, Hg²⁺ is always paired with at least two thiolates and thus is present at both the C-terminal and catalytic binding sites as a neutral complex. Prior to Hg²⁺ transfer, C141 is negatively charged. As Hg²⁺ is transferred into the catalytic site, a proton is transferred from C136 to C559′ while C558′ becomes negatively charged, resulting in the net transfer of a negative charge over a distance of ∼7.5 Å. Thus, the transport of this soft divalent cation is made energetically feasible by pairing a competition between multiple Cys thiols and/or thiolates for Hg²⁺ with a competition between the Hg²⁺ and protons for the thiolates.

Efficient isolation of high-quality RNA from lotus Nelumbo nucifera ssp nucifera tissues

Nelumbo nucifera is widely used as food, as an ornamental, in medicine, and as packing material; it is also reported to have anti-HIV effects and antioxidant capacity. We sought an improved method for extracting high-quality total RNA from different tissues of N. nucifera. Four methods for RNA extraction were assessed for their ability to recover high-quality RNA applicable for evaluation of polyphenol oxidase (PPO) gene expression profiles. The recovery and quality of the RNA obtained from five different tissues by the best CTAB-LiCl method were evaluated through UV light absorbance. Both A(260)/A(280) and A(260)/A(230) absorbance ratios were more than 2.0; the yield ranged from 59.87 to 163.75 μg/g fresh weight. The brightness of the 28S band was approximately twice that of 18S; the latter was also considered as high-quality RNA. The PPO gene fragment (606 bp) was successfully amplified by RT-PCR, demonstrating the integrity of the isolated RNA. The relative expression levels of the PPO gene based on RT-PCR in five tissues of lotus were: rhizome buds (2.66), young leaves (2.42), fresh cut rhizome (2.02), petals (1.80), and petiole (1.65), using housekeeping gene β-actin as an internal control. We concluded that the total RNA isolated by this protocol is of sufficient quality for molecular applications.

Time-dependent density functional theory assessment of UV absorption of benzoic acid derivatives

Benzoic acid (BA) derivatives of environmental relevance exhibit various photophysical and photochemical characteristics. Here, time-dependent density functional theory (TDDFT) is used to calculate photoexcitations of eight selected BAs and the results are compared with UV spectra determined experimentally. High-level gas-phase EOM-CCSD calculations and experimental aqueous-phase spectra were used as the references for the gas-phase and aqueous-phase TDDFT results, respectively. A cluster-continuum model was used in the aqueous-phase calculations. Among the 15 exchange-correlation (XC) functionals assessed, five functionals, including the meta-GGA hybrid M06-2X, double hybrid B2PLYPD, and range-separated functionals CAM-B3LYP, ωB97XD, and LC-ωPBE, were found to be in excellent agreement with the EOM-CCSD gas-phase calculations. These functionals furnished excitation energies consistent with the pH dependence of the experimental spectra with a standard deviation (STDEV) of ∼0.20 eV. A molecular orbital analysis revealed a πσ* feature of the low-lying transitions of the BAs. The CAM-B3LYP functional showed the best overall performance and therefore shows promise for TDDFT calculations of processes involving photoexcitations of benzoic acid derivatives.

Structure and conformational dynamics of the metalloregulator MerR upon binding of Hg(II)

The bacterial metalloregulator MerR is the index case of an eponymous family of regulatory proteins, which controls the transcription of a set of genes (the mer operon) conferring mercury resistance in many bacteria. Homodimeric MerR represses transcription in the absence of mercury and activates transcription upon Hg(II) binding. Here, the average structures of the apo and Hg(II)-bound forms of MerR in aqueous solution are examined using small-angle X-ray scattering, indicating an extended conformation of the metal-bound protein and revealing the existence of a novel compact conformation in the absence of Hg(II). Molecular dynamics (MD) simulations are performed to characterize the conformational dynamics of the Hg(II)-bound form. In both small-angle X-ray scattering and MD, the average torsional angle between DNA-binding domains is approximately 65 degrees. Furthermore, in MD, interdomain motions on a timescale of approximately 10 ns involving large-amplitude (approximately 20 A) domain opening-and-closing, coupled to approximately 40 degrees variations of interdomain torsional angle, are revealed. This correlated domain motion may propagate allosteric changes from the metal-binding site to the DNA-binding site while maintaining DNA contacts required to initiate DNA underwinding.

Mechanism of histone methylation catalyzed by protein lysine methyltransferase SET7/9 and origin of product specificity

Methylation of certain lysine residues in the N-terminal tails of core histone proteins in nucleosome is of fundamental importance in the regulation of chromatin structure and gene expression. Such histone modification is catalyzed by protein lysine methyltransferases (PKMTs). PKMTs contain a conserved SET domain in almost all of the cases and may transfer one to three methyl groups from S-adenosyl-L-methionine (AdoMet) to the epsilon-amino group of the target lysine residue. Here, quantum mechanical/molecular mechanical molecular dynamics and free-energy simulations are performed on human PKMT SET7/9 and its mutants to understand two outstanding questions for the reaction catalyzed by PKMTs: the mechanism for deprotonation of positively charged methyl lysine (lysine) and origin of product specificity. The results of the simulations suggest that Tyr-335 (an absolute conserved residue in PKMTs) may play the role as the general base for the deprotonation after dissociation of AdoHcy (S-adenosyl-L-homocysteine) and before binding of AdoMet. It is shown that conformational changes could bring Y335 to the target methyl lysine (lysine) for proton abstraction. This mechanism provides an explanation why methyl transfers could be catalyzed by PKMTs processively. The free-energy profiles for methyl transfers are reported and analyzed for wild type and certain mutants (Y305F and Y335F) and the active-site interactions that are of importance for the enzyme's function are discussed. The results of the simulations provide important insights into the catalytic process and lead to a better understanding of experimental observations concerning the origin of product specificity for PKMTs.

The QM/MM molecular dynamics and free energy simulations of the acylation reaction catalyzed by the serine-carboxyl peptidase kumamolisin-As

Quantum mechanical/molecular mechanical molecular dynamics and free energy simulations are performed to study the acylation reaction catalyzed by kumamolisin-As, a serine-carboxyl peptidase, and to elucidate the catalytic mechanism and the origin of substrate specificity. It is demonstrated that the nucleophilic attack by the serine residue on the substrate may not be the rate-limiting step for the acylation of the GPH*FF substrate. The present study also confirms the earlier suggestions that Asp164 acts as a general acid during the catalysis and that the electrostatic oxyanion hole interactions may not be sufficient to lead a stable tetrahedral intermediate along the reaction pathway. Moreover, Asp164 is found to act as a general base during the formation of the acyl-enzyme from the tetrahedral intermediate. The role of dynamic substrate assisted catalysis (DSAC) involving His at the P1 site of the substrate is examined for the acylation reaction. It is demonstrated that the bond-breaking and -making events at each stage of the reaction trigger a change of the position for the His side chain and lead to the formation of the alternative hydrogen bonds. The back and forth movements of the His side chain between the C=O group of Pro at P2 and Odelta2 of Asp164 in a ping-pong-like mechanism and the formation of the alternative hydrogen bonds effectively lower the free energy barriers for both the nucleophilic attack and the acyl-enzyme formation and may therefore contribute to the relatively high activity of kumamolisin-As toward the substrates with His at the P1 site.

Nonequilibrium Solid Phase Formation Studied by Lattice Dynamics Calculation and Ion Beam Mixing in an Immiscible Co−Ag System

For the equilibrium immiscible Co-Ag system, a proven realistic ab initio derived n-body potential is applied to study the nonequilibrium solid phase formation at three chemical stoichiometries of Co/Ag = 1:3, 1:1, and 3:1. To predict the structural stability, the elastic constants and the phonon spectra are calculated at the chosen stoichiometries with a total of eight hypothetical crystalline structures. The calculated results suggest that four compounds, that is, D0(3) CoAg3, B1 CoAg, B2 CoAg, and D0(3) Co3Ag, are unstable, as they all feature negative elastic constants as well as imaginary phonons, and that another four compounds of both fcc-type L1(2) and hcp-type D0(19) structures at chemical stoichiometries of Co/Ag = 1:3 and 3:1, respectively, may elastically be favored and therefore obtainable under some specific conditions. It is also found that all the calculated elastic constants and phonon spectra are coincident within the framework of the elastic theory. Moreover, the calculated elastic constants are in good agreement with those acquired directly from ab initio calculations, lending support to the validity of the ab initio derived n-body Co-Ag potential as well as its resultant elastic constants and the phonon spectra. Interestingly, some of the predicted nonequilibrium solid phases, that is, two hcp-type compounds at chemical stoichiometries of Co/Ag = 1:3 and 3:1, respectively, are indeed obtained in ion beam mixing experiments and their lattice constants determined by diffraction analysis are in good agreement with those from calculations.

N-acetylglucosaminyltransferase V modifies the signaling pathway of epidermal growth factor receptor

Transfection of sense cDNA of N-acetylglucosamyltransferase V (GnTV-S) into human H7721 hepatocarcinoma cells resulted in an increase in the N-acetylglucosaminebeta1,6mannosealpha1,3- branch (GnT-V product) on the N-glycans of epidermal growth factor (EGF) receptor (EGFR), and promotion of its EGF binding and tyrosine autophosphorylation, but showed little effect on the expression of EGFR protein. The phosphorylation at T308, S473 and tyrosine residue(s) and the activity of protein kinase B (Akt/PKB) as well as the phosphorylation of p42/44 mitogen-activated protein kinase (MAPK) and MAPK kinase (MEK) before and after EGF stimulation were concomitantly increased. Conversely, in the antisense GnT-V (GnTV-AS)-transfected H7721 cells, all the results were the reverse of those with GnTV-S-transfected cells. After the cells were treated with 1-deoxymannojirimycin, an inhibitor of N-glycan processing at high mannose, or antibody against the extracellular glycan domain of EGFR, the differences in PKB activity, p42/44 MAPK and MEK phosphorylation among GnTV-S-, GnTV-AS- and mock-transfected cells were significantly attenuated. These findings indicate that the altered expression of GnT-V will change the glycan structure and function of EGFR, which may modify downstream signal transduction.

Relationship between metastasis-associated phenotypes and N-glycan structure of surface glycoproteins in human hepatocarcinoma cells

PURPOSE

To study the relation of N-glycan structure on cell surface glycoproteins to the metastatic phenotypes.

METHODS

Two human hepatocarcinoma 7721 cell lines transfected with sense or antisense cDNA of GnT-V, named GnT-V/7721 and GnT-V-AS/7721, respectively, were adopted, because GnT-V can change the antennary number and the content of the beta 1,6 GlcNAc branch in N-glycans. The effects of over- and under-expression of GnT-V on the metastasis-associated phenotype of the transfected cells were investigated and compared with the cells mock-transfected with the plasmid vector.

RESULTS

In GnT-V/7721 cells, GnT-V activity was increased by 92% compared with the mock cells. HRP-labeled lectin staining of transfected cells showed elevated intensity with HRP-L-PHA and reduced intensity with HRP-ConA, suggesting the increased antennary number and content of the beta 1,6 GlcNAc branch in N-glycans. Analysis of the N-glycan structure of [3H]-labeled glycopeptides prepared from cell-surface [3H] glycoproteins using DSA-affinity chromatography also revealed the above change of the N-glycan structure in a more quantitative manner. GnT-V/7721 cells showed a suppressed cell attachment to fibronectin (Fn) or laminin (Ln), and increased cell migration and invasion through matrigel. In contrast, GnT-V-AS/7721 cells showed reduction of both GnT-V activity and content of the beta 1,6 branch in N-linked glycans, elevation of cell attachment to Fn or Ln, and decline of cell migration and invasion through matrigel. These changes were just the opposite to those in GnT-V/7721 cells.

CONCLUSIONS

The alteration of N-glycan structure in surface glycoproteins resulting from the activity change of GnT-V contributes to the alterations in metastasis-associated phenotypes. The product of GnT-V, the beta 1,6 GlcNAc branch in N-linked glycans, is a structural factor of adhesion inhibition and invasion promotion. GnT-V is, therefore, closely related to cancer metastasis and its over-expression is an important molecular mechanism of metastasis.

Down-regulation of N-acetylglucosaminyltransferase V by tumorigenesis- or metastasis-suppressor gene and its relation to metastatic potential of human hepatocarcinoma cells

The effects of transfection of the metastasis suppressor gene nm23-H1 and cell-cycle related tumor-suppressor gene p16 on the activity of N-acetylglucosaminyltransferase V (GnT-V) and their relations to cancer metastatic potential were investigated. After transfection of nm23-H1 into 7721 human hepatocarcinoma cells and A549 human lung cancer cells, the activities of GnT-V were decreased by 28%-42% in the cells. In contrast, when p16 was transfected into these two cell lines, the decrease of GnT-V activity was only observed in A549 cells. This was probably to be due to the obvious expression of p16 gene in parental 7721 cells and the deletion of p16 in A549 cells. The decrease of GnT-V mRNA was only observed in nm23-H1-transfected cells, but not in p16-transfected A549 cells, suggesting that these two genes regulated GnT-V via different mechanisms. Horseradish peroxidase (HRP)-lectin staining showed that the 7721 cells transfected with nm23-H1 or the A549 cells transfected with p16 displayed a decreased intensity with HRP-leucoagglutinating phytohemagglutinin and increased intensity with HRP-concanavalin A, indicating the decline of beta1,6 N-acetylglucosamine branching structure on the asparagine-linked glycans of cell-surface and intracellular glycoproteins. The nm23-H1 transfected 7721 cells also displayed some changes in metastasis-related phenotypes, including the increase in cell adhesion to fibronectin (Fn), the decline in cell adhesion to laminin (Ln), and the decreased cell migration and invasion through matrigel. Transfection of antisense GnT-V cDNA into 7721 cells resulted in a decrease of GnT-V activity, an increase of cell adhesion to Fn or Ln, and a decrease in cell migration and invasion through matrigel. These phenotypes bore similarity to those of the 7721 cells transfected with nm23-H1. Our findings indicate that the down-regulation of GnT-V by nm23-H1 contributes to the alterations in metastasis-related phenotypes, and is an important molecular mechanism of metastasis suppression mediated by nm23-H1.

Modulation of the basal activity of phosphatidylinositol-3-kinase/protein kinase B signaling pathway in human hepatocarcinoma cells

The modulation of GnT-V activity by signaling molecules in PI-3-K/PKB pathway in human hepatocarcinoma cell line 7721 was studied. GnT-V activity was determined after the transfection of sense or antisense cDNA of PKB into the cells, as well as the addition of activators, specific inhibitors, and the antibodies to the enzyme assay system or culture medium. It was found that the basal activity of GnT-V was up regulated by the sense and down regulated by the antisense cDNA of PKB transfected into 7721 cells. GnT-V was activated by PIP2, PIP3 or GTPgamma[S] added to the assay system, and the activation of PIP2 or GTPgamma[S] was abolished by LY2940002, a specific inhibitor of PI-3-K, but the activation of PIP3 was not attenuated by LY2940002. In addition, GnT-V activity in cultured parental or H-ras transfected cells was inhibited by the antibody against PKB or PI-3-K. These findings demonstrated the involvement of PI-3-K/PKB signaling pathway in the regulation of GnT-V. Moreover, ET18-OCH3, an inhibitor of Raf translocation and PI-PLC enzyme, which produces the activator of PKC, as well as the antibodies against Raf-1 or MEK also inhibited GnT-V activity in the parental and H-ras transfected cells. The inhibitory rates, however, were less in the transfected cells than those in the parental cells. These results reveal that in parental and H-ras transfected 7721 cells, the basal activity of GnT-V is also regulated by the Ras/Raf-1/MEK/MAPK cascade in addition to PI-3-K/PKB signaling pathway. The significance of these two pathways in the regulation of GnT-V and their relations to the activation of PKC previously reported by our laboratory (Ju TZ et al., 1995 Glyconjugate J 12, 767-772) was discussed.

Effects of H-ras and v-sis overexpression on N-acetylglucosaminyltransferase V and metastasis-related phenotypes in human hepatocarcinoma cells

Oncogenes and N-acetylglucosaminyltransferase (GnT-V) are both commonly associated with carcinogenesis and metastasis. In order to elucidate the relationship between oncogenes and GnT-V, two oncogenes, H-ras and v-sis/PDGF (platelet-derived growth factor), were selected, and the effects of their overexpression on GnT-V in 7721 human hepatocarcinoma cells were investigated. The results showed that the over expression of H-ras or v-sis/PDGF-B up-regulated the activities of GnT-V to various degrees in the transfected cells. In H-ras- and PDGF-B-overexpressing cells, the activity of GnT-V was up-regulated to double the normal value. The transient expression of v-sis, which produces a protein almost identical to PDGF-B, stimulated the GnT-V activity by 80.3%, and the effect was more pronounced (increased by 182.5%) in 7721 cells with stable expression of v-sis. The stimulating effect was entirely abolished by treatment with PDGF-B antibody. The staining of asparagine-linked glycans (N-glycans) in the H-ras- and v-sis-overexpressing 7721 cells was intensified when horseradish peroxidase-labeled leucoagglutinating phytohemogglutinin was used as a probe, indicating the increased content of beta1,6GlcNAc branching on the N-glycans. The enhancement of GnT-V mRNA expression was also observed in H-ras- and v-sis- overexpressing cells, indicating that H-ras and v-sis regulated GnT-V via the transcription of GnT-V mRNA and the synthesis of GnT-V protein. The cells overexpressing H-ras and v-sis displayed some changes in metastasis-related phenotypes, including acceleration of cell growth, decline of cell adhesion to fibronectin, and an increase of cell adhesion to laminin, as well as increased invasiveness through Matrigel. These results indicated that the alteration of cell adhesion and invasion induced by oncogenes is closely related to the up-regulation of GnT-V activity and its product, beta1,6GlcNAc branching in N-glycans on the cell surface.

Opposing changes in N-acetylglucosaminyltransferase-V and -III during the cell cycle and all-trans retinoic acid treatment of hepatocarcinoma cell line

The changes in N-acetylglucosaminyltransferase-V and -III (GnT-V, GnT-III) during the cell-cycle of synchronized 7721 human hepatocarcinoma cell line were investigated. Using an HPLC method to assay GnT and flow cytometry (FCM) for cell cycle analysis, it was found that GnT-V showed the highest activity, but GnT-III reached the lowest activity when G(2)/M cells were most abundant. In contrast, GnT-V declined to the minimum while GnT-III elevated to maximum when G(0)/G(1) cells were most predominant. The opposing changes were more obvious when the activities of GnT-V and GnT-III were expressed as relative activities (activity of GnT-V or GnT-III/the sum of activities of GnT-V plus GnT-IV plus GnT-III). These opposing changes of GnT-V and GnT-III during the cell cycle might result from the different regulatory mechanisms of GnT-V and GnT-III expression in the cell cycle. The alterations in the structures of cell surface N-glycans were compatible with the changes of the activities of GnTs. The results from immunocytochemistry and Northern blot showed that the protein and mRNA contents of GnT-V were not significantly changed during the cell cycle. The activity of a cell cycle regulating protein kinase, p34(cdc2) kinase, correlated to the activity of GnT-V. These findings suggested that the change of GnT-V activity in cell cycle was not the consequence of the alteration of gene transcription or enzyme protein synthesis, but might be caused by the post-translational regulation. The decrease in GnT-V and the corresponding increase in GnT-III activities were also found after the cells were treated with all-trans retinoic acid (ATRA), and the mechanism of this might be different from that in the cell cycle.

Increased susceptibility to apoptosis of human hepatocarcinoma cells transfected with antisense N-acetylglucosaminyltransferase V cDNA

The antisense cDNA of N-acetylglucosaminyltransferase V (GnT-V, EC 2. 4.1.155) was constructed as pcDNA3/GnT-V-AS plasmid and transfected into 7721 cells, a human hepatocarcinoma cell line. The transfection was confirmed with Northern blot. By using HPLC and HRP-lectin staining, it was found that the cells transfected with pcDNA3/GnT-V-AS (GnT-V-AS/7721) expressed less GnT-V activity and beta-1,6-GlcNAc branching in the cell glycoproteins compared with the cells mock-transfected with the vector pcDNA3 (pcDNA3/7721). The growth rate of GnT-V-AS/7721 was decreased in serum-containing medium, while the cell death was accelerated in serum-free medium. The GnT-V-AS/7721 cells were more susceptible to the apoptosis induced by ATRA than the mock-transfected cells. This was evidenced by the obvious appearance of a hypoploid sub-G(1) fraction in the DNA histogram using FCM analysis, the more condensed new moon-type nuclei under morphological observation, and the more intensive TUNEL reaction for assaying the fragmented DNA. At the same time as GnT-V down-regulation by GnT-V-AS, an increase of another N-aceylglusaminyltransferase, GnT-III (EC 2.4.1.144), was observed, and the biological significance of this finding was discussed.

[Changes with cell cycle of N-acetylglucosaminyl transferase III, IV and V in 7721 human hepatocarcinoma cells]

In order to investigate the changes of N-acetylglucosaminyl transferase (GlcNAc-T) III, IV and V in cell cycle, the synchronization of 7721 human hepatocellular carcinoma cells was performed using serum hunger method. The percentages of cells in different phases during cell cycle were measured by flow cytometry (FCM) and the cell cycle was checked by determining the activity of cellular p34cdc2 kinase. It was found that the activities of GlcNAc-T III increased in G0/G1 cell peak phase and had correlation with the cell percentage of G0/G1 phase (r = 0.760, P < 0.05), while GlcNAc-T V showed the highest activity when G2/M cells were most abundant and had an apparent correlation with the cell percentage of G2/M phase (r = 0.868, P < 0.001). The changes of GlcNAc-T IV activity seemed not related to the cell cycle. The changes in opposite directions of relative activities (percentage of total GlcNAc-T III, IV, V) of GlcNAc-T III and GlcNAc-T V were observed during cell cycle (r = -0.951, P < 0.001), suggesting that these two enzymes might be regulated differently and functioned oppositely in the cells: GlcNAc-T V may be related to the proliferation of 7721 cells, while GlcNAc-T III may be related to the non-mitotic silence phase of the cells, or, it may be a factor against proliferation. Immunohistochemical results showed that the protein content of GlcNAc-T V was not significantly changed during cell cycle, and had no correlation with the activity of GlcNAc-T V, suggesting that the changes of GlcNAc-T V activity in cell cycle might not be resulted from the alteration of enzyme protein synthesis. The correlation between the activities of GlcNAc-T V and p34cdc2 kinase (r = 0.752, P < 0.05) was observed in cell cycle, implicating that GlcNAc-T V might possibly be regulated by p34cdc2 kinase.

Clinical significance of serum S100 in metastatic malignant melanoma

The clinical significance of serum S100 was assessed in comparison to neuron-specific enolase (NSE) in 126 patients with malignant melanoma: 80 patients with clinical stage I/II, 23 patients with stage III and 23 patients with stage IV according to the criteria of the American Joint Committee on Cancer (AJCC). Using cut-off values of 0.15 microgram/l for S100 and 12.5 micrograms/l for NSE, the sensitivity was found to be 1.3% (1/80) for S100 and 8.75% (7/80) for NSE in patients with stage I/II, 8.7% (2/23) for S100 and 13% (8/23) for NSE in patients with stage III, and 73.9% (17/23) for S100 and 34.8% (8/23) for NSE in patients with stage IV disease (P < 0.05). In 6 patients with stage III/IV tumours, serial measurement of serum S100 and NSE was performed. A rise of serum S100 indicated progression of the disease; a decline indicated response to treatment. Our preliminary results support the value of serum S100 as an adjunct to the clinical staging and monitoring of metastatic malignant melanoma.

Clinical significance of serum S100 in metastatic malignant melanoma

The clinical significance of serum S100 was assessed in comparison to neuron-specific enolase (NSE) in 126 patients with malignant melanoma: 80 patients with clinical stage I/II, 23 patients with stage III and 23 patients with stage IV according to the criteria of the American Joint Committee on Cancer (AJCC). Using cut-off values of 0.15 micrograms/l for S100 and 12.5 micrograms/l for NSE, the sensitivity was found to be 1.3% (1/80) for S100 and 8.75% (7/80) for NSE in patients with stage I/II, 8.7% (2/23) for S100 and 13% (8/23) for NSE in patients with stage III, and 73.9% (17/23) for S100 and 34.8% (8/23) for NSE in patients with stage IV disease (P < 0.05). In 6 patients with stage III/IV tumours, serial measurement of serum S100 and NSE was performed. A rise of serum S100 indicated progression of the disease; a decline indicated response to treatment. Our preliminary results support the value of serum S100 as an adjunct to the clinical staging and monitoring of metastatic malignant melanoma.

[Leptotrombidium (L.) scutellare as the transmitting vector of tsutsugamushi disease of autumn-winter type in Jiangsu Province]

In 1986, epidemic of tsutsugamushi disease of autumn-winter type was found in Jiangsu province. To clarify the vector of this disease, we carried out a series of studies in 1986-1992. Leptotrombidium (L.) scutellare was found to be a dominant species of chigger mite on rats in the endemic areas and its seasonal distribution was correlated with the incidence of tsutsugamushi disease in inhabitants. This mite could naturally be infected by Rickettsia tsutsugamushi, and R. tsutsugamushi could be transmitted via biting and transovarial transmission. Specific antibodies could be detected in the sera of mice bitten by the mites or inoculated with the suspension of mites. Serological typing of the sera of mice was of Gilliam type. The above results demonstrate that L. (L.) scutellare can serve as transmitting vector of tsutsugamushi disease of autumn-winter type.