Comparative genomics and transcriptomics analysis of the bHLH gene family indicate their roles in regulating flavonoid biosynthesis in Sophora flavescens.

The basic helix-loop-helix (bHLH) transcription factors play crucial roles in various processes, such as plant development, secondary metabolism, and response to biotic/abiotic stresses. Sophora flavescens is a widely used traditional herbal medicine in clinical practice, known for its abundant flavonoids as the main active compounds. However, there has been no comprehensive analysis of S. flavescens bHLH (SfbHLH) gene family reported currently. In this study, we identified 167 SfbHLH genes and classified them into 23 subfamilies based on comparative genomics and phylogenetic analysis. Furthermore, widespread duplications significantly contributed to the expansion of SfbHLH family. Notably, SfbHLH042 was found to occupy a central position in the bHLH protein-protein interaction network. Transcriptome analysis of four tissues (leaf, stem, root and flower) revealed that most SfbHLH genes exhibited high expression levels exclusively in specific tissues of S. flavescens. The integrated analysis of transcriptomics and metabolomics during pod development stages revealed that SfbHLH042 may play a central role in connecting SfbHLH genes, flavonoids, and key enzymes involved in the biosynthesis pathway. Moreover, we also checked the expression of 8 SfbHLH genes using RT-qPCR analysis to realize the expression profiles of these genes among various tissues at different cultivated periods and root development. Our study would aid to understand the phylogeny and expression profile of SfbHLH family genes, and provide a promising candidate gene, SfbHLH042, for regulating the biosynthesis of flavonoids in S. flavescens.

Chromosome-scale genome assembly of Astragalus membranaceus using PacBio and Hi-C technologies.

Astragalus membranaceus (Fisch.) Bge (AM) is a medicinal herb plant belonging to the Leguminosae family. In this study, we present a chromosome-scale genome assembly of AM, aiming to enhance the molecular biology and functional studies of Astragali Radix. The genome size of AM is about 1.43 Gb, with a contig N50 value of 1.67 Mb. A total of 98.16% of the assembly anchored to 9 pseudochromosomes using Hi-C technology. The assembly completeness was estimated to be 97.27% using BUSCO with the long terminal repeat assembly index (LAI) of 16.22 and quality value (QV) of 48.58. Additionally, the genome contained 67.98% repetitive sequences. Genome annotation predicted 29,914 protein-coding genes, including 73 genes involved in the flavonoid biosynthetic pathway and 2,048 transcription factors. The high-quality genome assembly and gene annotation resources will greatly facilitate future functional genomic studies in Leguminosae species.

Physiological and transcriptomic profiles reveal key regulatory pathways involved in cold resistance in sunflower seedlings.

During sunflower growth, cold waves often occur and impede plant growth. Therefore, it is crucial to study the underlying mechanism of cold resistance in sunflowers. In this study, physiological analysis revealed that as cold stress increased, the levels of ROS, malondialdehyde, ascorbic acid, and dehydroascorbic acid and the activities of antioxidant enzymes increased. Transcriptomics further identified 10,903 DEGs between any two treatments. Clustering analysis demonstrated that the expression of MYB44a, MYB44b, MYB12, bZIP2 and bZIP4 continuously upregulated under cold stress. Cold stress can induce ROS accumulation, which interacts with hormone signals to activate cold-responsive transcription factors regulating target genes involved in antioxidant defense, secondary metabolite biosynthesis, starch and sucrose metabolism enhancement for improved cold resistance in sunflowers. Additionally, the response of sunflowers to cold stress may be independent of the CBF pathway. These findings enhance our understanding of cold stress resistance in sunflowers and provide a foundation for genetic breeding.

The alleviating effect on the growth, chlorophyll synthesis, and biochemical defense system in sunflowers under cadmium stress achieved through foliar application of humic acid.

BACKGROUND: With the progress of industrialization and urbanization, cadmium (Cd) pollution in farmland is increasingly severe, greatly affecting human health. Sunflowers possess high resistance to Cd stress and great potential for phytoremediation of Cd-contaminated soil. Previous studies have shown that humic acid (HA) effectively mitigates plant damage induced by Cd; however, its alleviating effects on sunflower plants under Cd stress remain largely unknown. RESULTS: We employed four different concentrations of HA (50, 100, 200, and 300 mg L-1) via foliar application to examine their ability to alleviate Cd stress on sunflower plants' growth, chlorophyll synthesis, and biochemical defense system. The results revealed that Cd stress not only reduced plant height, stem diameter, fresh and dry weight, and chlorophyll content in sunflower plants but also altered their chlorophyll fluorescence characteristics compared to the control group. After Cd stress, the photosynthetic structure was damaged and the number of PSII reactive centers per unit changed. Application of 200 mg L-1 HA promotes sunflower growth and increases chlorophyll content. HA significantly enhances antioxidant enzyme activities (SOD, POD, CAT, and APX) and reduces ROS content (O2 -, H2O2 and -OH). Totally, Application of 200 mg L-1 HA had the best effect than other concentrations to alleviate the Cd-induced stress in sunflower plants. CONCLUSIONS: The foliar application of certain HA concentration exhibited the most effective alleviation of Cd-induced stress on sunflower plants. It can enhance the light energy utilization and antioxidant enzyme activities, while reduce ROS contents in sunflower plants. These findings provide a theoretical basis for using HA to mitigate Cd stress in sunflowers.

Effects of micro/nano-ozone bubble nutrient solutions on growth promotion and rhizosphere microbial community diversity in soilless cultivated lettuces.

Due to its high efficacy as a wide-spectrum disinfectant and its potential for the degradation of pollutants and pesticides, ozone has broad application prospects in agricultural production. In this study, micro/nano bubble technology was applied to achieve a saturation state of bubble nutrient solution, including micro-nano oxygen (O2 group) and micro-nano ozone (O3 group) bubble nutrient solutions. The effects of these solutions on lettuce physiological indices as well as changes in the microbial community within the rhizosphere substrate were studied. The application of micro/nano (O2 and O3) bubble nutrient solutions to substrate-cultured lettuce plants increased the amount of dissolved oxygen in the nutrient solution, increased the lettuce yield, and elevated the net photosynthetic rate, conductance of H2O and intercellular carbon dioxide concentration of lettuce plants. Diversity analysis of the rhizosphere microbial community revealed that both the abundance and diversity of bacterial and fungal communities in the substrate increased after plant cultivation and decreased following treatment with micro/nanobubble nutrient solutions. RDA results showed that the microbial community in the S group was positively associated with EC, that in the CK and O2 groups exhibited a positive correlation with SC, and that in the O3 group displayed a positive correlation with CAT and POD. Overall, the implementation of micro/nanobubble generation technology in soilless substrates can effectively increase the lettuce growth and yield, and O3 had a more pronounced effect on lettuce yield and quality and the microbial community structure in the substrate than O2. Our study would provide a reference and theoretical basis for developing sustainable and green technology for promoting lettuce production and can be a promising alternative to conventional methods for improving crop yields.

Genome-wide identification and expression analysis of the Dof gene family reveals their involvement in hormone response and abiotic stresses in sunflower (Helianthus annuus L.).

DNA binding with one finger (Dof), plant-specific zinc finger transcription factors, can participate in various physiological and biochemical processes during the life of plants. As one of the most important oil crops in the world, sunflower (Helianthus annuus L.) has significant economic and ornamental value. However, a systematic analysis of H. annuus Dof (HaDof) members and their functions has not been extensively conducted. In this study, we identified 50 HaDof genes that are unevenly distributed on 17 chromosomes of sunflower. We present a comprehensive overview of the HaDof genes, including their chromosome locations, phylogenetic analysis, and expression profile characterization. Phylogenetic analysis classified the 366 Dof members identified from 11 species into four groups (further subdivided into nine subfamilies). Segmental duplications are predominantly contributed to the expansion of sunflower Dof genes, and all segmental duplicate gene pairs are under purifying selection due to strong evolutionary constraints. Furthermore, we observed differential expression patterns for HaDof genes in normal tissues as well as under hormone treatment or abiotic stress conditions by analyzing RNA-seq data from previous studies and RT-qPCR data in our current study. The expression of HaDof04 and HaDof43 were not detected in any samples, which implied that they may be gradually undergoing pseudogenization process. Some HaDof genes, such as HaDof25 and HaDof30, showed responsiveness to exogenous plant hormones, such as kinetin, brassinosteroid, auxin or strigolactone, while others like HaDof15 and HaDof35 may participate in abiotic stress resistance of sunflower seedling. Our study represents the initial step towards understanding the phylogeny and expression characterization of sunflower Dof family genes, which may provide valuable reference information for functional studies on hormone response, abiotic stress resistance, and molecular breeding in sunflower and other species.

Physiology and transcriptomics highlight the underlying mechanism of sunflower responses to drought stress and rehydration.

Drought can adversely influence the crop growth and production. Accordingly, sunflowers have strong adaptability to drought; hence, we conducted analyses for sunflower seedlings with drought stress and rehydration drought acclimation through physiological measurements and transcriptomics. It showed that drought can cause the accumulation of ROS and enhance the activity of antioxidant enzymes and the content of osmolytes. After rehydration, the contents of ROS and MDA were significantly reduced concomitant with increased antioxidant activity and osmotic adjustment. Totally, 2,589 DEGs were identified among treatments. Functional enrichment analysis showed that DEGs were mainly involved in plant hormone signal transduction, MAPK signaling, and biosynthesis of secondary metabolites. Comparison between differentially spliced genes and DEGs indicated that bHLH025, NAC53, and SINAT3 may be pivotal genes involved in sunflower drought resistance. Our results not only highlight the underlying mechanism of drought stress and rehydration in sunflower but also provide a theoretical basis for crop genetic breeding.

RBBP4 is an epigenetic barrier for the induced transition of pluripotent stem cells into totipotent 2C-like cells.

Cellular totipotency is critical for whole-organism generation, yet how totipotency is established remains poorly illustrated. Abundant transposable elements (TEs) are activated in totipotent cells, which is critical for embryonic totipotency. Here, we show that the histone chaperone RBBP4, but not its homolog RBBP7, is indispensable for maintaining the identity of mouse embryonic stem cells (mESCs). Auxin-induced degradation of RBBP4, but not RBBP7, reprograms mESCs to the totipotent 2C-like cells. Also, loss of RBBP4 enhances transition from mESCs to trophoblast cells. Mechanistically, RBBP4 binds to the endogenous retroviruses (ERVs) and functions as an upstream regulator by recruiting G9a to deposit H3K9me2 on ERVL elements, and recruiting KAP1 to deposit H3K9me3 on ERV1/ERVK elements, respectively. Moreover, RBBP4 facilitates the maintenance of nucleosome occupancy at the ERVK and ERVL sites within heterochromatin regions through the chromatin remodeler CHD4. RBBP4 depletion leads to the loss of the heterochromatin marks and activation of TEs and 2C genes. Together, our findings illustrate that RBBP4 is required for heterochromatin assembly and is a critical barrier for inducing cell fate transition from pluripotency to totipotency.

Draft Genome of White-blotched River Stingray Provides Novel Clues for Niche Adaptation and Skeleton Formation.

The white-blotched river stingray (Potamotrygon leopoldi) is a cartilaginous fish native to the Xingu River, a tributary of the Amazon River system. As a rare freshwater-dwelling cartilaginous fish in the Potamotrygonidae family in which no member has the genome sequencing information, P. leopoldi provides the evolutionary details in fish phylogeny, niche adaptation, and skeleton formation. In this study, we present its draft genome of 4.11 Gb comprised of 16,227 contigs and 13,238 scaffolds, with contig N50 of 3937 kb and scaffold N50 of 5675 kb in size. Our analysis shows that P. leopoldi is a slow-evolving fish that diverged from elephant sharks about 96 million years ago. Moreover, two gene families related to the immune system, immunoglobulin heavy constant delta genes and T-cell receptor alpha/delta variable genes, exhibit expantion in P. leopoldi only. We also identified the Hox gene clusters in P. leopoldi and discovered that seven Hox genes shared by five representative fish species are missing in P. leopoldi. The RNA sequencing data from P. leopoldi and other three fish species demonstrate that fishes have a more diversified tissue expression spectrum as compared with the corresponding mammalian data. Our functional studies suggest that the lack of the GC gene encoding vitamin D-binding protein in cartilaginous fishes (both P. leopoldi and Callorhinchus milii) could partly explain the absence of hard bone in their endoskeleton. Overall, this genome resource provides new insights into the niche adaptation, body plan, and skeleton formation of P. leopoldi as well as the genome evolution in cartilaginous fishes.

Phylogenetic analysis and expression profiles of jasmonate ZIM-domain gene family provide insight into abiotic stress resistance in sunflower.

Jasmonate ZIM-domain (JAZ) proteins act as inhibitory factors of the jasmonic acid (JA) pathway, which is involved in regulating plant development and defense responses. However, there are no extensive studies available on JAZ genes in sunflower (Helianthus annuus L.). In this study, the phylogenetic analysis of 139 putative JAZ genes from eight plants demonstrated that these JAZs could be divided into five groups (Groups I-V), and the 27 sunflower JAZs (HaJAZs) were classified into these five groups. All groups contained genes from both monocotyledons and dicotyledons, indicating that the emergence of JAZ genes predates the differentiation of monocotyledons and dicotyledons. Both segmental and tandem duplications contributed greatly to this gene family's expansion in sunflower, especially in Group II. Moreover, the expression profiles of HaJAZ genes under normal conditions, hormone treatments or abiotic stresses were analyzed based on RNA-seq data. HaJAZ2 may be undergoing pseudogenization as a nonfunctional gene because it was not expressed in any tissue. Many HaJAZ genes in roots upregulated their expression when involved in responding to exogenous hormones, especially methyl-jasmonate. The abiotic stress treatments of sunflower showed that HaJAZ5, HaJAZ15, HaJAZ17, HaJAZ20, and HaJAZ21 tend to be sensitive to certain abiotic stresses. HaJAZs from different groups may share similar functions but also exercise their unique functions when responding to abiotic stresses. We speculated that this gene family was conserved in sequence but varied in its expression among duplicated HaJAZ genes, which implies that they may confer neofunctionalization in the adaptation to abiotic stresses; this work provides insight into the resistance of sunflowers and their adaptation to diverse environmental conditions.

Insights into the plateau adaptation of Salvia castanea by comparative genomic and WGCNA analyses.

INTRODUCTION: Salvia castanea, a wild plant species is adapted to extreme Qinghai-Tibetan plateau (QTP) environments. It is also used for medicinal purposes due to high ingredient of tanshinone IIA (T-IIA). Despite its importance to Chinese medicinal industry, the mechanisms associated with secondary metabolites accumulation (i.e. T-IIA and rosmarinic acid (RA)) in this species have not been characterized. Also, the role of special underground tissues in QTP adaptation of S. castanea is still unknown. OBJECTIVES: We explored the phenomenon of periderm-like structure in underground stem center of S. castanea with an aim to unravel the molecular evolutionary mechanisms of QTP adaptation in this species. METHODS: Morphologic observation and full-length transcriptome of S. castanea plants were conducted. Comparative genomic analyses of S. castanea with other 14 representative species were used to reveal its phylogenetic position and molecular evolutionary mechanisms. RNA-seq and WGCNA analyses were applied to understand the mechanisms of high accumulations of T-IIA and RA in S. castanea tissues. RESULTS: Based on anatomical observations, we proposed a "trunk-branches" developmental model to explain periderm-like structure in the center of underground stem of S. castanea. Our study suggested that S. castanea branched off from cultivated Danshen around 16 million years ago. During the evolutionary process, significantly expanded orthologous gene groups, 24 species-specific and 64 positively selected genes contributed to morphogenesis and QTP adaptation in S. castanea. RNA-seq and WGCNA analyses unraveled underlying mechanisms of high accumulations of T-IIA and RA in S. castanea and identified NAC29 and TGA22 as key transcription factors. CONCLUSION: We proposed a "trunk-branches" developmental model for the underground stem in S. castanea. Adaptations to extreme QTP environment in S. castanea are associated with accumulations of high secondary metabolites in this species.

Comparative metabolomics study of Tartary (Fagopyrum tataricum (L.) Gaertn) and common (Fagopyrum esculentum Moench) buckwheat seeds.

Tartary buckwheat has higher health-promoting value than common buckwheat. However, the related metabolites information except flavonoids is largely deficient. Here, we compared the seed metabolomes of the two species using a UHPLC-QqQ-MS-based metabolomics approach. In total, 722 metabolites were obtained, of which 84 and 78 were identified as the key active ingredients of Traditional Chinese Medicines and the active pharmaceutical ingredients for six major diseases-resistance, respectively. Comparative analysis showed there were obviously difference in metabolic profiles between the two buckwheat species, and further found 61 flavonoids and 94 non-flavonoids metabolites displayed significantly higher contents (≥2 fold) in Tartary buckwheat than in common buckwheat. Our results suggest that Tartary and common buckwheat seeds are rich in metabolites beneficial to human health, and non-flavonoids metabolites also contributed to Tartary buckwheat's higher health-promoting value than common buckwheat. This study provides valuable information for the development of new functional foods of Tartary buckwheat.

Thermal degradation of cyanidin-3-O-glucoside: Mechanism and toxicity of products.

The thermal degradation behavior of cyanidin-3-O-gluoside (Cy3G) in nitrogen and air was studied using thermogravimetric analysis (TGA), thermogravimetry-Fourier transform infrared spectroscopy (TG-FTIR) and pyrolysis-gas chromatography/mass spectrometry (Py-GCMS). The results show that the thermal degradation of Cy3G in nitrogen and in air can be divided into three steps. The total degradation rate was 63.09% in nitrogen and 99.42% in air, and the total activation energy (Ea) was 65.85 and 80.98 kJ·mol-1, respectively. The TG-FTIR analysis showed that Cy3G is significantly decomposed at 200-300 °C. The Py-GCMS analysis shows that the first step in the thermal degradation of Cy3G in nitrogen is the cleavage of glycosidic bonds to give cyanidin and glucoside. The glucoside and cyanidin then degrade further to give mainly low molecular weight compounds, together with furan derivatives, pyran derivatives and aromatic compounds. The phenols and furans found in the pyrolysis products are known to have a degree of toxicity.

Genome-wide investigation of bHLH genes and expression analysis under different biotic and abiotic stresses in Helianthus annuus L.

The basic helix-loop-helix (bHLH) transcription factors play important roles in many processes such as plant growth, metabolism and response to biotic/abiotic stresses. Sunflower (Helianthus annuus) is a major oil crop, cultivated throughout the world. However, no systematic characterization of bHLH gene members in sunflower (HabHLH) and their functions involved in drought, cadmium tolerance and Orobanche cumana resistance has been reported yet. In this study, 183 HabHLH genes were identified and named according to their chromosomal locations. We classified these proteins into 21 subfamilies by phylogenetic tree analysis. Subsequently, DNA-binding patterns, sequence analysis, duplication analysis and gene structures were analyzed. All of the HabHLH genes were randomly distributed on 17 chromosomes, and 10 pairs of tandem duplicated genes and one pair of segmental duplicated genes were detected in the HabHLH family. Among the duplicated gene pairs, eight pairs of HabHLH genes suffer from positive selection. Moreover, qRT-PCR results revealed significant up-regulated expression of HabHLH024 gene in response to both abiotic (cadmium, drought) and biotic (Orobanche cumana) stresses, suggesting its important functions in response to different stresses. Therefore, HabHLH024 would be the potential candidate gene for the sunflower tolerance breeding.

Genome-wide investigation and expression analysis of membrane-bound fatty acid desaturase genes under different biotic and abiotic stresses in sunflower (Helianthus annuus L.).

Membrane-bound fatty acid desaturase (FAD) gene family plays crucial roles in regulation of fatty acid (FA) compositions in plants. Sunflower (Helianthus annuus L.) is an important oilseed crop in the world; however, no comprehensive study on exploring the role of FAD family in relation to stress tolerance in sunflower has been performed yet. In this study, we identified 40 putative FAD genes in H. annuus (HaFAD), which were unevenly distributed across 13 of the total 17 chromosomes. Phylogenetic analysis indicated that HaFAD genes were divided into four subfamilies, as supported by highly conserved gene structures and motifs. Collinearity analysis showed that tandem duplication events played a crucial role in the expansion of HaFAD gene family. In addition, tissue-specific expression showed that 32 HaFAD genes were widely expressed in various tissues or organs of sunflower. Furthermore, qRT-PCR results revealed significant expression changes of HaFAD genes in response to abiotic (cadmium, drought) and biotic (Orobanche cumana) stresses, suggesting their important functions in response to different stresses. Therefore, our results provide insights into HaFAD gene family in response to different stresses, and some specific up-regulated genes such as HaFAD3.2, HaADS8, HaFAD2.1, and HaADS9 would be the potential candidate genes for the sunflower tolerance breeding.

Diversity patterns of soil microbial communities in the Sophora flavescens rhizosphere in response to continuous monocropping.

BACKGROUND: Continuous monocropping can affect the physicochemical and biological characteristics of cultivated soil. Sophora flavescens is a valuable herbal medicine and sensitive to continuous monocropping. Currently, diversity patterns of soil microbial communities in soil continuous monocropping with S. flavescens have not been extensively elucidated. RESULTS: In this study, comparative 16S rDNA and internal transcribed spacer (ITS) MiSeq sequencing analyses were used to examine the taxonomic community structure and microbial diversity in nonrhizosphere soil (CK) and rhizosphere soils (SCC, TCC, and FCC) sampled from fields that had undergone two, three, and five years of continuous monocropping, respectively. Among the microbial communities, a decreased abundance of Acidobacteria and increased abundances of Proteobacteria and Bacteroidetes were found with the increase in monocropping years of S. flavescens. As the continuous monocropping time increased, the diversity of the bacterial community decreased, but that of fungi increased. Redundancy analysis also showed that among the properties of the rhizosphere soil, the available phosphorus, organic matter, total nitrogen, and sucrase had the greatest impacts on the diversity of the rhizosphere microbial community. Moreover, a biomarker for S. flavescens soil was also identified using the most differentially abundant bacteria and fungi in soil samples. CONCLUSIONS: Our study indicates that long-term monocropping exerted great impacts on microbial community distributions and soil physicochemical properties. The relationship between microbial community and physicochemical properties of rhizosphere soil would help clarify the side effects of continuous S. flavescens monocropping. Our study may aid in uncovering the theoretical basis underlying obstacles to continuous monocropping and provide better guidance for crop production.

Comparative transcriptome analysis provides insights into the molecular mechanism underlying double fertilization between self-crossed Solanum melongena and that hybridized with Solanum aethiopicum.

Wild relatives represent a source of variation for many traits of interest for eggplant (Solanum melongena) breeding, as well as for broadening its genetic base. However, interspecific hybridization with wild relatives has been barely used in eggplant breeding programs, and reproductive barriers have resulted in reduced seed numbers in interspecific combinations. The mechanism underlying this phenomenon remains unclear. We hybridized females of cultivated eggplant 177 (Solanum melongena) with males of wild relatives 53 and Y11 (Solanum aethiopicum). Self-crossed 177 was the control. The seed number per control fruit was significantly higher than that of the hybrids. Paraffin sections showed no significant difference between control and 177×53 and 177×Y11. Double fertilization began 4 days post-pollination. Sperm cells were fused with egg cells 6 days post-pollination. To understand the differences in molecular mechanisms underlying this process, transcriptomes of ovaries at 0, 4, and 6 days after self-crossing and hybridization were analyzed. We screened 22,311 differentially expressed genes (DEGs) between the control and hybrids 4 and 6 days post-pollination. A total of 497 DEGs were shared among all pollination combinations. These DEGs were enriched in plant hormone transduction, cell senescence, metabolism, and biosynthesis pathways. DEG clustering analysis indicated distinct expression patterns between the control and hybrids but not between the hybrids. The DEGs in hybrids involved secondary metabolic process, phenylpropanoid metabolic process, and carboxypeptidase activity, while those in the control involved xyloglucan metabolic process, auxin-activated signaling pathway, cell wall polysaccharide metabolic process, and xyloglucosyl transferase activity. Additionally, 1683 transcription factors, including members of the AP2-ERF, MYB, bHLH, and B3 families may play important roles in self-crossing and hybridization. Our results provide insights into the regulatory mechanisms underlying variations between ovaries of self-crossed and hybrid eggplants and a basis for future studies on crossbreeding Solanum and genetic mechanisms underlying double fertilization.

Genomic Survey of Tyrosine Kinases Repertoire in Electrophorus electricus With an Emphasis on Evolutionary Conservation and Diversification.

Tyrosine kinases (TKs) play key roles in the regulation of multicellularity in organisms and involved primarily in cell growth, differentiation, and cell-to-cell communication. Genome-wide characterization of TKs has been conducted in many metazoans; however, systematic information regarding this superfamily in Electrophorus electricus (electric eel) is still lacking. In this study, we identified 114 TK genes in the E electricus genome and investigated their evolution, molecular features, and domain architecture using phylogenetic profiling to gain a better understanding of their similarities and specificity. Our results suggested that the electric eel TK (EeTK) repertoire was shaped by whole-genome duplications (WGDs) and tandem duplication events. Compared with other vertebrate TKs, gene members in Jak, Src, and EGFR subfamily duplicated specifically, but with members lost in Eph, Axl, and Ack subfamily in electric eel. We also conducted an exhaustive survey of TK genes in genomic databases, identifying 1674 TK proteins in 31 representative species covering all the main metazoan lineages. Extensive evolutionary analysis indicated that TK repertoire in vertebrates tended to be remarkably conserved, but the gene members in each subfamily were very variable. Comparative expression profile analysis showed that electric organ tissues and muscle shared a similar pattern with specific highly expressed TKs (ie, epha7, musk, jak1, and pdgfra), suggesting that regulation of TKs might play an important role in specifying an electric organ identity from its muscle precursor. We further identified TK genes exhibiting tissue-specific expression patterns, indicating that members in TKs participated in subfunctionalization representing an evolutionary divergence required for the performance of different tissues. This work generates valuable information for further gene function analysis and identifying candidate TK genes reflecting their unique tissue-function specializations in electric eel.

Volatile Organic Compounds of Endophytic Burkholderia pyrrocinia Strain JK-SH007 Promote Disease Resistance in Poplar.

Volatile organic compounds (VOCs) play important roles in the regulation of plant growth and pathogen resistance. However, little is known about the influence of VOCs released from endophytic strains (Burkholderia pyrrocinia strain JK-SH007) on controlling pathogens or inducing systemic resistance in poplar. In this study, we found that VOCs produced by strain JK-SH007 inhibit three poplar canker pathogens (Cytospora chrysosperma, Phomopsis macrospora, and Fusicoccum aesculi) and promote defense enzyme activity and malondialdehyde (MDA) and total phenol (TP) accumulation. Thirteen kinds of VOC components were identified using the solid-phase microextraction combined with gas chromatography-mass spectrometry method. Dimethyl disulfide (DMDS) accounted for the largest proportion of these VOCs. Treatments of poplar seedlings with different volumes of VOC standards (DMDS, benzothiazole, dimethylthiomethane, and phenylacetone) showed that DMDS had the greatest effects on various defense enzyme activities and MDA and TP accumulation. We also found that the inhibitory effect of the VOCs on the three pathogens was gradually enhanced with increasing standard volume. Moreover, the treatment of samples with DMDS significantly reduced the severity and development of the disease caused by three poplar canker pathogens. Comparative transcriptomics analysis of poplar seedlings treated with DMDS showed that there were 1,586 differentially expressed genes in the leaves and stems, and quantitative PCR showed that the gene expression trends were highly consistent with the transcriptome sequencing results. The most significant transcriptomic changes induced by VOCs were related to hormone signal transduction, transcriptional regulation of plant-pathogen interactions, and energy metabolism. Moreover, 137 transcription factors, including members of the ethylene response factor, NAC, WRKY, G2-like, and basic helix-loop-helix protein families, were identified to be involved in the VOC-induced process. This study elucidates the resistance induced by Burkholderia pyrrocinia strain JK-SH007 to poplar canker at the molecular level and can make possible a new method for the comprehensive prevention and control of poplar disease.