Pathogenic and Comparative Genomic Analysis of Ralstonia pseudosolanacearum Isolated from Casuarina.

Casuarina equisetifolia is crucial in protecting coastal regions of China against typhoon attacks, but has faced a substantial challenge due to wilt disease caused by pathogens of the Ralstonia solanacearum species complex (RSSC). Although the initial outbreak of Casuarina wilt in 1970s was effectively controlled by disease-resistant C. equisetifolia varieties, the disease has recently re-emerged in coastal regions of Guangdong. In this study, we report the isolation, characterization, and comparative genomic analysis of 11 RSSC strains from diseased C. equisetifolia at various locations along the coast of Guangdong. Phylogenomic analysis showed that the strains were closely related and clustered with phylotype I strains previously isolated from peanuts. Single-gene based analysis further suggested these strains could be derived from strains present in Guangdong since the 1980s, indicating a historical context to their current pathogenicity. Casuarina-isolated strains exhibited notably higher virulence against C. equisetifolia and peanuts than representative RSSC strains GMI1000 and EP1, suggesting host-specific adaptations which possibly contributed to the recent outbreak. Comparative genomic analysis among RSSC strains revealed a largely conserved genome structure and high levels of conservation in gene clusters encoding extracellular polysaccharides biosynthesis, secretion systems, and quorum sensing regulatory systems. However, we also found a number of unique genes in the Casuarina-isolated strains that were absent in GMI1000 and EP1, and vice versa, pointing to potential genetic factors underpinning their differential virulence. These unique genes offer promising targets for future functional studies. Overall, our findings provide crucial insights into the RSSC pathogens causing Casuarina wilt in Guangdong, guiding future efforts in disease control and prevention.

Effects of copy number variations on longevity in late-onset Alzheimer's disease patients: insights from a causality network analysis.

Alzheimer's disease (AD), particularly late-onset Alzheimer's disease (LOAD), is a prevalent form of dementia that significantly affects patients' cognitive and behavioral capacities and longevity. Although approximately 70 genetic risk factors linked with AD have been identified, their influence on patient longevity remains unclear. Further, recent studies have associated copy number variations (CNVs) with the longevity of healthy individuals and immune-related pathways in AD patients. This study aims to investigate the role of CNVs on the longevity of AD patients by integrating the Whole Genome Sequencing (WGS) and transcriptomics data from the Religious Orders Study/Memory and Aging Project (ROSMAP) cohort through causality network inference. Our comprehensive analysis led to the construction of a CNV-Gene-Age of Death (AOD) causality network. We successfully identified three key CNVs (DEL5006, mCNV14192, and DUP42180) and seven AD-longevity causal genes (PLGRKT, TLR1, PLAU, CALB2, SYTL2, OTOF, and NT5DC1) impacting AD patient longevity, independent of disease severity. This outcome emphasizes the potential role of plasminogen activation and chemotaxis in longevity. We propose several hypotheses regarding the role of identified CNVs and the plasminogen system on patient longevity. However, experimental validation is required to further corroborate these findings and uncover precise mechanisms. Despite these limitations, our study offers promising insights into the genetic influence on AD patient longevity and contributes to paving the way for potential therapeutic interventions.

PhcX Is a LqsR-family response regulator that contributes to Ralstonia solanacearum virulence and regulates multiple virulence factors.

IMPORTANCE: The bacterial wilt caused by the soil-borne phytopathogen Ralstonia solanacearum is one of the most destructive crop diseases. To achieve a successful infection, R. solanacearum has evolved an intricate regulatory network to orchestrate the expression of an arsenal of virulence factors and fine-tune the allocation of energy. However, despite the wealth of knowledge gained in the past decades, many players and connections are still missing from the network. The importance of our study lies in the identification of PhcX, a novel conserved global regulator with critical roles in modulating the virulence and metabolism of R. solanacearum. PhcX affects many well-characterized regulators and exhibits contrasting modes of regulation from the central regulator PhcA on a variety of virulence-associated traits and genes. Our findings add a valuable piece to the puzzle of how the pathogen regulates its proliferation and infection, which is critical for understanding its pathogenesis and developing disease control strategies.

Pseudomonas chlororaphis L5 and Enterobacter asburiae L95 biocontrol Dickeya soft rot diseases by quenching virulence factor modulating quorum sensing signal.

Virulence factor modulating (VFM) is a quorum sensing (QS) signal shared by and specific to Dickeya bacteria, regulating the production of plant cell wall degrading enzymes (PCWDEs) and virulence of Dickeya. High polarity and trace of VFM signal increase the difficulty of signal separation and structure identification, and thus limit the development of quorum quenching strategy to biocontrol bacterial soft rot diseases caused by Dickeya. In order to high-throughput screen VFM quenching bacteria, a vfmE-gfp biosensor VR2 (VFM Reporter) sensitive to VFM signal was first constructed. Subsequently, two bacterial strains with high quenching efficiency were screened out by fluorescence intensity measurement and identified as Pseudomonas chlororaphis L5 and Enterobacter asburiae L95 using multilocus sequence analysis (MLSA). L5 and L95 supernatants reduced the expression of vfm genes, and both strains also decreased the production of PCWDEs of D. zeae MS2 and significantly reduced the virulence of D. oryzae EC1 on rice seedlings, D. zeae MS2 on banana seedlings, D. dadantii 3937 on potato and D. fangzhongdai CL3 on taro. Findings in this study provide a method to high-throughput screen VFM quenching bacteria and characterize novel functions of P. chlororaphis and E. asburiae in biocontrolling plant diseases through quenching VFM QS signal.

Pseudomonas forestsoilum sp. nov. and P. tohonis biocontrol bacterial wilt by quenching 3-hydroxypalmitic acid methyl ester.

Bacterial wilt caused by Ralstonia solanacearum ranks the second top important bacterial plant disease worldwide. It is also the most important bacterial disease threatening the healthy development of Casuarina equisetifolia protection forest. 3-hydroxypalmitic acid methyl ester (3-OH PAME) functions as an important quorum sensing (QS) signal regulating the expression of virulence genes in R. solanacearum, and has been regarded as an ideal target for disease prevention and control. To screen native microorganisms capable of degrading 3-OH PAME, samples of C. equisetifolia branches and forest soil were collected and cultured in the medium containing 3-OH PAME as the sole carbon source. Bacteria with over 85% degradation rates of 3-OH PAME after 7-day incubation were further separated and purified. As a result, strain Q1-7 isolated from forest soil and strain Q4-3 isolated from C. equisetifolia branches were obtained and identified as Pseudomonas novel species Pseudomonas forestsoilum sp. nov. and P. tohonis, respectively, according to whole genome sequencing results. The degradation efficiencies of 3-OH PAME of strains Q1-7 and Q4-3 were 95.80% and 100.00% at 48 h, respectively. Both strains showed high esterase activities and inhibited R. solanacearum exopolysaccharide (EPS) and cellulase production. Application of strains Q1-7 and Q4-3 effectively protects C. equisetifolia, peanut and tomato plants from infection by R. solanacearum. Findings in this study provide potential resources for the prevention and control of bacterial wilt caused by R. solanacearum, as well as valuable materials for the identification of downstream quenching genes and the research and development of quenching enzymes for disease control.

A novel genome sequence of Jasminum sambac helps uncover the molecular mechanism underlying the accumulation of jasmonates.

Jasminum sambac is a well-known plant for its attractive and exceptional fragrance, the flowers of which are used to produce scented tea. Jasmonate (JA), an important plant hormone was first identified in Jasminum species. Jasmine plants contain abundant JA naturally, of which the molecular mechanisms of synthesis and accumulation are not clearly understood. Here, we report a telomere-to-telomere consensus assembly of a double-petal J. sambac genome along with two haplotype-resolved genomes. We found that gain-and-loss, positive selection, and allelic specific expression of aromatic volatile-related genes contributed to the stronger flower fragrance in double-petal J. sambac compared with single- and multi-petal jasmines. Through comprehensive comparative genomic, transcriptomic, and metabolomic analyses of double-petal J. sambac, we revealed the genetic basis of the production of aromatic volatiles and salicylic acid (SA), and the accumulation of JA under non-stress conditions. We identified several key genes associated with JA biosynthesis, and their non-stress related activities lead to extraordinarily high concentrations of JA in tissues. High JA synthesis coupled with low degradation in J. sambac results in accumulation of high JA under typical environmental conditions, similar to the accumulation mechanism of SA. This study offers important insights into the biology of J. sambac, and provides valuable genomic resources for further utilization of natural products.

Genome-Wide Analysis and Characterization of Eggplant F-Box Gene Superfamily: Gene Evolution and Expression Analysis under Stress.

F-box genes play an important role in plant growth and resistance to abiotic and biotic stresses. To date, systematic analysis of F-box genes and functional annotation in eggplant (Solanum melongena) is still limited. Here, we identified 389 F-box candidate genes in eggplant. The domain study of F-box candidate genes showed that the F-box domain is conserved, whereas the C-terminal domain is diverse. There are 376 SmFBX candidate genes distributed on 12 chromosomes. A collinearity analysis within the eggplant genome suggested that tandem duplication is the dominant form of F-box gene replication in eggplant. The collinearity analysis between eggplant and the three other species (Arabidopsis thaliana, rice and tomato) provides insight into the evolutionary characteristics of F-box candidate genes. In addition, we analyzed the expression of SmFBX candidate genes in different tissues under high temperature and bacterial wilt stress. The results identified several F-box candidate genes that potentially participate in eggplant heat tolerance and bacterial wilt resistance. Moreover, the yeast two-hybrid assay showed that several representative F-box candidate proteins interacted with representative Skp1 proteins. Overexpression of SmFBX131 and SmFBX230 in tobacco increased resistance to bacterial wilt. Overall, these results provide critical insights into the functional analysis of the F-box gene superfamily in eggplant and provide potentially valuable targets for heat and bacterial resistance.

Comparative Pathogenomic Analysis of Two Banana Pathogenic Dickeya Strains Isolated from China and the Philippines.

Dickeya is a major and typical member of soft rot Pectobacteriaceae (SRP) with a wide range of plant hosts worldwide. Previous studies have identified D. zeae as the causal agent of banana soft rot disease in China. In 2017, we obtained banana soft rot pathogen strain FZ06 from the Philippines. Genome sequencing and analysis indicated that FZ06 can be classified as D. dadantii and represents a novel subspecies of D. dadantii, which we propose to name as subsp. paradisiaca. Compared with Chinese banana soft rot pathogenic strain D. zeae MS2, strain FZ06 has a similar host range but different virulence; FZ06 is significantly less virulent to banana and potato but more virulent to Chinese cabbage and onion. Characterization of virulence factors revealed obviously less production of pectate lyases (Pels), polygalacturonases (Pehs), proteases (Prts), and extrapolysaccharides (EPSs), as well as lower swimming and swarming motility and biofilm formation in strain FZ06. Genomic comparison of the two strains revealed five extra gene clusters in FZ06, including one Stt-type T2SS, three T4SSs, and one T4P. Expression of cell wall degrading enzyme (CWDE)-encoding genes is significantly lower in FZ06 than in MS2.

Machine Learning Models for Data-Driven Prediction of Diabetes by Lifestyle Type.

The prevalence of diabetes has been increasing in recent years, and previous research has found that machine-learning models are good diabetes prediction tools. The purpose of this study was to compare the efficacy of five different machine-learning models for diabetes prediction using lifestyle data from the National Health and Nutrition Examination Survey (NHANES) database. The 1999-2020 NHANES database yielded data on 17,833 individuals data based on demographic characteristics and lifestyle-related variables. To screen training data for machine models, the Akaike Information Criterion (AIC) forward propagation algorithm was utilized. For predicting diabetes, five machine-learning models (CATBoost, XGBoost, Random Forest (RF), Logistic Regression (LR), and Support Vector Machine (SVM)) were developed. Model performance was evaluated using accuracy, sensitivity, specificity, precision, F1 score, and receiver operating characteristic (ROC) curve. Among the five machine-learning models, the dietary intake levels of energy, carbohydrate, and fat, contributed the most to the prediction of diabetes patients. In terms of model performance, CATBoost ranks higher than RF, LG, XGBoost, and SVM. The best-performing machine-learning model among the five is CATBoost, which achieves an accuracy of 82.1% and an AUC of 0.83. Machine-learning models based on NHANES data can assist medical institutions in identifying diabetes patients.

The integration host factor regulates multiple virulence pathways in bacterial pathogen Dickeya zeae MS2.

Dickeya zeae is an aggressive bacterial phytopathogen that infects a wide range of host plants. It has been reported that integration host factor (IHF), a nucleoid-associated protein consisting of IHFα and IHFß subunits, regulates gene expression by influencing nucleoid structure and DNA bending. To define the role of IHF in the pathogenesis of D. zeae MS2, we deleted either and both of the IHF subunit encoding genes ihfA and ihfB, which significantly reduced the production of cell wall-degrading enzymes (CWDEs), an unknown novel phytotoxin and the virulence factor-modulating (VFM) quorum-sensing (QS) signal, cell motility, biofilm formation, and thereafter the infection ability towards both potato slices and banana seedlings. To characterize the regulatory pathways of IHF protein associated with virulence, IHF binding sites (consensus sequence 5'-WATCAANNNNTTR-3') were predicted and 272 binding sites were found throughout the genome. The expression of 110 tested genes was affected by IHF. Electrophoretic mobility shift assay (EMSA) showed direct interaction of IhfA protein with the promoters of vfmE, speA, pipR, fis, slyA, prtD, hrpL, hecB, hcp, indA, hdaA, flhD, pilT, gcpJ, arcA, arcB, and lysR. This study clarified the contribution of IHF in the pathogenic process of D. zeae by controlling the production of VFM and putrescine QS signals, phytotoxin, and indigoidine, the luxR-solo system, Fis, SlyA, and FlhD transcriptional regulators, and secretion systems from type I to type VI. Characterization of the regulatory networks of IHF in D. zeae provides a target for prevention and control of plant soft rot disease.

Isolation and Genome Analysis of Pectobacterium colocasium sp. nov. and Pectobacterium aroidearum, Two New Pathogens of Taro.

Bacterial soft rot is one of the most destructive diseases of taro (Colocasia esculenta) worldwide. In recent years, frequent outbreaks of soft rot disease have seriously affected taro production and became a major constraint to the development of taro planting in China. However, little is known about the causal agents of this disease, and the only reported pathogens are two Dickeya species and P. carotovorum. In this study, we report taro soft rot caused by two novel Pectobacterium strains, LJ1 and LJ2, isolated from taro corms in Ruyuan County, Shaoguan City, Guangdong Province, China. We showed that LJ1 and LJ2 fulfill Koch's postulates for taro soft rot. The two pathogens can infect taro both individually and simultaneously, and neither synergistic nor antagonistic interaction was observed between the two pathogens. Genome sequencing of the two strains indicated that LJ1 represents a novel species of the genus Pectobacterium, for which the name "Pectobacterium colocasium sp. nov." is proposed, while LJ2 belongs to Pectobacterium aroidearum. Pan-genome analysis revealed multiple pathogenicity-related differences between LJ1, LJ2, and other Pectobacterium species, including unique virulence factors, variation in the copy number and organization of Type III, IV, and VI secretion systems, and differential production of plant cell wall degrading enzymes. This study identifies two new soft rot Pectobacteriaceae (SRP) pathogens causing taro soft rot in China, reports a new case of co-infection of plant pathogens, and provides valuable resources for further investigation of the pathogenic mechanisms of SRP.

Microevolution of the mexT and lasR Reinforces the Bias of Quorum Sensing System in Laboratory Strains of Pseudomonas aeruginosa PAO1.

The Pseudomonas aeruginosa strain PAO1 has routinely been used as a laboratory model for quorum sensing (QS). However, the microevolution of P. aeruginosa laboratory strains resulting in genetic and phenotypic variations have caused inconsistencies in QS research. To investigate the underlying causes of these variations, we analyzed 5 Pseudomonas aeruginosa PAO1 sublines from our laboratory using a combination of phenotypic characterization, high throughput genome sequencing, and bioinformatic analysis. The major phenotypic variations among the sublines spanned across the levels of QS signals and virulence factors such as pyocyanin and elastase. Furthermore, the sublines exhibited distinct variations in motility and biofilm formation. Most of the phenotypic variations were mapped to mutations in the lasR and mexT, which are key components of the QS circuit. By introducing these mutations in the subline PAO1-E, which is devoid of such mutations, we confirmed their influence on QS, virulence, motility, and biofilm formation. The findings further highlight a possible divergent regulatory mechanism between the LasR and MexT in the P. aeruginosa. The results of our study reveal the effects of microevolution on the reproducibility of most research data from QS studies and further highlight mexT as a key component of the QS circuit of P. aeruginosa.

Genomic and Functional Dissections of Dickeya zeae Shed Light on the Role of Type III Secretion System and Cell Wall-Degrading Enzymes to Host Range and Virulence.

Dickeya zeae is a worldwide destructive pathogen that causes soft rot diseases on various hosts such as rice, maize, banana, and potato. The strain JZL7 we recently isolated from clivia represents the first monocot-specific D. zeae and also has reduced pathogenicity compared to that of other D. zeae strains (e.g., EC1 and MS2). To elucidate the molecular mechanisms underlying its more restricted host range and weakened pathogenicity, we sequenced the complete genome of JZL7 and performed comparative genomic and functional analyses of JZL7 and other D. zeae strains. We found that, while having the largest genome among D. zeae strains, JZL7 lost almost the entire type III secretion system (T3SS), which is a key component of the virulence suite of many bacterial pathogens. Importantly, the deletion of T3SS in MS2 substantially diminished the expression of most type III secreted effectors (T3SEs) and MS2's pathogenicity on both dicots and monocots. Moreover, although JZL7 and MS2 share almost the same repertoire of cell wall-degrading enzymes (CWDEs), we found broad reduction in the production of CWDEs and expression levels of CWDE genes in JZL7. The lower expression of CWDEs, pectin lyases in particular, would probably make it difficult for JZL7 to break down the cell wall of dicots, which is rich in pectin. Together, our results suggest that the loss of T3SS and reduced CWDE activity together might have contributed to the host specificity and virulence of JZL7. Our findings also shed light on the pathogenic mechanism of Dickeya and other soft rot Pectobacteriaceae species in general. IMPORTANCE Dickeya zeae is an important, aggressive bacterial phytopathogen that can cause severe diseases in many crops and ornamental plants, thus leading to substantial economic losses. Strains from different sources showed significant diversity in their natural hosts, suggesting complicated evolution history and pathogenic mechanisms. However, molecular mechanisms that cause the differences in the host range of D. zeae strains remain poorly understood. This study carried out genomic and functional dissections of JZL7, a D. zeae strain with restricted host range, and revealed type III secretion system (T3SS) and cell wall-degrading enzymes (CWDEs) as two major factors contributing to the host range and virulence of D. zeae, which will provide a valuable reference for the exploration of pathogenic mechanisms in other bacteria and present new insights for the control of bacterial soft rot diseases on crops.

Microbial Diversity Analysis and Genome Sequencing Identify Xanthomonas perforans as the Pathogen of Bacterial Leaf Canker of Water Spinach (Ipomoea aquatic).

Ipomoea aquatica is a leafy vegetable widely cultivated in tropical Asia, Africa, and Oceania. Bacterial leaf canker disease has been attacking the planting fields and seriously affecting the quality of I. aquatica in epidemic areas in China. This study examined the microbial composition of I. aquatica leaves with classical symptoms of spot disease. The results showed that Xanthomonas was overwhelmingly dominant in all four diseased leaf samples but rarely present in rhizospheric soil or irrigation water samples. In addition, Pantoea was also detected in two of the diseased leaf samples. Pathogen isolation, identification, and inoculation revealed that both Xanthomonas sp. TC2-1 and P. ananatis were pathogenic to the leaves of I. aquatic, causing crater-shaped ulcerative spots and yellowing with big brown rot lesions on leaves, respectively. We further sequenced the whole genome of strain TC2-1 and showed that it is a member of X. perforans. Overall, this study identified X. perforans as the causal pathogen of I. aquatica bacterial leaf canker, and P. ananatis as a companion pathogen causing yellowing and brown rot on leaves. The correct identification of the pathogens will provide important basis for future efforts to formulate targeted application strategy for bacterial disease control.

Isolation, Characterization, and Genomic Investigation of a Phytopathogenic Strain of Stenotrophomonas maltophilia.

Stenotrophomonas maltophilia is ubiquitous in diverse environmental habitats. It merits significant concern because of its increasing incidence of nosocomial and community-acquired infection in immunocompromised patients and multiple drug resistance. It is rarely reported as a phytopathogen except in causing white stripe disease of rice in India and postharvest fruit rot of Lanzhou lily. For this study, Dickeya zeae and S. maltophilia strains were simultaneously isolated from soft rot leaves of Clivia miniata in Guangzhou, China, and were both demonstrated to be pathogenic to the host. Compared with the D. zeae strains, S. maltophilia strains propagated faster for greater growth in lysogeny broth medium and produced no cellulases or polygalacturonases, but did produce more proteases and fewer extracellular polysaccharides. Furthermore, S. maltophilia strains swam and swarmed dramatically less on semisolid media, but formed a great many more biofilms. Both D. zeae and S. maltophilia strains isolated from clivia caused rot symptoms on other monocot hosts, but not on dicots. Similar to previously reported S. maltophilia strains isolated from other sources, the strain JZL8 survived under many antibiotic stresses. The complete genome sequence of S. maltophilia strain JZL8 consists of a chromosome of 4,635,432 bp without a plasmid. Pan-genome analysis of JZL8 and 180 other S. maltophilia strains identified 50 genes that are unique to JZL8, seven of which implicate JZL8 as the potential pathogen contributor in plants. JZL8 also contains three copies of Type I Secretion System machinery; this is likely responsible for its greater production of proteases. Findings from this study extend our knowledge on the host range of S. maltophilia and provide insight into the phenotypic and genetic features underlying the plant pathogenicity of JZL8.

Identification of an Embryonic Cell-Specific Region within the Pineapple SERK1 Promoter.

Plant tissue culture methods, such as somatic embryogenesis, are attractive alternatives to traditional breeding methods for plant propagation. However, they often suffer from limited efficiency. Somatic embryogenesis receptor kinase (SERK)1 is a marker gene of early somatic embryogenesis in several plants, including pineapple. It can be selectively induced and promotes a key step in somatic embryogenesis. We investigated the embryonic cell-specific transcriptional regulation of AcSERK1 by constructing a series of vectors carrying the GUS（Beta-glucuronidase） reporter gene under the control of different candidate cis-regulatory sequences. These vectors were transfected into both embryonic and non-embryonic callus, and three immature embryo stages and the embryonic-specific activity of the promoter fragments was analyzed. We found that the activity of the regulatory sequence of AcSERK1 lacking -983 nt ~-880 nt, which included the transcription initiation site, was significantly reduced in the embryonic callus of pineapple, accompanied by the loss of embryonic cell-specific promoter activity. Thus, this fragment is an essential functional segment with highly specific promoter activity for embryonic cells, and it is active only from the early stages of somatic embryo development to the globular embryo stage. This study lays the foundation for identifying mechanisms that enhance the efficiency of somatic embryogenesis in pineapple and other plants.

Genome-wide investigation of WRKY gene family in pineapple: evolution and expression profiles during development and stress.

BACKGROUND: WRKY proteins comprise a large family of transcription factors that play important roles in many aspects of physiological processes and adaption to environment. However, little information was available about the WRKY genes in pineapple (Ananas comosus), an important tropical fruits. The recent release of the whole-genome sequence of pineapple allowed us to perform a genome-wide investigation into the organization and expression profiling of pineapple WRKY genes. RESULTS: In the present study, 54 pineapple WRKY (AcWRKY) genes were identified and renamed on the basis of their respective chromosome distribution. According to their structural and phylogenetic features, the 54 AcWRKYs were further classified into three main groups with several subgroups. The segmental duplication events played a major role in the expansion of pineapple WRKY gene family. Synteny analysis and phylogenetic comparison of group III WRKY genes provided deep insight into the evolutionary characteristics of pineapple WRKY genes. Expression profiles derived from transcriptome data and real-time quantitative PCR analysis exhibited distinct expression patterns of AcWRKY genes in various tissues and in response to different abiotic stress and hormonal treatments. CONCLUSIONS: Fifty four WRKY genes were identified in pineapple and the structure of their encoded proteins, their evolutionary characteristics and expression patterns were examined in this study. This systematic analysis provided a foundation for further functional characterization of WRKY genes with an aim of pineapple crop improvement.

Genome-wide organization and expression profiling of the R2R3-MYB transcription factor family in pineapple (Ananas comosus).

BACKGROUND: The MYB proteins comprise one of the largest families of plant transcription factors, which are involved in various plant physiological and biochemical processes. Pineapple (Ananas comosus) is one of three most important tropical fruits worldwide. The completion of pineapple genome sequencing provides a great opportunity to investigate the organization and evolutionary traits of pineapple MYB genes at the genome-wide level. RESULTS: In the present study, a total of 94 pineapple R2R3-MYB genes were identified and further phylogenetically classified into 26 subfamilies, as supported by the conserved gene structures and motif composition. Collinearity analysis indicated that the segmental duplication events played a crucial role in the expansion of pineapple MYB gene family. Further comparative phylogenetic analysis suggested that there have been functional divergences of MYB gene family during plant evolution. RNA-seq data from different tissues and developmental stages revealed distinct temporal and spatial expression profiles of the AcMYB genes. Further quantitative expression analysis showed the specific expression patterns of the selected putative stress-related AcMYB genes in response to distinct abiotic stress and hormonal treatments. The comprehensive expression analysis of the pineapple MYB genes, especially the tissue-preferential and stress-responsive genes, could provide valuable clues for further function characterization. CONCLUSIONS: In this work, we systematically identified AcMYB genes by analyzing the pineapple genome sequence using a set of bioinformatics approaches. Our findings provide a global insight into the organization, phylogeny and expression patterns of the pineapple R2R3-MYB genes, and hence contribute to the greater understanding of their biological roles in pineapple.