New Genes Interacted with Recent Whole Genome Duplicates in the Fast Stem Growth of Bamboos

Coexpression Regulation of New and Ancient Genes in the Dynamic Transcriptome Landscape of Stem and Rhizome Development in "Bainianzhe"-An Ancient Chinese Sugarcane Variety Ratooned for Nearly 300 Years

The sucrose yield in sugarcane largely depends on stem morphology, including length, diameter and sugar content, making sugarcane stem a key trait in breeding. The "Bainianzhe" variety from Songxi County, Fujian Province, possesses both aerial stems and rhizomes, providing a unique model for studying stem development. We performed a spatiotemporal transcriptomic analysis of the base, middle and apical sections of both aerial stems and rhizomes. The analysis categorized transcriptomes by developmental stage-base, middle and apical-rather than environmental differences. Apical segments were enriched with genes related to cell proliferation, while base segments were linked to senescence and fibrosis. Gene regulatory networks revealed key TFs involved in stem development. Orphan genes may be involved in rhizome development through coexpression networks. Plant hormones, especially genes involved in ABA and GAs synthesis, were highly expressed in rhizomes. Thiamine-related genes were also more prevalent in rhizomes. Furthermore, the apical segments of rhizomes enriched in photosynthesis-related genes suggest adaptations to light exposure. Low average temperatures in Songxi have led to unique cold acclimation in Bainianzhe, with rhizomes showing higher expression of genes linked to unsaturated fatty acid synthesis and cold-responsive calcium signalling. This indicates that rhizomes may have enhanced cold tolerance, aiding in the plant's overwintering success.

De novo gene integration into regulatory networks via interaction with conserved genes in peach

De novo genes can evolve "from scratch" from noncoding sequences, acquiring novel functions in organisms and integrating into regulatory networks during evolution to drive innovations in important phenotypes and traits. However, identifying de novo genes is challenging, as it requires high-quality genomes from closely related species. According to the comparison with nine closely related Prunus genomes, we determined at least 178 de novo genes in P. persica "baifeng". The distinct differences were observed between de novo and conserved genes in gene characteristics and expression patterns. Gene ontology enrichment analysis suggested that Type I de novo genes originated from sequences related to plastid modification functions, while Type II genes were inferred to have derived from sequences related to reproductive functions. Finally, transcriptome sequencing across different tissues and developmental stages suggested that de novo genes have been evolutionarily recruited into existing regulatory networks, playing important roles in plant growth and development, which was also supported by WGCNA analysis and quantitative trait loci data. This study lays the groundwork for future research on the origins and functions of genes in Prunus and related taxa.

Spatiotemporal transcriptome atlas reveals gene regulatory patterns during the organogenesis of the rapid growing bamboo shoots

Bamboo with its remarkable growth rate and economic significance, offers an ideal system to investigate the molecular basis of organogenesis in rapidly growing plants, particular in monocots, where gene regulatory networks governing the maintenance and differentiation of shoot apical and intercalary meristems remain a subject of controversy. We employed both spatial and single-nucleus transcriptome sequencing on 10× platform to precisely dissect the gene functions in various tissues and early developmental stages of bamboo shoots. Our comprehensive analysis reveals distinct cell trajectories during shoot development, uncovering critical genes and pathways involved in procambium differentiation, intercalary meristem formation, and vascular tissue development. Spatial and temporal expression patterns of key regulatory genes, particularly those related to hormone signaling and lipid metabolism, strongly support the hypothesis that intercalary meristem origin from surrounded parenchyma cells. Specific gene expressions in intercalary meristem exhibit regular and dispersed distribution pattern, offering clues for understanding the intricate molecular mechanisms that drive the rapid growth of bamboo shoots. The single-nucleus and spatial transcriptome analysis reveal a comprehensive landscape of gene activity, enhancing the understanding of the molecular architecture of organogenesis and providing valuable resources for future genomic and genetic studies relying on identities of specific cell types.

Single-cell phylotranscriptomics of developmental and cell type evolution

Single-cell phylotranscriptomics is an emerging tool to reveal the molecular and cellular mechanisms of evolution. We summarize its utility in studying the hourglass pattern of ontogenetic evolution and for understanding the evolutionary history of cell types. The developmental hourglass model suggests that the mid-embryonic stage is the most conserved period of development across species, which is supported by morphological and molecular studies. Single-cell phylotranscriptomic analysis has revealed previously underappreciated heterogeneity in transcriptome ages among lineages and cell types throughout development, and has identified the lineages and tissues that drive the whole-organism hourglass pattern. Single-cell transcriptome age analyses also provide important insights into the origin of germ layers, the different selective forces on tissues during adaptation, and the evolutionary relationships between cell types.

Nitrogen fertilization practices alter microbial communities driven by clonal integration in Moso bamboo

Nitrogen (N) fertilization is crucial for maintaining plant productivity. Clonal plants can share resources between connected ramets through clonal integration influencing microbial communities and regulating soil biogeochemical cycling, especially in the rhizosphere. However, the effect of various N fertilization practices on microbial communities in the rhizosphere of clonal ramets remain unknown. In this study, clonal fragments of Moso bamboo (Phyllostachys edulis), consisting of a parent ramet, an offspring ramet, and an interconnecting rhizome, were established in the field. NH4NO3 solution was applied to the parent, offspring ramets or rhizomes to investigate the effect of fertilization practices on the structure and function of rhizosphere microbial communities. The differences in N availability, microbial biomass and community composition, and abundance of nitrifying genes among rhizosphere soils of ramets gradually decreased during the rapid growth of Moso bamboo, irrespective of fertilization practice. The soil N availability variation, particularly in NO3-, caused by fertilization practices altered the rhizosphere microbial community. Soil N availability and stable microbial biomass N in parent fertilization were the highest, being 9.0 % and 18.7 %, as well as 60.8 % and 90.4 % higher than rhizome and offspring fertilizations, respectively. The microbial network nodes and links in rhizome fertilization were 1.8 and 7.5 times higher than in parent and offspring fertilization, respectively. However, the diversity of bacterial community and abundance of nitrifying and denitrifying genes were the highest in offspring fertilization among three practices, which may be associated with increased N loss. Collectively, the rhizosphere microbial community characteristics depended on fertilization practices by altering the clonal integration of N in Moso bamboo. Parent and rhizome fertilization were favorable for N retention in plant-soil system and resulted in more stable microbial functions than offspring fertilization. Our findings provide new insights into precision fertilization for the sustainable Moso bamboo forest management.

Genome assemblies of 11 bamboo species highlight diversification induced by dynamic subgenome dominance

Polyploidy (genome duplication) is a pivotal force in evolution. However, the interactions between parental genomes in a polyploid nucleus, frequently involving subgenome dominance, are poorly understood. Here we showcase analyses of a bamboo system (Poaceae: Bambusoideae) comprising a series of lineages from diploid (herbaceous) to tetraploid and hexaploid (woody), with 11 chromosome-level de novo genome assemblies and 476 transcriptome samples. We find that woody bamboo subgenomes exhibit stunning karyotype stability, with parallel subgenome dominance in the two tetraploid clades and a gradual shift of dominance in the hexaploid clade. Allopolyploidization and subgenome dominance have shaped the evolution of tree-like lignified culms, rapid growth and synchronous flowering characteristic of woody bamboos as large grasses. Our work provides insights into genome dominance in a remarkable polyploid system, including its dependence on genomic context and its ability to switch which subgenomes are dominant over evolutionary time.

Expectations of duplicate gene retention under the gene duplicability hypothesis

BACKGROUND

Gene duplication is an important process in evolution. What causes some genes to be retained after duplication and others to be lost is a process not well understood. The most prevalent theory is the gene duplicability hypothesis, that something about the function and number of interacting partners (number of subunits of protein complex, etc.), determines whether copies have more opportunity to be retained for long evolutionary periods. Some genes are also more susceptible to dosage balance effects following WGD events, making them more likely to be retained for longer periods of time. One would expect these processes that affect the retention of duplicate copies to affect the conditional probability ratio after consecutive whole genome duplication events. The probability that a gene will be retained after a second whole genome duplication event (WGD2), given that it was retained after the first whole genome duplication event (WGD1) versus the probability a gene will be retained after WGD2, given it was lost after WGD1 defines the probability ratio that is calculated.

RESULTS

Since duplicate gene retention is a time heterogeneous process, the time between the events (t1) and the time since the most recent event (t2) are relevant factors in calculating the expectation for observation in any genome. Here, we use a survival analysis framework to predict the probability ratio for genomes with different values of t1 and t2 under the gene duplicability hypothesis, that some genes are more susceptible to selectable functional shifts, some more susceptible to dosage compensation, and others only drifting. We also predict the probability ratio with different values of t1 and t2 under the mutational opportunity hypothesis, that probability of retention for certain genes changes in subsequent events depending upon how they were previously retained. These models are nested such that the mutational opportunity model encompasses the gene duplicability model with shifting duplicability over time. Here we present a formalization of the gene duplicability and mutational opportunity hypotheses to characterize evolutionary dynamics and explanatory power in a recently developed statistical framework.

CONCLUSIONS

This work presents expectations of the gene duplicability and mutational opportunity hypotheses over time under different sets of assumptions. This expectation will enable formal testing of processes leading to duplicate gene retention.

How to utilize comparative transcriptomics to dissect morphological diversity in plants

Comparative transcriptomics has emerged as a powerful approach that allows us to unravel the genetic basis of organ morphogenesis and its diversification processes during evolution. However, the application of comparative transcriptomics in studying plant morphological diversity addresses challenges such as identifying homologous gene pairs, selecting appropriate developmental stages for comparison, and extracting biologically meaningful networks. Methods such as phylostratigraphy, clustering, and gene co-expression networks are explored to identify functionally equivalent genes, align developmental stages, and uncover gene regulatory relationships. In the current review, we highlight the importance of these approaches in overcoming the complexity of plant genomes, the impact of heterochrony on stage alignment, and the integration of gene networks with additional data for a comprehensive understanding of morphological evolution.

An Entropy-Balanced Orthogonal Learning Bamboo Forest Growth Optimization Algorithm with Quasi-Affine Transformation Evolutionary and Its Application in Capacitated Vehicle Routing Problem

The bamboo forest growth optimization (BFGO) algorithm combines the characteristics of the bamboo forest growth process with the optimization course of the algorithm. The algorithm performs well in dealing with optimization problems, but its exploitation ability is not outstanding. Therefore, a new heuristic algorithm named orthogonal learning quasi-affine transformation evolutionary bamboo forest growth optimization (OQBFGO) algorithm is proposed in this work. This algorithm combines the quasi-affine transformation evolution algorithm to expand the particle distribution range, a process of entropy increase that can significantly improve particle searchability. The algorithm also uses an orthogonal learning strategy to accurately aggregate particles from a chaotic state, which can be an entropy reduction process that can more accurately perform global development. OQBFGO algorithm, BFGO algorithm, quasi-affine transformation evolutionary bamboo growth optimization (QBFGO) algorithm, orthogonal learning bamboo growth optimization (OBFGO) algorithm, and three other mature algorithms are tested on the CEC2017 benchmark function. The experimental results show that the OQBFGO algorithm is superior to the above algorithms. Then, OQBFGO is used to solve the capacitated vehicle routing problem. The results show that OQBFGO can obtain better results than other algorithms.

Review: WRKY transcription factors: Understanding the functional divergence

WRKY transcription factors (TFs) play crucial roles in the growth and development of plants and their response to environmental changes. WRKY TFs have been detected in sequenced plant genomes. The functions and regulatory networks of many WRKY TFs, especially from Arabidopsis thaliana (AtWRKY TFs), have been revealed, and the origin of WRKY TFs in plants is clear. Nonetheless, the relationship between WRKY TFs function and classification is unclear. Furthermore, the functional divergence of homologous WRKY TFs in plants is unclear. In this review, WRKY TFs were explored based on WRKY-related literature published from 1994 to 2022. WRKY TFs were identified in 234 species at the genome and transcriptome levels. The biological functions of ∼ 71 % of AtWRKY TFs were uncovered. Although functional divergence occurred in homologous WRKY TFs, different WRKY TF groups had no preferential function.

Identifying candidate de novo genes expressed in the somatic female reproductive tract of Drosophila melanogaster

Most eukaryotic genes have been vertically transmitted to the present from distant ancestors. However, variable gene number across species indicates that gene gain and loss also occurs. While new genes typically originate as products of duplications and rearrangements of preexisting genes, putative de novo genes-genes born out of ancestrally nongenic sequence-have been identified. Previous studies of de novo genes in Drosophila have provided evidence that expression in male reproductive tissues is common. However, no studies have focused on female reproductive tissues. Here we begin addressing this gap in the literature by analyzing the transcriptomes of 3 female reproductive tract organs (spermatheca, seminal receptacle, and parovaria) in 3 species-our focal species, Drosophila melanogaster-and 2 closely related species, Drosophila simulans and Drosophila yakuba, with the goal of identifying putative D. melanogaster-specific de novo genes expressed in these tissues. We discovered several candidate genes, located in sequence annotated as intergenic. Consistent with the literature, these genes tend to be short, single exon, and lowly expressed. We also find evidence that some of these genes are expressed in other D. melanogaster tissues and both sexes. The relatively small number of intergenic candidate genes discovered here is similar to that observed in the accessory gland, but substantially fewer than that observed in the testis.

Identifying candidate de novo genes expressed in the somatic female reproductive tract of Drosophila melanogaster

Most eukaryotic genes have been vertically transmitted to the present from distant ancestors. However, variable gene number across species indicates that gene gain and loss also occurs. While new genes typically originate as products of duplications and rearrangements of pre-existing genes, putative de novo genes - genes born out of previously non-genic sequence - have been identified. Previous studies of de novo genes in Drosophila have provided evidence that expression in male reproductive tissues is common. However, no studies have focused on female reproductive tissues. Here we begin addressing this gap in the literature by analyzing the transcriptomes of three female reproductive tract organs (spermatheca, seminal receptacle, and parovaria) in three species - our focal species, D. melanogaster - and two closely related species, D. simulans and D. yakuba , with the goal of identifying putative D. melanogaster -specific de novo genes expressed in these tissues. We discovered several candidate genes, which, consistent with the literature, tend to be short, simple, and lowly expressed. We also find evidence that some of these genes are expressed in other D. melanogaster tissues and both sexes. The relatively small number of candidate genes discovered here is similar to that observed in the accessory gland, but substantially fewer than that observed in the testis.

Dosage balance acts as a time-dependent selective barrier to subfunctionalization

BACKGROUND

Gene duplication is an important process for genome expansion, sometimes allowing for new gene functions to develop. Duplicate genes can be retained through multiple processes, either for intermediate periods of time through processes such as dosage balance, or over extended periods of time through processes such as subfunctionalization and neofunctionalization.

RESULTS

Here, we built upon an existing subfunctionalization Markov model by incorporating dosage balance to describe the interplay between subfunctionalization and dosage balance to explore selective pressures on duplicate copies. Our model incorporates dosage balance using a biophysical framework that penalizes the fitness of genetic states with stoichiometrically imbalanced proteins. These imbalanced states cause increased concentrations of exposed hydrophobic surface areas, which cause deleterious mis-interactions. We draw comparison between our Subfunctionalization + Dosage-Balance Model (Sub + Dos) and the previous Subfunctionalization-Only (Sub-Only) Model. This comparison includes how the retention probabilities change over time, dependent upon the effective population size and the selective cost associated with spurious interaction of dosage-imbalanced partners. We show comparison between Sub-Only and Sub + Dos models for both whole-genome duplication and small-scale duplication events.

CONCLUSION

These comparisons show that following whole-genome duplication, dosage balance serves as a time-dependent selective barrier to the subfunctionalization process, by causing an overall delay but ultimately leading to a larger portion of the genome retained through subfunctionalization. This higher percentage of the genome that is ultimately retained is caused by the alternative competing process, nonfunctionalization, being selectively blocked to a greater extent. In small-scale duplication, the reverse pattern is seen, where dosage balance drives faster rates of subfunctionalization, but ultimately leads to a smaller portion of the genome retained as duplicates. This faster rate of subfunctionalization is because the dosage balance of interacting gene products is negatively affected immediately after duplication and the loss of a duplicate restores the stoichiometric balance. Our findings provide support that the subfunctionalization of genes that are susceptible to dosage balance effects, such as proteins involved in complexes, is not a purely neutral process. With stronger selection against stoichiometrically imbalanced gene partners, the rates of subfunctionalization and nonfunctionalization slow; however, this ultimately leads to a greater proportion of subfunctionalized gene pairs.

Experimental characterization of de novo proteins and their unevolved random-sequence counterparts

De novo gene emergence provides a route for new proteins to be formed from previously non-coding DNA. Proteins born in this way are considered random sequences and typically assumed to lack defined structure. While it remains unclear how likely a de novo protein is to assume a soluble and stable tertiary structure, intersecting evidence from random sequence and de novo-designed proteins suggests that native-like biophysical properties are abundant in sequence space. Taking putative de novo proteins identified in human and fly, we experimentally characterize a library of these sequences to assess their solubility and structure propensity. We compare this library to a set of synthetic random proteins with no evolutionary history. Bioinformatic prediction suggests that de novo proteins may have remarkably similar distributions of biophysical properties to unevolved random sequences of a given length and amino acid composition. However, upon expression in vitro, de novo proteins exhibit moderately higher solubility which is further induced by the DnaK chaperone system. We suggest that while synthetic random sequences are a useful proxy for de novo proteins in terms of structure propensity, de novo proteins may be better integrated in the cellular system than random expectation, given their higher solubility.

Binary Bamboo Forest Growth Optimization Algorithm for Feature Selection Problem

Inspired by the bamboo growth process, Chu et al. proposed the Bamboo Forest Growth Optimization (BFGO) algorithm. It incorporates bamboo whip extension and bamboo shoot growth into the optimization process. It can be applied very well to classical engineering problems. However, binary values can only take 0 or 1, and for some binary optimization problems, the standard BFGO is not applicable. This paper firstly proposes a binary version of BFGO, called BBFGO. By analyzing the search space of BFGO under binary conditions, the new curve V-shaped and Taper-shaped transfer function for converting continuous values into binary BFGO is proposed for the first time. A long-mutation strategy with a new mutation approach is presented to solve the algorithmic stagnation problem. Binary BFGO and the long-mutation strategy with a new mutation are tested on 23 benchmark test functions. The experimental results show that binary BFGO achieves better results in solving the optimal values and convergence speed, and the variation strategy can significantly enhance the algorithm's performance. In terms of application, 12 data sets derived from the UCI machine learning repository are selected for feature-selection implementation and compared with the transfer functions used by BGWO-a, BPSO-TVMS and BQUATRE, which demonstrates binary BFGO algorithm's potential to explore the attribute space and choose the most significant features for classification issues.

The nearly complete assembly of the Cercis chinensis genome and Fabaceae phylogenomic studies provide insights into new gene evolution

Fabaceae is a large family of angiosperms with high biodiversity that contains a variety of economically important crops and model plants for the study of biological nitrogen fixation. Polyploidization events have been extensively studied in some Fabaceae plants, but the occurrence of new genes is still concealed, owing to a lack of genomic information on certain species of the basal clade of Fabaceae. Cercis chinensis (Cercidoideae) is one such species; it diverged earliest from Fabaceae and is essential for phylogenomic studies and new gene predictions in Fabaceae. To facilitate genomic studies on Fabaceae, we performed genome sequencing of C. chinensis and obtained a 352.84 Mb genome, which was further assembled into seven pseudochromosomes with 30 612 predicted protein-coding genes. Compared with other legume genomes, that of C. chinensis exhibits no lineage-specific polyploidization event. Further phylogenomic analyses of 22 legumes and 11 other angiosperms revealed that many gene families are lineage specific before and after the diversification of Fabaceae. Among them, dozens of genes are candidates for new genes that have evolved from intergenic regions and are thus regarded as de novo-originated genes. They differ significantly from established genes in coding sequence length, exon number, guanine-cytosine content, and expression patterns among tissues. Functional analysis revealed that many new genes are related to asparagine metabolism. This study represents an important advance in understanding the evolutionary pattern of new genes in legumes and provides a valuable resource for plant phylogenomic studies.

Genome-wide identification and functional analysis of silicon transporter family genes in moso bamboo (Phyllostachys edulis)

Silicon (Si) has crucial effects on plant development and stress resistance. Silicon transporters regulate Si absorption, transport, and distribution in plants. In this study, we identified and characterized the Si transporter gene family of moso bamboo (Phyllostachys edulis) and cloned seven putative Si transporter genes. Moso bamboo Si transporters contain conserved functional domains that mediate the accumulation of considerable amounts of Si. The analysis of gene duplication patterns and divergence times suggested that the expansion of the moso bamboo Si transporter family was mainly due to segmental duplications. The expression of moso bamboo Si transporter genes, which varied among organs, was significantly modulated by Si treatments. The subcellular localization analysis showed that Si transporters are plasma membrane proteins. The Si content increased in transgenic Arabidopsis overexpressing PeLsi1-1 or PeLsi1-2, which affected vegetative and reproductive growth. Our single-particle tracking analysis revealed the four diffusion modes of PeLsi1-1 on the plasma membrane. Moreover, the particle velocity, dwell time, and motion range of PeLsi1-1 decreased in response to Si treatments. The results of this study will further clarify the molecular mechanisms underlying Si absorption and accumulation in bamboo plants.

WGDTree: a phylogenetic software tool to examine conditional probabilities of retention following whole genome duplication events

BACKGROUND

Multiple processes impact the probability of retention of individual genes following whole genome duplication (WGD) events. In analyzing two consecutive whole genome duplication events that occurred in the lineage leading to Atlantic salmon, a new phylogenetic statistical analysis was developed to examine the contingency of retention in one event based upon retention in a previous event. This analysis is intended to evaluate mechanisms of duplicate gene retention and to provide software to generate the test statistic for any genome with pairs of WGDs in its history.

RESULTS

Here a software package written in Python, 'WGDTree' for the analysis of duplicate gene retention following whole genome duplication events is presented. Using gene tree-species tree reconciliation to label gene duplicate nodes and differentiate between WGD and SSD duplicates, the tool calculates a statistic based upon the conditional probability of a gene duplicate being retained after a second whole genome duplication dependent upon the retention status after the first event. The package also contains methods for the simulation of gene trees with WGD events. After running simulations, the accuracy of the placement of events has been determined to be high. The conditional probability statistic has been calculated for Phalaenopsis equestris on a monocot species tree with a pair of consecutive WGD events on its lineage, showing the applicability of the method.

CONCLUSIONS

A new software tool has been created for the analysis of duplicate genes in examination of retention mechanisms. The software tool has been made available on the Python package index and the source code can be found on GitHub here: https://github.com/cnickh/wgdtree .

The genomic basis of evolutionary novelties in a leafhopper

Evolutionary innovations generate phenotypic and species diversity. Elucidating the genomic processes underlying such innovations is central to understanding biodiversity. In this study, we addressed the genomic basis of evolutionary novelties in the glassy-winged sharpshooter (Homalodisca vitripennis, GWSS), an agricultural pest. Prominent evolutionary innovations in leafhoppers include brochosomes, proteinaceous structures that are excreted and used to coat the body, and obligate symbiotic associations with two bacterial types that reside within cytoplasm of distinctive cell types. Using PacBio long-read sequencing and Dovetail Omni-C technology, we generated a chromosome-level genome assembly for the GWSS and then validated the assembly using flow cytometry and karyotyping. Additional transcriptomic and proteomic data were used to identify novel genes that underlie brochosome production. We found that brochosome-associated genes include novel gene families that have diversified through tandem duplications. We also identified the locations of genes involved in interactions with bacterial symbionts. Ancestors of the GWSS acquired bacterial genes through horizontal gene transfer (HGT), and these genes appear to contribute to symbiont support. Using a phylogenomics approach, we inferred HGT sources and timing. We found that some HGT events date to the common ancestor of the hemipteran suborder Auchenorrhyncha, representing some of the oldest known examples of HGT in animals. Overall, we show that evolutionary novelties in leafhoppers are generated by the combination of acquiring novel genes, produced both de novo and through tandem duplication, acquiring new symbiotic associations that enable use of novel diets and niches, and recruiting foreign genes to support symbionts and enhance herbivory.

The De Novo Genome Assembly of Olea europaea subsp. cuspidate, a Widely Distributed Olive Close Relative

The olive complex, comprising six subspecies, is a valuable plant for global trade, human health, and food safety. However, only one subspecies (Olea europaea subsp. europaea, OE) and its wild relative (Olea europaea subsp. europaea var. sylvestris, OS) have genomic references, hindering our understanding of the evolution of this species. Using a hybrid approach by incorporating Illumina, MGI, Nanopore, and Hi-C technologies, we obtained a 1.20-Gb genome assembly for the olive subspecies, Olea europaea subsp. cuspidate (OC), with contig and scaffold N50 values of 5.33 and 50.46 Mb, respectively. A total of 43,511 protein-coding genes were predicted from the genome. Interestingly, we observed a large region (37.5 Mb) of “gene-desert” also called “LTR-hotspot” on chromosome 17. The gene origination analyses revealed a substantial outburst (19.5%) of gene transposition events in the common ancestor of olive subspecies, suggesting the importance of olive speciation in shaping the new gene evolution of OC subspecies. The divergence time between OC and the last common ancestor of OE and OS was estimated to be 4.39 Mya (95% CI: 2.58–6.23 Mya). The pathways of positively selected genes of OC are related to the metabolism of cofactors and vitamins, indicating the potential medical and economic values of OC for further research and utilization. In summary, we constructed the de novo genome assembly and protein-coding gene pool for Olea europaea subsp. cuspidate (OC) in this study, which may facilitate breeding applications of improved olive varieties from this widely distributed olive close relative.

Single-base methylome analysis reveals dynamic changes of genome-wide DNA methylation associated with rapid stem growth of woody bamboos

CG and CHG methylation levels in the rapid shoot growth stages (ST2-ST4) of woody bamboos were obviously decreased, which might regulate the internode elongation during rapid shoot growth, while CHH methylation was strongly associated with shoot developmental time or age. DNA methylation plays a critical role in the regulation of plant growth and development. Woody bamboos have a unique trait of rapid stem growth resulted from internode elongation at the shooting period. However, it is still unclear whether DNA methylation significantly controls the bamboo rapid stem growth. Here we present whole-genome DNA methylation profiles of the paleotropical woody bamboo Bonia amplexicaulis at five newly defined stages of shoot growth, named ST1-ST5. We found that CG and CHG methylation levels in the rapid shoot growth stages (ST2-ST4) were significantly lower than in the incubation (ST1) and plateau stages (ST5). The changes in methylation levels mainly occurred in flanking regions of genes and gene body regions, and 23647 differentially methylated regions (DMRs) were identified between ST1 and rapid shoot growth stages (ST2-ST4). Combined with transcriptome analysis, we found that DMR-related genes enriched in the auxin and jasmonic acid (JA) signal transduction, and other pathways closely related to plant growth. Intriguingly, CHH methylation was not involved in the rapid shoot growth, but strongly associated with shoot developmental time by gradually accumulating in transposable elements (TEs) regions. Overall, our results reveal the importance of DNA methylation in regulating the bamboo rapid shoot growth and suggest a role of DNA methylation associated with development time or age in woody bamboos.

Anatomical Characteristics and Variation Mechanisms on the Thick-Walled and Dwarfed Culm of Shidu Bamboo (Phyllostachys nidularia f. farcta)

Stable culm variants are valuable and important material for the study of culm development in bamboo plants. However, to date, there are few reports on the mechanism of variation of these bamboo variants. Phyllostachys nidularia f. farcta (Shidu bamboo) is a bamboo variant with stable phenotypes such as a dwarf culm with a thickened wall. In this study, we systematically investigated the cytological characteristics and underlying mechanism of morphological variation in culms of this variant using anatomical, mathematical statistical, physiological, and genomic methods. The anatomical observation and statistical results showed that the lateral increase of ground tissue in the inner layer of culm wall and the enlargement of vascular bundles are the anatomical essence of the wall thickening of Shidu bamboo; the limited elongation of fiber cells and the decrease in the number of parenchyma cells longitudinally are probably the main causes of the shortening of its internodes. A number of genes involved in the gibberellin synthesis pathway and in the synthesis of cell wall components are differentially expressed between the variant and its prototype, Ph. nidularia, and may play an important role in determining the phenotype of internode shortening in Shidu bamboo. The decrease in gibberellin content and the content of the major chemical components of the cell wall of Shidu bamboo confirmed the results of the above transcriptome. In addition, the variation in culm morphology in Shidu bamboo had little effect on the volume of the culm wall of individual internodes, suggesting that the decrease in the total number of internodes and the decrease in dry matter content (lignin, cellulose, etc.) may be the main factor for the sharp decline in culm biomass of Shidu bamboo.