Genomic Resequencing Unravels the Genetic Basis of Domestication, Expansion, and Trait Improvement in Morus Atropurpurea.

Mulberry is an economically important plant in the sericulture industry and traditional medicine. However, the genetic and evolutionary history of mulberry remains largely unknown. Here, this work presents the chromosome-level genome assembly of Morus atropurpurea (M. atropurpurea), originating from south China. Population genomic analysis using 425 mulberry accessions reveal that cultivated mulberry is classified into two species, M. atropurpurea and M. alba, which may have originated from two different mulberry progenitors and have independent and parallel domestication in north and south China, respectively. Extensive gene flow is revealed between different mulberry populations, contributing to genetic diversity in modern hybrid cultivars. This work also identifies the genetic architecture of the flowering time and leaf size. In addition, the genomic structure and evolution of sex-determining regions are identified. This study significantly advances the understanding of the genetic basis and domestication history of mulberry in the north and south, and provides valuable molecular markers of desirable traits for mulberry breeding.

Role of long non-coding RNA in regulatory network response to Candidatus Liberibacter asiaticus in citrus.

Long non-coding RNAs (lncRNAs) serve as crucial regulators in plant response to various diseases, while none have been systematically identified and characterized in response to citrus Huanglongbing (HLB) caused by Candidatus Liberibacter asiaticus (CLas) bacteria. Here, we comprehensively investigated the transcriptional and regulatory dynamics of the lncRNAs in response to CLas. Samples were collected from leaf midribs of CLas- and mock-inoculated HLB-tolerant rough lemon (Citrus jambhiri) and HLB-sensitive sweet orange (C. sinensis) at week 0, 7, 17, and 34 following inoculation using CLas+ budwood of three biological replicates in the greenhouse. A total of 8,742 lncRNAs, including 2,529 novel lncRNAs, were identified from RNA-seq data with rRNA-removed from strand-specific libraries. Genomic variation analyses of conserved lncRNAs from 38 citrus accessions showed that 26 single nucleotide polymorphisms (SNPs) were significantly correlated with HLB. In addition, lncRNA-mRNA weighted gene co-expression network analysis (WGCNA) showed a significant module correlated with CLas-inoculation in rough lemon. Notably, the most significant LNC_28805 and multiple co-expressed genes related to plant defense in the module were targeted by miRNA5021, suggesting that LNC28805 might compete with endogenous miR5021 to maintain the homeostasis of immune gene expression levels. Candidate WRKY33 and SYP121 genes targeted by miRNA5021 were identified as two key hub genes interacting with bacteria pathogen response genes based on the prediction of protein-protein interaction (PPI) network. These two genes were also found within HLB-associated QTL in linkage group 6. Overall, our findings provide a reference for a better understanding of the role of lncRNAs involved in citrus HLB regulation.

Integration of Transcriptome and Metabolome Reveals the Formation Mechanism of Red Stem in Prunus mume.

Prunus mume var. purpurea, commonly known as "Red Bone", is a special variety with pink or purple-red xylem. It is famous due to gorgeous petals and delightful aromas, playing important roles in urban landscaping. The regulation mechanism of color formation in P. mume var. purpurea stem development is unclear. Here, we conducted a comprehensive analysis of transcriptome and metabolome in WYY ('Wuyuyu' accession, red stem) and FLE ('Fei Lve' accession, green stem), and found a total of 256 differential metabolites. At least 14 anthocyanins were detected in WYY, wherein cyanidin 3,5-O-diglucoside and peonidin3-O-glucoside were significantly accumulated through LC-MS/MS analysis. Transcriptome data showed that the genes related to flavonoid-anthocyanin biosynthesis pathways were significantly enriched in WYY. The ratio of dihydroflavonol 4-reductase (DFR) and flavonol synthase (FLS) expression levels may affect metabolic balance in WYY, suggesting a vital role in xylem color formation. In addition, several transcription factors were up-regulated, which may be the key factors contributing to transcriptional changes in anthocyanin synthesis. Overall, the results provide a reference for further research on the molecular mechanism of xylem color regulation in P. mume and lay a theoretical foundation for cultivating new varieties.

The chromosome-level genome provides insight into the molecular mechanism underlying the tortuous-branch phenotype of Prunus mume.

Plant with naturally twisted branches is referred to as a tortuous-branch plant, which have extremely high ornamental value due to their zigzag shape and the natural twisting of their branches. Prunus mume is an important woody ornamental plant. However, the molecular mechanism underlying this unique trait in Prunus genus is unknown. Here, we present a chromosome-level genome assembly of the cultivated P. mume var. tortuosa created using Oxford Nanopore combined with Hi-C scaffolding, which resulted in a 237.8 Mb genome assembly being anchored onto eight pseudochromosomes. Molecular dating indicated that P. mume is the most recently differentiated species in Prunus. Genes associated with cell division, development and plant hormones play essential roles in the formation of tortuous branch trait. A putative regulatory pathway for the tortuous branch trait was constructed based on gene expression levels. Furthermore, after transferring candidate PmCYCD genes into Arabidopsis thaliana, we found that seedlings overexpressing these genes exhibited curled rosette leaves. Our results provide insights into the evolutionary history of recently differentiated species in Prunus genus, the molecular basis of stem morphology, and the molecular mechanism underlying the tortuous branch trait and highlight the utility of multi-omics in deciphering the properties of P. mume plant architecture.

Phenotypic and Genetic Variation in Morphophysiological Traits in Huanglongbing-Affected Mandarin Hybrid Populations.

Huanglongbing (HLB) caused by 'Candidatus Liberibacter asiaticus' (CLas) is the most costly disease for the global citrus industry. Currently, no effective tools have been found to control HLB. Most commercial citrus varieties are susceptible to HLB, though some citrus hybrid cultivars have reduced sensitivity to the disease. Citrus breeding populations contain a large diversity of germplasm, with thousands of unique genotypes exhibiting a broad range of phenotypes. Understanding phenotypic variation and genetic inheritance in HLB-affected mandarin hybrid populations are crucial for breeding tolerant citrus varieties. In this study, we assessed 448 diverse mandarin hybrids coming from 30 crosses, and 45 additional accessions. For HLB tolerance, we measured HLB severity visual score and CLas titers by qPCR. We also measured seven morphophysiological traits indirectly related to HLB tolerance with leaf area index (LAI), leaf area (LA), leaf mass per area (LMA), photosystem II parameters (Fv/Fo, Fv/Fm), and photochemical performance index (PIabs). By estimating the genetic variation in five half-sib families, we estimated the heritability of phenotypic traits and found a significant genetic effect on HLB visual score and photosynthesis parameters, which indicates opportunities for the genetic improvement of HLB tolerance. In addition, although it is easy to identify infected trees based on HLB symptomatic leaves, visually phenotyping whole trees can be difficult and inconsistent due to the interpersonal subjectivity of characterization. We investigated their relationships and found that LAI was highly correlated with HLB score, with correlation coefficients of r = 0.70 and r = 0.77 for the whole population and five half-sib families, respectively. Photochemical parameters showed significant correlation with HLB severity and responded differentially with the side of the canopy. Our study suggests that LAI and photochemical parameters could be used as a rapid and cost-effective method to evaluate HLB tolerance and inheritance in citrus breeding programs.

Histone Acetyltransferases and Deacetylases Are Required for Virulence, Conidiation, DNA Damage Repair, and Multiple Stresses Resistance of Alternaria alternata.

Histone acetylation, which is critical for transcriptional regulation and various biological processes in eukaryotes, is a reversible dynamic process regulated by HATs and HDACs. This study determined the function of 6 histone acetyltransferases (HATs) (Gcn5, RTT109, Elp3, Sas3, Sas2, Nat3) and 6 histone deacetylases (HDACs) (Hos2, Rpd3, Hda1, Hos3, Hst2, Sir2) in the phytopathogenic fungus Alternaria alternata by analyzing targeted gene deletion mutants. Our data provide evidence that HATs and HDACs are both required for mycelium growth, cell development and pathogenicity as many gene deletion mutants (ΔGcn5, ΔRTT109, ΔElp3, ΔSas3, ΔNat3, ΔHos2, and ΔRpd3) displayed reduced growth, conidiation or virulence at varying degrees. In addition, HATs and HDACs are involved in the resistance to multiple stresses such as oxidative stress (Sas3, Gcn5, Elp3, RTT109, Hos2), osmotic stress (Sas3, Gcn5, RTT109, Hos2), cell wall-targeting agents (Sas3, Gcn5, Hos2), and fungicide (Gcn5, Hos2). ΔGcn5, ΔSas3, and ΔHos2 displayed severe growth defects on sole carbon source medium suggesting a vital role of HATs and HDACs in carbon source utilization. More SNPs were generated in ΔGcn5 in comparison to wild-type when they were exposed to ultraviolet ray. Moreover, ΔRTT109, ΔGcn5, and ΔHos2 showed severe defects in resistance to DNA-damaging agents, indicating the critical role of HATs and HDACs in DNA damage repair. These phenotypes correlated well with the differentially expressed genes in ΔGcn5 and ΔHos2 that are essential for carbon sources metabolism, DNA damage repair, ROS detoxification, and asexual development. Furthermore, Gcn5 is required for the acetylation of H3K4. Overall, our study provides genetic evidence to define the central role of HATs and HDACs in the pathological and biological functions of A. alternata.

Identification of the PmWEEP locus controlling weeping traits in Prunus mume through an integrated genome-wide association study and quantitative trait locus mapping.

Weeping Prunus mume (mei) has long been cultivated in East Asia for its specific ornamental value. However, little is known about the regulatory mechanism of the weeping trait in mei, which limits molecular breeding for the improvement of weeping-type cultivars. Here, we quantified the weeping trait in mei using nested phenotyping of 214 accessions and 342 F1 hybrids. Two major associated loci were identified from the genome-wide association study (GWAS), which was conducted using 3,014,409 single nucleotide polymorphisms (SNPs) derived from resequencing, and 8 QTLs and 55 epistatic loci were identified from QTL mapping using 7,545 specific lengths amplified fragment (SLAF) markers. Notably, an overlapping PmWEEP major QTL was fine mapped within a 0.29 Mb region on chromosome 7 (Pa7), and a core SNP locus closely associated with the weeping trait was screened and validated. Furthermore, a total of 22 genes in the PmWEEP QTL region were expressed in weeping or upright mei based on RNA-seq analysis. Among them, only a novel gene (Pm024213) containing a thioredoxin (Trx) domain was found to be close to the core SNP and specifically expressed in buds and branches of weeping mei. Co-expression analysis of Pm024213 showed that most of the related genes were involved in auxin and lignin biosynthesis. These findings provide insights into the regulatory mechanism of the weeping trait and effective molecular markers for molecular-assisted breeding in Prunus mume.

Genome-wide identification, characterisation, and evolution of ABF/AREB subfamily in nine Rosaceae species and expression analysis in mei (Prunus mume).

Rosaceae is an important family containing some of the highly evolved fruit and ornamental plants. Abiotic stress responses play key roles in the seasonal growth and development of plants. However, the molecular basis of stress responses remains largely unknown in Rosaceae. Abscisic acid (ABA) is a stress hormone involving abiotic stress response pathways. The ABRE-binding factor/ABA-responsive element-binding protein (ABF/AREB) is a subfamily of the basic domain/leucine zipper (bZIP) transcription factor family. It plays an important role in the ABA-mediated signaling pathway. Here, we analyzed the ABF/AREB subfamily genes in nine Rosaceae species. A total of 64 ABF/AREB genes were identified, including 18, 28, and 18 genes in the Rosoideae, Amygdaloideae, and Maloideae traditional subfamilies, respectively. The evolutionary relationship of the ABF/AREB subfamily genes was studied through the phylogenetic analysis, the gene structure and conserved motif composition, Ka/Ks values, and interspecies colinearity. These gene sets were clustered into four groups. In the Prunus ABF/AREB (PmABF) promoters, several cis-elements related to light, hormone, and abiotic stress response were predicted. PmABFs expressed in five different tissues, except PmABF5, which expressed only in buds. In the dormancy stages, PmABF1, 2, 5 and 7 showed differential expression. The expression of PmABF3, 4 and 6 was positively correlated with the ABA concentration. Except for PmABF5, all the PmABFs were sensitive to ABA. Several ABRE elements were contained in the promoters of PmABF1, 3, 6, 7. Based on the findings of our study, we speculate that PmABFs may play a role in flower bud dormancy in P. mume.

Comparative gene expression analysis reveals that multiple mechanisms regulate the weeping trait in Prunus mume.

Prunus mume (also known as Mei) is an important ornamental plant that is popular with Asians. The weeping trait in P. mume has attracted the attention of researchers for its high ornamental value. However, the formation of the weeping trait of woody plants is a complex process and the molecular basis of weeping stem development is unclear. Here, the morphological and histochemical characteristics and transcriptome profiles of upright and weeping stems from P. mume were studied. Significant alterations in the histochemical characteristics of upright and weeping stems were observed, and the absence of phloem fibres and less xylem in weeping stems might be responsible for their inability to resist gravity and to grow downward. Transcriptome analysis showed that differentially expressed genes (DEGs) were enriched in phenylpropanoid biosynthesis and phytohormone signal transduction pathways. To investigate the differential responses to hormones, upright and weeping stems were treated with IAA (auxin) and GA3 (gibberellin A3), respectively, and the results revealed that weeping stems had a weaker IAA response ability and reduced upward bending angles than upright stems. On the contrary, weeping stems had increased upward bending angles than upright stems with GA3 treatment. Compared to upright stems, interestingly, DEGs associated with diterpenoid biosynthesis and phenylpropanoid biosynthesis were significantly enriched after being treated with IAA, and expression levels of genes associated with phenylpropanoid biosynthesis, ABC transporters, glycosylphosphatidylinositol (GPI)-anchor biosynthesis were altered after being treated with GA3 in weeping stems. Those results reveal that multiple molecular mechanisms regulate the formation of weeping trait in P. mume, which lays a theoretical foundation for the cultivation of new varieties.

Transcriptome profiles reveal that gibberellin-related genes regulate weeping traits in crape myrtle.

Plant architecture includes vital traits that influence and benefit crops, and economically important trees. Different plant architectures provide natural beauty. Weeping ornamental plants are aesthetically appealing to people. The regulatory mechanism controlling the weeping trait is poorly understood in crape myrtle. To investigate the weeping trait mechanism, transcriptional profiling of different organs in weeping and upright crape myrtle was performed based on phenotype. Phenotypic and histological analyses demonstrated that endodermal cells were absent, and that new shoot phenotypes could be rescued by the GA3 treatment of weeping plants. The transcriptional analysis and coexpression network analysis (WGCNA) of differentially expressed genes indicated that GA synthesis and signal transduction pathways play a role in weeping traits. When the expression level of a negative element of GA signaling, LfiGRAS1, was reduced by virus-induced gene silencing (VIGS), new branches grew in infected plants in a negatively geotropic manner. An integrated analysis implied that GA had a strong influence on weeping crape myrtle by interacting with other factors. This study helps to elucidate the mechanism governing the weeping trait and can improve the efficiency of breeding in Lagerstroemia.

Corrigendum: A Computational Model for Inferring QTL Control Networks Underlying Developmental Covariation.

[This corrects the article DOI: 10.3389/fpls.2019.01557.].

Genome-wide identification, characterization and expression analysis of the HD-Zip gene family in the stem development of the woody plant Prunus mume.

The homeodomain-leucine zipper (HD-Zip) gene family, a group of plant-specific transcriptional factors (TFs), participates in regulating growth, development, and environmental responses. However, the characteristics and biological functions of HD-Zip genes in Prunus mume, which blooms in late winter or early spring, have not been reported. In this study, 32 HD-Zip genes, named PmHB1-PmHB32 based on their chromosomal positions, were identified in the genome of P. mume. These genes are distributed among seven chromosomes and are phylogenetically clustered into four major groups. Gene structure and motif composition were mostly conserved in each group. The Ka/Ks ratios showed that purifying selection has played a leading role in the long-term evolution of the genes, which maintained the function of this family. MicroRNA target site prediction indicated that the genes of the HD-Zip III subfamily may be regulated by miR165/166. Expression pattern analysis showed that the 32 genes were differentially expressed across five different tissues (leaf, flower bud, stem, fruit, and root) and at different stages of stem and leaf-bud development, suggesting that 10 of the genes may play important roles in stem development. Protein-protein interaction predictions showed that the subfamily III genes may regulate vascular development and shoot apical meristem (SAM) maintenance. Promoter analysis showed that the HD-Zip III genes might be involved in responses to light, hormones, and abiotic stressors and stem development. Taken together, our results provide an overview of the HD-Zip family in P. mume and lay the foundation for the molecular breeding of woody ornamental plants.

Identification and comparative analysis of the CIPK gene family and characterization of the cold stress response in the woody plant Prunus mume.

Prunus mume is an important ornamental woody plant that grows in tropical and subtropical regions. Freezing stress can adversely impact plant productivity and limit the expansion of geographical locations. Understanding cold-responsive genes could potentially bring about the development of new ways to enhance plant freezing tolerance. Members of the serine/threonine protein kinase (CIPK) gene family play important roles in abiotic stress. However, the function of CIPK genes in P. mume remains poorly defined. A total of 16 CIPK genes were first identified in P. mume. A systematic phylogenetic analysis was conducted in which 253 CIPK genes from 12 species were divided into three groups. Furthermore, we analysed the chromosomal locations, molecular structures, motifs and domains of CIPK genes in P. mume. All of the CIPK sequences had NAF domains and promoter regions containing cis-acting regulatory elements of the related stress response. Three PmCIPK genes were identified as Pmu-miR172/167-targeted sites. Transcriptome data showed that most PmCIPK genes presented tissue-specific and time-specific expression profiles. Nine genes were highly expressed in flower buds in December and January, and 12 genes were up-regulated in stems in winter. The expression levels of 12 PmCIPK genes were up-regulated during cold stress treatment confirmed by qRT-PCR. Our study improves understanding of the role of the PmCIPK gene family in the low temperature response in woody plants and provides key candidate genes and a theoretical basis for cold resistance molecular-assisted breeding technology in P. mume.

A Computational Model for Inferring QTL Control Networks Underlying Developmental Covariation.

How one trait developmentally varies as a function of others shapes a spectrum of biological phenomena. Despite its importance to trait dissection, the understanding of whether and how genes mediate such developmental covariation is poorly understood. We integrate developmental allometry equations into the functional mapping framework to map specific QTLs that govern the correlated development of different traits. Based on evolutionary game theory, we assemble and contextualize these QTLs into an intricate but organized network coded by bidirectional, signed, and weighted QTL-QTL interactions. We use this approach to map shoot height-diameter allometry QTLs in an ornamental woody species, mei (Prunus mume). We detect "pioneering" QTLs (piQTLs) and "maintaining" QTLs (miQTLs) that determine how shoot height varies with diameter and how shoot diameter varies with height, respectively. The QTL networks inferred can visualize how each piQTL regulates others to promote height growth at a cost of diameter growth, how miQTL regulates others to benefit radial growth at a cost of height growth, and how piQTLs and miQTLs regulate each other to form a pleiotropic web of primary and secondary growth in trees. Our approach provides a unique gateway to explore the genetic architecture of developmental covariation, a widespread phenomenon in nature.

Comparative Genomic and Transcriptomic Analyses of Family-1 UDP Glycosyltransferase in Prunus Mume.

Glycosylation mediated by Family-1 UDP-glycosyltransferases (UGTs) plays crucial roles in plant growth and adaptation to various stress conditions. Prunus mume is an ideal crop for analyzing flowering for its early spring flowering characteristics. Revealing the genomic and transcriptomic portfolio of the UGT family in P. mume, a species in which UGTs have not yet been investigated, is therefore important. In this study, 130 putative UGT genes were identified and phylogenetically clustered into 14 groups. These PmUGTs were distributed unevenly across eight chromosomes and 32 tandem duplication and 8 segmental duplication pairs were revealed. A highly conserved intron insertion event was revealed on the basis of intron/exon patterns within PmUGTs. According to RNA-seq data, these PmUGTs were specifically expressed in different tissues and during the bud dormancy process. In addition, we confirmed the differential expression of some representative genes in response to abscisic acid treatment. Our results will provide important information on the UGT family in P. mume that should aid further characterization of their biological roles in response to environmental stress.

Genome-Wide Analysis of the NAC Transcription Factor Gene Family Reveals Differential Expression Patterns and Cold-Stress Responses in the Woody Plant Prunus mume.

NAC transcription factors (TFs) participate in multiple biological processes, including biotic and abiotic stress responses, signal transduction and development. Cold stress can adversely impact plant growth and development, thereby limiting agricultural productivity. Prunus mume, an excellent horticultural crop, is widely cultivated in Asian countries. Its flower can tolerate freezing-stress in the early spring. To investigate the putative NAC genes responsible for cold-stress, we identified and analyzed 113 high-confidence PmNAC genes and characterized them by bioinformatics tools and expression profiles. These PmNACs were clustered into 14 sub-families and distributed on eight chromosomes and scaffolds, with the highest number located on chromosome 3. Duplicated events resulted in a large gene family; 15 and 8 pairs of PmNACs were the result of tandem and segmental duplicates, respectively. Moreover, three membrane-bound proteins (PmNAC59/66/73) and three miRNA-targeted genes (PmNAC40/41/83) were identified. Most PmNAC genes presented tissue-specific and time-specific expression patterns. Sixteen PmNACs (PmNAC11/19/20/23/41/48/58/74/75/76/78/79/85/86/103/111) exhibited down-regulation during flower bud opening and are, therefore, putative candidates for dormancy and cold-tolerance. Seventeen genes (PmNAC11/12/17/21/29/42/30/48/59/66/73/75/85/86/93/99/111) were highly expressed in stem during winter and are putative candidates for freezing resistance. The cold-stress response pattern of 15 putative PmNACs was observed under 4 °C at different treatment times. The expression of 10 genes (PmNAC11/20/23/40/42/48/57/60/66/86) was upregulated, while 5 genes (PmNAC59/61/82/85/107) were significantly inhibited. The putative candidates, thus identified, have the potential for breeding the cold-tolerant horticultural plants. This study increases our understanding of functions of the NAC gene family in cold tolerance, thereby potentially intensifying the molecular breeding programs of woody plants.

Transcriptome Profiles Reveal the Crucial Roles of Hormone and Sugar in the Bud Dormancy of Prunus mume.

Bud dormancy transition is a vital developmental process for perennial plant survival. The process is precisely regulated by diverse endogenous genetic factors and environmental cues, but the mechanisms are not yet fully understood. Prunus mume is an ideal crop for bud dormancy analysis because of its early spring-flowering characteristics and small sequenced genome. Here, we analyzed the transcriptome profiles at the three endodormancy stages and natural flush stage using RNA sequencing combined with phytohormone and sugar content measurements. Significant alterations in hormone contents and carbohydrate metabolism have been observed, and α-amylases, Glucan Hydrolase Family 17 and diphosphate-glycosyltransferase family might play crucial roles in the interactions between hormones and sugars. The following hypothetical model for understanding the molecular mechanism of bud dormancy in Prunus mume is proposed: low temperatures exposure induces the significant up-regulation of eight C-repeat binding factor genes, which directly promotes all six dormancy-associated MADS-box genes, resulting in dormancy establishment. The prolonged cold and/or subsequently increasing temperature then decreases the expression levels of these two gene families, which alleviates the inhibition of FLOWERING LOCUS T and reopens the growth-promoting pathway, resulting in dormancy release and the initiation of the bud break process.

Crosstalk of PmCBFs and PmDAMs Based on the Changes of Phytohormones under Seasonal Cold Stress in the Stem of Prunus mume.

Plants facing the seasonal variations always need a growth restraining mechanism when temperatures turn down. C-repeat binding factor (CBF) genes work essentially in the cold perception. Despite lots of researches on CBFs, the multiple crosstalk is still interesting on their interaction with hormones and dormancy-associated MADS (DAM) genes in the growth and dormancy control. Therefore, this study highlights roles of PmCBFs in cold-induced dormancy from different orgens. And a sense-response relationship between PmCBFs and PmDAMs is exhibited in this process, jointly regulated by six PmCBFs and PmDAM4-6. Meantime, GA3 and ABA showed negative and positive correlation with PmCBFs expression levels, respectively. We also find a high correlation between IAA and PmDAM1-3. Finally, we display the interaction mode of PmCBFs and PmDAMs, especially PmCBF1-PmDAM1. These results can disclose another view of molecular mechanism in plant growth between cold-response pathway and dormancy regulation together with genes and hormones.