In Silico Prediction of Skin Sensitization for Compounds via Flexible Evidence Combination Based on Machine Learning and Dempster-Shafer Theory.

Skin sensitization is increasingly becoming a significant concern in the development of drugs and cosmetics due to consumer safety and occupational health problems. In silico methods have emerged as alternatives to traditional in vivo animal testing due to ethical and economic considerations. In this study, machine learning methods were used to build quantitative structure-activity relationship (QSAR) models on five skin sensitization data sets (GPMT, LLNA, DPRA, KeratinoSens, and h-CLAT), achieving effective predictive accuracies (correct classification rates of 0.688-0.764 on test sets). To address the complex mechanisms of human skin sensitization, the Dempster-Shafer theory was applied to merge multiple QSAR models, resulting in an evidence-based integrated decision model. Various evidence combinations and combination rules were explored, with the self-defined Q3 rule showing superior balance. The combination of evidence such as GPMT and KeratinoSens and h-CLAT achieved a correct classification rate (CCR) of 0.880 and coverage of 0.893 while maintaining the competitiveness of other combinations. Additionally, the Shapley additive explanations (SHAP) method was used to interpret important features and substructures related to skin sensitization. A comparative analysis of an external human test set demonstrated the superior performance of the proposed method. Finally, to enhance accessibility, the workflow was implemented into a user-friendly software named HSkinSensDS.

Candidate Gene Identification and Transcriptome Analysis of Tomato male sterile-30 and Functional Marker Development for ms-30 and Its Alleles, ms-33, 7B-1, and stamenless-2.

Male sterility is a valuable trait for hybrid seed production in tomato (Solanum lycopersicum). The mutants male sterile-30 (ms-30) and ms-33 of tomato exhibit twisted stamens, exposed stigmas, and complete male sterility, thus holding potential for application in hybrid seed production. In this study, the ms-30 and ms-33 loci were fine-mapped to 53.3 kb and 111.2 kb intervals, respectively. Tomato PISTILLATA (TPI, syn. SlGLO2), a B-class MADS-box transcription factor gene, was identified as the most likely candidate gene for both loci. TPI is also the candidate gene of tomato male sterile mutant 7B-1 and sl-2. Allelism tests revealed that ms-30, ms-33, 7B-1, and sl-2 were allelic. Sequencing analysis showed sequence alterations in the TPI gene in all these mutants, with ms-30 exhibiting a transversion (G to T) that resulted in a missense mutation (S to I); ms-33 showing a transition (A to T) that led to alternative splicing, resulting in a loss of 46 amino acids in protein; and 7B-1 and sl-2 mutants showing the insertion of an approximately 4.8 kb retrotransposon. On the basis of these sequence alterations, a Kompetitive Allele Specific PCR marker, a sequencing marker, and an Insertion/Deletion marker were developed. Phenotypic analysis of the TPI gene-edited mutants and allelism tests indicated that the gene TPI is responsible for ms-30 and its alleles. Transcriptome analysis of ms-30 and quantitative RT-PCR revealed some differentially expressed genes associated with stamen and carpel development. These findings will aid in the marker-assisted selection for ms-30 and its alleles in tomato breeding and support the functional analysis of the TPI gene.

In silico prediction of ocular toxicity of compounds using explainable machine learning and deep learning approaches.

The accurate identification of chemicals with ocular toxicity is of paramount importance in health hazard assessment. In contemporary chemical toxicology, there is a growing emphasis on refining, reducing, and replacing animal testing in safety evaluations. Therefore, the development of robust computational tools is crucial for regulatory applications. The performance of predictive models is heavily reliant on the quality and quantity of data. In this investigation, we amalgamated the most extensive dataset (4901 compounds) sourced from governmental GHS-compliant databases and literature to develop binary classification models of chemical ocular toxicity. We employed 12 molecular representations in conjunction with six machine learning algorithms and two deep learning algorithms to create a series of binary classification models. The findings indicated that the deep learning method GCN outperformed the machine learning models in cross-validation, achieving an impressive AUC of 0.915. However, the top-performing machine learning model (RF-Descriptor) demonstrated excellent performance with an AUC of 0.869 on the test set and was therefore selected as the best model. To enhance model interpretability, we conducted the SHAP method and attention weights analysis. The two approaches offered visual depictions of the relevance of key descriptors and substructures in predicting ocular toxicity of chemicals. Thus, we successfully struck a delicate balance between data quality and model interpretability, rendering our model valuable for predicting and comprehending potential ocular-toxic compounds in the early stages of drug discovery.

In Silico Prediction of Chemical Acute Dermal Toxicity Using Explainable Machine Learning Methods.

The research on acute dermal toxicity has consistently been a crucial component in assessing the potential risks of human exposure to active ingredients in pesticides and related plant protection products. However, it is difficult to directly identify the acute dermal toxicity of potential compounds through animal experiments alone. In our study, we separately integrated 1735 experimental data based on rabbits and 1679 experimental data based on rats to construct acute dermal toxicity prediction models using machine learning and deep learning algorithms. The best models for the two animal species achieved AUC values of 78.0 and 82.0%, respectively, on 10-fold cross-validation. Additionally, we employed SARpy to extract structural alerts, and in conjunction with Shapley additive explanation and attentive FP heatmap, we identified important features and structural fragments associated with acute dermal toxicity. This approach offers valuable insights for the detection of positive compounds. Moreover, a standalone software tool was developed to make acute dermal toxicity prediction easier. In summary, our research would provide an effective tool for acute dermal toxicity evaluation of pesticides, cosmetics, and drug safety assessment.

Using a longitudinal network structure to subgroup depressive symptoms among adolescents.

BACKGROUND: Network modeling has been proposed as an effective approach to examine complex associations among antecedents, mediators and symptoms. This study aimed to investigate whether the severity of depressive symptoms affects the multivariate relationships among symptoms and mediating factors over a 2-year longitudinal follow-up. METHODS: We recruited a school-based cohort of 1480 primary and secondary school students over four semesters from January 2020 to December 2021. The participants (n = 1145) were assessed at four time points (ages 10-13 years old at baseline). Based on a cut-off score of 5 on the 9-item Patient Health Questionnaire at each time point, the participants were categorized into the non-depressive symptom (NDS) and depressive symptom (DS) groups. We conducted network analysis to investigate the symptom-to-symptom influences in these two groups over time. RESULTS: The global network metrics did not differ statistically between the NDS and DS groups at four time points. However, network connection strength varied with symptom severity. The edge weights between learning anxiety and social anxiety were prominently in the NDS group over time. The central factors for NDS and DS were oversensitivity and impulsivity (3 out of 4 time points), respectively. Moreover, both node strength and closeness were stable over time in both groups. CONCLUSIONS: Our study suggests that interrelationships among symptoms and contributing factors are generally stable in adolescents, but a higher severity of depressive symptoms may lead to increased stability in these relationships.

AttentiveSkin: To Predict Skin Corrosion/Irritation Potentials of Chemicals via Explainable Machine Learning Methods.

Skin Corrosion/Irritation (Corr./Irrit.) has long been a health hazard in the Globally Harmonized System (GHS). Several in silico models have been built to predict Skin Corr./Irrit. as an alternative to the increasingly restricted animal testing. However, current studies are limited by data amount/quality and model availability. To address these issues, we compiled a traceable consensus GHS data set comprising 731 Corr., 1283 Irrit., and 1205 negative (Neg.) samples from 6 governmental databases and 2 external data sets. Then, a series of binary classifiers were developed with five machine learning (ML) algorithms and six molecular representations. For 10-fold cross-validation, the best Corr. vs Neg. classifier achieved an Area Under the Receiver Operating Characteristic Curve (AUC) of 97.1%, while the best Irrit. vs Neg. classifier achieved an AUC of 84.7%. Compared with existing in silico tools on external validation, our Attentive FP classifiers showed the highest metrics on Corr. vs Neg. and the second highest accuracy on Irrit. vs Neg. The SHapley Additive exPlanation approach was further applied to figure out important molecular features, and the attention weights were visualized to perform interpretable prediction. Structural alerts associated with Skin Corr./Irrit. were also identified. The interpretable Attentive FP classifiers were integrated into the software AttentiveSkin at https://github.com/BeeBeeWong/AttentiveSkin. The conventional ML classifiers are also provided on our platform admetSAR at http://lmmd.ecust.edu.cn/admetsar2/. Considering the data deficiency and the limited model availability of Skin Corr./Irrit., we believe that our data set and models could facilitate chemical safety assessment and relevant studies.

Comparison of the carotid corrected flow time and tidal volume challenge for assessing fluid responsiveness in robot-assisted laparoscopic surgery.

We aimed to compare the ability of carotid corrected flow time assessed by ultrasound and the changes in dynamic preload indices induced by tidal volume challenge predicting fluid responsiveness in patients undergoing robot-assisted laparoscopic gynecological surgery in the modified head-down lithotomy position. This prospective single-center study included patients undergoing robot-assisted laparoscopic surgery in the modified head-down lithotomy position. Carotid Doppler parameters and hemodynamic data, including corrected flow time, pulse pressure variation, stroke volume variation, and stroke volume index at a tidal volume of 6 mL/kg predicted body weight and after increasing the tidal volume to 8 mL/kg predicted body weight (tidal volume challenge), respectively, were measured. Fluid responsiveness was defined as a stroke volume index ≥ 10% increase after volume expansion. Among the 52 patients included, 26 were classified as fluid responders and 26 as non-responders based on the stroke volume index. The area under the receiver operating characteristic curve measured to predict the fluid responsiveness to corrected flow time and changes in pulse pressure variation (ΔPPV6-8) after tidal volume challenge were 0.82 [95% confidence interval (CI) 0.71-0.94; P < 0.0001] and 0.85 (95% CI 0.74-0.96; P < 0.0001), respectively. The value for pulse pressure variation at a tidal volume of 8 mL/kg was 0.79 (95% CI 0.67-0.91; P = 0.0003). The optimal cut-off values for corrected flow time and ΔPPV6-8 were 357 ms and > 1%, respectively. Both the corrected flow time and Changes in pulse pressure variation after tidal volume challenge reliably predicted fluid responsiveness in patients undergoing robot-assisted laparoscopic gynecological surgery in the modified head-down lithotomy position. And pulse pressure variation at a tidal volume of 8 mL/kg maybe also a useful predictor.Trial registration: Chinese Clinical Trial Register (CHiCTR2200060573, Principal investigator: Hongliang Liu, Date of registration: 05/06/2022).

Fine Mapping of fw6.3, a Major-Effect Quantitative Trait Locus That Controls Fruit Weight in Tomato.

Tomato (Solanum lycopersicum) is a widely consumed vegetable, and the tomato fruit weight is a key yield component. Many quantitative trait loci (QTLs) controlling tomato fruit weight have been identified, and six of them have been fine-mapped and cloned. Here, four loci controlling tomato fruit weight were identified in an F2 population through QTL seq.; fruit weight 6.3 (fw6.3) was a major-effect QTL and its percentage of variation explanation (R2) was 0.118. This QTL was fine-mapped to a 62.6 kb interval on chromosome 6. According to the annotated tomato genome (version SL4.0, annotation ITAG4.0), this interval contained seven genes, including Solyc06g074350 (the SELF-PRUNING gene), which was likely the candidate gene underlying variation in fruit weight. The SELF-PRUNING gene contained a single-nucleotide polymorphism that resulted in an amino acid substitution in the protein sequence. The large-fruit allele of fw6.3 (fw6.3HG) was overdominant to the small-fruit allele fw6.3RG. The soluble solids content was also increased by fw6.3HG. These findings provide valuable information that will aid the cloning of the FW6.3 gene and ongoing efforts to breed tomato plants with higher yield and quality via molecular marker-assisted selection.

All-flesh fruit in tomato is controlled by reduced expression dosage of AFF through a structural variant mutation in the promoter.

The formation of locule gel is an important process in tomato and is a typical characteristic of berry fruit. In this study, we examined a natural tomato mutant that produces all-flesh fruit (AFF) in which the locule tissue remains in a solid state during fruit development. We constructed different genetic populations to fine-map the causal gene for this trait and identified SlMBP3 as the locus conferring the locule gel formation, which we rename as AFF. We determined the causal mutation as a 416-bp deletion in the promoter region of AFF, which reduces its expression dosage. Generally, this sequence is highly conserved among Solanaceae, as well as within the tomato germplasm. Using BC6 near-isogenic lines, we determined that the reduced expression dosage of AFF did not affect the normal development of seeds, whilst producing unique, non-liquefied locule tissue that was distinct from that of normal tomatoes in terms of metabolic components. Combined analysis using mRNA-seq and metabolomics indicated the importance of AFF in locule tissue liquefaction. Our findings provide insights into fruit-type differentiation in Solanaceae crops and also present the basis for future applications of AFF in tomato breeding programs.

OBV (obscure vein), a C2H2 zinc finger transcription factor, positively regulates chloroplast development and bundle sheath extension formation in tomato (Solanum lycopersicum) leaf veins.

Leaf veins play an important role in plant growth and development, and the bundle sheath (BS) is believed to greatly improve the photosynthetic efficiency of C4 plants. The OBV mutation in tomato (Solanum lycopersicum) results in dark veins and has been used widely in processing tomato varieties. However, physiological performance has difficulty explaining fitness in production. In this study, we confirmed that this mutation was caused by both the increased chlorophyll content and the absence of bundle sheath extension (BSE) in the veins. Using genome-wide association analysis and map-based cloning, we revealed that OBV encoded a C2H2L domain class transcription factor. It was localized in the nucleus and presented cell type-specific gene expression in the leaf veins. Furthermore, we verified the gene function by generating CRISPR/Cas9 knockout and overexpression mutants of the tomato gene. RNA sequencing analysis revealed that OBV was involved in regulating chloroplast development and photosynthesis, which greatly supported the change in chlorophyll content by mutation. Taken together, these findings demonstrated that OBV affected the growth and development of tomato by regulating chloroplast development in leaf veins. This study also provides a solid foundation to further decipher the mechanism of BSEs and to understand the evolution of photosynthesis in land plants.

Fine Mapping of the Ph-2 Gene Conferring Resistance to Late Blight (Phytophthora infestans) in Tomato.

Late blight is a devastating tomato disease. Breeding new varieties with multiple resistance (R) genes is highly effective for preventing late blight. The Ph-2 gene mediates resistance to Phytophthora infestans race T1 in tomato. In this study, we used an F2 population derived from a cross between Solanum lycopersicum Moboline (resistant) and LA3988 (susceptible) cultivars for the fine mapping of Ph-2. Two flanking markers, CAPS-1 and CC-Ase, mapped Ph-2 to a 141-kb genomic region containing 21 projected genes, 5 of which were identified as putative R genes. The Solyc10g085460 coding sequence varied significantly between the parents. The markers developed and candidate genes identified in this study shall be useful for the molecular breeding of tomato exhibiting increased late blight resistance and for the cloning of the Ph-2 gene.

SlGID1a Is a Putative Candidate Gene for qtph1.1, a Major-Effect Quantitative Trait Locus Controlling Tomato Plant Height.

Plant height is an important agronomic trait in crops. Several genes underlying tomato (Solanum lycopersicum) plant height mutants have been cloned. However, few quantitative trait genes for plant height have been identified in tomato. In this study, seven quantitative trait loci (QTLs) controlling plant height were identified in tomato. Of which, qtph1.1 (QTL for tomato plant height 1.1), qtph3.1 and qtph12.1 were major QTLs and explained 15, 16, and 12% of phenotypic variation (R2), respectively. The qtph1.1 was further mapped to an 18.9-kb interval on chromosome 1. Based on the annotated tomato genome (version SL2.50, annotation ITAG2.40), Solyc01g098390 encoding GA receptor SlGID1a was the putative candidate gene. The SlGID1a gene underlying the qtph1.1 locus contained a single nucleotide polymorphism (SNP) that resulted in an amino acid alteration in protein sequence. The near-isogenic line containing the qtph1.1 locus (NIL-qtph1.1) exhibited shorter internode length and cell length than the wild type (NIL-WT). The dwarf phenotype of NIL-qtph1.1 could not be rescued by exogenous GA3 treatment. Transcriptome analysis and real-time quantitative reverse transcription PCR (qPCR) showed that several genes related to biosynthesis and signaling of GA and auxin were differentially expressed in stems between NIL-qtph1.1 and NIL-WT. These findings might pave the road for understanding the molecular regulation mechanism of tomato plant height.

Exploration of resistance to Phelipanche aegyptiaca in tomato.

BACKGROUND: Cultivated tomatoes are highly susceptible to the destructive parasite Phelipanche aegyptiaca. Wild relatives show the potential resistance for genetic improvement. However, their genetic and molecular mechanisms are still unknown. RESULTS: Among 50 wild tomato accessions were evaluated for resistance to P. aegyptiaca, most of the wild relatives exhibited varying degrees of resistance compared to the cultivars. Solanum pennellii LA0716 performed the most promising and solid resistance with very low infection by the broomrape. The resistance involved in LA0716 was further confirmed by cytological analysis, and explored by employing a permanent introgression line (IL) population. Thirteen putative quantitative trait loci (QTLs) conferring the different resistance traits were identified. They are located on chromosomes 1, 2, 3, 4, 6, 8 and 9. The most attractive QTLs are positioned in IL6-2 and overlap with IL6-3. Specially, IL6-2 showed the highest and most consistent resistance for multiple traits and explained the major phenotypic variation of LA0716. Analysis of candidate genes involved in these regions showed that Beta (Solyc06g074240) and P450 (Solyc06g073570, Solyc06g074180 and Solyc06g074420) genes are substantially related to the strigolactone (SL) pathway. Transcript analysis further demonstrated that both Solyc06g073570 and Solyc06g074180 might play an important role in the reduction of P. aegyptiaca infection. CONCLUSION: Germplasms resistant to P. aegyptiaca were found in wild tomato species. QTLs conferring P. aegyptiaca tolerance in LA0716 were identified. IL6-2 is identified as a prospective line possessing the major QTLs. The candidate genes would provide the availability to assist the introgression of the resistance in future breeding programmes. © 2020 Society of Chemical Industry.

Anthocyanin Fruit encodes an R2R3-MYB transcription factor, SlAN2-like, activating the transcription of SlMYBATV to fine-tune anthocyanin content in tomato fruit.

Anthocyanin fruit (Aft) and atroviolacea (atv) were characterized in wild tomato and can enhance anthocyanin content in tomato fruit. However, the gene underlying the Aft locus and the mechanism by which Aft and atv act remain largely unknown. In this study, the Aft locus was fine-mapped to an approximately 145-kb interval on chromosome 10, excluding SlAN2 (Solyc10g086250), SlANT1 (Solyc10g086260) and SlANT1-like (Solyc10g086270), which have previously been suggested as candidates. Thus, the R2R3-MYB transcription factor SlAN2-like (Solyc10g086290) was considered the best candidate gene for Aft. The CRISPR/Cas9-mediated SlAN2-like mutants show a much lower accumulation of anthocyanins associated with the downregulation of multiple anthocyanin-related genes compared to the wild-type tomato, indicating that SlAN2-like is responsible for the Aft phenotype. The repressive function of SlMYBATV also was confirmed through the CRISPR/Cas9 approach. A yeast-two-hybrid assay revealed that SlMYBATV interacts with the bHLH protein SlJAF13. Furthermore, yeast-one-hybrid and dual-luciferase transient expression assays showed that Aft directly binds to the SlMYBATV promoter and activates its expression. The results herein provide candidate genes to enhance anthocyanin content in tomato fruit. This research also provides insight into a mechanism involving the Aft-SlMYBATV pathway that fine-tunes anthocyanin accumulation in tomato fruit.

A putative bHLH transcription factor is a candidate gene for male sterile 32, a locus affecting pollen and tapetum development in tomato.

The tomato (Solanum lycopersicum) male sterile 32 (ms32) mutant has been used in hybrid seed breeding programs largely because it produces no pollen and has exserted stigmas. In this study, histological examination of anthers revealed dysfunctional pollen and tapetum development in the ms32 mutant. The ms32 locus was fine mapped to a 28.5 kb interval that encoded four putative genes. Solyc01g081100, a homolog of Arabidopsis bHLH10/89/90 and rice EAT1, was proposed to be the candidate gene of MS32 because it contained a single nucleotide polymorphism (SNP) that led to the formation of a premature stop codon. A codominant derived cleaved amplified polymorphic sequence (dCAPS) marker, MS32D, was developed based on the SNP. Real-time quantitative reverse-transcription PCR showed that most of the genes, which were proposed to be involved in pollen and tapetum development in tomato, were downregulated in the ms32 mutant. These findings may aid in marker-assisted selection of ms32 in hybrid breeding programs and facilitate studies on the regulatory mechanisms of pollen and tapetum development in tomato.

Tomato locule number and fruit size controlled by natural alleles of lc and fas.

Improving yield by increasing the size of produce is an important selection criterion during the domestication of fruit and vegetable crops. Genes controlling meristem organization and organ formation work in concert to regulate the size of reproductive organs. In tomato, lc and fas control locule number, which often leads to enlarged fruits compared to the wild progenitors. LC is encoded by the tomato ortholog of WUSCHEL (WUS), whereas FAS is encoded by the tomato ortholog of CLAVATA3 (CLV3). The critical role of the WUS-CLV3 feedback loop in meristem organization has been demonstrated in several plant species. We show that mutant alleles for both loci in tomato led to an expansion of the SlWUS expression domain in young floral buds 2-3 days after initiation. Single and double mutant alleles of lc and fas maintain higher SlWUS expression during the development of the carpel primordia in the floral bud. This augmentation and altered spatial expression of SlWUS provided a mechanistic basis for the formation of multilocular and large fruits. Our results indicated that lc and fas are gain-of-function and partially loss-of-function alleles, respectively, while both mutations positively affect the size of tomato floral meristems. In addition, expression profiling showed that lc and fas affected the expression of several genes in biological processes including those involved in meristem/flower development, patterning, microtubule binding activity, and sterol biosynthesis. Several differentially expressed genes co-expressed with SlWUS have been identified, and they are enriched for functions in meristem regulation. Our results provide new insights into the transcriptional regulation of genes that modulate meristem maintenance and floral organ determinacy in tomato.

Identification and Characterization of EI (Elongated Internode) Gene in Tomato (Solanum lycopersicum).

Internode length is an important agronomic trait affecting plant architecture and crop yield. However, few genes for internode elongation have been identified in tomato. In this study, we characterized an elongated internode inbred line P502, which is a natural mutant of the tomato cultivar 05T606. The mutant P502 exhibits longer internode and higher bioactive GA concentration compared with wild-type 05T606. Genetic analysis suggested that the elongated internode trait is controlled by quantitative trait loci (QTL). Then, we identified a major QTL on chromosome 2 based on molecular markers and bulked segregant analysis (BSA). The locus was designated as EI (Elongated Internode), which explained 73.6% genetic variance. The EI was further mapped to a 75.8-kb region containing 10 genes in the reference Heinz 1706 genome. One single nucleotide polymorphism (SNP) in the coding region of solyc02g080120.1 was identified, which encodes gibberellin 2-beta-dioxygenase 7 (SlGA2ox7). SlGA2ox7, orthologous to AtGA2ox7 and AtGA2ox8, is involved in the regulation of GA degradation. Overexpression of the wild EI gene in mutant P502 caused a dwarf phenotype with a shortened internode. The difference of EI expression levels was not significant in the P502 and wild-type, but the expression levels of GA biosynthetic genes including CPS, KO, KAO, GA20ox1, GA20ox2, GA20ox4, GA3ox1, GA2ox1, GA2ox2, GA2ox4, and GA2ox5, were upregulated in mutant P502. Our results may provide a better understanding of the genetics underlying the internode elongation and valuable information to improve plant architecture of the tomato.

B-class MADS-box TM6 is a candidate gene for tomato male sterile-1526.

KEY MESSAGE: Tomato male sterile-1526 locus was fine-mapped to an interval of 44.6 kb, and a B-class MADS-box gene TM6 was identified as the candidate gene. Male sterile lines have been widely used for hybrid seed production in many crop plants. The tomato male sterile-1526 (ms-1526) mutant displays abnormal stamens and exerted stigmas and is suitable for practical use. In this study, the ms-1526 locus was fine-mapped to a 44.6 kb interval that contained four putative genes. Thereinto, Solyc02g084630 encodes tomato B-class MADS-box gene TM6 (syn. TDR6), which plays an important role in stamen development. Sequencing revealed that there was a 12.7 kb deletion in the ms-1526 region, where the promoter and first four exons of the TM6 gene were absent. ms-1547, an allele of ms-1526, also contained the same deletion in the TM6 gene. And the other allele ms-15 mutant contained a single-nucleotide polymorphism (SNP, C to A) in the coding region of the TM6 gene, which led to a missense mutation (G to W). The codominant insertion/deletion (InDel) marker MS26D and codominant derived cleaved amplified polymorphic sequence (dCAPS) marker MS15C were developed based on the deletion and SNP, respectively. A real-time quantitative reverse-transcription PCR showed that expression of the TM6 gene was barely detectable in the flowers of the ms-1526 and ms-1547 mutants. In addition, other floral organ identity genes, pollen development marker genes, and pistil marker genes were differentially expressed between wild type and mutant flowers. These findings may facilitate functional analysis of the TM6 gene and help in the marker-assisted selection of ms-15 and its alleles in tomato breeding.

Identification of Candidate HY5-Dependent and -Independent Regulators of Anthocyanin Biosynthesis in Tomato.

High quantities of anthocyanins in plants confer potential protective benefits against biotic and abiotic stressors. Studies have shown that the bZIP transcription factor HY5 plays a key role in controlling anthocyanin accumulation in response to light. However, in hy5 mutants, residual anthocyanins have been detected, indicating that other regulators exist to regulate anthocyanin biosynthesis in an HY5-independent manner. Here, we employed the CRISPR/Cas9 (clustered regularly interspersed short palindromic repeats/CRISPR-associated protein 9) system specifically to induce targeted mutagenesis of SlHY5 in the purple tomato cultivar 'Indigo Rose'. The T2 generation of tomato plants homozygous for the null allele of the SlHY5 frameshift mutated by a 1 bp insertion contained a lower anthocyanin content. Transcriptional analysis showed that most of the anthocyanin biosynthesis structural genes and several regulatory genes were down-regulated in the hy5 mutant lines. With transcriptome analyses of the various tissues from hy5 mutant lines, eight candidate transcription factors were identified that may regulate anthocyanin biosynthesis in an HY5-independent manner. These findings deepen our understanding of how light controls anthocyanin accumulation and facilitate the identification of the regulators of anthocyanin biosynthesis in an HY5-independent manner.

Fine mapping and molecular marker development of the Sm gene conferring resistance to gray leaf spot (Stemphylium spp.) in tomato.

KEY MESSAGE: The tomato gray leaf spot resistance gene Sm was fine-mapped in a 185-kb region through a map-based cloning strategy and genome-wide association study; a candidate gene was proved to be involved in Sm-mediated resistance through transient gene silencing. Gray leaf spot, caused by Stemphylium spp., is a warm weather foliar disease in tomato (Solanum lycopersicum L). Resistance against gray leaf spot is conferred by a single incompletely dominant gene (Sm) located on chromosome 11. This study aimed to map and identify molecular marker tightly linked to the Sm gene for the use of marker-assisted selection in breeding. Using an F2 population derived from a cross between the resistant line '9706' and the susceptible line 'Heinz 1706', the Sm gene was mapped to a 185-kb interval between two markers, InDel343 and InDel-FT-32 on chromosome 11, which was consistent with the result of a genome-wide association study using 289 diverse accessions. An ORF predicted in this region was proved to be involved in Sm-mediated resistance through transient gene silencing and seems to be a good candidate of the Sm locus. To clone the Sm gene, a bacterial artificial chromosome (BAC) library was screened and one BAC clone B80B15 containing the predicted ORF was identified. The analysis of sequence and structure characteristics demonstrated that the candidate gene was not a typical type resistance gene. Additionally, a co-dominant marker Sm-InDel, which produced a 122-bp or 140-bp fragment for resistant or susceptible alleles, respectively, was developed. This marker was validated in 289 germplasm and could be used in marker-assisted selection for gray leaf spot resistance.