Intestinal microbiota imbalance resulted by anti-Toxoplasma gondii immune responses aggravate gut and brain injury.

BACKGROUND: Toxoplasma gondii infection affects a significant portion of the global population, leading to severe toxoplasmosis and, in immunocompromised patients, even death. During T. gondii infection, disruption of gut microbiota further exacerbates the damage to intestinal and brain barriers. Therefore, identifying imbalanced probiotics during infection and restoring their equilibrium can regulate the balance of gut microbiota metabolites, thereby alleviating tissue damage. METHODS: Vimentin gene knockout (vim-/-) mice were employed as an immunocompromised model to evaluate the influence of host immune responses on gut microbiota balance during T. gondii infection. Behavioral experiments were performed to assess changes in cognitive levels and depressive tendencies between chronically infected vim-/- and wild-type (WT) mice. Fecal samples were subjected to 16S ribosomal RNA (rRNA) sequencing, and serum metabolites were analyzed to identify potential gut probiotics and their metabolites for the treatment of T. gondii infection. RESULTS: Compared to the immunocompetent WT sv129 mice, the immunocompromised mice exhibited lower levels of neuronal apoptosis and fewer neurobehavioral abnormalities during chronic infection. 16S rRNA sequencing revealed a significant decrease in the abundance of probiotics, including several species of Lactobacillus, in WT mice. Restoring this balance through the administration of Lactobacillus murinus and Lactobacillus gasseri significantly suppressed the T. gondii burden in the intestine, liver, and brain. Moreover, transplantation of these two Lactobacillus spp. significantly improved intestinal barrier damage and alleviated inflammation and neuronal apoptosis in the central nervous system. Metabolite detection studies revealed that the levels of various Lactobacillus-related metabolites, including indole-3-lactic acid (ILA) in serum, decreased significantly after T. gondii infection. We confirmed that L. gasseri secreted much more ILA than L. murinus. Notably, ILA can activate the aromatic hydrocarbon receptor signaling pathway in intestinal epithelial cells, promoting the activation of CD8+ T cells and the secretion of interferon-gamma. CONCLUSION: Our study revealed that host immune responses against T. gondii infection severely disrupted the balance of gut microbiota, resulting in intestinal and brain damage. Lactobacillus spp. play a crucial role in immune regulation, and the metabolite ILA is a promising therapeutic compound for efficient and safe treatment of T. gondii infection.

Comparative effects of different metals on the Japanese medaka embryos and larvae.

Rapid evaluation of the toxicity of metals using fish embryo acute toxicity is facilitative to ecological risk assessment of aquatic organisms. However, this approach has seldom been utilized for the comparative study on the effects of different metals to fish. In this study, acute and sub-chronic tests were used to compare the toxicity of Se(IV) and Cd in the embryos and larvae of Japanese medaka (Oryzias latipes). The embryos with different levels of dechorionation and/or pre-exposure were also exposed to Se(IV) and Cd at various concentrations. The results showed that the LC50-144 h of Cd was 1.3-5.2 folds higher than that of Se(IV) for the embryos. In contrast, LC50-96 h of Se(IV) were 200-400 folds higher than that of Cd for the larvae. Meanwhile, dechorionated embryos were more sensitive to both Se and Cd than the intact embryos. At elevated concentrations, both Se and Cd caused mortality and deformity in the embryos and larvae. In addition, pre-exposure to Cd at the embryonic stages enhanced the resistance to Cd in the larvae. However, pre-exposure to Se(IV) at the embryonic stages did not affect the toxicity of Se(IV) to the larvae. This study has distinguished the nuance differences in effects between Se(IV) and Cd after acute and sub-chronic exposures with/without chorion. The approach might have a potential in the comparative toxicology of metals (or other pollutants) and in the assessment of their risks to aquatic ecosystems.

Population-based cohort study of Toxoplasma gondii P22 antibody positivity correlation with anxiety.

BACKGROUND: Accumulating evidence suggests that latent infection with Toxoplasma gondii (T. gondii) is associated with a variety of neuropsychiatric and behavioral conditions. This research aims to explore the potential correlation between T. gondii antibody positivity and neuropsychiatric disorders through a comprehensive prospective cohort study. METHODS: The cohort study utilized the UK Biobank database to recruit 8814 individuals with no prior diagnosis of neuropsychiatric disorders. Cox proportional hazards models were employed to investigate the associations between T. gondii P22 antibody seropositivity (P22+) and the development of various types of neuropsychiatric disorders. RESULTS: Of the population, 14.65 % tested positive for T. gondii P22 antibody. The presence of T. gondii P22 antibody showed a slight inverse association with epilepsy (HR: 0.28; 95 % CI: 0.10-0.77), while it was positively associated with an increased risk of developing anxiety disorders (HR: 1.38; 95 % CI: 1.04-1.83). LIMITATIONS: The study sample consisted mostly of white British individuals aged 40 to 69 years old. Although we adjusted for potential confounders, there may be other unmeasured and residual confounding factors that could have influenced our reported associations. CONCLUSIONS: The findings suggested an increased risk of anxiety and potential evidence of epilepsy associated with T. gondii P22+. However, our analysis did not reveal an increased risk of several other neuropsychiatric conditions including Alzheimer's disease, dementia, substance abuse disorders, depression, and neurodegenerative disorders, associated with P22 antibody seropositivity.

The Idesia polycarpa genome provides insights into its evolution and oil biosynthesis.

The deciduous tree Idesia polycarpa can provide premium edible oil with high polyunsaturated fatty acid contents. Here, we generate its high-quality reference genome, which is ∼1.21 Gb, comprising 21 pseudochromosomes and 42,086 protein-coding genes. Phylogenetic and genomic synteny analyses show that it diverged with Populus trichocarpa about 16.28 million years ago. Notably, most fatty acid biosynthesis genes are not only increased in number in its genome but are also highly expressed in the fruits. Moreover, we identify, through genome-wide association analysis and RNA sequencing, the I. polycarpa SUGAR TRANSPORTER 5 (IpSTP5) gene as a positive regulator of high oil accumulation in the fruits. Silencing of IpSTP5 by virus-induced gene silencing causes a significant reduction of oil content in the fruits, suggesting it has the potential to be used as a molecular marker to breed the high-oil-content cultivars. Our results collectively lay the foundation for breeding the elite cultivars of I. polycarpa.

RNA-Seq-Based WGCNA and Association Analysis Reveal the Key Regulatory Module and Genes Responding to Salt Stress in Wheat Roots.

Soil salinization is the main abiotic stressor faced by crops. An improved understanding of the transcriptional response to salt stress in roots, the organ directly exposed to a high salinity environment, can inform breeding strategies to enhance tolerance and increase crop yield. Here, RNA-sequencing was performed on the roots of salt-tolerant wheat breeding line CH7034 at 0, 1, 6, 24, and 48 h after NaCl treatment. Based on transcriptome data, a weighted gene co-expression network analysis (WGCNA) was constructed, and five gene co-expression modules were obtained, of which the blue module was correlated with the time course of salt stress at 1 and 48 h. Two GO terms containing 249 differentially expressed genes (DEGs) related to osmotic stress response and salt-stress response were enriched in the blue module. These DEGs were subsequently used for association analysis with a set of wheat germplasm resources, and the results showed that four genes, namely a Walls Are Thin 1-related gene (TaWAT), an aquaporin gene (TaAQP), a glutathione S-transfer gene (TaGST), and a zinc finger gene (TaZFP), were associated with the root salt-tolerance phenotype. Using the four candidate genes as hub genes, a co-expression network was constructed with another 20 DEGs with edge weights greater than 0.6. The network showed that TaWAT and TaAQP were mainly co-expressed with fifteen interacting DEGs 1 h after salt treatment, while TaGST and TaZFP were mainly co-expressed with five interacting DEGs 48 h after salt treatment. This study provides key modules and candidate genes for understanding the salt-stress response mechanism in wheat roots.

Comparative study of rheological properties and Pickering emulsion stabilizing capacity of nonenyl succinic anhydride and octenyl succinic anhydride modified amaranth starches.

Functional properties and ability to stabilize Pickering emulsions of amaranth starch with the novel nonenyl succinic anhydride (NSA) modification and the widely used octenyl succinic anhydride (OSA) modification were compared. The NSA modification was more effective in altering the rheological properties of amaranth starches. NSA-modified amaranth starch showed significantly higher peak viscosity (7.13 Pa·s at DS of 0.02209) than the OSA-modified amaranth starch (6.10 Pa·s at DS of 0.03042). The gelatinization temperature, gelatinization enthalpy, and relative crystallinity of amaranth starch were more affected by the OSA than the NSA. The Pickering emulsions stabilized with NSA-modified starches had higher stability than those with the OSA-modified starches as characterized by particle size distribution, morphological, and rheological approaches. A lower degree of substitution by NSA than by OSA is needed to achieve a similar emulsification capacity. Thus, the NSA modification could be an efficient alternative to OSA modification in tailoring physicochemical and rheological functions, as well as stabilizing Pickering emulsions.

Toxoplasma gondii gra5 deletion mutant protects hosts against Toxoplasma gondii infection and breast tumors.

Toxoplasma gondii is the causative agent of toxoplasmosis, a zoonotic disease that poses a threat to human health and a considerable loss to livestock farming. At present, clinical therapeutic drugs mainly target T. gondii tachyzoites and fail to eradicate bradyzoites. Developing a safe and effective vaccine against toxoplasmosis is urgent and important. Breast cancer has become a major public health problem and the therapeutic method needs to be further explored. Many similarities exist between the immune responses caused by T. gondii infection and the immunotherapy for cancers. T. gondii dense granule organelles secrete immunogenic dense granule proteins (GRAs). GRA5 is localized to the parasitophorous vacuole membrane in the tachyzoite stage and the cyst wall in the bradyzoite stage. We found that T. gondii ME49 gra5 knockout strain (ME49Δgra5) was avirulent and failed to form cysts but stimulated antibodies, inflammatory cytokines, and leukocytes infiltration in mice. We next investigated the protective efficacy of ME49Δgra5 vaccination against T. gondii infection and tumor development. All the immunized mice survived the challenge infection of either wild-type RH, ME49, VEG tachyzoites, or ME49 cysts. Moreover, ME49Δgra5 tachyzoite inoculation in situ attenuated the growth of murine breast tumor (4T1) in mice and prevented 4T1's lung metastasis. ME49Δgra5 inoculation upregulated the levels of Th1 cytokines and tumor-infiltrating T cells in the tumor microenvironment and triggered anti-tumor responses by increasing the number of natural killer, B, and T cells, macrophages, and dendritic cells in the spleen. Collectively, these results suggested that ME49Δgra5 was a potent live attenuated vaccine against T. gondii infection and breast cancer.

Genome-wide identification of YABBY genes in three Cymbidium species and expression patterns in C. ensifolium (Orchidaceae).

Members of the YABBY gene family play significant roles in lamina development in cotyledons, floral organs, and other lateral organs. The Orchidaceae family is one of the largest angiosperm groups. Some YABBYs have been reported in Orchidaceae. However, the function of YABBY genes in Cymbidium is currently unknown. In this study, 24 YABBY genes were identified in Cymbidium ensifolium, C. goeringii, and C. sinense. We analyzed the conserved domains and motifs, the phylogenetic relationships, chromosome distribution, collinear correlation, and cis-elements of these three species. We also analyzed expression patterns of C. ensifolium and C. goeringii. Phylogenetic relationships analysis indicated that 24 YABBY genes were clustered in four groups, INO, CRC/DL, YAB2, and YAB3/FIL. For most YABBY genes, the zinc finger domain was located near the N-terminus and the helix-loop-helix domain (YABBY domain) near the C-terminus. Chromosomal location analysis results suggested that only C. goeringii YABBY has tandem repeat genes. Almost all the YABBY genes displayed corresponding one-to-one relationships in the syntenic relationships analysis. Cis-elements analysis indicated that most elements were clustered in light-responsive elements, followed by MeJA-responsive elements. Expression patterns showed that YAB2 genes have high expression in floral organs. RT-qPCR analysis showed high expression of CeYAB3 in lip, petal, and in the gynostemium. CeCRC and CeYAB2.2 were highly expressed in gynostemium. These findings provide valuable information of YABBY genes in Cymbidium species and the function in Orchidaceae.

Genome-wide identification and analysis of anthocyanin synthesis-related R2R3-MYB genes in Cymbidium goeringii.

The MYB gene family plays a vital regulatory role in plant metabolism, stress response, and floral color. The R2R3-MYB gene family of C. goeringii was identified, and its expression was analyzed using bioinformatics in this article. The R2R3-MYB genes of Arabidopsis thaliana were used as a reference to determine 104 CgMYB genes and categorize them into 22 subfamilies. Exon/intron organizations and conserved motif analysis revealed that the majority of CgMYB genes were highly conserved, and chromosome localization and collinearity analysis provided evidence of tandem duplication and segmental duplication events, indicating the phenomenon of gene family expansion and contraction. The function of CgMYB genes was analyzed by cis-acting element and gene ontology (GO) enrichment. In addition, we selected CgMYB91 and CgMYB32 for RT-qPCR, suggesting that CgMYB91 and CgMYB32 are associated with anthocyanin formation. In short, this study provides a comprehensive and specific function of the R2R3-MYB transcription factors (TFs) in orchids.

Toxoplasma gondii infection possibly reverses host immunosuppression to restrain tumor growth.

Tumor cells can successfully escape the host immune attack by inducing the production of immunosuppressive cells and molecules, leading to an ineffective tumor treatment and poor prognosis. Although immunotherapies have improved the survival rate of cancer patients in recent years, more effective drugs and therapies still need to be developed. As an intracellular parasite, Toxoplasma gondii can trigger a strong Th1 immune response in host cells, including upregulating the expression of interleukin-12 (IL-12) and interferon-γ (IFN-γ). Non-replicating uracil auxotrophic strains of T. gondii were used to safely reverse the immunosuppression manipulated by the tumor microenvironment. In addition to the whole lysate antigens, T. gondii-secreted effectors, including Toxoplasma profilin, rhoptry proteins (ROPs), and dense granule antigens (GRAs), are involved in arousing the host's antigen presentation system to suppress tumors. When T. gondii infection relieves immunosuppression, tumor-related myeloid cells, including macrophages and dendritic cells (DCs), are transformed into immunostimulatory phenotypes, showing a powerful Th1 immune response mediated by CD8+ T cells. Afterwards, they target and kill the tumor cells, and ultimately reduce the size and weight of tumor tissues. This article reviews the latest applications of T. gondii in tumor therapy, including the activation of cellular immunity and the related signal pathways, which will help us understand why T. gondii infection can restrain tumor growth.

Toxoplasma gondii DNA methyltransferases regulate parasitic energy metabolism.

The intracellular protozoan Toxoplasma gondii results in serious diseases such as encephalitis, and retinochoroiditis in immunocompromised patients. The interconversion between tachyzoites and bradyzoites under the host's immune pressure results in the interchange of acute infection and chronic infection. We previously reported two functional DNA methyltransferases (DNMT) in Toxoplasma gondii named TgDNMTa and TgDNMTb. In this research, proteomics analysis for T. gondii tachyzoites of ME49 WT, dnmta knockout (ME49-∆Tgdnmta), and dnmtb knockout (ME49-∆Tgdnmtb) strains, revealed 362 significantly regulated proteins for ME49-∆Tgdnmta, and 219 for ME49-∆Tgdnmtb, compared with the proteins of ME49 WT. TgDNMTa down regulated three glycolytic enzymes, one gluconeogenic enzyme and four pyruvate metabolic enzymes. Furthermore, TgDNMTb up regulated two proteins in the tricarboxylic acid (TCA) cycle. Glucose metabolic flux detection showed that TgDNMTa inhibited the glycolysis pathway, while TgDNMTb promoted the tricarboxylic acid (TCA) cycle so as to promote parasite's proliferation. These findings demonstrated that the functions of Toxoplasma gondii DNA methyltransferases extended beyond DNA methylation to the regulation of parasitic energy metabolism.

Genome-Wide Identification and Expression Analysis of Terpene Synthase Genes in Cymbidium faberi.

Terpene synthases (TPSs) are essential for forming terpenes, which play numerous functional roles in attracting pollinators, defending plants, and moderating the interaction between plants. TPSs have been reported in some orchids, but genome-wide identification of terpenes in Cymbidium faberi is still lacking. In this study, 32 putative TPS genes were classified in C. faberi and divided into three subfamilies (TPS-a, TPS-b, and TPS-e/f). Motif and gene structure analysis revealed that most CfTPS genes had the conserved aspartate-rich DDxxD motif. TPS genes in the TPS-a and TPS-b subfamilies had variations in the RRX8W motif. Most cis-elements of CfTPS genes were found in the phytohormone responsiveness category, and MYC contained most of the numbers associated with MeJA responsiveness. The Ka/Ks ratios of 12/13 CfTPS gene pairs were less than one, indicated that most CfTPS genes have undergone negative selection. The tissue-specific expression patterns showed that 28 genes were expressed in at least one tissue in C. faberi, and TPS genes were most highly expressed in flowers, followed by leaves and pseudobulbs. In addition, four CfTPS genes were selected for the real-time reverse transcription quantitative PCR (RT-qPCR) experiment. The results revealed that CfTPS12, CfTPS18, CfTPS23, and CfTPS28 were mainly expressed in the full flowering stage. CfTPS18 could convert GPP to ß-myrcene, geraniol, and α-pinene in vitro. These findings of CfTPS genes of C. faberi may provide valuable information for further studies on TPSs in orchids.

Toxoplasma gondii Type-I ROP18 Targeting Human E3 Ligase TRIM21 for Immune Escape.

Toxoplasma gondii is an intracellular pathogen that exerts its virulence through inhibiting host's innate immune responses, which is mainly related to the type II interferon (IFN-γ) response. IFN-γ inducible tripartite motif 21 (TRIM21), an E3 ligase, plays an important role in anti-infection responses against the intracellular pathogens including bacteria, virus, and parasite. We found that T. gondii virulence factor ROP18 of the type I RH strain (TgROP18I) interacted with human TRIM21, and promoted the latter's phosphorylation, which subsequently accelerated TRIM21 degradation through lysosomal pathway. Furthermore, TRIM21 protein level was found to be upregulated during RH and CEP strains of T. gondii infection. TRIM21 knocking down reduced the ubiquitin labeling on the parasitophorous vacuole membrane (PVM) [which led to parasitophorous vacuole (PV) acidification and death of CEP tachyzoites], and relieved the inhibition of CEP proliferation induced by IFN-γ in human foreskin fibroblast (HFF) cells which was consistent with the result of TRIM21 overexpression. On the other hand, TRIM21 overexpression enhanced the inhibition of CEP proliferation, and inhibited the binding of IκB-α with p65 to activate the IFN-γ-inducible NF-κB pathway, which might be resulted by TRIM21-IκB-α interaction. In brief, our research identified that in human cells, IFN-γ-inducible TRIM21 functioned in the innate immune responses against type III T. gondii infection; however, TgROP18I promoted TRIM21 phosphorylation, leading to TRIM21 degradation for immune escape in type I strain infection.

A Review on Dengue Vaccine Development.

Dengue virus (DENV) has become a global health threat with about half of the world's population at risk of infection. Although the disease caused by DENV is self-limiting in the first infection, the antibody-dependent enhancement (ADE) effect increases the mortality in the second infection with a heterotypic virus. Since there is no specific efficient medicine in treatment, it is urgent to develop vaccines to prevent infection and disease progression. Currently, only a live attenuated vaccine, chimeric yellow fever 17D-tetravalent dengue vaccine (CYD-TDV), has been licensed for clinical use in some countries, and many candidate vaccines are still under research and development. This review discusses the progress, strengths, and weaknesses of the five types of vaccines including live attenuated vaccine, inactivated virus vaccine, recombinant subunit vaccine, viral vectored vaccine, and DNA vaccine.

Application of the Scorpion Neurotoxin AaIT against Insect Pests.

Androctonus australis Hector insect toxin (AaIT), an insect-selective toxin, was identified in the venom of the scorpion Androctonus australis. The exclusive and specific target of the toxin is the voltage-gated sodium channels of the insect, resulting in fast excitatory paralysis and even death. Because of its strict toxic selectivity and high bioactivity, AaIT has been widely used in experiments exploring pest bio-control. Recombinant expression of AaIT in a baculovirus or a fungus can increase their virulence to insect pests and diseases vectors. Likewise, transgenic plants expressing AaIT have notable anti-insect activity. AaIT is an efficient toxin and has great potential to be used in the development of commercial insecticides.

Expression of Bacillus thuringiensis toxin Cyt2Ba in the entomopathogenic fungus Beauveria bassiana increases its virulence towards Aedes mosquitoes.

BACKGROUND: The entomopathogenic fungus Beauveria bassiana has been widely used to kill mosquito larvae and adults in the laboratory and field. However, its slow action of killing has hampered its widespread application. In our study, the B. bassiana fungus was genetically modified to express the Bacillus thuringiensis (Bt) toxin Cyt2Ba to improve its efficacy in killing mosquitoes. METHODOLOGY/PRINCIPAL FINDINGS: The efficacy of the wild type (WT) of B. bassiana and a transgenic strain expressing Cyt2Ba toxin (Bb-Cyt2Ba) was evaluated against larval and adult Aedes mosquitoes (Aedes aegypti and Aedes albopictus) using insect bioassays. The Bb-Cyt2Ba displayed increased virulence against larval and adult Aedes mosquitoes compared with the WT: for Ae. aegypti adults, the median lethal time (LT50) was decreased by 33% at the concentration of 1× 108 conidia/ml, 19% at 1× 107 conidia/ml and 47% at 1× 106 conidia/ml. The LT50 for Ae. albopictus adults was reduced by 20%, 23% and 29% at the same concentrations, respectively. The LT50 for Ae. aegypti larvae was decreased by 42% at 1× 107 conidia/ml and 25% at 1× 106 conidia/ml, and that for Ae. albopictus larvae was reduced by 33% and 31% at the same concentrations, respectively. In addition, infection with Bb-Cyt2Ba resulted in a dramatic reduction in the fecundity of Aedes mosquitoes. CONCLUSIONS/SIGNIFICANCE: In conclusion, our study demonstrated that the virulence of B. bassiana against mosquitoes can be significantly improved by introducing the Bt toxin gene Cyt2Ba into the genome to express the exogenous toxin in the fungus. The transgenic strain Bb-Cyt2Ba significantly reduced the survival and fecundity of Ae. aegypti and Ae. albopictus compared with the WT strain, which suggested that this recombinant B. bassiana has great potential for use in mosquito control.

[Larvicidal activity of recombinant Escherichia coli expressing scorpion neurotoxin AaIT or B.t.i toxin Cyt2Ba against mosquito larvae and formulations for enhancing the effects].

OBJECTIVE: To assess the larvicidal effects of recombinant Escherichia coli expressing scorpion neurotoxin AaIT or Bacillus thuringiensis subsp israelensis (B.t.i) toxin Cyt2Ba against the second instar larvae of Culex pipiensquinquefasciatus and Aedes albopictus and compare different formulations for their larvicidal effects. METHODS: The AaIT- or Cyt2Ba-coding sequences were cloned into pET28a(+) and the recombinant plasmids were transformed into E. coli BL21(DE3). After induction with IPTG, the recombinant proteins expressed by the recombinant E. coli were detected and identified by SDS-PAGE and Western blotting, respectively. The larvicidal activity of the bacterial suspension was tested at different concentrations against mosquitoes. The effective engineered bacteria were prepared into dry powder with different formulations, and their larvicidal activity was tested. RESULTS: AaIT and Cyt2Ba proteins were successfully expressed in E. coli. The recombinant AaIT protein showed no virulence to the mosquito larvae. The suspension of the recombinant E. coli expressing Cyt2Ba protein exhibited a stronger killing effect on Aedes albopictus larvae than on Culex pipiens quinquefasciatus larvae at 48 h (P<0.001) with LC50 of 3.00×106 cells/mL and 1.25×107 cells/mL, respectively. The dry powder of the engineered bacteria formulated with yeast extract, wheat flour or white pepper powder at the mass ratio of 1:1 showed the strongest killing effect on mosquito larvae (P=0.044), and the formulation with white pepper powder produced a stronger killing effect than formulations with yeast extract or wheat flour (P=0.002). CONCLUSION: The B.t.i Cyt2Ba protein expressed in E. coli BL21(DE3) shows a good larvicidal activity against mosquitoes, and appropriate formulations of the engineered bacteria can enhance its efficiency in mosquito control.

Towards Photo Watercolorization with Artistic Verisimilitude.

We present a novel artistic-verisimilitude driven system for watercolor rendering of images and photos. Our system achieves realistic simulation of a set of important characteristics of watercolor paintings that have not been well implemented before. Specifically, we designed several image filters to achieve: 1) watercolor-specified color transferring; 2) saliency-based level-of-detail drawing; 3) hand tremor effect due to human neural noise; and 4) an artistically controlled wet-in-wet effect in the border regions of different wet pigments. A user study indicates that our method can produce watercolor results of artistic verisimilitude better than previous filter-based or physical-based methods. Furthermore, our algorithm is efficient and can easily be parallelized, making it suitable for interactive image watercolorization.